Scipy integrate.quad循环报错：无法将表达式转换为浮点数

问题原因 #

你遇到的错误来自三个核心问题：

1. 闭包变量延迟绑定：循环内定义的函数f会捕获循环变量i，但quad是延迟调用f，当f实际执行时，i已经变成循环的最后一个值（2），而非循环时的当前i值。
2. 表达式类型不匹配：你提前用未定义的x（或符号变量x）定义了Ni_1D和dN1D数组，导致f返回的是符号表达式而非数值，而quad要求输入函数必须接受数值并返回数值。
3. 导数笔误：原代码中dN1D的第二个元素x-2*x是错误的——x-x²的导数应该是1-2x，这个笔误会导致计算结果完全错误。

修复后的代码方案 #

方案1：用默认参数绑定循环变量

import numpy as np
from scipy import integrate

Tn = [10.0, 10.0, 10.0]
Tt = [0, 0, 0]
y_coor = [0.0, 1.25, 2.5]
x_coor = [4.0, 4.0, 4.0]
fx = []

for i in range(3):
    # 提前计算当前i对应的系数
    coeff = Tn[i] * y_coor[i] + Tt[i] * x_coor[i]
    # 用默认参数绑定当前的i和coeff，避免闭包延迟问题
    def f(x, i=i, coeff=coeff):
        # 直接在函数内计算当前i对应的Ni_1D和dN1D值
        if i == 0:
            ni = (x/2)*(x-1)
            dni = x - 0.5
        elif i == 1:
            ni = x - x**2
            dni = 1 - 2*x  # 修正导数笔误
        else:
            ni = (x/2)*(x+1)
            dni = x + 0.5
        return ni * coeff * dni
    # 执行积分
    fx.append(integrate.quad(f, -1, 1))

print(fx)

方案2：用functools.partial绑定参数

如果你更倾向于简洁的写法，可以用functools.partial来绑定循环参数：

import numpy as np
from scipy import integrate
from functools import partial

def calc_integrand(x, i, coeff):
    # 按索引获取当前i对应的基函数和导数
    ni_list = [(x/2)*(x-1), x-x**2, (x/2)*(x+1)]
    dni_list = [x-0.5, 1-2*x, x+0.5]
    return ni_list[i] * coeff * dni_list[i]

Tn = [10.0, 10.0, 10.0]
Tt = [0, 0, 0]
y_coor = [0.0, 1.25, 2.5]
x_coor = [4.0, 4.0, 4.0]
fx = []

for i in range(3):
    coeff = Tn[i] * y_coor[i] + Tt[i] * x_coor[i]
    # 绑定i和coeff到积分函数
    bound_func = partial(calc_integrand, i=i, coeff=coeff)
    fx.append(integrate.quad(bound_func, -1, 1))

print(fx)

验证结果 #

运行上述代码后，你会得到每个i对应的积分结果（积分值+误差估计），比如i=1时的结果和你手动定义函数计算的结果一致。

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