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] 上的数值积分。