龙贝格求积法.ppt

龙贝格求积法.ppt

《龙贝格求积法.ppt》由会员分享，可在线阅读，更多相关《龙贝格求积法.ppt（24页珍藏版）》请在冰豆网上搜索。

4.4外推原理与外推原理与Romberg求积方法求积方法4.4.1外推原理外推原理在科学与工程计算中，很多算法与步长h有关，特别是数值积分、数值微分和微分方程数值解的问题。

对于这些算法，我们可以通过外推技巧提高计算精度。

例1计算的近似值。

由函数的Taylor展开式有若记，则有由此构造新的表达式可见计算的近似值的算法F（h）的截断误差是O（h2）,而算法F1（h）的截断误差是O（h4）。

外推一次，精度就提高了。

这就是外推法的基本思想。

若重复以上过程，不断外推，即不断折半步长h,得到计算的算法序列Fk（h）。

随着k的增加，算法的截断误差阶越来越高，计算精度越来越好。

可将上述外推思想推广到一般情况。

设F（h）是计算F（0）的一种近似算式，带截断误差的表达式为其中，与h无关。

如果我们用h和h/q（q1）两种步长分别计算F（h）和F（h/q），则有削去截断误差的主项，得新的算法我们称这个计算过程为Richardson外推法。

F1（h）逼近F（0）的截断误差是只要知道F（h）的更加完整的关于h幂的展开式，而无需知道展开式中各个系数的具体数值，就能重复使用Richardson外推法，直到截断误差达到容许误差。

用归纳法可以证明下面更一般的定理。

定理4.5假设F（h）逼近F（0）的余项为其中，是与h无关的非零常数。

则由定义的序列Fk（h）有无关，q1.Richardson外推法应用非常广泛且有效，下面介绍应用于数值积分的情形。

变步长梯形求积法变步长梯形求积法算法简单，但精度较差，收敛速度较慢，算法简单，但精度较差，收敛速度较慢，但可以利用梯形法算法简单的优点，形成一个新算法，这就是龙但可以利用梯形法算法简单的优点，形成一个新算法，这就是龙贝格求积公式。

贝格求积公式。

龙贝格公式龙贝格公式又称又称逐次分半加速法逐次分半加速法。

Romberg算法算法是在积分区间逐次分半的过程中，对用复化梯是在积分区间逐次分半的过程中，对用复化梯形法形法产生的近似值产生的近似值进行进行加权平均加权平均，以获得准确度较高的近似值的，以获得准确度较高的近似值的一种方法，具有公式简练，使用方便，结果较可靠的优点。

一种方法，具有公式简练，使用方便，结果较可靠的优点。

二、二、Romberg算法算法根据梯形法的误差公式，可知积分值根据梯形法的误差公式，可知积分值的截断误差大致与的截断误差大致与成正比，因此当步长二分后，截断误差将减至原有误差的成正比，因此当步长二分后，截断误差将减至原有误差的1/4，即有即有将上式移项整理，可得将上式移项整理，可得（4.4.3）由此可见，只要二分前后的两个积分值由此可见，只要二分前后的两个积分值与与相当接近，就可以保相当接近，就可以保证证计算结果的误差很小。

这样直接用计算结果来估计误差的方法计算结果的误差很小。

这样直接用计算结果来估计误差的方法通常称作通常称作误差的事后估计法误差的事后估计法。

按式按式（4.4.3）,积分近似值与积分近似值与的误差大致等于的误差大致等于因此如果因此如果用这个误差值作为用这个误差值作为的的一种补偿一种补偿，可以期望所得到的，可以期望所得到的（4.4.4）可能是更好的结果。

可能是更好的结果。

（4.4.5）有可能比有可能比更好地接近于积分更好地接近于积分的真值的真值I即即这就是说，用梯形法二分前后两个积分值这就是说，用梯形法二分前后两个积分值和和作线性组合作线性组合与与这表明在收敛缓慢的梯形数列这表明在收敛缓慢的梯形数列的基础上，若将的基础上，若将按按（4.4.6）作线性组合就可产生收敛速度较快的作线性组合就可产生收敛速度较快的Simpson序列：

序列：

、关于关于（4.4.3.6）的证明：

由的证明：

由（4.4.1）可知：

可知：

这种新近似值这种新近似值实质上又是什么呢？

可以验证：

实质上又是什么呢？

可以验证：

即：

即：

（4.4.6）故（故（4.4.6）式成立）式成立.（由上节梯形公式由上节梯形公式）（由上节（由上节Simpson公式）公式）同理，由上节近似式同理，由上节近似式类似推导可得：

类似推导可得：

（4.4.7）即将即将Simpson序列序列按按（4.4.7）作线性组合就可产生收作线性组合就可产生收敛敛速度更快的新序列速度更快的新序列-Cotes序列序列即在即在Cotes序列序列的基础上，产生了一个称为的基础上，产生了一个称为Romberg、.序列序列的新序列的新序列（4.4.8）由近似式，类似推导可得：

我们在变步长的过程中运用公式我们在变步长的过程中运用公式（4.4.4-4.4.8）（也称它们也称它们为加速公式为加速公式），就能将粗糙的梯形值，就能将粗糙的梯形值Tn逐步加工成精度较高的逐步加工成精度较高的辛普森值辛普森值Sn、柯特斯值柯特斯值Cn和龙贝格值和龙贝格值Rn.可以证明：

可以证明：

当当f（x）满足一定条件时，满足一定条件时，Romberg序列序列比比Cotes序列序列能更快地收敛到积分能更快地收敛到积分的真值的真值I。

越越精确的近似值精确的近似值也就是将收敛速度缓慢也就是将收敛速度缓慢综上可知：

在积分区间逐次分半过程中利用公式综上可知：

在积分区间逐次分半过程中利用公式可以将粗糙的近似值可以将粗糙的近似值逐步地逐步地“加工加工”成越来成越来的梯形序列的梯形序列逐步地逐步地“加工加工”成收敛速度越来越快的新成收敛速度越来越快的新新序列新序列（4.4.3-4.4.8）这种加速的方法称为这种加速的方法称为Romberg算法。

其算法。

其“加工加工”过程如下图，过程如下图，其中圆圈中号码表示计算顺序。

其中圆圈中号码表示计算顺序。

T1,T2、T4、T6由梯由梯形公式、复化梯形公式的形公式、复化梯形公式的递推化公式求得递推化公式求得020=1T1121=2T2S1222=4T4S2C1323=8T8S4C2R1424=16T16S8C4R2525=32T32S16C8R4区间等分数区间等分数梯形序列梯形序列辛普森序列辛普森序列柯特斯序列柯特斯序列龙贝格序列龙贝格序列龙贝格求积算法也可用下表来表示：

龙贝格求积算法也可用下表来表示：

例例11:

利用利用Romberg算法计算算法计算解解:

由题意由题意计算到计算到R1.例例2用用Romberg算法算法计算计算得到的梯形值，计得到的梯形值，计算结果见算结果见表表4-5（k代表二分次数代表二分次数）。

计算值的误差不超过。

计算值的误差不超过0.510-6.表表4-5我们看到，这里利用二分我们看到，这里利用二分3次的数据次的数据（它们的精度都很差，只有二三位它们的精度都很差，只有二三位是有效数字是有效数字），通过三次加速求得，通过三次加速求得=0.9460831，这个结果的每一位，这个结果的每一位数字都是有效数字，可见加速的效果是十分显著的。

三次外推后达到数字都是有效数字，可见加速的效果是十分显著的。

三次外推后达到6位有效数字。

位有效数字。

注意注意教材中介绍的教材中介绍的Richardson外推法外推法，为便于上机计算，引用记，为便于上机计算，引用记号号来表示各近似值，其中来表示各近似值，其中k仍代表积分区间的二分次数，而仍代表积分区间的二分次数，而下标下标m则指出了近似值则指出了近似值所在序列的性质。

如所在序列的性质。

如m=1在梯形序列中，在梯形序列中，m=2在在Simpson序列中，序列中，m=3在在Cotes序列中，序列中，引入上面记号引入上面记号后，后，Romberg算法所用到的各个计算公式可统一化为：

算法所用到的各个计算公式可统一化为：

三、三、Romberg算法计算公式的简化算法计算公式的简化由此可逐行构造由此可逐行构造出一个三角形数表出一个三角形数表-称为称为T数表数表ki0123Romberg算法中止准则，一般取同列或同行相邻两数值的误差算法中止准则，一般取同列或同行相邻两数值的误差绝对绝对实际计算中常常只计算到第实际计算中常常只计算到第4列列,只使用只使用3次理查逊外推法。

次理查逊外推法。

注：

注：

值小于事先给定的精度要求。

取最后一次的数值积分值作为积值小于事先给定的精度要求。

取最后一次的数值积分值作为积分的分的近似值。

近似值。

外推次数外推次数分分半半次次数数复化梯复化梯形序列形序列Simpson序列序列Cotes序列序列Romberg序列序列T1=T8=T4=T2=?

S1=R1=S2=C1=C2=S4=1、在上面、在上面“加工加工”过程中的系数过程中的系数和和，当，当m4时，时，而另一个而另一个系数则接近于系数则接近于1，也就是，也就是新新公式与原公式差别不大，公式与原公式差别不大，但工作量却大增。

但工作量却大增。

因此，在实际计算中常规定因此，在实际计算中常规定m3，即，即计算到出现计算到出现Romberg序列为止。

序列为止。

2、可用二维数组来存放并参加运算，也可用一维数组。

可用二维数组来存放并参加运算，也可用一维数组。

四、几点说明：

四、几点说明：

3、对于积分限为无穷的积分对于积分限为无穷的积分，可利用变量代换化成有限区，可利用变量代换化成有限区间的积分然后再进行计算。

例如：

间的积分然后再进行计算。

例如：

4、若被积函数有奇异点（间断点）存在于积分区间内，则、若被积函数有奇异点（间断点）存在于积分区间内，则可将积分可将积分区间分成小部分，使间断点在子区间的端点处。

区间分成小部分，使间断点在子区间的端点处。

也可用变量代换法处理。

也可用变量代换法处理。

令令则则代入得：

代入得：

例例3用龙贝格方法计算椭圆用龙贝格方法计算椭圆x2/4+y2l的周长，使结果的周长，使结果具有五位有效数字具有五位有效数字分析分析为便于计算，先将椭圆方程采用参数形式表示为便于计算，先将椭圆方程采用参数形式表示,再根再根据弧长公式将椭圆周长用积分形式表示由于计算结果要求具据弧长公式将椭圆周长用积分形式表示由于计算结果要求具有五位有效数字，因此需要估计所求积分值有几位整数，从而有五位有效数字，因此需要估计所求积分值有几位整数，从而确定所求积分值的绝对误差限最后再应用龙贝格方法计算积确定所求积分值的绝对误差限最后再应用龙贝格方法计算积分分解解令令x2cos，ysin，则椭圆的周长为则椭圆的周长为下表给出了用龙贝格方法计算积分下表给出了用龙贝格方法计算积分I=1+1+3sin2dx的过程的过程./20012.356194122.4199212.441163242.4221032.4228302.421608382.4221122.4221152.4220672.4220744162.4221122.4221122.4221122.4221130.0000395322.4221122.4221122.4221122.4221120.0000010.12510-4故积分故积分I2.422112,椭圆周长的近似值为椭圆周长的近似值为l=4I9.6884。