0

我试图找到矩阵的逆矩阵，但是，sympy 总是将最终答案乘以 4，从而使答案错误。当我使用因子函数时，也会发生这种乘法。

这是我要更改的矩阵


import sympy as sy
z = sy.symbols('z')

t = sy.Matrix([[z, -1], [sy.Rational(1,4), z + 1]])
t.inv()

这是它给出的“错误”答案 在此处输入图像描述

真正的答案是相同的，只是不乘以 4。为什么 sympy 会像这样乘以答案？

更新

为了隔离问题，我们可以尝试分解行列式。

如果我们考虑行列式，我希望答案是 (z+1/2)**2 但是这是我得到的答案，它除以 4 和 ，再乘以 2。预期的答案必须是 (z+ 1/2)^2 不是相乘版本的 sympy 输出。


d = t.det()
print("I want to factor this")
display(d)

print("After factoring")
display(sy.factor(d))

print(" \n The right answer I was expecting")
display((z+sy.Rational(1,2))**2)

在此处输入图像描述

我正在解决的问题演变了 Z 变换，并且因子的这种乘法使最终答案错误。因子或极点必须与它们完全相同，无需任何乘法，即答案必须完全为 (z+1/2)^2

有没有办法让 sympy 给我确切的因素？

4

1 回答 1

1

小心你的话：你得到的答案没有任何“错误”。您希望得到不同形式的答案，这很好，但给出的答案是正确的。

如果你考虑高斯有理数，你可以得到你想要的形式：

In [41]: e = z**2 + z + Rational(1, 4)                                                                                                         

In [42]: e                                                                                                                                     
Out[42]: 
 2       1
z  + z + ─
         4

In [43]: factor(e)                                                                                                                             
Out[43]: 
         2
(2⋅z + 1) 
──────────
    4     

In [44]: factor(e, gaussian=True)                                                                                                                
Out[44]: 
         2
(z + 1/2) 

In [45]: t.inv()                                                                                                                               
Out[45]: 
⎡    4⋅z + 4              -1        ⎤
⎢ ──────────────   ─────────────────⎥
⎢    2               ⎛   2         ⎞⎥
⎢ 4⋅z  + 4⋅z + 1     ⎜  z    z   1 ⎟⎥
⎢                  4⋅⎜- ── - ─ - ──⎟⎥
⎢                    ⎝  4    4   16⎠⎥
⎢                                   ⎥
⎢       1                4⋅z        ⎥
⎢────────────────   ──────────────  ⎥
⎢  ⎛   2       1⎞      2            ⎥
⎢4⋅⎜- z  - z - ─⎟   4⋅z  + 4⋅z + 1  ⎥
⎣  ⎝           4⎠                   ⎦

In [46]: t.inv().applyfunc(lambda e: factor(e, gaussian=True))                                                                                   
Out[46]: 
⎡   z + 1          1     ⎤
⎢ ──────────   ──────────⎥
⎢          2            2⎥
⎢ (z + 1/2)    (z + 1/2) ⎥
⎢                        ⎥
⎢    -1            z     ⎥
⎢────────────  ──────────⎥
⎢           2           2⎥
⎣4⋅(z + 1/2)   (z + 1/2) ⎦
于 2020-10-16T16:02:03.667 回答