数值分析(研究生)第三章数值积分及数值微分二.ppt

三、利用数值积分求导 微分是积分的逆运算，因此可利用数值积分的方法来计算数 值微分．设 f (x)是一个充分光滑的函数，设 则有： 对上式右边积分采用不同的求积公式就可得到不同的数值微分公式．例如，对 用中矩形公式 则得 从而得到中点微分公式 若对(7.6)右端积分用辛普森求积公式，则有 上式略去余项，并记 普森数值微分公式 的近似值为mk，则得到辛 这是关于m0，m1，…，mn这 n+1 个未知量的 n-1 个方程组， 可用解方程的方法求解. 四、三次样条求导 三次样条函数 S(x) 作为 f (x)的近似，不但函数值很接近，导数值也很接近，并有 (见第2章定理4)，因此利用三次样条函数S(x)直接得到 根据第2章(6.6) 可求得 式中 f [xk，xk+1]为一阶均差 . 五、数值微分的外推算法 (6.8) (6.8) 5 3 - 1 (6.8) 3 - 1 * 上页 下页 返回 第三章 数值积分与数值微分（二） 第五节 Romberg求积算法 第六节 Gauss求积公式 第七节 数值微分 一、梯形公式的递推公式及事后估计法 上一节介绍的复化求积方法对提高精度是行之有效的，但在使用求积公式之前必须给出合适的步长，步长取得太大精度难以保证，步长太小则会导致计算量的增加，而事先给出一个恰当的步长又往往是困难的． 实际计算中常常采用变步长的计算方案，即在步长逐次分半(即步长二分)的过程中，反复利用复化求积公式进行计算，直至所求得的积分值满足精度要求为止． 设将求积区间[a，b]分成n等分，则一共有n+1个分点，按梯形公式计算积分值Tn，需要提供n+1个函数值．如果将求积区间再二分一次，则分点增至2n+1个，我们将二分前后两个积分值联系起来加以考察． §5 龙贝格求积公式 注意到每个子区间[xk，xk+1]经过二分只增加了一个分点xk+1/2＝( xk+xk+1)/2，用复化梯形公式求得该子区间上的积分值为 二、龙贝格算法 根据复化梯形公式的余项表达式 可见，利用两种步长计算的结果能估计截断误差.这种利用计算结果估计误差的方法称为事后估计法. 若将该截断误差加到计算结果中，就可得出“改进梯形求积公式”： 改进梯形求积公式的右边实际是 这就是说用梯形法二分前后的两个积分值Tn与T2n的线性组合的结果得到复化辛普森法求积公式 类似的情况，用辛普森法二分前后的两个积分值Sn与S2n的线性组合的结果可得到复化柯特斯求积公式 重复同样的手续，用柯特斯法二分前后的两个积分值Cn与C2n的线性组合的结果可得到龙贝格(Romberg)求积公式 我们在变步长的过程中运用加速公式(5.1)、(5.2)、(5.3)，就能将粗糙的梯形值Tn逐步加工成精度较高的辛普森值Sn 、柯特斯值Cn和龙贝格值Rn . 龙贝格求积算法可用下表来表示： 例2 用龙贝格方法计算椭圆 x2/4 + y2 ＝ l 的周长，使结果具有五位有效数字． 分析 为便于计算，先将椭圆方程采用参数形式表示，再根据弧长公式将椭圆周长用积分形式表示．由于计算结果要求具有五位有效数字，因此需要估计所求积分值有几位整数，从而确定所求积分值的绝对误差限．最后再应用龙贝格方法计算积分． 解 令 x ＝2cosq，y ＝sinq 则椭圆的周长为 三、理查森(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 ，即 机械求积公式