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

《地球物理计算方法》

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