实习题答案p50,p239(数值分析(第五版)).pdf

该【实习题答案p50,p239(数值分析(第五版)) 】是由【1781111****】上传分享，beplayapp体育下载一共【6】页，该beplayapp体育下载可以免费在线阅读，需要了解更多关于【实习题答案p50,p239(数值分析(第五版)) 】的内容，可以使用beplayapp体育下载的站内搜索功能，选择自己适合的beplayapp体育下载，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此beplayapp体育下载到您的设备，方便您编辑和打印。:..一1求牛顿插值多项式、差商、插值及其误差估计的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=;p=;q1=;c1=;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(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);,然后在MATLAB工作窗口输入程序>>symsM,X=[,,,,];Y=[,,,,];x=;[y,R,A,C,P]=newdscg(X,Y,x,M)运行后输出插值y=f()及其误差限公式R,三阶牛顿插值多项式P及其系数向量C,差商的矩阵A如下:y==(7*M)/800000A=----------=--:..P=-(25*x^4)/48+(5*x^3)/6-(53*x^2)/48+(23*x)/120+49/502:用拉格朗日插值法求其4次插值多项式的过程如下:先编写一个M文件来求其拉格朗日差值多项式,其程序如下:functionf=Language(x,y,x0)symstl;if(length(x)==length(y))n=length(x);elsedisp('x和y的维数不相等!');return;%检错endh=sym(0);for(i=1:n)l=sym(y(i));for(j=1:i-1)l=l*(t-x(j))/(x(i)-x(j));end;for(j=i+1:n)l=l*(t-x(j))/(x(i)-x(j));end;h=h+l;endsimplify(h);if(nargin==3)f=subs(h,'t',x0);%计算插值点的函数值?elsef=collect(h);f=vpa(f,6);%将插值多项式的系数化成6位精度的小数?,然后在MATLAB工作窗口输入程序:>>x=[];>>y=[];>>f=Language(x,y,)其运行结果为:f=>>f=Language(x,y)其运行结果为:f=:..-*t^4+*t^3-*t^2+*t+>>plot(x,y)::functionfun=dfun(x)%求原方程的导数f′(x)%UNTITLED4Summaryofthisfunctiongoeshere%Detailedexplanationgoesheresymstf=diff(fun0(t));fun=subs(f,t,x);end:..=fun0(x)%UNTITLED3Summaryofthisfunctiongoeshere%Detailedexplanationgoesherefun=x^2-3*x+2-exp(x);%[y,k]=newton1(x0,n,derta)%y为迭%代序列,k为迭代次数,x0初始位置,n为最大迭代次数,derta为%容许误差%UNTITLED5Summaryofthisfunctiongoeshere%Detailedexplanationgoesherek=1;y=[x0];t=x0-fun0(x0)./dfun(x0);whileabs(t-x0)>=dertay=[y,t];x0=t;k=k+1;t=x0-fun0(x0)./dfun(x0);if(k-1)>nerror('nisfull');:>>[y,k]=newton1(1,40,)y==:程序如下:f=inline('(x^2+2-exp(x))/3');x=feval(f,);disp(x);Eps=1E-8;i=1;while1:..x0=x;i=i+1;x=feval(f,x);disp(x);ifabs(x-x0)