第五章插值型数值微分与数值积分精选文档.ppt

《第五章插值型数值微分与数值积分精选文档.ppt》由会员分享，可在线阅读，更多相关《第五章插值型数值微分与数值积分精选文档.ppt（39页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第五章插值型数值微分与数值积分本讲稿第一页，共三十九页5.1 插值型数值微分公式插值型数值微分公式 故一般限于对节点上的导数值采用插值多项式的相应导数值进行故一般限于对节点上的导数值采用插值多项式的相应导数值进行近似计算，以便估计误差。近似计算，以便估计误差。这类公式称为这类公式称为插值型数值微分公式插值型数值微分公式。当当 x 为插值节点为插值节点xi 时，上式简化为时，上式简化为 一般地一般地 本讲稿第二页，共三十九页 5.1.1 常用的数值微分公式常用的数值微分公式 即 1.两点公式两点公式(n=1)这称为这称为两点公式两点公式。本讲稿第三页，共三十九页两点公式的截断误差为两点公式的截断

2、误差为 这里这里 本讲稿第四页，共三十九页2.三点公式三点公式(n=2)这称为这称为三点公式三点公式，其中（，其中（54b）又称为）又称为中点公式中点公式。本讲稿第五页，共三十九页三点公式的截断误差为三点公式的截断误差为 这里这里 本讲稿第六页，共三十九页进一步由进一步由 可得计算公式可得计算公式 本讲稿第七页，共三十九页为估计二阶导数数值微分公式的误差，可设为估计二阶导数数值微分公式的误差，可设 f(x)四阶连续可微，四阶连续可微，故得故得 从而得到误差估计式从而得到误差估计式 二阶导数的截断误差二阶导数的截断误差本讲稿第八页，共三十九页 例例1：已知列表：已知列表x 2.5 2.55 2.

3、60 2.65 2.70y 1.58114 1.59687 1.61245 1.62788 1.64317解解:h=0.05本讲稿第九页，共三十九页 例例 5.1 为计算为计算 在在 x=2 处的一阶导数值，我们可选用中处的一阶导数值，我们可选用中点公式点公式当计算保留四位小数时，得到计算结果如表当计算保留四位小数时，得到计算结果如表5-1(书书103页）。页）。而精确值为而精确值为 ，可见当，可见当 h=0.1时近似结果最好，时近似结果最好，步长太大或太小计算效果均不好。步长太大或太小计算效果均不好。本讲稿第十页，共三十九页5.2 插值型数值积分插值型数值积分 x0 x1-xi-1xixi+

4、1-xnf(x0)f(x1)-f(xi-1)f(xi)f(xi+1)-f(xn)2.由下列列表函数求由下列列表函数求L-插值多项式插值多项式称为称为插值型求积公式插值型求积公式，称为称为求积求积节点节点，称为称为求积系数求积系数，其和 本讲稿第十一页，共三十九页 5.2.1 Newton-CotesNewton-Cotes公式公式 则则 ，考虑等距节点的情形考虑等距节点的情形考虑等距节点的情形考虑等距节点的情形 牛顿牛顿-柯特斯公式柯特斯公式CotesCotes系数系数 本讲稿第十二页，共三十九页n=1,2,4n=1,2,4的的的的N-CN-C公式公式公式公式 这称为这称为梯形公式梯形公式；几

5、何意义：用梯形面积几何意义：用梯形面积代替代替f(x)作为曲边的曲边作为曲边的曲边梯形面积。梯形面积。图图1 梯形公式梯形公式 ab本讲稿第十三页，共三十九页这称为这称为Simpsion公式公式图图2 Simpson公式公式 ab几何意义：用抛物线几何意义：用抛物线 作曲边的曲边作曲边的曲边梯形面积代替梯形面积代替f(x)作作为曲边的曲边梯形面积。为曲边的曲边梯形面积。本讲稿第十四页，共三十九页对应于对应于 情形的情形的Cotes系数见表系数见表5-2(书书106页页)。这称为这称为Cotes公式公式。求积公式的稳定性分析求积公式的稳定性分析求积公式的稳定性分析求积公式的稳定性分析 本讲稿第十

6、五页，共三十九页5.2.2 复合求积公式复合求积公式 本讲稿第十六页，共三十九页当取当取 m=1 时，称为时，称为复合梯形公式复合梯形公式，简记为，简记为Tn1.复合梯形公式复合梯形公式=1为本讲稿第十七页，共三十九页当取当取 m=2 时，称为时，称为复合复合Simpson公式公式，简记为，简记为Sn2.复合复合Simpson公式公式本讲稿第十八页，共三十九页当取当取 m=4 时，称为时，称为复合复合Cotes公式公式，简记为，简记为Cn(公式见书公式见书107页）页）3.复合复合Cotes公式公式本讲稿第十九页，共三十九页 例例 5.2 试利用表试利用表5-3的函数表，分别用复合梯形公式、复

7、合的函数表，分别用复合梯形公式、复合Simpson公式和复合公式和复合Cotes公式计算定积分公式计算定积分 解：解：1.写出公式写出公式本讲稿第二十页，共三十九页2.确定确定h本讲稿第二十一页，共三十九页3.列表列表 k xk f(xk)T8 S4 C2 0 0 0.000 000 1 1 7 1 1/8 0.110 312 2 4 32 2 1/4 0.194 700 2 2 12 3 3/8 0.257 733 2 4 32 4 1/2 0.303 265 2 2 14 5 5/8 0.334 538 2 4 32 6 3/4 0.354 275 2 2 12 7 7/8 0.364 7

8、54 2 4 32 8 1 0.367 879 1 1 7 4.20704 6.341712 47.56338 本讲稿第二十二页，共三十九页本讲稿第二十三页，共三十九页5.2.3 插值型求积公式的误差分析与步长减半算法插值型求积公式的误差分析与步长减半算法 1、求积公式的误差、求积公式的误差 记记 为采用插值型求积公式进行积分近似的为采用插值型求积公式进行积分近似的截断误差，则由多项式插值公式的误差估计式（截断误差，则由多项式插值公式的误差估计式（5-1）得）得 因此，当因此，当 f(x)为次数不超过为次数不超过 n 次的多项式时，插值型求积公式精确次的多项式时，插值型求积公式精确成立。成立。

9、本讲稿第二十四页，共三十九页特别地，特别地，本讲稿第二十五页，共三十九页从而可得从而可得 本讲稿第二十六页，共三十九页 为便于估计误差，实际计算时常常采用步长逐次减半的算法，为便于估计误差，实际计算时常常采用步长逐次减半的算法，下面介绍其思想。下面介绍其思想。由由(5-17a)得得所以所以 2、变步长法则、变步长法则本讲稿第二十七页，共三十九页 类似地，可对类似地，可对 Simpson 公式和公式和 Cotes 公式分别利用公式分别利用（5-18b）和（）和（5-18c）进行事后误差估计，建立步长逐次减）进行事后误差估计，建立步长逐次减半的算法。半的算法。因此，可先用因此，可先用 计算出计算出

10、T1，并把步长，并把步长减半算出减半算出T2，若，若 则则T2 即为即为所求的近似值，否则再把步长减半，算出所求的近似值，否则再把步长减半，算出T4，根据式根据式(5-18a)进行进行事后误差估计事后误差估计 ，如此递推计算，直如此递推计算，直到某个到某个n 满足满足 为止为止，取，取 为所求的近似为所求的近似值，这就是梯形公式的步长逐次减半算法。值，这就是梯形公式的步长逐次减半算法。本讲稿第二十八页，共三十九页为减少计算量，需建立递推公式，现对复合梯形公式推导之。为减少计算量，需建立递推公式，现对复合梯形公式推导之。这里这里 对应于新的步长，对应于新的步长，对对应于新分点。应于新分点。本讲稿

11、第二十九页，共三十九页因此可建立梯形公式的步长逐次减半递推公式：因此可建立梯形公式的步长逐次减半递推公式：本讲稿第三十页，共三十九页解解:本讲稿第三十一页，共三十九页计算结果见下表计算结果见下表本讲稿第三十二页，共三十九页 例例 5.3 试用梯形公式的步长逐次减半算法计算定积分试用梯形公式的步长逐次减半算法计算定积分 使误差小于使误差小于 。解解 一般的计算结果见表一般的计算结果见表5-4(书书112页页)。本讲稿第三十三页，共三十九页 5.2.4 龙贝格积分法龙贝格积分法 这说明收敛较快的这说明收敛较快的 Simpson 步长减半序列步长减半序列 可由梯形公式的可由梯形公式的步长减半序列步长

12、减半序列 构造生成。构造生成。本讲稿第三十四页，共三十九页类似地，类似地，(5-20c)称为称为龙贝格（龙贝格（Romberg）积分公式）积分公式。按以上方法可继续。按以上方法可继续外推下去，建立如下收敛较快的外推算法外推下去，建立如下收敛较快的外推算法龙贝格积分法龙贝格积分法（书（书114页）。页）。本讲稿第三十五页，共三十九页 其计算公式为其计算公式为 注：这样的计算格式注：这样的计算格式可根据精度自动停机。只可根据精度自动停机。只要竖线上相邻两结果之差要竖线上相邻两结果之差不超过给定精度为止。计不超过给定精度为止。计算过程实质是将区间逐次算过程实质是将区间逐次分半计算分半计算 ，然后利用，然后利用加速公式，故又叫逐次分加速公式，故又叫逐次分半加速法。半加速法。本讲稿第三十六页，共三十九页解解:本讲稿第三十七页，共三十九页本讲稿第三十八页，共三十九页本讲稿第三十九页，共三十九页