《《MATLAB+7.X程序设计》习题答案汇总.doc》由会员分享，可在线阅读，更多相关《《MATLAB+7.X程序设计》习题答案汇总.doc（37页珍藏版）》请在三一文库上搜索。
1、 习题答案 37 习题答案第1章 1-2 (1) 使用Windows“开始”菜单。(2) 运行MATLAB系统启动程序matlab.exe。(3) 利用快捷方式。1-3 MATLAB系统主要由开发环境、MATLAB函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。1-4 窗口:M文件编辑/调试器、历史指令窗、当前目录浏览器、工作空间浏览器、内存数组编辑器、交互界面分类目录窗、及帮助导航/浏览器和图形窗。1-5 在操作桌面上选择建立新文件或打开文件操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。MATLAB的开发环境中包括了专门的M文件
2、编辑器，该编辑器不但提供了M文件的编辑功能，同时还与MATLAB的开发环境一起实现了MATLAB命令和函数文件的运行与调试。1-6 Copy将所选历史命令拷贝到剪贴板EvaluateSelection运行所选历史命令CreateM-File打开编辑器，将所选历史命令复制到编辑器DeleteSelection删除所选历史命令(不确认是否删除)DeletetoSelection从头删除到所选历史命令(不确认是否删除)DeleteEntireHistory清除全部历史命令(将弹出对话框确认是否清除)1-7 当用户在MATLAB命令窗口输入一条命令后，MATLAB按照一定次序寻找相关的文件。基本的搜索
3、过程是：（1）检查该命令是不是一个变量。(2)检查该命令是不是一个内部函数。(3)检查该命令是否当前目录下的M文件。(4)检查该命令是否MATLAB搜索路径中其他目录下的M文件。用户可以将自己的工作目录列入MATLAB搜索路径，从而将用户目录纳入MATLAB系统统一管理。设置搜索路径的方法有：(1) 用path命令设置搜索路径。例如，将用户目录c:mydir加到搜索路径下，可在命令窗口输入命令： path(path,d:my)(2) 用对话框设置搜索路径1-8 在MATLAB中有多种获得帮助的途径：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help
4、菜单项可以打开帮助浏览器；（2）help命令：在命令窗口键入help 命令可以列出帮助主题，键入help 函数名可以得到指定函数的在线帮助信息；（3）lookfor命令：在命令窗口键入lookfor 关键词可以搜索出一系列与给定关键词相关的命令和函数；注意：lookfor和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help命令显示详细信息。第2章 2-1 FORMAT SHORT 默认设置，以两位整数和4位小数形式输出FORMAT LONG 以16 位十进制数形式输出FORMAT SHORT E以5位十进制数加指数形式输出FORMAT LONG E 以16位十进制数加指数形
5、式输出FORMAT SHORT G 从format short和format short e中自动选择最佳输出形式FORMAT LONG G 从format long和format long e中自动选择最佳输出形式FORMAT HEX 以16位十六进制数形式输出FORMAT + 以正号、负号和零形式输出FORMAT BANK 以两位小数形式输出FORMAT RAT 以近似分数形式输出2-2 x=4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i; Rx=real(x)Ix=imag(x)Ax=abs(x)Cx=conj(x)CCx=x2-3 x
6、=-5:2:5 %建立行向量length(x) %向量x的长度x(1) %x的第一个元素x(end) % x的最后一个元素x % 向量转置为列向量y=linspace(-5,5,8) %建立行向量y % 向量转置为列向量2-4 执行结果为：ans = 3 3ans = 2ans = 9a12 = 1ah1 = 8 1 6al1 = 8 3 4ah1_3 = 8 1 6 4 9 2al1_3 = 8 6 3 7 4 2a = 8 1 6 0 3 5 7 0 4 9 2 0 0 0 0 9a = 8 1 6 0 3 5 7 0 4 9 2 0a = 8 1 6 3 5 7 4 9 22-5 执行结
7、果为：ans = 6ans = 1 5 1ans = -3 1 5ans = 7 3 -2 6 -3 5ans = 6 12 -15ans = -6 6 6 -12 12 12 -10 10 10ans = 14 7 12 6 -12 14 10 -19 4ans = 7 2 3 4 10 12 -9 -3 5ans = 38 15 42 -42 29 54 -38 7 22ans = 49 1 9 4 25 36 9 1 25ans = 0.0976 0.4634 -0.4146 -0.1707 0.4390 -0.0244 0.5732 -0.9024 0.9390ans = 0.1429
8、 2.0000 0.3333 1.0000 0.4000 0.3333 -1.0000 -3.0000 0.20002-6执行结果为：ans = 1 -2 3ans = 1 2 3ans = 1 3 -2 -4 3 5ans = 1 2 3 3 4 5ans = 5 6 7 8ans = 8ans = 6 7 8 10 11 122-7已知角度，求x的正弦、余弦、正切和余切。在三角函数运算中，自变量可以是角度也可以是弧度。x=30 45 60; x1=x/180*pi;sin(x1),cos(x1),tan(x1),cot(x1)% 自变量是弧度sind(x),cosd(x),tand(x),
9、cotd(x) %自变量是角度2-8 x=10:99;length(x(rem(10:99,5)=0)2-9 x=10:20;log10(x)/log10(3)2-10 a=1.55 -1.45 1.45 -1.55;floor(a),ceil(a),round(a),fix(a)执行结果为：ans = 1 -2 1 -2ans = 2 -1 2 -1ans = 2 -1 1 -2ans = 1 -1 1 -12-11 s1=upper(China);s2=china;s3=lower(Chinese);c1=strcmp(s1,s2)c2=strcmpi(s1,s2)c3=strncmp(s
10、1,s3,4)c4=strncmpi(s1,s3,2)执行结果为：c1 = 0c2 = 1c3 = 0c4 = 12-12 执行结果为：ans = 5 5 5 5ans = -3 -1 1 3ans = 20ans = 4 6 6 4ans = 1 8 9 4ans = 1 4 9 16ans = 16 8 4 2ans = 2 4 8 16 16 8 4 22-13 执行结果为：ans = 2h = 2l = 3ans = 6ans = 2x2 doubleans = 1x27 char 2x2 double 2x2 doubleans = 2x2 double ans = 2x2 doub
11、le 2x2 double 3x3 double ans = 0 -3 1 0ans = 0ans = -3ans = 0 1c = 1x27 char 2x2 double 2x2 double 2x2 double 3x3 double 5c = 1x27 char 2x2 double 2x2 double 2x2 double 3x3 double 2-14 t.college=computert.name=Wang Ming; t.number=19990101; t.sex=man; t.corse.name=MATLAB; t.corse.class=2004;t.corse(2
12、).name=maple; t.corse(2).class=2005;t2.college=computert2.name=zhang Ming; t2.number=19990101; t2.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=1 5 8;3 0 -4;b=6 -7 3;5 0 2;aba=baba6; %产生与A同型的01
13、逻辑值数组x=find(a6);a(L) % 把L中逻辑值1对应的A元素取出x %大于6的元素在原数组中的序号3-4 执行结果为：x = 0 0 0 0 0 0 0 1 1y = 1 1 1 1 0 0 0 0 0L = 1 1 1 1 0 0 0 1 1aL = 1.0000 2.1000 3.2000 4.3000 8.7000 9.8000xh =1 2 3 4 8 93-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
14、(lb)any_lab=any(la,lb)all_lab=all(la,lb)执行结果为：islogical_la = 1islogical_lb = 1any_lab = 1 1 1 1 1 1 1 1 1 1all_lab =1 0 1 1 1 1 1 1 1 13-6a=magic(3);b=magic(4);P1=rem(a,3)=0;P2=rem(b,3)=0;a(P1)b(P2)3-7 执行结果为：all_a = 1 0 1 1 1any_a = 1 1 1 1 1isnan_a = 0 0 0 0 0 0 0 0 0 0isinf_a = 0 0 0 0 0 0 0 0 0 0
15、isfinite_a = 1 1 1 1 1 1 1 1 1 13-8 a=input(Please input a matrix? );if isreal(a)=1 sign(a)else disp(The matrix isnot real matrix)end3-9x=input(请输入x的值:);if isreal(x)=0 disp(No defination)elseif x=-1 & x1|x-1 y=1+x else y=x.2+1 end case 0 disp(No defination)end 3-10x=input(please enter x value:);if x
16、=0&x=2&x=3&x10 y=x.2-5*x+6else y=x.2-x-1end3-11 执行结果：a = 2 -1 0 0 0 0 -1 2 -1 0 0 0 0 -1 2 -1 0 0 0 0 -1 2 -1 0 0 0 0 -1 2 -13-12s= 108 3-13 程序1：y=0;for n=1:100 f=1/n/n/n; y=y+f; n=n+1;endy程序2：y=0;n=1;while n x4_4(1 -6 3 2;3 -5 4 0;-1 -11 2 4)no inversepinv_a = -0.1000 0.1860 -0.0300 0.5000 -0.2500
17、-0.2500 0.7000 -0.2020 -0.2900 1.0000 -0.5400 -0.3000no detrank_a = 3 x4_4(1 2 1 0;6 2 4 1;0 2 1 0;3 1 4 1)det_a = 1inv_a = 1.0000 0 -1.0000 0 -3.0000 1.0000 3.0000 -1.0000 6.0000 -2.0000 -5.0000 2.0000 -24.0000 7.0000 20.0000 -6.0000rank_a =44-5 a=1 1 3 2;-1 1 -1 3;5 -2 8 9;-1 3 1 7;rank(a)ans =3线性
18、相关4-6 a=2 2 -1;1 -2 4;5 8 2; X,D=eig(a)X = 0.0000 -0.7160 0.2760 -0.4472 0.4200 -0.7472 -0.8944 0.5576 0.6045D = 6.0000 0 0 0 1.6056 0 0 0 -5.60564-7 A=6 5 -2 5;9 -1 4 -1;3 4 2 -2;3 -9 2 0;b=-4 13 1 11; inv(A)*bA=6 5 -2 5;9 -1 4 -1;3 4 2 -2;3 -9 2 0;b=-4 13 1 11; Abans = -8.3333 2.0000 27.0000 18.00
19、00A=6 5 -2 5;9 -1 4 -1;3 4 2 -2;3 -9 2 0;b=-4 13 1 11;Ab=A,b; rref(Ab)ans = 1.0000 0 0 0 -8.3333 0 1.0000 0 0 2.0000 0 0 1.0000 0 27.0000 0 0 0 1.0000 18.0000 4-9x=1 1 3;4 3 2;1 2 5*inv(1 1 -1;2 1 0;1 1 1)x = -1 0 2 0 1 2 -1 -1 44-10 p=1 0 -2 -5;r=roots(p);for ii=1:length(r) if r(ii)=0&r(ii)=3 r(ii)
20、 endendans = 2.09464-11 （1）、的根；p1=3 -1 2 1 0 3;p2=1/3 1 -3 -1;roots(p1),roots(p2)ans = 0.7330 + 0.7416i 0.7330 - 0.7416i -0.8952 -0.1188 + 1.0066i -0.1188 - 1.0066ians = -4.7588 2.0642 -0.3054（2），和；p1+0 0 p2ans = 3.0000 -1.0000 2.3333 2.0000 -3.0000 2.0000 p1p1 = 3 -1 2 1 0 3conv(p1,p2)ans =1.0000 2
21、.6667 -9.3333 2.3333 -4.0000 -4.0000 2.0000 -9.0000 -3.0000r p k=residue(p1,p2)r = -792.2450 22.1035 -0.8584p = -4.7588 2.0642 -0.3054k = 9 -30 177（3）、的导数。polyder(p1)ans = 15 -4 6 2 0polyder(p2)ans = 1 2 -3P Q=polyder(p1,p2)P = 2.0000 8.6667 -38.0000 -4.3333 -8.0000 -12.0000 -8.0000 9.0000Q =0.1111
22、0.6667 -1.0000 -6.6667 7.0000 6.0000 4-12 （1）时求的值。P=3 2 1;polyval(p,5 7 9)ans = -108.9538 -219.0359 -367.1882（2）时，求的值。P=3 2 1;polyvalm(p,2 5 ;7 9 )ans = -181.7691 -251.4115 -351.9761 -533.74524-13 p=poly(6 3 8)p = 1 -17 90 -144roots(p)ans = 8.0000 6.00003.00004-14 poly(-0.5 -0.3+0.4*i -0.3-0.4*i)ans
23、 = 1.0000 1.1000 0.5500 0.12504-15 p=4 -12 -14 5 9; diff(p)ans = -16 -2 19 4p=4 -12 -14 5 9; polyint(p)ans = 0.8000 -3.0000 -4.6667 2.5000 9.0000 04-16 x=0:0.01:1;y=sin(x).*exp(-x/10);interp1(x,y,0.25)ans = 0.24134-17 x=1:10;y=log(x)/log(3);interp1(x,y,3.25)ans = 1.06554-18 x=-1 2 2.5 3 4;y=0.2 3 3.
24、5 3.6 2.5;f=polyfit(x,y,2)polyval(f,3.5)4-19 a=1 1 1;polyint(a)quad(x.*log(x),1,exp(1)ans =2.0973第5章5-1 （1）f=3*x2+5*x+2表示在给定x时，将3*x2+5*x+2的数值运算结果赋值给变量f，如果没有给定x则指示错误信息。（2）f=3*x2+5*x+2表示将字符串3*x2+5*x+2赋值给字符变量f，没有任何计算含义，因此也不对字符串中的内容做任何分析。（3）x=sym(x) f=3*x2+5*x+2表示x是一个符号变量， f=3*x2+5*x+2就具有了符号函数的意义，f也自然成为
25、符号变量了。5-2 syms phi1 phi2;y=simple(sin(phi1)*cos(phi2)-cos(phi1)*sin(phi2) 5-3 syms x;syms y;f=x2+exp(x+y)-y*log(x)-3;subs(f,x,y,2,4)5-4 syms a11 a12 a21 a22;A=a11,a12,a13;a21,a22,a23;a31,a32,a33AD=det(A) % 行列式AI=inv(A) % 逆AE=eig(A) % 特征值5-5 （1）求f,g的反函数和复合函数，（2）分别求f,g的微分，（3）分别求f,g的积分。syms x y t;f=sym
26、(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-6 syms x;f=x4-5*x3+5*x2+5*x-6;factor(f)horner(f)5-7 syms x;f=(x+1)3+(x-1)2+5*x-6collect(f,x)5-8 syms x;expand(x+1)6)5-9 syms(x,y);n,d=numden(1/(x3-1)+1/(x2+y
27、+1)+1/(x+y+1)+8)5-10求下列极限。 (1) x=sym(x);limit(1/(x+1)-3/(x3+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(x2*exp(1/x/x),x,0)(5) x=sym(x);limit(5*x2)/(1-x2)+2(1/x),x,inf) (6) x=sym(x);limit(x2-2*x+1)/(x3-x),x,1) 5-11求下列幂级数的和函数。（1） n=sym