中南大学matlab课后练习习题.docxVIP

第二章 1·求下列表达式的值。 1） w=sqrt(2)*(1+*10^-6) w = 2） a=;b=5;c=; x=(2*pi*a+(c+b)/(pi+a*b*c)-exp(2))/(tan(b+c)+a); x x = 3） a=;b=; y=2*pi*a^2*[(1-pi/4)*b-4)*a]; y y = 4） t=[2,1-3i;5,]; z=1/2*exp(2*t)*log(t+sqrt(1+t^2)); z = +004 * - - 2，已知 a， b，求下列表达式的值。 a=[-1,5,-4;0,7,8;3,61,7];b=[8,3,-1;2,5,3;-3,2,0]; 1） a+6*b ans = 47 23 -10 12 37 26 -15 73 7 a^2-b+eye(3) ans = -18 -217 17 22 533 109 21 867 526 2） a*b ans = 14 14 16 -10 51 21 125 328 180 a.*b ans = -8 15 4 0 35 24 -9 122 0 b*a ans = -11 0 -15 7 228 53 3 -1 28 3） a/b ans = b\a ans = 4） [a,b] ans = -1 5 -4 8 3 -1 0 7 8 2 5 3 3 61 7 -3 2 0 [a([1,3],:);b^2] ans = -1 5 -4 3 61 7 73 37 1 17 37 13 -20 1 9 已知 a，完成下列操作。 a=[23,10,,0;41,-45,65,5;32,5,0,32;6,,54,]; 1） 输出 a 在 [10,25] 范围内的全部元素。 k=find(a10a25) a(k) k = 1 ans = 23 (2) 取出 a 前 3 行构成矩阵 b，前两列构成矩阵 c，右下角 3*2 子矩阵构成矩阵 d， b 与 c 的乘积构成矩阵 e。 b=a(1:3,:) b = 0 0 c=a(:,1:2) c = d=a(2:4,3:4) d = 0 e=b*c = +003 * (3) 分别求表达式额 ed,ed,e|d 和 ~e|~d 的值。 ed ans = 1 0 1 ed ans = 1 1 1 e|d ans = 1 1 1 ~e|~d ans = 0 0 0 4, 产生 5 阶希尔伯特矩阵 H 和 5 阶帕斯卡矩阵 P，且求其行列式的值 Hh和 Hp 以及他们的条件数 Th 和 Tp，判断哪个矩阵性能更好，为什么 format rat H=hilb(5) H = 1 1/2 1/3 1/4 1/5 1/2 1/3 1/4 1/5 1/6 1/3 1/4 1/5 1/6 1/7 1/4 1/5 1/6 1/7 1/8 1/5 1/6 1/7 1/8 1/9 P=pascal(5) P = 1 1 1 1 1 1 2 3 4 5 1 3 6 10 15 1 4 10 20 35 1 5 15 35 70 Hh=det(H) Hh = 1/0000 Hp=det(P) Hp = 1 Th=cond(H) Th = 476607 Tp=cond(P) Tp = 178868/21 P 矩阵性能更好，因为 Tp 更接近 1. 已知 A，求 A 的特征值及特征向量，并分析其数学意义。 A=[-29,6,18;20,5,12;-8,8,5] A = -29 6 18 20 5 12 -8 8 5 [V,D]=eig(A,nobalance) V = 1 361/1013 493/1574 -1227/1438 -1 1 290/593 551/788 331/713 D = -3595/142 0 0 0 -3755/357 0 0 0 4697/279 第三章 1． 从键盘输入一个 3 位整数，将它反向输出。如输入 639，输出为 936. clear m=input( 请输入一个三位数： ); m1=fix(m/100); m2=rem(fix(m/10),10); m3=rem(m,10); m=m1+m2*10+m3*100; disp(m) ； 请输入一个三位数： 453 354 2．输入一个百分制成绩，要求输出成绩等级 A、 B、 C、 D、 E。其中 90~100 分为 A， 80~89 分为 B， 70~79 分为 C,60~69 分为 D， 60 分以下为 E。要求： 1）分别用 if 语句和 switch 语句实现。 2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 c=input( 请输入一个百分制分数： ); if c=90c=100 % 用 if 语句