文档介绍：该【MATLAB基础及其应用教程-周开利-邓春晖课后答案 】是由【幸福人生】上传分享，文档一共【38】页，该文档可以免费在线阅读，需要了解更多关于【MATLAB基础及其应用教程-周开利-邓春晖课后答案 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。MATLAB基础及其应用教程-周开利-邓春晖课后答案第三章****题及参考答案解答:>>p=[1-1-1];>>roots(p)ans= -:取n=5,m=61>>x=linspace(0,2*pi,5);y=sin(x);>>xi=linspace(0,2*pi,61);>>y0=sin(xi);>>y1=interp1(x,y,xi);>>y2=interp1(x,y,xi,'spline');>>plot(xi,y0,'o',xi,y1,xi,y2,'-.');>>plot(xi,y0,'o',xi,y1,xi,y2,'-.');>>subplot(2,1,1);plot(xi,y1-y0);gridon>>subplot(2,1,2);plot(xi,y2-y0);gridon分段线性和三次样条插值方法与精确值之差解答:>>x=[0,300,600,1000,1500,2000];>>y=[,,,,,];>>xi=0:100:2000;>>y0=*exp(-(xi+500)/7756);>>y1=interp1(x,y,xi,'spline');>>p3=polyfit(x,y,3);>>y3=polyval(p3,xi);>>subplot(2,1,1);plot(xi,y0,'o',xi,y1,xi,y3,'-.');>>subplot(2,1,2);plot(xi,y1-y0,xi,y3-y0);gridon插值和拟合方法相比较,都合理,误差也相近。解答:梯形法积分>>x=-3::3;>>y=exp(-x.^2/2);>>z=trapz(x,y)/(2*pi)z=>>z=quad('exp(-x.^2/2)',-3,3)/(2*pi)z=-5~5:梯形法积分>>x=-5::5;>>y=exp(-x.^2/2);>>z=trapz(x,y)/(2*pi)z=>>z=quad('exp(-x.^2/2)',-5,5)/(2*pi)z=,两种积分的结果依然相同。梯形积分中改变x的维数为2维数组>>x(1,:)=-5::5>>x(2,:)=-5::5>>y=exp(-x.^2/2);>>z=trapz(x,y)/(2*pi)???Errorusing==>trapzLENGTH(X)mustequalthelengthofthefirstnon-。解答:>>x=linspace(0,1,4);>>y=x./(x.^2+4);>>t=cumsum(y)*(1-0)/(4-1);>>z1=t(end)>>z2=trapz(x,y)>>z3=quad('x./(x.^2+4)',0,1)>>z4=quadl('x./(x.^2+4)',0,1)z1====:>>A=[5121;2511;12102;12210];>>b=[991515]';>>tol=*10^-6;>>imax=5;>>x0=zeros(1,4);>>tx=jacobi(A,b,imax,x0,tol);>>forj=1:size(tx,1)fprintf('%4d%%%%\n',...j-1,tx(j,1),tx(j,2),tx(j,3),tx(j,4)):>>A=[5121;2511;12102;12210];>>b=[991515]';>>tol=*10^-6;>>imax=5;>>x0=zeros(1,4);>>w=;>>tx=sor(A,b,imax,x0,tol,w);