如何用Sympy求解矩陣方程（矩陣= 0）？

Question

如果系統被寫爲X矩陣，如何用Sympy求解方程組？

我有這樣的代碼

from sympy import * 
x1, x2 = symbols('x1 x2') 
u1 = 0.0005 u2 = 0.0005 
K = sqrt(2)*10^-3 
A = 0.2 
X = Matrix([40*u1/(A*x2) + 10*u2/(A*x2) - (x1*K*sqrt(x2))/(A*x2), u1/A + u2/A - K*sqrt(x2)/A]) 

的X是一個矩陣，你可以看到。當X矩陣= 0時，我需要找到x1和x2的值。我將對此矩陣進行線性化。

是的。我想要X矩陣在這個表單上。不是Ax = b形式。

Answer 1

solve(X) 

回報[{x1: 25.0000000000000, x2: 0.500000000000000}]

等同於0隱含理解solve，如果你傳遞了一堆符號表達式。

請注意，您的系統是而不是線性：除以x2和x2的平方根。