如何用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}]

如果你传递一堆符号表达式，则在solve隐含地理解等于0。

请注意，您的系统不是线性的：有x2的除法和x2的平方根。