Matlab程设计教程(第二版)课后参考答案.docVIP

Matlab程序设计教程（第二版）课后参考答案 第二章 实验 1.(1) w=sqrt(2)*(1+0.34245*10^(-6)) w = 1.4142 (2) a=3.5; b=5; c=-9.8; x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a) x =0.9829 (3) a=3.32; b=-7.9; y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a) y = -128.4271 (4) t=[2,1-3i;5,-0.65]; z=0.5*exp(2*t)*log(t+sqrt(1+t*t)) z = 1.0e+004 * 0.0048 + 0.0002i 0.0048 - 0.0034i 1.5899 2.0090 - 1.3580i 2．实验二 第1题 3．实验二 第2题 4． H=hilb(5); P=pascal(5); Hh=det(H) Hh = 3.7493e-012 Hp=det(P) Hp = 1 Th=cond(H) Th = 4.7661e+005 Tp=cond(P) Tp = 8.5175e+003 条件数越趋近于1，矩阵的性能越好，所以帕斯卡矩阵性能更好。 5． A=[-29,6,18;20,5,12;-8,8,5] A = -29 6 18 20 5 12 -8 8 5 [V,D]=eig(A) V = 0.7130 0.2803 0.2733 -0.6084 -0.7867 0.8725 0.3487 0.5501 0.4050 D = -25.3169 0 0 0 -10.5182 0 0 0 16.8351 V为A的特征向量，D为A的特征值。它们之间满足A*V=V*D 思考练习 在MATLAB中6+7i是一个复数常量，6+7*i则是一个表达式。i是虚数单位，而I是单位向量。 A.*B表示A和B单个元素之间对应相乘，A*B是按矩阵乘法的规则乘。 A./B表示A矩阵除以B矩阵的对应元素，B.\A等价于A./B。 A/B表示A *inv（B）而B\A表示inv(B)* A，对于矩阵运算一般两式不等。 3.(1) A(7)=[] (2) abs（x） (3) reshape(x,3,4) (4) k=find(t==0); t(k)=0 (5) ones(size(A)) (6) diag(A) B= diag(diag(A)) 4. 3+randn(25,20) 5.实验二 第3题 6. A=[34,NaN,Inf,-Inf,-pi,eps,0]; all(A) ans = 0 any(A) ans = 1 isnan(A) ans = 0 1 0 0 0 0 0 isinf(A) ans = 0 0 1 1 0 0 0 isfinite(A) ans = 1 0 0 0 1 1 1 7.结构矩阵建立P46 结构矩阵名.成员名=表达式 单元矩阵建立 P47 矩阵元素用大括号括起来。 8．用稀疏存储方式可以可以节省内存空间; 运算规则上，稀疏矩阵与普通矩阵一样，只是矩阵的存储方式不同。在运算过程中，稀疏存储矩阵可以直接参与运算。当参与运算的对象不全是稀疏存储矩阵时，所得结果一般是完全存储形式。 第三章 实验 1.提示： 先根据输入的3位整数分别求个、十、百位上的数字a、b、c，反向输出整数：100*a+10*b+c 2（1） m=input(请输入一个百分制的成绩： ); if m=0m=100 switch fix(m/10) case {9,10} disp(这个成绩所属的成绩等级:A); case {8} disp(这个成绩所属的成绩等级:B); case {7} disp(这个成绩所属的成绩等级:C); case {6} disp(这个成绩所属的成绩等级:D); otherwise disp(这个成绩所属的成绩等级:E); end else disp(错误:输入的成绩不合理); end 3 a=[1