《复化数值积分.doc》由会员分享，可在线阅读，更多相关《复化数值积分.doc（13页珍藏版）》请在金锄头文库上搜索。
1、复化数值积分由插值的龙格现象可知，高阶牛顿-柯特斯积分不能保证等距数值积分系列的收敛性，同时可证（略）高阶牛顿-柯特斯积分的计算是不稳定的。因此，实际计算中常用低阶复化梯形等积分公式。7.3.1 复化梯形积分把积分区间分割成若干小区间，在每个小区间上用梯形积分公式，再将这些小区间上的数值积分累加起来，称为复化梯形公式。复化梯形公式用若干个小梯形面积逼近积分比用一个大梯形公式效果显然更好，如图7.4所示。这种作法使我们想起定积分定义，即它为被积函数无限分割的代数和。这也正是计算定积分最朴素的算法。图7.4 复化梯形公式积分视图复化梯形积分计算公式对作等距分割，有，于是在上，则有记等分的复化梯形公
2、式为，有（7.11）复化梯形公式截断误差由，根据均值定理，当时，存在，有，于是 （7.12） 由此看到复化梯形公式的截断误差按照或者的速度下降，事实上，可以证明，只要在上有界并黎曼可积，当分点无限增多时，复化梯形公式收敛到积分。记，则有对于任给的误差控制小量，有 或 就有，式中表示取其最大整数。7.3.2 复化辛普森积分把积分区间分成偶数等分，记，其中是节点总数，是积分子区间的总数。记，在每个区间上用辛普森数值积分公式计算，则得到复化辛普森公式，记为。复化辛普森积分计算公式而，称（7.13）为复化辛普森积分公式，它是在上采用辛普森积分公式叠加而得。下面用图7.5显示复化辛普森积分计算公式中节点
3、与系数的关系，取，在每个积分区间上提出因子后，三个节点的系数分别是1，4，1；将4个积分区间的系数按节点的位置累加，可以清楚地看到，首尾节点的系数是1，奇数点的系数是4，偶数点的系数是2。141141141141142424241图7.5 复化辛普森积分系数复化辛普森公式的截断误差设，在上的误差为因此，即 （7.14）与复化梯形公式类似，误差的截断误差按照或者的速度下降。可以证明，只要在上有界并黎曼可积，当分点无限增多时，复化辛普森公式收敛到积分。记，则有对任给的误差控制小量，只要 或 就有。例7.5 求，计算中要求有5位有效数字。用复化梯形和复化辛普森求积公式的分点应取多少？解： 由复化梯形
4、误差公式得到：计算出，复化梯形公式至少要在.00等分n = 68。由复化辛普森误差公式，有在复化辛普森公式中取 或。 7.3.3 复化积分的自动控制误差算法复化积分的误差公式表明，截断误差随分点的增大而减小，对于给定的误差量，用估计函数导数的界的方法可计算出。用误差公式计算满足精度的分点数，像是在做一道计算导数上界的微积分习题（如例7.5所示）。但是在实际运算中，一般难以估计出函数的各阶导数界，也就无法确定分点数。在计算中常用误差的事后估计方法，即用估计误差。T2n (f )的计算公式对定积分，取分点，计算得取分点，计算得这里，。可以看到，的值是与新增分点的组合。取分点，计算得这里，。同理，计
5、算时只要在的基础上计算新增分点，的值再做组合，如图7.6所示。图7.6 与一般地，每次总对前一次的小区间分半，分点加密一倍，并可充分利用老分点上的函数值，每次只需计算新增分点的和。对上等分， ，则有记上的中点为则 （7.15）其中。或 其中。类似地，可得积分节点为，的辛普森求积公式的关系式：（7.16）其中：由误差公式：由于，分别为及个点上的均值，可视，于是上式表明的误差大约是误差的4倍。或（7.17）由此得到启发，对任给的误差控制量，要，只需即可，而用作为控制手段简单直接，序列在计算机上也不难实现。复化积分的算法描述从数值积分的误差公式可以看到，截断误差随分点的增长而减少，控制计算的精度也就
6、是确定分点数。在计算中不用数值积分的误差公式确定分点数的理论模式，而用作为控制，通过增加分点自动满足精度的方法称为数值积分公式的自动积分法。即在计算中构造序列，直到或时停止计算，由分点数自动控制积分值的误差，并取。下面描述复化数值积分公式的自动控制误差算法，详细程序和算例请看本章7.6节。1输入：误差控制精度e = eps；初始分点值 。2计算分点的复化梯形积分T1=T2+100 /迭代计算中T1和T2分别表示和3while | T1T2|eT1 = T2H=Hn /计算新增节点的值T2= (T1+H )/2H = h/2，n =2n /将区间一分为二end while4输出积分值T2。在自动
7、控制误差算法中初始分点值不宜过小，以防假收敛。7.3.4 龙贝格（Romberg）积分由前面得到的关系式（7.17），可以将 作为的修正值补充到，即（7.18）其结果是将梯形求积公式组合成辛普森求积公式，截断误差由提高到。这种手段称为外推算法。外推算法在不增加计算量的前题下提高了误差的精度，是计算方法中一种常用手法。不妨对再做一次组合。由得到 （7.19）复化辛普森公式组成复化柯特斯公式，其截断误差是。同理对柯特斯公式进行组合：得到具有7次代数精度和截断误差是 的龙贝格公式：还可以继续对做上去。为了便于在计算机上实现龙贝格算法，将统一用表示，列标 分别表示梯形、辛普森、柯特斯积分，行标表示分点数或步长j。龙贝格计算公式：对每一个从2做到，一直做到小于给定控制精度停止计算。龙贝格算法龙贝格算法按表7.2元素的行序进行运算， 在计算中每个元素只用到上一行和本行的元素。对上面的算法进一步优化，对每k行可将计算定义在两行元素之间，令；在每计算一行元素后，要将。表7.2 龙贝格算法计算元素顺序表 1输入区间端点 ，精度控制值，循环次数，定义函数，取；2；3for =2 to 4输出。
