数值计算方法课件CH4 数值积分43 Romberg算法.ppt

华长生制作 * 第四章 数值积分 4.3 Romberg 算法 4.3 Romberg算法 由前几节的内容可知,梯形公式、Simpson 公式、Cotes 公式的代数精度分别为： 1 次、3 次和 5 次 复合梯形、复合Simpson、复合Cotes 公式的收敛阶分别为: 2 阶、4 阶和 6 阶 无论从代数精度还是收敛速度, 复合梯形公式都是较差的 问：有没有办法改善梯形公式呢? 龙贝格(Romberg)算法：在复合梯形公式误差估计的基础上，应用线形外推的方法构造出的一种加速算法.这种方法利用误差补偿，将复合梯形公式加工成具有高收敛阶，高精确度的积分近似值. 龙贝格(Romberg)算法又称作逐次分半加速法. -------- (1) -------- (3) 各节点为： 复合梯形公式： -------- (2) 将定积分 的积分区间 分割为n 等份, 一、复合梯形公式的递推化 新步长 新增节点 步长折半后, 递推形式的复合梯形公式只需计算新增节点处的函数值, 节省一半的计算量! 若 , n为区间等分数,经若干次步长折半所得,必为偶数！ 记 ， 则 , 因此(1)(2)(3)式可化为如下的 递推梯形公式 注: 递推梯形公式加上一个控制精度, 即可成为自动选取 步长的复合梯形公式. 递推梯形公式的前后两步之差: 可视为当前步骤的计算误差. 步长折半后的新步长 步长折半后的新增节点 (4)------- 例1 用递推的梯形公式计算积分 要求误差不超过 解 0 1 2 3 0 1 1.000 000 00 0.841 470 98 0.920 735 49 0.939 793 28 0.958 851 08 0.5 0.25 0.75 0.125 0.375 0.625 0.875 0.989 615 84 0.908 851 68 0.997 397 88 0.976 726 74 0.936 155 64 0.945 690 87 0.944 513 52 0.945 690 86 6.352?10-3 1.573?10-3 3.924?10-4 取 对于上例,若要求 ,则需要对分10次,涉及到1000多个点处的函数值,得结果0.9460831. 问: 是否存在其他的方法, 能够加速积分的收敛？ 下面利用复合梯形余项公式、复合 Simpson 余项公式和复合 Cotes 余项公式对上述递推公式进行加速——— Romberg算法 复合梯形公式的余项: 二、外推加速公式 令 可以证明： 复合 Simpson 公式的余项: 可得 令 即 自然有 即 同样, 得 令 自然有 由复合 Cotes 公式的余项: 外推加速公式 Romberg 算法 将上述结果综合后得: 递推梯形公式 Romberg 算法求解步骤 k 0123 m 方向 k 方 向 外推加速公式的一般形式: --------(9) Romberg算法控制计算停止的条件: 对于给定的误差限ε，当 时，取 ，停止计算。 例2 用Romberg算法计算积分 要求误差不超过 * * * *