【推荐】数值计算方法：-数值微分与数值积分.pptVIP

当取 m =1 时，称为 复化梯形公式 ，简记为 T n （ 1 ）复化梯形公式 31 当取 m =2 时，称为 复化 Simpson 公式 ，简记为 S n （ 2 ）复化 Simpson 公式 32 复合求积公式的误差（余项） 33 解： 使用复化梯形公式 使用复化 Simpson 公式： 例 5. 3 分别用复化梯形公式、复化 Simpson 公式计算定积分 34 为便于估计误差，实际计算时常常采用步长逐次减半的算 法，下面介绍其思想。 由 得 5.2.3 Romberg 积分法 （ 1 ）梯形逐步减半算法 35 所以得 所以 36 第五章 数值微分与数值积分 ? 5.1 数值微分 ? 5.2 数值积分 1 微积分中，关于导数的定义如下： 自然，而又简单的方法就是， 取极限的近似值，即差商 . 5.1 数值微分 5.1.1 差商型求导公式 2 由 Taylor 展开 因此，有误差 向前差商 5.1 数值微分 5.1.1 差商型求导公式 3 由 Taylor 展开 因此，有误差 向后差商 5.1 数值微分 5.1.1 差商型求导公式 4 由 Taylor 展开 因此，有误差 中心差商 5.1 数值微分 5.1.1 差商型求导公式 5 由误差表达式， h 越小，误差越小，但同时舍入误差增大， 所以，有个 最佳步长 我们可以用事后误差估计的方法来确定 设 D(h),D(h/2) 分别为步长为 h,h/2 的差商公式。则 时的步长 h/2 就是合适的步长 6 f(x)=exp(x) h f(1.15) R(x) h f(1.15) R(x) 0.10 3.1630 -0.0048 0.05 3.1590 -0.0008 0.09 3.1622 -0.0040 0.04 3.1588 -0.0006 0.08 3.1613 -0.0031 0.03 3.1583 -0.0001 0.07 3.1607 -0.0025 0.02 3.1575 -0.0007 0.06 3.1600 -0.0018 0.01 3.1550 -0.0032 例： 7 插值是建立逼近函数的手段，用以研究原函数的性质。 因此，可以用插值函数的导数近似为原函数的导数。 误差 5.1 数值微分 5.1.2 插值型求导公式 注意 ： 为了便于估计误差，限定只能对节点上的导数值采用插值 多项式的相应导数进行近似。 8 1 、两点公式 5.1 数值微分 5.1.2 插值型求导公式 给定两点上的函数值 这称为 两点公式 。 9 截断误差 : 10 若给定三点上的函数值 则由 这称为 三点公式 ，其中（ 5.1.9 ）又称为 中点公式 。 2 、 三点公式 11 例 1 ：已知列表 X 2.5 2.55 2.60 2.65 2.70 Y 1.58114 1.59687 1.61245 1.62788 1.64317 12 解 : h=0.05 13 例 5.1 为计算 在 x =2 处的一阶导数值，我们 可选用中点公式 而精确值为 ， 经计算，发现：当 步长 h =0.1 时近似结果最好，步长太大或 太小计算效果均不好。 14 二阶导数公式及误差 对其求二阶导数得 由 Taylar 展开可得误差估计式 15 5.2 数值积分 5.2.1 插值型求积公式 （ 1 ）插值型求积公式 （ 2 ） Newton-Cotes 型求积公式 （ 3 ）梯形公式、 Simpson 公式和 Cotes 公式 5.2.2 复化求积公式 （ 1 ）复化梯形公式 （ 2 ）复化 Simpson 公式 5.2.3 Romberg 积分