數(shù)據(jù)和函數(shù)的可視化

第五講數(shù)據(jù)和函數(shù)的可視化——matlab語言豐富的圖形表現(xiàn)方法，使得數(shù)學(xué)計(jì)算結(jié)果可以方便地、多樣性地實(shí)現(xiàn)可視化，這是其它語言所不能比擬的。12/12/20221數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第1頁!5.1引導(dǎo)任何二元實(shí)數(shù)標(biāo)量對（x,y）可用平面上的一個(gè)點(diǎn)表式；任何二元實(shí)數(shù)向量對（x,y）可用平面上的一組點(diǎn)表示。對于離散實(shí)函數(shù)yn=f(xn)，當(dāng)xn以遞增（或遞減）次序取值x=[x1,x2,…,xN]T時(shí)，根據(jù)函數(shù)關(guān)系可求得同樣數(shù)目的yn，y=[y1,y2,…,yN]T。當(dāng)把這向量對用直角坐標(biāo)中的點(diǎn)序列圖示時(shí)，就實(shí)現(xiàn)了離散函數(shù)的可視化。5.1.1離散數(shù)據(jù)和離散函數(shù)的可視化12/12/20222數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第2頁!例5.1-1：離散函數(shù)y=|n|的可視化n=(-10:10)';y=abs(n);plot(y,'r.')figure(2),plot(n,y,'b')12/12/20223數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第3頁!5.1.2連續(xù)函數(shù)的可視化與離散函數(shù)可視化一樣，進(jìn)行連續(xù)函數(shù)可視化也必須先在一組離散自變量上計(jì)算相應(yīng)的函數(shù)值，并把這一組“數(shù)據(jù)對”用點(diǎn)圖示。但這些離散的點(diǎn)不能表現(xiàn)函數(shù)的連續(xù)性。為了進(jìn)一步表示離散點(diǎn)之間的函數(shù)情況，有兩種常用處理方法：對區(qū)間進(jìn)行更細(xì)的分割，計(jì)算更多的點(diǎn)，去近似表現(xiàn)函數(shù)的連續(xù)變化；把兩點(diǎn)用直線連接，近似表現(xiàn)兩點(diǎn)間的（一般是非線性的）函數(shù)性狀。在MATLAB中，以上兩種表現(xiàn)方法都可以采用。注意：倘若自變量的采樣點(diǎn)數(shù)不足夠多，則無論哪種方法都不能真實(shí)地反映原函數(shù)。12/12/20224數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第4頁!N=9; t=0:2*pi/N:2*pi; x=sin(t);y=cos(t); tt=reshape(t,2,(N+1)/2); tt=flipud(tt); tt=tt(:); xx=sin(tt);yy=cos(tt);subplot(1,2,1),plot(x,y) title('(1)正常排序圖形'),axisequaloff,shg subplot(1,2,2),plot(xx,yy) title('(2)非正常排序圖形'),axisequaloff,shg

例5.1-3：繪制奇數(shù)正多邊形及圓。開始點(diǎn)開始點(diǎn)tt=01.39632.79254.18885.58510.69812.09443.49074.88696.2832tt=0.69812.09443.49074.88696.283201.39632.79254.18885.5851tt=0.698102.09441.39633.49072.79254.88694.18886.28325.585112/12/20225數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第5頁!1.基本調(diào)用格式plot(x,y,‘s')(1)一維數(shù)組(x,y)指定采樣點(diǎn)的橫坐標(biāo)和縱坐標(biāo)。(2)’

s’是字符串，用來設(shè)定“離散點(diǎn)形”或/和“連續(xù)線型”，也指定“點(diǎn)線色彩”。顏色字符串用英文單詞的前1~3個(gè)字母，如yellow—yel（或y或ye）表示等。(3)若沒有第三個(gè)輸入量，將采用默認(rèn)設(shè)置：藍(lán)色細(xì)實(shí)線。輸入量(x,y,‘s’)稱為平面繪線三元組。12/12/20226數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第6頁!表5.2-4連續(xù)線型允許設(shè)置值符號-:-.--含義實(shí)線虛線點(diǎn)劃線雙劃線表5.2-4點(diǎn)線色彩允許設(shè)置值符號bgrcmykw含義藍(lán)綠紅青品紅黃黑白如：s=‘r:’表示紅色虛線t=0:10;plot(t,'pr:','MarkerSize',18)12/12/20227數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第7頁!plot(x1,y1,’s1’,x2,y2,’s2’)在此格式中，每個(gè)繪線“三元組”(X,Y,’s’)的結(jié)構(gòu)和作用，與plot(X,Y,’s’)相同。不同的“三元組”之間沒有約束關(guān)系。（3）單輸入量繪線plot(Y)（2）多三元組繪制多條曲線Y是實(shí)向量時(shí)，以該向量元素的下標(biāo)為橫坐標(biāo)、元素值為縱坐標(biāo)畫出一條連續(xù)曲線。

Y是實(shí)矩陣時(shí)，則按列繪制每列元素值相對其下標(biāo)的曲線。圖中曲線數(shù)等于Y陣列數(shù)。

Y是復(fù)數(shù)矩陣時(shí)，則按列分別以元素實(shí)部和虛部為橫、縱坐標(biāo)繪制多條曲線12/12/20228數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第8頁!例5.2-1：二維曲線繪圖指令演示。Clft=(0:pi/50:2*pi)';k=0.4:0.1:1; Y=cos(t)*k; %是矩陣subplot(1,2,1),plot(t,Y,'LineWidth',1.5)title('Byplot(t,Y)'),xlabel('t')subplot(1,2,2),plot(Y,'LineWidth',1.5)title('Byplot(Y)'),xlabel('rowsubscriptofY')12/12/20229數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第9頁!5.2.2坐標(biāo)控制和圖形標(biāo)識表5.2-6常用的坐標(biāo)控制指令坐標(biāo)軸控制方式、取向和范圍坐標(biāo)軸的高寬比指令含義指令含義axisauto使用缺省設(shè)置axisequal縱、橫軸采用等長刻度axismanual使當(dāng)前坐標(biāo)范圍不變axisfill在manual方式下起作用，使坐標(biāo)充滿整個(gè)繪圖區(qū)axisoff取消軸背景axisimage縱、橫軸采用等長刻度,且坐標(biāo)框緊貼數(shù)據(jù)范圍axison使用軸背景axisnormal缺省矩形坐標(biāo)系axisij矩陣式坐標(biāo)，原點(diǎn)在左上方axissquare產(chǎn)生正方形坐標(biāo)系axisxy普通直角坐標(biāo)，原點(diǎn)在左下方axistight把數(shù)據(jù)范圍直接設(shè)為坐標(biāo)范圍Axis(V)V=[x1,x2,y1,y2];V=[x1,x2,y1,y2,z1,z2];人工設(shè)定坐標(biāo)范圍。社定植：二維，4個(gè)；三維，6個(gè)axisvis3d保持高寬比不變，用于三維旋轉(zhuǎn)時(shí)避免圖形大小變化說明：坐標(biāo)范圍設(shè)定向量V中的元素必須服從：x1<x2,y1<y2,z1<z2。V的元素值允許取inf或-inf，那意味著上限或下限是自動(dòng)產(chǎn)生的，即坐標(biāo)范圍半自動(dòng)確定。1.坐標(biāo)軸的控制12/12/202210數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第10頁!2.分格線和坐標(biāo)框grid 是否畫分格線的雙向切換指令（使當(dāng)前分格線狀態(tài)翻轉(zhuǎn)）gridon畫出分格線gridoff不畫分格線box坐標(biāo)形式在封閉式和開啟式之間切換指令boxon使當(dāng)前坐標(biāo)呈封閉形式boxoff使當(dāng)前坐標(biāo)呈開啟形式默認(rèn)設(shè)置：不畫出分格線，所畫坐標(biāo)呈封閉形式。12/12/202211數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第11頁!4.標(biāo)識指令中字符的精細(xì)控制表5.2-7圖形標(biāo)識用的希臘字母指令字符指令字符指令字符指令字符\alphaα\thetaθ\XiЕ\phiφ\betaβ\ThetaΘ\piπ\(zhòng)PhiΦ\gammaγ\iotaι\PiΠ\chiχ\GammaГ\kappaκ\rhoρ\psiψ\deltaδ\lambdaλ\sigmaσ\PsiΨ\DeltaΔ\LambdaΛ\SigmaΣ\omegaω\epsilonε\muμ\tauτ\OmegaΩ\zetaζ\Nuν\upsilonυ\etaη\xiξ\UpsilonY指令效果指令效果指令效果'sin\bezt'sinβ'\zeta\omega'ξω'\itA{\in}R^{m\timesn}'A

Rm×n12/12/202212數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第12頁!4.標(biāo)識指令中字符的精細(xì)控制若想設(shè)置上下標(biāo)，對字體和字體大小進(jìn)行控制，須在被控制字符前使用5.2-9,5.2-10的指令表5.2-9上下標(biāo)的控制指令

分類指令arg取值舉例示例指令效果上標(biāo)^{arg}任何合法字符'\ite^{-t}sint'e-tsint下標(biāo)_{arg}任何合法字符'\x～{\chi}_{alpha}^{2}(3)'x

12/12/202213數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第13頁!表5.2-2離散數(shù)據(jù)點(diǎn)形允許設(shè)置值符號含義符號含義.實(shí)心點(diǎn)d菱形符diamond+十字符h六角星hexagram*米字符o空心圓圈^朝上三角符p五角星pentagram<朝坐三角符s方塊符square>朝右三角符x叉字符v朝下三角符1.基本調(diào)用格式plot(x,y,‘s')平面繪線三元組表5.2-4連續(xù)線型允許設(shè)置值符號-:-.--含義實(shí)線虛線點(diǎn)劃線雙劃線表5.2-4點(diǎn)線色彩允許設(shè)置值符號bgrcmykw含義藍(lán)綠紅青品紅黃黑白Review12/12/202214數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第14頁!例5.2-4：標(biāo)識實(shí)驗(yàn)。clf;t=0:pi/50:2*pi;y=sin(t);plot(t,y);axis([0,2*pi,-1.2,1.2])text(pi/2,1,'\fontsize{16}\leftarrow\itsin(t)\fontname{隸書}極大值')title('y=sin(t)')xlabel('t')ylabel('y')

12/12/202215數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第15頁!5.2.3多次疊繪、雙縱坐標(biāo)和多子圖前面已介紹過plot在同一次調(diào)用中畫多條曲線的功能。實(shí)際應(yīng)用中，還會(huì)遇到在已經(jīng)存在的圖上再繪制一條或多條曲線的情況。為此，MATLAB提供了以下指令：1.多次疊繪holdon使當(dāng)前軸及圖形保持而不被刷新，準(zhǔn)備接受此后將繪制的新曲線holdoff 使當(dāng)前軸及圖形不再具備不刷新的性質(zhì)hold 當(dāng)前圖形是否具備刷新性質(zhì)的雙向切換開關(guān)12/12/202216數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第16頁!2.雙縱坐標(biāo)圖在實(shí)際應(yīng)用中常常提出這樣一種需求：把同一自變量的兩個(gè)不同量綱、不同數(shù)量級的函數(shù)量的變化繪制在同一張圖上。為此，MATLAB5.x新增了以下指令。plotyy(X1,Y1,X2,Y2) %以左、右不同縱軸繪制X1-Y1，X2-Y2兩條曲線plotyy(X1,Y1,X2,Y2,FUN)%以左、右不同縱軸把X1-Y1，X2-Y2繪制成繪圖函數(shù)名FUN指定形式的兩條曲線。plotyy(X1,Y1,X2,Y2,FUN1,FUN2)

%以左,右不同縱軸把X1-Y1，X2-Y2繪制成繪圖函數(shù)名FUN1，F(xiàn)UN2指定的不同形式的兩條曲線12/12/202217數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第17頁!3.多子圖MATLAB允許用戶在同一圖形窗中布置幾幅獨(dú)立的子圖，具體指令是：subplot(m,n,k) 使(m×n)幅子圖中的第k幅成為當(dāng)前圖。subplot(‘position’,[leftbottomwidthheight])

在指定位置上開辟子圖，并成為當(dāng)前圖。12/12/202218數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第18頁!例5.2-8演示subplot指令對圖形窗的分割。clf;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])subplot('position',[0.2,0.1,0.6,0.40])

plot(t,y12,'b-',t,[y1,-y1],'r:')axis([0,pi,-1,1])

12/12/202219數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第19頁!clfx=0:0.01:1;y=(x+2).^x-2;plot(x,y);gridon

例5.2-9采用圖解法求的解。（1）繪制的曲線。[x,y]=ginput(1);（2）放大圖形，用鼠標(biāo)獲取曲線與橫軸交點(diǎn)坐標(biāo)值。formatlongx,y x=0.698286930399556y=-3.381823762854459e-00912/12/202220數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第20頁!gtextgtext(arg),

用鼠標(biāo)把字符串或字符串異質(zhì)(元胞)數(shù)組arg放置到圖形上。該指令對二維、三維圖形都適用。具體操作方法：指令運(yùn)行后，當(dāng)前圖形窗自動(dòng)由后臺轉(zhuǎn)為前臺，鼠標(biāo)光標(biāo)變?yōu)槭植?；移?dòng)鼠標(biāo)，使十字叉移到希望的位置；點(diǎn)動(dòng)鼠標(biāo)右鍵，arg所承載的字符將被放在緊靠十字叉點(diǎn)的“象限”位置。A={‘My',‘God';‘It’s',‘here!'}

gtext(A)12/12/202221數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第21頁!5.3三維曲線和曲面圖plot3(X,Y,Z’s’)plot3(X1,Y1,Z1,'s1',X2,Y2,Z2,'s2',…)5.3.1三維線圖指令plot3X,Y,Z是同維向量時(shí),則繪制以X,Y,Z元素為x,y,z坐標(biāo)的三維曲線．X,Y,Z是同維矩陣時(shí)，則以X,Y,Z對應(yīng)列元素為x,y,z坐標(biāo)分別繪制曲線．曲線條數(shù)等于矩陣列數(shù)s,s1,s2的意義與二維情況相同，是主要用于指定線型，色彩，數(shù)據(jù)點(diǎn)形的選項(xiàng)字符串。12/12/202222數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第22頁!例5.3-1三維曲線繪圖t=(0:0.02:2)*pi;x=sin(t);y=cos(t);z=cos(2*t);plot3(x,y,z,'r-',x,y,z,'bd'),view([-82,58]),boxon,legend('鏈','寶石')12/12/202223數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第23頁!三維網(wǎng)線圖曲面圖作圖要領(lǐng)x——1×m數(shù)組,y——1×n數(shù)組生成X,Y格點(diǎn)坐標(biāo)——[X,Y]=meshgrid(x,y)X由x（1×m）按行放置，再縱向擴(kuò)展n行Y由y（1×n）按列放置，再橫向擴(kuò)展m行結(jié)果：X和Y均成為n×m矩陣表達(dá)式點(diǎn)運(yùn)算——Z=X.^2+Y.^2X,Y,Z是n×m的數(shù)組，維數(shù)可任定默認(rèn)方位角：-37.5o，俯角30o。12/12/202224數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第24頁!clf,x=-4:4;y=x;[X,Y]=meshgrid(x,y);%生成x-y坐標(biāo)"格點(diǎn)"矩陣Z=X.^2+Y.^2;surf(X,Y,Z);%繪曲面圖holdon;colormap(hot)stem3(X,Y,Z,‘bo’);%表現(xiàn)在格點(diǎn)上計(jì)算函數(shù)值例5.3-2：用曲面圖表現(xiàn)函數(shù)figure(2),mesh(X,Y,Z)%繪網(wǎng)線圖holdon;colormap(hot)12/12/202225數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第25頁!5.3.3曲面/網(wǎng)線圖的精細(xì)修飾Colormap(CM)%設(shè)置當(dāng)前圖形窗的著色色圖為CM。2.色彩控制色圖為(m×3)矩陣，它的每一行是RGB三元組。常用色圖矩陣見表5.3-1。這些預(yù)定義色圖矩陣是由[0,1]區(qū)間數(shù)據(jù)組成的(64×3)矩陣.Jet是默認(rèn)色圖。12/12/202226數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第26頁!【例5.3-3】三種濃淡處理方式比較。clf x=-4:4;y=x;[X,Y]=meshgrid(x,y);Z=X.^2+Y.^2;surf(X,Y,Z)colormap(jet)subplot(1,3,1),surf(Z),axisoff%默認(rèn)shadingfacetedsubplot(1,3,2),surf(Z),axisoff,shadingflatsubplot(1,3,3),surf(Z),axisoff,shadinginterpset(gcf,'Color','w') 12/12/202227數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第27頁!【例5.3-4】半透明的表面圖。clfsurf(peaks)shadinginterpalpha(0.7) colormap(summer)

alpha(v)：0<v<1用于曲面圖形，0表示完全透明，1表示不透明。4.透明控制12/12/202228數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第28頁!6.照明模式lighting該指令只有在light指令執(zhí)行后才起作用。lightingoptions指令功能flat入射光均勻的灑在圖形對象的每個(gè)面上（缺省模式）gouraud先對頂點(diǎn)顏色進(jìn)行插補(bǔ)，再對頂點(diǎn)勾畫的面色進(jìn)行插補(bǔ)，用于曲面表現(xiàn)。phong對頂點(diǎn)處法線進(jìn)行插值，再計(jì)算各像素的反光，表現(xiàn)效果最好，但費(fèi)時(shí)。none關(guān)閉所有的光源。12/12/202229數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第29頁!clf;[X,Y,Z]=sphere(40);colormap(jet)subplot(1,2,1),surf(X,Y,Z),axisequaloff,shadinginterplight('position',[0-101.5],'style','infinite')lightingphongmaterialshinysubplot(1,2,2),surf(X,Y,Z,-Z),axisequaloff,shadingflatlight;lightingflatlight('position',[-1,-1,-2],'color','y')light('position',[-1,0.5,1],'style','local','color','w')set(gcf,'Color','w') 【例5.3-5】燈光、照明、材質(zhì)指令所表現(xiàn)的圖形。12/12/202230數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第30頁!例5.3-6透視演示[X0,Y0,Z0]=sphere(30);X=2*X0;Y=2*Y0;Z=2*Z0;surf(X0,Y0,Z0);shadinginterp;holdon;mesh(X,Y,Z),colormap(hot),holdoffhiddenoff;axisequal,axisoff12/12/202231數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第31頁!3.裁切由NaN處理的圖形不會(huì)產(chǎn)生切面，為了看清圖形而需要表現(xiàn)切面，那么應(yīng)該把被切部分強(qiáng)制為零．例5.3-8表現(xiàn)切面x=[-8:0.1:8];y=x;[X,Y]=meshgrid(x,y);ZZ=X.^2+Y.^2-40;ii=find(abs(X)>6|abs(Y)>6);%確定超出[-6,6]范圍的格點(diǎn)下標(biāo)ZZ(ii)=zeros(size(ii));%強(qiáng)制為0surf(X,Y,ZZ),shadinginterp;colormap(copper)light('position',[0,-15,1]);lightingphongmaterial([0.8,0.8,0.5,10,0.5])12/12/202232數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第32頁!closeall;clear;[X,Y,Z]=peaks(40);n=6;%等高線分級數(shù)subplot(1,2,1),pcolor(X,Y,Z)shadinginterpzmax=max(max(Z));zmin=min(min(Z));caxis([zmin,zmax]) ColorbarholdonC=contour(X,Y,Z,n,'k:'); clabel(C)holdoffsubplot(1,2,2)[C,h]=contourf(X,Y,Z,n,'k:'); clabel(C,h)colormap(cool)set(gcf,'Color','w')

例5.4-1

二維半圖指令pcolor,contour,contourf應(yīng)用。12/12/202233數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第33頁!5.4.2四維表現(xiàn)clfx=3*pi*(-1:1/15:1);y=x;[X,Y]=meshgrid(x,y);R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;[dzdx,dzdy]=gradient(Z);dzdr=sqrt(dzdx.^2+dzdy.^2); dz2=del2(Z);subplot(1,2,1),surf(X,Y,Z,abs(dzdr)) shadingfaceted;colorbar('SouthOutside') brighten(0.6); colormaphsvtitle('No.1surf(X,Y,Z,abs(dzdr))')subplot(1,2,2);surf(X,Y,Z,abs(dz2)) shadingfaceted colorbar('NorthOutside') title('No.2surf(X,Y,Z,abs(dz2))')

例5.4-2

用顏色表現(xiàn)z=f(x,y)

函數(shù)的其他特征（如梯度，曲率等）。1.準(zhǔn)四維表現(xiàn)12/12/202234數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第34頁!5.4.3動(dòng)態(tài)圖形1.

彗星狀軌跡圖et(x,y,p)二維彗星軌線et(x,y,z,p)三維彗星軌線p指定彗星體的長度p*length(y)，缺省的p值為0.1。三維圖形中彗星體的長度p*length(z)shg;n=50;t=n*pi*(0:0.000005:1);x=sin(t);y=cos(t);plot(x,y,'g');axissquareholdonet(x,y,0.001)holdoff

例5.4-4簡單二維示例。12/12/202235數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第35頁!2.色圖的變幻spinmap(t,inc)>分別用t,inc(缺省為2)控制色圖旋轉(zhuǎn)的時(shí)間和快慢ezsurf('x*y','circ');shadingflat;view([-18,28])C=summer;CC=[C;flipud(C)];colormap(CC)spinmap(30,4)

例5.4-5指令spinmap的使用。12/12/202236數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第36頁!例5.4-6三維圖形的影片動(dòng)畫。clf,x=3*pi*(-1:0.05:1);y=x;[X,Y]=meshgrid(x,y);R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;h=surf(X,Y,Z);colormap(jet);axisoff;n=12;mmm=moviein(n);fori=1:nrotate(h,[001],25);%使圖形繞Z軸旋轉(zhuǎn)25°/每次mmm(:,i)=getframe;%捕獲畫面endcloseshg;axisoff;movie(mmm,20,10)%以每秒10幀速度重復(fù)播放20次12/12/202237數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第37頁!例模擬布朗運(yùn)動(dòng)。程序如下：n=30;%指定布朗運(yùn)動(dòng)的點(diǎn)數(shù)s=0.02;%指定溫度或速率%產(chǎn)生n個(gè)隨機(jī)點(diǎn)(x,y)，處于-0.5到0.5之間x=rand(n,1)-0.5;y=rand(n,1)-0.5;h=plot(x,y,'.');%繪制隨機(jī)點(diǎn)axis([-11-11]);axissquaregridoffset(h,'EraseMode','Xor','MarkerSize',20);%設(shè)置擦除模式為Xor%循環(huán)5000次，產(chǎn)生動(dòng)畫效果fori=linspace(1,10,5000)drawnowx=x+s*randn(n,1);%在坐標(biāo)點(diǎn)附近添加隨機(jī)噪聲y=y+s*randn(n,1);set(h,'XData',x,'YData',y);%通過改變數(shù)據(jù)屬性來重新繪圖end12/12/202238數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第38頁!f=anim_zzy1(2,450);image(f.cdata),axisoff

12/12/202239數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第39頁!5.5圖形窗功能簡介圖5.5-2個(gè)性化的圖形（1）t=(0:pi/100:pi)'; y1=sin(t)*[1,-1]; y2=sin(t).*sin(9*t); t3=pi*(0:9)/9;y3=zeros(size(t3));plot(t,y1,t,y2,t3,y3)

（2）【例5.5-1】（12/12/202240數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第40頁!5.5圖形窗功能簡介圖5.5-2個(gè)性化的圖形【例5.5-1】（5）圖5.5-6調(diào)制曲線的設(shè)置示例（6）12/12/202241數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第41頁!第四章習(xí)題作業(yè)存在的問題十八題問題出現(xiàn)在加載數(shù)據(jù)x,y上，用Matlab6.5打不開prob_data401.mat，需經(jīng)過轉(zhuǎn)換。加載指令loadprob_data418.mat,loadprob_data418,檢查文件名?？捎媚M數(shù)據(jù)，注意x,y數(shù)據(jù)長度必須相等，指令polyfit，不是ployfit。polyfit是數(shù)值運(yùn)算指令，不是符號運(yùn)算指令。再有就是繪圖時(shí)區(qū)域要與x取值一致,x=-1:0.1:4,有同學(xué)取xx=0:0.01:1。xx=-1:0.01:4,yy=polyval(P,xx),plot(xx,yy,x,y,'*r')12/12/202242數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第42頁!例5.1-2：用圖形表示連續(xù)調(diào)制波形y=sin(t)sin(9t)t1=(0:11)/11*pi; t2=(0:400)/400*pi;

t3=(0:50)/50*pi;

y1=sin(t1).*sin(9*t1); y2=sin(t2).*sin(9*t2);y3=sin(t3).*sin(9*t3);subplot(2,2,1),plot(t1,y1,'r.') axis([0,pi,-1,1]),title('(1)點(diǎn)過少的離散圖形')subplot(2,2,2),plot(t1,y1,t1,y1,‘r.’)%紅色部分畫折線，其余畫點(diǎn)axis([0,pi,-1,1]),title('(2)點(diǎn)過少的連續(xù)圖形')subplot(2,2,3),plot(t2,y2,'r.') axis([0,pi,-1,1]),title('(3)點(diǎn)密集的離散圖形')subplot(2,2,4),plot(t3,y3) axis([0,pi,-1,1]),title('(4)點(diǎn)足夠的連續(xù)圖形')12/12/202243數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第43頁!5.2二維曲線和圖形plot——最基本的二維圖形指令plot的功能：plot命令自動(dòng)打開一個(gè)圖形窗口Figure用直線連接相鄰兩數(shù)據(jù)點(diǎn)來繪制圖形根據(jù)圖形坐標(biāo)大小自動(dòng)縮擴(kuò)坐標(biāo)軸,自動(dòng)標(biāo)注數(shù)據(jù)標(biāo)尺及單位標(biāo)注5.2.1二維曲線繪制的基本指令plot12/12/202244數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第44頁!表5.2-2離散數(shù)據(jù)點(diǎn)形允許設(shè)置值符號含義符號含義.實(shí)心點(diǎn)d菱形符diamond+十字符h六角星hexagram*米字符o空心圓圈^朝上三角符p五角星pentagram<朝坐三角符s方塊符square>朝右三角符x叉字符v朝下三角符演示數(shù)據(jù)點(diǎn)型的變化t=0:10;plot(t,'p','MarkerSize',18)12/12/202245數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第45頁!2.plot的衍生調(diào)用格式plot(X,Y)指令采用細(xì)實(shí)線繪制多條彩色曲線，各曲線的顏色順序見表5.2-4.（1）單色或多色繪制多條曲線plot(X,Y,’s’)X、Y是均為（m×n）矩陣時(shí)，則以X、Y對應(yīng)列元素為橫、縱坐標(biāo)分別繪制n條曲線。X、Y之一是一維數(shù)組，且數(shù)組長度與另一個(gè)矩陣輸入量的“行數(shù)”（或“列數(shù)”）相等時(shí)，將繪制出“列數(shù)”（或“行數(shù)”）條曲線。’s’

用來指定多條曲線用同一色彩繪制。X、Y均是一維數(shù)組時(shí)，就是基本調(diào)用格式。plot(X,Y,'p-','MarkerSize',18)X=[114;225;336;]Y=[147;258;369;]12/12/202246數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第46頁!3.

plot的屬性可控調(diào)用格式plot(x,y,'s','PropertyName',PropertyValue,…)含義屬性名屬性值說明點(diǎn)、線顏色Color[Vr,Vg,Vb],取[0,1]默認(rèn)為'b'線型LineStyle'-',':','-.','-'默認(rèn)為實(shí)線線寬LineWidth正實(shí)數(shù)默認(rèn)為0.5點(diǎn)形狀Marker'd','+',…可通過s設(shè)置點(diǎn)大小MarkerSize正實(shí)數(shù)默認(rèn)為6.0plot(x,y,'LineStyle',':','LineWidth',8.5)12/12/202247數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第47頁!例5.2-2：用圖形表示連續(xù)調(diào)制波形y=sin(t)sin(9t)極其包絡(luò)線。t=(0:pi/100:pi)';y1=sin(t)*[1,-1];y2=sin(t).*sin(9*t);

t3=pi*(0:9)/9;y3=sin(t3).*sin(9*t3);plot(t,y1,'r:',t,y2,'-bo')holdon plot(t3,y3,'s','MarkerSize',10,'MarkerEdgeColor',[0,1,0],'MarkerFaceColor',[1,0.8,0]) axis([0,pi,-1,1])holdoff 12/12/202248數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第48頁!例5.2-3觀察各種軸控制指令的影響。演示采用長軸為3.25，短軸為1.15的橢圓。注意：采用多子圖表現(xiàn)時(shí)，圖形形狀不僅受“控制指令”影響，而且受整個(gè)圖面“寬高比”及“子圖數(shù)目”的影響。想準(zhǔn)確體會(huì)控制指令的影響，可在全圖狀態(tài)下進(jìn)行觀察。t=0:2*pi/99:2*pi;x=1.15*cos(t);y=3.25*sin(t); subplot(2,3,1),plot(x,y),axisnormal,gridon,title('NormalandGridon')subplot(2,3,2),plot(x,y),axisequal,gridon,title('Equal')subplot(2,3,3),plot(x,y),axissquare,gridon,title('Square')subplot(2,3,4),plot(x,y),axisimage,boxoff,title('ImageandBoxoff')subplot(2,3,5),plot(x,y),axisimagefill,boxofftitle('ImageandFill')subplot(2,3,6),plot(x,y),axistight,boxoff,title('Tight')12/12/202249數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第49頁!3.圖形標(biāo)識圖形標(biāo)識包括：圖名（Title）、坐標(biāo)軸名（Label）、圖形注釋（Text）和圖例（Legend）。標(biāo)識指令的最簡捷使用格式如下：title(S) 書寫圖名xlabel(S) 橫坐標(biāo)軸名ylabel(S) 縱坐標(biāo)軸名legend(S1,S2,…) 繪制曲線所用線型、色彩或數(shù)據(jù)點(diǎn)形圖例text(xt,yt,S) 在圖面（xt,yt）坐標(biāo)處書寫字符注釋12/12/202250數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第50頁!表5.2-8圖形標(biāo)識用的其他特殊字符指令字符指令字符指令字符指令字符指令字符\approx≈\propto∝\exists\cap∩\downarrow↓\cong\sim～\forall\cup∪\leftarrow←\div÷\times×\in\subset\leftrightarrow\equiv≡\oplus+\infty∞\subseteq\rightarrow→\geq≥\oslash\perp⊥\supset\uparrow↑\leq≤\otimes×\prime’\supseteq\circ○\neq≠\int\cdot·\Im\bullet●\pm±\partial\ldots…\Re\copyrightc4.標(biāo)識指令中字符的精細(xì)控制12/12/202251數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第51頁!4.標(biāo)識指令中字符的精細(xì)控制表4字體式樣設(shè)置規(guī)則

x字樣指令arg取值舉例示例指令效果名稱\fontname{arg}arial;courier;roman;宋體；隸書；黑體'\fontname{courier}Example1''\fontname{隸書}范例2'Example1范例2風(fēng)格\argbf（黑體）it（斜體1）bf（斜體2）bf（正體）'\bfExample3''\itExample4'Example3Example4大小\fontsize{arg}任何合法字符'\x～{\chi}_{alpha}^{2}(3)'x

12/12/202252數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第52頁!2坐標(biāo)控制和圖形標(biāo)識axisequal;axissquare;axisoff;axison(1).坐標(biāo)軸的控制grid 是否畫分格線的雙向切換指令（使當(dāng)前分格線狀態(tài)翻轉(zhuǎn)）gridon畫出分格線gridoff不畫分格線(2).分格線(3).圖形標(biāo)識title(S) 書寫圖名xlabel(S) 橫坐標(biāo)軸名ylabel(S) 縱坐標(biāo)軸名legend(S1,S2,…) 繪制曲線所用線型、色彩或數(shù)據(jù)點(diǎn)形圖例text(xt,yt,S) 在圖面（xt,yt）坐標(biāo)處書寫字符注釋Review12/12/202253數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第53頁!例5.2-5通過繪制二階系統(tǒng)階躍響應(yīng)，綜合演示圖形標(biāo)識。clf;t=6*pi*(0:100)/100;y=1-exp(-0.3*t).*cos(0.7*t);plot(t,y,'r-','LineWidth',3);holdontt=t(find(abs(y-1)>0.05));ts=max(tt);plot(ts,0.95,'bo','MarkerSize',10);holdoff;axis([-inf,6*pi,0.6,inf])set(gca,'Xtick',[2*pi,4*pi,6*pi],'Ytick',[0.95,1,1.05,max(y)]) set(gca,'XtickLabel',{'2*pi';'4*pi';'6*pi'})set(gca,'YtickLabel',{'0.95';'1';'1.05';'max(y)'});gridontext(13.5,1.2,'\fontsize{12}{\alpha}=0.3')text(13.5,1.1,'\fontsize{12}{\omega}=0.7')cell_string{1}='\fontsize{12}\uparrow'; cell_string{2}='\fontsize{16}\fontname{隸書}鎮(zhèn)定時(shí)間';cell_string{3}='\fontsize{6}';cell_string{4}=['\fontsize{14}\rmt_{s}='num2str(ts)];text(ts,0.85,cell_string,'Color','b','HorizontalAlignment','Center')title('\fontsize{14}\ity=1-e^{-\alphat}cos{\omegat}') xlabel('\fontsize{14}\bft\rightarrow') ylabel('\fontsize{14}\bfy\rightarrow')

12/12/202254數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第54頁!例5.2-6利用hold繪制離散信號通過零階保持器后產(chǎn)生的波形。t=2*pi*(0:20)/20;y=cos(t).*exp(-0.4*t);stem(t,y,'g','Color','k');holdonstairs(t,y,':r','LineWidth',3)holdofflegend('\fontsize{14}\itstem','\fontsize{14}\itstairs')boxon12/12/202255數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第55頁!clf;dx=0.1;x=0:dx:4;y=x.*sin(x);s=cumtrapz(y)*dx; a=plotyy(x,y,x,s,'stem','plot');text(0.5,1.5,'\fontsize{14}\ity=xsinx')sint='{\fontsize{16}\int_{\fontsize{8}0}^{x}}';ss=['\fontsize{14}\its=',sint,'\fontsize{14}\itxsinxdx'];text(2.5,3.5,ss)set(get(a(1),'Ylabel'),'String','被積函數(shù)\ity=xsinx')set(get(a(2),'Ylabel'),'String',ss) xlabel('x')例5.2-7用雙縱坐標(biāo)畫出函數(shù)和積分在[0,4]區(qū)間上的曲線。12/12/202256數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第56頁!subplot('position',[leftbottomwidthheight])例：subplot('position',[0.50.50.450.45])subplot('position',[0.050.50.40.45])subplot('position',[0.20.050.60.4])12/12/202257數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第57頁!5.2.4獲取二維圖形數(shù)據(jù)的指令ginput[x,y]=ginput(n)

用鼠標(biāo)從二維圖形上獲取n個(gè)點(diǎn)的數(shù)據(jù)坐標(biāo)(x,y)。該指令僅適用于二維圖形。具體操作方法：指令運(yùn)行后，會(huì)把當(dāng)前圖形從后臺調(diào)到前臺，同時(shí)鼠標(biāo)光標(biāo)變?yōu)槭植妫挥脩艨梢苿?dòng)鼠標(biāo)，使十字叉移到待取坐標(biāo)點(diǎn)；點(diǎn)動(dòng)鼠標(biāo)左鍵，便獲得該點(diǎn)數(shù)據(jù)；此后，用同樣的方法，獲取其余點(diǎn)的數(shù)據(jù)；當(dāng)n個(gè)點(diǎn)的數(shù)據(jù)全部取到后，圖形窗便退回后臺，機(jī)器回到ginput執(zhí)行前的環(huán)境。12/12/202258數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第58頁!交互式圖形指令在MATLAB中還有若干與鼠標(biāo)有關(guān)的圖形操作指令。它們是：ginputgtextlegend12/12/202259數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第59頁!legendlegend(arg,pos) 在指定位置建立圖例legendoff 擦除當(dāng)前圖上的圖例說明：輸入宗量arg是圖例中的說明文字注釋。輸入宗量pos是圖例在圖上位置的指定符，它可取下表中的值。Legend在圖形窗中產(chǎn)生后，可用鼠標(biāo)對其進(jìn)行拖拉，即把鼠標(biāo)光標(biāo)移到圖例上，按住鼠標(biāo)左鍵；圖例將隨鼠標(biāo)移動(dòng)，直到滿意位置；放開按鍵便完成操作。pos取值01234-1圖例位置自動(dòng)取最佳位右上角(缺省值)左上角左下角右下角圖右側(cè)12/12/202260數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第60頁!二維圖形的所有基本特性對三維圖形全都適用。定義三維坐標(biāo)軸大小:

axis([xminxmaxyminymaxzminzmax])gridon(off)繪制三維網(wǎng)格text(x,y,z,‘string’)三維圖形標(biāo)注子圖和多窗口也可以用到三維圖形中12/12/202261數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第61頁!5.3.2三維曲面/網(wǎng)線圖需要繪圖數(shù)據(jù)的準(zhǔn)備，可實(shí)現(xiàn)三維圖形的色彩，明暗，光照和視點(diǎn)處理．1.三維圖形的數(shù)據(jù)準(zhǔn)備1)確定自變量x,y的取值范圍和取值間距．如x=x1:dx:x2;y=y1:dy:y22)構(gòu)成xy平面上的自變量格點(diǎn)矩陣:[X,Y]=meshgrid(x,y);3)計(jì)算在自變量采樣格點(diǎn)上的函數(shù)值，即Z=f(X,Y)12/12/202262數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第62頁!2.繪制曲面/網(wǎng)線圖基本指令格式mesh(Z)>以Z矩陣列，行下標(biāo)為x,y軸自變量，畫網(wǎng)線圖．mesh(X,Y,Z)>最常用的網(wǎng)線圖調(diào)用格式mesh(X,Y,Z,C)>最完整的格式，畫由C指定用色的網(wǎng)線圖surf(Z)>以Z矩陣列，行下標(biāo)為x,y軸自變量，畫曲面圖．surf(X,Y,Z)>最常用的曲面圖調(diào)用格式surf(X,Y,Z,C)>最完整調(diào)用格式，畫由C指定用色曲面圖．●在最完整的調(diào)用格式中,四個(gè)輸入宗量都是維數(shù)相同的矩陣。X,Y是自變量格點(diǎn)矩陣,Z是格點(diǎn)上的函數(shù)矩陣。Ｃ指定點(diǎn)的用色,可以缺省,此時(shí)取C=Z?！駟屋斎胱诹扛袷嚼L圖時(shí),把Z矩陣的列下標(biāo)當(dāng)作x坐標(biāo)軸的自變量,把Z的行下標(biāo)當(dāng)作y坐標(biāo)軸的自變量。12/12/202263數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第63頁!5.3.3曲面/網(wǎng)線圖的精細(xì)修飾1.視角控制view改變觀察點(diǎn)是獲得較好三維視覺效果的重要途徑．一般而言，三維圖形的觀測角度是由方位角(Azimuth)及仰角(Elevation)來決定對二維圖形而言，默認(rèn)值為az=0°,el=90°;對三維圖形而言，默認(rèn)值為az=-37.5°,el=30°.ElevationAzimuth觀察點(diǎn)原點(diǎn)xzyview([az,el])>通過方位角，俯視角設(shè)置視點(diǎn).view([vx,vy,vz])>通過直角坐標(biāo)設(shè)置視點(diǎn).12/12/202264數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第64頁!指令功能shadingfaceted瓷磚”的顏色是均勻一致的，而且同時(shí)顯色“瓷磚”交接的“邊”。（此為缺省值）shadingflat“瓷磚”的顏色是均勻一致的。shadinginterp使用雙線性插值來使“瓷磚”根據(jù)四頂點(diǎn)的顏色產(chǎn)生連續(xù)的變化shadingoptions%options可取以下方式：3.濃淡處理12/12/202265數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第65頁!hiddenoff:用于網(wǎng)格圖形alpha(v)：0<v<1用于曲面圖形，0表示完全透明，1表示不透明。根據(jù)3個(gè)數(shù)據(jù)m×n矩陣X，Y和Z所繪制得到的曲面為例，Matlab有三種透明度的處理方式標(biāo)量：使所有數(shù)據(jù)點(diǎn)都設(shè)置相同的透明度線性數(shù)據(jù)：使曲面的數(shù)據(jù)點(diǎn)的透明度按照某個(gè)指定的維數(shù)的方向線性變化矩陣：使每個(gè)數(shù)據(jù)點(diǎn)選取不同的透明度4.透明控制12/12/202266數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第66頁!5.燈光設(shè)置lightcolor屬性表示光的顏色；可取相應(yīng)的顏色字符或RGB三元組。style屬性表示光源形式；選項(xiàng)有'infinite'和'local'兩個(gè)取值，分別表示無限遠(yuǎn)的平行光源和近光源。position屬性表示光源位置。position屬性的選項(xiàng)取三維坐標(biāo)點(diǎn)組成的向量形式[x，y，z]。默認(rèn)設(shè)置：白光，無窮遠(yuǎn)，穿過[1,0,1]射向坐標(biāo)原點(diǎn)。light('Color',option1,'Style',option2,'Position',option3)12/12/202267數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第67頁!7.控制光反射的材質(zhì)指令materialmaterialoptions其有如下一些選項(xiàng)：指令功能shiny使對象比較明亮，鏡面反射額較大，反射光的顏色取決于光源的顏色。dull使對象比較暗淡，漫反射額較大，反射光的顏色取決于光源的顏色。metal使對象有金屬光澤，反射光的顏色取決于光源的顏色和土星表面的顏色（此為缺省值）default返回到Matlab缺省值12/12/202268數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第68頁!5.3.4透視，鏤空和裁切1.圖形的透視matlab在采用缺省設(shè)置畫mesh圖形時(shí)，對疊壓在后面的圖形采取了消隱措施，若需要透視效果，可以通過相關(guān)的指令來實(shí)現(xiàn)．hiddenoff>透視被疊壓的圖形hiddenon>消隱被疊壓的圖形12/12/202269數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第69頁!2.圖形的鏤空例5.3-7演示如何利用非數(shù)NaN對圖形進(jìn)行鏤空處理P=peaks(30);surfc(P);%完整的圖形P(18:20,9:15)=NaN;%進(jìn)行鏤空處理surfc(P);%鏤空后的圖形colormap(hot);light('position',[50,-10,5]);lightingflatmaterial([0.9,0.9,0.6,15,0.4])12/12/202270數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第70頁!偽彩圖pcolor,等位線指令contour、contourf5.4高維可視化5.4.1二維半圖指令pcolor,contour,contourfcloseall;clear;[X,Y,Z]=peaks(40);subplot(1,2,1),pcolor(X,Y,Z)%shadinginterpsubplot(1,2,2),surf(X,Y,Z)%shadinginterp例5.4-1

二維半圖指令pcolor,contour,contourf應(yīng)用。12/12/202271數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第71頁!5.4高維可視化5.4.1二維半圖指令pcolor,contour,contourf偽彩圖pcolor；等位線指令contour、contourf；等位線標(biāo)高指令clabel的配合使用和區(qū)別。練習(xí)本例時(shí)注意：（1）本例等位線指令中的第4輸入宗量n設(shè)定高度的等級數(shù)，第5輸入宗量設(shè)定等位線的線型、色彩。（2）左右兩圖的標(biāo)高方法不同。左圖的標(biāo)識以“+”引導(dǎo)，水平放置。右圖沿線布置。這是由clabel的調(diào)用格式不同產(chǎn)生的。（3）左右兩圖色彩的形成方法不同，色彩效果也不同。（4）在左圖中，colorbar畫出一根垂直色標(biāo)尺，而caxis決定該色標(biāo)尺的刻度。例5.4-1

二維半圖指令pcolor,contour,contourf應(yīng)用。12/12/202272數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第72頁!clf;[x,y,z]=meshgrid(-2:.2:2,-2:.25:2,-2:.16:2);v=x.*exp(-x.^2-y.^2-z.^2);xs=[-0.7,0.7];ys=0;zs=0;slice(x,y,z,v,xs,ys,zs) colorbar;shadinginterpcolormaphsvxlabel('x'),ylabel('y'),zlabel('z')title('Thecolor-to-v(x,y,z)mapping')view([-22,39]);alpha(0.3)

[X,Y,Z]=meshgrid(x,y,z)slice(X,Y,Z,V,sx,sy,sz)2.切片圖例5.4-3

圖形表現(xiàn)。12/12/202273數(shù)據(jù)和函數(shù)的可視化共81頁，您現(xiàn)在瀏覽的是第73頁!2.色圖的變幻matlab為顏色的動(dòng)態(tài)變化提供了一個(gè)指令spinmap。它的功能是使當(dāng)前圖形的色圖做循環(huán)變化，以產(chǎn)生動(dòng)畫效果．該指令不涉及圖形對象特性的操作，只限于對色圖的操作。該指令只對256色設(shè)置有效。spinmap--->使色圖周期旋轉(zhuǎn)3sspinmap(t