求积公式
一般形式
插值型求积公式
Newton-Cotes公式
形式
Cotes 系数表
n \ k | 0 | 1 | 2 | 3 | 4 |
---|---|---|---|---|---|
1 | 12 | 12 | |||
2 | 16 | 46 | 16 | ||
3 | 18 | 38 | 38 | 18 | |
4 | 790 | 1645 | 215 | 1645 | 790 |
梯形求积公式
当
n=1
时,Newton-Cotes公式就成了梯形求积公式,如下:
说明: M2 代表 f(x) 的二阶导数的最大值。
Simpson 求积公式
当
n=2
时,Newton-Cotes公式就成了
Simpson
求积公式,如下:
说明: M4 代表 f(x) 的四阶导数的最大值。
Cotes 求积公式
当
n=4
时,Newton-Cotes公式就成了
Cotes
求积公式,如下:
说明: M6 代表 f(x) 的六阶导数的最大值。
复化求积公式
把积分区间划分为
n
个小区间,则每小段区间长度
复化梯形求积公式
对梯形求积公式进行复化,可以得到:
其中, xk=a+kh , M2 代表 f(x) 的二阶导数的最大值。
复化 Simpson 求积公式
对
Simpson
求积公式进行复化,可以得到:
其中, xk=a+kh , xk−12=a+(k−12)h , M4 代表 f(x) 的四阶导数的最大值。
复化 Cotes 求积公式
对
Cotes
求积公式进行复化,可以得到:
其中, xk=a+kh , xk−34=a+(k−34)h , xk−12=a+(k−12)h , xk−14=a+(k−14)h , M6 代表 f(x) 的六阶导数的最大值。
Gauss 型求积公式
定义6.1 代数精度
若求积公式
则称此公式 具有 m 次代数精度.
可见, 若公式具有
定理6.1 Gauss型求积公式的一般理论
区间
[a,b]
上权函数为
ρ(x)
的具有
n
个节点的数值积分公式代数精度不超过
证明
若记 x1,x2,…xn 为求积节点,则积分公式为
定义6.2 Gauss 型求积公式
使求积公式具有 2n−1 次代数精度的节点 x1,x2,…xn 称为 Gauss 点,此时的插值型求积公式称为 Gauss 型求积公式.
定理6.2 Gauss 点定理
对区间
[a,b]
上权函数为
ρ(x)
的积分
区间 [a,b] 上权函数为 ρ(x) 的正交多项式 pn(x) 的 n 个零点
构造方法
1.求
[a,b]
上权函数为
ρ(x)
的
n
次正交多项式
2.求
\pn(x)=0
的
n
个零点,作为
3.计算系数
Ai=∫baρ(x)li(x)dx,i=1,2,...,n.
;
4.给出结果
5.计算余项步骤
1>构造 2n−1 次 Hermite 插值多项式 H2n−1(x) 使 H2n−1(xi)=f(xi) , H′2n−1(xi)=f′(xi) .则 R2n−1(x)=f(2n)(ξx)(2n)!ω2n(x) .
因为 Gauss 型求积公式具有 2n−1 次代数精度,所以有
2>