如何让SymPy计算含参数的收敛无穷级数？

我之前也碰到过类似的情况，这其实是SymPy在处理符号参数的无穷级数时，对参数的约束条件有更严格的要求导致的。

问题根源 #

你定义的a仅指定了real=True，但SymPy需要明确这个参数是非零且非整数的实数，才能触发对应的收敛级数求和公式——毕竟如果a是整数，级数里会出现分母为0的项直接发散；如果a=0，级数就变成调和级数的平方，同样发散。没有这些明确约束，SymPy不敢贸然应用求和规则，所以只会返回未计算的级数表达式。

而当你换成exp(1)这类具体值时，SymPy能直接判断它是非零非整数的实数，满足收敛条件，所以自动调用了对应的求和逻辑，返回了结果。

解决办法 #

给符号a添加更精确的约束条件即可，有两种常用方式：

方式1：定义符号时直接指定约束

from sympy import Symbol, Sum, pi, coth

# 定义a为非整数的实数，同时隐含a≠0（整数包含0）
a = Symbol('a', real=True, noninteger=True)
result = Sum(1/(n**2 + a**2), (n, 1, oo)).doit()
print(result)
# 输出：(π*coth(π*a))/(2*a) - 1/(2*a**2)

方式2：用refine函数事后补充约束

如果需要保留符号的灵活性，也可以先求和，再用refine函数明确参数满足的条件：

from sympy import Symbol, Sum, pi, coth, refine, Q

a = Symbol('a', real=True)
# 先尝试求和，此时会返回未计算的级数或带条件的表达式
raw_result = Sum(1/(n**2 + a**2), (n, 1, oo)).doit()
# 明确a是非零且非整数的实数
refined_result = refine(raw_result, Q.nonzero(a) & Q.noninteger(a))
print(refined_result)
# 同样得到：(π*coth(π*a))/(2*a) - 1/(2*a**2)

验证结果 #

你可以替换具体值验证正确性，比如把a=2代入结果：

print(refined_result.subs(a, 2))
# 输出：(π*coth(2*π))/4 - 1/8

直接用具体值求和也会得到相同结果，说明这个方法是可靠的。

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