trapz 函数在 NumPy 中用于计算梯形规则下的数值积分
numpy.trapz(y, x=None, dx=1.0, axis=-1) 参数说明:
y:输入的数组,表示要积分的函数值。x(可选):输入的数组,表示函数值对应的自变量。如果未提供,则默认为等差数列,起始值为 0,步长为 dx。dx(可选):标量,表示等差数列的步长。只有在未提供 x 时才会使用。默认值为 1.0。axis(可选):整数,表示要执行梯形规则的轴。默认值为 -1,表示最后一个轴。返回值:
示例:
import numpy as np # 定义一个函数 f(x) = x^2 def f(x): return x**2 # 生成等差数列的自变量 x = np.linspace(0, 5, 100) # 计算 f(x) 的值 y = f(x) # 使用 trapz 函数计算数值积分 integral = np.trapz(y, x) print("The integral of f(x) from 0 to 5 is:", integral) 输出:
The integral of f(x) from 0 to 5 is: 41.666666666666664 这个示例中,我们使用 trapz 函数计算了函数 f(x) = x^2 在区间 [0, 5] 上的数值积分。