第4章 插值法2.ppt

* §3.牛顿(Newton)插值 3.1差商及其性质 一.差商（均差）定义 拉格朗日插值公式可看作直线方程两点式的 推广，若从直线方程点斜式 出发，将它推广到具有n+1个插值点的情况，可 把插值多项式表示为 其中 为待定系数，可由插值条件 确定。 当 依次可得到 。为写出系数 的一般表达式，现引入差商（均差）定义。 定义：称 为函数 关于节点 的一阶差商，记为 。 一阶差商 的差商 称为 关于节点 的二阶差商，记 为 。 递归地用k-1阶差商来定义k阶差商， 称为 关于k+1个节点 的k阶 差商。 性质1:k阶差商可以表示成k+1个函数值 的线性组合,即 可用归纳法证明。 例: 二. 差商（均差）的性质 这个性质也表明差商与节点的排列顺序无关(差商的对 称性)。即 性质2： k阶差商定义: 依对称性，对调定义公式左端k阶差商中 与 的位置， 再将各差商中的节点按原来次序排列。 性质3:若 是 的n次多项式,则一阶差商 是 的n-1次多项式,二阶差商 是 的n-2次多项式; 一般地,函数 的k阶差商 是 的n-k次多项式 ,而 时,k阶 差商为零。 若 是 的n次多项式，则 也是n次多项式，且 。于是 可分 解为 其中 为n-1次多项式。所以 为n-1次多项式 三.利用差商表计算差商 利用差商的递推定义,可以用递推来计算差商。如下表： 如要计算四阶差商,再增加一个节点,表中还要增加一行。 一阶差商 二阶差商 三阶差商 ? ? ? ? ? ? ┇ ┇ ┇ ┇ ┇ 例:已知如下,计算三阶差商 。 解：列表计算 1 3 4 7 0 2 15 12 一阶差商 二阶差商 三阶差商 1 0 ? ? ? 3 2 1 ? ? 4 15 13 4 ? 7 12 -1 -3.5 -1.25 3.2 牛顿插值公式 根据差商定义，把 看成 上的一点，可得 ┅ 只要把后一式代入前一式，得: 最后一项中, 差商部分含有 ,为余项部分,记作 而前面n+1项中, 差商部分都不含有 ,因而前面n+1项 是关于 的n次多项式,记作 这就是牛顿插值公式。于是,上式记为 由牛顿插值公式与 比较知： 例如:当n=1时， 其中, 这就是牛顿一次插值多项式，也就是点斜 式直线方程。 当n=2时, 这就是牛顿二次插值多项式。 显然, 即 满足二次插值条件。 若 是 上n次连续可微函数， 并且 是 中不同的点，则在 中存在一点 ，使得 设 是函数 在节点 上次数 至多为n-1次的插值多项式，则 例:已知 求满足以上插值条件的牛顿型插值多项式。 解:在上例中,我们已计算出 则牛顿三次插值多项式为 1 3 4 7 0 2 15 12 例3:已知 在六个点的函数值如下表，运用牛顿型插值多项式求 的近似值。 一阶 二阶 三阶 四阶 五阶 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 解： 欲求 ,只需在 之后再加一项: 故 。 截断误差： *