第六章 插值.pptVIP

华长生制作 1. 线性插值 2. 抛物线插值 3. n 次Lagrange插值 Lagrange插值公式的标准型公式: Lagrange插值算法特点局限性 例：对于如下数据，试求牛顿插值多项式 解：先计算各阶均差见下表 故 85 40 15 4 1 yi 4 3 2 1 0 xi 0 1 10 45 85 4 4 1 7 25 40 3 3 4 11 15 2 2 3 4 1 1 1 0 0 四阶均差 三阶均差 二阶均差 一阶均差 yi xi 序 此题如用拉格朗日插值，则有 结果与牛顿插值一样，但较为繁琐。 ? 等距节点公式 向前差分 i i i y y y - = ? + 1 i k i k i k i k y y y y 1 1 1 1 ) ( - + - - ? - ? = ? ? = ? 向后差分 1 1 1 - - - ? - ? = ? i k i k i k y y y i?1 i i y y y - = ? 当节点等距分布时: y0 y1 y2 y3 ? y0 y1 y2 ?2y0 ?2y1 ?3y0 ? 节点等距情况下： 步长h=x i+1-x i 一般地： 可见： 在节点等距情况下，牛顿插值公式中的差商可换成相应的 差分，其形式可得到进一步简化。 牛顿公式 6.3 Hermite插值法 , , , , , , , ) ( 1 0 1 0 n n y y y b x x x a x f L L 处的函数值为 在节点 设 ￡ ￡ 值函数 上的具有一阶导数的插 的在区间 为 设 ] , [ ) ( ) ( b a x f x P ) ( ] , [ ) ( ) 1 ( 一阶光滑度 上具有一阶导数 在 若要求 b a x P 埃尔米特插值 /* Hermite Interpolation */ 不仅要求函数值重合，而且要求若干阶导数也重合。 即：要求插值函数 P (x) 满足P(xi) = f (xi), P’ (xi) = f ’ (xi), …, P(m) (xi) = f (m) (xi). 注：? N 个条件可以确定 阶多项式。 N ? 1 ?要求在1个节点 x0 处直到m 阶导数都重合的插值多项式即为Taylor多项式 其余项为 ?一般只考虑 f 与f ’的值。 两点三次Hermite插值 先考虑只有两个节点的插值问题 希望插值系数与Lagrange插值一样简单 重新假设 其中 可知 由 可得 Lagrange 插值基函数 类似可得 即 将以上结果代入 得两个节点的三次Hermite插值公式 两点三次Hermite插值的余项 两点三次Hermite插值的误差为 构造辅助函数 均是 二重根 连续使用4次Rolle定理，可得， 使得 即 所以,两点三次Hermite插值的余项为 以上分析都能成立吗？ 例. 解: 作为多项式插值,三次已是较高的次数，次数再高就有 可能发生Runge现象,因此，对有n+1节点的插值问题， 我们可以使用分段两点三次Hermite插值 * * § 6.1 Lagrange插值多项式 第六章 函数近似计算的插值法 若通过求解线性方程组(1)来求解插值多项式 系数 , 不但计算工作量较大, 且难于得到 的简单表达式. 一、 代数多项式的构造: 可通过找插值基函数的方法,得到插值多项式! 十八世纪法国数学家Lagrange对以往的插值算法进 行研究与整理，提出了易于掌握和计算的统一公式， 称为Lagrange插值公式。 它的特例是线性插值公式和抛物线插值公式。 Lagrange插值多项式 已知两个插值点及其函数值： f1 f0 f(x) x1 x0 x 插值节点 对应的函数值 求一次多项式 使得 由于方程组的系数行列式 所以，按Gramer法则，有唯一解 于是 或 (B-1） 容易验证，过点（x0，f0）与（x1，f1）直线方程就是式（B-1），如图5-3所示。 y x x0 x1 P1(x) f(x) P1(x) f(x) 误差 图5-3 已知三个插值节点及其函数值： f2 f1 f0 f(x) x2 x1 x0 x 求一个二次多项式 使得 由于该方程组的系数行列式 所以，有唯一解。即满足这样条件的二次多项式是唯一确定的。 满足上述条件，所以它就是所求的二次多项式。 容易看出 （B-2） 容易验证，P2(x)是过点(x0, f0)、(x1, f1)与(x2, f2)三点的抛物线，如图5-4所示。 y x x1 x0 x2 P2