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
計算這些點上的數值積分。