《Matlab基础与应用教程答案(人民邮电出版社).docx》由会员分享，可在线阅读，更多相关《Matlab基础与应用教程答案(人民邮电出版社).docx（18页珍藏版）》请在课桌文档上搜索。
1、第1章一、思索题4. (1)B=A(2:5,1:2:5)(2) A(7)=(3) A=A+30(4) size（八）(5) t(find(t=O)=eps(6) t=reshape(x,3,4)(7) abs(matlab,)(8) char(93)54z7；9,3,0;5. A=97,67,34,10;-78,75,65,5;32,5,-23,-59;0,-12,(1) B=A(I:3,:)C=A(:,1:2)D=A(2:4,3:4)E=B*C(2)A(find(A=50&A120gz=dj*120+1.15*dj*(h-120);elseifh60gz=dj*h-700;elsegz=dj
2、*h;endformatbank;display(ghzgz)3.循环结构n=input(,inputn:);s=0;fork=l:ns=s+lk2;enddisplay(sqrt(s*6)向量运算n=input(inputn:);k=l:n;display(sqrt(sum(1./k.2)*6)4.y=0;k=0;whiley=le-5&k500x=x;x=a(b+x);k=k+l;enddisplay(krx);display(-b+sqrt(b2+4*a)/2,(-b-sqrt(b2+4*a)2);6.y=fun(40)/(fun(30)+fun(20)(1)函数文件fun.mfunct
3、ionf=fun(n)f=n+log(n2+5);(2)函数文件fun.mfunctionf=fun(n)a=l:n;f=sum(a.*(a+l);第3章一、思索题4.t=0:0.01:1.5;y=sqrt(3)2*exp(-4*t).*sin(4*sqrt(3)*t+pi3);5.x=-10:0.01:10;y=linspace(-6,6,size(xz2)z=x.A3+3*x.*y.a2;plot3(x,y,z)6.x=100:100:400;y=100:100:400;z=636z697z624z478;698z712z630z478;680z674,598z412;662r626z55
4、2z334;X,Y=meshgrid(xry);mesh(X,Yzz)二、试验题1.(1) x=-10:0.01:10;plot(xx-x.a332)(2) plot(xzexp(-x.*x2)2pi)(3) x=-8:0.01:8;plot(x,sqrt(64-x.*x)/2)(4) t=0:pi/100:6*pi;plot(t.*sin(t)zt.*cos(t)2. (1)Xl=Iinspace(0,1z100);yl=2*xl-0.5;t=linsace(0ziz100);x=sin(3*t).*cos(t);y=sin(3*t).*sin(t);plot(xlzyl,r-zx,yzb:
5、,);text(0.4zl,y=2x-0.5,);text(-0.4,0.5,x=sin(3t)cos(t),);text(-0.40.3z,y=sin(3t)sin(t),);(2)subplot(1,2,1);scatter(xlryl,10)title(,y=2x-0.5,);subplot(1,2,2);scatter(x,y,10)3.subplot(l,2zl);x=l:1:100;y=sin(1.x);plot(x,y)subplot(1,2z2);fplot(,sin(lx),zlz100)4.subplot(2,2zl);bar(t,y);subplot(2z2z2);sta
6、irs(tzy);subplot(2,2,3);stem(t,y)subplot(2,2,4);Semilogy(t,y);5.theta=linspace(-pizpi,100);ro=5.*cos(theta)+4;polar(thetazro);fi=linspace(0r2*pi,100);a=lr=a.*(l+cos(fi);polar(fizr);6.t=0:pi/20:2*pi;x=exp(-t./20).*cos(t);y=exp(-t./20).*sin(t);z=t;Plot3(,yfz);t=0:0.01:1;x=t;y=t.a2;z=t.a3;plot3(x,yrz);
7、7.x=-30:0.1:0;y=0:0.1:30;xzy=meshgrid(x,y);z=10.*sin(sqrt(x.2+y.a2).sqrt(l+x.A2+y.a2);meshe(,yfz);8.X=Iinspace(-3,3z100);y=linspace(-3z3z100);xy=meshgrid(x,y);fxy=-5./(1+x.A2+y.a2);i=find(abs(x)=0.8&abs(y)0.5)/size（八）2.h=466r715,950z1422z1635;w=7.04z4.28,3.40z2.52,2.13;hh=500z900z1500;ww=interpl(hzw
8、zhh,spline)3.X=Iinspace(1z10z50);y=log(x);f=polyfit(xzy,5);yy=polyval(fzx);Plot(X,y,r-,x,yy,g.)4.x=3*exp(-t) ; %求各采样点样本值XN=64;T=5;t=linsace(0zTzN);dt=t(2)-t(l);f=ldt;X=fft(x);F=X(l:N/2+l);f=f*(0:N/2)/N;plot(fzabs(F)-*)xlabel(Frequency);ylabel(F(k),)%采样点数%采样时间终点%给出N个采样时间ti(I=1:N)%采样周期%采样频率(HZ)%计算X的快速
9、傅立叶变换X%F(k)=X(k)(k=l:N/2+l)%使频率轴f从零起先%绘制振幅-频率图(1) pl=l2007;p2=l-2;p3=l0051;pl2=conv(plzp2);p=pl2+zeros(lzsize(pl2z2)-size(p3z2),p3;roots(p)(2) A=-1z4z3;2,1,5;0,5,6;Px=polyval(pzA)Pxm=polyvalm(p,A)6.(1)z=fzero(,3*x-sin(x)+1,0)(2)建立函数文件myfun.mfunctionF=myfun(X)x=X(l)；Y=X(2);F(l)=x*x+y*y-9;F(2)=x+y-l;在
10、叮嘱窗口中输入以下叮嘱：x=fsolve(Qmyfunz3z0,zOptimset(Display,off)第6章一、思索题2.fx=inline(l./(l+x.A2),);I,n=quad(fxz-100000z100OOO,le-3);Izn=quadl(fx,-100000z100000rle-3);x=-100000:0.01:100000;y=l./(l+x.*x);trapz(xzy);3.(1) fx=inline(,-2*y+2*x*x+2*x,);t,y=ode23(fx,0z0.5zl)(2) fx=inline(,y-exp(x)*cos(x),);t,y=ode23(
11、fxz03,1)二、试验题1. forx=l:3fx=xzxA2,xA3;l,2*x,3*x;0,2,6*x;diff(fx)end2.(1) x=0:0.01:l;y=x.a10+10.x+1./IoglO(x);dy=diff(y)/0.01;(2) x=0:0.01:l;y=log(l+x);dy=diff(y,2)0.01;plot(x(l:99)zdy)3.(1) fx=inline(,x.a2.*sqrt(2*x.*x+3),);quad(fxz1,5)(2) fx=inline(x.sin(x).a2);quad(fxzpi4zpi3)(3) fx=inline(,abs(cos
12、(x+y),);dblquad(fx,0zpiz0,pi)(4) symsxy;fx=x*y;int(int(fxzy2zy+2)z-lz2)X的积分区间为10,2时fx=inline(x.*y);dblquad(fxz0z2z-lz2)4.x=0.3:0.2:1.5;y=0.3895,0.6598z0.9147,1.1611z13971z1.6212z1.8325;trapz(xzy)5.(1)yp=inline(-(1.2+sin(10*x)*y,);tzy=ode23(ypz0z5zl);(2)令E=y,/=y,3=y”,则可写出原方程的状态方程形式：X；=%2M=X3Xy=cosr+5
13、cos27彳372(/+D23+sinr王，X2，X3010,0001X2+01-15cos2rCOS/3+sinr(/+D2L3J建立函数文件ztfun.mfunctiony=ztfun(tzx)b=0;0;cos(t);z-lz5*cos(2*t)(t+l)a2*+b;y=0,lz00,0,l-l(3+sin(t)解微分方程组t,y=ode23(ztfun,0,5z1;0;2);6.建立函数文件ztfun.mfunctionyy=ztfun(tzy)yy=y(2)*y(3);-y(l)*y(3);-0.51*y(1)*y(2);解微分方程组t,y=ode23(ztfun,0,5,0;1;1
14、)第7章一、思索题3.(1)数值积分fx=inline(exp(x).*(l+exp(x).a2,);quad(fx,0zlog(2)符号积分f=sym(,exp(x)*(l+exp(x)a2,);v=int(fz0zlog(2);eval(v)(2)略二、试验题A=sym(,l,2z3;xzy,z;3z2,l,)rank（八）inv（八）det（八）2.(1) y=sym(sqrt(x+sqrt(x+sqrt(x),);yl=diff(y)y2=diff(y,x,f2)(2) symsxy;fxy=sin(x2*y)*exp(-x2-y);diff(diff(fxy,x),y)3.(1) s
15、ymsXint(l(l+x4)(2) symsxtint(-2*x*x+l)/(2*x*x-3*x+l)2zxzcos(t),exp(2*t)4.symsnxsymsum(l/(2*n+l)/(2*x+l)(2*n+l),nz0zinf)symsum(l/(2*n+l)/(2*x+l)A(2*n+l)zn,0z5)5.(1) symsxtaylor(exp(x)+exp(-x)/2,5,0)(2) symsaxtaylor(exp(-5*x)*sin(3*x+pi3),5,a)(l)x=solve(sym(xA3+a*x+l=0)(2)xy=solve(sym(sqrt(xa2+ya2)-10
16、0=0z3*x+5*y-8=0,)7.方程转化为:y=yit2y+y=1符号解yl,yll=dsolvCDy=yl,Dyl+y=l-t*2pi,y(-2)=5,yl(-2)=5,t,)数值解编写函数文件ZtfUn.mfunctionyy=ztfun(tzy)yy=y(2)l-t2pi-y(l);在叮嘱窗口输入以下叮嗝t,y=ode45(ztfun,-2z7,-5;5);t=linspace(-2,7,49)V2=V8.xy=dsolve(,Dx=3*x+4*y,Dy=-4*x+3*y,x(0)=0zy(0)=1,)第9章二、试验题1.(1)新建一个BIankGUE(2)在窗口上添加一个Stat
17、iCText、一个EditText,一个PUShBotton,两个CheCkBox,一组RadiOButtono各控件的属性设置如表9.1所示。表9.1界面控件属性控件属性名属性值静态文本框String试验一:编辑框String按钮String设置字体大小单选按钮ISiringred单选按钮2Stringyellow单选按钮3Stringblue单选按钮4StringblackValue1.0复选框IString倾斜复选框2String加粗(3)在按钮的CaHbaCk函数中添加以下代码:s=eval(get(handles.editl,String,);set(handles.textl,Fo
18、ntSize,s);(4)在复选框1的Callback函数中添加以下代码：if(get(hbject,Value)=get(hbject,Max,)set(handles.textl,FontAngle,italic);elseset(handles.textl,FontAngle,normal);end在复选框2的Callback函数中添加以下代码：if(get(hbject,Value)=get(hbject,Max)set(handles.textl,FontWeight,bold*);elseset(handles.textl,FontAngle,normal,);end(5)在单选按
19、钮组的SeIeCtiOnChangeFCn函数中添加以下代码：set(handles.textl,ForegroundColor,get(eventdata.NewValue,string,);2.(1)新建一个BIankGUI。(2)用菜单编辑器编辑菜单，各菜单项属性如表9.2.表9.2菜单项属性菜单项Label属性值Tag属性值MenuLineLabelm_lineHelpHelpm_helpLineMenuItemLineStyleLineStylem-lincstylcLineWidthLineWidthminewidthLineColorLineColormlinecolorHelp
20、MenuItemAboutAboutm_about(3)绘图。在窗口的OPeningFCn函数添加以下代码:globalh;t=0:0.01:6;h=plot(t.*sin(t),t.*cos(t)(4)在LineSIyIe菜单项的Callback函数中添加以下代码:globalh;set(h,LineStyle在LineWidth菜单项的Callback函数中添加以下代码：globalh;set(h,LineWidth,2);在LineCOk)I菜单项的Callback函数中添加以下代码：globalh;set(h,color,red,);在About菜单项的Callback函数中添加以下代
21、码：msgbox(,HelpItem*);3.t=0:pi/20:6*pi;a=2;b=3;r=a+b*t;h_draw=polar(t,r,);h_menu=uicontextmenu;%Aoiv5Y2EP%-A2EYIiuimenu(h_menu,Label,red,Call,set(h_draw,Color,red,),);uimenu(h_menu,Label,yellow,Call,set(h_draw,Color,yellow*,)；,)；uimenu(h_menu,Label,blue,Call,set(h_draw,Color,blue,),);set(h_draw,UICon
22、textMenu,h_menu);%=0&u=2&u=5&u=,a&textail(n)=,ztextail(n)=char(textali(n)-(,a-,A,);endm=tetall(n)-,A,+1;ch(m)=ch(m)+1;endend4.(1)编辑C源程序exercisell_4.c,代码如下：#include#include#include#include#include*Forstrcm()*/*ForEXIT_FAILURE,EXIT_SUCCESS*/,mat.h,#defineBUFSIZE256#definekl39!definek271intmain()MATFil
23、e*pmat;mxArray*pal,*pa2,*pa3;doubledataklk2;constchar*file=,matll_4.mat;charstrBUFSIZE;intstatus;intnl,n2;for(nl=O;nlkl;nl+)for(n2=0;n2k2;n2+)datanln2=pow(nl+1.0)*(112+1),1/3.0);printf(nCreatingfile%s.nn,file);pmat=matOpen(file,w);if(pmat=NULL)printf(Errorcreatingfile%sn,file);printf(Doyouhavewritep
24、ermissioninthisdirectory?)n);return(EXIT_FAILURE);pal=InxCreateDoubleMatrix(kl,k2,mxREAL);if(pal=NULL)printf(,%s:Outofmemoryonline%dn,printf(,UnabletocreatemxArray.n);return(EXIT_FAILURE);pa2=mxCreateDoubleMatrix(kl,k2zmxREAL);if(pa2=NULL)printf(,%s:Outofmemoryonline%dnu,_FILE_,_LINE_);printf(,Unabl
25、etocreatemxArray.nu);return(EXIT_FAILURE);)memcpy(void*)(mxGetPr(pa2),(void*)data,sizeof(data);pa3=mxCreateString(MATLAB:thelanguageoftechnicalcomputing*);if(pa3=NULL)printf(,%s:Outofmemoryonline%dn,_FILE_,_LINE_);printf(,UnabletocreatestringmxArray.nn);return(EXIT_FAILURE);status=matPutVariable(pma
26、t,nLocalDouble,pal);if(status!=O)printf(,%s:ErrorusingmatPutVariableonline%dn/_FILE_,_LINE_);return(EXIT_FAILURE);status=InatputVariableAsGlobal(pmat,nGlobalDoublen,pa2);if(status!=O)printf(ErrorusingmatPutVariableAsGlobaln,);return(EXITFAILURE);status=InatPutVariable(pmat,Localstring,pa3);if(status
27、!=0)printf(s:ErrorusingmatPutVariableonline%dn*/FILE,_LINE_);return(EXIT_FAILURE);memcpy(void*)(mxGetPr(pal),(void*)data,sizeof(data);status=InatPutVariable(pmat,LocaIDouble,pal);if(status!=O)printf(,%s:ErrorusingmatPutVariableonline%dn,_FILE_,_LINE_);return(EXIT_FAILURE);mxDestroyArray(pal);InxDest
28、royArray(pa2);mxDestroyArray(pa3);if(matclose(pmat)!=O)printf(,Errorclosingfile%sn,file);return(EXIT_FAILURE);pmat=matpen(file,rn);if(pmat=NULL)printf(,Errorreopeningfile%sn,zfile);return(EXIT_FAILURE);pal=matGetVariable(pmat,LocalDoublen);if(pal=NULL)printf(,ErrorreadingexistingmatrixLocalDoublen,);return(EXIT_FAILURE);)if(InxGetNumberOfDimensions(pal)!=2)printf(Errorsavingmatrix:resultdoesnothavetwodi
