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