【正文】 lab在采用缺省設(shè)置畫 mesh圖形時，對疊壓在后面的圖形采取了消隱措施，若需要透視效果，可以通過相關(guān)的指令來實(shí)現(xiàn)． hidden off透視被疊壓的圖形 hidden on消隱被疊壓的圖形 MATLAB 2022教程 第 4章 數(shù)值計(jì)算 2022年 5月 27日 第 58頁 例 5. 36 透視演示 [X0,Y0,Z0]=sphere(30)。 X=2*X0。Y=2*Y0。Z=2*Z0。 surf(X0,Y0,Z0)。 shading interp。hold on。 mesh(X,Y,Z), colormap(hot), hold off hidden off。 axis equal, axis off MATLAB 2022教程 第 4章 數(shù)值計(jì)算 2022年 5月 27日 第 59頁 2. 圖形的鏤空 例 演示如何利用非數(shù) NaN對圖形進(jìn)行鏤空處理 P=peaks(30)。 surfc(P)。 % 完整的圖形 P(18:20,9:15)=NaN。 % 進(jìn)行鏤空處理 surfc(P)。 % 鏤空后的圖形 colormap(hot )。 light( 39。position39。 , [50,10,5])。 lighting flat material([,15,]) MATLAB 2022教程 第 4章 數(shù)值計(jì)算 2022年 5月 27日 第 60頁 3. 裁切 由 NaN處理的圖形不會產(chǎn)生切面，為了看清圖形而需要表現(xiàn)切面，那么應(yīng)該把被切部分強(qiáng)制為零． 例 表現(xiàn)切面 x=[8::8]。 y=x。 [X,Y]=meshgrid(x,y)。 ZZ=X.^2+Y.^240。 ii=find(abs(X)6|abs(Y)6)。 %確定超出 [6,6]范圍的格點(diǎn)下標(biāo) ZZ(ii)=zeros(size(ii))。 % 強(qiáng)制為 0 surf(X,Y,ZZ),shading interp。colormap(copper) light(39。position39。,[0,15,1])。lighting phong material([,10,]) MATLAB 2022教程 第 4章 數(shù)值計(jì)算 2022年 5月 27日 第 61頁 偽彩圖 pcolor, 等位線指令 contour 、 contourf 高維可視化 二維半圖指令 pcolor, contour, contourf close all。 clear。 [X,Y,Z]=peaks(40)。 subplot(1,2,1),pcolor(X,Y,Z) % shading interp subplot(1,2,2), surf(X,Y,Z) % shading interp 例 二維半圖指令 pcolor, contour, contourf 應(yīng)用。 MATLAB 2022教程 第 4章 數(shù)值計(jì)算 2022年 5月 27日 第 62頁 close all。 clear。[X,Y,Z]=peaks(40)。n=6。 %等高線分級數(shù) subplot(1,2,1),pcolor(X,Y,Z) shading interp zmax=max(max(Z))。zmin=min(min(Z))。 caxis([zmin,zmax]) Colorbar hold on C=contour(X,Y,Z,n,39。k:39。)。 clabel(C) hold off subplot(1,2,2) [C,h]=contourf(X,Y,Z,n,39。k:39。)。 clabel(C,h) colormap(cool) set(gcf,39。Color39。,39。w39。) 例 二維半圖指令 pcolor, contour, contourf 應(yīng)用。 MATLAB 2022教程 第 4章 數(shù)值計(jì)算 2022年 5月 27日 第 63頁 高維可視化 二維半圖指令 pcolor, contour, contourf 偽彩圖 pcolor ；等位線指令 contour 、 contourf ；等位線標(biāo)高指令 clabel 的配合使用和區(qū)別。練習(xí)本例時注意：（ 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)尺的刻度。 例 二維半圖指令 pcolor, contour, contourf 應(yīng)用。 MATLAB 2022教程 第 4章 數(shù)值計(jì)算 2022年 5月 27日 第 64頁 四維表現(xiàn) clf x=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)) shading faceted。 colorbar(39。SouthOutside39。) brighten()。 colormap hsv title(39。No. 1 surf(X,Y,Z,abs(dzdr))39。) subplot(1,2,2)。surf(X,Y,Z,abs(dz2)) shading faceted colorbar(39。NorthOutside39。) title(39。No. 2 surf(X,Y,Z,abs(dz2))39。) 例 用顏色表現(xiàn) z=f(x,y) 函數(shù)的其他特征（如梯度，曲率等） 。 1. 準(zhǔn)四維表現(xiàn) MATLAB 2022教程 第 4章 數(shù)值計(jì)算 2022年 5月 27日 第 65頁 clf。[x,y,z] = meshgrid(2:.2:2,2:.25:2,2:.16:2)。 v =x.*exp(x.^2y.^2z.^2)。 xs = [,]。 ys = 0。 zs =0。 slice(x,y,z,v,xs,ys,zs) colorbar。 shading interp colormap hsv xlabel(39。x39。),ylabel(39。y39。),zlabel(39。z39。) title(39。The colortov(x,y,z) mapping39。) view([22,39])。 alpha() [X,Y,Z] = meshgrid(x,y,z) slice(X,Y,Z, V, sx, sy, sz) 2. 切片圖 例 圖形表現(xiàn) 。 2 2 2x y zv x e ? ? ??MATLAB 2022教程 第 4章 數(shù)值計(jì)算 2022年 5月 27日 第 66頁 動態(tài)圖形 1. 彗星狀軌跡圖 et(x,y,p) 二維彗星軌線 et(x,y,z,p) 三維彗星軌線 p 指定彗星體的長度 p*length(y)，缺省的 p值為 。三維圖形中彗星體的長度 p*length(z) shg。n=50。t=n*pi*(0::1)。x=sin(t)。y=cos(t)。 plot(x,y,39。g39。)。axis square hold on et(x,y,) hold off 例 簡單二維示例。 MATLAB 2022教程 第 4章 數(shù)值計(jì)算 2022年 5月 27日 第 67頁 2. 色圖的變幻 ?matlab為顏色的動態(tài)變化提供了一個指令 spinmap。它的功能是使當(dāng)前圖形的色圖做循環(huán)變化，以產(chǎn)生動畫效果．該指令不涉及圖形對象特性的操作，只限于對色圖的操作。該指令只對 256色設(shè)置有效。 ?spinmap使色圖周期旋轉(zhuǎn) 3s ?spinmap(t)使色圖周期旋轉(zhuǎn) ts ?spinmap(inf)使色圖 無限制旋轉(zhuǎn) 下去，用 Ctrl+C鍵終止旋轉(zhuǎn) ?spinmap(t,inc)分別用 t,inc(缺省為 2)控制色圖旋轉(zhuǎn)的時間和快慢 MATLAB 2022教程 第 4章 數(shù)值計(jì)算 2022年 5月 27日 第 68頁 2. 色圖的變幻 ?spinmap(t,inc)分別用 t,inc(缺省為 2)控制色圖旋轉(zhuǎn)的時間和快慢 ezsurf(39。x*y39。,39。circ39。)。shading flat。view([18,28]) C=summer。 CC=[C。 flipud( C )]。 colormap(CC) spinmap(30,4) 例 指令 spinmap的使用。 MATLAB 2022教程 第 4章 數(shù)值計(jì)算 2022年 5月 27日 第 69頁 3. 影片動畫 matlab支持影片動畫 (movie):先把一組二維或三維圖形存儲起來，然后再把這組圖形回放。由于“視覺殘留”，于是產(chǎn)生動畫效果。 M(i)=getframe對當(dāng)前圖形拍照后產(chǎn)生的數(shù)據(jù)向量依次存放于畫面結(jié)構(gòu)數(shù)組中。 movie(M,k)以不超過每秒 12幀的速度把 M中的畫面播放 k次。 ● 由 M(i)=getframe所產(chǎn)生的 M(i)是一個結(jié)構(gòu)數(shù)組，它有兩個域： M(i).cdata和 M(i).colormap ● 動畫的幾個典型的產(chǎn)生方法 : (1)改變某參數(shù) ,獲得一組畫面 ,如駐波 ,行波的產(chǎn)生 . (2)對產(chǎn)生的某三維圖形 ,改變觀察角 ,獲得一組畫面 . (3)對產(chǎn)生的某三維圖形 ,運(yùn)用 rotate旋轉(zhuǎn)指令 ,獲得一組畫面． MATLAB 2022教程 第 4章 數(shù)值計(jì)算 2022年 5月 27日 第 70頁 例 三維圖形的影片動畫。 clf,x=3*pi*(1::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)。axis off。 n=12。 mmm=moviein(n)。 for i=1:n rotate(h,[0 0 1],25)。 % 使圖形繞 Z軸 旋轉(zhuǎn) 25176。 /每次 mmm(:,i)=getframe。 % 捕獲畫面 end close shg。 axis off。 movie(mmm,20,10) % 以每秒 10幀 速度重復(fù)播放 20次 MATLAB 2022教程 第 4章 數(shù)值計(jì)算 2022年 5月 27日 第 71頁 4．實(shí)時動畫制作 ? 制作實(shí)時動畫的基本方法是，先畫出初始圖形，再計(jì)算活動對象的新位置，并在新位置上把它顯示出來，最后擦除原位置上原有的對象，刷新屏幕。重復(fù)操作即可產(chǎn)生動畫效果。 MATLAB 2022教程 第 4章 數(shù)值計(jì)算 2022年 5月 27日 第 72頁 例 模擬布朗運(yùn)動。程序如下：