中学课件数值分析牛顿插值.pptVIP

?牛顿插值 (Newton’s Interpolation) Lagrange 插值多项式具有结构紧凑、便于编程计算等特点，但在计算过程中若需要改变节点或增加节点时，全部基函数 li(x) 都需要重新计算，这在实际计算中很不方便，为了克服这一缺点，介绍一种灵活增加节点的插值多项式。 本讲主要内容: ? Newton插值多项式的构造 ? 差商的定义及性质 ? 差分的定义及性质 ? 等距节点Newton插值公式 ?差商 (亦称均差 divided difference ) 定义 设实数 x0,x1, ?, xn 互异 , 一般情况用归纳法证之. 性质1 表明 ：差商具有对称性. 即, 差商的值与结点排列顺序无关. 牛顿插值多项式 Nn(x) 与 Lagrange 插值多项式 Ln(x) 都是过给定节点 x0,x1, ?, xn 的n 次多项式. 过 n ?1个点的 n 次多项式是唯一性的, 可知 Nn(x) ? Ln(x)， 从而, 它们的余项也相同，即 N2(x)=0.78846 + 0.43505(x ? 2.20) ? 0.087375(x ? 2.20) (x?2.40) N4(x)= 0.78846 + 0.43505(x ? 2.20） ? 0.087375(x ? 2.20)(x ? 2.40) +0.0225(x ? 2.20)(x ? 2.40)(x ? 2.60) ? 0.00755(x?2.20)(x?2.40)(x?2.60)(x?2.80) 练习：已知 f(0)=1，f(1)=2，f(2)=4，求 f(x)的二次插值多项式并求f(1.5). 插值函数的近似精度 利用插值函数pn(x)来近似计算被插函数 f(x)在区间[a,b]上任意点的函数值, 理论说明高次插值不可取。 二元函数插值： 设实值函数 f (x, y) 定义在矩形区域 D ={a ? x ? b, c ? y ? d }, 插值节点集: Z ={(xi, yj)| a ? x0 ? x1 ??? xn? b, c ? y0? y1??? ym? d }. 取在 Z 上线性无关的函数组 ??kr?x, y? | k=0,1, ?, n; r=0,1, ?, m?. 其中， ?kr?x, y? 是次数关于 x 不高于 n 次、 关于 y 不高于 m 次的二元多项式。 在函数空间 D ?Span??00, ?, ?0m, ?, ?n0, ?, ?nm ? 上寻找二元插值多项式 使其满足插值条件 此问题就是二元函数的代数插值问题。 如同一元的情况, 满足插值条件(5.11)的二元插值函数是唯一存在的。 Lagrange 插值曲面 取插值基函数 式(5.13)叫做 Lagrange 形式的插值曲面. 近似式 ?(x, y) ? pnm(x, y) (5.14) 叫做二元函数的Lagrange插值公式. 式 (5.14) 的余项 或 截断误差为 Rnm(x, y) ? ?(x, y) ? pnm(x, y) 例2 试利用f (x,y)的函数表 建立x为二次、y为一次的二元插值多项式p21(x,y)，用以计算f (0.3,0.8)的近似值。 数学符号 x0,x1, ?, xn, y0,y1, ?, ym, ??0?x?,?1?x?, ?, ?n?x?? k=0,1, ?, n ?(x, y) ? pnm(x, y) ??? ?? ? ?? ? ? ?????? ??? ??????? ΓΔΘΛΞΦΨΩ ????? ????? ????? ????? ????? ????? ????? ??????? ? ????? ??? ? ? ? ? ? ? ? ? ????? ? ? ? ? ? ?? ? ? ? ? ? ? ?? ?? ????????? ????? ????? ?????? ????????? ? ? ? ??????? ? ? ? ? ? ? ? ? ? ? ? ? ①②③④⑤⑥⑦⑧⑨ ?注： 实际使用时给定的函数表常常是等距节点的情形，这时只需考察函数值之差 。于是均差变成了差分，相应的Newton插值变成Newton前插与后插公式，当插值节点 a≤x0 , x1 , … , xn≤b 由小到大排列得到的是前插公式，反之，插值点由大到小排列得到的是后插公式，而利用插值计算f(x)的值时如果只用到函数表中的部分值。那么计算x0附近点