Matlab試驗(yàn)三答案

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

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

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

4、作，也就是只能訪問函數(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)語句實(shí)現(xiàn)求1100的和編程：sum=0;fori=1:100sum=sum+i;endsum運(yùn)行結(jié)果：>>ljiasum=50503、用for循環(huán)語句實(shí)現(xiàn)編寫一個(gè)求階乘的函數(shù)文件。編程：functiony=jcheng(n)y=1;fori=1:ny=y*i;endy運(yùn)行結(jié)果

5、：假設(shè)n=5>>jcheng(5)y=120ans=1204、找到一個(gè)n!>10100的值(利用上題的階乘函數(shù)文件)編程：i=1;while(1)y=jcheng(i);ify>10A100break;endi=i+1;endi運(yùn)行結(jié)果：>>zhaony=24y=120y=720y=5040y=40320y=3628803628800y=39916800y=479001600y=6.2270e+009y=8.7178e+010y=1.3077e+012y=2.0923e+013y=3.5569e+014y=1.2165e+017y=2.4329e+018y=

6、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+0401.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.586

7、2e+059y=1.2414e+061y=6.0828e+062y=3.0414e+064y=1.5511e+066y=8.0658e+067y=4.2749e+069y=1.2696e+073y=7.1100e+074y=4.0527e+076y=2.3506e+078y=1.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=5、編寫一個(gè)乘法表（9X9）編程：fo

8、ri=1:9A=；forj=1:i;A=A,sprintf('%dx%d=%d',j,i,i*j);enddisp(A)end運(yùn)行結(jié)果：>>cfb1x1=11x2=22x2=41x3=32x3=63x3=91x4=42x4=83x4=124x4=161x5=52x5=103x5=154x5=205x5=251x6=62x6=123x6=184x6=245x6=306x6=361x7=72x7=143x7=214x7=285x7=356x7=427x7=491x8=82x8=163x8=244x8=325x8=406x8=487x8=568x8=641x9=92x9=

9、183x9=274x9=365x9=456x9=547x9=638x9=729x9=816、判斷y向量的元素屬于3n、3n+1或3n+2,設(shè)y=3,5,7,9,11,16。編程：假設(shè)：將3n放進(jìn)數(shù)組A,將3n+1放進(jìn)數(shù)組B,將3n+2放進(jìn)數(shù)組C。y=35791116;A=;B=;C=;j=1;k=1;m=1;fori=1:6x=rem(y(i),3);ifx=0A(j)=y(i);j=j+1;endifx=1B(k)=y(i);k=k+1;endifx=2C(m)=y(i);m=m+1;endenddisp(A)disp(B)disp(C)運(yùn)行結(jié)果：>>yush16117、已知一維

10、數(shù)組A=2,4,5,8,10、B=4,9,6,7,4,用for循環(huán)語句實(shí)現(xiàn)AiBni1。求i1和函數(shù)可用sum().編程：A=245810;B=49674;C=;y=0;fori=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：functiony=ywjz(A)x=0;n=numel(A);fori=1:nx=x+A(i);endv=x/n;運(yùn)行結(jié)果：>>y=ywjz(12345)y=3編程2：functiony=ewjz(A

11、)x=0;n=numel(A);fori=1:nx=x+A(i);endy=x/n;運(yùn)行結(jié)果：>>y=ewjz(23;45;67;)y=調(diào)用驗(yàn)證：A=133;557;y=ewjz(A);y>>aay=49、編寫驗(yàn)證魔方矩陣的函數(shù)文件，輸出要求如下：( 1) 如果輸入矩陣的維數(shù)小于3，輸出顯示error( 2) 如果輸入矩陣的不是方陣，輸出顯示thesizeofmatrixXmustbeN-by-Nmatrix( 3) 顯示行、列和及其對(duì)角線求和后的值，并判斷其和是否相同，若不同，顯示No,相同顯示Yes編程：functionyanzh(x)r,c=size(x);ifr

12、=cdisp('thesizeofmatrixXmustbeN-by-Nmatrix')elseifr<3disp('eroor')endendS=sum(x,1),sum(x,2)',sum(diag(x)',sum(diag(rot90(x)'ifS(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(816;357;492;)yes> >yan

13、zh(816;357;)thesizeofmatrixXmustbeN-by-Nmatrixno> >yanzh(86;35;)eroorno> >yanzh(826;357;492;)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);forii=1:size(A,1)forjj=1:size(B,2)fork=1:size(A,2)C(ii,jj尸C(ii,j

14、j)+A(ii,k)*B(k,jj);endendendC輸入矩陣A:12;34輸入矩陣B:56;7819224350A.*B>>A=input('請(qǐng)輸入矩陣A:');B=input('請(qǐng)輸入矩陣B:');ifsize(A)=size(B)m,n=size(A);fori=1:mforj=1:nC(i,j尸A(i,j)*B(i,j);endenddisp('輸出A.*B:')Celsedisp('無效')end請(qǐng)輸入矩陣A:123請(qǐng)輸入矩陣B:456輸出A.*B:C=41018(2)編寫程序?qū)崿F(xiàn)求矩陣的轉(zhuǎn)置>>A=input(