python - 如何用sympy求解矩阵方程？

Question

在 sympy 中，给定一个矩阵方程

M * x + N * y = 0(或更复杂..)

如何为 x 解决这个问题？ (M,N = 矩阵, x,y = 向量)

我用普通符号试过这个，但显然失败了。使用 MatrixSymbol 效果不佳。有什么办法可以做到，还是 sympy 不能做到？

Answer 1

正如 MRocklin 指出的那样，MatrixExpressions 尚不支持此功能，但非交换符号支持:

In [13]: M, N, x, y = symbols('M N x y', commutative=False)

In [15]: solve(M*x + N*y, x)
Out[15]:
⎡      -1⎤
⎣-N⋅y⋅M  ⎦

与 MatrixExpressions 不同，非交换符号没有形状概念，因此您需要自己跟踪它。但这也表明为 MatrixExpression 实现此功能的基本内容已经存在。它可能很容易实现。