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

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

用法:

scipy.integrate.quadrature(func, a, b, args=(), tol=1.49e-08, rtol=1.49e-08, maxiter=50, vec_func=True, miniter=1)#

使用fixed-tolerance 高斯求積計算定積分。

使用具有絕對容差 tol 的高斯正交積分 func 從 a 到 b。

參數 ：：

func： 函數

要集成的 Python 函數或方法。

a： 浮點數

積分的下限。

b： 浮點數

積分上限。

args： 元組，可選

傳遞給函數的額外參數。

tol, rtol： 浮點數，可選

當最後兩次迭代之間的誤差小於 tol 或相對變化小於 rtol 時，迭代停止。

maxiter： 整數，可選

高斯求積的最大階數。

vec_func： 布爾型，可選

如果 func 將數組作為參數處理(是 “vector” 函數)，則為 True 或 False。默認為真。

miniter： 整數，可選

高斯求積的最小階。

返回 ：：

val： 浮點數

積分的高斯正交近似(在容差範圍內)。

err： 浮點數

積分的最後兩個估計之間的差異。

例子：

>>> from scipy import integrate
>>> import numpy as np
>>> f = lambda x: x**8
>>> integrate.quadrature(f, 0.0, 1.0)
(0.11111111111111106, 4.163336342344337e-17)
>>> print(1/9.0)  # analytical result
0.1111111111111111

>>> integrate.quadrature(np.cos, 0.0, np.pi/2)
(0.9999999999999536, 3.9611425250996035e-11)
>>> np.sin(np.pi/2)-np.sin(0)  # analytical result
1.0