《MATLAB+7.X程序设计》习题答案汇总.doc 37页

'习题答案　　　　　　　37习题答案第1章1-2(1)使用Windows“开始”菜单。(2)运行MATLAB系统启动程序matlab.exe。(3)利用快捷方式。1-3MATLAB系统主要由开发环境、MATLAB函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。1-4窗口:M文件编辑/调试器、历史指令窗、当前目录浏览器、工作空间浏览器、内存数组编辑器、交互界面分类目录窗、及帮助导航/浏览器和图形窗。1-5在操作桌面上选择建立新文件或打开文件操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。MATLAB的开发环境中包括了专门的M文件编辑器，该编辑器不但提供了M文件的编辑功能，同时还与MATLAB的开发环境一起实现了MATLAB命令和函数文件的运行与调试。1-6Copy将所选历史命令拷贝到剪贴板EvaluateSelection运行所选历史命令CreateM-File打开编辑器，将所选历史命令复制到编辑器DeleteSelection删除所选历史命令(不确认是否删除)DeletetoSelection从头删除到所选历史命令(不确认是否删除)DeleteEntireHistory清除全部历史命令(将弹出对话框确认是否清除)1-7当用户在MATLAB命令窗口输入一条命令后，MATLAB按照一定次序寻找相关的文件。基本的搜索过程是：（1）检查该命令是不是一个变量。(2)检查该命令是不是一个内部函数。(3)检查该命令是否当前目录下的M文件。(4)检查该命令是否MATLAB搜索路径中其他目录下的M文件。用户可以将自己的工作目录列入MATLAB搜索路径，从而将用户目录纳入MATLAB系统统一管理。设置搜索路径的方法有：(1)用path命令设置搜索路径。例如，将用户目录c:mydir加到搜索路径下，可在命令窗口输入命令：37 习题答案　　　　　　　37path(path,’d:my’)(2)用对话框设置搜索路径1-8在MATLAB中有多种获得帮助的途径：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLABHelp菜单项可以打开帮助浏览器；（2）help命令：在命令窗口键入help命令可以列出帮助主题，键入help函数名可以得到指定函数的在线帮助信息；（3）lookfor命令：在命令窗口键入lookfor关键词可以搜索出一系列与给定关键词相关的命令和函数；注意：lookfor和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help命令显示详细信息。第2章2-1FORMATSHORT默认设置，以两位整数和4位小数形式输出FORMATLONG以16位十进制数形式输出FORMATSHORTE以5位十进制数加指数形式输出FORMATLONGE以16位十进制数加指数形式输出FORMATSHORTG从formatshort和formatshorte中自动选择最佳输出形式FORMATLONGG从formatlong和formatlonge中自动选择最佳输出形式FORMATHEX以16位十六进制数形式输出FORMAT+以正号、负号和零形式输出FORMATBANK以两位小数形式输出FORMATRAT以近似分数形式输出2-2x=[4+8i3+5i2-7i1+4i7-5i;3+2i7-6i9+4i3-9i4+4i];Rx=real(x)Ix=imag(x)Ax=abs(x)Cx=conj(x)CCx=x"2-3x=-5:2:5%建立行向量length(x)%向量x的长度x(1)%x的第一个元素37 习题答案　　　　　　　37x(end)%x的最后一个元素x"%向量转置为列向量y=linspace(-5,5,8)%建立行向量y"%向量转置为列向量2-4执行结果为：ans=33ans=2ans=9a12=1ah1=816al1=834ah1_3=816492al1_3=863742a=8160357049200009a=8160357037 习题答案　　　　　　　374920a=8163574922-5执行结果为：ans=6ans=151ans=-315ans=73-26-35ans=612-15ans=-666-121212-101010ans=147126-121410-194ans=72341012-9-3537 习题答案　　　　　　　37ans=381542-422954-38722ans=4919425369125ans=0.09760.4634-0.4146-0.17070.4390-0.02440.5732-0.90240.9390ans=0.14292.00000.33331.00000.40000.3333-1.0000-3.00000.20002-6执行结果为：ans=1-23ans=123ans=13-2-435ans=123345ans=5678ans=837 习题答案　　　　　　　37ans=6781011122-7已知角度，求x的正弦、余弦、正切和余切。在三角函数运算中，自变量可以是角度也可以是弧度。x=[304560];x1=x/180*pi;sin(x1),cos(x1),tan(x1),cot(x1)%自变量是弧度sind(x),cosd(x),tand(x),cotd(x)%自变量是角度2-8x=10:99;length(x((rem(10:99,5)==0)))2-9x=10:20;log10(x)/log10(3)2-10a=[1.55-1.451.45-1.55];floor(a),ceil(a),round(a),fix(a)执行结果为：ans=1-21-2ans=2-12-1ans=2-11-2ans=1-11-12-11s1=upper("China");s2="china";s3=lower("Chinese");c1=strcmp(s1,s2)c2=strcmpi(s1,s2)c3=strncmp(s1,s3,4)c4=strncmpi(s1,s3,2)执行结果为：c1=0c2=137 习题答案　　　　　　　37c3=0c4=12-12执行结果为：ans=5555ans=-3-113ans=20ans=4664ans=1894ans=14916ans=16842ans=24816168422-13执行结果为：ans=2h=2l=3ans=6ans=[2x2double]ans=37 习题答案　　　　　　　37[1x27char][2x2double][2x2double]ans=[2x2double][]ans=[2x2double][2x2double][3x3double][]ans=0-310ans=0ans=-3ans=01c=[1x27char][2x2double][2x2double][2x2double][3x3double][][][][5]c=[1x27char][2x2double][2x2double][2x2double][3x3double][]2-14t.college="computer"t.name="WangMing";t.number="19990101";t.sex="man";t.corse.name="MATLAB";t.corse.class="2004";t.corse(2).name="maple";t.corse(2).class="2005";t2.college="computer"t2.name="zhangMing";t2.number="19990101";37 习题答案　　　　　　　37t2.sex="man";t2.corse.name="math";t2.corse.class="2006";t2.corse(2).name="MATLAB";t2.corse(2).class="2005";t2.corse(3).name="math";t2.corse(3).class="2006";t第3章3-1a=[158;30-4];b=[6-73;502];a>ba>=ba6;%产生与A同型的0—1逻辑值数组x=find(a>6);a(L)%把L中逻辑值1对应的A元素取出x%大于6的元素在原数组中的序号3-4执行结果为：x=000000011y=111100000L=111100011aL=1.00002.10003.20004.30008.70009.8000xh=1234893-5a=reshape(-6:2:12,2,5);b=reshape(-8:3:19,2,5);la=a&b;lb=a|b;islogical_la=islogical(la)islogical_lb=islogical(lb)any_lab=any([la,lb])all_lab=all([la,lb])37 习题答案　　　　　　　37执行结果为：islogical_la=1islogical_lb=1any_lab=1111111111all_lab=10111111113-6a=magic(3);b=magic(4);P1=rem(a,3)==0;P2=rem(b,3)==0;a(P1)b(P2)3-7执行结果为：all_a=10111any_a=11111isnan_a=0000000000isinf_a=0000000000isfinite_a=11111111113-8a=input("Pleaseinputamatrix?");ifisreal(a)==1sign(a)else37 习题答案　　　　　　　37disp("Thematrixisnotrealmatrix")end3-9x=input("请输入x的值:");ifisreal(x)==0disp("Nodefination")elseifx>=-1&x<=1y=x.^2+1elsey=x+1end程序2：x=input("请输入x的值:");isreal_x=isreal(x);switchisreal_xcase1ifx>1|x<-1y=1+xelsey=x.^2+1endcase0disp("Nodefination")end3-10x=input("pleaseenterxvalue:");ifx<0&x~=-3y=x.^2+x-6elseifx>=0&x~=2&x~=3&x<10y=x.^2-5*x+6elsey=x.^2-x-1end3-11执行结果：37 习题答案　　　　　　　37a=2-10000-12-10000-12-10000-12-10000-12-13-12s=1083-13程序1：y=0;forn=1:100f=1/n/n/n;y=y+f;n=n+1;endy程序2：y=0;n=1;whilen<=100f=1/n/n/n;y=y+f;n=n+1;endy3-14n=input("pleaseenternvalue:");s=0;fori=1:2:ns=s+i;ends3-15程序：a=0;b=3;n=1000;h=(b-a)/n;37 习题答案　　　　　　　37x=a:h:b;f=x./(x.^2-1);fori=1:ns(i)=(f(i)+f(i+1))*h/2;ends=sum(s)3-16程序1：n=input("Pleaseentern:");y1=0;fori=1:n;y1=y1+1/i/i;y=sqrt(6*(y1));endy程序2：n=input("Pleaseenterthenumbern:");i=1:n;f=1./i.^2;y1=sum(f);y=sqrt(6*(y1))3-17m=input("矩阵行数：m=");n=input("矩阵列数：n=");fori=1:mforj=1:ndisp(["输入第",num2str(i),"行，第",num2str(j),"列元素"])A(i,j)=input(‘’)endend3-18程序1：functiony=f2(x)y=x.^3+2*x.^2+1;37 习题答案　　　　　　　37y=f2([-5-3-1135])程序2：f=inline("x.^3+2*x.^2+1");f([-5-3-1135])3-19formonth=1:12;switchmonthcase{3,4,5}season="spring"case{6,7,8}season="summer"case{9,10,11}season="autumn"otherwiseseason="winter"endend3-203-21定义函数文件：functionz=ff(x,y)z=x^2+exp(x+y)-y*log(x)-3;调用：z=ff(2,4)第4章4-1结果为：Incotdat16y=00.50000.86601.00000.86600.500001.00000.86600.50000-0.5000-0.8660-1.000000.57741.7321Inf-1.7321-0.57740Inf1.73210.57740-0.5774-1.7321Infy_2_max=37 习题答案　　　　　　　3711InfInfy_max=Infy_2_min=0-1.0000-1.7321-1.7321y_min=-1.73214-2cx=cov(x),cy=cov(y)cxy=cov(x,y)4-3执行结果为：ans=9-3630-36192-18030-180180ans=9192180ans=30-180180-36192-1809-3630ans=30-369-180192-36180-18030ans=30-18018037 习题答案　　　　　　　37-36192-1809-3630ans=9-3630-36192-18030-180180ans=900-36192030-180180ans=9-36300192-180001804-4定义函数文件：functionx4_4(a)%amatrix[mn]=size(a);ifm==ndet_a=det(a)ifdet_a~=0inv_a=inv(a)endelsedisp("noinverse")pinv_a=pinv(a)disp("nodet")endrank_a=rank(a)>>x4_4([1-632;3-540;-1-1124])noinversepinv_a=-0.10000.1860-0.03000.5000-0.2500-0.25000.7000-0.2020-0.29001.0000-0.5400-0.300037 习题答案　　　　　　　37nodetrank_a=3>>x4_4([1210;6241;0210;3141])det_a=1inv_a=1.00000-1.00000-3.00001.00003.0000-1.00006.0000-2.0000-5.00002.0000-24.00007.000020.0000-6.0000rank_a=44-5a=[1132;-11-13;5-289;-1317];rank(a)ans=3线性相关4-6a=[22-1;1-24;582];[X,D]=eig(a)X=0.0000-0.71600.2760-0.44720.4200-0.7472-0.89440.55760.6045D=6.00000001.6056000-5.60564-7A=[65-25;9-14-1;342-2;3-920];b=[-413111]";inv(A)*bA=[65-25;9-14-1;342-2;3-920];b=[-413111]";Abans=-8.33332.000027.000037 习题答案　　　　　　　3718.0000A=[65-25;9-14-1;342-2;3-920];b=[-413111]";Ab=[A,b];rref(Ab)ans=1.0000000-8.333301.0000002.0000001.0000027.00000001.000018.00004-9x=[113;432;125]*inv([11-1;210;111])x=-102012-1-144-10p=[10-2-5];r=roots(p);forii=1:length(r)ifr(ii)>=0&r(ii)<=3r(ii)endendans=2.09464-11（1）、的根；p1=[3-12103];p2=[1/31-3-1];roots(p1),roots(p2)ans=0.7330+0.7416i0.7330-0.7416i-0.8952-0.1188+1.0066i-0.1188-1.0066ians=-4.75882.064237 习题答案　　　　　　　37-0.3054（2），和；p1+[00p2]ans=3.0000-1.00002.33332.0000-3.00002.0000p1p1=3-12103conv(p1,p2)ans=1.00002.6667-9.33332.3333-4.0000-4.00002.0000-9.0000-3.0000[rpk]=residue(p1,p2)r=-792.245022.1035-0.8584p=-4.75882.0642-0.3054k=9-30177（3）、、的导数。polyder(p1)ans=15-4620polyder(p2)ans=12-3[PQ]=polyder(p1,p2)P=2.00008.6667-38.0000-4.3333-8.0000-12.0000-8.00009.0000Q=37 习题答案　　　　　　　370.11110.6667-1.0000-6.66677.00006.00004-12（1）时求的值。P=[321];polyval(p,[579])ans=-108.9538-219.0359-367.1882（2）时，求的值。P=[321];polyvalm(p,[25;79])ans=-181.7691-251.4115-351.9761-533.74524-13p=poly([638])p=1-1790-144roots(p)ans=8.00006.00003.00004-14poly([-0.5-0.3+0.4*i-0.3-0.4*i])ans=1.00001.10000.55000.12504-15p=[4-12-1459];diff(p)ans=-16-2194p=[4-12-1459];polyint(p)ans=0.8000-3.0000-4.66672.50009.000004-1637 习题答案　　　　　　　37x=0:0.01:1;y=sin(x).*exp(-x/10);interp1(x,y,0.25)ans=0.24134-17x=1:10;y=log(x)/log(3);interp1(x,y,3.25)ans=1.06554-18x=[-122.534];y=[0.233.53.62.5];f=polyfit(x,y,2)polyval(f,3.5)4-19a=[111];polyint(a)quad("x.*log(x)",1,exp(1))ans=2.0973第5章5-1（1）f=3*x^2+5*x+2表示在给定x时，将3*x^2+5*x+2的数值运算结果赋值给变量f，如果没有给定x则指示错误信息。（2）f="3*x^2+5*x+2"表示将字符串"3*x^2+5*x+2"赋值给字符变量f，没有任何计算含义，因此也不对字符串中的内容做任何分析。（3）x=sym("x")f=3*x^2+5*x+2表示x是一个符号变量，f=3*x^2+5*x+2就具有了符号函数的意义，f也自然成为符号变量了。5-2symsphi1phi2;y=simple(sin(phi1)*cos(phi2)-cos(phi1)*sin(phi2))…37 习题答案　　　　　　　375-3symsx;symsy;f=x^2+exp(x+y)-y*log(x)-3;subs(f,{x,y},{2,4})5-4symsa11a12a21a22;A=[a11,a12,a13;a21,a22,a23;a31,a32,a33]AD=det(A)%行列式AI=inv(A)%逆AE=eig(A)%特征值5-5（1）求f,g的反函数和复合函数，（2）分别求f,g的微分，（3）分别求f,g的积分。symsxyt;f=sym(sin(x));g=sym(y/exp(-2*t));finverse(f)finverse(g)finverse(g,"t")compose(f,g)compose(f,g,"z")diff(f),diff(g,"t")int(f)int(g)int(g,"t")int(g,"t",0,10)5-6symsx;f=x^4-5*x^3+5*x^2+5*x-6;factor(f)horner(f)5-7symsx;f=(x+1)^3+(x-1)^2+5*x-6collect(f,x)37 习题答案　　　　　　　375-8symsx;expand((x+1)^6)5-9syms("x","y");[n,d]=numden(1/(x^3-1)+1/(x^2+y+1)+1/(x+y+1)+8)5-10求下列极限。(1)x=sym("x");limit(1/(x+1)-3/(x^3+1),x,-1)(2)x=sym("x");limit((1+cos(x))^(3*sec(x)),x,pi/2)(3)x=sym("x");limit(log(sin(x))/((pi-2*x)^2),x,pi/2)(4)x=sym("x");limit(x^2*exp(1/x/x),x,0)(5)x=sym("x");limit(((5*x^2)/(1-x^2)+2^(1/x)),x,inf)(6)x=sym("x");limit((x^2-2*x+1)/(x^3-x),x,1)5-11求下列幂级数的和函数。（1）n=sym("n");symsum((-1)^n*(x^n)/sqrt(n*2-n),"n",2,inf)（2）n=sym("n");symsum((-1)^n*(sin(pi/(2^n)*(x^n))),"n",1,inf)（3）n=sym("n");symsum(x^n/n,"n",1,inf)（4）n=sym("n");symsum(n*x^n,n,1,inf)（5）n=sym("n");symsum((-1)^n*(x^n)/sqrt(n^2-n),"n",1,inf)（6）n=sym("n");symsum((-1)^n*sin(pi*(x^n)/(2^n)),"n",1,inf)5-12求下列函数的导数。(1)x=sym("x");diff(x^10+10^x+log(10)/log(x))(2)x=sym("x");subs(diff((1-sqrt(x))/(1+sqrt(x))),4)（3）x=sym("x");a=diff([sqrt(x^2-2*x+5)，cos(x^2)+2*cos(2*x)，4^(sin(x))，log(log(x))])5-13(1)symst;dx_dt=diff(t^4,t);dy_dt=diff(4*t,t);dy_dx=dy_dt/dx_dt(2)37 习题答案　　　　　　　37symst;dx_dt=diff(log(1+t^2),t);dy_dt=diff(t-atan(t),t);dy_dx=dy_dt/dx_dt5-14symsx;dy_dx=diff(exp(x)*cos(x),x,4)5-15symsx;symsy;d2y_d2x=diff(exp(x)*sin(x),x,2);dy_dx=diff(exp(x)*sin(x));f=d2y_d2x-2*dy_dx+2*y5-16%symsxy;diff(x*log(x+y),x,2)%symsxy;diff(x*log(x+y),y,2)%symsxy;diff(diff(x*log(x+y),x),y)5-17symsax;f=[a,x^2,1/x;exp(a*x),log(x),sin(x)];df=diff(f)5-181~3symsx;symsareal-y=[sqrt(a^2-x^2),(x-1)/(3*x-1)^(1/3),x^2*asin(x)];int(y,x)4symsx;int(x^2/(x+1),x)37 习题答案　　　　　　　375symsx;int(sin(2*x)/sqrt(1+sin(x)^2),x)6symsx;int(1/sqrt(x^2+5),x)5-19计算下列定积分。(1)symsx;int(x*log(x),1,2.18)(2)symsx;int(x/sin(x)/sin(x),pi/4,pi/3)(3)symsx;int(sin(log(x)),-1,2.718)(4)symsx;int(x^3*(sin(x)^2)/(x^4+x^2+1),-1,1)（5）symsx;int(sqrt(sin(x)^3-sin(x)^5),0,pi)(6)symsx;int(exp(-x*x/2),0,1)5-20解下列方程：（1）symsxsolve(x^5-5*x-1)（2）symsxsolve(x*2^x-1)（3）symsxsolve(x-3*sin(x)-1)（4）symsxsymspqsolve(x^3+p*x+q,"x")5-21（1）x=sym("x");y=sym("y");z=sym("z");r=sym("r");s=solve("6*x+5*y-2*z+5*r+4","9*x-y+4*z-r-13","3*x+4*y+2*z-2*r-1","3*x-9*y+2*z-11")s.xs.ys.zs.r(2)x=sym("x");y=sym("y");z=sym("z");r=sym("r");37 习题答案　　　　　　　37s=solve("2*x+2*y-z+r-4","4*x+3*y-z+2*r-6","8*x+3*y-3*z+4*r-12","3*x+3*y-2*z-2*r-6")s.xs.ys.zs.r(3)A=sym([1-11-1;1-1-11;1-1-22]);b=sym([1;0;-0.5]);X=Ab%求一个特解：最少非零元素的最小二乘解symsk;null_A=null(A)[mn]=size(null(A));forii=1:nXx=X+k*null_A(:,ii)%构成通解end（4）x=sym("x");y=sym("y");z=sym("z");r=sym("r");s=solve("x-y+4*z-2*r","x-y-z+2*r","3*x+y+7*z-2*r","x-3*y-12*z+6*r")s.xs.ys.zs.r5-221dsolve("D2y+3*Dy+exp(x)=0","x")2dsolve("D2y-exp(2*y)*Dy=0","x")5-231dsolve("(x^2-1)*Dy+2*x*y-cos(x)=0","y(0)=1","x")2.dsolve("x^2+2*x*y-y^2+(y^2+2*x*y-x^2)*Dy=0","y(1)=1","x")3.symsn;symsV0dsolve("D2x+2*n*Dx+a^2*x=0","x(0)=x0","Dx(0)=V0","t")37 习题答案　　　　　　　375-24（1）。taylor(x^5-x^3+2*x+1,-1,7)(2)taylor(x*exp(x),0,6)第6章6-1x=-2:0.01:2;plot(x.^2.*sin(x.^2-x-2))fplot("x.^2*sin(x.^2-x-2)",[-22])ezplot("x.^2*sin(x.^2-x-2)",-2,2)6-2subplot(2,2,1);ezplot("cos(x)",[0,2*pi])subplot(2,2,2);ezplot("sin(x-pi/2)",[0,2*pi])subplot(2,2,3);ezplot("(x^2)*cos(x-pi)",[0,2*pi])subplot(2,2,4);ezplot("exp(sin(x))",[0,2*pi])6-3(1)symsx;symsy;y=1+log(x+2);g=finverse(y);x=-2:0.001:2;y=subs(y,x);plot(x,y)holdon;g=subs(g,y);plot(y,g)(2)symsx;symsy;y=(x^2+1)^(1/3);g=finverse(y);x=-2:0.001:2;y=subs(y,x);plot(x,y)holdon;g=subs(g,y);37 习题答案　　　　　　　37plot(y,g)6-4(1)symsx;f=1/(x+1)-3/(x^3+1);y=limit(f,x,-1)ezplot(f);holdon;vy=double(y)plot(-1,vy,"r.")(2)symsx;f=(1+2*x)^(1/x);y=limit(f,x,0);ezplot(f);holdon;vy=double(y);plot(0,vy,"r.")6-5x=-10:0.01:10;y=2*exp(-0.5*x);%y=2*a^(-0.5*x)%y=abs(sin(x))%y=sin(abs(x))%y=2*sin(3*x)plot(x,y)figuresubplot(2,2,1);bar(x,y,"g")subplot(2,2,2);fill(x,y,"r")subplot(2,2,3);stairs(x,y,"b")subplot(2,2,4);stem(x,y,"k")6-6symsx;f=x^2+x+1;g=log(x);37 习题答案　　　　　　　37fg=compose(f,g)gf=compose(g,f)x=[10:0.1:20];plot(x,subs(fg,x),x,subs(gf,x))6-7x=0:0.1:10;y=2*x.^3+10*x.^2+5*x+6;figure(1)subplot(2,2,1)plot(x,y);title("LinearPlot");xlabel("x");ylabel("y");gridon;subplot(2,2,2)semilogx(x,y);title("SemilogxPlot");xlabel("x");ylabel("y");gridon;subplot(2,2,3)semilogy(x,y);title("SemilogyPlot");xlabel("x");ylabel("y");gridon;subplot(2,2,4)loglog(x,y);title("LoglogPlot");xlabel("x");ylabel("y");gridon;6-8（1）t=-1:0.001:1;plot(t.^4,4*t)（2）t=-10:0.001:10;plot(log(1+t.^2),t-atan(t))37 习题答案　　　　　　　376-9symsx;symsy;f=y^2-x;g=y-x^2;[x,y]=solve(f,g,x,y);x=x(1):0.001:x(2);plot(x,sqrt(x),x,x.^2)s=int(sqrt(x)-x.^2,x(1),x(2))6-10%只画实部t=0:pi/20:4*pi;y=log(t.^2+1).*cos(t)+i*sin(2*t);subplot(2,2,1),plot(t,y);title("plotofcomplexfunctionvstime");xlabel("t");ylabel("y(t)");%实部与虚部分别画t=0:pi/20:4*pi;y=log(t.^2+1).*cos(t)+i*sin(2*t);subplot(2,2,2),,plot(t,real(y),"b-");holdon;plot(t,imag(y),"r--");title("plotofcomplexfunctionvstime");xlabel("t");ylabel("y(t)")legend("real","imaginary");holdoff;%同时画实部与虚部t=0:pi/20:4*pi;y=log(t.^2+1).*cos(t)+i*sin(2*t);subplot(2,2,3),plot(y,"b-");title("plotofcomplexfunction");xlabel("frealpart");ylabel("imagainarypart");%极坐标t=0:pi/20:4*pi;37 习题答案　　　　　　　37y=log(t.^2+1).*cos(t)+i*sin(2*t);subplot(2,2,4),,polar(angle(y),abs(y));title("plotofcomplexfunction");6-11theta=-2*pi:0.01:2*pi;polar(theta,sin(theta))6-12t=-2*pi:0.01:2*pi;plot3(exp(-0.2*t).*cos(2*t),exp(-0.2*t).*sin(2*t),t)6-13%三维曲线x=-1:0.01:1;y=x;z=sin(pi*sqrt(x.^2+y.^2));plot3(x,y,z)%三维网格x=-1:0.01:1;y=x;[x,y]=meshgrid(x,y);z=sin(pi*sqrt(x.^2+y.^2));mesh(x,y,z)%三维网面x=-1:0.01:1;y=x;[x,y]=meshgrid(x,y);z=sin(pi*sqrt(x.^2+y.^2));surf(x,y,z)6-14%三维曲线u=0:0.01:2*pi;v=0:0.01:2*pi;x=(1+cos(u)).*cos(v);y=(1+cos(u)).*sin(v);z=sin(u);plot3(x,y,z)第7章7-1Simulink的基本模块包括16个子模块库。7-237 习题答案　　　　　　　377-37-47-57-6方法1：37 习题答案　　　　　　　37方法2：37 习题答案　　　　　　　377-77-87-937 习题答案　　　　　　　37第8章8-1菜单项组合键功能DefineInputCellAlt+I定义输入单元DefineAutoInitCellAlt+A定义自动初始化单元DefineCalcZoneAlt+Z定义计算区UndefineCellsAlt+U将单元转换为文本PurgeOutputCellsAlt+P清除输出单元GroupCellsAlt+G定义单元组UngroupCellsAlt+p将单元组转换为单个单元Hide/ShowCellMarkersAlt+C隐藏/显示单元标志ToggleGraphOutputforcell为每个单元锁定图形输出EvaluateCellCtrl+Enter运行当前单元或单元组EvaluateCalcZoneAlt+Enter运行当前计算区EvaluateM-bookAlt+M运行M-book中所有单元EvaluateLoopAlt+L循环运行单元BringMATLABtofront将MATLAB置于屏幕之前NotebookOptions…Alt+O定义输出显示选项8-2定义输入单元的方法是：首先选中所需命令，然后在Notebook菜单项中选择DefineInputCell命令，于是被选中的MATLAB命令成为输入单元。定义输入单元也可以在选中所需命令后，直接按组合键Alt+I。为了执行输入单元，应选择Notebook菜单项中的EvaluateCell命令或直接按组合键Ctrl+Enter。8-3把已有的多个独立输入单元或自初始化单元同时选中，然后在Notebook菜单中选择GroupCells，于是，便获得以第一个独立单元的性质组合而成的单元组。37 习题答案　　　　　　　378-4输入单元执行后产生输出单元，输出单元包含MATLAB的输出结果，包括数据、图形和出错信息。如果输入单元经修改后重新执行，那么新的输出单元将替换原有的输出单元。输出单元的输出格式可以重新设置，设置的方法为选择菜单Notebook→NotebookOptions…，则出现设置对话框。8-5方法是：先选定包含普通Word文本、输入单元和输出单元的一个连续区，然后选择Notebook菜单中的DefineCalcZone命令。一旦计算区被定义后，不管光标在计算区的什么位置，只要选择Notebook菜单中的EvaluateCalcZone命令即可执行计算区中的全部输入单元，且在每个输入单元后面以输出单元形式给出相应的计算结果。8-61单元的循环运行Notebook提供了循环运行单元的命令，先选定需要循环运行的输入单元，然后选择菜单Notebook→EvaluateLoop，就会出现设置对话框。2整个M-book文件的运行菜单Notebook→EvaluateM-book是运行整个M-book文件，即把文档中所有输入单元送到MATLAB中去运行。37'