无法完成多运算符求和：NumPy ufunc 'add'类型转换错误求助

你遇到的错误核心是初始累加变量的类型不匹配：用it = 0（Python原生整数）作为累加起点，第一次循环时会和numpy数组（np.cos(theta)的结果）相加，此时it会被强制转换成object dtype的数组。当m增大、循环次数变多后，后续计算的数组无法和object dtype的数组完成同类型相加，触发类型转换错误。

修复步骤 #

把累加变量it初始化为和theta同类型的numpy数组，而非原生整数：

from scipy.special import factorial2
import numpy as np  # 补充缺失的numpy导入

m = 5

theta = np.linspace(0, 1.5, 100)

def I(n):
    it = np.zeros_like(theta)  # 初始化为和theta同形状、同类型的0数组
    for t in range(n+1):
        term = ((2**t)*(m**t)*(np.cos(theta)**(2*t)))/factorial2(2*t+1)
        it += term
    return it

Itheta= I(20)

额外说明 #

1. 原代码遗漏了numpy的导入语句，运行时会触发NameError，已补充完整。
2. 当m较大、循环次数n过高时，(2*m)**t会指数级增长，可能出现数值溢出变为inf，这属于数值稳定性问题，而非当前的类型错误。若要优化，可考虑用对数计算再转换回原值，或改用np.float128类型提升精度。

内容的提问来源于stack exchange，提问作者augusto