大连理工大学 计算机科学计算 第五章3.pdfVIP

第 5 章 插值函数的应用 5.3 外推加速原理与Romberg算法 5.3 外推加速原理与Romberg算法 1 1 5.3 外推加速原理与Romberg算法 由复合求积公式的余项表达式看到，数值求积的精度与 步长h有关。 如 4 b −a b −a h 2 ⎛ ⎞ (4) I −T − h f (η ) , I −Sn − ⎜ ⎟ f (η ) n 12 180 2 ⎝ ⎠ 因此，在实际计算中应选择比较小的步长，使其能计算 出好的数值结果。本节介绍便于在计算机上实现的递推算法。 5.3.1 逐次分半递推算法 设将区间[a,b]分为n等分，共有n+1个分点，如果将求积区 间再二分一次，则分点增至2n+1个，我们将二分前后两个积分 值联系起来加以考察。 每个子区间[xk, xk+ 1]经过二分只增加了一个分点： a L xk xk +1 xk +1 L b 2 2 2 则在区间[xk ,xk+ 1]上利用复化梯形公式可导出： h ⎡ ( ) +2 + ( )⎤ f x f x 1 f x k ( ) k +1 4 ⎢⎣ k +2 ⎥⎦ b a − 这里h 代表二分前的步长。 n 再将每个子区间上的积分值相加得 h n−1 h n−1 f x f x + f x , T [ ( ) + ( )] 1 ∑ k k +1 ∑ ( ) 2n k + 4 k 0 2 k 0 2 进一步有 T n 1 h n −1 h n−1 ⎛ ⎞ [f (x ) +f (x ) ] + f x 1 , T