
当使用求解来计算二次方程的根时,SymPy返回可以简化的表达式,但我无法简化它们.最小的例子如下:
from sympy import *sqrt(-24-70*I)在这里,SymPy只返回sqrt(-24-70 * I),而Mathematica或Maple将回复相当于5-7 * I.
我知道有两个平方根,但这种行为需要SymPy,例如,将返回相当复杂的解决方案
z = symbols("z")solve(z ** 2 + (1 + I) * z + (6 + 18 * I),z)同时,Maple和Mathematica将很高兴地给我两个高斯整数来解决这个等式.
有没有选项或我缺少的东西?
最佳答案找到z的平方根在逻辑上与求解方程(x I * y)** 2 = z相同.所以你可以做到这一点:from sympy import *z = -24-70*Ix,y = symbols('x y',real=True)result = solve((x+I*y)**2 - z,(x,y))结果是[(-5,7),(5,-7)]
为方便起见,这可以包装为一个功能:
def my_sqrt(z): x,real=True) sol = solve((x+I*y)**2 - z,y)) return sol[0][0] + sol[0][1]*I现在你可以使用my_sqrt(-24-70 * I)并得到-5 7 * I.
在您的示例中,相同的策略有助于使用二次方程:
x,real=True)z = x + I*ysolve(z ** 2 + (1 + I) * z + (6 + 18 * I),y))输出:[( – 3,3),(2,-4)] 总结
以上是内存溢出为你收集整理的python – SymPy和复数的平方根全部内容,希望文章能够帮你解决python – SymPy和复数的平方根所遇到的程序开发问题。
如果觉得内存溢出网站内容还不错,欢迎将内存溢出网站推荐给程序员好友。
欢迎分享,转载请注明来源:内存溢出
微信扫一扫
支付宝扫一扫
评论列表(0条)