Как решить матричное уравнение (matrix = 0) с помощью Sympy?

Question

Как мне решить систему уравнений с помощью Sympy, если система записана в виде X-матрицы?

У меня есть этот код

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 - это матрица. Мне нужно найти значения для x1 и x2, когда матрица X = 0. Я собираюсь линеаризовать эту матрицу.

Да. Я хочу, чтобы матрица Х была в этой форме. Не Ax = b форма.

Answer 1

solve(X)

возвращается [{x1: 25.0000000000000, x2: 0.500000000000000}]

Приравнивание к 0 неявно понимается в solve, если вы передадите ему кучу символических выражений.

Обратите внимание, что ваша система не является линейной: есть деление на x2 и квадратный корень из x2.