数值分析 第4章 数值积分与数值微分.pptVIP

* * * * * * * * * * * * * * * 我们看到，中点公式具有如下形式 (8.7) 式中的系数均与步长无关。 若将步长二分，则有 (8.8) 则可消去误差主项 ，得 取(8.7)与(8.8)加权平均 (8.9) 重复同样的手续，再导出下列加速公式 这种加速过程还可继续下去。这种加速方法通常称作Richardson外推加速法。 若令 (8.10) 则进一步消去D1(h)误差主项 , 有 例8.2 运用加速公式和加工例8.1的结果。 表8.2 Richardson外推加速法计算结果 h D(h) 0.8 0.4 0.2 0.1 3.01765 2.79135 2.73644 2.72281 2.715917 2.718137 2.718267 2.718285 2.718276 2.71828 解 计算结果见表8.2。这里，加速的效果同样是相当显著的。 注: 中点心公式及其加速方法适合用表达式表示的函数。对于列表函数，则宜使用插值方法等导出数值求导公式。 最小,步长h不宜太大,也不宜太小. 其最优步长应为 它表明h 越小, 舍入误差 越大, 故它是病态的. 用中点 公式计算 的误差上界为 要使误差E(h) 则计算 当 f (a+h) 及 f (a-h) 分别有舍入误差 ?1 及 ?2 时，若令 的舍入误差上界为 二、插值型的求导公式 x x0 x1 x2 … xn y y0 y1 y2 … yn 对于列表函数 y = f (x)： 插值多项式 y = Pn(x)作为它的近似, 我们取 统称插值型的求导公式． 的近似值， 作为 建立的数值公式 依据插值余项定理，求导公式(8.11)的余项为 式中 (8.11) 我们限定:求某个节点 xk 上的导数值，上面的第二项变为零，这时有余项公式 (8.12) 1．两点公式 已给两节点 x0, x1 上的函数值 f (x0), f (x1)，做线性插值 记 x1 – x0 = h，对上式两端求导，有 于是有下列求导公式： 下面我们仅仅考察节点处的导数值．为简化讨论, 假定所给的节点是等距的． 而利用余项公式(8. 12)知，带余项的两点公式是： 设已给出三节点x0, xl=x0+h, x2=x0+2h上的函数值,做二次插值 令 x = x0+th, 则 2．三点公式 (8.13) 这里撇号(’)表示对x求导．上式分别取 t = 0,1,2, 得到三种三点公式： (8.14) 而带余项的三点求导公式如下： 公式(8.6)是我们所熟悉的中点公式．在三点公式中，它由于少用了一个函数值 f(x1) 而引人注目． 例如，将式(6.5)再对t求导一次，有 于是有 用插值多项式Pn(x)作为 f (x)的近似函数，还可以建立高阶数值微分公式： 而带余项的二阶三点公式如下： (8.15) 例 已知函数y= ex 的下列数值: