matlab課后習(xí)題答案

1、習(xí)題二1 .如何理解“矩陣是MATLAB 最基本的數(shù)據(jù)對象”？答： 因?yàn)橄蛄靠梢钥闯墒莾H有一行或一列的矩陣，單個(gè)數(shù)據(jù)(標(biāo)量)可以看成是僅含一個(gè)元素的矩陣，故向量和單個(gè)數(shù)據(jù)都可以作為矩陣的特例來處理。因此，矩陣是MATLAB 最基本、最重要的數(shù)據(jù)對象。2 .設(shè) A 和 B 是兩個(gè)同維同大小的矩陣，問：(1)A*B 和 A.*B 的值是否相等？答： 不相等。(2)A./B 和 B.A 的值是否相等？答：相等。(3)A/B 和 BA 的值是否相等？答： 不相等。(4)A/B 和 BA 所代表的數(shù)學(xué)含義是什么？答： A/B 等效于 B 的逆右乘 A 矩陣，即 A*inv(B) ，而 BA 等效于 B

2、矩陣的逆左乘 A 矩陣，即 inv(B)*A 。3 .寫出完成下列操作的命令。(1)將矩陣 A 第 25 行中第 1, 3, 5 列元素賦給矩陣B 。答： B=A(2:5,1:2:5); 或 B=A(2:5,1 3 5)(2)刪除矩陣 A 的第 7 號元素。答： A(7)=(3)將矩陣A的每個(gè)元素值加30。答： A=A+30;(4)求矩陣A的大小和維數(shù)。答： size(A);ndims(A);(5)將向量t的 0 元素用機(jī)器零來代替。答： t(find(t=0)=eps;(6)將含有 12 個(gè)元素的向量x 轉(zhuǎn)換成 3 4 矩陣。答： reshape(x,3,4);(7)求一個(gè)字符串的 ASCI

3、I 碼。答： abs(123); 或 double(123);(8)求一個(gè)ASCII 碼所對應(yīng)的字符。答： char(49);4 .下列命令執(zhí)行后，L1 、 L2 、 L3 、 L4 的值分別是多少？A=1:9;B=10-A;.L1=A=B;L2=A3&A3&A7);答： L1 的值為0, 0,0,0,1, 0,0, 0,0L2 的值為1, 1,1,1,1, 0,0, 0,0L3 的值為0, 0,0,1,1, 1,0, 0,0L4 的值為4, 5,65 .已知23100.77804145655A3250326 9.54543.14完成下列操作：(1)取出 A 的前 3 行構(gòu)成矩陣B， 前兩列構(gòu)

4、成矩陣C， 右下角 3 2 子矩陣構(gòu)成矩陣 D， B 與 C 的乘積構(gòu)成矩陣E。答： B=A(1:3,:);C=A(:,1:2);D=A(2:4,3:4);E=B*C;(2)分別求 E=10&A25)。1001 ， E|D= 1 001111答： E=10&A=0.5 & x=1.5 & x=3.5 & x5.5y=log(abs(b+c/x);enddisp(y=,num2str(y);(2) 用 switch 語句實(shí)現(xiàn)計(jì)算：a=input(請輸入a的值：)；b=input(請輸入b的值：)；c=input(請輸入c的值：)； x=input(請輸入x的值：);switch fix(x/0.

5、5)case 1,2y=a*xA2+b*x+c;case num2cell(3:6)y=a*(sin(b)Ac)+x;case num2cell(7:10)y=log(abs(b+c/x);enddisp(y=,num2str(y);3 . 產(chǎn)生 20 個(gè)兩位隨機(jī)整數(shù)，輸出其中小于平均值的偶數(shù)。答：A=fix(10+89*rand(1,20);sum=0;for i=1:20sum=sum+A(i);endB=A(find(A v_max v_max=x;end;if x v_minv_min=x;end;enddisp(最大數(shù)為：,num2str(v_max);disp(最小數(shù)為：,num2

6、str(v_min);(2)用max函數(shù)、min函數(shù)實(shí)現(xiàn)：for i=1:5A(i)=input(請輸入第,num2str(i),數(shù)：);enddisp(最大數(shù)為：,num2str(max(A);disp(最小數(shù)為：,num2str(min(A);5 .已知：s 1 2 22 23 L 263 ,分別用循環(huán)結(jié)構(gòu)和調(diào)用MATLAB的sum函 數(shù)求s的值。答：(1)用循環(huán)結(jié)構(gòu)實(shí)現(xiàn)：s=0;fo門=0:63s=s+2Ai; end s(2)調(diào)用sum函數(shù)實(shí)現(xiàn)：s=0:63;s=2.As; sum(s)6 .當(dāng)n分別取100、1000、10000時(shí)，求下列各式的值。(1) 1 1 1 1 L ( 1)

7、n 11 L ( ln2)2 3 4n(2) 1( 7)i16i6414nL (2 2 4 4 6 6 L (2n) (2n) T3 3-5 5-5(2n 1) (2n 1)要求分別用循環(huán)結(jié)構(gòu)和向量運(yùn)算(使用sum或prod函數(shù))來實(shí)現(xiàn)。(1)用循環(huán)結(jié)構(gòu)實(shí)現(xiàn)：sum=0;for k=1:100sum=sum+(-1)A(k+1)/k;endsum使用sum函數(shù)：x=；for k=1:10000x=x, (-1)A(k+1)/k； endsum(x)(2)用循環(huán)結(jié)構(gòu)實(shí)現(xiàn)：sum=0;for k=1:100sum=sum+(-1)A(k+1)/(2*k-1);endsum使用sum函數(shù)：x=；fo

8、r k=1:100x=x, (-1)A(k+1)/(2*k-1);endsum(x)(3)用循環(huán)結(jié)構(gòu)實(shí)現(xiàn):sum=0;for k=1:100sum=sum+1/(4Ak); end sum 使用 sum 函數(shù)實(shí)現(xiàn)： x=;for k=1:100x=x, 1/(4Ak); end sum(x)(4) 用循環(huán)結(jié)構(gòu)實(shí)現(xiàn)：t=1;for k=1:100t=t*(2*k)*(2*k)/(2*k-1)*(2*k+1); endt使用 prod 函數(shù)實(shí)現(xiàn)：x=;for k=1:100x=x, (2*k)*(2*k)/(2*k-1)*(2*k+1); endprod(x)7.編寫一個(gè)函數(shù)文件，求小于任意自然數(shù)

9、n的斐波那契(Fibnacci)數(shù)列各項(xiàng)斐波那契數(shù)列定義如下：f11,n1f21,n2fnfn1fn2 ,n2答：function x=fibnacci(n)for i=1:nif i=2x(i)=1;elsex(i)=x(i-1)+x(i-2);endend編寫一個(gè)函數(shù)文件，用于求兩個(gè)矩陣的乘積和點(diǎn)乘，然后在命令文件中調(diào)用該 函數(shù)。答：函數(shù)文件myfnc.m：function x, y= myfnc(A, B)tryx=A*B;catchx=；endy=A.*B;命令文件myexe.m：A=input(請輸入矩陣A :);B=input(請輸入矩陣B:);x, y=myfnc(A, B);i

10、f length(x)=0display(兩矩陣的維數(shù)不匹配，無法進(jìn)行乘積運(yùn)算！);elsedisp(矩陣A和矩陣B的乘積為：);xenddisp(矩P$ A和矩陣B的點(diǎn)乘為：);yn先用函數(shù)的遞歸調(diào)用定義一個(gè)函數(shù)文件求im,然后調(diào)用該函數(shù)文件求i 11005010 4kk2Lk 1 k 1 k 1 k答：函數(shù)文件myfnc.m：function sum=myfnc(n, m)if n=1sum=1;elsesum= myfnc (n-1, m)+nAm;end1005010 .在命令窗口中調(diào)用myfnc.m文件，計(jì)算 kk2-：k 1 k 1 k 1 ksum=myfnc(100, 1)+ myfnc(50, 2)+m