1
で行列の配列は:一覧/私はsympyに次の式を作成したいsympy
M
、U
とA^(k)
はサイズm
の正方行列です。
はこれまでのところ、私は与えられたp
のためにこれを行うために管理:
import sympy
m = sympy.symbols('m', integer=True)
p = 3
A = [sympy.MatrixSymbol('A^({})'.format(k), m, m) for k in range(p)]
M = sympy.MatrixSymbol('M', m, m)
U = sympy.MatrixSymbol('U', m, m)
expr = M*A[0]*U
for k in range(1, p):
expr += M*A[k]*U
expr
しかし、私は固定数にp
を設定する必要はありません。代わりにp
をシンボルにする必要があります。ちょうどm
と同じです。これはsympyでどのように達成できますか?