第三章_数值积分和数值微分.pptVIP

(3.5) 考察 与n等份辛卜生公式 之间的关系。将 复化梯形公式 梯形变步长公式 代入(3.5) 表达式得 故 这就是说，用梯形法二分前后两个积分值 和 作线性组合，结果却得到复化辛卜生公式计算得到 的积分值 。 再考察辛卜生法。其截断误差与 成正比，因此，如果将步长折半，则误差减至 ，即有 由此可得 可以验证,上式右端的值其实等于Cn，就是说，用辛卜生公式二等份前后的两个积分值Sn和S2n 作线性组合后，可得到柯特斯公式求得的积分值Cn，即有 (3.6) 用同样的方法，根据柯特斯公式的误差公式，可进一步导出龙贝格公式 (3.7) 在变步长的过程中运用(3.6)和(3.7)，就能将粗糙的梯形值Tn逐步加工成精度较高的辛卜生值Sn、柯特斯值Cn和龙贝格值Rn或者说，将收敛缓慢的梯形值序列Tn加工成收敛迅速的龙贝格值序列Rn，这种加速方法称为龙贝格算法（龙贝格公式），见下表所示。 T1 T2 S1 T4 S2 C1 T8 S4 C2 R1 T16 S8 C4 R2 T32 S16 C8 R4. ................................................ 上面是Romberg的计算表若 则计算停止 3.3.4 龙贝格求积法算法实现 （1） 龙贝格求积法计算步骤 用梯形公式计算积分近似值 按变步长梯形公式计算积分近似值 将区间逐次分半,令区间长度 计算 ③ 按加速公式求加速值 梯形加速公式： 辛卜生加速公式： 龙贝格求积公式： ④ 精度控制；直到相邻两次积分值 （其中ε为允许的误差限）则终止计算并取Rn 作为积分 的近似值，否则将区间再对分，重复 ②，③，④ 的计算，直到满足精度要求为止。 例 用龙贝格算法计算定积分 要求相邻两次龙贝格值的偏差不超过 解:由题意 由于 ，于是有 3.7 高斯（Gauss）型求积公式* 3.7.1 高斯积分问题的提出 在前面建立牛顿-柯特斯公式时，为了简化计算，对插值公式中的节点限定为等分的节点，然后再定求积系数，这种方法虽然简便，但求积公式的精度受到限制。我们已经知道，过n+1个节点的插值形求积公式至少具有n次代数精度，我们不仅要问，是否存在具有最高代数精度的求积公式呢？若有，最高代数精度能达到多少呢？让我们先看一个例子： 在构造形如 的两点公式时,如果限定求积节点, 那么所得插值求积公式 （3.8） 的代数精度仅为1。但是,如果对式(3.8)中的系数 和 节点都不加限制，那么就可适当选取 和 ,使所得公式的代数精度 。事实上，若要使求积公式(3.8)对函数 都准确成立，只要 和 满足方程组 解之得 代入(3.8)即得 (3.9) 可以验证，所得公式（3.9）是具有3次代数精度的插值型求积公式。 这个例子告诉我们，只要适当选择求积节点，可使插值型求积公式的代数精度达到最高。这就是本节要介绍的高斯求积公式。 同理，对于一般的插值求积公式 （3.10） 只要适当地选取其2n+2个待定参数 xk 和 ,就可使它的代数精度达到2n+1次。 定义3.3 若插值求积公式（3.10）具有2n+1次代数精度，则称之为高斯求积公式，并称相应的求积节点 为高斯点。 可以证明，n个节点的高斯求积公式具有最高不超过2n+1次的代数精度，这就是我们所要讨论的具有最高代数精度的插值型求积公式。 3.4几个低阶求积公式 在牛顿-柯特斯求积公式中n=1,2,4时，就分别 得到下面的梯形公式、辛卜生公式和柯特斯公式。(1)?梯形公式 当n=1时，牛顿-柯特斯公式就是梯形公式 定理3.2 （梯形公式的误差）设f(x)在[a,b]上具有连续的二阶导数，则梯形公式的误差（余项）为 证:由插值型求积公式的余项 其中 可知梯形公式的误差为 由于(x-a)(x-b)在[a,b]中不变号, 在[a,b]上连续,根据积分中值定理 ,在[a,b]上存在一点η，使 因此 （2） 辛卜生公式 当n=2时，牛顿-柯特斯公式就是辛卜生公式（或 称抛物线公式） 定理3.3（辛卜生公式的误差）设在[a,b]上具有连续的四阶导