数值分析20140930(广东工业大学课件).ppt

* * * * * * 2.5.2 一般埃尔米特（Hermite)插值多项式 2.三点三次带一个导数值的插值多项式 P34 2.5.3 利用重节点的差商构造埃尔米特（Hermite)插值多项式 书上P36 例2.6 问：这种方法和前面方法得到的插值多项式是否相同？ §2.6 分段低次插值法 一、高次插值的病态性质 龙格(Runge)现象…. 二、分段线性插值 所谓分段线性插值就是用通过插值点的折线段逼近f(x). 二、分段三次埃尔米特插值 分段线性插值函数导数间断，若已知节点上函数值和导数，可构造一个导数连续的插值函数Ih(x)，满足 * * * * * * * * * * * * * * * * * 第二章 插值 §3 牛顿插值 /* Newton’s Interpolation */ Lagrange 插值虽然易算，但若要增加一个节点时， 全部基函数 li(x) 都需重新算过。 将Ln (x) 改写成 的形式，希望每加一个节点时， 只附加一项上去即可。 差商(亦称均差) /* divided difference */ f (xi) - f (x j) xi - x j f[xi , x j]= 1阶差商 /* the 1st divided difference of f w.r.t. xi and xj */ f[xi , xj , xk]= 2阶差商 f[xi , xj] - f[xj , xk] xi - xk f[x0,x1,...,xk] - f[x1,...,xk,xk+1] x0 - xk+1 f[x0,...,xk - 1,xk] - f[x0,...,xk - 1,xk+1] xk - xk+1 f[x0,...,xk+1]= = §3 Newton’s Interpolation (k+1)阶差商： 注： k 阶差商必须由 k + 1个节点构成， 个节点是构造不 k 出 k 阶差商的。 为统一 起见,补充定义函数 f(x0)为零阶差商。 差商的值与 xi 的顺序无关 ！ n次多项式，则其 k阶差商 f [x 0 1L , x k -1, x] 2).差商具有对称性，即 f[x0,x1,L,xk]= f[x1,x0,x2L,xk]=L= f[x1,L,xk,x0] 4).若 f(x)是 , x §3 Newton’s Interpolation 差商性质： 1). k阶差商可表为函数值 f ( x 0 ), f ( x 1 ),L , f ( x k )的线性组合，即 由上式得差商表: k xk f(xk) 一阶差商 二阶差商 三阶差商 … 0 1 2 3 4 ┆ x0 x1 x2 x3 x4 ┆ f(x0) f(x1) f(x2) f(x3) f(x4) ┆ f[x0, x1] f[x1, x2] f[x0,x1,x2] f[x2, x3] f[x1,x2,x3] f[x0,x1,x2,x3] f[x3, x4] f[x2, x3,x4] f[x1,x2,x3,x4] … ┆ ┆ ┆ 规定函数值为零阶差商 作为Newton插值的以 x0，x1，L，xn为节点的基函数，而次数 其中 c0，c1，L，cn 是待定系数，由插值条件(2.1)决定。 Newton插值是通过选取特殊的基函数来实现的，这时，取 不超过n的多项式N n(x)可表示为 牛顿插值 /* Newton’s Interpolation */ §3 Newton’s Interpolation 因此就得到下列的满足插值条件(2.1)的 n次插值多项式 N n( x) = f ( x 0 ) + f [x 0, x1]( x -x 0 ) + f [x 0, x1, x 2]( x - x 0 )( x - x1) + L + f[x0, x1,L, xn](x - x0)(x - x1)L(x - xn - 1) §3 Newton’s Interpolation 通过插值条件运用数学归纳法可以求得 c k = f [ x 0 , x 1 ,L , x k ] §3 Newton’s Interpolation Nn(x)=c0+c1(x - x0)+c2(x - x0)(x -