如何用sympy求解矩阵方程？

Question

同情,给出一个矩阵方程

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

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

我用普通符号尝试了这个,但显然失败了.使用MatrixSymbol也不行.有没有办法做到这一点,或者是否有能力做到这一点？

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实现这一点的基本内容已经存在.它可能很容易实现.