SciPy - integration.quad() 方法
SciPy integrate.quad() 方法用于执行定积分任务。它通常被称为四边形(两点的集合)。
语法
以下是 SciPy integrate.quad() 方法的语法 −
quad(custom_func, upper_lim, lower_lim)
参数
此方法接受以下参数 −
- custom_func:此参数设置用户定义的积分计算函数。
- upper_lim:此参数用于设置上限。
- lower_lim:此参数用于设置下限。
返回值
此方法返回数值整数或绝对误差作为结果。
示例 1
以下是说明 SciPy integrate.quad() 方法用法的基本示例。
from scipy import integration # 定义函数 def fun(x): return x**2 # 执行集成 res, err = integration.quad(fun, 0, 1) print("结果为 ", res) print("估计误差为 ", err)
输出
上述代码产生以下输出 −
结果为 0.33333333333333337 估计误差为3.700743415417189e-15
示例 2
此程序定义函数 f(x) = sin(x) 在 0 到 pi 区间内的定积分。
import scipy.integrate as sp import numpy as np # 定义函数 def fun(x): return np.sin(x) # 执行积分 res, err = sp.quad(fun, 0, np.pi) print("结果为 ", res) print("估计误差为 ", err)
输出
上述代码产生以下输出 −
结果为 2.0 估计误差为 2.220446049250313e-14
示例 3
在这里,我们在用户定义函数中说明指数函数(e-x2),并使用 quad() 执行结果积分。这里,范围间隔设置在 -inf 和 inf 之间。
import scipy.integrate as sp import numpy as np # 定义函数 def fun(x): return np.exp(-x**2) # 执行积分 res, err = sp.quad(fun, -np.inf, np.inf) print("结果为 ", res) print("估计误差为 ", err)
输出
上述代码产生以下输出 −
结果为 1.7724538509055159 估计误差为1.4202636780944923e-08