MATLAB-軟件使用簡介

1、MATLAB 軟件使用簡介MATLAB 是一個(gè)功能強(qiáng)大的常用數(shù)學(xué)軟件, 它不但可以解決數(shù)學(xué)中的數(shù)值計(jì)算問題, 還可以解決符號演算問題, 并且能夠方便地繪出各種函數(shù)圖形。MATLAB自1984年由美國的MathWorks公司推向市場以來，歷經(jīng)十幾年的開展和競爭，現(xiàn)已成為國際最優(yōu)秀的科技應(yīng)用軟件之一。這里主要以適用于Windows操作系統(tǒng)的MATLAB5.3版本向讀者介紹MATLAB的使用命令和內(nèi)容。一、MATLAB 的進(jìn)入/退出 MATLAB 的安裝成功后, 系統(tǒng)會(huì)在Windows【開始】菜單的【程序】子菜單中參加啟動(dòng)MATLAB命令的圖標(biāo), 用鼠標(biāo)單擊它就可以啟動(dòng)MATLAB系統(tǒng),見圖2.1。

2、圖2.1 啟動(dòng)MATLAB啟動(dòng)MATLAB后, 屏幕上出現(xiàn)MATLAB命令窗口:圖2.2 MATLAB命令窗口 圖2.2的空白區(qū)域是MATLAB 的工作區(qū)命令輸入?yún)^(qū), 在此可輸入和執(zhí)行命令。退出MATLAB系統(tǒng)像關(guān)閉Word文件一樣, 只要用鼠標(biāo)點(diǎn)擊MATLAB系統(tǒng)集成界面右上角的關(guān)閉按鈕即可。二、 MATLAB 操作的本卷須知l 在MATLAB工作區(qū)輸入MATLAB命令后, 還須按下Enter鍵, MATLAB才能執(zhí)行你輸入的MATLAB命令, 否那么MATLAB不執(zhí)行你的命令。l MATLAB 是區(qū)分字母大小寫的。l 一般，每輸入一個(gè)命令并按下Enter鍵, 電腦就會(huì)顯示此次輸入的執(zhí)行結(jié)果

3、。(以下用表示回車)。如果用戶不想電腦顯示此次輸入的結(jié)果，只要在所輸入命令的后面再加上一個(gè)分號“；即可以到達(dá)目的。如： x= 2 + 3 x=5 x = 2 + 3 ; 不顯示結(jié)果5l 在MATLAB工作區(qū)如果一個(gè)表達(dá)式一行寫不下，可以用在此行結(jié)尾處鍵入三個(gè)英文句號的方法到達(dá)換行的目的。如： q=56+sin(pi)+exp(3)+(1+2+3+4+5)/sin(x) -5x+1/2-567/(x+y)l MATLAB 可以輸入字母、漢字，但是標(biāo)點(diǎn)符號必須在英文狀態(tài)下書寫。l MATLAB 中不需要專門定義變量的類型，系統(tǒng)可以自動(dòng)根據(jù)表達(dá)式的值或輸入的值來確定變量的數(shù)據(jù)類型。l 命令行與M文件

4、中的百分號“%標(biāo)明注釋。在語句行中百分號后面的語句被忽略而不被執(zhí)行，在M文件中百分號后面的語句可以用Help命令打印出來。三、MATLAB的變量與表達(dá)式l MATLAB的變量名MATLAB的變量名是用一個(gè)字母打頭，后面最多跟19個(gè)字母或數(shù)字來定義的。如x,y,ae3,d3er45都是合法的變量名。應(yīng)該注意不要用MATLAB中的內(nèi)部函數(shù)或命令名作為變量名。MATLAB中的變量名是區(qū)分大小寫字母的。如在MATLAB中，ab與 Ab表示兩個(gè)不同的變量。列出當(dāng)前工作空間中的變量命令為Who 將內(nèi)存中的當(dāng)前變量以簡單形式列出；Whos 列出當(dāng)前內(nèi)存變量的名稱、大小、類型等信息；Clear 去除內(nèi)存中的所

5、有變量與函數(shù)。l MATLAB的運(yùn)算符 數(shù)學(xué)運(yùn)算符：+加號，-減號，*乘號, 左除, / (右除) , (乘冪) 關(guān)系運(yùn)算符：< (小于), > (大于), <= (小于等于), >= (大于等于), = =(等于), = (不等于) 邏輯運(yùn)算符：&(邏輯與運(yùn)算), |( 邏輯或運(yùn)算), ( 邏輯非運(yùn)算)l MATLAB的表達(dá)式及語句表達(dá)式由運(yùn)算符、函數(shù)、變量名和數(shù)字組成的式子。MATLAB語句由變量、表達(dá)式及MATLAB命令組成，用戶輸入的語句由MATLAB 系統(tǒng)解釋運(yùn)行。MATLAB 語句的2種最常見的形式為：形式1：表達(dá)式形式2：變量=表達(dá)式 在第一種形式

6、中，表達(dá)式運(yùn)算后產(chǎn)生的結(jié)果如果為數(shù)值類型，系統(tǒng)自動(dòng)賦值給變量ans，并顯示在屏幕上。例1：用兩種形式計(jì)算算術(shù)運(yùn)算結(jié)果。解：Matlab命令為形式1：56+sin(pi)+exp(3) ans =1.5645e+004形式2：a=56+sin(pi)+exp(3) a =1.5645e+004如果在表達(dá)式的后面加“；，有 a=56+sin(pi)+exp(3)；執(zhí)行后不顯示運(yùn)算結(jié)果。例2：矩陣，對它們做簡單的關(guān)系與邏輯運(yùn)算解：Matlab命令為A=1,2;1,2; B=1,1;2,2; C=(A<B)&(A=B) C= 0 0 0 0四、MATLAB的數(shù)據(jù)顯示格式雖然在MATLAB

7、系統(tǒng)中數(shù)據(jù)的存儲(chǔ)和計(jì)算都是雙精度進(jìn)行的，但MATLAB可以利用菜單或format命令來調(diào)整數(shù)據(jù)的顯示格式。Format命令的格式和作用如下：l Format|format short 5位定點(diǎn)表示l Format long 15位定點(diǎn)表示l format short e 5位浮點(diǎn)表示l Format long e 15位浮點(diǎn)表示l Format short g 系統(tǒng)選擇5位定點(diǎn)和5位浮點(diǎn)中更好的表示l Format long g 系統(tǒng)選擇15位定點(diǎn)和15位浮點(diǎn)中更好的表示l Format rat 近似的有理數(shù)的表示l Format hex 十六進(jìn)制的表示l Format bank 用元 角分美

8、制定點(diǎn)表示l Format compact 變量之間沒有空行l(wèi) Format loose 變量之間有空行 例3：對數(shù)用五位定點(diǎn)、十五位定點(diǎn)以及有理數(shù)形式表示出來。解：Matlab命令為a=5+sin(7)format short , a a = 5.6570format rat,aa = 3117/551 format long,aa = 5.65698659871879 五、MATLAB 中的常用函數(shù)MATLAB的常用內(nèi)部函數(shù)有：表2.1 常用的三角函數(shù)函數(shù)名稱函數(shù)功能sinx函數(shù)名稱函數(shù)功能sin(x)正弦函數(shù)cosxasin(x)反正弦函數(shù)asinxcos(x)余弦函數(shù)tanxacos(

9、x)反余弦函數(shù)acosxtan(x)正切函數(shù)cotxatan(x)反正切函數(shù)atanxcot(x)余切函數(shù)cotxacot(x)反余切函數(shù)acotxsec(x)正割函數(shù)secxasec(x)反正割函數(shù)asecxsinh(x)雙曲函數(shù)sinhxasinh(x)反雙曲函數(shù)asinhx表2.2 常用的計(jì)算函數(shù)函數(shù)名稱函數(shù)功能abs(x)求變量x絕對值|x|angle(x)復(fù)數(shù)x的相角sqrt(x)求變量x的算術(shù)平方根real(x)求復(fù)數(shù)x的實(shí)部image(x)求復(fù)數(shù)x的虛部conj(x)求復(fù)數(shù)x的共軛復(fù)數(shù)round(x)四舍五入至最近整數(shù)fix(x)無論正負(fù)，舍去小數(shù)至最近整數(shù)ceil(x)參加正小

10、數(shù)至最近整數(shù)floor(x)舍去正小數(shù)至最近整數(shù)rat(x)將實(shí)數(shù)化為分?jǐn)?shù)表示rats(x)將實(shí)數(shù)化為多項(xiàng)分?jǐn)?shù)表示sign(x)符號函數(shù)rem(x，y)求x除以y的余數(shù)gcd(x,y)整數(shù)x和y的最大公因數(shù)lcm(x,y)整數(shù)x和y的最小公倍數(shù)exp(x)自然指數(shù)pow2(x)2的指數(shù)log(x)自然對數(shù)lnxlog2(x)以2為底的對數(shù)log10(x)以10為底的對數(shù)六、矩陣的操作 MATLAB的根本單位是矩陣，它是的MATLAB精髓，掌握矩陣的輸入、各種數(shù)值運(yùn)算以及矩陣函數(shù)的使用是以后能否學(xué)好MATLAB的關(guān)鍵。l 矩陣的輸入 I. 直接輸入創(chuàng)立矩陣 輸入方法是先鍵入左方括弧“，然后按行直

11、接鍵入矩陣的所有元素，最后鍵入右方括弧“。注意：整個(gè)矩陣以“和“ 作為首尾，同行的元素用“，或空格隔開，不同行的元素用“；或按Enter鍵來分隔；矩陣的元素可以為數(shù)字也可以為表達(dá)式，如果進(jìn)行的是數(shù)值計(jì)算，表達(dá)式中不可包含未知的變量。例4：直接輸入創(chuàng)立矩陣解：Matlab命令為A = 1,2,3; 4,15,60; 7,8,9 A = 1 2 3 4 15 60 7 8 9 或用Matlab命令A(yù)=1,2,34,15,667,8,9 A = 1 2 3 4 15 60 7 8 9 I. 用矩陣函數(shù)來生成矩陣MATLAB 提供了大量的函數(shù)來創(chuàng)立一些特殊的矩陣，表2.3給出MATLAB常用的矩陣函數(shù)

12、。表2.3 常用的矩陣函數(shù)函數(shù)名稱函數(shù)功能函數(shù)名稱函數(shù)功能zero(m,n)m行n列的零矩陣eig(A)求矩陣A的特征值eye(n)n階方矩陣poly(A)求矩陣A的特征多項(xiàng)式ones(m,n)m行n列的元素為1的矩陣trace(A)求矩陣A的跡rand(m,n)m行n列的隨機(jī)矩陣cond(A)求矩陣A的條件數(shù)randn(m,n)m行n列的正態(tài)隨機(jī)矩陣rref(A)求矩陣A的行最簡形magic(n)n階魔方矩陣inv(A)求矩陣A的逆矩陣hess(A)hess 矩陣det(A)求矩陣A的行列式sqrtm(A)求矩陣A的平方根expm(A)求矩陣A的指數(shù)值funm(A)按矩陣計(jì)算的函數(shù)值logm

13、(A)求矩陣A的對數(shù)值rank(A)求矩陣A的秩morm(A,1)求矩陣A的范數(shù)例5：輸入矩陣。解：Matlab命令為ones(3) %生成元素都為1的3階方陣ans = 1 1 1 1 1 1 1 1 1例6：輸入矩陣解：Matlab命令為zeros(2,5) %生成元素都為0的2行5列零矩陣ans = 0 0 0 0 0 0 0 0 0 0例7：生成3階魔方矩陣。解：Matlab命令為magic(3) ans = 8 1 6 3 5 7 4 9 2l 操作符“：的說明j:k 表示步長為1的等差數(shù)列構(gòu)成的數(shù)組：j, j+1, j+2, kj:i:k 表示步長為i的等差數(shù)列構(gòu)成的數(shù)組：j,j+

14、i,j+2*i, kA(i:j) 表示A(i)，A(i+1)，A(j)例8：操作符冒號:的應(yīng)用解：Matlab命令為1:5 %步長為1的等差數(shù)列。Ans = 1 2 3 4 51:2:7 %步長為2的等差數(shù)列。Ans = 1 3 5 78:-2:0 %步長為-2的等差、遞減數(shù)列。Ans = 8 6 4 2 0l 對矩陣元素的操作設(shè)A是一個(gè)矩陣,那么在MATLAB中有如下符號表示它的元素:A(i,j) 表示矩陣A的第i行第j列元素。A(:,j) 表示矩陣A的第j列。A(i,:) 表示矩陣A的第i行。A(:,:) 表示A的所有元素構(gòu)造2維矩陣A(:) 表示以矩陣A的所有元素按列做成的一個(gè)列矩陣。A

15、(i) 表示矩陣A(:)的第i個(gè)元素。 表示空矩陣I. 元素的抽取與賦值 例9：矩陣，抽取與修改矩陣A的一些元素.解：Matlab命令為A=1 23 56;sin(3) 7 9;log(2) 6 1 %輸入矩陣A。A = 1.0000 23.0000 56.0000 0.1411 7.0000 9.0000 0.6931 6.0000 1.0000A(2,3) %求矩陣A的第二行第三列元素。ans = 9A(4) %求矩陣A的第四個(gè)元素。ans = 23A(2:4) %取矩陣A的A(2),A(3),A(4)。ans = 0.1411 0.6931 23.0000A(1,:) %取矩陣A的第一行

16、。ans = 1 23 56A(:,3) %取矩陣A的第三列ans = 56 9 1a=A(1,3) %把矩陣A的第一行第三列元素賦值給變量a。a =56 A(2,1)=100 %把矩陣A的第二行第一列元素修改為100。A = 1.0000 23.0000 56.0000 100.0000 7.0000 9.0000 0.6931 6.0000 1.0000II. 矩陣的擴(kuò)充例10：矩陣，利用A與B生成矩陣，。解：Matlab命令為A=1,3;6,9; %輸入矩陣AC= A C(1,3)=100; %把矩陣A擴(kuò)充為1行3列矩陣C C = 1 3 1006 9 0B=1,5;0,8; %輸入矩陣

17、BD=A,B %由矩陣A與B合成矩陣DD = 1 3 1 5 6 9 0 8AA=A,zeros(2);zeros(2),B %由矩陣A與B合成分塊矩陣AAAA = 1 3 0 0 6 9 0 0 0 0 1 5 0 0 0 8III. 矩陣的局部刪除例11：矩陣,刪除矩陣A的第一行。解：Matlab命令為A=1 23 56;sin(3) 7 9;log(2) 6 1; A(1,:)= %刪除矩陣A的第一行A = 0.1411 7.0000 9.00000.6931 6.0000 1.0000l 矩陣的運(yùn)算A+B: 矩陣加法A-B: 矩陣減法A*B: 矩陣乘法AB: 矩陣的左除A/B: 矩陣的

18、右除 transpose(A)或A：A的轉(zhuǎn)置：數(shù)k乘以Adet(A)： A的行列式：rank(A)：A的秩七、數(shù)組在MATLAB中數(shù)組就是一行或者一列的矩陣，前邊介紹的對矩陣輸入、修改、保存都適用于數(shù)組，同時(shí)MATLAB還提供了一些創(chuàng)立數(shù)組的特殊指令。l 特殊數(shù)組的創(chuàng)立linspace(a,b,n) 給出區(qū)間a,b的n個(gè)等分點(diǎn)數(shù)據(jù)logspace(a,b,n) 給出區(qū)間的n個(gè)等比點(diǎn)數(shù)據(jù)，公比為。例12：linspace(0,1,6) %給出區(qū)間0,1 的6個(gè)等分點(diǎn)數(shù)據(jù)ans = 0 0.2000 0.4000 0.6000 0.8000 1.0000logspace(0,1,6) %給出區(qū)間的

19、6個(gè)等比點(diǎn)數(shù)據(jù)，公比為ans = 1.0000 1.5849 2.5119 3.9811 6.3096 10.0000l 數(shù)組運(yùn)算數(shù)組的運(yùn)算除了作為1×n的矩陣應(yīng)遵循矩陣的運(yùn)算規(guī)那么外，MATLAB中還為數(shù)組提供了一些特殊的運(yùn)算： 乘法為：.* ,左除為：. ,右除為：./ ,乘冪為：. 。設(shè)數(shù)組,，那么對應(yīng)的運(yùn)算具體為: 例13：數(shù)組運(yùn)算例題 a=1:5 %定義數(shù)組aa = 1 2 3 4 5b=3:2:11 %定義數(shù)組bb = 3 5 7 9 11a.2 %數(shù)組a的每一個(gè)元素求平方ans = 1 4 9 16 25a.*b %數(shù)組a的每一個(gè)元素乘以對應(yīng)的數(shù)組b的元素ans = 3

20、 10 21 36 55例14：計(jì)算的值。解：Matlab命令為x=-pi:pi/2:pi; %定義自變量xy=sin(x) %求自變量x的每一個(gè)元素對應(yīng)的正弦值y =-0.0000 -1.0000 0 1.0000 0.0000八、M文件M文件有兩種形式：命令文件和M函數(shù)文件。它們都是由假設(shè)干MATLAB語句或命令組成的文件。兩種文件的擴(kuò)展名都是.m。要注意的是M文件名一定以字母開頭，而且最好不要與內(nèi)置函數(shù)重名。在M文件中，當(dāng)表達(dá)式后面接分號時(shí)，表達(dá)式的計(jì)算結(jié)果雖不顯示但中間結(jié)果仍保存在內(nèi)存中。假設(shè)程序?yàn)槊钗募?，那么程序?zhí)行完以后，中間變量仍予以保存；假設(shè)程序?yàn)楹瘮?shù)文件，那么程序執(zhí)行完以后

21、，中間變量被全部刪除。l 文件的操作 為表達(dá)方便，用記號“主菜單名|子菜單名|.來指示子菜單。例如File|set path 表示單擊file主菜單后再選擇其中的子菜單set path。 MATLAB 對文件的翻開、關(guān)閉和保存等操作與Word完全類似，在此不再說明。在MATLAB中新建M文件的操作是在命令窗口中選擇File|New|M-File 見圖2.3，然后用鼠標(biāo)單擊M-File，可以翻開MATLAB自帶的“M函數(shù)與M文件編輯器見圖2.4，用戶就可以在此編輯窗口來編輯一個(gè)新的M文件了。MATLAB自帶的M函數(shù)與M文件編輯器還可以用來對已經(jīng)存在的M文件進(jìn)行編輯、存儲(chǔ)、修改和讀取。 圖2.3

22、新建M文件圖2. M函數(shù)與M文件編輯器(編輯窗口)l 命令文件 命令文件的一般形式為： <M文件名>.m如a1.m, pp.m等都是合法的M文件名。M文件有兩種運(yùn)行方式：一是在命令窗口直接寫文件名，按Enter鍵; 二是在編輯窗口翻開菜單Tools，再單擊Run。M文件保存的路徑一定要在搜索路徑上,否那么M文件不能運(yùn)行。以下例題中如果不做特別說明，都是以第一種方式運(yùn)行的。例15：用M命令文件畫出衰減振蕩曲線及其它的包絡(luò)線。的取值范圍是。解：步驟： 1.翻開MATLAB命令窗口，單擊File|New|Mfile(見圖2-3)翻開編輯窗口； 2.在編輯窗口逐行寫以下語句；t=0:pi/

23、50:4*pi;y0=exp(-t/3);y=exp(-t/3).*sin(3*t);plot(t,y,'-r',t,y0,':b',t,-y0,':b')3.保存M文件，并且保存在搜索路徑上，文件名為a1.m;4.運(yùn)行M文件。在命令窗口寫a1,并按Enter鍵，或者在編輯窗口翻開菜單Tools，在選擇Run。圖 2.5 衰減振蕩曲線與包絡(luò)l M函數(shù)文件 M函數(shù)文件的一般形式為： function <因變量>=<函數(shù)名>(<自變量>) M函數(shù)文件可以有多個(gè)因變量和多個(gè)自變量，當(dāng)有多個(gè)因變量時(shí)用括起來。為了更好的

24、理解函數(shù)文件，請看下例：例16：設(shè)可逆方陣為A，編寫同時(shí)求 的M函數(shù)文件。解：步驟： 1.翻開MATLAB命令窗口，單擊File|New|Mfile(見圖2-3)翻開編輯窗口； 2.在編輯窗口逐行寫以下語句； function da,a2,inva,traa=comp4(x)%M函數(shù)文件comp4.m同時(shí)求矩陣x的四個(gè)值%da為矩陣x的行列式%a2為矩陣x的平方%inva為矩陣x的逆矩陣%traa為矩陣x的轉(zhuǎn)置da=det(x)a2=x2inva=inv(x)traa=x'3.保存M函數(shù)文件，并且保存在搜索路徑上，文件名為comp4.m;4 命令窗口執(zhí)行以下語句：A=1,2;5,8；

25、%輸入矩陣A。comp4(A) %調(diào)用comp4.m函數(shù)計(jì)算矩陣A的。da = -2a2 = 11 18 45 74inva = -4.0000 1.0000 2.5000 -0.5000traa = 1 5 2 8九、程序設(shè)計(jì)語句l for循環(huán)for循環(huán)的語句為： for i=表達(dá)式可執(zhí)行語句1.可執(zhí)行語句n end說明：1表達(dá)式是一個(gè)向量，可以是m:n，m:s:n，也可以是字符串、字符串矩陣等。(2) for循環(huán)的循環(huán)體中，可以屢次嵌套for和其它的結(jié)構(gòu)體。例17：利用for循環(huán)求1100的整數(shù)之和。解:(1). 建立命令文件exam1.m。%利用for循環(huán)求1100的整數(shù)之和sum=0;

26、for i=1:100 sum=sum+i; endsum(2) 執(zhí)行命令文件exam1.m。exam1sum = 5050例18：利用for循環(huán)找出100200 之間的所有素?cái)?shù)。解: (1).建立命令文件exam2.m。 %利用for循環(huán)找出100200 之間的所有素?cái)?shù) disp('100200 之間的所有素?cái)?shù)為：') for m=100:200 k=fix(sqrt(m); %求m的算術(shù)平方根然后取整. for i=2:k+1 if rem(m,i)=0 %求整數(shù)m與i的余數(shù) break; end end if i>=k+1 disp(int2str(m) %以字符串

27、的形式顯示素?cái)?shù). end end(2) 執(zhí)行命令文件exam2.m。exam2101 103 107 109 113 127 131 137 139 149 151 157 163 167 173 179 181 191 193 197 199說明:break語句能在for循環(huán)和while循環(huán)中退出循環(huán)，繼續(xù)執(zhí)行循環(huán)后面的命令。l while循環(huán)while循環(huán)的語句為：while 表達(dá)式 循環(huán)體語句end說明：表達(dá)式一般是由邏輯運(yùn)算和關(guān)系運(yùn)算以及一般的運(yùn)算組成的表達(dá)式，以判斷循環(huán)要繼續(xù)進(jìn)行還是要停止循環(huán)。只要表達(dá)式的值非零，即為邏輯為“真“，程序就繼續(xù)循環(huán)，只要表達(dá)式的值為零就停止循環(huán)。例19

28、：利用while循環(huán)來計(jì)算1!+2!+¼+50！的值。解: (1).建立命令文件exam3.m %利用while循環(huán)來計(jì)算1!+2!+.+50！的值sum=0;i=1;while i<51 prd=1; j=1; while j<=i prd=prd*j; j=j+1; end sum=sum+prd; i=i+1;enddisp(1!+2!+.+50！的和為：)sum(2) 執(zhí)行命令文件exam3.m。 exam31!+2!+.+50！的和為：sum = 2.5613e+018l if-else-end分支此分支結(jié)構(gòu)有3種形式：（1） if 表達(dá)式執(zhí)行語句 end 功能

29、:如果表達(dá)式的值為真,就執(zhí)行語句,否那么執(zhí)行end后面的語句.（2） if 表達(dá)式執(zhí)行語句1else 執(zhí)行語句2 end 功能:如果表達(dá)式的值為真,就執(zhí)行語句1,否那么執(zhí)行語句2.（3） if 表達(dá)式1執(zhí)行語句1elseif 表達(dá)式2 執(zhí)行語句2else 語句nend 功能:如果表達(dá)式1的值為真,就執(zhí)行語句1,然后跳出if執(zhí)行語句;否那么判斷表達(dá)式2,如果表達(dá)式2的值為真,就執(zhí)行語句2,然后跳出if執(zhí)行語句.否那么依此類推,一直進(jìn)行下去.如果所有的表達(dá)式的值都為假,就執(zhí)行end后面的語句.例19：編一函數(shù)計(jì)算函數(shù)值： 解:(1).建立M函數(shù)文件yx.m。function y=yx(x)if x

30、<1 y=xelseif x>=1 & x<=10 y=2*x-1elseif x>10 & x<=30 y=3*x-11else y=sin(x)+log(x)end(2).調(diào)用M函數(shù)文件計(jì)算。result=yx(0.2),yx(2),yx(30),yx(10*pi)result = 0.2000 3.0000 79.0000 3.4473l switch-case-end分支Switch語句的形式為： switch 表達(dá)式 case 常量表達(dá)式1 語句塊1 case 常量表達(dá)式2 語句塊2 case 常量表達(dá)式n，常量表達(dá)式n+1， 語句塊no

31、therwise語句塊n+1end 功能：switch語句后面的表達(dá)式可以為任何類型；每個(gè)case后面的常量表達(dá)式可以是多個(gè)，也可以是不同類型；與if語句不同的是，各個(gè)case 和otherwise 語句出現(xiàn)的先后順序不會(huì)影響程序運(yùn)行的結(jié)果。例20: 編一個(gè)轉(zhuǎn)換成績等級的函數(shù)文件，其中成績等級轉(zhuǎn)換標(biāo)準(zhǔn)為考試成績分?jǐn)?shù)在90，100分顯示優(yōu)秀；在80，90)分顯示良好；在60，80)分顯示及格；在0，60)分顯示不及格。解:(1).建立M函數(shù)文件ff.m function result=ff(x) n=fix(x/10); switch n case 9,10 disp('優(yōu)秀')

32、 case 8 disp('良好')case 6,7 disp('及格')otherwise disp('不及格')end (2).調(diào)用M函數(shù)文件判斷99分,56分,72分各屬于哪個(gè)范圍.。ff(99) 優(yōu)秀ff(56) 不及格ff(72) 及格十、 Matlab繪圖 1plot(y)功能: 畫一條或多條折線圖。其中y是數(shù)值向量或數(shù)值矩陣。說明：當(dāng)y是數(shù)值向量時(shí)，plot(y)在坐標(biāo)系中順序的用直線段連接頂點(diǎn)i,y(i)畫出一條折線圖;當(dāng)y是數(shù)值矩陣時(shí)，Matlab為矩陣的每一列畫出一條折線，繪圖時(shí)，以矩陣y每列元素的相應(yīng)行下標(biāo)值為橫坐標(biāo)，以y的

33、元素為縱坐標(biāo)繪制的連線圖。例21：畫出向量1,3,2,9,0.5折線圖。解:MATLAB命令為 y=1,3,2,9,0.5; plot(y) 圖2.6向量式圖形2. plot(x,y)功能:畫一條或多條折線圖。其中x可以是長度為n的數(shù)值向量或是n´m的數(shù)值矩陣，y 也可以是長度為n的數(shù)值向量或是n´m的數(shù)值矩陣。說明：¬當(dāng)x ,y 都是長度為n的數(shù)值向量時(shí)，plot(x,y)在坐標(biāo)系中順序的用直線段連接頂點(diǎn)x(i),y(i)畫出一條折線圖;­當(dāng)x 是長度為n的數(shù)值向量且y是n´m的數(shù)值矩陣時(shí)，plot(x,y)用向量x分別與矩陣y的每一列匹配,

34、 在同一坐標(biāo)系中繪出m條不同顏色的折線圖;®當(dāng)x 和y都是n´m的數(shù)值矩陣時(shí)，plot(x,y)分別用矩陣x的第i列與矩陣y的第i列匹配,在同一坐標(biāo)系中繪出m條不同顏色的折線圖。注: plot(x,y)命令可以用來畫通常的函數(shù)f(x)圖形, 此時(shí)向量x常用命令x=a:h:b的形式獲得f(x)函數(shù)在繪圖區(qū)間a,b上的自變量點(diǎn)向量數(shù)據(jù),對應(yīng)的函數(shù)向量值取為y= f(x)。步長h可以任意選取，一般,步長越小,曲線越光滑,但是步長太小,會(huì)增加計(jì)算量,運(yùn)算速度要降低。通常步長h取為0.1可以到達(dá)較好的繪圖效果。如果想在圖形中標(biāo)出網(wǎng)格線，用命令：plot(x,y)，grid on即可。

35、例22：畫出函數(shù)y = sin x2 在-5 £ x £ 5 的圖形。 解: Matlab 命令: x=-5:.1:5; %取繪圖橫坐標(biāo)向量點(diǎn)xy=sin(x.2); plot(x,y),grid on 圖2.7曲線y=sinx2例23：畫出橢圓的曲線圖。解:對于這種情形，首先把它寫成參數(shù)方程 。Matlab 命令: t=0:pi/50:2*pi; x=5*cos(t); y=2*sin(t); plot(x,y) ,grid on圖2.8 橢圓3. plot(x1,y1,x2,y2,x3,y3)功能:在同一圖形窗口畫出多條不同顏色曲線，曲線關(guān)系為 。例24：在同一圖形窗口

36、畫出三個(gè)函數(shù)y = cos 2x ,y = x 2 ,y = x的圖形,-2£ x £ 2 。解:Matlab命令：x=-2:.1:2; plot(x,cos(2*x),x,x.2,x,x) legend('cos(2x)','x2','x') 圖2.9 例5的繪圖結(jié)果4ezplot(F,xmin,xmax)功能:畫出符號函數(shù)F在區(qū)間xmin,xmax內(nèi)的圖像說明: F是只含有一個(gè)變量的函數(shù)。如果區(qū)間xmin,xmax缺省,默認(rèn)區(qū)間為-2pi,2pi例25：繪制在間的圖形。Matlab命令：syms tezplot(2/3*exp(-t/2)*cos(3/2*t),0,4*pi)圖2.10 符號函數(shù)的圖形5二維特殊圖形除了plot指令外，Matlab還提供了許多其它的二維繪圖指令，這些指令大大擴(kuò)充了Matlab的曲線