第7讲(2014-11-12)概念.docVIP

PAGE  PAGE 53 数值积分与数值微分 上次课要点： §1 数值积分的基本概念与基本方法 一、求积公式的一般形式 其中仅依赖于节点而与无关。 更一般的形式：由节点处的函数值、一阶导数值（甚至高阶导数值）的线性组合。 二、有关求积公式的几个概念 1．代数精度 定义1 如果某个求积公式对于次数不超过的多项式都能精确成立，但对于次多项式不精确成立，则称该求积公式具有次代数精度。 2．求积公式的收敛性 定义2 在求积公式 中，若 其中，则称求积公式是收敛的。 3．求积公式的稳定性 一个求积公式是数值稳定的，是指当输入数据的误差很小时，输出数据的误差也很小。 三、插值型求积公式 插值型求积公式: 式中。 公式余项: 定理1 具有个插值节点的插值型求积公式 其中，为Lagrange插值基函数，至少具有次代数精度。 定理2 若插值型求积公式 中系数，则此求积公式是数值稳定的。 §2 牛顿——科特斯公式 ：等分， 步长：， 节点： 牛顿——科特斯公式: 其中称为科特斯系数。 可以证明，当为偶数时，牛顿——科特斯公式具有次代数精度。 （1）梯形公式 （2）辛普森公式 （3）科特斯公式 这里。 当时，科特斯系数有正有负，数值稳定性没有保证。 三、几种低阶求积公式的余项 1、梯形公式的余项 2、辛普森公式的余项 3、科特斯公式的余项 。 §3 复化求积公式 一、复化梯形公式 余项： 。 二、复化辛普森公式 余项： 三、埃尔米特积分法 余项： 。 本次课继续。。。。。。 §4 龙贝格求积公式 (Romberg Integration) 实践中，固定复化求积公式中的步长会导致： （1）步长取得太大，精度难以保证； （2）步长太小，则计算量太大，有时不必要。 一种可行的解决思路：采用变步长方法。即在步长逐步分半的过程中，反复使用复化求积公式，直到所得的积分值满足精度要求为止。 一、梯形法的递推化 将区间进行等分，设节点为，，则梯形法求积公式为 再将每一个子区间二等分，记相应的中点为，则此时梯形法求积公式为 于是有递推公式 。 注意：上述递推公式中不变。 二、龙贝格算法的思路 上述算法过于依赖梯形公式，效率较低。 问题：能否利用和的某种组合，得到精度更高的求积公式？ 理查德外推法：用精度较低的近似公式，组合成精度较高的近似公式，在数值计算的许多问题中都有应用。 理查德在1927年发现的外推法，当时并没有用来做数值积分问题。而数学家龙贝格在1955年将它应用到数值积分上，取得很好的成效。 下面考察梯形法、辛普森法和科特斯法，当步长减半时误差的变化情况。 可以看出，当减半时， 梯形法：误差大致减至原有误差的； 辛普森法：误差大致减至原有误差的； 科特斯法：误差大致减至原有误差的。 由上述讨论，可得 移项整理可得 说明当很小时，也很小。因此，实践中常给定计算精度，用作为停止计算的准则。 进一步，由于 可以期望，用 计算积分值可能得到更好的结果，而不需要增加计算量。 定理4 梯形法二分后的两个积分值的线性组合 等于辛普森法积分值，即 . 用同样方法，依据辛普森公式可导出科特斯公式 依据科特斯公式可导出龙贝格公式 在变步长过程中，应用上述公式，就能将粗糙的梯形公式逐步加工成精度较高的辛普森公式、科特斯公式和龙贝格公式。 三、理查德外推法与龙贝格算法 将区间等分后，设子区间的长度为，用梯形法计算的积分值记为。 定理5 设，则有 （1） 其中，系数与无关。 利用公式（1），可得 （2） 记 ，则 （消去了项） （3） 其中也与无关。 由（3）得 （4） 记，则又可消去项，得 如此进行下去，每加速一次，误差的量级便提高2阶。 一般地，记 ， 则有递推公式 上述方法称为理查德外推加速法。 以表示二分次后的梯形值，表示的次加速值，则有 龙贝格算法： 步1 初始化：计算，置（记录区间的二分次数），； 步2（二分） 计算积分值： ； 步3（加速） 求加速值：计算 ，； 步4 精度检验：对指定的精度，若 ，则终止计算，并取作为所求的结果；否则置，，转步2。 计算次序： 第1次循环 二分： 加速： 第2次循环 二分： 加速： 第3次循环 二分： 加速： 第4次循环 二分： 加速： …… §5 自适应积分方法 一般而言，用复合求积方法求积时： 1）被积函数变化剧烈 ——将区间分得细一些； 2）被积函数变化平缓 ——不需要分得很细。 问题： 如果被积函数在积分区间上，有的部分变化剧烈，有的部分变化平缓，则将整个积分区间一起细分是非理性的。 解决方法： 根据各个子区间满足精度的情况，采用不同的步长——这就是自适应积分方法。 下面以复合辛普森公式为例，介绍自适应积分方法。 设给定精度要求，计算积分 的近似值。