Matlab程序設(shè)計(jì)教程(第二版)課后參考答案

1、Matlab程序設(shè)計(jì)教程（第二版）課后參考答案第一章實(shí)驗(yàn)1.實(shí)驗(yàn)一第1題2.自己驗(yàn)證總結(jié)Matlab的主要優(yōu)點(diǎn)3.實(shí)驗(yàn)一第2題4.實(shí)驗(yàn)一第3題5.網(wǎng)站思考練習(xí)1. 啟動(dòng)見書P5 退出見書P62. Matlab主要功能見書P2 4種功能3. 分行輸入行末尾加續(xù)行符，即三個(gè)點(diǎn)“”4. 見書P115. 直接在命令窗口輸入fac第二章實(shí)驗(yàn)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

2、.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實(shí)驗(yàn)二 第1題3實(shí)驗(yàn)二 第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

3、= 8.5175e+003條件數(shù)越趨近于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思考練習(xí)1. 在MATLAB中6+7i是一個(gè)復(fù)數(shù)常量，6+7*i則是一個(gè)表達(dá)式。i是虛數(shù)單位，而I是單位向量。

4、2. A.*B表示A和B單個(gè)元素之間對(duì)應(yīng)相乘，A*B是按矩陣乘法的規(guī)則乘。A./B表示A矩陣除以B矩陣的對(duì)應(yīng)元素，B.A等價(jià)于A./B。A/B表示A *inv（B）而BA表示inv(B)* A，對(duì)于矩陣運(yùn)算一般兩式不等。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.實(shí)驗(yàn)二 第3題6. A=34,NaN,Inf,-Inf,-pi,eps,0;all(A)ans = 0any(A)ans = 1isn

5、an(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.結(jié)構(gòu)矩陣建立P46 結(jié)構(gòu)矩陣名.成員名=表達(dá)式單元矩陣建立 P47 矩陣元素用大括號(hào)括起來(lái)。8用稀疏存儲(chǔ)方式可以可以節(jié)省內(nèi)存空間; 運(yùn)算規(guī)則上，稀疏矩陣與普通矩陣一樣，只是矩陣的存儲(chǔ)方式不同。在運(yùn)算過(guò)程中，稀疏存儲(chǔ)矩陣可以直接參與運(yùn)算。當(dāng)參與運(yùn)算的對(duì)象不全是稀疏存儲(chǔ)矩陣時(shí)，所得結(jié)果一般是完全存儲(chǔ)形式。第三章實(shí)驗(yàn)1.提示： 先根據(jù)輸入的3位整數(shù)分別求個(gè)、十、百位上的數(shù)字a、b、c，反向輸出整數(shù)：100*a+10*b+c2（1）m=i

6、nput('請(qǐng)輸入一個(gè)百分制的成績(jī)： ');if m>=0&m<=100switch fix(m/10) case 9,10 disp('這個(gè)成績(jī)所屬的成績(jī)等級(jí):A'); case 8 disp('這個(gè)成績(jī)所屬的成績(jī)等級(jí):B'); case 7 disp('這個(gè)成績(jī)所屬的成績(jī)等級(jí):C'); case 6 disp('這個(gè)成績(jī)所屬的成績(jī)等級(jí):D'); otherwise disp('這個(gè)成績(jī)所屬的成績(jī)等級(jí):E'); endelse disp('錯(cuò)誤:輸入的成績(jī)不合理'

7、;);end3a=1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20b=max(a)c=min(a)b=0c=0for i=1:20 if a(i)>b b=a(i); end if a(i)<c c=a(i); endendbc4a=-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:61y=(exp(0.3*a(i)-exp(-0.3*a(i)/2*sin(a(i)+0.3)+log(0.3+a(i)/2)end5（1

8、）A=0;n=100for i=1:nA=1/n/n+A;end（2）A=1;n=10000;for i=1:nA=(2*n*n*2)/(2*n-1)*(2*n+1)*A;end6A=rand(5,6)n=input('ÇëÊäÈëÒ»¸ön');if n>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

9、*(i+1);end然后在命令文件p3_11.m中調(diào)用f1.m和f2.mp3_11.m代碼如下:y1=f1(40)/(f1(30)+f1(20);y2=f2(40)/(f2(30)+f2(20);y1y28函數(shù)文件factor1.m代碼如下:function f=factor1(k)if k=1 f=1;else f=factor1(k-1)+k;end函數(shù)文件factor2.m代碼如下:function f=factor2(k)if k=1 f=1;else f=factor2(k-1)+2;end函數(shù)文件factor3.m代碼如下:function f=factor3(k)if k=1 f

10、=1;else f=factor3(k-1)+1/k;end思考練習(xí)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=c

11、hengji(A,B)m=A.*B;n=A*B;然后在命令文件p3_12.m中調(diào)用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);6 自己獨(dú)立完成第四章實(shí)驗(yà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(

12、'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)bar(t,y)subplot(2,2,2)stairs(t,y)subplot(2,2,3)stem(t,y)subplot(2,2,4)loglog(t,y, 'y')3(1)x=-2*pi:0.1:2*piy=5*cos(x)+4polar(x,y)

13、 (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;思考練習(xí)1 務(wù)必自己看課本完成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

14、(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)<1e-4);x1=x(k);y3=2*x1-0.5;plot(x1,y3,'rp')4x=-pi*6:0.1:pi*6y=sin(1./x)plot(x,y)fplot('sin(1./x)',-pi*6,pi*6,1e-4)5(1)x=1:0.1:10y=12./xpolar(x,y)(2)x=-pi/6:0.1:pi/6y=3*si

15、n(x).*cos(x)./(sin(x).3+cos(x).3)polar(x,y)6u,v=meshgrid(-pi:pi/100:pi);x=3*u.*sin(v);y=2*u.*cos(v);z=4*u.*u;mesh(x,y,z);(2) x,y=meshgrid(-3:6/100:3);z=-5./(1+x.*x+y.*y);mesh(x,y,z)第五章實(shí)驗(yàn)1 a=randn(10,5);b=mean(a)b=std(a,0,1)b=max(a);c=max(b,2)b=min(a);c=min(b,2)b=sum(a)c=sum(b,2)sort(a)sort(a,2,'descend&#