numpy.trapz
是 NumPy 库中的一个函数,用于计算梯形法则下的数值积分
y
:要进行数值积分的输入数组。通常,这是一个表示函数在一系列点上的值的数组。x
(可选):一个与y
同样长度的数组,表示对应的 x 坐标。如果没有提供x
,则默认为y
的索引,即range(len(y))
。dx
(可选):一个标量,表示两个相邻点之间的间距。如果提供了x
,则不能使用dx
,因为间距将由x
数组确定。axis
(可选):一个整数,指定要沿其执行梯形法则的轴。默认为 -1,表示最后一个轴。
以下是一个使用 numpy.trapz
的示例:
import numpy as np # 定义一个函数及其在一系列点上的值 def f(x): return x * np.exp(-x) x = np.linspace(0, 5, 100) y = f(x) # 使用梯形法则计算数值积分 integral = np.trapz(y, x) print("Numerical integral:", integral)
在这个示例中,我们定义了一个函数 f(x)
,并在 0 到 5 的区间上生成了 100 个点的值。然后,我们使用 numpy.trapz
计算这些点上的数值积分。