也许像这样?
s = Sum(Indexed('x',i),(i,1,3))
f = lambda x: Subs(s.doit(), [s.function.subs(s.variables[0], j)
for j in range(s.limits[0][1], s.limits[0][2] + 1)], x).doit()
>>> f((30,10,2))
42
Run Code Online (Sandbox Code Playgroud)
您可以使用lambdaify。只需确保总和的限制与 numpy 数组的可迭代项匹配即可。
from sympy import Sum, symbols, Indexed, lambdify
import numpy as np
x, i = symbols("x i")
s = Sum(Indexed('x',i),(i,0,3))
f = lambdify(x, s)
b = np.array([1, 2, 3, 4])
f(b)
Run Code Online (Sandbox Code Playgroud)