第三章-Matlab圖形繪制

Matlab圖形繪制一、二維曲線的繪制繪制二維圖形最常用的函數(shù)是plot，調(diào)用格式有以下四種：plot(Y)如果Y為實(shí)向量，則以Y的索引坐標(biāo)作為橫坐標(biāo)，以Y本身的元素作為縱坐標(biāo)。如果Y為復(fù)數(shù)向量，則以該向量實(shí)部為橫坐標(biāo)，虛部為縱坐標(biāo)。plot(X,Y)此時(shí)以X為橫坐標(biāo)，Y為縱坐標(biāo)。X和Y要求維數(shù)必須相同。當(dāng)變量X和Y是同階矩陣時(shí)，將按照矩陣的行或列進(jìn)行操作。plot(X,Y,s)第三個(gè)變量用于設(shè)置圖形顯示屬性。設(shè)置圖形的線型、顏色、標(biāo)記等。plot(X1,Y1,s1,X2,Y2,s2,…..)每三項(xiàng)為一組繪制多條曲線。y=rand(100,1);plot(y)x=rand(100,1);z=x+y.*i;plot(z)示例：例子：多條曲線繪制t=0:0.1:2*pi;y=[sin(t);cos(t)];plot(t,y)（5）在MATLAB中繪圖時(shí)，當(dāng)線條多于一條時(shí)，若用戶沒(méi)有指定各線條的顏色時(shí)，MATLAB會(huì)自動(dòng)以不同的顏色區(qū)分不同的線條，可以查看默認(rèn)的顏色順序：其顏色分別是：藍(lán)色、深綠色、紅色、淺橄欖綠色、紫紅色、土黃色或灰色。注意：當(dāng)指定了數(shù)據(jù)點(diǎn)標(biāo)志類型，但不指定線型時(shí)，則表示只標(biāo)記數(shù)據(jù)點(diǎn)，而不進(jìn)行連線繪圖。

x=linspace(0,2*pi,1000);y1=0.2*exp(-0.5*x).*cos(4*pi*x);y2=2*exp(-0.5*x).*cos(pi*x);plotyy(x,y1,x,y2)二、圖形處理坐標(biāo)軸的控制Grid命令Box命令圖形數(shù)據(jù)取點(diǎn)圖形保持Hold特殊坐標(biāo)系的圖形函數(shù)在基本的繪圖函數(shù)使用時(shí)，它們的坐標(biāo)軸刻度均為線性刻度。但是當(dāng)實(shí)際的數(shù)據(jù)出現(xiàn)指數(shù)變化時(shí)，指數(shù)變化就不能直觀的從圖形上體現(xiàn)出來(lái)。又假如當(dāng)函數(shù)以極坐標(biāo)等表示時(shí)，我們就不能用線性刻度來(lái)表示。為了解決這些問(wèn)題，matlab提供了相應(yīng)的繪圖函數(shù)。這里所謂的特殊坐標(biāo)系是區(qū)別于均勻直角坐標(biāo)系而言，具體來(lái)說(shuō)就是對(duì)數(shù)坐標(biāo)系、極坐標(biāo)系、柱坐標(biāo)和球坐標(biāo)等。對(duì)數(shù)坐標(biāo)系MATLAB語(yǔ)言提供了繪制不同形式的對(duì)數(shù)坐標(biāo)曲線的功能，具體實(shí)現(xiàn)該功能的函數(shù)是semilogx、semilogy

和

loglog。semilogx（x軸對(duì)數(shù)刻度坐標(biāo)圖）用該函數(shù)繪制圖形時(shí)x軸采用對(duì)數(shù)坐標(biāo)。

semilogx（y）：對(duì)x軸的刻度求常用對(duì)數(shù)（以10為底），而y為線性刻度。對(duì)數(shù)坐標(biāo)系>>x=0.001:0.01*pi:2*pi;>>y=log10(x);>>semilogx(x,y,'-*');（x軸對(duì)數(shù)刻度坐標(biāo)圖）>>plot(x,y)（均勻直角坐標(biāo)系）對(duì)數(shù)坐標(biāo)系

Semilogx繪制的y=log10(x)plot繪制的y=log10(x)

如圖所示，x=（0.001~2Π）以0.01Π為間隔。semilogx函數(shù)中x軸的相鄰坐標(biāo)值之間的刻度是不等距的，它是是按照相等的指數(shù)變化來(lái)增加的

。而plot函數(shù)中的相鄰坐標(biāo)值之間的刻度是等距的。對(duì)數(shù)坐標(biāo)系semilogy（y軸對(duì)數(shù)刻度坐標(biāo)圖）用該函數(shù)繪制圖形時(shí)y軸采用對(duì)數(shù)坐標(biāo)。調(diào)用格式與semilogx基本相同。>>x=0.001:0.1*pi:2*pi;>>y=10.^x;>>semilogy(x,y,'-*')>>plot(x,y)semilogy函數(shù)plot函數(shù)對(duì)數(shù)坐標(biāo)系loglog

用該函數(shù)繪制圖形時(shí)x和y軸均采用對(duì)數(shù)坐標(biāo)。調(diào)用格式與semilogx基本相同。>>m=1:1*pi:20*pi;>>x=10.^m;>>y=log10(m);>>loglog(x,y,'-*');>>plot(x,y,'-*')loglog函數(shù)plot函數(shù)極坐標(biāo)系polar命令該命令用于畫(huà)極坐標(biāo)圖，它接受極坐標(biāo)形式的函數(shù)rho=f(θ)。

polar的調(diào)用格式：

polar(theta,rho,LineSpec)

用極角theta和極徑rho畫(huà)出極坐標(biāo)圖形。極角theta為從x軸到半徑的單位為弧度的向量，極徑rho為各數(shù)據(jù)點(diǎn)到極點(diǎn)的半徑向量,LineSpec指定極坐標(biāo)圖中線條的線型、標(biāo)記符號(hào)和顏色等。

極坐標(biāo)系>>rho0=1;theta=0:pi/20:4*pi;rho=rho0+theta*rho0;polar(theta,rho,'r')polar的調(diào)用格式：

polar(theta,rho)用極角theta和極徑rho畫(huà)出極坐標(biāo)圖形。極角theta為從x軸到半徑的單位為弧度的向量，極徑rho為各數(shù)據(jù)點(diǎn)到極點(diǎn)的半徑向量。Matlab實(shí)現(xiàn)頻數(shù)直方圖----hist的使用

“hist”是“Histogram”(直方圖、柱狀圖)的簡(jiǎn)稱。

1、N=hist(Y)功能：將向量Y的元素平均到十個(gè)等間隔的容器中，并且返回每個(gè)容器的元素個(gè)數(shù)；如果Y是一個(gè)矩陣，hist指令逐列元素操作。例1：執(zhí)行指令Y=[1:10];Hist(Y)得到：10個(gè)藍(lán)色方條，每個(gè)方條對(duì)應(yīng)一個(gè)容器，其長(zhǎng)度代表容器中數(shù)據(jù)的多少。由圖知，容器中的數(shù)據(jù)量均為1。例2：執(zhí)行指令Y=[1,2,2,5,6,6,8,11]Hist(Y)得到：長(zhǎng)條圖之繪製Bar(y)為。範(fàn)例5-1：bar01.m x=[13452]; bar(x);Fig.5-15-1長(zhǎng)條圖之繪製(cont.)bar指令也可接受矩陣輸入，它會(huì)將同一橫列的資料聚集在一起。範(fàn)例5-2：bar02.mx=[23457;12321];bar(x);Fig.5-25-1長(zhǎng)條圖之繪製(cont.)bar及barh指令還有一項(xiàng)特異功能，就是可以將同一橫列的資料以堆疊（Stack）方式來(lái)顯示。範(fàn)例5-3：bar03.mx=[23457;12321];bar(x,'stack')Fig.5-35-1長(zhǎng)條圖之繪製(cont.)除了平面長(zhǎng)條圖之外，MATLAB亦可使用bar3指令來(lái)畫(huà)出立體長(zhǎng)條圖。範(fàn)例5-4：bar04.m x=[23457;12321]; bar3(x)Fig.5-45-1長(zhǎng)條圖之繪製(cont.)bar3指令還可以使用群組（Group）方式來(lái)呈現(xiàn)長(zhǎng)條圖範(fàn)例5-5：bar05.m x=[23457;12321]; bar3(x,'group')Fig.5-55-1長(zhǎng)條圖之繪製(cont.)長(zhǎng)條圖的指令和類別:垂直長(zhǎng)條圖水平長(zhǎng)條圖平面barbarh立體bar3bar3h5-1長(zhǎng)條圖之繪製(cont.)若要指定長(zhǎng)條圖的x座標(biāo)，可使用兩個(gè)輸入向量給bar指令。假設(shè)新竹的月平均溫度如下：範(fàn)例5-6：bar06.mx=1:6; %月份 y=35*rand(1,6); %溫度值（假設(shè)是介於0～35的亂數(shù)）bar(x,y); xlabel('月份'); %x軸的說(shuō)明文字ylabel('平均溫度(^{o}c)'); %y軸的說(shuō)明文字%下列指令將x軸的數(shù)字改成月數(shù)set(gca,'xticklabel',{'一月','二月','三月','四月','五月','六月'});5-1長(zhǎng)條圖之繪製(cont.)Fig.5-65-2面積圖之繪製面積圖（AreaGraphs）和以堆疊方式呈現(xiàn)的長(zhǎng)條圖很類似，特別適用於具有疊加關(guān)係的資料。舉例來(lái)說(shuō)，若要顯示清華大學(xué)在過(guò)去10年來(lái)的人數(shù)（含大學(xué)部，研究生，及教職員）變化情況，可用面積圖顯示。範(fàn)例5-7：area01.m y=rand(10,3)*100; x=1:10; area(x,y); xlabel('Year'); ylabel('Count')Fig.5-75-3扇形圖之繪製使用pie指令，可畫(huà)出平面扇形圖（PieCharts），並可加上說(shuō)明。範(fàn)例5-8：pie01.m x=[2354]; label={'東','南','西','北'}; pie(x,label);Fig.5-85-3扇形圖之繪製(cont.)pie指令直接將x元素視為面積百分比，因此可畫(huà)出不完全的扇形圖。範(fàn)例5-9：pie02.m x=[0.21,0.14,0.38]; pie(x);Fig.5-95-3扇形圖之繪製(cont.)pie指令還有一特異功能，可將某個(gè)或數(shù)個(gè)扇形圖向外拖出，以強(qiáng)調(diào)部份資料。範(fàn)例5-10：pie03.m x=[2354]; explode=[1100]; pie(x,explode);Fig.5-10其中指令explode中非零的元素即代表要向外拖出的扇形。5-3扇形圖之繪製(cont.)欲畫(huà)出立體扇形圖，可用pie3指令。範(fàn)例5-11:pie301.m x=[2354]; explode=[1100]; label={'春','夏','秋','冬'}; pie3(x,explode,label);Fig.5-115-4針頭圖之繪製顧名思義，針頭圖（StemPlots）就是以一個(gè)大頭針來(lái)表示某一點(diǎn)資料，其指令為stem。範(fàn)例5-12：stem01.m t=0:0.2:4*pi; y=cos(t).*exp(-t/5); stem(t,y)Fig.5-125-4針頭圖之繪製(cont.)針頭圖特別適用於表示「數(shù)位訊號(hào)處理」（DSP，DigitalSignalProcessing）中的數(shù)位訊號(hào)。若要畫(huà)出實(shí)心的針頭圖，可加“fill”選項(xiàng)。範(fàn)例5-13：stem02.m t=0:0.2:4*pi; y=cos(t).*exp(-t/5); stem(t,y,'fill');Fig.5-135-4針頭圖之繪製(cont.)欲畫(huà)出立體的針頭圖， 可用stem3指令。範(fàn)例5-14：stem301.m theta=-pi:0.05:pi; x=cos(theta); y=sin(theta); z=abs(cos(3*theta)).*exp(-abs(theta/3)); stem3(x,y,z);Fig.5-145-5階梯圖之繪製使用stairs指令，可畫(huà)出階梯圖（StairstepPlots），其精神和針頭圖很相近，只是將目前資料點(diǎn)的高度向右水平畫(huà)至下一點(diǎn)為止。（在數(shù)位訊號(hào)處理，此種作法稱為Zero-orderHold。）範(fàn)例5-15：stairs01.m t=0:0.4:4*pi; y=cos(t).*exp(-t/5); stairs(t,y);Fig.5-155-5階梯圖之繪製(cont.)若再加上針頭圖，則可見(jiàn)兩 者相似之處。範(fàn)例5-16：stairs02.m t=0:0.4:4*pi; y=cos(t).*exp(-t/5); stairs(t,y); holdon %保留舊圖形 stem(t,y); %疊上針頭圖 holdoffFig.5-165-6實(shí)心圖之繪製MATLAB指令fill將資料點(diǎn)視為多邊形頂點(diǎn)，並將此多邊形塗上顏色，呈現(xiàn)出實(shí)心圖（FilledPlots）的結(jié)果。範(fàn)例5-17:fill01.m t=0:0.4:4*pi; y=sin(t).*exp(-t/5); fill(t,y,'b'); %'b'為藍(lán)色Fig.5-175-6實(shí)心圖之繪製(cont.)若與stem合用，則可創(chuàng)造出 一些不同的視覺(jué)效果。範(fàn)例5-18:fill02.m t=0:0.4:4*pi; y=sin(t).*exp(-t/5); fill(t,y,'y'); %'y'為黃色 holdon %保留舊圖形 stem(t,y,'b'); %疊上藍(lán)色針頭圖 holdoffFig.5-185-6實(shí)心圖之繪製(cont.)fill3可用於三維的實(shí)心圖。範(fàn)例5-19:fill301.m X=[0011]; Y=[0110]; Z=[0110]; C=[00.30.60.9]'; fill3(X,Y,Z,C);Fig.5-195-6實(shí)心圖之繪製(cont.)使用fill3指令，我們亦可以 畫(huà)出各種酷酷的圖形。範(fàn)例5-20:fill302.m t=(1/16:1/8:1)'*2*pi; x=sin(t); y=cos(t); c=linspace(0,1,length(t)); fill3(x,y/sqrt(2),y/sqrt(2),c,x/sqrt(2),y,x/sqrt(2),c); axistightFig.5-205-7向量場(chǎng)圖之繪製使用quiver指令可畫(huà)出平面 上的向量場(chǎng)圖（QuiverPlots） ，特別適用於表示分布於平面 的向量場(chǎng)（VectorFields）， 例如平面上的電場(chǎng)分布，或是流 速分布。範(fàn)例5-21：quiver01.m [x,y,z]=peaks(20); [u,v]=gradient(z); contour(x,y,z,10); holdon,quiver(x,y,u,v);holdoff axisimageFig.5-215-7向量場(chǎng)圖之繪製(cont.)欲畫(huà)出空間中的向量場(chǎng)圖， 可用quiver3指令。範(fàn)例5-22：quiver301.m [x,y]=meshgrid(-2:0.2:2,-1:0.1:1); z=x.*exp(-x.^2-y.^2); [u,v,w]=surfnorm(x,y,z); quiver3(x,y,z,u,v,w); holdon,surf(x,y,z);holdoff axisequalFig.5-225-8等高線圖之繪製我們可用contour指令來(lái)畫(huà)出「等高線圖」（ContourPlots）。範(fàn)例5-23：contour01.m z=peaks; contour(z,30); %畫(huà)出30條等高線Fig.5-235-8等高線圖之繪製(cont.)若要畫(huà)出特定高度的等高線，可執(zhí)行如下：範(fàn)例5-24：contour02.m z=peaks; contour(z,[025]);Fig.5-245-8等高線圖之繪製(cont.)欲標(biāo)明等高線的高度， 可用clabel指令。範(fàn)例5-25：contour03.m z=peaks; [c,handle]=contour(z,10); clabel(c,handle);Fig.5-255-8等高線圖之繪製(cont.)若欲在等高線之間填入顏色，可用contourf 指令。範(fàn)例5-26：contour04.m z=peaks; contourf(z);Fig.5-265-8等高線圖之繪製(cont.)若要使畫(huà)出的等高線對(duì) 應(yīng)至正確的x及y 座標(biāo)，則可執(zhí)行如下：範(fàn)例5-27：contour05.m [x,y,z]=peaks; contour(x,y,z);%使用三個(gè)輸入Fig.5-275-8等高線圖之繪製(cont.)contourf亦可接受x、y、z輸入引數(shù)。若要將等高線畫(huà)在曲面的正下方，可用surfc或meshc指令。範(fàn)例5-28：contour06.m [x,y,z]=peaks; meshc(x,y,z); axistightFig.5-285-8等高線圖之繪製(cont.)若要畫(huà)出三度空間中的等高線，可用contour3指令。範(fàn)例5-29：contour301.m [x,y,z]=peaks; contour3(x,y,z,30); axistighFig.5-295-8等高線圖之繪製(cont.)使用contour指令亦可畫(huà)出極座標(biāo)中的等高線，但過(guò)程較為複雜，以下列複數(shù)函數(shù)為例： 其中z代表複數(shù)平面中的任一點(diǎn)複數(shù)，如果我們要畫(huà)出此函數(shù)的等高線，可見(jiàn)下列範(fàn)例：範(fàn)例5-30：contour07.m t=linspace(0,2*pi,61); %角度的格子點(diǎn) r=0:0.05:1; %長(zhǎng)度的格子點(diǎn) [tt,rr]=meshgrid(t,r); %產(chǎn)生二維的格子點(diǎn) [xx,yy]=pol2cart(tt,rr); %將極座標(biāo)轉(zhuǎn)換至直角座標(biāo) zz=xx+sqrt(-1)*yy; %複數(shù)表示 ff=abs(zz.^3-1); %曲面的函數(shù) contour(xx,yy,ff,50); %畫(huà)出等高線 axisimage5-8等高線圖之繪製(cont.)Fig.5-305-8等高線圖之繪製(cont.)在上例中，座標(biāo)的標(biāo)示仍為直 角座標(biāo)。欲將等高線顯示於極 座標(biāo)上，需先用polar指令產(chǎn) 生一個(gè)極座標(biāo)圖，再移除圖形 ，留下圖軸，然後再進(jìn)行作圖 。範(fàn)例5-31：contour08.m h=polar([02*pi],[01]); delete(h); holdon contour(x,y,abs(f),30); holdoff%產(chǎn)生在極座標(biāo)上的一條直線%移除上述圖形，但留下極座標(biāo)圖軸Fig.5-31＊要注意的是，你必須先執(zhí)行contour07.m，然後再執(zhí)行contour08.m，才能得到上述的及座標(biāo)等高線的效果。5-8等高線圖之繪製(cont.)我們也可以同時(shí)畫(huà)出複數(shù)函數(shù)的曲面和等高線圖，例如，下列範(fàn)例可以畫(huà)出複數(shù)函數(shù):範(fàn)例5-32：contour09.m t=linspace(0,2*pi,61); %角度的格子點(diǎn) r=0:0.05:1; %長(zhǎng)度的格子點(diǎn) [tt,rr]=meshgrid(t,r); %產(chǎn)生二維的格子點(diǎn) [xx,yy]=pol2cart(tt,rr); %將極座標(biāo)轉(zhuǎn)換至直角座標(biāo) zz=xx+sqrt(-1)*yy; %複數(shù)表示 ff=abs(zz.^3-1); %曲面的函數(shù)值 h=polar([02*pi],[01]); %產(chǎn)生在極座標(biāo)上的一條直線 delete(h); %移除上述圖形，但留下極座標(biāo)圖軸 holdon contour(xx,yy,ff,20); %等高線 surf(xx,yy,ff); %曲面圖 holdoff view(-19,22); %設(shè)定觀測(cè)角度5-8等高線圖之繪製(cont.)Fig.5-325-8等高線圖之繪製(cont.)指令說(shuō)明bar,barh,bar3,bar3h長(zhǎng)條圖Area面積圖pie,pie3扇形圖stem,stem3針頭圖stairs階梯圖fill,fill3實(shí)心圖quiver,quiver3向量場(chǎng)圖contour,contourf,contour3等高線圖特殊繪圖函數(shù):5-9 其他進(jìn)階繪圖功能MATLAB在5.3版後，開(kāi)始支援「容積目視法」（VolumeVisualization）、因此能夠畫(huà)出在三度空間中的流線圖、向量場(chǎng)圖、等高面圖（Isosurfaces）、切面圖（Slices）等，相關(guān)指令可列表如下頁(yè)：5-9 其他進(jìn)階繪圖功能(cont.)指令說(shuō)明coneplot以圓錐瓶畫(huà)出三度空間的向量場(chǎng)圖contourslice在三度空間的切面上畫(huà)出等高線isosurface從容積資料中算出等高面資料isocaps計(jì)算等高面在端點(diǎn)切片的等高資訊isonormals計(jì)算等高面的法向量slice在三度空間的切片streamline從2–D或3–D的流線資料來(lái)畫(huà)出流線圖isocolors計(jì)算等高區(qū)面頂點(diǎn)的顏色（第六版才支援）divergence計(jì)算3-D向量場(chǎng)的亂度（Divergence）（第六版才支援）curl計(jì)算3-D向量場(chǎng)的curl及垂直方向的角速度（第六版才支援）※

MATLAB有關(guān)於「容積目視法」的指令:

5-9 其他進(jìn)階繪圖功能(cont.)streamtube由向量資料畫(huà)出流線管（StreamTubes）（第六版才支援）streamribbon由向量資料畫(huà)出流線緞帶（StreamRibbons）（第六版才支援）streamslice由向量資料畫(huà)出間隔分明的流線(第六版才支援)streamparticles由向量資料畫(huà)出流線粒子（StreamParticles）（第六版才支援）interpstreamspeed由速度對(duì)流線頂點(diǎn)做內(nèi)差（Interpolation）（第六版才支援）volumebounds傳回容積資料的座標(biāo)及顏色極限值（第六版才支援）※

MATLAB有關(guān)於「容積目視法」的指令:

5-9 其他進(jìn)階繪圖功能(cont.)這些指令的用法較為繁複，由於篇幅有限，在此不詳細(xì)說(shuō)明，讀者可查閱相關(guān)的線上支援。若要一睹這些「容積目視」指令所能創(chuàng)造的繽紛效果，您可在MATLAB指令視窗下輸入「volvec」，以開(kāi)啟展示視窗，在點(diǎn)選「Multiple」之後，產(chǎn)生圖形如右：②極坐標(biāo)系函數(shù)polar，調(diào)用形式為：polar(theta,rho)或polar(theta,rho,s)③雙縱坐標(biāo)（雙y軸坐標(biāo)系）函數(shù)plotyy，調(diào)用形式為：plotyy(X1,Y1,X2,Y2)plotyy(X1,Y1,X2,Y2,fun)fun可以是plot、semilogx、semilogy或loglogplotyy(X1,Y1,X2,Y2,fun1,fun2)fun1繪制(X1,Y1)，fun2繪制(X2,Y2)注：雙坐標(biāo)繪制圖形的調(diào)用過(guò)程中，不能夠像前面的plot函數(shù)那樣對(duì)曲線屬性進(jìn)行設(shè)置，需要使用句柄圖形控制完成。x=0:0.01:20;y1=200*exp(-0.05*x).*sin(x);y2=0.8*exp(-0.5*x).*sin(10*x);[AX,H1,H2]=plotyy(x,y1,x,y2,'plot');AX是對(duì)應(yīng)兩個(gè)坐標(biāo)軸的句柄，H1和H2分別是對(duì)應(yīng)plot的句柄。（2）二維特殊函數(shù)圖形函數(shù)名說(shuō)明函數(shù)名說(shuō)明函數(shù)名說(shuō)明area填充繪圖feather矢量圖fplot函數(shù)圖繪制bar條形圖stem離散序列餅狀圖hist直方圖barh水平條形圖fill多邊形填充paretoPareto圖comet彗星圖stairs階梯圖pie餅圖errorbar誤差帶圖contour等高線圖plotmatrix分散矩陣?yán)L制ezplot簡(jiǎn)單繪制函數(shù)圖contourf填充的等高線圖ribbon三維圖的二維條狀顯示ezpolar簡(jiǎn)單繪制極坐標(biāo)圖scatter散射圖Polar繪制極坐標(biāo)示例：theta=linspace(0,2*pi);r=cos(2*theta);polar(theta,r)Hist用來(lái)顯示資料的分段情況和統(tǒng)計(jì)特性，適合于大量數(shù)據(jù)的情況示例：x=randn(9999,1);hist(x,50)Rose與hist接近，將資料的大小視為角度，資料的個(gè)數(shù)作為距離，采用極坐標(biāo)繪圖。示例：x=randn(9999,1);rose(x,50)Bar函數(shù)適合于數(shù)據(jù)點(diǎn)不多時(shí)使用示例：x=1:15;y=rand(size(x));bar(x,y);如果已知數(shù)據(jù)的誤差量，可用errorbar表示。示例：x=linspace(0,2*pi,20);y=sin(x);z=std(y)*ones(size(x));errorbar(x,y,z);fplot函數(shù)可以繪制較精確的圖形，適合于變化劇烈的函數(shù)。示例：fplot(‘sin(1/x)’,[0.01,0.3])Stem產(chǎn)生針狀圖，常被用來(lái)繪制數(shù)位訊號(hào)。x=linspace(0,10,100);y=sin(x).*exp(-x/4);stem(x,y);Fill將資料點(diǎn)視為多邊形頂點(diǎn)，并將此多邊形涂上顏色。x=linspace(0,10,100);y=sin(x).*exp(-x/4);fill(x,y,’c’);Feather函數(shù)將每一個(gè)數(shù)據(jù)點(diǎn)視為復(fù)數(shù)，并以箭頭畫(huà)出。theta=linspace(0,2*pi,40);z=cos(theta)+i*sin(theta);feather(z);Compass函數(shù)與feather函數(shù)接近，只是每個(gè)箭頭的起點(diǎn)都在圓內(nèi)theta=linspace(0,2*pi,40);z=cos(theta)+i*sin(theta);compass(z);Quiver用來(lái)繪制向量場(chǎng)圖。[x,y]=meshgrid(-2:0.2:2);z=x.*exp(-x.^2-y.^2);[Dx,Dy]=gradient(z,0.2,0.2);Contour(x,y,z);holdon;quiver(x,y,Dx,Dy);colormaphsv;holdoff3、顯函數(shù)、隱函數(shù)和符號(hào)函數(shù)作圖(1)fplot(fun,lims)繪制由字符串fun指定函數(shù)名的函數(shù)在x軸區(qū)間為lims=[xmin,xmax]的函數(shù)圖。若lims=[xmin,xmax,ymin,ymax]，則y軸也被限制。例題：解

1）首先用M文件fun1.m定義函數(shù)f(x)如下functiony=Afun1(x);ifx<1y=x+1;elsey=1+1./x;end在matlab命令窗口輸入fplot('Afun1',[-3,3])就可畫(huà)出函數(shù)f(x)的圖形。2）這里也可以使用匿名函數(shù)，編寫(xiě)程序如下fun2=@(x)(x+1)*(x<1)+(1+1/x)*(x>=1);fplot(fun2,[-3,3])(2)ezplot(f)可用來(lái)繪制顯函數(shù)、隱函數(shù)和參數(shù)方程的圖形，x軸的默認(rèn)范圍為[-2pi,2pi]。1)對(duì)于顯函數(shù)y=f(x)，調(diào)用格式為：ezplot(f,[xmin,xmax])2)對(duì)于隱函數(shù)f(x,y)=0，調(diào)用格式為：ezplot(f,[xmin,xmax],[ymin,ymax])ezplot(f,[a,b]),x與y在（a,b）區(qū)間之內(nèi)。3)對(duì)于參數(shù)方程x=x(t),y=y(t)，調(diào)用格式為：ezplot(x,y)或者ezplot(f,[tmin,tmax])示例1：ezplot(‘x^2*sin(x+y^2)+y^2*exp(x+y)+5*cos(x^2+y)’,[-1010]);示例2：ezplot('cot(x)')示例3：ezplot('x^2+y^2/4=1')示例4：ezplot(‘x^2+y^2-1’,[-1,1],[-1,1])示例5：ezplot(‘x^3+y^3-5*x*y’)示例6：ezplot(‘2*cos(t)’,’sin(t)’,[0,2*pi]);3.2三維曲線的繪制繪制三維曲線的函數(shù)是plot3，與plot函數(shù)的調(diào)用格式相同，也有以下四種：plot3(x,y,z)plot3(X,Y,Z)X,Y,Z是相同階數(shù)的矩陣，繪出三矩陣的列向量的曲線。plot3(X,Y,Z,s)plot3(x1,y1,z1,s1,x2,y2,z2,s2,....)組合圖形以第二種形式為例：[x,y]=meshgrid(-2:0.1:2,-2:0.1:2);z=x.*exp(-x.^2-y.^2);plot3(x,y,z)函數(shù)名說(shuō)明函數(shù)名說(shuō)明mesh三維網(wǎng)格圖meshc將網(wǎng)格與等高線結(jié)合meshz屏蔽的網(wǎng)格圖meshgrid生成網(wǎng)格點(diǎn)MATLAB對(duì)于網(wǎng)格的處理方法是：將xy平面按指定方式分隔成平面網(wǎng)格，然后根據(jù)程序中給定的方式計(jì)算第三維變量的值，即z軸的值，與對(duì)應(yīng)的xy平面的坐標(biāo)構(gòu)成三維點(diǎn)元素，根據(jù)由此得到的(x,z)和(y,z)計(jì)算各平面的曲線，彼此相連就構(gòu)成了網(wǎng)格圖。3.3三維網(wǎng)線圖和曲面圖的繪制（1）網(wǎng)線圖mesh函數(shù)[X,Y]=meshgrid(x,y)x和y是給定的向量，可以定義網(wǎng)格的劃分區(qū)域，也可以定義網(wǎng)格的劃分方法，X和Y是網(wǎng)格劃分后的數(shù)值矩陣。mesh(X,Y,Z,C)繪制四個(gè)矩陣變量的彩色網(wǎng)格面圖形，觀測(cè)點(diǎn)由view定義，坐標(biāo)軸由axis定義，C設(shè)置顏色，由colormap實(shí)現(xiàn)。mesh(X,Y,Z)使用C=Z，即顏色正比于圖高。mesh(x,y,Z,C)由兩個(gè)向量代替矩陣，要求length(x)=n,length(y)=m,且[m,n]=size(Z)示例1：x=-2:0.1:2;y=-2:0.1:2;[xx,yy]=meshgrid(x,y);zz=xx.^2+yy.^2;%ii=find(xx.^2+yy.^2>4);zz(ii)=NaN;surf(xx,yy,zz);x=-10:0.5:10;y=-8:0.5:8;[x,y]=meshgrid(x,y);z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2);mesh(x,y,z);示例2：（2）著色函數(shù)surfsurf(X,Y,Z,C)輸入?yún)?shù)的設(shè)置和mesh相同，不同的是mesh繪制的是網(wǎng)格圖，而surf繪制的著色的三維表面。兩者對(duì)比如下圖所示函數(shù)名說(shuō)明函數(shù)名說(shuō)明函數(shù)名說(shuō)明bar3三維條形圖surfc著色圖+等高線comet3三維彗星圖trisurf三角形表面圖ezgraph3函數(shù)控制繪制圖waterfall瀑布圖pie3三維餅狀圖scatter3三維散射圖cylinder柱面圖stem3三維離散數(shù)據(jù)圖sphere球面圖trimesh三角形網(wǎng)格圖（3）特殊的三維繪圖命令（4）旋轉(zhuǎn)曲面圖cylinder[X,Y,Z]=cylinder(r)這里的r表示構(gòu)成旋轉(zhuǎn)曲面的曲線。方法一：

Matlab程序如下。x=0:10:600;[X,Y,Z]=cylinder(30*exp(-x/400).*sin((x+25*pi)/100)+130);surf(X,Y,Z)方法二：將旋轉(zhuǎn)曲面用參數(shù)方程表示。

解因?yàn)檫@里的函數(shù)是隱函數(shù)，化成顯函數(shù)后有兩支，必須使用參數(shù)方程，旋轉(zhuǎn)面的參數(shù)方程為畫(huà)圖的Matlab程序如下：alpha=[0:0.1:2*pi]';beta=0:0.1:2*pi;x=4*cos(alpha)*ones(size(beta));y=(5+4*sin(alpha))*cos(beta);z=(5+4*sin(alpha))*sin(beta);surf(x,y,z)畫(huà)圖的Matlab程序也可以寫(xiě)成x=@(alpha,beta)4*cos(alpha);y=@(alpha,beta)(5+4*sin(alpha))*cos(beta);z=@(alpha,beta)(5+4*sin(alpha))*sin(beta);ezsurf(x,y,z)（5）繪制柱面圖柱面平行于某個(gè)坐標(biāo)軸，方程中不出現(xiàn)某個(gè)坐標(biāo)軸的變量，方程表示為F(x,y)=0或者F(x,z)=0或者F(y,z)=0示例一：畫(huà)出方程

表示的柱面解：方程中的x是自變量矩陣，z是因變量，則另一個(gè)自變量矩陣為y，自變量平面是xoy面，x軸是真正的自變量，y軸是柱面方向。Matlab程序如下：u=linspace(-5,5,10)';%設(shè)定參數(shù)列向量uv=linspace(-5,10,10);%設(shè)定參數(shù)行向量vX=u*ones(size(v));%構(gòu)成自變量矩陣XY=ones(size(u))*v;%構(gòu)成自變量矩陣YZ=2-X.^2;%求因變量Zmesh(X,Y,Z)（6）其他二次曲面Matlab中使用繪圖命令ezmesh或ezsurf也很方便，只需要把曲面方程寫(xiě)成兩個(gè)變量的顯函數(shù)方程或參數(shù)方程即可。對(duì)于旋轉(zhuǎn)面，如果母線的方程可以表示成關(guān)于旋轉(zhuǎn)軸變量的顯式函數(shù)，則可以直接使用Matlab工具箱中的命令cylinder，否則必須把旋轉(zhuǎn)面化成參數(shù)方程，然后使用ezmesh或ezsurf命令繪圖。對(duì)于其它的二次曲面，如果可以寫(xiě)成顯函數(shù)直接使用命令ezmesh或ezsurf，否則必須先化成參數(shù)方程。幾個(gè)示例：（1）x=@(s,t)3*tan(s)*cos(t);y=@(s,t)2*tan(s)*sin(t);z=@(s,t)tan(s);ezsurf(x,y,z)（2）ellipsoid(0,0,0,3,2,sqrt(6))（3）ezsurf(@(x,y)x*y)（4）x=@(s,t)3*cos(s);y=@(s,t)2*sin(s);z=@(s,t)t;ezmesh(x,y,z)（7)Matlab三維散點(diǎn)插值曲面圖的繪制1、scatter3函數(shù)，調(diào)用方法：scatter3(X,Y,Z)示例1：[x,y,z]=sphere(16);X=[x(:)*.5x(:)*.75x(:)];Y=[y(:)*.5y(:)*.75y(:)];Z=[z(:)*.5z(:)*.75z(:)];S=repmat([1.75.5]*10,numel(x),1);C=repmat([123],numel(x),1);scatter3(X(:),Y(:),Z(:),S(:),C(:),'filled'),view(-60,60)示例2：x=[4229042.63

4230585.02

4231384.96

4231773.63

4233028.58

4233296.71

4235869.68

4236288.29];

y=[431695.4

441585.8

432745.6436933.7

428734.4

431946.3428705.0

432999.5];

z=[1.0191.023

1.011

1.022

1.020

1.022

1.022

1.023];

scatter3(x,y,z)2、griddata函數(shù)，用來(lái)產(chǎn)生經(jīng)插值后均勻間隔數(shù)據(jù)作圖。常用的調(diào)用方法是：[XI,YI,ZI]=griddata(x,y,z,XI,YI,’method’)，其中x,y,z來(lái)自關(guān)系式z=f(x,y)，也可以是離散的數(shù)據(jù)點(diǎn)，ZI是有XI和YI插值后得到的對(duì)應(yīng)點(diǎn)。其中method可以是linear，cubic和nearest。V1=[110;

-110;

-1-10;

1-10];x=V1(:,1);y=V1(:,2);z=V1(:,3);xmin=min(x);xmax=max(x);ymin=min(y);ymax=max(y);x0=linspace(xmin,xmax,20);y0=linspace(ymin,ymax,20);[x0,y0]=meshgrid(x0,y0);z0=griddata(x,y,z,x0,y0,'v4');figure;holdon;view(3);plot3(x,y,z,'r.','MarkerSize',20);plot3(x0,y0,z0,'k.');col=zeros(size(z0));col(:,end)=1;mesh(x0,y0,z0,col)示例1：x=10*(0.5-rand(1000,1));y=10*(0.5-rand(1000,1));z=sqrt(x.*x+y.*y).*(9.995+rand(1000,1)/100);clf;plot3(x,y,z,'rp','markerSize',4);holdon;[xx,yy]=meshgrid(linspace(min(x),max(x)),linspace(min(y),max(y)));zz=griddata(x,y,z,xx,yy,'cubic');mesh(xx,yy,zz);shadinginterp;colormapgray;holdon;示例2：示例3：x=rand(100,1)*4-2;y=rand(100,1)*4-2;z=x.*exp(-x.^2-y.^2);ti=-2:.25:2;[xi,yi]=meshgrid(ti,ti);zi=griddata(x,y,z,xi,yi);mesh(xi,yi,zi),holdon,plot3(x,y,z,'o'),holdoff3、TriScatteredInterp函數(shù)，該函數(shù)可以看作是griddata的升級(jí)版，可以取代griddata對(duì)離散數(shù)據(jù)點(diǎn)進(jìn)行插值。主要用于不規(guī)則數(shù)據(jù)點(diǎn)的插值。主要調(diào)用方法：F=TriScatteredInterp(x,y,z)，相當(dāng)于定義了一個(gè)關(guān)于x,y,z的插值函數(shù)。示例1：x=10*(0.5-rand(1000,1));y=10*(0.5-rand(1000,1));%生成隨機(jī)數(shù)據(jù)z=sqrt(x.*x+y.*y).*(9.995+rand(1000,1)/100);clf;plot3(x,y,z,‘rp’,‘markerSize’,4);%繪制三維散點(diǎn)圖holdon;F=TriScatteredInterp(x,y,z);%可以理解為定義了一個(gè)函數(shù)F(.)[xx,yy]=meshgrid(linspace(min(x),max(x)),linspace(min(y),max(y)));zz=F(xx,yy);%把x，y的坐標(biāo)輸入其中，可以計(jì)算出插值的z維度數(shù)值mesh(xx,yy,zz);%然后可以采用常規(guī)的方法繪制曲面shadinginterp;colormapgray;holdon;注：1、對(duì)F還可以設(shè)置方法F.method=‘nature’/’linear’/’nearest’,前兩個(gè)針對(duì)連續(xù)曲線。2、類似的一些地形圖的繪制也是通過(guò)這種插值的方法進(jìn)行繪制。示例2：z=x.^2+y.^2;F=TriScatteredInterp(x,y,z);[xi,yi]=meshgrid(-0.3:.02:0.3,-0.0688:0.01:0.0688);zi=F(xi,yi);mesh(xi,yi,zi)xlabel('Interpolatedsurface','fontweight','b');figure;zi=xi.^2+yi.^2;mesh(xi,yi,zi)（8）四維表現(xiàn)圖針對(duì)v=v(x,y,z)類型的函數(shù)，引入了三維實(shí)體的四維切片色圖，由函數(shù)slice實(shí)現(xiàn)。slice(X,Y,Z,V,Sx,Sy,Sz)繪制向量Sx,Sy,Sz中的點(diǎn)沿x,y,z方向切片圖，數(shù)組X，Y，Z用來(lái)定義V的坐標(biāo)，V是M×N×P的矩陣。slice(X,Y,Z,V,XI,YI,ZI)繪制沿XI，YI，ZI數(shù)組定義的曲面的通過(guò)容量V的切片圖。slice(V,Sx,Sy,Sz)或slice(V,XI,YI,ZI)默認(rèn)X=1:N,y=1:M,Z=1:Pslice(...,’method’)由method指定使用的插值方法，可以為linear，cubic或nearestH=slice(...)返回處理surface對(duì)象的向量。3.4

Matlab圖形處理的基本技術(shù)1、圖形的控制圖形處理的基本技術(shù)包括圖形控制、圖形標(biāo)注、圖形保持以及子圖的繪制等，前三項(xiàng)除了使用函數(shù)命令之外，還可以直接在生成的圖形界面內(nèi)進(jìn)行對(duì)應(yīng)修改。坐標(biāo)軸控制函數(shù)axis

兩種調(diào)用形式：axis(V)，V是一個(gè)數(shù)組，存儲(chǔ)坐標(biāo)軸的范圍，V=[Xmin,Xmax,Ymin,Ymax,Zmin,Zmax]。axis‘控制字符串’用戶可以通過(guò)選擇不同的控制字符串完成對(duì)坐標(biāo)軸的操作。坐標(biāo)軸縮放函數(shù)zoom

調(diào)用形式：zoom‘控制字符串’對(duì)圖形的縮放不會(huì)改變圖形的基本結(jié)構(gòu)。平面的坐標(biāo)網(wǎng)圖函數(shù)grid調(diào)用形式：gridon/off圖形中繪制/取消坐標(biāo)網(wǎng)格坐標(biāo)軸封閉函數(shù)box調(diào)用形式：boxon/off圖形四周都顯示/常規(guī)顯示坐標(biāo)軸2、圖形的標(biāo)注坐標(biāo)軸標(biāo)注主要函數(shù)有title,xlabel,ylabel,調(diào)用形式為：xlabel(‘標(biāo)注’,’屬性1’,屬性值1,’屬性2’,屬性值2,...)屬性包括字體大小，字體名，字體粗細(xì)等。例如：title(‘正弦函數(shù)’,’FontSize’,12,’FontWeight’,’bold’,’FontName’,’隸書(shū)’)標(biāo)注中的特殊符號(hào)有相應(yīng)的轉(zhuǎn)換字符，例如\alpha表示α，\beta表示β，\gamma表示γ，\delta表示δ等等。顯示控制\bf:黑體，\it:斜體，\rm:標(biāo)準(zhǔn)形式，\fontname{fontname}:字體名\fontsize(fontsize):字體大小。其他符號(hào)如^,{等只需使用”\”來(lái)引導(dǎo)即可，有分?jǐn)?shù)或指數(shù)的形式，相應(yīng)的指數(shù)或分母放到大括號(hào)內(nèi)即可。例如：title(‘\bfy=e^{-x^{2}/2}’);文本標(biāo)注函數(shù)texttext(x,y,’標(biāo)注文本或控制字符串’)x和y表示標(biāo)注的位置，標(biāo)注中可以實(shí)時(shí)調(diào)用返回值為字符串的函數(shù)，如char,num2str等

交互式文本標(biāo)注legend用于解決同一張圖中繪制多條曲線的情況。legend（’標(biāo)注1’,’標(biāo)注2’,....,’定位代號(hào)’）標(biāo)注1，標(biāo)注2等分別對(duì)應(yīng)繪圖過(guò)程按繪制先后順序所生成的曲線。定位代號(hào)為0-4和-1，分別針對(duì)不同的位置，默認(rèn)為1，置于圖形的右上角。③圖形的保持和子圖圖形的保持函數(shù)holdon/off子圖命令subplot(m,n,p)3.5動(dòng)畫(huà)技術(shù)Matlab使科學(xué)計(jì)算余動(dòng)畫(huà)自然結(jié)合，實(shí)現(xiàn)完美的效果，Matlab中有很多動(dòng)畫(huà)的應(yīng)用程序，如travel,truss,lorenz等。Matlab中產(chǎn)生動(dòng)畫(huà)的兩種形式：1、以影片方式產(chǎn)生動(dòng)畫(huà)。這種方式是以圖像的方式預(yù)存多個(gè)畫(huà)面，再將這些畫(huà)面助陣播放，即可得到動(dòng)畫(huà)的效果。這種方法需要事先準(zhǔn)備圖片，無(wú)法進(jìn)行實(shí)時(shí)成像，占內(nèi)存大。產(chǎn)生動(dòng)畫(huà)的方法是使用getframe命令來(lái)抓取圖形作為動(dòng)畫(huà)，每個(gè)畫(huà)面都是以一個(gè)列向量的方式置于存放整個(gè)電影的矩陣M中。然后使用movie(M,k)命令來(lái)播放電影。示例：clear;s=0.2;x1=0;nframes=50;fork=1:nframesx1=x1+s;x=0:0.01:x1;y=sin(x);plot(x,y);axis([02*pi-11]);gridoff;M(k)=getframe;//將當(dāng)前圖形存入Mendmovie(M,3);2、以對(duì)象方式產(chǎn)生動(dòng)畫(huà)。這種方式保持圖形窗口中的大部分對(duì)象，即整個(gè)背景不變，而只更新部分運(yùn)動(dòng)的對(duì)象，以便加快整幅圖像的實(shí)時(shí)生成速度。缺點(diǎn)是無(wú)法產(chǎn)生復(fù)雜的動(dòng)畫(huà)。產(chǎn)生移動(dòng)的動(dòng)畫(huà)效果需要先計(jì)算對(duì)象的新位置，并在新位置上顯示出對(duì)象，然后擦除原位置上的舊對(duì)象，并刷新屏幕。需要用到擦除屬性drasemode(取值有normal,xor,background,none)，一般設(shè)置為xor。然后用drawnow命令刷新屏幕，有時(shí)會(huì)通過(guò)pause函數(shù)進(jìn)行延遲運(yùn)行。a=45;b=2;c=-300;s=0;t=0;forn=1:5000x=t-sign(s)*(sqrt(abs(b*s-c)));y=a-s;plot(x,y,'o');holdon;pause(0.1)s=x;t=y;end示例1：示例2：t=0:pi/20:2*pi;y=exp(sin(t));h=plot(t,y,'YDataSource','y');fork=1:.1:10y=exp(sin(t.*k));refreshdata(h,‘caller’);drawnow;pause(.1)end3.6

Matlab圖形用戶界面設(shè)計(jì)圖形窗口本身是由不同的圖形對(duì)象組成，如坐標(biāo)軸，曲線，曲面或者文字等，Matlab給每一個(gè)圖形對(duì)象都分配了一個(gè)標(biāo)識(shí)符，成為句柄(handle)，可以通過(guò)該句柄對(duì)對(duì)象的屬性進(jìn)行編程控制。直接對(duì)圖形句柄進(jìn)行操作的繪圖方法稱為底層繪圖操作。GUI就是由窗口、菜單、對(duì)話框等各種圖形元素組成的用戶界面。1、創(chuàng)建圖形對(duì)象的底層函數(shù)。（1）h=figure(n)，創(chuàng)建圖形窗口，n為窗口序號(hào)；（2）h=uicontrol(‘property’,value),圖形界面控制，property和value確定控制類型，如BackGroundColor，CallBack，F(xiàn)ontName等（3）h=uimenu(‘property’,value)，創(chuàng)建用戶界面菜單（4）h=axes(‘property’,RECT),創(chuàng)建軸對(duì)象，RECT用于定義軸對(duì)象的位置和大小。（5）h=line(x,y,z)，創(chuàng)建線對(duì)象，二維用x和y。（6）h=patch(x,y,z,c),創(chuàng)建塊對(duì)象，x,y,z定義多邊形，c填充顏色。（7）h=surface(x,y,z,c)創(chuàng)建面對(duì)象，x,y,z定義三維曲面坐標(biāo)，c定義顏色矩陣。（8）h=image(x)顯示圖像，x為圖像矩陣（9）h=text(x,y,’string’)標(biāo)注文字示例1：xdata=[22025;

28245;

88248];ydata=[44420;

84622;

40400];cdata=[1504610;

12579;

23083];hp=patch(xdata,ydata,cdata,'Marker','o','MarkerFaceColor','flat','FaceColor','none')set(hp,‘EdgeColor’,‘g’)或者set(hp,'EdgeColor','none')示例2：t=0:pi/20:2*pi;hline1=plot(t,sin(t),'k');hline2=line(t+.06,sin(t),'LineWidth',4,'Color',[.8.8.8]);2、圖形對(duì)象屬性的操作（1）對(duì)象屬性的獲取。主要命令有：get:獲得句柄圖形對(duì)象的屬性和返回某些對(duì)象的句柄值。get(h)，get(gcf)set:改變圖形對(duì)象的屬性。gcf:當(dāng)前窗口對(duì)象的句柄，對(duì)應(yīng)Figure。h=gcf:返回當(dāng)前窗口對(duì)象的句柄。gca:當(dāng)前軸對(duì)象的句柄，對(duì)應(yīng)Axes。delete(gcf):刪除當(dāng)前窗口的屬性。示例1：h=figure(1);h1=line(2:36,3:37);h3