《地球物理计算方法》第2章 数值积分——Romberg加速算法.pptxVIP

《地球物理计算方法》

第2章数值积分——Romberg加速算法

2复化积分可以提高精度，h-0时，I(h)=I.1、h-0,,计算量太大。2、逐步二分，如果前后偏差满足精度范围就可以了。同时继承前面的积分。1、梯形法的递推（变步长）

建立T2n与Tn之间的递推关系：得到递推关系：

变步长梯形法的流程图计算流程图

例用变步长梯形求积法计算定积分，．解:由得由于,故有二分一次

6新分点上的函数值二分二次

2、龙贝格公式办法：可以采用进行外推算法，加工梯形值提高精度。梯形法的算法简单，但精度低，收敛慢。根据梯形法的误差公式，积分值Tn的误差与h2成正比，因此步长减半后误差将减至1/4，即有1、梯形法的加速

移项整理，知这是T2n的事后估计法，可以用这个误差来补偿数值结果；这个值正好是辛普生公式的结果。也可写为

2、辛普生法的再加速根据Simpson法的误差公式，积分值Sn的误差与h4成正比，因此步长减半后误差将减至1/16，即有移项整理，知这个值正好是科特斯公式的结果。也可写为

103、Romberg公式重复同样的手续，依据柯特斯法的误差公式可进一步导出下列龙贝格公式：这个公式也可写为

11一般递推公式：

用梯形公式计算积分近似值按变步长梯形公式计算积分