科學(xué)計(jì)算之MATLAB

1、 在在歐美各高等學(xué)校歐美各高等學(xué)校，MatlabMatlab成為線性代數(shù)、自動(dòng)控制理論、成為線性代數(shù)、自動(dòng)控制理論、數(shù)字信號(hào)處理、時(shí)間序列分析、動(dòng)態(tài)系統(tǒng)仿真、圖像處理等數(shù)字信號(hào)處理、時(shí)間序列分析、動(dòng)態(tài)系統(tǒng)仿真、圖像處理等諸多課程的基本教學(xué)工具，成為本科生、碩士生和博士生的諸多課程的基本教學(xué)工具，成為本科生、碩士生和博士生的必須掌握的基本技能。必須掌握的基本技能。 在在設(shè)計(jì)研究單位和工業(yè)部門(mén)設(shè)計(jì)研究單位和工業(yè)部門(mén)，MatlabMatlab已被廣泛地用于研究和已被廣泛地用于研究和解決各種具體的工程問(wèn)題。解決各種具體的工程問(wèn)題。 可以預(yù)見(jiàn)，可以預(yù)見(jiàn)，MatlabMatlab將在我國(guó)將在我國(guó)科學(xué)研究和

2、工程應(yīng)用科學(xué)研究和工程應(yīng)用中發(fā)揮越來(lái)中發(fā)揮越來(lái)越大的作用。越大的作用。1.1 Matlab語(yǔ)言的作用1.2 MATLAB的發(fā)展歷史 70年代后期，Cleve Moler擔(dān)任美國(guó)墨西哥大學(xué)計(jì)算機(jī)系系主任，為了讓學(xué)生方便的調(diào)用EISPACK和LINPACK，設(shè)計(jì)了接口程序，取名MATLAB (MATrix LABoratory)，即Matrix和Laboratory的組合。 MATLAB的第一個(gè)商業(yè)化的版本是1984年推出的是3.0的DOS版本 1992年MathWorks 公司于推出了4.0版本 1994年的4.2版本擴(kuò)充了4.0版本的功能，尤其在圖形界面設(shè)計(jì)方面更提供了新的方法。 1997年推

3、出的5.0版允許了更多的數(shù)據(jù)結(jié)構(gòu)，如單元數(shù)據(jù)、多維矩陣、對(duì)象與類等，使其成為一種更方便編程的語(yǔ)言。 1999年推出的 MATLAB 5.3版在很多方面又進(jìn)一步改進(jìn)了 MATLAB 語(yǔ)言的功能。 2000年10月底推出了其全新的 MATLAB 6.0正式版(Release 12)，在核心數(shù)值算法、界面設(shè)計(jì)、外部接口、應(yīng)用桌面等諸多方面有了極大的改進(jìn)。 現(xiàn)在的MATLAB支持各種操作系統(tǒng)，它可以運(yùn)行在十幾個(gè)操作平臺(tái)上，其中比較常見(jiàn)的有基于Windows 9X/NT、OS/2、Macintosh、Sun、Unix、Linux等平臺(tái)的系統(tǒng)?，F(xiàn)在的MATLAB再也不是一個(gè)簡(jiǎn)單的矩陣實(shí)驗(yàn)室了，它已經(jīng)演變成

4、為一種具有廣泛應(yīng)用前景的全新的計(jì)算機(jī)高級(jí)編程語(yǔ)言了。其功能也越來(lái)越強(qiáng)大，會(huì)不斷根據(jù)科研需求提出新的解決方法。1.3 Matlab 特點(diǎn)1、語(yǔ)言簡(jiǎn)潔緊湊，庫(kù)函數(shù)豐富，使用方便、語(yǔ)言簡(jiǎn)潔緊湊，庫(kù)函數(shù)豐富，使用方便2、運(yùn)算符豐富、運(yùn)算符豐富 由于Matlab是用C語(yǔ)言編寫(xiě)的，Matlab提供了和C語(yǔ)言幾乎一樣多的運(yùn)算符3、結(jié)構(gòu)化編程，面向?qū)ο筇匦?、結(jié)構(gòu)化編程，面向?qū)ο筇匦?、程序設(shè)計(jì)限制寬松，自由度大、程序設(shè)計(jì)限制寬松，自由度大 無(wú)需對(duì)矩陣預(yù)定義就可以使用5、程序的可移植性很好、程序的可移植性很好 基本上不做修改就可以在各種型號(hào)的計(jì)算機(jī)和操作系統(tǒng)上運(yùn)行6、圖形功能強(qiáng)大、圖形功能強(qiáng)大 數(shù)據(jù)的可視化非常

5、簡(jiǎn)單7、功能強(qiáng)大的工具箱、功能強(qiáng)大的工具箱8、源程序的開(kāi)放性、源程序的開(kāi)放性 所有Matlab的核心文件和工具文件都是可讀可改的源文件9、缺點(diǎn)：和其他高級(jí)程序相比，程序的執(zhí)行速度較慢、缺點(diǎn)：和其他高級(jí)程序相比，程序的執(zhí)行速度較慢1.4MATLAB 工具箱列表2 MATLAB的桌面環(huán)境及入門(mén)知識(shí) 啟動(dòng)與退出MATLAB MATLAB缺省桌面 命令窗口的使用 工作空間 歷史窗口 當(dāng)前目錄窗口和搜索路徑2.1 啟動(dòng)與退出MATLAB啟動(dòng)啟動(dòng)MATLABMATLAB直接用鼠標(biāo)雙擊桌面上直接用鼠標(biāo)雙擊桌面上MATLAB7.1MATLAB7.1圖標(biāo)圖標(biāo)或或WindowsWindows桌面的桌面的“開(kāi)始開(kāi)始

6、”“所有程序所有程序”“MATLAB7.1” “MATLAB7.1” “MATLAB7.1”“MATLAB7.1”。退出退出MATLABMATLAB關(guān)閉關(guān)閉MATLABMATLAB桌面桌面在命令窗口執(zhí)行在命令窗口執(zhí)行quitquit或或exitexit命令命令 菜單欄菜單欄工具欄工具欄文件編輯窗口文件編輯窗口工作空間工作空間歷史命令歷史命令窗口窗口Start Start 菜單菜單命令窗口命令窗口命令提示符命令提示符2.2 MATLAB缺省桌面2.2.1 命令窗口的使用激活命令窗口。激活命令窗口。“” ” 與閃爍的光標(biāo)一起表明系統(tǒng)就緒，等待輸入與閃爍的光標(biāo)一起表明系統(tǒng)就緒，等待輸入。n簡(jiǎn)單簡(jiǎn)單計(jì)

7、算計(jì)算【例例2.2-12.2-1】計(jì)算計(jì)算（1 1）在）在MATLABMATLAB命令窗口輸入命令窗口輸入 以下內(nèi)容：以下內(nèi)容：(12+2(12+2* *(7-4)/32(7-4)/32（2 2）按）按【EnterEnter】鍵，指令執(zhí)行。鍵，指令執(zhí)行。（3 3）返回的計(jì)算結(jié)果：）返回的計(jì)算結(jié)果：ans= 2 常見(jiàn)通用命令命令命令 含義含義clc 清除命令窗口的顯示內(nèi)容clear 清除Matlab工作空間中保存的變量who或whos 顯示Matlab工作空間中的變量信息dir 顯示當(dāng)前工作目錄的文件和子目錄清單cd 顯示或設(shè)置當(dāng)前工作目錄type 顯示指定m文件的內(nèi)容help或doc 獲取在線

8、幫助quit或exit 關(guān)閉/推出MATALB2.2.2 工作空間查看工作空間內(nèi)存變量，可以由查看工作空間內(nèi)存變量，可以由whowho、whoswhos 。命名新變量。命名新變量。修改變量名修改變量名刪除變量刪除變量繪圖繪圖保存變量數(shù)據(jù)保存變量數(shù)據(jù)裝入數(shù)據(jù)裝入數(shù)據(jù) 2.2.3 歷史窗口 歷史窗口：歷史窗口：n首先記錄首先記錄每次啟動(dòng)時(shí)間每次啟動(dòng)時(shí)間n并記錄在命令窗口并記錄在命令窗口輸入命令輸入命令，此次運(yùn)行期間此次運(yùn)行期間，輸入的所有命令被，輸入的所有命令被記錄為一組，并以此次啟動(dòng)時(shí)間為標(biāo)志。記錄為一組，并以此次啟動(dòng)時(shí)間為標(biāo)志。 使用歷史窗口：使用歷史窗口：n可以查看命令窗口可以查看命令窗口輸

9、入過(guò)的命令或語(yǔ)句輸入過(guò)的命令或語(yǔ)句n可以選擇一條或多條命令執(zhí)行拷貝、執(zhí)行、創(chuàng)建可以選擇一條或多條命令執(zhí)行拷貝、執(zhí)行、創(chuàng)建MM文件等。文件等。 要清除歷史記錄，可以選擇要清除歷史記錄，可以選擇Edit菜單中的菜單中的Clear Command History 命命令令2.2.4 當(dāng)前目錄窗口和搜索路徑當(dāng)前目錄窗口當(dāng)前目錄窗口：指指Matlab運(yùn)行時(shí)的工作目錄。運(yùn)行時(shí)的工作目錄。l只有在當(dāng)前目錄和搜索路徑下的文件、函數(shù)才可以被運(yùn)行和調(diào)用。l如果沒(méi)有特殊指明，數(shù)據(jù)文件也將存放在當(dāng)前目錄下；l用戶可以將自己的工作目錄設(shè)置成當(dāng)前目錄，從而使得所有操作都在當(dāng)前目錄中進(jìn)行。搜索路徑搜索路徑：指：指Matla

10、b執(zhí)行過(guò)程中對(duì)變量、函數(shù)和文件進(jìn)行搜索的路徑。執(zhí)行過(guò)程中對(duì)變量、函數(shù)和文件進(jìn)行搜索的路徑。l在在File菜單中選擇菜單中選擇Set Path命令命令或在命令窗口輸入或在命令窗口輸入pathtool命令，出現(xiàn)搜索路命令，出現(xiàn)搜索路徑設(shè)置對(duì)話框：徑設(shè)置對(duì)話框：！修改完搜索路徑后，需要！修改完搜索路徑后，需要進(jìn)行保存。進(jìn)行保存。3 MATLAB 數(shù)值計(jì)算MATLAB中的變量與表達(dá)式數(shù)組構(gòu)造矩陣運(yùn)算關(guān)系邏輯運(yùn)算數(shù)據(jù)分析3.1 MATLAB中的變量與表達(dá)式變量變量MATLAB 中變量名是以字母開(kāi)頭，可以由字母、數(shù)字或下劃線組成MATLAB 區(qū)分大小寫(xiě)，他的命令和函數(shù)全是小寫(xiě)的MATLAB 變量名不超過(guò)1

11、9個(gè)字符表達(dá)式表達(dá)式MATLAB中的表達(dá)式由運(yùn)算符、函數(shù)、變量名和數(shù)字組成MATLAB語(yǔ)句由變量和表達(dá)式組成，由兩種常見(jiàn)形式：形式1：表達(dá)式。形式2：變量=表達(dá)式。若不想顯示形式1或2的運(yùn)算結(jié)果，可在表達(dá)式后加分號(hào)“；”。特殊變量3.2 數(shù)組構(gòu)造以左方括號(hào)開(kāi)始，以空格或逗號(hào)為間隔輸入元素，最后以右方括號(hào)結(jié)束X=0 1 3 5 7 9 10X=0,1,3,5,7,9,10X=0 1 3 5 7 9 10a : 步長(zhǎng) ：b a:b (默認(rèn) 步長(zhǎng)=1)x=0:0.5:2y=0:5x=0 0.5000 1.0000 1.5000 2.0000 y=0 1.0000 2.0000 3.0000 4.00

12、00 5.0000以空格或逗號(hào)分隔，指定不同的列以分號(hào)或回車分隔，指定不同的行A=2, 4, 6, 8; 1 3 5 7; 0 0 0 0; 1, 0, 1, 03.3 矩陣的運(yùn)算矩陣矩陣與標(biāo)量間的四則運(yùn)算與標(biāo)量間的四則運(yùn)算x=1 3 4; 2,6,5; 3 2,4a=2*x-2矩陣矩陣間的四則運(yùn)算間的四則運(yùn)算必須具有相同的維數(shù)，對(duì)應(yīng)元素進(jìn)行運(yùn)算 + - .* ./ .乘（*）：要求第一個(gè)數(shù)組的列數(shù)等于第二個(gè)數(shù)組的行數(shù)。a=1 3 4;2,6,5;3 2,4b=2 3 1; 4 1 2;4 5 3c=a+bd=a.*be=a./bf=a*ba./b=b.a矩陣矩陣的的冪運(yùn)算冪運(yùn)算 . 乘方（）

13、:在指數(shù)n為自然數(shù)時(shí)相當(dāng)于n次自乘。a=1 3 4;2,6,5;3 2,4b=2 3 1; 4 1 2;4 5 3c=a.2d=a2e=a.b特殊的矩陣特殊的矩陣3.4 關(guān)系邏輯運(yùn)算關(guān)系操作符關(guān)系操作符 小于 大于= 大于等于= 等于= 不等于邏輯操作符邏輯操作符& 與| 或 非運(yùn)算法則運(yùn)算法則：若關(guān)系式成立，結(jié)果為1； 若關(guān)系式不成立，結(jié)果為0.運(yùn)算法則運(yùn)算法則：若邏輯真，結(jié)果為1； 若邏輯假，結(jié)果為0.3.5數(shù)據(jù)分析常用函數(shù)4 MATLAB數(shù)據(jù)可視化 二維曲線繪圖的基本二維曲線繪圖的基本操作操作 多次疊繪和多多次疊繪和多子圖子圖 繪制圖形的輔助操作繪制圖形的輔助操作 坐標(biāo)、刻度和分

14、格坐標(biāo)、刻度和分格線控制線控制 刻度、分格線和坐標(biāo)框刻度、分格線和坐標(biāo)框 4.1 二二維曲線繪圖的基本操作維曲線繪圖的基本操作nplot指令的基本調(diào)用格式指令的基本調(diào)用格式 （1）plot(x) nx為向量時(shí)，以該元素的下標(biāo)為橫坐標(biāo)、元素值為縱坐標(biāo)繪出曲線為向量時(shí)，以該元素的下標(biāo)為橫坐標(biāo)、元素值為縱坐標(biāo)繪出曲線nx為實(shí)數(shù)二維數(shù)組時(shí)，則按列繪制每列元素值相對(duì)其下標(biāo)的曲線，曲線數(shù)等于為實(shí)數(shù)二維數(shù)組時(shí)，則按列繪制每列元素值相對(duì)其下標(biāo)的曲線，曲線數(shù)等于x數(shù)組的列數(shù)。數(shù)組的列數(shù)。nx為復(fù)數(shù)二維數(shù)組時(shí)，則按列分別以數(shù)組的實(shí)部和虛部為橫、縱坐標(biāo)繪制多條曲線為復(fù)數(shù)二維數(shù)組時(shí)，則按列分別以數(shù)組的實(shí)部和虛部為橫、

15、縱坐標(biāo)繪制多條曲線（2）plot(x, y)nx、y為同維數(shù)組時(shí)，繪制以為同維數(shù)組時(shí)，繪制以x、y元素為橫縱坐標(biāo)的曲線元素為橫縱坐標(biāo)的曲線nx為向量，為向量，y為二維數(shù)組、且其列數(shù)或行數(shù)等于為二維數(shù)組、且其列數(shù)或行數(shù)等于x的元素?cái)?shù)時(shí)，繪制多條不同顏色的曲線的元素?cái)?shù)時(shí)，繪制多條不同顏色的曲線nx為二維數(shù)組，為二維數(shù)組，y為向量時(shí)，情況與上相同，只是為向量時(shí)，情況與上相同，只是y仍為縱坐標(biāo)。仍為縱坐標(biāo)。（3）plot(x1, y1, x2, y2, )n繪制以繪制以x1為橫坐標(biāo)、為橫坐標(biāo)、y1為縱坐標(biāo)的曲線為縱坐標(biāo)的曲線1，以，以x2為橫坐標(biāo)、為橫坐標(biāo)、y2為縱坐標(biāo)的曲線為縱坐標(biāo)的曲線2，等等。，

16、等等。n其中其中x為橫坐標(biāo)，為橫坐標(biāo)，y為縱坐標(biāo)，繪制為縱坐標(biāo)，繪制y=f(x)函數(shù)曲線函數(shù)曲線。 在在0，2區(qū)間內(nèi)，繪制曲線區(qū)間內(nèi)，繪制曲線x = 0:pi/100:2*pi;y = 2*exp(-0.5*x).*sin(2*pi*x);plot(x,y)0.52sin(2)xyex繪制曲線繪制曲線t = -pi:pi/100:pi;x = t.*cos(3*t);y = t.*sin(t).2;plot(x,y)2cos(3 ),sinxtttytt 4.2多次多次疊疊繪繪和和多多子圖子圖多次疊繪多次疊繪多次調(diào)用多次調(diào)用plot命令命令在一幅圖上在一幅圖上繪制多條曲線，需要繪制多條曲線，需

17、要hold指令的配合。指令的配合。 hold on 保持當(dāng)前坐標(biāo)軸和圖形，并可以接受下一次繪制。保持當(dāng)前坐標(biāo)軸和圖形，并可以接受下一次繪制。hold off 取消當(dāng)前坐標(biāo)軸和圖形保持，這種狀態(tài)下，調(diào)用取消當(dāng)前坐標(biāo)軸和圖形保持，這種狀態(tài)下，調(diào)用plot繪制完全新繪制完全新的圖形，不保留以前的坐標(biāo)格式、曲線。的圖形，不保留以前的坐標(biāo)格式、曲線。 重繪波形重繪波形 及其包絡(luò)線。及其包絡(luò)線。t=(0:pi/100:pi);%長(zhǎng)度為長(zhǎng)度為101的時(shí)間采樣列向量的時(shí)間采樣列向量y1=sin(t)*1,-1;%包絡(luò)線函數(shù)值，是（包絡(luò)線函數(shù)值，是（101x2）的矩陣）的矩陣y2=sin(t).*sin(9*t

18、);%長(zhǎng)度為長(zhǎng)度為101的調(diào)制波列向量的調(diào)制波列向量t3=pi*(0:9)/9;y3=sin(t3).*sin(9*t3);plot(t,y1,r:)hold onplot(t,y2,b)plot(t3,y3,bo)axis(0,pi,-1,1)%控制軸的范圍控制軸的范圍 hold off)9sin()sin(tty MATLAB允許在同一圖形窗口布置幾幅獨(dú)立的子圖。具體指令：允許在同一圖形窗口布置幾幅獨(dú)立的子圖。具體指令：subplot(m, n, k) 使（使（mxn）幅子圖中第）幅子圖中第k個(gè)子圖成為當(dāng)前圖個(gè)子圖成為當(dāng)前圖subplot(postion, left, bottom, wi

19、dth, height) 在指定的位置上開(kāi)辟子圖，并成為當(dāng)前圖在指定的位置上開(kāi)辟子圖，并成為當(dāng)前圖說(shuō)明：說(shuō)明：subplot(m, n, k)的含義的含義：圖形窗口包含（：圖形窗口包含（mxn）個(gè)子圖，）個(gè)子圖，k為要指定的當(dāng)前子圖的編號(hào)。為要指定的當(dāng)前子圖的編號(hào)。其編號(hào)其編號(hào)原則：原則： 左上方為第左上方為第1子圖，然后向右向下依次排序。該指令按缺省值分割子圖區(qū)域。子圖，然后向右向下依次排序。該指令按缺省值分割子圖區(qū)域。subplot(postion, left, bottom, width, height)用于手工指定子圖位置用于手工指定子圖位置，指定位置的四元組采用，指定位置的四元組采用

20、歸一化的標(biāo)稱單位，即認(rèn)為整個(gè)圖形窗口繪圖區(qū)域的高、寬的取值范圍都是歸一化的標(biāo)稱單位，即認(rèn)為整個(gè)圖形窗口繪圖區(qū)域的高、寬的取值范圍都是0， 1，而左下角，而左下角為（為（0,0）坐標(biāo)。）坐標(biāo)。產(chǎn)生的子圖彼此獨(dú)立。所有的繪圖指令均可以在子圖中使用。產(chǎn)生的子圖彼此獨(dú)立。所有的繪圖指令均可以在子圖中使用。 t=(pi*(0:1000)/1000);y1=sin(t);y2=sin(10*t);y12=sin(t).*sin(10*t);subplot(2,2,1),plot(t,y1);axis(0,pi,-1,1)subplot(2,2,2),plot(t,y2);axis(0,pi,-1,1)su

21、bplot(position,0.2,0.05,0.6,0.45)plot(t,y12,b-,t,y1,-y1,r:);4.3繪制繪制圖形的輔助操作圖形的輔助操作 曲線線形控制符曲線線形控制符 曲線顏色控制符曲線顏色控制符 符號(hào) - : -.-含義 實(shí)線 虛線 點(diǎn)劃線 雙劃線 符號(hào) b g rcmykw含義蘭綠紅青品紅黃黑白曲線的曲線的線形控制符、顏色控制符、數(shù)據(jù)點(diǎn)形控制符線形控制符、顏色控制符、數(shù)據(jù)點(diǎn)形控制符可以組合使用可以組合使用其先后次序不影響繪圖結(jié)果其先后次序不影響繪圖結(jié)果也可以單獨(dú)使用。也可以單獨(dú)使用。 例子例子Specify Line properties t=(0:15)*2*p

22、i/15;y=sin(t);subplot(3,2,1), plot(t, y); title(Lins style is default)subplot(3,2,2), plot(t, y, o); title(Lins style is o)subplot(3,2,3), plot(t, y, k:); title(Lins style is k:)subplot(3,2,4), plot(t, y, k:*); title(Lins style is k:*)subplot(3,2,5), plot(t, y, m-d); title(Lins style is m-d)subplot(

23、3,2,6), plot(t, y, r-p); title(Lins style is r-p)4.4坐標(biāo)坐標(biāo)、刻度和分格線控制、刻度和分格線控制 例子例子繪制繪制橢圓，長(zhǎng)軸為橢圓，長(zhǎng)軸為3.25，短軸為，短軸為1.15t=0:2*pi/99:2*pi;x=1.15*cos(t); y=3.25*sin(t); % y為長(zhǎng)軸，為長(zhǎng)軸，x為短軸為短軸subplot(2,2,1); plot(x, y);axis offtitle(axis off); subplot(2,2,2); plot(x,y);axis image;title(axis image); subplot(2,2,3);

24、plot(x,y);axis equal;title(axis Equal); subplot(2,2,4); plot(x,y);axis square;title(axis Square);axis off-101-202axis image-4-2024-202axis Equal-202-4-2024axis Square4.5刻度刻度、分格線和坐標(biāo)框、分格線和坐標(biāo)框 n分格線與分格線與grid指令指令 ngrid on畫(huà)出分格線畫(huà)出分格線ngrid off不畫(huà)分格線不畫(huà)分格線nMATLAB的的缺省設(shè)置是不畫(huà)分格線；缺省設(shè)置是不畫(huà)分格線；分格線的疏密取決于坐標(biāo)刻度（改變坐標(biāo)刻度，可改分

25、格線的疏密取決于坐標(biāo)刻度（改變坐標(biāo)刻度，可改變分格線的疏密）。變分格線的疏密）。n坐標(biāo)框坐標(biāo)框nbox on控制加邊框線控制加邊框線nbox off 控制不加邊框線控制不加邊框線n刻度設(shè)置刻度設(shè)置n指令及格式：指令及格式：set(gca, xtick, xs, ytick, ys)nxs、ys可以使任何合法的實(shí)數(shù)向量，用于分別設(shè)置可以使任何合法的實(shí)數(shù)向量，用于分別設(shè)置x、y軸的刻度。軸的刻度。繪制繪制y=1-exp(-0.3*t).*cos(0.7*t) t=6*pi*(0:100)/100;y=1-exp(-0.3*t).*cos(0.7*t);tt=t(find(abs(y-1)0.05)

26、; ts=max(tt); plot(t,y,r-);grid on; axis(0,6*pi,0.6,max(y); title(y=1-exp(-alpha*t)*cos(omega*t); hold on;plot(ts,0.95,bo);hold off; set(gca,xtick,2*pi,4*pi,6*pi,ytick,0.95,1,1.05,max(y); grid on;5 Matlab程序設(shè)計(jì) Matlab命令的執(zhí)行方式 程序控制結(jié)構(gòu) 順序結(jié)構(gòu) 條件結(jié)構(gòu) 循環(huán)結(jié)構(gòu)5.1 Matlab命令的執(zhí)行方式交互式命令執(zhí)行方式（命令窗口） 逐條輸入，逐條執(zhí)行，操作簡(jiǎn)單、直觀，但速度慢，

27、執(zhí)行過(guò)程不能保留。M文件的程序執(zhí)行方式 將命令編成程序存儲(chǔ)在一個(gè)文件中（M文件），依次運(yùn)行文件中的命令，可以重復(fù)進(jìn)行。例：建立文件將變量a、b的值互換。命令文件：clear;a = 1:10;b = 11,12,13,14;15,16,17,18;c = a; a = b; b = c;ab將文件保存為exch，并在命令窗口執(zhí)行。執(zhí)行結(jié)果：a = 11 12 13 14 15 16 17 18b = 1 2 3 4 5 6 7 8 9 105.2 程序控制結(jié)構(gòu)n 順序結(jié)構(gòu)n 選擇結(jié)構(gòu)n 循環(huán)結(jié)構(gòu) 任何復(fù)雜的程序都可以由這3種基本結(jié)構(gòu)構(gòu)成。順序結(jié)構(gòu)例5.2 求一元二次方程 的根。由于Matlab

28、能進(jìn)行復(fù)數(shù)運(yùn)算，所以不需要判斷方程的判別式，可直接根據(jù)求根公式求根。程序如下：a = input(a=?);b = input(b=?);c = input(c=?);d = b*b-4*a*c;x = (-b+sqrt(d)/(2*a),(-b-sqrt(d)/(2*a);disp(x1=,num2str(x(1),x2=,num2str(x(2);程序輸出為：a=?4b=?78c=?54x1=-0.7188,x2=-18.781220axbxc選擇選擇結(jié)構(gòu)結(jié)構(gòu)根據(jù)給定的條件成立或不成立，分別執(zhí)行不同的語(yǔ)句。（if語(yǔ)句，switch語(yǔ)句和try語(yǔ)句）1. if語(yǔ)句在Matlab中，if語(yǔ)句有3種格式。(2)雙分支雙分支if語(yǔ)句語(yǔ)句語(yǔ)句格式：if 條件 語(yǔ)句組 1else 語(yǔ)句組 2end當(dāng)條件成立時(shí)，執(zhí)行語(yǔ)