相关疑难解决方法(0)

我有2个矩阵,其中第一个是稀疏的整数系数.

import sympy
A = sympy.eye(2)
A.row_op(1, lambda v, j: v + 2*A[0, j])

第二个是象征性的,我在它们之间执行操作:

M = MatrixSymbol('M', 2, 1)
X = A * M + A.col(1)

现在,我想要的是获得元素方程式:

X_{0,0} = A_{0,0}
X_{0,1} = 2*A_{0,0} + A_{0,1}

一种方法是指定一个矩阵sympy,每个元素都是一个单独的符号:

rows = []
for i in range(shape[0]):
    col = []
    for j in range(shape[1]):
        col.append(Symbol('%s_{%s,%d}' % (name,i,j)))
    rows.append(col)
M = sympy.Matrix(rows)

有没有办法用MatrixSymbol上面的方法来做,然后得到由此产生的元素方程式？