[工学]数值计算PPTcxj42.pptVIP

若对(6.8)右端积分用辛普森求积公式，则有 上式略去余项，并记 普森数值微分公式 的近似值为mk，则得到辛 这是关于m0，m1，…，mn这 n+1 个未知量的 n-1 个方程组， 可用解方程的方法求解. * 上页 下页 返回 第四章 数值积分与数值微分（二） 第四节 龙贝格求积公式 第五节 高斯求积公式 第六节 数值微分 一、梯形法的递推化 上一节介绍的复化求积方法对提高精度是行之有效的，但在使用求积公式之前必须给出合适的步长，步长取得太大精度难以保证，步长太小则会导致计算量的增加，而事先给出一个恰当的步长又往往是困难的． 实际计算中常常采用变步长的计算方案，即在步长逐次分半(即步长二分)的过程中，反复利用复化求积公式进行计算，直至所求得的积分值满足精度要求为止． 设将求积区间[a，b]分成n等分，则一共有n+1个分点，按梯形公式计算积分值Tn，需要提供n+1个函数值．如果将求积区间再二分一次，则分点增至2n+1个，我们将二分前后两个积分值联系起来加以考察． §4 龙贝格求积公式 注意到每个子区间[xk，xk+1]经过二分只增加了一个分点xk+1/2＝( xk+xk+1)/2，用复化梯形公式求得该子区间上的积分值为 二、龙贝格算法 根据复化梯形公式的余项表达式 可见，利用两种步长计算的结果能估计截断误差.若将该截断误差加到计算结果中，就可得出“改进梯形求积公式”： 改进梯形求积公式的右边实际是 这就是说用梯形法二分前后的两个积分值Tn与T2n的线性组合的结果得到复化辛普森法求积公式 类似的情况，用辛普森法二分前后的两个积分值Sn与S2n的线性组合的结果可得到复化柯特斯求积公式 重复同样的手续，用柯特斯法二分前后的两个积分值Cn与C2n的线性组合的结果可得到龙贝格(Romberg)求积公式 我们在变步长的过程中运用加速公式(1)、(2)、(3)，就能将粗糙的梯形值Tn逐步加工成精度较高的辛普森值Sn 、柯特斯值Cn和龙贝格值Rn . 龙贝格求积算法可用下表来表示： 三、理查森(Richardson)外推加速法 上面讨论说明由梯形公式出发，将区间[a，b]逐次二分可提高求积公式的精度，上述加速过程还可继续下去，其理论依据是梯形公式的余项展开，即 若记Tn = T(h)，当区间[a，b]分为2n等分时，有 ，则 可见I = T(h)的误差为O(h2)阶.若记 ，则 显然T1(h)与积分值 I 近似的阶为O(h4) .这样构造的 就是辛普森公式序列Sn，S2n，… .若令 ， 则又可进一步从余项中消去 h4 项，这样构造出的 ，其实就是柯特斯公式序列，它与积分值 I 的逼近阶为 O(h6) . 如此继续下去，每加速一次，误差的量级便提高2阶，一般地，若记T0(h) = T(h)，经过m (m = 1，2，…)次加速后，则有 可以证明，如果 f (x) 充分光滑，那么T 数表每一列的元素及对角线元素均收敛到所求的积分值 I ，即 机械求积公式 含有2n+2个待定参数xk、Ak(k＝0，1，…，n)．当 xk 为等距节点时得到的插值求积公式的代数精度至少为n次，如果适当选取 xk (k＝0，1，…，n)，有可能使求积公式具有 2n+1 次代数精度，这类求积公式称为高斯(Gauss)求积公式. 为使问题更具一般性，我们研究带权积分 ，这里r (x)为权函数，类似机械求积公式，它的求积公式为 Ak(k＝0，1，…，n)为不依赖于f (x)的求积系数，xk (k＝0，1，…，n)为求积节点，可适当选取 xk 及 Ak (k＝0，1，…，n)使积分(1)具有2n+1次代数精度． §5 高斯求积公式 一、高斯点 定义4 如果求积公式(5.1)具有2n+1次代数精度，则称其节点 xk (k＝0，1，…，n)为高斯点，相应公式(5.1)称为高斯求积公式 . 根据定义要使(5．1)具有2n+1次代数精度，只要取f(x)＝xm， 对m＝0，1，…，2n+1，(5.1)精确成立，则得 当给定权函数r (x)，求出右端积分，则可由(5.2)解得 Ak 及 xk (k ＝