脑动力：MATLAB函数功能速查效率手册

lili456

7.3  三维图形
: t& u+ }( m$ \) U( c, R: r- ~+ ~三维图形表现的是三个变量之间的依赖关系，使用户能够直观看出二元函数数值的变化趋势。下面的小节将介绍三维图形的绘制函数。
7.3.1  plot3——绘制三维曲线
) G; r, t: i: p【功能简介】绘制三维曲线。
$ y! v# L9 W+ ^4 W& |+ F【语法格式】
1.plot3(X1,Y1,Z1,…)
X1、Y1、Z1是同型的向量或矩阵，函数在三维空间中绘制出一条或多条曲线，绘制多条曲线时，曲线的条数等于矩阵的列数。
$ X1 Q& r; ^+ M# F% N格式变体：
        plot3(X1,Y1,Z1,LineSpec)：用参数LineSpec指定画线的线型、标记符号和颜色。
6 Y' \: ~& `) F% H4 U% N: _/ R; x2.plot3(…,'PropertyName','PropertyValue')
设置所绘图形指定属性的属性值。
1 N0 S( P3 A% q) {7 j5 S【实例7.27】绘制三维螺旋图。
>> t=0:pi/50:10*pi;
, ]% r- v0 [$ a4 Y  F* t; u>> plot3(sin(t),cos(t),t);        %三维螺旋图
3 ~/ G; ]4 y' K, `>> grid on;
* M5 ~+ k! Z+ F( f+ M; q' v1 V2 e7 E7 p>> axis square;
执行结果如图7-28所示。
' l0 R$ v3 E: Z2 x  s
图7-28  三维螺旋图
【实例分析】三维螺旋图的参数方程形式为x=sin(t),y=cos(t),z=t。

lili456

7.3.2  mesh——绘制三维网格图
: v' K) Y$ U9 v! b* n0 W1 l, B【功能简介】绘制三维网格图。
【语法格式】
1.mesh(X,Y,Z)
生成由X、Y和Z定义的网格图。X和Y如果分别是长度为m、n的向量，且(n,m)=size(Z)，则生成的网格线的交叉点为[X(j),Y(i),Z(i,j)]。如果X、Y分别为矩阵，则生成网格线的交叉点为[X(i,j),Y(i,j),Z(i,j)]。网格线的颜色由Z定义， 颜色与高度成比例。
2.mesh(Z)
X与Y自动生成。[m,n]=size(Z)，则X=1:n，Y=1:m。颜色由高度决定。
3.mesh(…,C)
' L$ W3 H5 }: {( a& r" u/ ]. r图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
6 r) w9 E5 O7 V- Z3 y【实例7.28】绘制三维网格图。
>> [X,Y]=meshgrid(-8:.5:8);        %构造X、Y矩阵
, G2 Q2 ^" i" ]>> R=sqrt(X.^2+Y.^2)+eps;
5 Y/ i+ ?- Q1 C& J1 L8 b% Q1 x  D>> Z=sin(R)./R;
>> mesh(X,Y,Z);                        %绘制三维网格图
执行结果如图7-29所示。


3 \" \2 y) n6 c: v( A图7-29  三维网格图
【实例分析】不指定颜色时，网格颜色由高度决定。

lili456

7.3.3  surf——三维曲面图
【功能简介】绘制三维曲面图。
【语法格式】
/ k7 g9 @9 Q/ L( m; C5 @; X+ u1.surf(X,Y,Z)
生成由X、Y和Z定义的曲面图。如果X和Y分别是长度为m、n的向量，且(n,m)=size(Z)，则生成的曲面中的交叉点为[X(j),Y(i),Z(i,j)]。如果X、Y分别为矩阵，则生成曲面中的交叉点为[X(i,j),Y(i,j),Z(i,j)]。线条之间的区域用颜色填充，surf函数调用格式与mesh函数一致。
; T; b# x' E. Y& K' u2.surf (Z)
6 c+ q+ @/ v! t5 lX与Y自动生成。[m,n]=size(Z)，则X=1:n，Y=1:m。颜色由高度决定。
# f- y. |$ u+ ]9 c3.surf (…,C)
3 }, {1 P& Q5 b0 L) _/ T8 p图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.29】绘制三维曲面图。
>> [X,Y]=meshgrid(-8:.5:8);        %生成X、Y矩阵
>> R=sqrt(X.^2+Y.^2)+eps;
) n" ~9 l8 L8 B+ n6 {9 G4 F3 H3 V>> Z=sin(R)./R;
>> surf(X,Y,Z);                                %绘制三维曲面图
% k; j4 ?+ P1 ~" X9 Y2 j4 E执行结果如图7-30所示。
3 d  c4 d' b! G# w5 N0 m1 x0 q
图7-30  三维曲面图
8 l* E9 w+ H" \0 T+ z, x& J【实例分析】surf函数与mesh函数调用格式是相同的。

lili456

7.3.4  contour3——三维等高线绘制
! m+ w) ^  Z7 E【功能简介】绘制三维等高线。
2 F. p2 O* a3 d【语法格式】
1.contour3(Z)
画出矩阵Z的等高线图，矩阵Z至少为2×2大小，系统按矩阵Z中元素距离XY平面的高度，自动选择等高线的条数进行绘制。X和Y的值是自动确定的，[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
3 T3 E, g, w( c+ W' x- i' }6 c格式变体：
        contour3(Z,n)：画出矩阵Z的n条等高线。
        contour3(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
2.contour3(X,Y,Z)
画出矩阵Z的等高线图，X轴和Y轴的范围由参数X、Y指定。如果X为矩阵，则使用X(1,来定义X轴的范围，如果Y为矩阵，则使用Y(:,1)来定义Y轴的范围。如果X与Y均为矩阵，则两者必须同型。
格式变体：
" X" Z9 C! P/ g        contour3(X,Y,Z,n)：画出矩阵Z的n条等高线。
0 f* a' y/ V5 k3 y' S        contour3(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
3.contour3(…,LineSpec)
2 {. W+ F0 C* O# b8 h参数LineSpec指定线型、标记符号和颜色。
【实例7.30】绘制等高线。
>> [X,Y]=meshgrid([-2:.25:2]);
# e9 X4 B( U# y* ?( P/ K' o>> Z=X.*exp(-X.^2-Y.^2);
# g: }. _2 Q. T8 d- b( i& p5 X) h: a>> contour3(X,Y,Z,30);        %绘制等高线
+ A" t5 v9 p" ~: ~( y' Z4 |, u>> colormap cool;
执行结果如图7-31所示。
& c/ h' j! ~& d' t- [0 V% O
图7-31  contour3函数绘制等高线
【实例分析】等高线绘制有contour3和contour两个函数均可完成，contour函数绘制的是二维的等高线。

lili456

7.3.5  contour——曲面的等高线
【功能简介】绘制曲面的等高线。
【语法格式】
1.contour(Z)
/ i/ f( R' h6 B  U8 h) M画出矩阵Z的二维等高线图，矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
8 `- \- I* @1 m4 v格式变体：
* x8 [, C3 f  j# e' ^        contour(Z,n)：画出矩阵Z的n条等高线。
        contour(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
2.contour(X,Y,Z)
画出矩阵Z的二维等高线图，X轴和Y轴的范围由参数X、Y指定。如果X与Y均为矩阵，则两者必须同型且单调递增。
- ]- L! m6 R% l, G) O) K格式变体：
$ R4 v6 M. E; t3 O) i2 C( H        contour(X,Y,Z,n)：画出矩阵Z的n条等高线。
- `2 y$ B  y0 A$ y* x        contour(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
3.contour(…,LineSpec)
参数LineSpec指定线型、标记符号和颜色。
- L8 P9 \* a, O2 w. D【实例7.31】画出曲面的等高线图。
8 ?3 y+ Q- D' ]; v7 w9 H$ l  E% o>> [X,Y]=meshgrid([-2:.25:2]);
7 \! m$ O* f5 i; Y>> Z=X.*exp(-X.^2-Y.^2);
>> contour(X,Y,Z,30);        %画出曲面的二维等高线
2 S7 D9 \5 t& Q; H0 ]0 V5 I3 b执行结果如图7-32所示。

图7-32  曲面的等高线
【实例分析】contour函数画出的是二维图形，contour3函数画出的是三维图形。
+ S' V1 v9 U8 b! U8 U' X, K) J

lili456

7.3.6  clabel——等高线高度标签
4 ?- F( L. ]6 |1 M3 A8 [【功能简介】在二维等高线图中绘制等高线高度标签。
【语法格式】
1.clabel(C,h)
C为标签矩阵，h为等高线句柄。函数旋转标签到适当角度，在空间允许的情况下插入标签。
格式变体：
        clabel(C,h,v)：在高度v处插入标签C。
- x& {* K2 d6 ~" J! h# x        clabel(C,h,'manual')：手动添加标签到句柄h指定的等高线图中，用鼠标左键或空格键在最接近位置上放置标签，用回车键结束操作。
1 N5 k  T  x, Z# b/ X+ N3 ?2.clabel(C)
添加标签到当前等高线图中，随机选择标签位置。
格式变体：
        clabel(C,v)：在v指定的高度处添加标签。
( ?+ e# S' b7 P* U2 R        clabel(C,'manual')：手动添加标签到当前等高线图中。
【实例7.32】给等高线做标注。
>> [x,y,z]=peaks;
>> C=contour(x,y,z);        %绘制等高线
" L; V$ l2 u* J, Y>> clabel(C);                %给等高线添加标注
) f$ Z: Z+ R4 R3 W执行结果如图7-33所示。

图7-33  给等高线添加标注
【实例分析】peaks是MATLAB自带函数。
. Q; E6 w# T% E  W% J; ]2 ?7.3.7  contourc——等高线图形计算
【功能简介】计算等高线矩阵C，用函数contour、contour3和contourf来显示。
! V2 d) a! P, J/ a! O& p【语法格式】
1.C=contourc(Z)
0 h. D$ {2 h8 R3 ?从矩阵Z中计算等高矩阵C，Z可视为XY平面的高度矩阵。等高线的数量和高度值是系统自动确定的。Z至少为2×2矩阵，至少包含两个不同的值，X和Y的范围是1:n与1:m，其中[m,n]=size(Z)。
# |* ]/ B: n! O' n3 R8 {格式变体：
        C=contour(Z,n)：确定等高线条数为n，返回等高矩阵。
        C=contour(Z,v)：在高度v处计算等高线。
) \5 m2 ?: _! H+ w; Q2.C=contourc(X,Y,Z)
在X和Y上计算Z的等高矩阵C，Z可视为XY平面的高度矩阵。等高线的数量和高度值是系统自动确定的。Z至少为2×2矩阵，至少包含两个不同的值。
* t7 }$ o2 I3 P6 l  q0 F, p格式变体：
* o* ]' j3 f" W+ Q. U* c  \2 O        C=contour(X,Y,Z,n)：确定等高线条数为n，返回等高矩阵。
        C=contour(X,Y,Z,v)：在高度v处计算等高线。
3 d# Q, M* I7 ~8 L【实例7.33】计算peaks函数的等高矩阵。
7 N( a# O' ]9 B* |3 U  G0 u4 C0 X>> a=peaks;
( G  z5 S$ C7 \, m) D7 L! D6 k/ j>> c=contourc(a,10);                %计算10条等高线的等高矩阵
>> s=size(c)                        %矩阵维度
s =
     2   800
>> contour(a,c)                        %画出等高线
执行结果如图7-34所示。
7 C4 \. l0 L8 C# o$ e
# T! x/ y0 N# l% p4 [1 q% z
图7-34  计算等高矩阵再绘制等高线
* O1 k. \: L9 i2 J% ?【实例分析】矩阵C是一个2×m矩阵，m的值取决于数据矩阵及等高线的条数。

lili456

7.3.8  fill3——填充三维图
【功能简介】填充三维多边形。
【语法格式】
1.fill3(X,Y,Z,C)
参数X、Y和Z定义多边形的顶点，如果X、Y和Z是矩阵，fill3将创建n个多边形，n为矩阵的列数。函数将在必要时闭合最后一个顶点与第一个顶点来封闭多边形。参数C指定颜色，如果C是一个行向量，则length(C)等于size(X,2)和size(Y,2)，如果C是一个列向量，则length(C)必须等于size(X,1)和size(Y,1)。
( H4 b; I' K7 [6 j格式变体：
        fill3(X1,Y1,Z1,C1,X2,Y2,Z2,C2,…)：绘制多个三维多边形。
2.fill3(X,Y,Z,ColorSpec)
用参数ColorSpec指定的颜色对三维多边形进行填充。
【实例7.34】用渐变的颜色来对三维多边形进行着色。
>> X=[0,1,1,2;1,1,2,2;0,0,1,1];
4 Y8 f) B2 J" @) Y>> Y=[1,1,1,1;1,0,1,0;0,0,0,0];
>> Z=[1,1,1,1;1,0,1,0;0,0,0,0];
& Q5 [- W4 a% X" j4 i>> C=[.5,1,1,.5;1,.5,.5,.1667;.3333,.3333,.5,.5]        %指定颜色
& `1 E3 l% l* v6 dC =
    0.5000    1.0000    1.0000    0.5000
9 j8 {8 H- o6 }; x2 B    1.0000    0.5000    0.5000    0.1667
    0.3333    0.3333    0.5000    0.5000
4 P; u! R+ Y- x5 k; i+ ?>> fill3(X,Y,Z,C)                                                %填充三维多边形
4 M& |" |# ~) Y. w执行结果如图7-35所示。

5 }5 O% v: C4 n; J图7-35  填充三维多边形
【实例分析】X、Y、Z决定三维多边形的形状，C确定颜色。
7.3.9  sphere——绘制球体
【功能简介】绘制球体。
【语法格式】
1.sphere(n)
在三维直角坐标系中绘制由n×n个面组成的单位球体。
; |+ c( O7 @4 E格式变体：
        sphere：默认n=20，绘制单位球体。
2.[X,Y,Z]=sphere(n)
% f! h# M  A) D2 g$ m返回三个阶数为(n+1)×(n+1)的坐标矩阵。该命令不画图，只是返回矩阵，要画图可以调用mesh(X,Y,Z)或surf(X,Y,Z)。
; Z& ]7 C$ Q# w9 p* X【实例7.35】绘制多个球体。
>> [x,y,z]=sphere;
>> surf(x,y,z);                %画出第一个球体
>> hold on;
>> surf(x+3,y+2,z);        %画出第二个球体
- V/ t" y, W6 S7 y# b- y6 w>> surf(x,y-1,z+2);        %画出第三个球体
>> daspect([1 1 1])
" Q/ h1 b3 |) V& T# y+ {执行结果如图7-36所示。
8 V; B/ o4 y4 _+ h( {9 j
) I% }+ K5 ^. _图7-36  绘制多个球体
1 N& e5 t* ]( R* z2 f! W2 g! H【实例分析】sphere绘制单位球体，半径为1。
6 `" ]# F  o, l. x4 Y) s( e

lili456

7.3.10  contourf——填充二维等高线
- M  A9 l" `. [0 D7 v. t【功能简介】填充二维等高线。
【语法格式】
1.contourf(Z)
画出矩阵Z的二维等高线图，再对等高线之间的区域进行填充，填充的颜色使用当前窗口的颜色映射表。矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
格式变体：
* Q; x: ?# M: N2 s' f" o+ ~& c: a        contourf(Z,n)：画出矩阵Z的n条等高线并对等高线间的区域进行填充。
        contourf(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
% i+ Z7 A9 C1 p2.contourf(X,Y,Z)
- L3 O7 G$ K! r; x% p画出矩阵Z的等高线图并进行填充，X轴和Y轴的范围由参数X、Y指定。如果X与Y均为矩阵，则两者必须同型且单调递增。
格式变体：
& O$ L  `" V0 @" }2 u& ]        contourf(X,Y,Z,n)：画出矩阵Z的n条等高线并进行填充。
        contourf(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
【实例7.36】画出等高线并进行填充。
: h3 W) B+ F0 _5 k>> a=peaks;
>> contourf(a);                %画出用颜色填充的等高线
& R6 t+ r6 t4 p; T) g/ ?; i, u. Y! c执行结果如图7-37所示。
- l# D% l9 ?# x0 }$ K3 G
- O; d/ `  v( _" _  P) B& b4 B4 Z图7-37  用颜色填充的等高线图
【实例分析】用颜色对等高线间的区域进行填充，使图形更直观。
9 m. l# J5 a( H8 z7 a) U1 M7.3.11  pie3——三维饼图
) r6 R+ u. b) L7 S* g  y. T【功能简介】绘制三维饼图。
6 ]& X) o6 c& _. z; H1 ~【语法格式】
+ x+ z4 S- n: m2 ]! l, Y1.pie3(X)
使用X中的数据绘制一个三维饼图，X中的每一个元素都是饼图的一个部分。元素所占扇形面积的大小由X(i)/sum(X)决定，如果sum(X)<1，则只绘制不完整的三维饼形图。
% A" N9 |( p$ X# {3 h- Y/ V3 ^2.pie3(X,explode)
指定三维饼图中的每一部分是否分离出来。explode是与X同型的数组，用非零值表示相应位置的X中的元素分离出来显示。
【实例7.37】绘制三维饼图。
>> x=[1,3,.7,2.5,2];
7 a8 u* }1 P( z/ ~: q>> explode=[0,1,0,0,1];        %设置第二个元素和最后一个元素对应的区域分离出来独立显示
4 h/ D% F! Z& F6 t, T' q1 p>> pie3(x,explode);                %绘制三维饼图
>> colormap hsv
执行结果如图7-38所示。

1 p" W0 a/ e& I+ D; A# J2 P8 u! \图7-38  三维饼图
% d- g; C4 O( g1 b【实例分析】pie3的调用格式与pie类似。
: s5 ~7 W# J% Q! e! u3 g" T+ X

lili456

7.3.12  comet3——三维彗星图
【功能简介】绘制三维彗星图。
【语法格式】
( S0 U- Z2 e4 w1.comet3 (z)
显示向量z的三维彗星图，X与Y轴范围由系统自动确定。
0 ^: F7 _6 R- {- ~% Z* n2.comet3(x,y,z)
5 ~) j# i8 d4 V( n$ }- a% q6 ]! _/ B显示由x、y、z确定的三维曲线的彗星图。
3.comet3(x,y,z,p)
& P" ~  e$ g1 H2 @* y- |) O" l指定彗星体的长度为p*length(y)。
【实例7.38】绘制一个三维彗星图。
1 v# p  s& b6 V/ U# D, p5 o9 O>> t=-10*pi:pi/250:10*pi;
>> comet3((cos(2*t).^2).*sin(t),(sin(2*t).^2).*cos(t),t);                %绘制彗星图
8 E- F9 D; S5 \' j3 y执行结果如图7-39所示。

图7-39  三维彗星图
【实例分析】彗星图的绘制是一个动态的过程，用户可按实例中的代码自行运行。
7.3.13  cylinder——生成圆柱图形
【功能简介】绘制圆柱图形。
1.[X,Y,Z]=cylinder(r,n)
返回一个半径为r、高度为1的圆柱体的X、Y、Z轴坐标值，所绘制的圆柱体的圆周有指定的n个距离相同的点。用户可用surf或mesh函数和X、Y、Z参数画出圆柱。
/ z! b8 F& n) v/ L" n7 E格式变体：
        [X,Y,Z]=cylinder(r)：所绘制的圆柱体的圆周有指定的20个距离相同的点。
5 F6 ~$ D" z' Y) ~  d2 I        [X,Y,Z]= cylinder：半径采用默认值1。
9 \* G! Q0 {8 v5 z$ O0 T* r) R; W2.cylinder(…)
没有输出参量时，直接画出圆柱体。
【实例7.39】绘制一个两头粗、中间细的圆柱。
% k5 h/ W& v/ W>> t=0:pi/10:2*pi;
>> [X,Y,Z]=cylinder(2+cos(t));                %用cylinder函数算出X、Y、Z参数
9 n$ }/ s% H9 B+ }; f& z>> surf(X,Y,Z);                                        %画出圆柱体
>> axis square                                        %调整坐标轴
执行结果如图7-40所示。

图7-40  绘制圆柱体
【实例分析】也可不带返回值，直接用cylinder(2+cos(t))画出圆柱。
7 s' b' r& q& U+ G, |9 @

lili456

7.3.14  surfc——绘制阴影图及等高线
【功能简介】绘制阴影图及等高线。先画出三维曲面图，再在下方画出二维等高线。
5 p3 ]! o. L+ p3 T8 J; r, i) k【语法格式】
+ _; J, r8 W& M: X1.surfc(X,Y,Z)
生成由X、Y和Z定义的曲面图，并在下方绘制二维等高线。X和Y如果分别是长度为m、n的向量，且[n,m]=size(Z)，则生成的曲面中的交叉点为[X(j),Y(i),Z(i,j)]。如果X、Y分别为矩阵，则生成曲面中的交叉点为[X(i,j),Y(i,j),Z(i,j)]。线条之间的区域用颜色填充，surfc函数调用格式与surf函数一致。
' I% z) K2 g! |2 i* k4 \# s2.surfc(Z)
( A. ^6 n1 H! Q( c! pX与Y自动生成。[m,n]=size(Z)，则X=1:n，Y=1:m。颜色由高度决定。
6 j4 H5 x: n$ W) v5 v9 F3.surfc(…,C)
& k0 U$ |; R9 g! q  U图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.40】绘制三维曲面及其二维等高线        。
  Y& z" u: F8 B  c>> [X,Y]=meshgrid(-8:.5:8);
, K! g8 b( i. g  J7 w  ]7 |) `>> R=sqrt(X.^2+Y.^2)+eps;
>> Z=sin(R)./R;
1 m' a! P' F/ [& A9 s>> surfc (X,Y,Z);                        %绘制三维曲面及其二维等高线               
2 o6 a6 B$ [; w3 ]3 ^6 J执行结果如图7-41所示。
8 V0 `- A5 F; T# g
3 X4 h) U2 y6 _5 i, ]" H* @' T图7-41  绘制三维曲面及其二维等高线
# P# Z8 ~4 G; K; ?' F7 |【实例分析】suefc函数调用格式与surf相同，只是多画了一个二维等高线。
: |( j" W3 o6 |; `7.3.15  surfl——带光照模式的曲面图
$ Y: u" v7 q: p) ~. g' T【功能简介】绘制带光照模式的曲面图。
【语法格式】
1.surfl(X,Y,Z)
生成由X、Y和Z定义的曲面图，光照方向和系数采用默认值。
6 A4 Q5 j0 `6 L1 @格式变体：
        surfl (Z)：X、Y的值由系统自动生成，等于Z中元素的索引。
2.surfl(…,'light')
+ k1 K$ Y8 n0 U; F2 Q9 D1 k7 ?# y用MATLAB光照对象生成一个带光照的曲面。
1 ^* b" E, H/ z5 r, t1 L$ v% A3 H3.surfc(…,s)
指定光源的方向s。s是长度为2或3的向量，即s=[azimuth,elevation]或[sx,sy,sz]。默认光源方向是从当前视角开始逆时针45度方向。
( T' O0 ~- E% O; K! K) T( a7 K7 c4.surfc(…,s,k)
2 P$ S1 k% k% _4 U) Q* h指定反射系数k，k是一个常量。k为长度为4的向量[ka,kd,ks,shine]，四个元素分别表示环境光系数、漫反射系数、镜面反射系数和镜面反射亮度，默认值为[0.55,0.6,0.4,10]。
9 K) v& F8 Q8 Q8 m1 T. d3 U【实例7.41】对peaks函数表示的曲面用surfl函数进行描绘。
>> [x,y]=meshgrid(-3:1/8:3);
>> z=peaks(x,y);                %peaks函数
/ M, P) Z; J* n>> surfl(x,y,z);                        %绘制带光照的曲面图
>> shading interp;
>> colormap gray
$ N$ C3 ^0 A- G* s& b: H执行结果如图7-42所示。

' h" @- H! P( B4 {6 G  T9 D$ d图7-42  绘制带光照的曲面图
) S) e5 m7 A7 d- e2 p8 V" W【实例分析】命令colormap gray将色图指定为灰度色图。
5 B$ D8 b- r0 @: N. I