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

《地球物理计算方法》

第2章数值积分——Romberg加速算法*复化积分可以提高精度，h-0时，I(h)=I.1、h-0,,计算量太大。2、逐步二分，如果前后偏差满足精度范围就可以了。同时继承前面的积分。1、梯形法的递推（变步长）建立T2n与Tn之间的递推关系：得到递推关系：变步长梯形法的流程图计算流程图例用变步长梯形求积法计算定积分，．解:由得由于,故有二分一次*新分点上的函数值二分二次2、龙贝格公式办法：可以采用进行外推算法，加工梯形值提高精度。梯形法的算法简单，但精度低，收敛慢。根据梯形法的误差公式，积分值Tn的误差与h2成正比，因此步长减半后误差将减至1/4，即有1、梯形法的加速移项整理，知这是T2n的事后估计法，可以用这个误差来补偿数值结果；这个值正好是辛普生公式的结果。也可写为2、辛普生法的再加速根据Simpson法的误差公式，积分值Sn的误差与h4成正比，因此步长减半后误差将减至1/16，即有移项整理，知这个值正好是科特斯公式的结果。也可写为*3、Romberg公式重复同样的手续，依据柯特斯法的误差公式可进一步导出下列龙贝格公式：这个公式也可写为*一般递推公式：