【正文】 用 MATLAB提供的函數來定義色圖矩陣。 3．三維表面圖形的著色 三維表面圖實際上就是在網格圖的每一個網格片上涂上顏色。 surf函數用缺省的著色方式對網格片著色。除此之外，還可以用 shading命令來改變著色方式。 shading faceted命令將每個網格片用其高度對應的顏色進行著色，但網格線仍保留著，其顏色是黑色。這是系統(tǒng)的缺省著色方式。 shading flat命令將每個網格片用同一個顏色進行著色，且網格線也用相應的顏色，從而使得圖形表面顯得更加光滑。 shading interp命令在網格片內采用顏色插值處理，得出的表面圖顯得最光滑。 例 223 3種圖形著色方式的效果展示。 程序如下： [x,y,z]=sphere(20)。 colormap(copper)。 subplot(1,3,1)。 surf(x,y,z)。 axis equal subplot(1,3,2)。 surf(x,y,z)。shading flat。 axis equal subplot(1,3,3)。 surf(x,y,z)。shading interp。 axis equal 101101101101101101101101101 圖形的裁剪處理 例 225 繪制三維曲面圖，并進行插值著色處理，裁掉圖中 x和 y都小于 0部分。 程序如下： [x,y]=meshgrid(5::5)。 z=cos(x).*cos(y).*exp(sqrt(x.^2+y.^2)/4)。 surf(x,y,z)。shading interp。 pause %程序暫停 i=find(x=0amp。y=0)。 z1=z。z1(i)=NaN。 surf(x,y,z1)。shading interp。 為了展示裁剪效果，第一個曲面繪制完成后暫停，然后顯示裁剪后的曲面。 505505 0 . 500 . 51 動畫制作 MATLAB提供 getframe、 moviein和 movie函數進行動畫制作。 1． getframe函數 getframe函數可截取一幅畫面信息 (稱為動畫中的一幀 )，一幅畫面信息形成一個很大的列向量。顯然，保存 n幅圖面就需一個大矩陣。 2． moviein函數 moviein(n)函數用來建立一個足夠大的 n列矩陣。該矩陣用來保存 n幅畫面的數據，以備播放。之所以要事先建立一個大矩陣，是為了提高程序運行速度。 3． movie函數 movie(m,n)函數播放由矩陣 m所定義的畫面 n次，缺省時播放一次。 例 227 繪制了 peaks函數曲面并且將它繞 z軸旋轉。 程序如下 [X,Y,Z]=peaks(30)。 surf(X,Y,Z) axis([3,3,3,3,10,10]) axis off。 shading interp。 colormap(hot)。 m=moviein(20)。 %建立一個 20列大矩陣 for i=1:20 view(+24*(i1),30) %改變視點 m(:,i)=getframe。 %將圖形保存到 m矩陣 end movie(m,2)。 %播放畫面 2次 應用舉例 clear。 t=1:9。 d1=[ ]。 d2=[ ]。 d3=[ ]。 plot(t,d1,39。39。,t,d2,39。gx:39。,t,d3,39。m*.39。)。 title(39。穩(wěn)定性變化規(guī)律 39。)。 xlabel(39。時間 39。)。 ylabel(39。穩(wěn)定性 39。)。 axis([0 10 0 100])。 text(,39。\leftarrow樣品 a39。)。 text(3,39。樣品 b\rightarrow39。)。 text(,39。\leftarrow樣品 c39。)。 0 1 2 3 4 5 6 7 8 9 100102030405060708090100穩(wěn)定性變化規(guī)律時間穩(wěn)定性? 樣品 a樣品 b ?? 樣品 c作業(yè) : 1在同一坐標內，分別用不同線型和顏色繪制曲線 y1=(4πx) 和 y2=(πx)，標記兩曲 線交叉點。 2 繪制三維曲面圖 z=sin(x+sin(y))x/10。 ? x=linspace(0,2*pi,1000)。 ? y1=*exp(*x).*cos(4*pi*x)。 ? y2=2*exp(*x).*cos(pi*x)。 ? k=find(abs(y1y2)1e2)。 % ? x1=x(k)。 % ? y3=*exp(*x1).*cos(4*pi*x1)。 % ? plot(x,y1,x,y2,39。k:39。,x1,y3,39。bp39。)。 繪制山區(qū)地貌圖 要在某山區(qū)方圓大約 27km2范圍內修建一條公路，從山腳出發(fā)經過一個居民區(qū)，再到達一個礦區(qū) .橫向縱向分別每隔 400m測量一次，得到一些地點的高程： (平面區(qū)域0≤x≤ 5600,0≤y≤4800)，需作出 該山區(qū)的地貌圖和等高線圖 . 返回 要求 x0,y0單調； x， y可取 為矩陣，或 x取行向量， y取為列向量， x,y的值分別不能超出x0,y0的范圍。 z=interp2(x0,y0,z0,x,y,’method’) 被插值點 插值方法 用 MATLAB作網格節(jié)點數據的插值 插值節(jié)點 被插值點的函數值 ‘nearest’ 最鄰近插值 ‘ linear’ 雙線性插值 ‘ cubic’ 雙三次插值 缺省時 , 雙線性插值