MATLAB m文件和函數(shù)課件

MATLAB提供了強(qiáng)大的圖形繪制功能，用戶只需制定繪圖方式，并提供充足的繪圖數(shù)據(jù)，即可以得到所需的圖形。五、繪圖功能與技巧

Matlab與系統(tǒng)仿真1MATLABm文件和函數(shù)

（一）二維圖形繪制（二）三維圖形繪制（三）圖形可視編輯工具

主要內(nèi)容Matlab與系統(tǒng)仿真2MATLABm文件和函數(shù)（1）簡(jiǎn)單二維圖形繪制（2）特殊二維圖形繪制（3）二維圖形的修飾函數(shù)（4）圖形窗口的控制與表現(xiàn)（5）二維圖形繪制流程（總結(jié)）

（一）二維圖形繪制Matlab與系統(tǒng)仿真3MATLABm文件和函數(shù)1.Plot(y)

當(dāng)y為一向量，以y的序號(hào)作為X軸，按向量y得值繪制曲線。（1）繪制簡(jiǎn)單的二維圖形2.Plot(x,y)x，y均為向量，以x為X軸，y為Y軸繪制曲線（x,y等長(zhǎng)）。3.Plot(x,y,’option’)曲線的屬性（顏色、粗細(xì)。。。）由’option’確定。Matlab與系統(tǒng)仿真4MATLABm文件和函數(shù)’option’選項(xiàng)：線的類型、標(biāo)記符號(hào)、顏色等。曲線線型曲線顏色標(biāo)記符號(hào)選項(xiàng)意義選項(xiàng)意義選項(xiàng)意義選項(xiàng)意義選項(xiàng)意義'–''––'':''–.'實(shí)線虛線點(diǎn)線點(diǎn)劃

'b''g''m''w'

藍(lán)色綠色紅紫白色'c''k''r''y'藍(lán)綠黑色紅色黃色'*''.''x''o'星號(hào)點(diǎn)號(hào)叉號(hào)圓圈'pentagram''square''diamond''hexagram'五角星正方塊菱形塊六角星線寬LineWidth例：plot(t,y2,'k--','LineWidth',5)Matlab與系統(tǒng)仿真5MATLABm文件和函數(shù)例1：

>>t=0:0.1:pi；>>plot(t)Matlab與系統(tǒng)仿真6MATLABm文件和函數(shù)例2：繪制出一個(gè)周期內(nèi)的正弦曲線*語(yǔ)句命令及繪圖結(jié)果：>>t=0:0.1:2*pi;>>y=sin(t);>>plot(t,y)Matlab與系統(tǒng)仿真7MATLABm文件和函數(shù)例3：>>plot(t,y2,':g')>>plot(t,y3,'xk')>>t=0:0.1:2*pi;>>y1=sin(t);>>y2=cos(t);>>y3=sin(t).*cos(t);>>plot(t,y1,'-')Matlab與系統(tǒng)仿真8MATLABm文件和函數(shù)1.常用調(diào)用格式（2）繪制特殊二維圖形函數(shù)名和常用調(diào)用格式意義comet(x,y)compass(x,y)feather(x,y)fill(x,y,c)hist(y,n)loglog(x,y)polar(x,y)quiver(x,y)stairs(x,y)stem(x,y)彗星狀軌跡圖羅盤(pán)圖羽毛狀圖二維填充函數(shù)直方圖對(duì)數(shù)圖極坐標(biāo)圖磁力線圖階梯圖形火柴桿圖Matlab與系統(tǒng)仿真9MATLABm文件和函數(shù)例：>>t=-pi:0.3:pi;>>y=1./(1+exp(-t));>>plot(t,y);嘗試：>>stem(t,y);>>polar(t,y);>>stairs(t,y);……Matlab與系統(tǒng)仿真10MATLABm文件和函數(shù)Matlab與系統(tǒng)仿真11MATLABm文件和函數(shù)給定函數(shù)的曲線繪制Plot()函數(shù)可以容易地繪制顯式表達(dá)式的函數(shù)，但隱函數(shù)曲線如何繪制？按常規(guī)：>>x=-sqrt(5):0.01:sqrt(5);>>y1=sqrt((5-x.^2)/3);>>y2=-sqrt((5-x.^2)/3);>>plot(x,y1,x,y2)如：x2+3y2=5

Matlab與系統(tǒng)仿真12MATLABm文件和函數(shù)ezplot()函數(shù)>>ezplot('x^2+3*y^2-5')

Matlab與系統(tǒng)仿真13MATLABm文件和函數(shù)MATLAB提供圖形函數(shù)，專門(mén)用于對(duì)由plot命令所畫(huà)圖形進(jìn)行修飾。（3）二維圖形的修飾函數(shù)函數(shù)并非圖形修飾唯一方法——菜單、右鍵快捷菜單、屬性設(shè)置頁(yè)也可以。Matlab與系統(tǒng)仿真14MATLABm文件和函數(shù)1.坐標(biāo)軸范圍的設(shè)定：axis([Xmin,Xmax,Ymin,Ymax])2.加坐標(biāo)軸名稱：xlabel(‘string’)，ylabel(‘string’)3.給圖形加標(biāo)題：title(‘string’)4.加圖例標(biāo)注：legend(‘string1’,’string2’,…)5.加網(wǎng)格：gridon,gridoff6.加普通標(biāo)注：gtext(‘string’),在彈出的圖形窗口中，選擇地點(diǎn)，單擊鼠標(biāo)。Matlab與系統(tǒng)仿真15MATLABm文件和函數(shù)Matlab與系統(tǒng)仿真16MATLABm文件和函數(shù)例：原圖>>t=0:0.1:2*pi;>>y1=sin(t);>>y2=cos(t);>>plot(t,y1,'r:',t,y2,'g*')Matlab與系統(tǒng)仿真17MATLABm文件和函數(shù)將圖形的x軸的大小范圍限定在[0,2pi]之間，y軸的大小范圍限定在[-2,2]之間；2)x、y坐標(biāo)軸分別標(biāo)注為“弧度值”、“函數(shù)值”；3)圖形標(biāo)題標(biāo)注為“正弦曲線和余弦曲線”；4)添加圖例標(biāo)注，標(biāo)注字符分別為“y1”、“y2”；5)給圖形添加網(wǎng)格線；6)在兩條曲線上分別標(biāo)注文本“y1=sin(t)”、

“y2=cos(t)”。要求：進(jìn)一步修飾例2.16的圖形，實(shí)現(xiàn)以下要求：

Matlab與系統(tǒng)仿真18MATLABm文件和函數(shù)解：

>>helplineline([起點(diǎn)橫坐標(biāo)，終點(diǎn)橫坐標(biāo)],[起點(diǎn)縱坐標(biāo)，終點(diǎn)縱坐標(biāo)])，例:line([1,2],[3,4])將畫(huà)出（1,3）到(2,4)的一條直線，而不是(1,2)到(3,4)。Matlab與系統(tǒng)仿真19MATLABm文件和函數(shù)最后結(jié)果（注意gtext指令

）：

要產(chǎn)生這個(gè)箭頭，一.可以>>inset>>arrow二.可以通過(guò)函數(shù)實(shí)現(xiàn)，詳見(jiàn)補(bǔ)充知識(shí)的文件夾，如下圖請(qǐng)大家看同學(xué)備注Matlab與系統(tǒng)仿真20MATLABm文件和函數(shù)1.創(chuàng)建或打開(kāi)圖形窗口（4）圖形窗口的控制與表現(xiàn)在命令窗口運(yùn)行繪圖指令，自動(dòng)創(chuàng)建（當(dāng)前窗口，唯一）；Figure命令，創(chuàng)建另一窗口（當(dāng)前），保留原窗口；菜單建立新窗口；Clf命令，清除當(dāng)前窗口所有內(nèi)容Matlab與系統(tǒng)仿真21MATLABm文件和函數(shù)圖形重疊（同一坐標(biāo)系中繪制多幅圖形）（1）Plot命令會(huì)清屏，但可在一個(gè)plot命令中實(shí)現(xiàn)：Plot(x,y1,’option1’,x,y2,’option2’,…)x為公共X軸，分別以y1,y2為Y軸繪制多條曲線，每條曲線的屬性由相應(yīng)的’option’確定。Matlab與系統(tǒng)仿真22MATLABm文件和函數(shù)例：>>t=0:0.1:2*pi;

>>y1=sin(t);

>>y2=cos(t);>>y3=sin(t).*cos(t);>>plot(t,y1,'-',t,y2,':g',t,y3,'xk')Matlab與系統(tǒng)仿真23MATLABm文件和函數(shù)思考：但很多實(shí)驗(yàn)情況下，實(shí)際上是對(duì)同一個(gè)變量進(jìn)行測(cè)量，并將數(shù)據(jù)進(jìn)行繪圖比較。Matlab與系統(tǒng)仿真24MATLABm文件和函數(shù)（2）holdon、holdoff指令

holdon：保留當(dāng)前窗口圖形；

holdoff：解除holdon；例：t=0:0.1:2*pi;y1=sin(t);y2=cos(t);y3=sin(t).*cos(t);plot(t,y1,'-')

holdonplot(t,y2,':g')plot(t,y3,'xk')line([0,8],[0,0])

holdoffMatlab與系統(tǒng)仿真25MATLABm文件和函數(shù)3.圖形窗口分割利用subplot()函數(shù)將整個(gè)圖形窗口分割成若干個(gè)子圖形部分，然后在每個(gè)部分用不同函數(shù)繪制不同曲線。subplot(m,n,p)：將當(dāng)前窗口分割成m行、n列，并且當(dāng)前在其中的第p個(gè)區(qū)域繪圖。繪圖區(qū)域編號(hào)原則：“從左到右，從上到下”。允許每個(gè)子窗口不同的坐標(biāo)系。Matlab與系統(tǒng)仿真26MATLABm文件和函數(shù)subplot(3,2,2)t=0:0.1:2*pi;y1=sin(t);plot(t,y1,'-')subplot(3,2,4)t=-pi:0.1:pi;y2=sin(t);plot(t,y2,'xk')例：Matlab與系統(tǒng)仿真27MATLABm文件和函數(shù)MATLAB圖形窗口控制函數(shù)（匯總）函數(shù)功能說(shuō)明figurefigure(n)每調(diào)用一次就打開(kāi)一個(gè)新的圖形窗口創(chuàng)建或打開(kāi)第n個(gè)圖形窗口，使之成為當(dāng)前窗口clf清除當(dāng)前圖形窗holdonholdoff保留當(dāng)前窗口的圖形不被后繼圖形覆蓋，可實(shí)現(xiàn)在同一坐標(biāo)系中多幅圖形的重疊解除holdon命令，一般與holdon成對(duì)使用subplot(m,n,p)將當(dāng)前繪圖窗口分割成m行、n列，并在第p個(gè)區(qū)域繪圖closecloseall關(guān)閉當(dāng)前圖形窗口關(guān)閉所有圖形窗口Matlab與系統(tǒng)仿真28MATLABm文件和函數(shù)（5）二維圖形繪制流程圖輸入圖形顯示的相關(guān)函數(shù)選擇圖形顯示區(qū)域調(diào)用相關(guān)的繪圖指令設(shè)置圖形的屬性圖形的修飾保存打印輸出Matlab與系統(tǒng)仿真29MATLABm文件和函數(shù)>>x=0:0.05:pi;>>y1=5*sin(2*x);>>y2=4*cos(4*x);>>subplot(3,1,1);>>plot(x,y1);>>axis([0pi-5.55.5]);>>gridon;>>xlabel('Time');>>ylabel('Amplitude');>>title('y=5sin(2x)');相關(guān)函數(shù)圖形顯示區(qū)域繪圖指令、屬性設(shè)置圖形修飾舉例：Matlab與系統(tǒng)仿真30MATLABm文件和函數(shù)>>subplot(3,1,2);>>plot(x,y2);>>axis([0pi-5.55.5]);>>gridon;>>xlabel('Time');>>ylabel('Amplitude');>>title('y=4cos(4x)');>>subplot(3,1,3);>>plot(x,y1,'-.or',x,y2,'-b+');>>xlabel('Time');>>ylabel('Amplitude');>>title('characteristiccurve');Matlab與系統(tǒng)仿真31MATLABm文件和函數(shù)例：可以用subplot()函數(shù)將整個(gè)圖形窗口分割成若干個(gè)子圖形部分，然后在每個(gè)部分用不同函數(shù)繪制不同曲線。>>t=-pi:0.3:pi;y=1./(1+exp(-t));>>subplot(221),plot(t,y);>>title('plot(t,y)')>>subplot(222),stem(t,y);>>title('stem(t,y)')>>subplot(223),polar(t,y);>>title('polar(t,y)')>>subplot(224),stairs(t,y);>>title('stairs(t,y)')Matlab與系統(tǒng)仿真32MATLABm文件和函數(shù)

（二）三維圖形繪制

（1）三維曲線的基本繪圖（2）三維網(wǎng)格曲面的繪制（3）圖像顯示與處理Matlab與系統(tǒng)仿真33MATLABm文件和函數(shù)調(diào)用格式

Plot3(x1,y1,z1,‘option1’,x2,y2,z2,‘option2’,…)（1）三維曲線的基本繪圖參數(shù)意義與二維繪圖命令相似。以逐點(diǎn)連折線的方式繪制三維折線圖，當(dāng)各個(gè)數(shù)據(jù)點(diǎn)的間距較小時(shí)，繪制的即三維曲線。也可用各種特殊繪圖函數(shù)。也可用各種圖形修飾函數(shù)，但要求多考慮一個(gè)Z軸方向。Matlab與系統(tǒng)仿真34MATLABm文件和函數(shù)t=0:pi/50:8*pi;x=sin(t);y=cos(t);plot3(x,y,t);stem3(x,y,t)例：三維螺旋線繪制。Matlab與系統(tǒng)仿真35MATLABm文件和函數(shù)三維網(wǎng)格曲面圖所構(gòu)成的網(wǎng)格狀表面由X-Y平面上的矩形柵格及相應(yīng)的Z坐標(biāo)構(gòu)成，相鄰點(diǎn)之間用直線連接。（2）三維網(wǎng)格曲面的繪制等高圖密度圖建筑立面。。。對(duì)于顯示大型數(shù)據(jù)矩陣或雙變量的函數(shù)是很有用。Matlab與系統(tǒng)仿真36MATLABm文件和函數(shù)即定義四邊形各個(gè)頂點(diǎn)的三維坐標(biāo)值；1.柵格數(shù)據(jù)點(diǎn)的產(chǎn)生一般先知道（x,y)，再利用函數(shù)公式計(jì)算z的坐標(biāo)；二維坐標(biāo)（x,y)是種柵格形的數(shù)據(jù)點(diǎn)，由meshgrid命令產(chǎn)生。格式

meshgrid(x,y):x,y是向量，表示坐標(biāo)范圍。Matlab與系統(tǒng)仿真37MATLABm文件和函數(shù)x=-2:0.1:2;y=x;[X,Y]=meshgrid(x,y);例Matlab與系統(tǒng)仿真38MATLABm文件和函數(shù)2.繪制三維網(wǎng)格曲面函數(shù)調(diào)用格式說(shuō)明[X,Y]=meshgrid(x,y)根據(jù)(x,y)二維坐標(biāo)數(shù)據(jù)生成x-y網(wǎng)格點(diǎn)坐標(biāo)數(shù)據(jù)其中：x,y是向量；X,Y是矩陣mesh(X,Y,Z)繪制三維網(wǎng)格曲面：通過(guò)直線連接相鄰的點(diǎn)構(gòu)成的三維曲面surf(X,Y,Z)繪制三維陰影曲面：通過(guò)小平面連接相鄰的點(diǎn)構(gòu)成的三維曲面Matlab與系統(tǒng)仿真39MATLABm文件和函數(shù)Z=X.*exp(-X.^2-Y.^2);mesh(X,Y,Z)gridon

Matlab與系統(tǒng)仿真40MATLABm文件和函數(shù)例：繪制由函數(shù)表示的曲面圖形。Matlab與系統(tǒng)仿真41MATLABm文件和函數(shù)Matlab與系統(tǒng)仿真42MATLABm文件和函數(shù)

imread()％讀取放在當(dāng)前工作目錄中的圖像文件函數(shù)，

％保存為3維數(shù)組數(shù)據(jù)；

image()％讀取保存的3維圖像數(shù)據(jù)并顯示圖像。（3）圖像顯示與處理例：W=imread(‘pic.jpg’);％圖形文件存當(dāng)前目錄

image(W)

Matlab與系統(tǒng)仿真43MATLABm文件和函數(shù)W是三維數(shù)據(jù)，由三層矩陣疊放而成，分別存放表示紅、綠、藍(lán)三原色分量。嘗試：改變圖片的色調(diào)？

1）將照片變?yōu)辄S顏色基調(diào)；2）將照片變?yōu)榈灼Ч?。Matlab與系統(tǒng)仿真44MATLABm文件和函數(shù)>>W(:,:,3)=0;>>image(W)

>>image(uint8(255-double(W)))

Matlab與系統(tǒng)仿真45MATLABm文件和函數(shù)（三）MATLAB的圖形可視化編輯工具圖形管理窗口除了簡(jiǎn)單的顯示圖形功能外，本身就是一個(gè)功能強(qiáng)大的圖形可視編輯工具，可實(shí)現(xiàn)的功能主要有：1)通用的圖形文件管理功能，如保存、打開(kāi)、新建圖形文件等；2)通用的圖形效果編輯功能，如圖形放大、縮小、旋轉(zhuǎn)、對(duì)齊等；3)圖形對(duì)象插入功能，如插入坐標(biāo)軸名稱、圖形標(biāo)題、圖例標(biāo)注、線段、文字等；4)獨(dú)立地設(shè)置窗口中各圖形對(duì)象屬性功能，如線段的類型、顏色、粗細(xì)等。

Matlab與系統(tǒng)仿真46MATLABm文件和函數(shù)（1）選擇菜單[insert]三種編輯圖形方法：t=0:0.1:2*pi;y1=sin(t);y2=cos(t);y3=sin(t).*cos(t);plot(t,y1,'-')holdonplot(t,y2,':g')plot(t,y3,'xk')holdoffMatlab與系統(tǒng)仿真47MATLABm文件和函數(shù)單擊圖形窗口工具欄中的編輯圖標(biāo)（箭頭）即可進(jìn)入圖形編輯狀態(tài)。（2）右鍵快捷菜單用戶任意選擇特性對(duì)象，出現(xiàn)選中標(biāo)志，單擊右鍵，可進(jìn)行選項(xiàng)設(shè)置。2012aMatlab與系統(tǒng)仿真48MATLABm文件和函數(shù)在方法（2）的快捷菜單中，選擇“屬性”項(xiàng)，可進(jìn)行如下設(shè)置環(huán)境：（3）“屬性”菜單項(xiàng)鼠標(biāo)在圖形中非圖形對(duì)象區(qū)域單擊，將顯示圖形修飾對(duì)話框：Matlab與系統(tǒng)仿真49MATLABm文件和函數(shù)Matlab與系統(tǒng)仿真50MATLABm文件和函數(shù)Matlab與系統(tǒng)仿真51MATLABm文件和函數(shù)第二章回顧2.1MATLAB的工作環(huán)境MATLAB程序主界面：

菜單、工具欄、命令窗口、當(dāng)前路徑瀏覽器、工作空間瀏覽器、命令歷史瀏覽器Matlab與系統(tǒng)仿真52MATLABm文件和函數(shù)MATLAB兩種編程方式比較交互式的指令行操作方式M文件編程工作方式工作過(guò)程用戶在工作窗口中按MATLAB語(yǔ)法規(guī)則輸入命令行后回車(chē)確認(rèn)，系統(tǒng)將執(zhí)行該命令并給出運(yùn)算結(jié)果。當(dāng)用戶在工作窗口中輸入M文件名并回車(chē)確認(rèn)后，系統(tǒng)將自動(dòng)搜索該文件。若該文件存在，系統(tǒng)將按M文件中語(yǔ)句所規(guī)定的計(jì)算任務(wù)以解釋方式逐一執(zhí)行語(yǔ)句，并返回運(yùn)算結(jié)果。優(yōu)點(diǎn)簡(jiǎn)便易行，交互性強(qiáng)輸入、編輯和調(diào)試、保存簡(jiǎn)便缺點(diǎn)當(dāng)要解決的問(wèn)題變得復(fù)雜后，輸入、編輯和調(diào)試?yán)щy需要在文本編輯器下編輯并保存文件，過(guò)程較復(fù)雜適用情況非常適合于對(duì)簡(jiǎn)單問(wèn)題的數(shù)學(xué)演算、結(jié)果分析及測(cè)試非常適合于大型或復(fù)雜問(wèn)題的解決文本編輯窗口：Matlab與系統(tǒng)仿真53MATLABm文件和函數(shù)2.2MATLAB語(yǔ)言的基本元素

變量：變量的命名應(yīng)遵循如下規(guī)則

賦值語(yǔ)句

矩陣的表示：一般、等差輸入、特殊、矩陣元素函數(shù)功能ones(n,m)產(chǎn)生n行m列的全1矩陣zeros(n,m)產(chǎn)生n行m列的全0矩陣rand(n,m)產(chǎn)生n行m列的在[0，1]區(qū)間均勻分布的隨機(jī)矩陣randn(n,m)產(chǎn)生n行m列的正態(tài)分布的隨機(jī)矩陣eye(n)產(chǎn)生(n×n)維的單位矩陣Matlab與系統(tǒng)仿真54MATLABm文件和函數(shù)2.3矩陣的運(yùn)算運(yùn)算符名稱指令示例說(shuō)明+加A+B若A、B為同維矩陣，則表示A與B對(duì)應(yīng)元素相加；若其中一個(gè)矩陣為標(biāo)量則表示另一矩陣的所有元素加上該標(biāo)量。－減A-B若A、B為同維矩陣，則表示A與B對(duì)應(yīng)元素相減；若其中一個(gè)矩陣為標(biāo)量則表示另一矩陣的所有元素減去該標(biāo)量*矩陣乘A*B矩陣A與B相乘，A和B均可為向量或標(biāo)量，但A和B的維數(shù)必須符合矩陣乘法的定義\矩陣左除A\B方程A*X=B

的解X/矩陣右除A/B方程X*A=B

的解X^矩陣乘方A^BA、B均為標(biāo)量時(shí)，表示A的B次方冪；A為方陣，B為正整數(shù)時(shí)，表示矩陣A的B次乘積；當(dāng)A、B均為矩陣時(shí)，無(wú)定義代數(shù)運(yùn)算

Matlab與系統(tǒng)仿真55MATLABm文件和函數(shù)函數(shù)功能d=size(A)[m,n]=size(A)將矩陣A的行數(shù)和列數(shù)賦值給變向量b將矩陣A的行數(shù)和列數(shù)分別賦值給變量m和nA'計(jì)算矩陣A的轉(zhuǎn)置矩陣inv(A)計(jì)算矩陣A的逆矩陣length(A)計(jì)算矩陣A的長(zhǎng)度（列數(shù)）sum(A)若A為向量，計(jì)算A所有元素之和；若A為矩陣，產(chǎn)生一行向量，其元素分別為矩陣A各列元素之和max(A)若A為向量，求出A所有元素的最大值；若A為矩陣，產(chǎn)生一行向量，其元素分別為矩陣A各列元素的最大值函數(shù)運(yùn)算

prod（A）min（A）Matlab與系統(tǒng)仿真56MATLABm文件和函數(shù)運(yùn)算符名稱指令示例說(shuō)明.*元素群乘A.*B矩陣A與B對(duì)應(yīng)元素相乘，A和B必須為同維矩陣或其中之一為標(biāo)量.\元素群左除A.\B矩陣B的元素除以矩陣A的對(duì)應(yīng)元素，A、B必須為同維矩陣或其中之一為標(biāo)量./元素群右除A./B矩陣A的元素除以矩陣B的對(duì)應(yīng)元素，A、B必須為同維矩陣或其中之一為標(biāo)量.^元素群乘方A.^B矩陣A的各元素與矩陣B的對(duì)應(yīng)元素的乘方運(yùn)算，即[A(i,j)^B(i,j)]，A、B必須為同維矩陣數(shù)組（元素群）運(yùn)算

Matlab與系統(tǒng)仿真57MATLABm文件和函數(shù)符號(hào)意義符號(hào)意義>大于>=大于或等于<小于<=小于或等于==等于~=不等于關(guān)系運(yùn)算

Matlab與系統(tǒng)仿真58MATLABm文件和函數(shù)符號(hào)意義&與邏輯|或邏輯~非邏輯邏輯運(yùn)算

Matlab與系統(tǒng)仿真59MATLABm文件和函數(shù)2.4MATLAB的程序流程控制

循環(huán)控制結(jié)構(gòu)

條件轉(zhuǎn)移結(jié)構(gòu)

Matlab與系統(tǒng)仿真60MATLABm文件和函數(shù)2.5M文件的編寫(xiě)

命令文件

函數(shù)文件

function[返回參數(shù)1,返回參數(shù)2,…]=函數(shù)名(輸入?yún)?shù)1,輸入?yún)?shù)2,…)%注釋說(shuō)明語(yǔ)句段，由％引導(dǎo)輸入、返回變量格式的檢測(cè)函數(shù)體語(yǔ)句函數(shù)文件具有特點(diǎn)

Matlab與系統(tǒng)仿真61MATLABm文件和函數(shù)例：求n!（MATLAB函數(shù)的可遞歸調(diào)用性）

functionk=my_fact(n)ifnargin~=1,error(‘輸入變量個(gè)數(shù)錯(cuò)誤，只能一個(gè)’)；endifnargout>1,error(‘輸出變量個(gè)數(shù)過(guò)多’)；endifabs(n-floor(n))>eps|n<0%判定n是否為非負(fù)整數(shù)

error(‘n應(yīng)該為非負(fù)整數(shù)’);endifn>1

%如果n>1,進(jìn)行遞歸調(diào)用

k=n*my_fact(n-1);elseifany([01]==n)%0!=1!=1為已知

k=1;end>>my_fact(11)

ans=39916800Matlab與系統(tǒng)仿真62MATLABm文件和函數(shù)2.6MATLAB的圖形繪制

二維圖形的繪制

函數(shù)調(diào)用格式說(shuō)明plot(y)y為向量，以y的序號(hào)作為X軸，按y的值繪制曲線plot(x,y,'option')x,y均為向量，以x作為x軸、y作為y軸繪制曲線，曲線的屬性由選項(xiàng)“option”來(lái)確定plot(x,y1,'option1',x,y2,'option2',…)以公共的向量x作為x軸，分別以向量y1，y2，…為y軸繪制多條曲線，每條曲線的屬性由相應(yīng)的選項(xiàng)“option”來(lái)確定plot(x1,y1,'option1',x1,y2,'option2',…)分別以向量xl，x2，…作為x軸，以y1，y2，…為y軸繪制多條曲線，每條曲線的屬性由相應(yīng)的選項(xiàng)“'option”來(lái)確定Matlab與系統(tǒng)仿真63MATLABm文件和函數(shù)符號(hào)顏色符號(hào)線型符號(hào)標(biāo)記'b'藍(lán)色'-'實(shí)線'v'▽'g'綠色'--'虛線'^'△'m'洋紅色':'點(diǎn)線'。'圓圈'w'白色'-.'點(diǎn)劃線'*'星號(hào)'c'青色'x'