龙贝格（Romberg）求积法.pptxVIP

; 龙贝格（Romberg）求积法 复化求积方法对于提高计算精度是行之有效的方法，但复化公式的一个主要缺点在于要先估计出步长。若步长太大，则难以保证计算精度，若步长太小，则计算量太大，并且积累误差也会增大。在实际计算中通常采用变步长的方法，即把步长逐次分半，直至达到某种精度为止。 1.1变步长的梯形公式 变步长复化求积法的基本思想是在求积过程中，通过对计算结果精度的不断估计，逐步改变步长（逐次分半），直至满足精度要求为止。即按照给定的精度实现步长的自动选取。 ; 设将积分区间[a,b]n等分，即分成n个子区间，一共有n+1个节点，即x=a+kh, k=0,1,…，n，步长 。对于某个子区间 ,利用梯形公式计算积分近似值有 ;; 当把积分区间分成n等份，用复化梯形 公式计算积分I的近似值 时，截断误差为 ; 可见,当步长二分后误差将减至 ,将 上式移项整理，可得验后误差估计式 ; 1.2 变步长的梯形求积算法实现 （1）变步长的梯形求积法的计算步骤 ① 变步长梯形求积法。它是以梯形求积公式为基础，逐步减少步长，按如下递推公式求二分后的梯形值;（2）变步长梯形公式的流程图 ;(3)程序实现(见附录A A-13 变步长梯形求积法) ;;1.3 龙贝格求积公式 变步长梯形求积法算法简单，但精度较差，收敛速度较慢，但可以利用梯形法算法简单的优点，形成一个新算法，这就是龙贝格求积公式。龙贝格公式又称逐次分半加速法。 根据积分区间分成n等份和2n等份时的误差估计式(6.8)可得 ;;再考察辛卜生法。其截断误差与 成正比，因此，如果将步长折半，则误差减至 ，即有 ;用同样的方法，根据柯特斯公式的误差公式，可进一步导出龙贝格公式 ;1.4 龙贝格求积法算法实现 （1） 龙贝格求积法计算步骤 用梯形公式计算积分近似值 按变步长梯形公式计算积分近似值 将区间逐次分半,令区间长度 ;④ 精度控制；直到相邻两次积分值 ;例6.16 用龙贝格算法计算定积分 要求相邻两次龙贝格值的偏差不超过 解:由题意 ;;;