Python中sympy.solve()函数的使用

1、函数的简单使用
　　用 sympy.solve() 函数来解单一方程的调用方式可以表示为 sympy.solve(表达式,符号) 。下面的代码演示了怎样解exp ⁡ ( x ) − 3 = 0 \exp(x)-3=0exp(x)−3=0这个方程。

import sympy
from sympy import *
d = sympy.Symbol("d")
b = sympy.solve(sympy.exp(d)-3,d)
print(b)

可以得到log3的精确表达式。

2、遇见的问题

　　问题1 TypeError: can’t convert expression to float报错

　　该问题是将 math 库中的运算符用在了表达式中，如log ⁡ \loglog,sin ⁡ \sinsin,exp ⁡ \expexp等。应该用 sympy 库中的运算符，可以在方程表达式中用 sympy.log() 语句指定，或者去掉开头的 import math 语句。

　　问题2 No algorithms are implemented to solve equation报错

如下面的代码所示，我们尝试解  x + exp ⁡ ( x ) + sin ⁡ ( x ) − 10 = 0 这个方程。

import sympy
from sympy import *
d = sympy.Symbol("d")
b = sympy.solve(d+sympy.exp(d)+sympy.sin(d)-10,d)
print(b)

结果是报错，没有算法来解这个方程，换种说法就是不能得到解析解。我们可以用数值计算函数 sympy.nsolve() 来近似计算的到数值解。具体用法如下：

import sympy
from sympy import *
d = sympy.Symbol("d")
b=sympy.nsolve(d+sympy.exp(d)+sympy.sin(d)-10,0)
print(b)

运行程序可以得到结果约为1.962，十分精确。

备注：

原文链接：https://blog.csdn.net/weixin_44892951/article/details/106611169