Python | Scipy integration.quad() 方法

Scipy是一个面向科学计算的Python库，它提供了许多高效的数学算法和函数。其中，integration.quad()是一个在众多积分方法中的一种，用于数值积分。

语法

scipy.integrate.quad(func, a, b, args=(), **kwargs)

其中，参数说明如下：

• func：被积函数。
• a：积分下限。
• b：积分上限。
• args：被积函数的其他参数（可选）。
• **kwargs：其他可选参数。

返回值

返回一个元组 (y, abserr)，其中：

• y：数值积分的结果。
• abserr：基于达到所需精度而估计的误差的估计值。

示例

下面是一个简单的示例，计算正态分布的概率密度函数在区间 [0, 1] 上的积分值：

from scipy.integrate import quad
from scipy.stats import norm

def integrand(x):
    return norm.pdf(x)

result, error = quad(integrand, 0, 1)
print(result)

输出结果为：

0.34134474606854304

注意事项

• 被积函数必须是实数函数。
• 积分区间必须是有限区间。
• args 参数应该是一个元组，否则可能导致出错。
• quad() 会根据被积函数的特点和给定的积分区间自动选择合适的积分方法和步长大小，因此可以达到较高的精度。