python 解方程

最新推荐文章于 2025-07-04 01:30:00 发布

yichudu

最新推荐文章于 2025-07-04 01:30:00 发布

阅读量3.5k

点赞数

CC 4.0 BY-SA版权

分类专栏： python 文章标签：方程

天天开心

本文链接：https://blog.csdn.net/chuchus/article/details/75356014

python 专栏收录该内容

54 篇文章

订阅专栏

本文介绍使用Sympy与Scipy两个Python库进行方程求解的方法。Sympy适用于符号计算，能高效解决线性和非线性方程；Scipy则通过迭代法提供数值解，适合解决复杂的多元方程组。文章提供了具体的代码实例。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

我用到的有两个库. sympy 与 scipy.

sympy

Symbolic Computation,符号计算.
一个符号计算库. anaconda自带. 解方程的功能比 scipy 强大.

线性方程

solve the equation $-3x^2+2x=0$ . The answer is $x_1=0,x_2=2/3$ .

from sympy import *
x=symbols('x')  # indicate x is an symbol
eq=Eq(-3*x**2+2*x,0) # type the equation
print(solve(eq,x))  # [0, 2/3]

非线性方程

def get_point_of_contact(O: Point, A: Point, r: float)->str:
    """
    过圆外一点A, 求圆O的切点
    :param O: 圆心
    :param A: 圆外一点
    :param r: 圆的半径
    :return: 切点坐标
    """
    from sympy import symbols,nonlinsolve

    x,y=symbols('x y')  # indicate x is an symbol
    # 二元二次方程组:
    # * 切线与垂线向量乘法
    # * 圆的方程代入
    result=nonlinsolve([(A.x - x) * (x - O.x) + (A.y - y) * (y - O.y),
            (x - O.x) * (x - O.x) + (y - O.y) * (y - O.y) - r * r],[x,y])
    return result
print(get_point_of_contact(Point(0, 0), Point(2, 0), 1))
"""{(1/2, -sqrt(3)/2), (1/2, sqrt(3)/2)}"""