计算方法-数值积分b.ppt

5.3 复合求积 /* Composite Quadrature */ 高次插值有Runge 现象，故采用分段低次插值 ? 分段低次合成的 Newton-Cotes 复合求积公式。 ? 复合梯形公式 在每个 上用梯形公式： = Tm 将 [a, b]区间m等分，步长h=(b-a)/m，分点xk=a+kh, k=0,1,…, m。用低阶牛顿—柯特斯公式求子区间[xk, xk+1]上的积分值，再累加得到积分的近似值。 思路 复合公式的余项？ 5.3 Composite Quadrature ? 复化 Simpson 公式 4 4 4 4 4 = Sm ? 复化 Cotes 公式 Cm 若f(x)在积分区间[a, b]上分别具有二阶、四阶和六阶连续导数，则复化积分公式的余项分别是 定理 其中，ξ∈[a,b]，且当h充分小时，又有 5.3 Composite Quadrature 证 这里考察复化梯形公式的余项公式，其余由同学们自己完成。 因此 中值定理　　　 　　由余项公式可以看出，只要相关的各阶导数在积分区间上连续，则当m→∞(即h→0)时，Tm，Sm，Cm都收敛于积分真值，而且收敛速度一个比一个快。 另有 由于 在[a, b]上连续，故每个小区间上的积分使用梯形公式时，有误差为 5.3