[]newch5插值型数值微分与数值积分.pptVIP

第五章 插值型数值微分与数值积分 5.1 插值型数值微分公式 5.2 插值型数值积分　　 5.1 插值型数值微分公式 当 x 为插值节点 时，上式简化为 故一般限于对节点上的导数值采用插值多项式的相应导数值进行近似计算，f以便估计误差。 一般地 这类公式称为插值型数值微分公式。 5.1.1 常用的数值微分公式 1.两点公式(n=1) 这称为两点公式。 即 截断 误差 2.两点公式(n=2) 即 二阶导数（不要记忆） 例1：已知列表 X 2.5 2.55 2.60 2.65 2.70 Y 1.58114 1.59687 1.61245 1.62788 1.64317 解: h=0.05 例 5.1 为计算 在 x=2 处的一阶导数值，我们可选用中点公式 当计算保留四位小数时，得到计算结果如表5-1(书103页）。 而精确值为 ，可见当 h=0.1时近似结果最好，步长太大或太小计算效果均不好。 为估计二阶导数数值微分公式的误差，可设 f (x) 四阶连续可微，故得 从而得到误差估计式 5.2 插值型数值积分 插值型数值积分的思想是： 若已知 则利用拉格朗日插值多项式建立近似计算公式 这里 称为插值型求积公式， 称为求积节点， 称为求积系数，其和 下面求求积系数 ，设等距节点情形，即 ——牛顿-柯特斯公式 Cotes系数 特别地 这称为梯形公式； 几何意义：用梯形面积代替f(x)作为曲边的曲边梯形面积。 图1梯形公式 这称为Simpsion公式 图2 Simpson公式 几何意义：用抛物线替 作曲边的曲边梯形 面积代替f(x)作为曲边的曲边梯形面积。 这称为Cotes公式。 对应于 情形的Cotes系数见表5-2 (书106页)。 5.2.2 复合求积公式 求积公式的稳定性分析： 等距节点的插值求积公式，当n较大（n7）时，系数 中出现负数，而且有正有负，这将使舍入误差增大并难于 估计，因此实际计算时一般不用n较大的公式，而是将积 分区间(a,b)分成n个小区间，在每个小区间上用低阶New- Cotes公式计算积分的近似值，然后对这些近似值求和，从 而得到所求积分的近似值，由此得到一些有实际意义的求 积公式，称为复合求积公式。 1.复合梯形公式(n=1,简记为Tn ) 3.复合Cotes公式(n=4，简记为Cn )(公式见书107页） 2.复合Simpson公式(n=3,简记为Sn ) 2.确定h 解：1写出公式 例 1 计算 ，求 4.由表格计算结果 0 1 4 1 1 1 1 1.0625 3.7641 2 4 2 1.25 3.200 2 4 2 3 1.5625 2.5600 2 4 4 2 2 1 1 1 25.0494 18.8 37.6988 3.列表 例 2 试利用表5-3的函数表，分别用复合梯形公式、复合Simpson公式和复合Cotes公式计算定积分 解 三、求积公式的误差： 1.梯形公式误差： 大区间上的误差记为： 2.Simpson公式误差 不难推出 3.Cotes公式误差 四、变步长法则 1.基本思想 上面介绍的复化求积公式对提高进度是有效的，但是在使用求积公式之前，必须给出适当的步长。如果事先给出精度要求，在使用复化求积公式时，由于误差估计式中 含有