Python SciPy integrate.cumulative_trapezoid用法及代碼示例

本文簡要介紹 python 語言中 scipy.integrate.cumulative_trapezoid 的用法。

用法:

scipy.integrate.cumulative_trapezoid(y, x=None, dx=1.0, axis=-1, initial=None)#

使用複合梯形規則對 y(x) 進行累積積分。

參數 ：：

y： array_like

要整合的值。

x： 數組，可選

要整合的坐標。如果無(默認)，則在 y 中的連續元素之間使用間距 dx。

dx： 浮點數，可選

y 元素之間的間距。僅在 x 為無時使用。

axis： 整數，可選

指定要累積的軸。默認值為 -1(最後一個軸)。

initial： 標量，可選

如果給定，則將此值插入到返回結果的開頭。 0 或 None 是唯一接受的值。默認值為 None，這意味著 res 沿積分軸有一個比 y 少的元素。

返回 ：：

res： ndarray

y 沿軸累積積分的結果。如果 initial 為 None，則形狀使得積分軸的值比 y 小一個。如果給定初始值，則形狀等於 y 的形狀。

例子：

>>> from scipy import integrate
>>> import numpy as np
>>> import matplotlib.pyplot as plt

>>> x = np.linspace(-2, 2, num=20)
>>> y = x
>>> y_int = integrate.cumulative_trapezoid(y, x, initial=0)
>>> plt.plot(x, y_int, 'ro', x, y[0] + 0.5 * x**2, 'b-')
>>> plt.show()