2022年Matlab程序设计教程课后参考答案 .pdf

《2022年Matlab程序设计教程课后参考答案 .pdf》由会员分享，可在线阅读，更多相关《2022年Matlab程序设计教程课后参考答案 .pdf（10页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、1 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.003

2、4i 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

3、 0.5501 0.4050 名师资料总结-精品资料欢迎下载-名师精心整理-第 1 页，共 10 页 -2 D=-25.3169 0 0 0 -10.5182 0 0 0 16.8351 V为 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，对于矩阵运算

4、一般两式不等。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

5、矩阵元素用大括号括起来。8用稀疏存储方式可以可以节省内存空间;运算规则上，稀疏矩阵与普通矩阵一样，只是矩阵的存储方式不同。在运算过程中，稀疏存储矩阵可以直接参与运算。当参与运算的对象不全是稀疏存储矩阵时，所得结果一般是完全存储形式。第三章实验1.提示：先根据输入的3 位整数分别求个、十、百位上的数字a、b、c，反向输出整数：100*a+10*b+c2（1）名师资料总结-精品资料欢迎下载-名师精心整理-第 2 页，共 10 页 -3 m=input(请输入一个百分制的成绩：);if m=0&mb b=a(i);endif a(i)5 n=5 A(n,n)elseA(n)end7 f1.m 代码如

6、下: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.m p3_11.m代码如下:y1=f1(40)/(f1(30)+f1(20);y2=f2(40)/(f2(30)+f2(20);y1 y2 8 1001kk函数文件 factor1.m代码如下:名师资料总结-精品资料欢迎下载-名师精心整理-第 4 页，共 10 页 -5 function f=factor1(k)if k=1 f=1;else f=factor1

7、(k-1)+k;end5012kk函数文件 factor2.m代码如下:function f=factor2(k)if k=1 f=1;else f=factor2(k-1)+2;end101/1kk函数文件 factor3.m代码如下:function f=factor3(k)if k=1 f=1;else f=factor3(k-1)+1/k;end 思考练习1见课本3.1 2 N=1,2,3,4,5 A=N.*2 B=N./4 C=1./N D=1./(N./2)./(N./2)/4 3 B=0 for i=1:20 A(i)=fix(99-10+1)*rand(1)+10);B=A(i

8、)+B;endB=B/20 for i=1:20 if A(i)B&rem(A(i),2)=0 A(i)名师资料总结-精品资料欢迎下载-名师精心整理-第 5 页，共 10 页 -6 endend4(1)A=1 for i=2:100 A=(-1)(i+1)*1/(2*n-1)+A;endA (2)A=0 for i=1:1000 A=1/(4)i)+A;endA 5 function 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 9 b=-1 0 1;1-1 0;0 1 1 m,n=ch

9、engji(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 equal title(x2+2y2=64);2 t=-pi:0.1:pi;y=1./(1+exp(-t);名师资料总结-精品资料欢迎下载-名师精心整理-第 6 页，共 10 页 -7 plot(t,y);subplo

10、t(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*pi y=5*cos(x)+4 polar(x,y)(2)x=-pi/3:.01:pi/3 a=sin(x)b=cos(x)y=a.*a.*5./b polar(x,y)4 t=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)

11、;surf(x,y,z);shading interp;名师资料总结-精品资料欢迎下载-名师精心整理-第 7 页，共 10 页 -8 思考练习1 务必自己看课本完成2(1)x=1:0.1:10 a=exp(x.2)/2)y=a.*1/2/pi plot(x,y)(2)t=-pi:0.1:2*pi X=t.*sin(t)Y=t.*cos(t)plot(X,Y)3 t=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 on k=find(abs(y2-y1)0.5 c=c+1;end end b=mean(a)b=std(a,0,2)b=max(a)b=min(a)c/30000 5 a=fix(100-0+1)*rand(5,100);b,c=max(a,2)b,c=min(a,2)mean(a,2)std(a,0,2)b,c=sort(a,2,descend)zcj=sum(a,1)xsxh=1:100 名师资料总结-精品资料欢迎下载-名师精心整理-第 10 页，共 10 页 -