ch23龙贝格.pptVIP

内江师范学院数学与信息科学学院 吴开腾 制作 第三节 龙贝格算法（Romberg） 复合求积方法对提高精度是一种行之有效的方法，但是在使用求积公式之前必须先给出步长。步长取得太大精度难以保证，步长太小则会导致计算量的增加，而事先给出一个合适的步长往往是困难的。为了解决这个问题，较好的方法就是采用变步长的计算方法，即在步长逐次减半的过程中，反复利用复合求积公式进行计算，直到所求得的积分值满足精度要求为止。 一、梯形法的递推公式 目的：探索梯形方法的计算规律 1、前提：设将求积区间 分为 等分，则共有 等分点 ， 表示复合梯形求积的积分值。 2、考察子区间 ，其中点为 ， 于是有 于是有递推公式为 步长等分前 步长等分后 有 3、框图 准备初值 求变步长后的值 精度判定 修改步长 4、例题分析 方法：步长逐次分半，精度逐步改进。 问题：计算积分 取n=4用复化辛普生公式，可以求得 取n=8用梯形公式，可以求得 有2位有效数字 有6位有效数字 准确值 5、讨论 梯形方法简单，但精度较低，收敛速度缓慢。收敛速度为2阶。 二、龙贝格公式 逐次分半算出新近似值 ，然后，判定误差限是否满足，否则继续逐次分半，直到满足要求为止。 由于梯形方法的收敛速度为2阶，于是有： 即： 按照事后误差估计方法，积分值 的误差大致为 ，采用事后误差估计方法，用这个误差值作为的一种补偿，可以期望： 这就是改进之后更好的结果。这个公式的实质就是辛浦生公式的积分值 。 0.9460596 5 09.9460831 10 0.9459850 4 0.9460830 9 0.9456909 3 0.9460827 8 0.9445135 2 0.9460815 7 0.9397933 1 0.9460769 6 0.9207355 0 Tn k Tn k 内江师范学院数学与信息科学学院 吴开腾 制作