matlab基礎(chǔ)與應(yīng)用部分習(xí)習(xí)題答案

1、作業(yè)一4、寫出完成下列操作的命令。（1）將矩陣A第25行中第1,3,5列元素賦給矩陣B。>> M=0:1:48;>> A=reshape(M,7,7)A = 0 7 14 21 28 35 42 1 8 15 22 29 36 43 2 9 16 23 30 37 44 3 10 17 24 31 38 45 4 11 18 25 32 39 46 5 12 19 26 33 40 47 6 13 20 27 34 41 48>> B=A(2:5,1:2:5)B = 1 15 29 2 16 30 3 17 31 4 18 32（2）刪除矩陣A的第七行元素。

2、>> A(7,:)=A = 0 7 14 21 28 35 42 1 8 15 22 29 36 43 2 9 16 23 30 37 44 3 10 17 24 31 38 45 4 11 18 25 32 39 46 5 12 19 26 33 40 47(3)將矩陣A的每個(gè)元素值加30。>> A=A+30A = 30 37 44 51 58 65 72 31 38 45 52 59 66 73 32 39 46 53 60 67 74 33 40 47 54 61 68 75 34 41 48 55 62 69 7635 42 49 56 63 70 77(4求矩

3、陣A的大小和維素。sizeA = size(A) dA = ndims(A)sizeA = 6 7dA = 2（5）將向量t的0元素用機(jī)器零來代替。>> t=1 2 3 4 0 5;t = 1 2 3 4 0 5>> find(t=0)ans = 5>> t(5)=epst = （6）將含有12個(gè)元素的向量x轉(zhuǎn)換成3*4矩陣。>> x=0:11x = 0 1 2 3 4 5 6 7 8 9 10 11>> y=reshape(x,3,4)y = 0 3 6 9 1 4 7 10 2 5 8 11（7）求一個(gè)字符串的ASCII。>

4、> abs('a')ans =97（8）求一個(gè)ASCII所對(duì)應(yīng)的字符。>> char(99)ans =c5、已知A97 67 34 10-78 75 65 532 5 -23 -590 = -12 54 7寫出完成下列操作的命令。（1） 取出A的前3列構(gòu)成矩陣B，前兩列構(gòu)成矩陣C，右下角3*2子矩陣構(gòu)成矩陣D，B與C的乘積構(gòu)成矩陣E。>> B=A(1:3,:)B = 97 67 34 10 -78 75 65 5 32 5 -23 -59>> C=A(:,1:2)C = 97 67 -78 75 32 5 0 -12>> D

5、=A(2:4,3:4)D = 65 5 -23 -5954 7>> E=B*CE = 5271 11574 -11336 664 1978 3112（2） 輸出【50,100】范圍的全部元素。>> find(A>=50 & A<=100)ans = 1 5 6 1012>> output = A(1) A(5) A(6) A(10) A(12)output = 97 67 75 65 54二、實(shí)驗(yàn)題>> x = -74*pi/180;y=-27*pi/180; sin(x*x+y*y)/(sqrt(tan(abs(x+y)+p

6、i)ans = - >> a = :, exp.*a).*sin(a+ >> x=2 4; 5; *log(x+sqrt(1+x.*x)ans = >> A = 3 54 2;34 -45 7;87 90 15, B = 1 -2 67;2 8 74;9 3 0>> A*B A.*Bans = 129 432 4197 7 -407 -1052 402 591 12489ans = 3 -108 134 68 -360 518 783 270 0>> A3 A.3>> A/B BA>> A,B A(1,3,:)

7、;B2>> a=1+2i,b=3+4i,c=exp(pi/6*i) c+a*b/(a+b)作業(yè)二1、 從鍵盤輸入一個(gè)4位整數(shù)，按如下規(guī)則加密后輸出。加密規(guī)則：每位數(shù)字都加上7，然后用和除以10的余數(shù)取代該數(shù)字；然后將第1位數(shù)和第3位數(shù)交換，第2位數(shù)和第4位數(shù)交換。M 文件名為：程序：x=input('輸入一個(gè)四位整數(shù)：');y=fix(x/1000),mod(fix(x/100),10),mod(fix(x/10),10),mod(x,10)z=mod(y+7),10)x=z(3)*1000+z(4)*100+z(1)*10+z(2)結(jié)果：輸入一個(gè)四位整數(shù)：1234

8、y =1 2 3 4z =8 9 0 1x =1892、 硅谷公司員工的工資計(jì)算方法如下。（1） 工作時(shí)數(shù)超過120h者，超過部分加發(fā)15%。（2） 工作時(shí)數(shù)低于60h者，扣發(fā)700元。（3） 其余按每小時(shí)84元計(jì)發(fā)。試編程按輸入的工號(hào)和該號(hào)員工的工時(shí)數(shù)，計(jì)算應(yīng)發(fā)工資。程序：gonghao=input('輸入工號(hào)');hour=input('輸入工時(shí)');every=84;if hour>120gongzi=every*120+* every *(hour-120);elseif hour<60gongzi=every*hour-700;elsego

9、ngzi=devery*hour;endformat bank;display(gonghao,gongzi)（1）結(jié)果：輸入工號(hào)12輸入工時(shí)130ans = （2）結(jié)果：輸入工號(hào)12輸入工時(shí)56ans = （3）結(jié)果：輸入工號(hào)12輸入工時(shí)80ans = 3、 根據(jù)，求的近似值。當(dāng)n分別取100、1000、10000時(shí)，結(jié)果是多少要求：分別用循環(huán)結(jié)構(gòu)和向量運(yùn)算（使用sun函數(shù)）來實(shí)現(xiàn)。（1）程序：（循環(huán)結(jié)構(gòu)）：n=input('input n:');s=0;for i=1:ns=s+1/i2;pai= sqrt(s*6);enddisplay(pai)input n:100pa

10、i = （2）程序：（向量結(jié)構(gòu)）：n=input('input n:');i=1:n;pai= sqrt(sum(1./i.2)*6);display(pai);結(jié)果：input n:100pai = 4、 根據(jù)，求：（1） y3時(shí)的最大n值；（2） 與（1）的n值對(duì)應(yīng)的y值。y=0;k=0;while y<3 k=k+1; y=y+1/(2*k-1);enddisplay(k-1,y-1/(2*k-1)結(jié)果：ans = 5、 考慮一下迭代公式：其中，a、b為正的常數(shù)。（1） 編寫程序求迭代的結(jié)果，迭代的終止條件為|，迭代初值，迭代次數(shù)不超過500次。程序：x0=;a=i

11、nput('a=');b=input('b=');x= a/(b+x0); m=1;while abs(x-x0)>1e-5 & m<500x0=x;x=a/(b+x0);m=m+1;enddisplay(m,x);display(-b+sqrt(b2+4*a)/2,(-b-sqrt(b2+4*a)/2); 結(jié)果：a=5b=5ans = ans = （2） 如果迭代過程收斂于r，那么r的準(zhǔn)確值是，當(dāng)（a，b）的值取（1，1）、（8,3）、（10,）時(shí)，分別對(duì)迭代結(jié)果和準(zhǔn)確值進(jìn)行比較。a=1b=1ans = ans = a=8b=3ans =

12、ans = a=10b=ans = ans =3.11276、已知（1） 當(dāng)時(shí)，y的值是多少程序：該題使用了function的調(diào)用。函數(shù)文件function b=f(n);b=n+10*log(n2+5);結(jié)果：y=f(40)/(f(30)+f(20)y = （2） 當(dāng)時(shí)，y的值是多少程序：函數(shù)文件function f=fun(n)f=sum(1:n.*(1:n+1);結(jié)果：y=fun(40)/(fun(30)+fun(20)y =作業(yè)三1、 繪制下列曲線。（1） （2）（3） （4）解：（1）程序：x=-5:5; y=x-x.3/(3*2*1); plot(x,y);(2)程序：x=-10:

13、10; y= exp(-x.*x/2)/(2*pi);>> plot(x,y);（3） 程序：f='x.2+2*y.2-64'ezplot(f,-10,10,-10,10);ezplot('x.2+2*y.2-64',-10,10,-10,10);（4）程序：t=0:8*pi;>> x=t.*sin(t);>> y=t.*cos(t);>> plot(x,y);已知，完成下列操作：（1） 在同一坐標(biāo)系下用不同的顏色和線型繪制兩條曲線，給曲線添加文字說明；解：程序：x1=linspace(0,1,100);>&

14、gt; y1=2*;>> t=linspace(0,pi,100);>> x=sin(3*t).*cos(t);>> y=sin(3*t).*sin(t);>> plot(x1,y1,'k:',x,y,'bp');>> text,1,'y=');>> text,'x=sin(3t)cos(t)');>> text,'x=sin(3t)sin(t)');（2） 以子圖形式，分別用條形圖和散點(diǎn)圖繪制兩條曲線，并為各子圖添加函數(shù)標(biāo)題。程序

15、：subplot(1,2,1);bar(x1,'group');title('y1=');subplot(1,2,2);barh(t,'stack');title('stack');subplot(1,2,1);>> scatter(x1,y1,10);>> title('y=');>> subplot(1,2,2);>> scatter(x,y,10);>> title('x=sin(3t)cos(t),y=sin(3t)cos(t)')

16、;2、 分別用plot和fplot函數(shù)繪制函數(shù)的曲線，分別分析兩條曲線的差別。解：程序：subplot(1,2,1);x=1:1:100;y=sin(1./x);plot(x,y);subplot(1,2,2);>> fplot('sin(1/x)',1,100);3、 設(shè)，在同一圖形窗口采用子圖的形式繪制條形圖、階梯圖、桿圖和對(duì)數(shù)坐標(biāo)圖。解：程序：t=-pi:pi;>> y=1./(1+exp(-t);>> bar(t,y);>> subplot(2,2,2);>> stairs(t,y);>> subp

17、lot(2,2,3);>> stem(t,y);>> subplot(2,2,4);>> semilogy(t,y);結(jié)果：4、 繪制下列極坐標(biāo)圖。（1） （2），解：（1）程序：theta=linspace(-pi,pi,100);ro=5.*cos(theta)+4;polar(theta,ro);（2）程序：fi=linspace(0,2*pi,100);a=1r=a.*(1+cos(fi);polar(fi,r);6、繪制下列三維圖線。（1） （2）解：（1）程序：t=0:pi/20:2*pi;x=exp(-t./20).*cos(t);y=exp(

18、-t./20).*sin(t);z=t;plot3(x,y,z);結(jié)果：（3） 程序：t=0:1;x=t;y=t.2;z=t.3;plot3(x,y,z);結(jié)果：7、已知，繪制其在范圍內(nèi)的曲面圖和等高線。解：程序：x=-30:0;y=0:30;x,y=meshgrid(x,y);z=10.*sin(sqrt(x.2+y.2)./sqrt(1+x.2+y.2);meshc(x,y,z);結(jié)果：8、已知，繪制其曲面圖，并將部分鏤空。解：程序：x=linspace(-3,3,100);y=linspace(-3,3,100);x y=meshgrid(x,y);fxy=-5./(1+x.2+y.2)

19、;i=find(abs(x)<= & abs(y)<=; fxy(i)=NaN;surf(x,y,fxy)9、繪制曲面圖形，應(yīng)用插值著色處理，并設(shè)置光照效果。解：程序：u=linspace(1,100);v=linspace(-pi,pi,1000);u,v=meshgrid(u,v);x=3.*u.*sin(v);y=2.*u.*cos(v);z=4*u.2;x=3*u.*sin(v);y=2*u.*cos(v);z=4*u.2;surf(x,y,z);shading interp;light('position',1，1，1);10、設(shè)計(jì)一個(gè)藍(lán)色球體沿正

20、弦曲線運(yùn)動(dòng)的動(dòng)畫。解：程序：t=0:pi/250:10*pi;y=sin(t);plot(t,y,'m');h=line('color',0,0,1,'marker','.','markersize',50);for i=1:length(t); set(h,'xdata',t(i),'ydata',y(i); M(i)=getframe;endmovie(M); 第四次作業(yè)P103 思考題5、寫出完成下列操作的命令。（1）建立3階單位矩陣A。A=eye(3)A = 1 0 0 0

21、1 0 0 0 1（2）建立5*6隨機(jī)矩陣A，其元素為100，200范圍內(nèi)的隨機(jī)整數(shù)。a=100; b=200; A=fix(a+(b-a)*rand(5,6)A = 195 176 161 140 105 120 123 145 179 193 135 119 160 101 192 191 181 160 148 182 173 141 100 127 189 144 117 189 113 119（3）產(chǎn)生均值為1、方差為的500個(gè)正態(tài)分布的隨機(jī)數(shù)。y=1+sqrt*randn(25,20)（4）產(chǎn)生和A同樣大小的幺矩陣。B=ones(size(A)（5）將矩陣A的對(duì)角線的元素加30。A

22、=1,2,3,4,5,6;5,6,2,1,3,5;7,8,9,4,5,6A = 1 2 3 4 5 6 5 6 2 1 3 5 7 8 9 4 5 6>> A=A+30*eye(size(A)A = 31 2 3 4 5 6 5 36 2 1 3 5 7 8 39 4 5 6（6）從矩陣A提取主對(duì)角線元素，并以這些元素構(gòu)成對(duì)角陣B。C=diag(diag(A)C = 31 0 0 0 36 0 0 0 39實(shí)驗(yàn)題：4、分別用矩陣求逆、左除運(yùn)算符、矩陣分解等方法求下列對(duì)角方程組，并分析程序的執(zhí)行時(shí)間。其中A = 2 -1 0 0 0 -1 2 -1 0 0 0 -1 2 -1 0 0

23、 0 -1 2 -1 0 0 0 -1 2>> b=1;0;0;0;0b = 1 0 0 0 0解：程序：%矩陣求逆ticx1=inv(A)*b toc%左除運(yùn)算ticx2=Ab toc%矩陣分解%Q、R分解ticQ,R=qr(A);x4=R(Qb)；toc%L、U分解ticL,U=lu(A);x3=U(Lb)toc%Cholesky分解ticR=chol(A);x5=R(R'b)toc結(jié)果：P121 第五章實(shí)驗(yàn)題1、 利用，MATLAB提供的rand函數(shù)生成30000個(gè)符合均勻的隨機(jī)數(shù)，然后檢驗(yàn)隨機(jī)數(shù)的性質(zhì)。（1） 均值和標(biāo)準(zhǔn)方差。解：A=rand(1,30000);&g

24、t;> mean(A)std(A)（2） 最大元素和最小元素。max(A)min(A)（3） 大于的隨機(jī)數(shù)個(gè)數(shù)占總數(shù)的百分比。size(find(A>)/size(A)ans =2、在某處測(cè)的海洋不同深部處水溫度數(shù)據(jù)如表。用插值法求出水深500m，900m和1500m處的水溫（）。水深（m）4667159504221635水溫（）解：程序：m=466,715,950,1422,1635;w=,;h=500,900,1500;wi=interp1(沒,w,h,'spline')3、用5次多項(xiàng)式p（x）在區(qū)間【1,10】區(qū)間內(nèi)逼近函數(shù)lgx，并繪制出lgx和p（x）在【

25、1,101】區(qū)間的函數(shù)曲線。解：x=1:101;y=log(x);f=polyfit(x,y,5);yi=polyval(f,x);plot(x,y,'bp',x,yi,'k-p')4、已知，取N=64，對(duì)t從05s采樣，用fft作快速傅立葉變換，并繪制相應(yīng)的振幅-頻率圖。解：N=64;T=5;t=linspace(0,T,N);x=exp(-t);dt=t(2)-t(1);f=1/dt;X=fft(x);F=X(1:N/2+1);f=f*(0:N/2)/N;plot(f,abs(F),'r-+')xlabel('Frequency

26、9;)ylabel('|F(k)|')5、有3個(gè)多項(xiàng)式，試進(jìn)行下列操作。（1）求。解、程序：p1=1,2,0,7;p2=0,0,1,-2;p3=1,0,5,1;p12=conv(p1,p2);p=p12+zeros(1,size(p12,2)-size(p3,2),p3結(jié)果：p = 0 0 1 1 -4 12 -13（2）求p（x）的根。程序：roots(p)結(jié)果：ans = + - （3）設(shè)A矩陣，當(dāng)以矩陣A的每一元素為自變量時(shí)，求的值，當(dāng)以矩陣A為自變量時(shí)，求的根。程序：A=-1,4,3;2,1,5;0,5,6;px=polyval(p,A)pxm=polyvalm(p,A

27、)結(jié)果：px = -29 291 95 19 -3 697 -13 697 1427pxm = 391 2084 3273 502 2693 4207 720 3775 58926、求非線性方程的數(shù)值解。（1）、求在附近的根。解：程序：z=fzero('3*x-sin(x)+1',0)（2）、，初值。建立函數(shù)文件：function F=fun(X)x=X(1);y=X(2);F(1)=x*x+y*y-9;F(2)=x+y-1;在命令窗口輸入：x=fsolve(fun,3,0',optimset('Display','off')第7章符號(hào)與

28、計(jì)算課后實(shí)驗(yàn)題1、已知A=求矩陣的秩、逆矩陣和行列式的值。（這里的A是符號(hào)矩陣）解：程序及結(jié)果：A=sym('1,2,3;x,y,z;3,2,1');>> rank(A) %矩陣秩的值ans = 3>> inv(A) %逆矩陣的值 ans = 1/4*(y-2*z)/(-2*y+z+x), 1/(-2*y+z+x), -1/4*(-2*z+3*y)/(-2*y+z+x) -1/4*(x-3*z)/(-2*y+z+x), -2/(-2*y+z+x), 1/4*(-z+3*x)/(-2*y+z+x) 1/4*(2*x-3*y)/(-2*y+z+x), 1/(

29、-2*y+z+x), -1/4*(-y+2*x)/(-2*y+z+x)>> det(A) %矩陣的行列值ans = -8*y+4*z+4*x2、求函數(shù)的符號(hào)導(dǎo)數(shù)。（1）求。解：程序及結(jié)果：syms x;>> y=sqrt(x+sqrt(x+sqrt(x);>> diff(y) %一階導(dǎo)數(shù) ans =1/2/(x+(x+x(1/2)(1/2)(1/2)*(1+1/2/(x+x(1/2)(1/2)*(1+1/2/x(1/2)>> diff(y,x,2) %二階導(dǎo)數(shù)ans =-1/4/(x+(x+x(1/2)(1/2)(3/2)*(1+1/2/(x+x

30、(1/2)(1/2)*(1+1/2/x(1/2)2+1/2/(x+(x+x(1/2)(1/2)(1/2)*(-1/4/(x+x(1/2)(3/2)*(1+1/2/x(1/2)2-1/8/(x+x(1/2)(1/2)/x(3/2)（2）已知，求。解：程序：syms x y;>> fxy=sin(x2*y)*exp(-x2-y);>> c=diff(diff(fxy,x),y) 結(jié)果：c = -2*sin(x2*y)*x3*y*exp(-x2-y)+2*cos(x2*y)*x*exp(-x2-y)-2*cos(x2*y)*x*y*exp(-x2-y)-2*cos(x2*y)

31、*x3*exp(-x2-y)+2*sin(x2*y)*x*exp(-x2-y)3、求積分。（1） （2）（1）解：程序：f=sym('1/(x4+1)');>> s1=int(f) （2）解：程序：syms x t;>> f=sym('(-2*x2+1)/(2*x2-3*x+1)2');>> int(f,cos(t),exp(2*t)結(jié)果：ans = -(-2*exp(2*t)*cos(t)2-exp(2*t)+cos(t)+2*exp(2*t)2*cos(t)/(2*exp(2*t)-1)/(exp(2*t)-1)/(2*c

32、os(t)-1)/(cos(t)-1)4、求級(jí)數(shù)的和函數(shù)，并求之和。解：程序：syms x n;>> s1=symsum(2/(2*n+1)*(2*x+1)(2*n+1),n,0,inf)結(jié)果：s1 = 1/(2*x+1)*(4*x2+4*x+1)(1/2)*log(1+1/(4*x2+4*x+1)(1/2)/(1-1/(4*x2+4*x+1)(1/2)s2=symsum(2/(2*n+1)*(2*x+1)(2*n+1),n,1,5) s2 = 2/3/(2*x+1)3+2/5/(2*x+1)5+2/7/(2*x+1)7+2/9/(2*x+1)9+2/11/(2*x+1)11 &g

33、t;> s=s1+s2 s = 1/(2*x+1)*(4*x2+4*x+1)(1/2)*log(1+1/(4*x2+4*x+1)(1/2)/(1-1/(4*x2+4*x+1)(1/2)+2/3/(2*x+1)3+2/5/(2*x+1)5+2/7/(2*x+1)7+2/9/(2*x+1)9+2/11/(2*x+1)115、求函數(shù)在的泰勒展開式。(1) (2)（1）解：程序：sym ('x');>> y=(exp(x)+exp(-x)/2;>> taylor(y,x,5,0)結(jié)果：ans = 1+1/2*x2+1/24*x4（2）解：程序：>&g

34、t; syms a x>> s=taylor(exp(-5*x)*sin(3*x+pi/3),5,a)結(jié)果：s = exp(-5*a)*sin(3*a+1/3*pi)+(3*exp(-5*a)*cos(3*a+1/3*pi)-5*exp(-5*a)*sin(3*a+1/3*pi)*(x-a)+(8*exp(-5*a)*sin(3*a+1/3*pi)-15*exp(-5*a)*cos(3*a+1/3*pi)*(x-a)2+(33*exp(-5*a)*cos(3*a+1/3*pi)+5/3*exp(-5*a)*sin(3*a+1/3*pi)*(x-a)3+(-161/6*exp(-5*a)*sin(3*a+1/3*pi)-40*exp(-5*a)*cos(3*a+1/3*pi)*(x-a)46、求非線性方程的符號(hào)解。（1） （2）（1）解：程序：s=solve(sym('x3+a