Matlab實(shí)驗(yàn)三答案

1、實(shí)驗(yàn)三MATLA晤言的程序設(shè)計(jì)一、實(shí)驗(yàn)?zāi)康模?、熟悉 MATLAB 程序編輯與設(shè)計(jì)環(huán)境2、掌握各種編程語(yǔ)句語(yǔ)法規(guī)則及程序設(shè)計(jì)方法3、函數(shù)文件的編寫和設(shè)計(jì)4、了解和熟悉跨空間變量傳遞和賦值二、實(shí)驗(yàn)基本知識(shí)：1、程序流程控制語(yǔ)句for 循環(huán)結(jié)構(gòu)語(yǔ)法：的門=初值：增量：終值語(yǔ)句 1語(yǔ)句 nend說(shuō)明：1. i=初值：終值，則增量為1。2初值、增量、終值可正可負(fù)，可以是整數(shù)，也可以是小數(shù)，只須符合數(shù)學(xué)邏輯。while 循環(huán)結(jié)構(gòu)語(yǔ)法： while 邏輯表達(dá)式循環(huán)體語(yǔ)句end說(shuō)明：1、while結(jié)構(gòu)依據(jù)邏輯表達(dá)式的值判斷是否執(zhí)行循環(huán)體語(yǔ)勾。若表達(dá)式的值為真，執(zhí)行循環(huán)體語(yǔ)句一次、在反復(fù)執(zhí)行時(shí)，每次都要進(jìn)行判

2、斷。若表達(dá)式的值為假，則程序執(zhí)行end之后的語(yǔ)句。2、為了避免因邏輯上的失誤，而陷入死循環(huán)，建議在循環(huán)體語(yǔ)句的適當(dāng)位置加break語(yǔ)句、以便程序能正常執(zhí)行。 （執(zhí)行循環(huán)體的次數(shù)不確定； 每一次執(zhí)行循環(huán)體后， 一定會(huì)改變while后面所跟關(guān)系式的值。）3、 while 循環(huán)也可以嵌套、其結(jié)構(gòu)如下：while邏輯表達(dá)式1循環(huán)體語(yǔ)句 1while邏輯表達(dá)式2循環(huán)體語(yǔ)句2end循環(huán)體語(yǔ)句3endif-else-end分支結(jié)構(gòu)if 表達(dá)式 1語(yǔ)句 1else if 表達(dá)式2（可選）語(yǔ)句 2else（可選）語(yǔ)句 3endend說(shuō)明：結(jié)構(gòu)是一個(gè)條件分支語(yǔ)句，若滿足表達(dá)式的條件，則往下執(zhí)行；若不滿足，則跳出

3、if 結(jié)構(gòu)。2 . else if表達(dá)式2與else為可選項(xiàng)，這兩條語(yǔ)句可依據(jù)具體情況取舍。3 .注意：每一個(gè)if都對(duì)應(yīng)一個(gè)end,即有幾個(gè)if,記就應(yīng)有幾個(gè)end。switch-case 結(jié)構(gòu)語(yǔ)法：switch表達(dá)式case量表達(dá)式1語(yǔ)句組1case量表達(dá)式2語(yǔ)句組2 otherwise語(yǔ)句組 nend說(shuō)明：1. switch后面的表達(dá)式可以是任何類型，如數(shù)字、字符串等。2.當(dāng)表達(dá)式的值與cas詬面常量表達(dá)式的值相等時(shí)，就執(zhí)行這個(gè)case后面的語(yǔ)句組如果所有的常量表達(dá)式的值都與這個(gè)表達(dá)式的值不相等時(shí)，則執(zhí)行 otherwise后的執(zhí)行語(yǔ)句。 程序流程控制指令 break、 return 、

4、pause 2、函數(shù)文件的結(jié)構(gòu)和編寫方法函數(shù)只能在函數(shù)體內(nèi)對(duì)變量進(jìn)行操作，也就是只能訪問(wèn)函數(shù)本身工作空間中的變量。M 函數(shù)文件的結(jié)構(gòu)（ 1） 函數(shù)定義行（ function ）（ 2） H1 行（函數(shù)幫助文本的第一行）（ 3） 函數(shù)幫助文本（ 4） 函數(shù)體（ 5） 注釋函數(shù)文件編寫后， 保存時(shí)， 其文件名必須與函數(shù)名相同。 注意： 函數(shù)名不要與 MATLAB 自身的函數(shù)命令相同。三 、實(shí)驗(yàn)內(nèi)容：2、用for循環(huán)語(yǔ)句實(shí)現(xiàn)求1100的和編程：sum=0;for i=1:100sum=sum+i;end sum運(yùn)行結(jié)果：>> ljiasum =50503、用for循環(huán)語(yǔ)句實(shí)現(xiàn)編寫一個(gè)求階

5、乘的函數(shù)文件。編程：function y=jcheng(n)y=1;for i=1:ny=y*i;endy運(yùn)行結(jié)果：假設(shè) n=5>> jcheng(5)y =120ans =1204、找到一個(gè)n!> 10100的值(利用上題的階乘函數(shù)文件)編程：i=1;while(1)y=jcheng(i);if y>10A100break;endi=i+1;endi運(yùn)行結(jié)果：>> zhaony =y =24y =120y =720y =5040y =40320y =3628803628800y =39916800y =479001600y =6.2270e+009y =8

6、.7178e+010y =1.3077e+012y =2.0923e+013y =3.5569e+014y =6.4024e+015y =1.2165e+017y =2.4329e+018y =5.1091e+019y =1.1240e+021y =2.5852e+022y =6.2045e+023y =1.5511e+025y =4.0329e+0263.0489e+029y =8.8418e+030y =2.6525e+032y =8.2228e+033y =2.6313e+035y =8.6833e+036y =2.9523e+038y =1.0333e+040y =3.7199e+04

7、11.1962e+0561.3764e+043y =5.2302e+044y =2.0398e+046y =8.1592e+047y =3.3453e+049y =1.4050e+051y =6.0415e+052y =2.6583e+054y =5.5026e+057y =2.5862e+059y =1.2414e+061y =6.0828e+062y =3.0414e+064y =1.5511e+066y =8.0658e+067y =4.2749e+069y =y =y =1.2696e+073y =7.1100e+074y =4.0527e+076y =2.3506e+078y =1.

8、3868e+080y =8.3210e+081y =5.0758e+083y =1.9826e+087y =1.2689e+089y =8.2477e+090y =5.4434e+092y =3.6471e+094y =2.4800e+096y =1.7112e+098y =1.1979e+100i =705、編寫一個(gè)乘法表（9X9）編程：for i=1:9A=；for j=1:i;A=A,sprintf('%dx%d=%d ',j,i,i*j);enddisp(A)end運(yùn)行結(jié)果：>> cfb1x1=11x2=2 2x2=41x3=3 2x3=6 3x3=91x4

9、=4 2x4=8 3x4=12 4x4=161x5=5 2x5=10 3x5=15 4x5=20 5x5=251x6=6 2x6=12 3x6=18 4x6=24 5x6=30 6x6=361x7=7 2x7=14 3x7=21 4x7=28 5x7=35 6x7=42 7x7=491x8=8 2x8=16 3x8=24 4x8=32 5x8=40 6x8=48 7x8=56 8x8=641x9=9 2x9=18 3x9=27 4x9=36 5x9=45 6x9=54 7x9=63 8x9=72 9x9=816、判斷y向量的元素屬于3n、3n+1或3n +2,設(shè)y = 3,5, 7, 9,11

10、,16。 編程：假設(shè)：將3n放進(jìn)數(shù)組A,將3n+1放進(jìn)數(shù)組B,將3n+2放進(jìn)數(shù)組C。y=3 5 7 9 11 16;A=;B=;C=;j=1;k=1;m=1;for i=1:6x=rem(y(i),3);if x=0A(j)=y(i);j=j+1;endif x=1B(k)=y(i);k=k+1;endif x=2C(m)=y(i);m=m+1;endenddisp(A)disp(B)disp(C)運(yùn)行結(jié)果：>> yush16117、已知一維數(shù)組 A= 2,4,5,8,10、B=4, 9,6,7, 4,用 for循環(huán)語(yǔ)句實(shí)現(xiàn)Ai Bn i 1。求i1和函數(shù)可用 sum().編程：A

11、=2 4 5 8 10;B=4 9 6 7 4;C=;y=0;for i=1:5C(i)=(A(i)*B(5-i+1); y=y+C(i);endy運(yùn)行結(jié)果：>> timu71788、(1 )編寫求解一維數(shù)組平均值的函數(shù)文件。( 2 )編寫求解二維數(shù)組平均值的函數(shù)文件。編程 1：function y=ywjz(A)x=0;n=numel(A);for i=1:nx=x+A(i);endv=x/n;運(yùn)行結(jié)果：>> y=ywjz(1 2 3 4 5) y=3編程 2：function y=ewjz(A)x=0;n=numel(A);for i=1:nx=x+A(i);end

12、y=x/n;運(yùn)行結(jié)果：>> y=ewjz(2 3;4 5;6 7;)y =調(diào)用驗(yàn)證：A=1 3 3;5 5 7;y=ewjz(A);y>> aay =49、編寫驗(yàn)證魔方矩陣的函數(shù)文件，輸出要求如下：( 1) 如果輸入矩陣的維數(shù)小于3，輸出顯示 error ( 2) 如果輸入矩陣的不是方陣，輸出顯示 the size of matrix X must be N-by-N matrix ( 3) 顯示行、列和及其對(duì)角線求和后的值，并判斷其和是否相同， 若不同，顯示No,相同顯示Yes 編程：function yanzh(x)r,c=size(x);if r=cdisp(&#

13、39;the size of matrix X must be N-by-N matrix')else if r<3disp('eroor')endendS=sum(x,1),sum(x,2)',sum(diag(x)',sum(diag(rot90(x)'if S(1)=S(2)&S(2)=S(3)&S(3)=S(4)&S(1)=S(4)disp('yes')elsedisp('no')end驗(yàn)證運(yùn)行：> > yanzh(8 1 6;3 5 7;4 9 2;)yes>

14、; > yanzh(8 1 6;3 5 7;)the size of matrix X must be N-by-N matrix no> > yanzh(8 6;3 5;)eroorno> > yanzh(8 2 6;3 5 7;4 9 2;)no.實(shí)驗(yàn)五循環(huán)結(jié)構(gòu)程序設(shè)計(jì)實(shí)驗(yàn)內(nèi)容(1)編寫程序?qū)崿F(xiàn)求兩矩陣A*B,A.*B.A*B>> A=input(' 輸入矩陣 A:');B=input('輸入矩陣 B:');C=zeros(size(A,1),size(B,2);for ii=1:size(A,1)for jj=1:

15、size(B,2)for k=1:size(A,2)C(ii,jj尸C(ii,jj)+A(ii,k)*B(k,jj);endendendC輸入矩陣A:1 2;3 4輸入矩陣B:5 6;7 819 2243 50A.*B>> A=input('請(qǐng)輸入矩陣 A:'); B=input('請(qǐng)輸入矩陣B:'); if size(A)=size(B) m,n=size(A);for i=1:mfor j=1:nC(i,j尸A(i,j)*B(i,j);endenddisp(' 輸出 A.*B:') Celsedisp(' 無(wú)效')end請(qǐng)輸入矩陣A:1 2 3請(qǐng)輸入矩陣B:4 5 6輸出A.*B:C =4 10 18(2)編寫程序?qū)崿F(xiàn)求矩陣的轉(zhuǎn)置>> A=i