您好, 欢迎来到 !    登录 | 注册 | | 设为首页 | 收藏本站

python – SymPy和复数的平方根

5b51 2022/1/14 8:21:48 python 字数 1865 阅读 521 来源 www.jb51.cc/python

当使用求解来计算二次方程的根时,SymPy返回可以简化的表达式,但我无法简化它们.最小的例子如下:from sympy import * sqrt(-24-70*I) 在这里,SymPy只返回sqrt(-24-70 * I),而Mathematica或Maple将回复相当于5-7 * I.我知道有两个平方根,但这种行为需要SymPy,例如,将返回相当复杂的解

概述

当使用求解来计算二次方程的根时,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将很高兴地给我两个高斯整数来解决这个等式.

有没有选项或我缺少的东西?

from sympy import *
z = -24-70*I
x,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*y
solve(z ** 2 + (1 + I) * z + (6 + 18 * I),y))

输出:[( – 3,3),(2,-4)]

总结

以上是编程之家为你收集整理的python – SymPy和复数的平方根全部内容,希望文章能够帮你解决python – SymPy和复数的平方根所遇到的程序开发问题。


如果您也喜欢它,动动您的小指点个赞吧

除非注明,文章均由 laddyq.com 整理发布,欢迎转载。

转载请注明:
链接:http://laddyq.com
来源:laddyq.com
著作权归作者所有。商业转载请联系作者获得授权,非商业转载请注明出处。


联系我
置顶