首页 文章

在sympy中用符号索引矩阵

提问于
浏览
0

我尝试在这样的求和中索引矩阵

from sympy import *
vx1,vx2,vx3,vx4,vx5, vy1,vy2,vy3,vy4,vy5,  = symbols('vx1 vx2 vx3 vx4 vx5 vy1 vy2 vy3 vy4 vy5')
vx=Matrix([vx1,vx2,vx3,vx4,vx5])
vy=Matrix([vy1,vy2,vy3,vy4,vy5])
p, n = symbols('p n', integer=True)
vx[0]
vx[1]
vx[2]
vx[3]
summation(p, (p, 0, 4))
summation(vx[p], (p, 0, 4))

但似乎同情不能这样做:

NameError: IndexError: Invalid index a[p]

有办法吗?

2 回答

  • 0

    如果要将符号索引转换为Matrix,请使用MatrixSymbol:

    In [15]: vx = MatrixSymbol('vx', 1, 4)
    
    In [16]: summation(vx[(0, p)], (p, 0, 4)).doit()
    Out[16]: vx₀₀ + vx₀₁ + vx₀₂ + vx₀₃ + vx₀₄
    
  • 4

    以下怎么样?

    >>> sum(vx[p] for p in range(5))
    vx1 + vx2 + vx3 + vx4 + vx5
    

相关问题