阅读量:0
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] 上的数值积分。