[所有分类]第6章 函数插值-2.pptVIP

? 二、差商(亦称均差) /* divided difference */ 此时 若 则 例: 已知 y = f (x) 的观测数据如下 求分段线性插值和分段二次插值 P1(x) 和 P2(x) . -1 -5 3 1 2 4 8 f (x) 0 1 2 3 4 5 6 x §6. 6 Hermite 插值 分段线性插值简单易操作, 但插值曲线不光滑, 即在内节点处一节导数不连续, 这种情况往往不能满足实际应用的需要. 为了克服这一缺陷, 通常添加一阶导数作为插值条件 . 一、 两个节点的情形： 设 x0 , x1为插值节点， x0 x1，且已知 在区间[x0 , x1]上求多项式 H (x)，使得满足插值条件 由于有4个条件，所以H (x)应为次数不超过3次的多项式，称为 Hermite三次插值。 定理1：设 ，则在区间[x0 , x1]上满足插值条件 的不超过3次的多项式H (x)存在且唯一，并可构造如下： 其中插值基函数 为三次多项式 如果 ，则插值余项为 其中 在 x0 与 x1 之间。 插值基函数的性质: 插值函数的唯一性： 设H (x)和 都是满足插值条件的不超过3次的Hermite插值多项式，则 是不超过3次的多项式，且满足 这说明 x0 和 x1 都是P (x)的二重根，从而P (x)为4次多项式，这是不可能的。 * §6.3 牛顿插值(Newton’s Interpolation) Lagrange 插值虽然易算，但若要增加一个节点时，全部基函数 l i ( x ) 都需要重新计算。也就是说，Lagrange 插值不具有继承性。 能否重新在Pn中寻找新的基函数 ？ 希望每加一个节点时，只在原有插值的基础上附加部分计算量（或者说添加一项）即可。 本讲主要内容: ● Newton插值多项式的构造 ● 差商的定义及性质 ● 差分的定义及性质 ● 等距节点Newton插值公式 问题1 求作 n 次多项式 (2) 为了使 的形式得到简化,引入如下记号 (3) (1) 一、 基函数 使满足 De f 1: 由式(3)定义的n+1个多项式 称为Newton插值的以x0 , x1,…, xn 为节点的基函数, 即 可以证明,这样选取的基函数是线性无关的, 由此得出的问题1的解便于求值, 而且新增加一个节点 xn+1时 只需加一个新项 即 而 依据条件(2),可以依次确定系数 c0 , c1,…, cn . 如: 取 x = x0 得 取 x = x2 得 为了得到计算系数 ci 的一般方法,下面引进差商的概念. 称为f (x)在x0 , x1处的1阶差商 称为f (x)在x0 , x1 , x2 处的2阶差商 一般地，n 阶差商： De f 2: 称为f (x)在x1 , x2 处的1阶差商 给定[a , b]中互不相同的点 x0, x1, x2 ,…, 以及 f (x)在这些点处相应的函数值 f (x0) , f (x1) , f (x2) ,…,则 差商的性质: 性质1 (差商与函数值的关系) 证: 归纳法. 当k=1时, 有 结论成立. 设 k = m -1 时, 结论成立. 则有 由差商的定义, 当k = m 时, 有 所以k = m 时结论成立, 由归纳假设知性质成立. 性质2 (对称性): 差商的值与结点排列顺序无关. 在n阶差商中任意调整节点的顺序, 差商的值不变. 性质3 (差商与导数的关系) 提示: 反复应用Rolle定理. 性质4: (重节点差商) 若 f (x)在 xi 处具有k 阶导数, 则有 性质5: (线性性) 若 为常数, 则有 规定：一个点 x0 的零阶差商为 f (x0) . 例：设 求 和 差商表 xk f (xk) 一阶差商 二阶差商 三阶差商 ……