matlab課后習題答案1到6章

1、習題二1. 如何理解“矩陣是 MATLABR 根本的數據對象答：由于向量可以看成是僅有一行 或一列的矩陣,單個數據(標量) 可以看成是僅含一個元素的矩陣, 故向量和單個數據都可以作為矩陣 的特例來處理.因此,矩陣是MATLA撮根本、 最重要的數據對象.2. 設A和B是兩個同維同大小的 矩陣,問：(1) A*B和A.*B的值是否相 等答：不相等.(2) A./B和 B.A 的值是否 相等答：相等.(3) A/B和BA的值是否相等答：不相等.(4) A/B和BA所代表的數 學含義是什么答：A/B等效于B的逆右乘A 矩陣,即 A*inv(B),而 BA 等效于B矩陣的逆左乘A矩 陣,即 inv(B)

2、*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的每個元素值 加30.答：A=A+30;(4) 求矩陣A的大小和維 數.答：size(A);ndims(A);(5) 將向量t的0元素用E=B*C;機器零來代替.(2) 分別求 E<D E&D E|D、答：t(find(t=0)=eps;E|D和(6) 將含有12個元素的向find(A>=10&A<25).量x轉換成3 4矩 陣.答

3、：reshape(x,3,4);0 1答：E<D= 00,(7)求一個字符串的 ASCII碼.答：abs( '123'); 或 double( '123');0111E&D=01,11(8) 求一個ASCII碼所對應的字符.答：char(49);11E|D=11,4.以下命令執(zhí)行后,L1、L2、L3、11L4的值分別是多少A=1:9;B=10-A;.L1=A=B;L2=A<=5;00E|D= 10L3=A>3&A<7;00L4=find(A>3&A<7);find(A>=10&A<2

4、5)=1;答：L1 的值為0, 0, 0, 0, 1,0, 0, 0, 05.L2 的值為1, 1, 1, 1,6. 當 A=34, NaN, Inf, -Inf,1, 0, 0, 0, 0-pi, eps, 0 時,分析以下函L3 的值為0, 0, 0, 1,數的執(zhí)行結果：all(A),1, 1, 0, 0, 0any(A) , isnan(A) , isinf(A),L4的值為4, 5,切isfinite(A) .5.答：all(A)的值為0any(A)的值為123100.7780isnan(A) 的值為0, 1, 0,A41456550, 0, 0, 0325032isinf(A) 的值

5、為0, 0, 1,69.54543.141, 0, 0, 0isfinite(A)的值為1, 0,完成以下操作：0, 0, 1, 1, 1(1) 取出A的前3行構成矩7.用結構體矩陣來存儲5名學生陣B,前兩列構成矩陣的根本情況數據,每名學生的c,右下角3 2子矩數據包括學號、姓名、專業(yè)和陣構成矩陣D, B與C6門課程的成績.的乘積構成矩陣E答：student(1).id='000T;答：B=A(1:3,:);student(1).name='Tom'C=A(:,1:2);student(1).major='computD=A(2:4,3:4);er'st

6、udent(1).grade=89,78,6產生均值為1,方差為A*B,由計算可知B*A=A*B,7,90,86,85;的500個正態(tài)分布的隨A,A-1= A-1 - A是互逆.8.建立單元矩陣B并答復有關問機數.題.答:1+sqrt*randn(5,100);4.求下面線性方程組的解.1.B1,1=1;(4)產生和A同樣大小的幺B1,2='Brenden'B2,1=reshape(1:9,3,3);B2,2=12,34,2;54,21,3;4,23,67;(1)答:(2)矩陣.4x12x2 X3答:(5)ones(size(A);將矩陣A對角線的元素加30osize(B)和

7、ndims(B)的值分別是多少size(B)的值為 2, 2 .ndims(B)的值為 2.B(2)和B(4)的值分別是多少2.答:(6)答:A+eye(size(A)*30;從矩陣A提取主對角線元素,并以這些元素構成對角陣B.B=diag(diag(A);使用函數,實現(xiàn)方陣左旋90°答:A=4,2,-1;3,-1,;12,3,0;b=2;10;8;x=inv(A)*b方程組的3x112為x2 2x33x2 82101 4 7答：B(2)=25836 9B(4)=123425421342367或右旋90°的功能.例如,原矩陣為A, A左旋后得到旋后得到CB(3)=和 B3=

8、執(zhí)行后,B的值分別是多少答：當執(zhí)行B(3)=后,B=1, 1, 4, 7; 2, 5,8; 3, 6, 9, 12, 34, 2; 54,21, 3; 4, 23, 67當執(zhí)行B3=后,B=1,;1, 4, 7; 2,5, 8; 3, 6, 9, 12, 34, 2;54, 21, 3; 4, 23, 67習題三寫出完成以下操作的命令.(1)建立3階單位矩陣Ao答:A=eye(3);3.(2)建立5乂 6隨機矩陣A,其元素為100,200范圍內的隨機整數.round(100+(200-100)*rand(5,6);B,右x=6.000026.666727.3333答:1471025811369

9、125.求以下矩陣的主對角線元素、上三角陣、下三角陣、秩、范數、(1)條件數和跡.10111211237895 A1424563052123111509321654987121110B=rot90(A);C=rot90(A,-1);建立一個方陣A,求A的逆矩陣和A的行列式的值,并驗證A與A1是互逆的.答：A=rand(3)*10;B=inv(A);C=det(A);先計算 B*A,再計算(2)答:0.43 438.9 421取主對角線元素:diag(A);上三角陣:triu(A)下三角陣:tril(A);秩:rank(A);范數:sum=sum+A(i);norm(A,1); 或norm(A)

10、; 或norm(A,inf);條件數：cond(A,1); 或 2ax bx c,asincb x,end0.5 x 1.5B=A(find(A<(sum/20);1.5 x 3.5C=B(find(rem(B,2)=0);ln b - ,3.5disrxC);5.5xcond(A,2);或cond(A,inf)跡：trace(A);(2)【請參考】6. 求矩陣A的特征值和相應的特答：(1)用if語句實現(xiàn)計算：a=input('請輸入a的值：'); b=input('請輸入b的值：'); c=input('請輸入c的值：'); x=inpu

11、t('請輸入x的值：');if x>= & x<4. 輸入20個數,求其中最大數和最小數.要求分別用循環(huán)結構和調用MATLAB max函數、min函數來實現(xiàn).答：(1)用循環(huán)結構實現(xiàn)：v_max=0;0.50.252算:v_min=0;for i=1:20x=input('請輸入第',num2str(i),' 數：');if x> v_maxv_max=x;end;if x< v_minv_min=x;end;enddisp(' 最大數為：', num2str(v_max);disp(' 最

12、小數為：',num2str(v_min);(2) 用max函數、min函數實現(xiàn)：for i=1:5A(i)=input(' 請輸入第',num2str(i),'數:');enddisp(' 最大數為：',num2str(max(A);disp(' 最小數為：', num2str(min(A);5. 已知：s 1 2 22 23 L263征向量.11A 110.5 0.25答：V,D=eig(A);習題四1. 從鍵盤輸入一個4位整數,按 如下規(guī)那么加密后輸出.加密規(guī) 那么：每位數字都加上 7,然后 用和除以10的余數取代該數

13、 字；再把第一位與第三位交換,第二位與第四位交換.答：a=input('請輸入4位整數：');A=a/1000,a/100,a/10,a; A=fix(rem(A,10);A=rem(A+7,10);b=A(3)*1000+A(4)*100+A(1 )*10+A(2);disp(' 加密后的值為： ',num2str(b);2. 分別用if語句和switch語句 實現(xiàn)以下計算,其中 a、b、c 的值從鍵盤輸入.y=a*xA2+b*x+c;endif x>= & x<y=a*(sin(b)Ac)+x;endif x>= & x&l

14、t;y=log(abs(b+c/x);enddisp('y=',num2str(y);(2)用switch 語句實現(xiàn)計a=input('請輸入a的值:'); b=input('請輸入b的值:'); c=input('請輸入c的值:'); x=input('請輸入x的值:'); switch fix(x/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('

15、;y=',num2str(y);3. 產生20個兩位隨機整數,輸 出其中小于平均值的偶數.答：A=fix(10+89*rand(1,20);sum=0;for i=1:20分別用循環(huán)結構和調用MATLAB勺sum函數求s的值.答:x=x,6.運算現(xiàn).答：(1)用循環(huán)結構實現(xiàn):s=0;endfor i=0:63s=s+2Ai;end調用sum函數實現(xiàn):s=0:63;s=2.As;sum(s)當 n 分別取 100、1000、10000時,求以下各式的值1);11 1 L(1)n1123 4n111357sum(x)(2)用循環(huán)結構實現(xiàn):sum=0;for k=1:100sum=sum+(

16、-1)A(k+1)/(2*k-endsum使用sum函數:x=;for k=1:100(噱(-1)A(k+1)/(2*k-1);endsum(x)(3)用循環(huán)結構實現(xiàn):for k=1:10011111sum=sum+1/(4Ak);4 16 64LL (-)4n3 endsum(4)使用sum函數實現(xiàn)：x=;2 2 446J L(2以=1L1 3355 7(2n 1) x2P 1/(41);sum=0;end要求分別用循環(huán)結構和向量(使用sum或prod函數)來實(1)用循環(huán)結構實現(xiàn):sum=0;for k=1:100sum=sum+(-1)A(k+1)/k;endsum使用sum函數:x=;

17、for k=1:10000sum(x)(4)用循環(huán)結構實現(xiàn):t=1;for k=1:100t=t*(2*k)*(2*k)/(2*k-1)*(2*k+1);end使用prod函數實現(xiàn):x=;for k=1:100x=x,(2*k)*(2*k)/(2*k-1)*(2*k+1);endprod(x)7. 編寫一個函數文件,求小于任意 自然數n的斐波那契(Fibnacci ) 數列各項.斐波那契數列定義如下：f11,n1f21,n2fnfn 1fn 2, n 2答：function x=fibnacci(n)for i=1:nif i<=2x(i)=1;elsex(i)=x(i-1)+x(i-2

18、);endend8. 編寫一個函數文件,用于求兩 個矩陣的乘積和點乘,然后在 命令文件中調用該函數.答：函數文件：function x, y= myfnc(A,B)tryx=A*B;catchx=;endy=A.*B;命令文件：A=input('請輸入矩陣A：');B=input('請輸入矩陣B：');x, y=myfnc(A, B);if length(x)=0display(' 兩矩陣的維數不匹配,無法進行乘積運算！');disp('矩陣A和矩陣B的乘endrho=pi/3.*theta.A2;積為：)；spolar(theta,rh

19、o)x答：執(zhí)行結果為ends=1083.disp('矩陣A和矩陣B的點命令文件執(zhí)行后的結果(1)乘為：)；為：t=0:pi/100:2*pi;yx =x=cos(t);41220y=sin(t);9,先用函數的遞歸調用定義一y=z=t;n246plot3(x,y,z)個函數文件求im ,然后第五章i 11.(1) x=-10:10;(2)調用該函數文件求y=100./(1+x.A2);u=0:pi/100:2*pi;1005010 a21k k2.k 1k 1k 1 kPlot(x,y)v=0:pi/100:2*pi;x=(1+cos(u).*cos(v);(2) x=-10:10;y

20、=(1+cos(u).*sin(v);答：y=1/(2*pi)*exp(-x.A2/2);z=sin(u);函數文件：Plot(x,y)plot3(x,y,z)function sum=myfnc(n, m)if n<=1(3) ezplot('xA2+yA2=1')(3)sum=1;(4)else(4)5.sum= myfnc (n-1,t=-10:10;plot函數：m)+nAm;x=t.A2;>> x=linspace(-10,10,200);endy=5*t.A3;>> y=;在命令窗口中調用文件,計plot(x,y)>> fo

21、r x0=x1005010 1算 kk2-：k 1k 1k 1 k2.(1)if x0>0y=y,x0.A2+(1+x0).A(1/4)+5;elseif x0=0sum=myfnc(100,1)+theta=0:2*pi;y=y,0;myfnc(50, 2)+myfnc(10,-1)rho=5*cos(theta)+4;elseif x0<0polar(theta,rho)y=y,x0.A3+sqrt(1-x0)-5;10.寫出以下程序的輸出結果.ends=0;(2)enda=12,13,14;15,16,17theta=:2*pi;>> plot(x,y);18,1

22、9,20;21,22,23;rho=12./sqrt(theta);for k=apolar(theta,rho)fplot 函數：for j=1:4fplot('(x<0).*(x.A3+sqrt(1-x)-if(3) theta=:2*pi;5)+(x=0).*0+(x>0).*(x.A2+(1+xrem(k(j),2)=0rho=5./cos(theta)-7;).a(1/4)+5)',-10,10)polar(theta,rho)s=s+k(j);第六章end(4)1.elseendtheta=:2*pi;A=randn(10,5)f=inline('

23、;sqrt(x.A2+1)');(1 ) mean(A);均值(3)dx=diff(f(x,4)std(A); 標準方差i=0:10;(2) max(max(A);最大元素xi=*i;結果：min(min(A); 最小元素polyval(P,xi)dx =(3) B=sum(A,2) ;A 每行元素ans =的和 sum(B) ;A全部元素之7. (1)和5.f=inline('sin(x).A5.*sin(5*x)'(4) sort(A) ;A的每列元素按升(1);序排列建立函數文件：quad(f,0,pi)function f=fxy(u)ans =sort(A,2

24、, ' descend' ) ;A的x=u(1);y=u(2);每行元素按將序排列仁 3.*x.A2+2*x.*y+y.A2(2)2.在命令窗口中輸入以下命令：f=inline('(1+x.A2)./(1+x.A4)')(1)U,fmin=fminsearch('fxy',1,1;(2)quad(f,-1,1)X=1 4 9 16 25 36 49 64 81 100;結果：ans =Y=1:10;U =X1 = 1:100;*(3)Y1=interp1(X,Y,X1,'cubic')f=inline('x.*sin(x)

25、./(1+cos(x)3.fmin =.')；x=165 123 150 123 141;quad(f,0,pi)y=187 126 172 125 148;(2)ans =P=polyfit(x,y,3)仁 inline('-sin(x)-cos(x<2)');P =fmax=fminbnd(f,0,pi)(4)+003 *fmax =f=inline('abs(cos(x+y)');dblquad(f,0,pi,0,pi)所以它的線性擬合曲線為：p(x)=6.ans =+8433(1) x=pi/6 pi/4 pi/3;4.仁 inline(&

26、#39;sin(x).A2+cos(x)<2'8.(1)P1=0 3 2;P2=5 -1 2;P3=1);N=64;%采樣點數0 ;dx=diff(f(x,5*pi/12)/(pi/12)T=5;%采樣時間終點P=conv(conv(P1,P2),P3)可參見第157頁例題t=linspace(0,T,N); % 給出 N個采P =dx =樣時間ti (i=1:N )000y=exp(-t); %求各采樣點樣0x=pi/2時單獨計算：本值y所以 P(x)=15x5+x=pi/2;dt=t(2)-t(1);%采樣(2)仁 inline('sin(x).A2+cos(x)&l

27、t;2'周期roots(P);f=1/dt;%采樣頻ans =diff(f(x,pi)/(pi/2)率ans =Y=fft(y);%計算 y 的+0快速傅里葉變換Y-(2)F=Y(1:N/2+1);%F(k)=Y(kx=1:3;)f=f*(0:N/2)/N;%軸f從0開始plot(f,abs(F)%繪制振幅-頻率圖9.(1)矩陣求逆法：A=2 3 5;3 7 4;1 -7 1;b=10；3;5;x=inv(A)*bx =使頻率else仁 inline('(sin(x)<2).*e矩陣除法法：A=2 3 5;3 7 4;1 -7 1;b=10;3;5;x=Abx =矩陣分解法：A=2 3 5;3 7 4;1 -7 1;b=10;3;5;L,U=lu(A);x=U(Lb)x =