牛顿插值法的MATLAB综合程序.doc

《牛顿插值法的MATLAB综合程序.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、6.3.5牛顿插值法的MATLAB综合程序求牛顿插值多项式、差商、插值及其误差估计的MATLAB主程序function[y,R,A,C,L]=newdscg(X,Y,x,M)n=length(X);m=length(x);fort=1:mz=x(t);A=zeros(n,n);A(:,1)=Y';s=0.0;p=1.0;q1=1.0;c1=1.0;forj=2:nfori=j:nA(i,j)=(A(i,j-1)-A(i-1,j-1))/(X(i)-X(i-j+1));endq1=abs(q1*(z-X(j-1)));c1=c1*j;endC=A(

2、n,n);q1=abs(q1*(z-X(n)));fork=(n-1):-1:1C=conv(C,poly(X(k)));d=length(C);C(d)=C(d)+A(k,k);endy(k)=polyval(C,z);endR=M*q1/c1;L(k,:)=poly2sym(C);例6.3.6给出节点数据，，，，作三阶牛顿插值多项式，计算，并估计其误差.解首先将名为newdscg.m的程序保存为M文件，然后在MATLAB工作窗口输入程序>>symsM,X=[-4,0,1,2];Y=[27,1,2,17];x=-2.345;[y,R,A,C,P

3、]=newdscg(X,Y,x,M)运行后输出插值y及其误差限公式R，三阶牛顿插值多项式P及其系数向量C，差商的矩阵A如下y=22.3211R=65133/1312*M（即R=2.3503*M）A=27.00000001.0000-6.5000002.00001.00001.5000017.000015.00007.00000.9167C=0.91674.2500-4.16671.0000P=11/12*x^3+17/4*x^2-25/6*x+1