資源描述：
《牛頓插值法的MATLAB綜合程序.doc》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。
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給出節(jié)點數(shù)據(jù)，，，，作三階牛頓插值多項式，計算，并估計其誤差.解首先將名為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ùn)行后輸出插值y及其誤差限公式R，三階牛頓插值多項式P及其系數(shù)向量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
