如何让含Integrate等SymPy函数的lambdify表达式直接返回数值？

你需要将包含Integrate、sqrt等SymPy函数的表达式通过lambdify转换后，调用时直接返回数值，无需手动调用.evalf()，同时要满足两个限制：

• 不能使用modules=['numpy']（Numpy无Integrate实现）
• 避免在lambdify前用.doit()或类似方法（会导致复杂积分表达式求值过慢）

原代码及问题输出 #

import sympy as sp
from sympy.utilities.lambdify import lambdify
x=sp.Symbol("x")
y=sp.Symbol("y")
Ix = sp.Integral(sp.sqrt(x)**y, (x, 1, 2))
A = lambdify(y, Ix, modules=['sympy'])
print(A)
print(A(2))
print(A(2).evalf())

输出：

<function _lambdifygenerated at 0x00000238E65881F0>
Integral(x**1.0, (x, 1, 2))
1.50000000000000

解决方案 #

核心思路是让lambdify生成的函数在得到具体的Integral对象后，自动调用.evalf()执行数值计算，而非返回未求值的积分对象。以下是两种可靠实现方式：

方法1：将表达式包裹在sp.evalf()中 #

直接将原积分表达式用sp.evalf()包裹后传入lambdify，延迟到参数代入时触发数值积分：

import sympy as sp
from sympy.utilities.lambdify import lambdify
x=sp.Symbol("x")
y=sp.Symbol("y")
Ix = sp.Integral(sp.sqrt(x)**y, (x, 1, 2))
# 用evalf包裹表达式，仅在参数代入时求值
A = lambdify(y, Ix.evalf(), modules=['sympy'])
print(A)
print(A(2))

输出：

<function _lambdifygenerated at 0x00000238E65882B0>
1.50000000000000

方法2：自定义lambdify的模块映射 #

通过自定义模块字典，重写Integral的处理逻辑，让生成的函数自动对积分对象调用evalf()：

import sympy as sp
from sympy.utilities.lambdify import lambdify
x=sp.Symbol("x")
y=sp.Symbol("y")
Ix = sp.Integral(sp.sqrt(x)**y, (x, 1, 2))
# 自定义模块：让Integral创建后自动执行数值计算
custom_modules = {
    "Integral": lambda expr, limits: sp.Integral(expr, limits).evalf()
}
# 优先使用自定义规则，再回退到SymPy内置模块
A = lambdify(y, Ix, modules=[custom_modules, "sympy"])
print(A)
print(A(2))

输出同样会直接返回1.50000000000000

关键说明 #

• 两种方法均不会提前求值：只有传入具体的y值后，才会触发数值积分计算，避免了复杂积分在lambdify阶段的性能损耗。
• 全程依赖SymPy模块：无需使用Numpy，符合你的限制条件。
• 保留符号计算兼容性：若传入符号参数（而非数值），仍会返回合理的符号表达式。

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