matlab習題及解答

1、第1章 MATLAB概論1.1 與其他計算機語言相比較，MATLAB語言突出的特點是什么？MATLAB具有功能強大、使用方便、輸入簡捷、庫函數(shù)豐富、開放性強等特點。1.2 MATLAB系統(tǒng)由那些部分組成？MATLAB系統(tǒng)主要由開發(fā)環(huán)境、MATLAB數(shù)學函數(shù)庫、MATLAB語言、圖形功能和應用程序接口五個部分組成。1.4 MATLAB操作桌面有幾個窗口？如何使某個窗口脫離桌面成為獨立窗口？又如何將脫離出去的窗口重新放置到桌面上？在MATLAB操作桌面上有五個窗口，在每個窗口的右上角有兩個小按鈕，一個是關閉窗口的Close按鈕，一個是可以使窗口成為獨立窗口的Undock按鈕，點擊Undock按鈕就

2、可以使該窗口脫離桌面成為獨立窗口，在獨立窗口的view菜單中選擇Dock 菜單項就可以將獨立的窗口重新防止的桌面上。1.5 如何啟動M文件編輯/調試器？在操作桌面上選擇“建立新文件”或“打開文件”操作時，M文件編輯/調試器將被啟動。在命令窗口中鍵入edit命令時也可以啟動M文件編輯/調試器。1.6 存儲在工作空間中的數(shù)組能編輯嗎？如何操作？存儲在工作空間的數(shù)組可以通過數(shù)組編輯器進行編輯：在工作空間瀏覽器中雙擊要編輯的數(shù)組名打開數(shù)組編輯器，再選中要修改的數(shù)據(jù)單元，輸入修改內容即可。1.7 命令歷史窗口除了可以觀察前面鍵入的命令外，還有什么用途？命令歷史窗口除了用于查詢以前鍵入的命令外，還可以直接

3、執(zhí)行命令歷史窗口中選定的內容、將選定的內容拷貝到剪貼板中、將選定內容直接拷貝到M文件中。1.8 如何設置當前目錄和搜索路徑，在當前目錄上的文件和在搜索路徑上的文件有什么區(qū)別？當前目錄可以在當前目錄瀏覽器窗口左上方的輸入欄中設置，搜索路徑可以通過選擇操作桌面的file菜單中的Set Path菜單項來完成。在沒有特別說明的情況下，只有當前目錄和搜索路徑上的函數(shù)和文件能夠被MATLAB運行和調用，如果在當前目錄上有與搜索路徑上相同文件名的文件時則優(yōu)先執(zhí)行當前目錄上的文件，如果沒有特別說明，數(shù)據(jù)文件將存儲在當前目錄上。1.9 在MATLAB中有幾種獲得幫助的途徑？在MATLAB中有多種獲得幫助的途徑：

4、（1）幫助瀏覽器：選擇view菜單中的Help菜單項或選擇Help菜單中的MATLAB Help菜單項可以打開幫助瀏覽器；（2）help命令：在命令窗口鍵入“help” 命令可以列出幫助主題，鍵入“help 函數(shù)名”可以得到指定函數(shù)的在線幫助信息；（3）lookfor命令：在命令窗口鍵入“l(fā)ookfor 關鍵詞”可以搜索出一系列與給定關鍵詞相關的命令和函數(shù)（4）模糊查詢：輸入命令的前幾個字母，然后按Tab鍵，就可以列出所有以這幾個字母開始的命令和函數(shù)。注意：lookfor和模糊查詢查到的不是詳細信息，通常還需要在確定了具體函數(shù)名稱后用help命令顯示詳細信息。第2章 MATLAB矩陣運算基礎2

5、.1 在MATLAB中如何建立矩陣，并將其賦予變量a？>> a=5 7 3;4 9 12.2 有幾種建立矩陣的方法？各有什么優(yōu)點？可以用四種方法建立矩陣：直接輸入法，如a=2 5 7 3，優(yōu)點是輸入方法方便簡捷；通過M文件建立矩陣，該方法適用于建立尺寸較大的矩陣，并且易于修改；由函數(shù)建立，如y=sin(x)，可以由MATLAB的內部函數(shù)建立一些特殊矩陣；通過數(shù)據(jù)文件建立，該方法可以調用由其他軟件產生數(shù)據(jù)。2.3 在進行算術運算時，數(shù)組運算和矩陣運算各有什么要求？進行數(shù)組運算的兩個數(shù)組必須有相同的尺寸。進行矩陣運算的兩個矩陣必須滿足矩陣運算規(guī)則，如矩陣a與b相乘（a*b）時必須滿足a

6、的列數(shù)等于b的行數(shù)。2.4 數(shù)組運算和矩陣運算的運算符有什么區(qū)別？在加、減運算時數(shù)組運算與矩陣運算的運算符相同，乘、除和乘方運算時，在矩陣運算的運算符前加一個點即為數(shù)組運算，如a*b為矩陣乘，a.*b為數(shù)組乘。2.5 計算矩陣與之和。>> a=5 3 5;3 7 4;7 9 8;>> b=2 4 2;6 7 9;8 3 6;>> a+bans = 7 7 7 9 14 13 15 12 142.6 求的共軛轉置。>> x=4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i; >> xan

7、s = 4.0000 - 8.0000i 3.0000 - 2.0000i 3.0000 - 5.0000i 7.0000 + 6.0000i 2.0000 + 7.0000i 9.0000 - 4.0000i 1.0000 - 4.0000i 3.0000 + 9.0000i 7.0000 + 5.0000i 4.0000 - 4.0000i2.7 計算與的數(shù)組乘積。>> a=6 9 3;2 7 5;>> b=2 4 1;4 6 8;>> a.*bans = 12 36 3 8 42 402.8 “左除”與“右除”有什么區(qū)別？在通常情況下，左除x=ab是a

8、*x=b的解，右除x=b/a是x*a=b的解，一般情況下，ab¹b/a。2.9 對于，如果，求解X。>> A=4 9 2;7 6 4;3 5 7;>> B=37 26 28;>> X=ABX = -0.5118 4.0427 1.33182.10 已知：，分別計算a的數(shù)組平方和矩陣平方，并觀察其結果。>> a=1 2 3;4 5 6;7 8 9;>> a.2ans = 1 4 9 16 25 36 49 64 81>> a2ans = 30 36 42 66 81 96 102 126 1502.11 ，觀察a與

9、b之間的六種關系運算的結果。>> a=1 2 3;4 5 6;>> b=8 7 4;3 6 2;>> a>bans = 0 1 0 1 0 1>> a>=bans = 0 1 0 1 0 1>> a<bans = 1 0 1 0 1 0>> a<=bans = 1 0 1 0 1 0>> a=bans = 0 0 0 0 0 0>> a=bans = 1 1 1 1 1 12.12 ，在進行邏輯運算時，a相當于什么樣的邏輯量。相當于a=1 1 0 1 1。2.13 在sin(

10、x)運算中，x是角度還是弧度？在sin(x)運算中，x是弧度，MATLAB規(guī)定所有的三角函數(shù)運算都是按弧度進行運算。2.14 角度，求x的正弦、余弦、正切和余切。>> x=30 45 60;>> x1=x/180*pi;>> sin(x1)ans = 0.5000 0.7071 0.8660>> cos(x1)ans = 0.8660 0.7071 0.5000>> tan(x1)ans = 0.5774 1.0000 1.7321>> cot(x1)ans = 1.7321 1.0000 0.57742.15 用四舍五入

11、的方法將數(shù)組2.4568 6.3982 3.9375 8.5042取整。>> b=2.4568 6.3982 3.9375 8.5042;>> round(b)ans = 2 6 4 92.16 矩陣，分別對a進行特征值分解、奇異值分解、LU分解、QR分解及Chollesky 分解。>> v,d=eig(a,b)v = -0.4330 -0.2543 -0.1744 -0.5657 0.9660 -0.6091 -0.7018 0.0472 0.7736d = 13.5482 0 0 0 4.8303 0 0 0 3.6216>> a=9 1 2

12、;5 6 3;8 2 7;>> u,s,v=svd(a)u = -0.5601 0.5320 -0.6350 -0.4762 -0.8340 -0.2788 -0.6779 0.1462 0.7204s = 15.5234 0 0 0 4.5648 0 0 0 3.3446v = -0.8275 0.3917 -0.4023 -0.3075 -0.9156 -0.2592 -0.4699 -0.0907 0.8781>> l,u=lu(a)l = 1.0000 0 0 0.5556 1.0000 0 0.8889 0.2041 1.0000u = 9.0000 1.00

13、00 2.0000 0 5.4444 1.8889 0 0 4.8367>> q,r=qr(a)q = -0.6903 0.3969 -0.6050 -0.3835 -0.9097 -0.1592 -0.6136 0.1221 0.7801r = -13.0384 -4.2183 -6.8260 0 -4.8172 -1.0807 0 0 3.7733>> c=chol(a)c = 3.0000 0.3333 0.6667 0 2.4267 1.1447 0 0 2.29032.17 將矩陣、和組合成兩個新矩陣：（1）組合成一個4´3的矩陣，第一列為按列順序排

14、列的a矩陣元素，第二列為按列順序排列的b矩陣元素，第三列為按列順序排列的c矩陣元素，即 （2）按照a、b、c的列順序組合成一個行矢量，即 >> a=4 2;5 7; >> b=7 1;8 3; >> c=5 9;6 2;% (1) >> d=a(:) b(:) c(:) d = 4 7 5 5 8 6 2 1 9 7 3 2% (2) >> e=a(:);b(:);c(:)' e = 4 5 2 7 7 8 1 3 5 6 9 2 或利用（1）中產生的d >> e=reshape(d,1,12) ans = 4 5

15、 2 7 7 8 1 3 5 6 9 2第3章 數(shù)值計算基礎3.2 求解多項式x3-7x2+2x+40=0的根。 >> r=1 -7 2 40; >> p=roots(r); -0.2151 0.4459 0.7949 0.27073.4 計算多項式乘法(x2+2x+2)(x2+5x+4)。 >> c=conv(1 2 2,1 5 4) c = 1 7 16 18 83.5 計算多項式除法(3x3+13x2+6x+8)/(x+4)。 >> d=deconv(3 13 6 8,1 4) d = 3 1 23.7 計算多項式的微分和積分。 >&

16、gt; p=4 12 14 5; >> pder=polyder(p); >> pders=poly2sym(pder) >> pint=polyint(p); >> pints=poly2sym(pint) pders = 12*x2-24*x-14 pints = x4-4*x3-7*x2+5*x3.8 解方程組。 >> a=2 9 0;3 4 11;2 2 6; >> b=13 6 6' >> x=ab x = 7.4000 -0.2000 -1.40003.11 矩陣，計算a的行列式和逆矩陣。

17、>> a=4 2 -6;7 5 4 ;3 4 9; >> ad=det(a) >> ai=inv(a) ad = -64 ai = -0.4531 0.6562 -0.5937 0.7969 -0.8437 0.9062 -0.2031 0.1562 -0.09373.12 y=sin(x)，x從0到2p，Dx=0.02p，求y的最大值、最小值、均值和標準差。 >> x=0:0.02*pi:2*pi; >> y=sin(x); >> ymax=max(y) >> ymin=min(y) >> yme

18、an=mean(y) >> ystd=std(y) ymax = 1 ymin = -1 ymean = 2.2995e-017 ystd = 0.70713.13 ，計算x的協(xié)方差、y的協(xié)方差、x與y的互協(xié)方差。 >> x=1 2 3 4 5; >> y=2 4 6 8 10; >> cx=cov(x) >> cy=cov(y) >> cxy=cov(x,y) cx = 2.5000 cy = 10 cxy = 2.5000 5.0000 5.0000 10.00003.15 有一正弦衰減數(shù)據(jù)y=sin(x).*exp(

19、-x/10)，其中x=0:pi/5:4*pi，用三次樣條法進行插值。 >> x0=0:pi/5:4*pi; >> y0=sin(x0).*exp(-x0/10); >> x=0:pi/20:4*pi; >> y=spline(x0,y0,x); >> plot(x0,y0,'or',x,y,'b')第4章 符號數(shù)學基礎4.1 創(chuàng)建符號變量有幾種方法？MATLAB提供了兩種創(chuàng)建符號變量和表達式的函數(shù)：sym和syms。sym用于創(chuàng)建一個符號變量或表達式，用法如x=sym(x) 及 f=sym(x+y+z)

20、，syms用于創(chuàng)建多個符號變量，用法如syms x y z。 f=sym(x+y+z)相當于 syms x y z f= x+y+z4.2 下面三種表示方法有什么不同的含義？（1）f=3*x2+5*x+2（2）f='3*x2+5*x+2'（3）x=sym('x') f=3*x2+5*x+2（1）f=3*x2+5*x+2表示在給定x時，將3*x2+5*x+2的數(shù)值運算結果賦值給變量f，如果沒有給定x則指示錯誤信息。（2）f='3*x2+5*x+2'表示將字符串'3*x2+5*x+2'賦值給字符變量f，沒有任何計算含義，因此也不對字符

21、串中的內容做任何分析。（3）x=sym('x') f=3*x2+5*x+2表示x是一個符號變量，因此算式f=3*x2+5*x+2就具有了符號函數(shù)的意義，f也自然成為符號變量了。4.3 用符號函數(shù)法求解方程at2+b*t+c=0。 >> r=solve('a*t2+b*t+c=0','t') r = 1/2/a*(-b+(b2-4*a*c)(1/2) 1/2/a*(-b-(b2-4*a*c)(1/2)4.5 求矩陣的行列式值、逆和特征根。 >> syms a11 a12 a21 a22; >> A=a11,a12

22、;a21,a22 >> AD=det(A) % 行列式 >> AI=inv(A) % 逆 >> AE=eig(A) % 特征值 A = a11, a12 a21, a22 AD = a11*a22-a12*a21 AI = -a22/(-a11*a22+a12*a21), a12/(-a11*a22+a12*a21) a21/(-a11*a22+a12*a21), -a11/(-a11*a22+a12*a21) AE = 1/2*a11+1/2*a22+1/2*(a112-2*a11*a22+a222+4*a12*a21)(1/2) 1/2*a11+1/2*

23、a22-1/2*(a112-2*a11*a22+a222+4*a12*a21)(1/2)4.6 因式分解： >> syms x; >> f=x4-5*x3+5*x2+5*x-6; >> factor(f) ans = (x-1)*(x-2)*(x-3)*(x+1)4.7 ，用符號微分求df/dx。 >> syms a x; >> f=a, x2, 1/x; exp(a*x), log(x), sin(x); >> df=diff(f) df = 0, 2*x, -1/x2 a*exp(a*x), 1/x, cos(x)第5

24、章 基本圖形處理功能5.1 繪制曲線，x的取值范圍為-5,5。 >> x=-5:0.2:5; >> y=x.3+x+1; >> plot(x,y)5.2 有一組測量數(shù)據(jù)滿足，t的變化范圍為010，用不同的線型和標記點畫出a=0.1、a=0.2和a=0.5三種情況下的曲線。 >> t=0:0.5:10; >> y1=exp(-0.1*t); >> y2=exp(-0.2*t); >> y3=exp(-0.5*t); >> plot(t,y1,'-ob',t,y2,':*r

25、9;,t,y3,'-.g')5.3 在5.1題結果圖中添加標題，并用箭頭線標識出各曲線a的取值。 >> title('ityrm=e-itat') >> title('ityrm=e-itat','FontSize',12) >> text(t(6),y1(6),'leftarrowitarm=0.1','FontSize',11) >> text(t(6),y2(6),'leftarrowitarm=0.2','FontSiz

26、e',11) >> text(t(6),y3(6),'leftarrowitarm=0.5','FontSize',11)5.4 在5.1題結果圖中添加標題和圖例框。 >> title('ityrm=e-itat','FontSize',12) >> legend('a=0.1','a=0.2','a=0.5')5.5表中列出了4個觀測點的6次測量數(shù)據(jù)，將數(shù)據(jù)繪制成為分組形式和堆疊形式的條形圖。第1次第2次第3次第4次第5次第6次觀測點13

27、67428觀測點2673247觀測點3972584觀測點4643274>> y=3 6 9 6;6 7 7 4;7 3 2 3;4 2 5 2;2 4 8 7;8 7 4 4;>> bar(y)>> bar(y,stack)5.6 x= 66 49 71 56 38，繪制餅圖，并將第五個切塊分離出來。 >> x=66 49 71 56 38; >> L=0 0 0 0 1; >> pie(x,L)5.7 ，當x和y的取值范圍均為-2到2時，用建立子窗口的方法在同一個圖形窗口中繪制出三維線圖、網(wǎng)線圖、表面圖和帶渲染效果的表面

28、圖。 >> x,y=meshgrid(-2:.2:2); >> z=x.*exp(-x.2-y.2); >> mesh(x,y,z) >> subplot(2,2,1), plot3(x,y,z) >> title('plot3 (x,y,z)') >> subplot(2,2,2), mesh(x,y,z) >> title('mesh (x,y,z)') >> subplot(2,2,3), surf(x,y,z) >> title('surf

29、 (x,y,z)') >> subplot(2,2,4), surf(x,y,z), shading interp >> title('surf (x,y,z), shading interp')5.8 繪制peaks函數(shù)的表面圖，用colormap函數(shù)改變預置的色圖，觀察色彩的分布情況。 >> surf(peaks(30);>> colormap(hot)>> colormap(cool)>> colormap(lines) 5.9 用sphere函數(shù)產生球表面坐標，繪制不通明網(wǎng)線圖、透明網(wǎng)線圖、

30、表面圖和帶剪孔的表面圖。 >> x,y,z=sphere(30); >> mesh(x,y,z) >> mesh(x,y,z),hidden off >> surf(x,y,z) >> z(18:30,1:5)=NaN*ones(13,5); >> surf(x,y,z)5.10 將5.9題中的帶剪孔的球形表面圖的坐標改變?yōu)檎叫?，以使球面看起來是圓的而不是橢圓的，然后關閉坐標軸的顯示。 axis square axis off第6章 M文件編程6.1 命令文件與函數(shù)文件的主要區(qū)別是什么？（1）命令文件是一系列命令的組合，

31、函數(shù)文件的第一行必須用function說明；（2）命令文件沒有輸入?yún)?shù)，也不用返回參數(shù)，函數(shù)文件可以接受輸入?yún)?shù)，也可以返回參數(shù)；（3）命令文件處理的變量為工作空間變量，函數(shù)文件處理的變量為函數(shù)內部的局部變量，也可以處理全局變量。6.2 如何定義全局變量？用關鍵字global可以把一個變量定義為全局變量，在M文件中定義全局變量時，如果在當前工作空間已經存在了相同的變量，系統(tǒng)將會給出警告，說明由于將該變量定義為全局變量，可能會使變量的值發(fā)生改變，為避免發(fā)生這種情況，應該在使用變量前先將其定義為全局變量。6.4 if語句有幾種表現(xiàn)形式？（1）if 表達式 語句體 end（2）if 表達式 語句體1

32、 else 語句體2 end（3）if 表達式1 語句體1 elseif 表達式2 語句體2 end（4）if 表達式1 語句體1 elseif 表達式2 語句體2 else 語句體3 end6.5 說明break語句和return語句的用法。break語句用于終止for循環(huán)和while循環(huán)的執(zhí)行。如果遇到break語句，則退出循環(huán)體，繼續(xù)執(zhí)行循環(huán)體外的下一行語句。在嵌套循環(huán)中，break只存在于最內層的循環(huán)中。return語句用于終止當前的命令序列，并返回到調用的函數(shù)或鍵盤，也用于終止keyboad方式，在MATLAB中，被調用的函數(shù)運行結束后會自動返回調用函數(shù)，使用return語句時將re

33、turn插入被調用函數(shù)的某一位置，根據(jù)某種條件迫使被調用函數(shù)提前結束并返回調用函數(shù)。6.7 編制一個解數(shù)論問題的函數(shù)文件：取任意整數(shù)，若是偶數(shù)，則用2除，否則乘3加1，重復此過程，直到整數(shù)變?yōu)?。 function c=collatz(n) % collatz % Classic “3n+1” Ploblem from number theory c=n; while n>1 if rem(n,2)=0 n=n/2; else n=3*n+1; end c=c n; end6.8 有一組學生的考試成績（見表），根據(jù)規(guī)定，成績在100分時為滿分，成績在9099之間時為優(yōu)秀，成績在8089分

34、之間時為良好，成績在6079分之間為及格，成績在60分以下時為不及格，編制一個根據(jù)成績劃分等級的程序。學生姓名王張劉李陳楊于黃郭趙成 績728356941008896685465 Name='王','張','劉','李','陳','楊','于','黃','郭','趙' Marks=72,83,56,94,100,88,96,68,54,65; % 劃分區(qū)域：滿分(100)，優(yōu)秀(90-99)，良好(80-89)，及格(60-79)，不及

35、格(<60)。 n=length(Marks); for i=1:n ai=89+i; bi=79+i; ci=69+i; di=59+i; end; c=d,c; % 根據(jù)學生的分數(shù)，求出相應的等級。 for i=1:n switch Marks(i) case 100 %得分為100時 Rank(i,:)=' 滿分' case a %得分在9099之間 Rank(i,:)=' 優(yōu)秀' case b %得分在8089之間 Rank(i,:)=' 良好' case c %得分在6079之間 Rank(i,:)=' 及格' o

36、therwise %得分低于60。 Rank(i,:)='不及格' end end % 將學生姓名，得分，級等信息打印出來。 disp(' ') disp('學生姓名 ',' 得分 ',' 等級'); disp('-') for i=1:10; disp(' ',Name(i),' ',num2str(Marks(i),' ',Rank(i,:); end 學生姓名 得分 等級 - 王 72 及格 張 83 良好 劉 56 不及格 李 94 優(yōu)秀 陳

37、100 滿分 楊 88 良好 于 96 優(yōu)秀 黃 68 及格 郭 54 不及格 趙 65 及格6.10 keyboard命令的作用是什么？當在函數(shù)中使用keyboard命令時，是否可以在工作空間瀏覽器中觀察和修改函數(shù)中的局部變量？退出keyboard狀態(tài)的命令是什么？編一個函數(shù)文件，在其中設置keyboard，觀察變量情況。當運行中的程序遇到keyboard命令時會暫時停止運行，將控制權交給鍵盤，這時命令窗口的提示符變?yōu)镵>>，可以在命令窗口或工作空間瀏覽器中觀察和修改變量內容，如果keyboard命令設置在函數(shù)中，則可以觀察和修改函數(shù)中的局部變量，檢查或修改完成后，鍵入retur

38、n命令，退出keyboard狀態(tài)，控制權交還給系統(tǒng)，程序將繼續(xù)運行。書上第四章習題4.1function X1,X2=matlab401(a,b,c)d=b2-4*a*c;if d>0 X1=(-b+sqrt(d)/(2*a); X2=(-b-sqrt(d)/(2*a);elseif d=0 X1=(-b)/(2*a); X2=X1;else r=(-b)/(2*a); k=(sqrt(abs(d)/(2*a); X1=r+k*i; X2=r-k*i;enddisp(X1,X2) 測試數(shù)據(jù)(1)>>matlab401(1,2,1) %兩個相等的實根(2)>>mat

39、lab401(1,-1,-2) %兩個不等的實根(3)>>matlab401(1,4,6.25) %兩個不等的虛根4.2function fenS=matlab402(n)fenZ=zeros(1,n);fenM=ones(1,n);fenZ(1)=2;fenZ(2)=3;fenM(1)=1;fenM(2)=2;fenSum=fenZ(1)/fenM(1)+fenZ(2)/fenM(2);for i=3:n fenZ(i)=fenZ(i-1)+fenZ(i-2); fenM(i)=fenM(i-1)+fenM(i-2); fenS(i)=fenZ(i)/fenM(i); fenSu

40、m=fenSum+fenS(i);enddisp(fenSum)測試數(shù)據(jù)(1)>>matlab402(3) %前3項的和(2)>>matlab402(15) %前15項的和4.3clear;clc;close all;funy=zeros(1,601);i=1;for x=-3:0.01:3 if (x>-1)&(x<=-3) y=(-x.2-4*x-3)/2; funy(i)=y; i=i+1; elseif (x>=-1)&(x<1) y=-x.2+1; funy(i)=y; i=i+1; else y=(-x.2+4*x-3

41、)/2; funy(i)=y; i=i+1; end endx=-3:0.01:3;plot(x,funy,'r*-')運行>>matlab4034.4clear;clc;close all;A=-4 -2 0 2 4;-3 -1 1 3 5;B=abs(A);C=find(B>3);%ÕÒ³öAÖоø¶ÔÖµ´óÓÚ3µÄÔªËØË

42、;ùÔÚµÄλÖÃD=;%´æ·ÅAÖоø¶ÔÖµ´óÓÚ3µÄÔªËØfor i=1:length(C) D(length(D)+1)=A(C(i);enddisp(C)disp(D)運行>>matlab404結果C=1 9 10 D=-4 4 54.5clear;clc

43、;close all;jc=1;i=1;while jc<1000 i=i+1; jc=jc*i; endn=i-1 運行>>matlab405結果n=64.6function matlab406(m,n)% ÕÒ³ö¡¾m,n¡¿·¶Î§ÄÚ²»Äܱ»3Õû³ýµÄÊý²¢&#

44、202;ä³ö%num=zeros(1,n-m);j=1;for i=m:n if mod(i,3)=0 num(j)=i; j=j+1; endenddisp(num)運行>>matlab406(100,200)4.7clear;clc;close allnumb=0;%ͳ¼Æ0µÄ¸öÊý±äÁ¿for i=1:1999 temp=num2str(i); numb=numb+length(strfind(temp,'0');enddisp('The numbers is:',num2str(numb);%Êä³&