52牛顿插值.ppt

* 24 * 24 * 24 * 24 第二节、牛顿插值 Newton插值多项式的构造 差商的定义及性质 Newton插值方法 算法设计 Lagrange 插值虽然易算，但若要增加一个节点时，全部基函数 li(x) 都需要重新计算。 能否重新在Pn中寻找新的基函数 ？ 希望每加一个节点时，只附加一项上去即可。 能否构成Pn的一组基函数 利用插值条件，代入上式，得关于Ak的线性代数方程组 引例 设一个点　　　　唯一确实0次多项式为 代入可得 设二个点　　　　　　　多项式为 第一个点代入，后一项为零，即为一点式 第二个点代入，得 1、2个点的牛顿多项式 三个点： 前两个点代入得，第三项恒为零，即为二点式 第3个点代入得 3点牛顿多项式 General Form 差商(亦称均差) /* divided difference */ 0阶差商 1阶差商 2阶差商 K 阶差商 General Form where Rewriting 利用插值条件和差商,可求出Nn(x)的系数 bi 差商表 xk f ( xk ) 一阶差商 二阶差商 三阶差商 n 阶差商 差商表 例 x 1 2 3 4 5 6 y -3 0 15 48 105 192 根据下表数据，求牛顿插值多项式 xk f [xk] 一阶 二阶 三阶 四阶 五阶 1 2 0 3 15 4 48 5 105 6 192 3 15 33 57 87 6 9 12 15 1 1 1 0 0 0 解：构造差商表 -3 例1 :给定f(x)=lnx的数据表 xi 2.20 2.40 2.60 2.80 3.00 f(xi) 0.78846 0.87547 0.95551 1.02962 1.09861 解:差商表 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) Largran Code X=[1 2 3 4 5 6]; %X is a vector that contains a list of abscissas Y=X.^3-4*X; % Y is a vector that contains a list of ordinates n=length(X); D=zeros(n,n); D(:,1)=Y; %Use fozmula (20) to form the divided-difference table for j=2:n for k=j:n D(k,j)=(D(k,j-1)-D(k-1,j-1))/(X(k)-X(k-j+1)); end end %Determine the coefficients of the Newton interpolating polynomial C=D(n,n); for k=(n-1):-1:1 C=conv(C,poly(X(k))); m=length(C); C(m)=C(m)+D(k,k); end * 24 * 24 * 24 * 24