如何用Scipy求解复合积分？解决TypeError: must be real number报错

问题描述 #

需要计算两个复合积分函数：

f(x) = ∫(px) dx（积分区间0-1）
g(x) = ∫(alpha(ln(f(x)))) dx（积分区间0-10）

编写的代码运行报错：

TypeError: must be real number, not function

原代码：

import numpy as np
import math
from scipy.integrate import quad

def firstMethod(x, p1):
    return p1 * x

def firstIntegral():
    return quad(firstMethod, 0, 1, args=(4))[0]

def secondMethod(x, alpha):
    return alpha*math.log(firstIntegral)

def secondIntegral():
    return quad(secondMethod, 0, 10, args=(0.5))[0]

print(secondIntegral())

错误原因 #

1. 函数调用错误：在secondMethod中，firstIntegral是函数对象，而非函数的返回值。直接传递函数对象给math.log会触发类型错误，必须调用firstIntegral()获取积分结果。
2. 元组参数格式错误：quad的args参数要求传入元组，args=(4)实际是整数类型，正确写法应为args=(4,)（末尾加逗号），否则在多参数场景下会出现异常。

修正后的代码 #

import numpy as np
import math
from scipy.integrate import quad

def firstMethod(x, p1):
    return p1 * x

def firstIntegral():
    # 修正args为元组格式
    return quad(firstMethod, 0, 1, args=(4,))[0]

def secondMethod(x, alpha):
    # 调用firstIntegral()获取积分结果，而非传递函数对象
    return alpha * math.log(firstIntegral())

def secondIntegral():
    # 修正args为元组格式
    return quad(secondMethod, 0, 10, args=(0.5,))[0]

print(secondIntegral())

额外说明 #

如果你的函数定义存在符号混淆（比如f(x)实际是依赖第二个积分变量x的函数，而非常数），需要调整firstIntegral接收参数x，例如：

def firstIntegral(x):
    # 假设p为固定值，第一个积分的被积函数是p*t，积分变量t从0到x
    return quad(lambda t: 4 * t, 0, x)[0]

def secondMethod(x, alpha):
    return alpha * math.log(firstIntegral(x))

这种情况下，f(x)会随第二个积分的变量x动态变化，最终积分结果也会不同。

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