如何解决Sympy 解决功能为我提供了比预期更多的解决方案
我正在尝试使用 SymPy 求解矩阵方程。我有一个 4x4 矩阵,通过手工计算我知道我应该得到什么解决方案。我的代码中的 Matrix 's' 等于零,所以我正在解决使用行列式找到 z 的问题。当我尝试在 SymPy 中执行此操作时,我得到 6 个解决方案而不是 4 个,并且解决方案中也包含正确的值。
我也遇到了显示为 P(z)/Q(z) 的表达式的问题,其中 P 和 Q 是 z 的函数。
我的代码如下。
from sympy import *
import time
init_printing()
z = Symbol('z',real=True)
s = Matrix([[10322.4 - (230.2*z),-5161.2,0],[-5161.2,10322.4-(230.2*z),[0,-5161.2],5161.2-(220.5*z)]])
a = s.det()
print(a)
start = time.time()
z_solutions = solve(a,z)
end = time.time()
print("the solutions to the equation are:")
print()
print(z_solutions)
我的输出是:
(70889915769.4571*z**8 - 23910733144618.3*z**7 + 3.37005688026419e+15*z**6 - 2.57290735158039e+17*z**5 + 1.15080765413656e+19*z**4 - 3.0326418528654e+20*z**3 + 4.44856872382415e+21*z**2 - 3.07669487225344e+22*z + 5.67055263945047e+22)/(26.3548111557343*z**4 - 4727.10517244052*z**3 + 304704.23*z**2 - 8316757.68*z + 79913956.32)
the solutions to the equation are:
[2.75414367052697,22.4205039097983,22.7444509610069,53.0547635060136,67.2615117276319,79.3764680453814]
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。