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

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*a2*(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.3580i2实验二 第1题3实验二 第2题4 H=hilb(5); P=pascal(5); Hh=det(H)Hh = 3.7493e-012Hp=det(P)Hp = 1Th=cond(H)Th = 4.7661e+005Tp=cond(P)Tp = 8.5175e+003条件数越趋近于1，矩阵的性能越好，所以帕斯卡矩阵性能更好。5 A=-29,6,18;20,5,12;-8,8,5A = -29 6 18 20 5 12 -8 8 5V,D=eig(A)V = 0.7130 0.2803 0.2733 -0.6084 -0.7867 0.8725 0.3487 0.5501 0.4050D = -25.3169 0 0 0 -10.5182 0 0 0 16.8351V为A的特征向量，D为A的特征值。它们之间满足A*V=V*D思考练习1. 在MATLAB中6+7i是一个复数常量，6+7*i则是一个表达式。i是虚数单位，而I是单位向量。2. A.*B表示A和B单个元素之间对应相乘，A*B是按矩阵乘法的规则乘。A./B表示A矩阵除以B矩阵的对应元素，B.A等价于A./B。A/B表示A *inv（B）而BA表示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 = 0any(A)ans = 1isnan(A)ans = 0 1 0 0 0 0 0isinf(A)ans = 0 0 1 1 0 0 0isfinite(A)ans = 1 0 0 0 1 1 17.结构矩阵建立P46 结构矩阵名.成员名=表达式单元矩阵建立 P47 矩阵元素用大括号括起来。8用稀疏存储方式可以可以节省内存空间; 运算规则上，稀疏矩阵与普通矩阵一样，只是矩阵的存储方式不同。在运算过程中，稀疏存储矩阵可以直接参与运算。当参与运算的对象不全是稀疏存储矩阵时，所得结果一般是完全存储形式。第三章实验1.提示： 先根据输入的3位整数分别求个、十、百位上的数字a、b、c，反向输出整数：100*a+10*b+c2（1）m=input(请输入一个百分制的成绩： );if m=0&mb b=a(i); end if a(i)5n=5A(n,n)else A(n)end7f1.m代码如下:function m=f1(n)m=n+10*log(n2+5);f2.m 代码如下:function m=f2(n)m=0;for i=1:n m=m+i*(i+1);end然后在命令文件p3_11.m中调用f1.m和f2.mp3_11.m代码如下:y1=f1(40)/(f1(30)+f1(20);y2=f2(40)/(f2(30)+f2(20);y1y28函数文件factor1.m代码如下:function f=factor1(k)if k=1 f=1;else f=factor1(k-1)+k;end函数文件factor2.m代码如下:function f=factor2(k)if k=1 f=1;else f=factor2(k-1)+2;end函数文件factor3.m代码如下:function f=factor3(k)if k=1 f=1;else f=factor3(k-1)+1/k;end思考练习1见课本3.12N=1,2,3,4,5A=N.*2B=N./4C=1./ND=1./(N./2)./(N./2)/43B=0for i=1:20A(i)=fix(99-10+1)*rand(1)+10);B=A(i)+B;endB=B/20for i=1:20if A(i)B&rem(A(i),2)=0A(i)endend4(1)A=1for i=2:100A=(-1)(i+1)*1/(2*n-1)+A;endA (2)A=0for i=1:1000A=1/(4)i)+A;endA5function m,n=chengji(A,B)m=A.*B;n=A*B;然后在命令文件p3_12.m中调用chengji.m文件a=1 2 3;4 5 6;7 8 9b=-1 0 1;1 -1 0;0 1 1m,n=chengji(a,b);disp(a.*b=);disp(m);disp(a*b= );disp(n);第四章实验1（1）x=-12:0.001:12;y=x-x.3/6;plot(x,y);xlabel(x轴);ylabel(y轴);title(Y=X-X3/3!);（2）ezplot(x2+2*y2-64,-8,8,-8,8)xlabel(x轴);ylabel(y轴);axis equaltitle(x2+2y2=64);2t=-pi:0.1:pi;y=1./(1+exp(-t);plot(t,y);subplot(2,2,1)stairs(t,y)subplot(2,2,2)bar(t,y)subplot(2,2,3)stem(t,y)subplot(2,2,4)fill(t,y, y)3(1)x=-2*pi:0.1:2*piy=5*cos(x)+4polar(x,y) (2)x=-pi/3:.01:pi/3a=sin(x)b=cos(x)y=a.*a.*5./bpolar(x,y)4t=0:pi/100:2*pi;x=exp(-t/20).*cos(t);y=exp(-t/20).*sin(t);z=t;plot3(x,y,z)(2) x,y=meshgrid(-5:5);z=ones(11);surf(x,y,z);shading interp;思考练习1 务必自己看课本完成2(1)x=1:0.1:10a=exp(x.2)/2)y=a.*1/2/piplot(x,y) (2)t=-pi:0.1:2*piX=t.*sin(t)Y=t.*cos(t)plot(X,Y)3t=0:0.00005:pi;x=sin(3*t).*cos(t);y1=2*x-0.5;y2=sin(3*t).*sin(t);plot(x,y2,x,y1);hold onk=find(abs(y2-y1)0.5 c=c+1; endendb