第7讲(2014-11-12)

《第7讲(2014-11-12)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第4章数值积分与数值微分上次课要点：§1数值积分的基本概念与基本方法一、求积公式的一般形式其中仅依赖于节点而与无关。更一般的形式：由节点处的函数值、一阶导数值（甚至高阶导数值）的线性组合。二、有关求积公式的几个概念1．代数精度52定义1如果某个求积公式对于次数不超过的多项式都能精确成立，但对于次多项式不精确成立，则称该求积公式具有次代数精度。2．求积公式的收敛性定义2在求积公式中，若其中，则称求积公式是收敛的。523．求积公式的稳定性一个求积公式是数值稳定的，是指当输入数据的误差很小时，输出数据的误差也很小。三、插值型求积公式插值型求积公式:式

2、中。52公式余项:定理1具有个插值节点的插值型求积公式其中，为Lagrange插值基函数，至少具有次代数精度。定理2若插值型求积公式52中系数，则此求积公式是数值稳定的。§2牛顿——科特斯公式：等分，步长：，节点：牛顿——科特斯公式:其中称为科特斯系数。52可以证明，当为偶数时，牛顿——科特斯公式具有次代数精度。（1）梯形公式（2）辛普森公式（3）科特斯公式52这里。当时，科特斯系数有正有负，数值稳定性没有保证。三、几种低阶求积公式的余项1、梯形公式的余项2、辛普森公式的余项523、科特斯公式的余项。§3复化求积公式一、复化梯形公式余项：。二、

3、复化辛普森公式52余项：三、埃尔米特积分法余项：。本次课继续。。。。。。§4龙贝格求积公式(RombergIntegration)52实践中，固定复化求积公式中的步长会导致：（1）步长取得太大，精度难以保证；（2）步长太小，则计算量太大，有时不必要。一种可行的解决思路：采用变步长方法。即在步长逐步分半的过程中，反复使用复化求积公式，直到所得的积分值满足精度要求为止。一、梯形法的递推化将区间进行等分，设节点为，，则梯形法求积公式为52再将每一个子区间二等分，记相应的中点为，则此时梯形法求积公式为于是有递推公式。52注意：上述递推公式中不变。二、龙

4、贝格算法的思路上述算法过于依赖梯形公式，效率较低。问题：能否利用和的某种组合，得到精度更高的求积公式？理查德外推法：用精度较低的近似公式，组合成精度较高的近似公式，在数值计算的许多问题中都有应用。理查德在1927年发现的外推法，当时并没有用来做数值积分问题。而数学家龙贝格在1955年将它应用到数值积分上，取得很好的成效。52下面考察梯形法、辛普森法和科特斯法，当步长减半时误差的变化情况。可以看出，当减半时，梯形法：误差大致减至原有误差的；辛普森法：误差大致减至原有误差的；52科特斯法：误差大致减至原有误差的。由上述讨论，可得移项整理可得说明当很

5、小时，也很小。因此，实践中常给定计算精度，用作为停止计算的准则。进一步，由于52可以期望，用计算积分值可能得到更好的结果，而不需要增加计算量。定理4梯形法二分后的两个积分值的线性组合等于辛普森法积分值，即.52用同样方法，依据辛普森公式可导出科特斯公式依据科特斯公式可导出龙贝格公式在变步长过程中，应用上述公式，就能将粗糙的梯形公式逐步加工成精度较高的辛普森公式、科特斯公式和龙贝格公式。三、理查德外推法与龙贝格算法将区间等分后，设子区间的长度为52，用梯形法计算的积分值记为。定理5设，则有（1）其中，系数与无关。利用公式（1），可得（2）记，则（

6、消去了项）（3）其中也与无关。52由（3）得（4）记，则又可消去项，得如此进行下去，每加速一次，误差的量级便提高2阶。一般地，记，则有递推公式52上述方法称为理查德外推加速法。以表示二分次后的梯形值，表示的次加速值，则有龙贝格算法：步1初始化：计算，置（记录区间的二分次数），；步2（二分）计算积分值：52；步3（加速）求加速值：计算，；步4精度检验：对指定的精度，若，则终止计算，并取作为所求的结果；否则置，，转步2。计算次序：第1次循环52二分：加速：第2次循环二分：加速：第3次循环二分：加速：第4次循环二分：加速：……§5自适应积分方法一般而

7、言，用复合求积方法求积时：521）被积函数变化剧烈——将区间分得细一些；2）被积函数变化平缓——不需要分得很细。问题：如果被积函数在积分区间上，有的部分变化剧烈，有的部分变化平缓，则将整个积分区间一起细分是非理性的。解决方法：根据各个子区间满足精度的情况，采用不同的步长——这就是自适应积分方法。52下面以复合辛普森公式为例，介绍自适应积分方法。设给定精度要求，计算积分的近似值。先取步长，记应用辛普森公式，有52（1）将区间二等分，步长，在每一个小区间上用辛普森公式，得，52记则有即（2）如果在区间上变化不大，可假设，此时，由（1）、（2）两式可

8、得52由（2）式，得比较上述两式，得这里。如果有52（3）则可期望得到此时可取作为的近似。考虑到前面假设不一定合理，（3）式改为（4）成立时，认为也成