2.3牛顿插值程序.pptVIP

* §2.3 均差与牛顿插值多项式 拉格朗日插值多项式结构对称，使用方便。但由于是用基函数构成的插值，这样要增加一个节点时，所有的基函数必须全部重新计算，不具备承袭性，还造成计算量的浪费。这就启发我们去构造一种具有承袭性的插值多项式来克服这个缺点，也就是说，每增加一个节点时，只需增加相应的一项即可。这就是牛顿插值多项式。 3.1 差商（均差）及性质 1 差商（均差） 已知y = 函数表 则 在 上平均变化率分别为： 即有定义： 定义为f(x) 的差商 §2.3 差商与牛顿插值多项式 . , ) ( ) ( ) ( ] , [ 2 0 0 0 0 的一阶（均差） 关于点 为函数 称 定义 k k k k x x x f x x x f x f x x f - - = 特别地 均差及其性质 2 基本性质 （2）k 阶差商 关于节点 是对称的，或说 均差与节点顺序无关，即 例如： 共6个 的线性组合，即 的k阶差商 是函数值 （1） 分析 ： 当k =1时, (1)可用归纳法证明。(2)利用(1)很容易得到。只证(1) 证明： （1）当k =1时, 一阶差商 二阶差商 k 阶差商 三阶差商 (0 阶差商) 表2.1 3 差商表 计算顺序:每次用前一列同行的差商与前一列 上一行的差商再作差商。 3.2 牛顿插值多项式 已知 函数表（3.1）, 由差商定义及对称性，得 1 牛顿插值多项式的推导 将(b)式两边同乘以, 抵消 抵消 抵消 (d)式两边同乘以 ,把所有式子相加,得 ,(c)式两边同乘以 记 --- 牛顿插值多项式 --- 牛顿插值余项 可以验证 ，即 满足插值条件, 因此 可得以下结论。 定理 则满足插值条件 的插值多项式为： （牛顿插值多项式） 其中， --- 牛顿插值多项式 --- 牛顿插值余项 2 n +1阶差商函数与导数的关系 由n次插值多项式的唯一性，则有 , 牛顿插值多项式 与拉格朗日插值多项式 都是次数小于或等于n的多项式, 只是表达方式不同. ? 因为 而 的基函数可为: 已知 函数表 牛顿插值多项式系数 牛顿插值多项式系数 牛顿插值多项式系数 阶导数存在时，由插值多项式的唯一性有余项公式 n+1阶差商函数 导数 其中 且 为包含 区间. 依赖于 则n 阶差商与导数 的关系为 其中 n +1阶差商函数与导数的关系 定理 。 Newton插值的优点是：每增加一个节点，插值多项式只增加一项，即： 因此便于递推运算。而且Newton插值的计算量小于Lagrange插值。 例1 解析 一阶均差 二阶均差 三阶均差 四阶均差 五阶均差 1.11600 1.18600 0.28000 1.27573 0.35893 0.19733 1.38410 0.43348 0.21300 0.03134 1.51533 0.52492 0.22863 0.03126 -0.00012 一阶均差 二阶均差 三阶均差 四阶均差 五阶均差 0.40 0.41075 0.55 0.57815 0.65 0.69675 0.80 0.88811 0.90 1.02652 1.05 0.25382 于是 ( ) ( ) 63192 . 0 596 . 0 596 . 0 4 = ? N f 3 3 0 1 -1 -1/3 -2 -3/2 -1/6 1/24 1 4 7 8 6 1 2 3 4 5 0 1 2 3 4 一阶差商 二阶差商 三阶差商 四阶差商 f(xk) xk k 定义 3.3 差 分 向前差分算子， 差 分 二阶中心差分为 差 分 性质1 各阶差分可用函数值表示 性质2 差 分 性质3 差 分 差分与导数的关系： 差 分 差 分 3.4 等距节点插值公式 此公式为牛顿向前插值公式，其余项为 等距节点插值公式 计算各阶差分，可列表进行 等距节点插值公式 类似有牛顿向后插值公式 等距节点插值公式 差分表： 等距节点插值公式 等距节点插值公式:牛顿向前插值公式、牛顿向后插值公式。 等距节点插值公式 已知f(x)=sinx的函