Sympy积分特定参数正态分布失败及绘图报错问题咨询

结论 #

这是SymPy符号积分模块的内置启发式匹配逻辑导致的功能限制，不属于使用错误。

差异原因 #

SymPy的sympy.integrate对正态分布这类有已知闭合解的初等函数积分，会优先尝试对输入的参数做简化匹配：

• 第一组参数10.1、0.333333333可被快速识别为简单有理数近似值（101/10、1/3），命中了正态分布积分的预定义求解规则，因此返回含erf函数的闭合解析结果。
• 第二组参数8.655555555555557、0.5212875796916135的浮点精度更高，无法被快速映射为低复杂度有理数，没有触发对应的积分匹配规则，因此直接返回未计算的Integral表达式对象。

报错原因 #

SymPy的绘图接口会自动将输入的符号表达式转换为可执行的数值计算函数，未计算的Integral对象在转换时会直接将字面量写入生成的代码，运行时自然找不到对应的Integral定义，因此抛出NameError。

修复方案 #

方案1：参数预简化 #

积分前调用sympy.nsimplify对浮点参数做简化处理，触发积分匹配规则：

import sympy as sp
x = sp.symbols('x')
# 预简化参数
mu = sp.nsimplify(8.655555555555557)
sigma = sp.nsimplify(0.5212875796916135)
pdf = 1/(sigma * sp.sqrt(2*sp.pi)) * sp.exp(-(x - mu)**2 / (2 * sigma**2))
# 可正常得到含erf的解析结果
integral_res = sp.integrate(pdf, x)

方案2：直接使用数值积分 #

如果不需要解析解仅用于绘图，可直接对Integral对象调用evalf做数值计算，或者改用matplotlib手动生成数值点绘图。

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