Matlab程序设计教程(第二版)课后参考答案(10页).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.0034i1.58992.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 =-296 18205 12-885 [V,D]=eig(A) V =0.7130 0.2803 0.2733-0.6084 -0.7867 0.87250.3487 0.5501 0.4050 D = -25.3169000 -10.5182000 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 = 0100000 isinf(A) ans = 0011000 isfinite(A) ans = 1000111 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);otherwisedisp(这个成绩所属的成绩等级:E);end elsedisp(错误:输入的成绩不合理); end 3 a=[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20] b=max(a) c=min(a) b=0 c=0 for i=1:20if a(i)bb=a(i);end if a(i)cc=a(i);end end b c 4 a=[-3:0.1:3]; b=((exp(0.3.*a)-exp(-0.3.*a))./2) c=sin(a+0.3) y=b.*c+log((0.3+a)./2)for i=1:61 y=(exp(0.3*a(i))-exp(-0.3*a(i)))/2*sin(a(i)+0.3)+log((0.3+a(i))/2) end 5 （1） A=0; n=