Skip to content
Advertisement

Sympy multiplies answer changing it and making it wrong

I am trying to find the inverse of a matrix, however, sympy always multiplies the final answer by 4, thereby making the answer wrong. This multiplication also happens when I use the factor function.

Here is the matrix I want to change


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

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


And here is the “wrong” answer it gives enter image description here

The true answer is the same only NOT multiplied by 4. Why does sympy multiply the answer like this?

Update

To isolate the problem we can try to factor the determinant.

If we factor the determinant I expect the answer to be(z+1/2)**2 However this is the answer I get, it is divided by 4 and , multiplied by 2. The expected answer has to be exactly (z+1/2)^2 not the multiplied version sympy outputs.


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)


enter image description here

The problem I am solving evolves the Z-transfrom and this multiplication of the factor is making the final answer wrong. The factors or poles has to be exactly what they are without any multiplication, ie the answer has to be exactly (z+1/2)^2

Is there a way to make sympy give me the exact factors?

Advertisement

Answer

Be careful with your words: there is nothing “wrong” about the answer that you get. You would prefer to have the answer in a different form and that’s fine but the answer given is correct.

You can get the form that you want if you factor over the Gaussian rationals:

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) ⎦

User contributions licensed under: CC BY-SA
6 People found this is helpful
Advertisement