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

lili456

7.3  三维图形
$ C, t* P  T- v% o' a三维图形表现的是三个变量之间的依赖关系，使用户能够直观看出二元函数数值的变化趋势。下面的小节将介绍三维图形的绘制函数。
; h6 z$ J, t6 y% }2 ^5 t: r7.3.1  plot3——绘制三维曲线
1 C/ G* a. G8 n% `: s【功能简介】绘制三维曲线。
【语法格式】
" L& n9 i1 I1 j% D4 X/ |% c1.plot3(X1,Y1,Z1,…)
# f. {: _$ V/ ~( i: sX1、Y1、Z1是同型的向量或矩阵，函数在三维空间中绘制出一条或多条曲线，绘制多条曲线时，曲线的条数等于矩阵的列数。
+ ]" P$ F9 u9 h9 b5 `; E格式变体：
        plot3(X1,Y1,Z1,LineSpec)：用参数LineSpec指定画线的线型、标记符号和颜色。
2.plot3(…,'PropertyName','PropertyValue')
5 ?. K* z" ]3 R设置所绘图形指定属性的属性值。
【实例7.27】绘制三维螺旋图。
>> t=0:pi/50:10*pi;
4 a4 {4 H7 a5 \- Z8 {. A2 T>> plot3(sin(t),cos(t),t);        %三维螺旋图
>> grid on;
>> axis square;
执行结果如图7-28所示。
0 _2 d% e5 p9 ^7 W2 f8 X2 j- m+ Z
图7-28  三维螺旋图
【实例分析】三维螺旋图的参数方程形式为x=sin(t),y=cos(t),z=t。
. \( V2 b# _8 b* \4 K: K" z

lili456

7.3.2  mesh——绘制三维网格图
$ H& E! q; I6 t! F【功能简介】绘制三维网格图。
【语法格式】
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定义， 颜色与高度成比例。
9 f3 D  C3 y1 h# G- F3 O2.mesh(Z)
' z+ ?/ e/ n! I* T+ ~' C  N9 y! AX与Y自动生成。[m,n]=size(Z)，则X=1:n，Y=1:m。颜色由高度决定。
, W, H2 ~4 R& Z6 j& F" F3.mesh(…,C)
5 _* Q( j$ p3 x' I图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
! H& f0 l9 G, |6 C6 A【实例7.28】绘制三维网格图。
+ }3 y' D8 `" m; b4 O>> [X,Y]=meshgrid(-8:.5:8);        %构造X、Y矩阵
>> R=sqrt(X.^2+Y.^2)+eps;
" r' V9 E8 F9 j( L7 F* i>> Z=sin(R)./R;
8 _6 D" o6 v- W2 q4 k( \>> mesh(X,Y,Z);                        %绘制三维网格图
3 b1 c* }0 U1 `0 Q执行结果如图7-29所示。
" v8 C- q: c; g+ R) E5 ]6 U
; s: H0 g# \% ?
图7-29  三维网格图
【实例分析】不指定颜色时，网格颜色由高度决定。
- i. p7 p; l/ T* s) n& J

lili456

7.3.3  surf——三维曲面图
【功能简介】绘制三维曲面图。
【语法格式】
1.surf(X,Y,Z)
& i4 u4 V. \) b+ ?* u生成由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函数一致。
! y% n7 e+ _8 P2 _# o* J2.surf (Z)
8 q7 c9 O/ j3 aX与Y自动生成。[m,n]=size(Z)，则X=1:n，Y=1:m。颜色由高度决定。
. s7 K( I1 o! e: t3.surf (…,C)
- T" p2 E' Y5 B2 T# j2 z3 f图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.29】绘制三维曲面图。
: R% q. _$ X5 L% ~  {>> [X,Y]=meshgrid(-8:.5:8);        %生成X、Y矩阵
7 U8 k) T: V3 |& I! f5 |>> R=sqrt(X.^2+Y.^2)+eps;
>> Z=sin(R)./R;
2 ?; z( a, V0 B, w6 T>> surf(X,Y,Z);                                %绘制三维曲面图
执行结果如图7-30所示。

" l4 P' e5 f4 G7 ~2 n' n图7-30  三维曲面图
8 a/ M. J* _/ w. A6 h( I7 X8 S! `【实例分析】surf函数与mesh函数调用格式是相同的。
2 r% Y- q4 G  b, w* K! Q& P

lili456

7.3.4  contour3——三维等高线绘制
' X- S; b0 G( ~# v  s* m4 W【功能简介】绘制三维等高线。
& M. G% F; g# ^# }5 f6 G6 A【语法格式】
1.contour3(Z)
) Y) e% p: a1 }7 g9 S# F画出矩阵Z的等高线图，矩阵Z至少为2×2大小，系统按矩阵Z中元素距离XY平面的高度，自动选择等高线的条数进行绘制。X和Y的值是自动确定的，[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
格式变体：
        contour3(Z,n)：画出矩阵Z的n条等高线。
* M4 B7 H! K: `7 b' T/ U3 m! s        contour3(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
) |! I9 w4 I+ {0 V; p1 j, L" ~2.contour3(X,Y,Z)
% Z# E$ r/ c6 E2 F3 P' m, J  |画出矩阵Z的等高线图，X轴和Y轴的范围由参数X、Y指定。如果X为矩阵，则使用X(1,来定义X轴的范围，如果Y为矩阵，则使用Y(:,1)来定义Y轴的范围。如果X与Y均为矩阵，则两者必须同型。
格式变体：
        contour3(X,Y,Z,n)：画出矩阵Z的n条等高线。
        contour3(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
3.contour3(…,LineSpec)
6 e, e5 }- J7 i  ?/ o9 Q8 s* {参数LineSpec指定线型、标记符号和颜色。
【实例7.30】绘制等高线。
>> [X,Y]=meshgrid([-2:.25:2]);
( [3 l* w# S! w7 i( k>> Z=X.*exp(-X.^2-Y.^2);
>> contour3(X,Y,Z,30);        %绘制等高线
>> colormap cool;
执行结果如图7-31所示。
) M4 f5 t# z) _3 ^* r- _
+ h! s( D- h) ?% b; y* F 图7-31  contour3函数绘制等高线
% u4 }/ I1 U# l) e, J0 l【实例分析】等高线绘制有contour3和contour两个函数均可完成，contour函数绘制的是二维的等高线。
$ n2 D, `7 Z" m' n0 a: l3 u

lili456

7.3.5  contour——曲面的等高线
【功能简介】绘制曲面的等高线。
( {2 B# j, B- }& E" H$ h9 R; J, f【语法格式】
1.contour(Z)
画出矩阵Z的二维等高线图，矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
2 P6 J9 S" R2 _0 }& J格式变体：
        contour(Z,n)：画出矩阵Z的n条等高线。
        contour(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
2.contour(X,Y,Z)
2 q2 g& `! d. ^  E/ P7 Z% _! h0 y画出矩阵Z的二维等高线图，X轴和Y轴的范围由参数X、Y指定。如果X与Y均为矩阵，则两者必须同型且单调递增。
格式变体：
        contour(X,Y,Z,n)：画出矩阵Z的n条等高线。
: x2 ^9 T$ ~" U' h3 R! ~" A        contour(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
3 G* t, c' P6 h/ {3.contour(…,LineSpec)
" t) R* u) y$ }5 D! b  m& V参数LineSpec指定线型、标记符号和颜色。
【实例7.31】画出曲面的等高线图。
7 @' h' |6 B. f2 x5 [/ F>> [X,Y]=meshgrid([-2:.25:2]);
& G. N# {. T; Z1 t$ ~0 P) B& p>> Z=X.*exp(-X.^2-Y.^2);
>> contour(X,Y,Z,30);        %画出曲面的二维等高线
执行结果如图7-32所示。
& @. j% b4 ?# i& i: Y: N+ R
图7-32  曲面的等高线
. V. @: B6 W% \2 i) ?3 g, v【实例分析】contour函数画出的是二维图形，contour3函数画出的是三维图形。

lili456

7.3.6  clabel——等高线高度标签
【功能简介】在二维等高线图中绘制等高线高度标签。
【语法格式】
# w9 [! t/ G7 }" E( U1.clabel(C,h)
C为标签矩阵，h为等高线句柄。函数旋转标签到适当角度，在空间允许的情况下插入标签。
格式变体：
3 ~1 I3 Q5 I$ k, ]6 C+ ^4 N  V        clabel(C,h,v)：在高度v处插入标签C。
        clabel(C,h,'manual')：手动添加标签到句柄h指定的等高线图中，用鼠标左键或空格键在最接近位置上放置标签，用回车键结束操作。
- ^% ^. G6 r# Q6 B0 _2.clabel(C)
" ?4 h* [# Y) Q: n' D添加标签到当前等高线图中，随机选择标签位置。
格式变体：
        clabel(C,v)：在v指定的高度处添加标签。
% ~8 E# x4 y( i2 n6 U1 Q        clabel(C,'manual')：手动添加标签到当前等高线图中。
# j7 H- {7 z( e) ?【实例7.32】给等高线做标注。
>> [x,y,z]=peaks;
>> C=contour(x,y,z);        %绘制等高线
>> clabel(C);                %给等高线添加标注
执行结果如图7-33所示。

8 t! S9 V4 D3 n7 O+ H+ ^图7-33  给等高线添加标注
. w6 O/ L0 S/ F【实例分析】peaks是MATLAB自带函数。
7.3.7  contourc——等高线图形计算
/ f, U2 i6 a% {" P9 ~3 F4 b9 R【功能简介】计算等高线矩阵C，用函数contour、contour3和contourf来显示。
【语法格式】
* B# a/ T( ^# j& c6 Y( W2 ^9 M1.C=contourc(Z)
$ g. m& F# t4 K/ P# j9 R3 ?0 d: c从矩阵Z中计算等高矩阵C，Z可视为XY平面的高度矩阵。等高线的数量和高度值是系统自动确定的。Z至少为2×2矩阵，至少包含两个不同的值，X和Y的范围是1:n与1:m，其中[m,n]=size(Z)。
格式变体：
        C=contour(Z,n)：确定等高线条数为n，返回等高矩阵。
0 @! n  U. ]- Y8 r' H        C=contour(Z,v)：在高度v处计算等高线。
2.C=contourc(X,Y,Z)
6 E+ J9 @: h: r3 }) q4 q在X和Y上计算Z的等高矩阵C，Z可视为XY平面的高度矩阵。等高线的数量和高度值是系统自动确定的。Z至少为2×2矩阵，至少包含两个不同的值。
2 @% ]4 ?" G, u9 q" R格式变体：
: R! F7 x! W1 N        C=contour(X,Y,Z,n)：确定等高线条数为n，返回等高矩阵。
, E1 `: s4 _, E. N/ o8 L7 Z        C=contour(X,Y,Z,v)：在高度v处计算等高线。
【实例7.33】计算peaks函数的等高矩阵。
>> a=peaks;
% |. B9 E3 j& x* V! M>> c=contourc(a,10);                %计算10条等高线的等高矩阵
>> s=size(c)                        %矩阵维度
) T! r; c0 W& a4 ^s =
     2   800
>> contour(a,c)                        %画出等高线
' y, b) o9 q. m. f3 O* i( z执行结果如图7-34所示。
1 N  t' b9 B5 }$ l! h$ C9 Z
  D3 E" n' i$ d+ h
图7-34  计算等高矩阵再绘制等高线
【实例分析】矩阵C是一个2×m矩阵，m的值取决于数据矩阵及等高线的条数。

lili456

7.3.8  fill3——填充三维图
. P' K, f& d, Z5 K' u4 W3 W& X7 @/ r【功能简介】填充三维多边形。
" S, h. T2 @+ {# X% r; d【语法格式】
8 z% t! \5 Q9 I; F0 B5 V& ^- \4 R! f1.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)。
格式变体：
, f0 q; v6 _- i8 V9 f% {        fill3(X1,Y1,Z1,C1,X2,Y2,Z2,C2,…)：绘制多个三维多边形。
2.fill3(X,Y,Z,ColorSpec)
用参数ColorSpec指定的颜色对三维多边形进行填充。
8 M# E( f& I  d9 d3 w- J【实例7.34】用渐变的颜色来对三维多边形进行着色。
>> X=[0,1,1,2;1,1,2,2;0,0,1,1];
>> 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];
>> C=[.5,1,1,.5;1,.5,.5,.1667;.3333,.3333,.5,.5]        %指定颜色
C =
' R# [: M$ D# s$ A3 N    0.5000    1.0000    1.0000    0.5000
. }) Z- [1 s, H% {4 r. r$ Q5 p    1.0000    0.5000    0.5000    0.1667
" G+ I  q( O* Q1 J    0.3333    0.3333    0.5000    0.5000
>> fill3(X,Y,Z,C)                                                %填充三维多边形
执行结果如图7-35所示。

& o; T" `& u$ D9 d6 b4 g& z图7-35  填充三维多边形
【实例分析】X、Y、Z决定三维多边形的形状，C确定颜色。
7.3.9  sphere——绘制球体
【功能简介】绘制球体。
【语法格式】
* F& i2 G6 x. y0 O) x1.sphere(n)
) e. S7 p- B. @$ x; n7 R- ?' K% I在三维直角坐标系中绘制由n×n个面组成的单位球体。
格式变体：
, P5 {$ v; `3 K) a7 |        sphere：默认n=20，绘制单位球体。
( a3 T$ m/ H4 [1 X5 i2.[X,Y,Z]=sphere(n)
返回三个阶数为(n+1)×(n+1)的坐标矩阵。该命令不画图，只是返回矩阵，要画图可以调用mesh(X,Y,Z)或surf(X,Y,Z)。
【实例7.35】绘制多个球体。
' D0 j3 `" F  O' [>> [x,y,z]=sphere;
>> surf(x,y,z);                %画出第一个球体
7 u2 B7 W1 t  {/ h! l: c>> hold on;
>> surf(x+3,y+2,z);        %画出第二个球体
+ i) r/ Z$ I2 K( x, Z>> surf(x,y-1,z+2);        %画出第三个球体
' T; d  q! {$ ]. X0 H( S>> daspect([1 1 1])
% O7 y) o. C* ?' q. r' _4 N执行结果如图7-36所示。
0 s  u) I9 F8 ^. A: ^4 G
5 g5 `$ c0 r! H+ E8 p, @8 g图7-36  绘制多个球体
【实例分析】sphere绘制单位球体，半径为1。

lili456

7.3.10  contourf——填充二维等高线
【功能简介】填充二维等高线。
【语法格式】
6 Q) |3 _  Q! ~1.contourf(Z)
- ^! H6 r$ N0 g& I$ ~; ]( i) A画出矩阵Z的二维等高线图，再对等高线之间的区域进行填充，填充的颜色使用当前窗口的颜色映射表。矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
$ Z1 N4 U5 w1 \- n$ |格式变体：
8 e* m3 P: v3 l1 A. W$ ^        contourf(Z,n)：画出矩阵Z的n条等高线并对等高线间的区域进行填充。
        contourf(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
2.contourf(X,Y,Z)
; K, P9 h; t" |0 I画出矩阵Z的等高线图并进行填充，X轴和Y轴的范围由参数X、Y指定。如果X与Y均为矩阵，则两者必须同型且单调递增。
格式变体：
( m3 T! k. ]( j# T5 R        contourf(X,Y,Z,n)：画出矩阵Z的n条等高线并进行填充。
- T4 d# \3 [8 L- e" k; I        contourf(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
【实例7.36】画出等高线并进行填充。
' J" b5 h. @# n3 O2 l" M>> a=peaks;
/ H- k! }. [. k/ u0 f/ s$ j>> contourf(a);                %画出用颜色填充的等高线
执行结果如图7-37所示。

图7-37  用颜色填充的等高线图
9 o. ~  k$ f( ?【实例分析】用颜色对等高线间的区域进行填充，使图形更直观。
7.3.11  pie3——三维饼图
【功能简介】绘制三维饼图。
【语法格式】
* }% [! l' Y  p, S" b! c% Y! Z6 @! `1.pie3(X)
使用X中的数据绘制一个三维饼图，X中的每一个元素都是饼图的一个部分。元素所占扇形面积的大小由X(i)/sum(X)决定，如果sum(X)<1，则只绘制不完整的三维饼形图。
2.pie3(X,explode)
1 q0 W) M2 w! T) T8 T- Q指定三维饼图中的每一部分是否分离出来。explode是与X同型的数组，用非零值表示相应位置的X中的元素分离出来显示。
# j6 f5 Y4 u2 H; f【实例7.37】绘制三维饼图。
$ ~+ @; y: p( I) B5 F* a* E>> x=[1,3,.7,2.5,2];
>> explode=[0,1,0,0,1];        %设置第二个元素和最后一个元素对应的区域分离出来独立显示
>> pie3(x,explode);                %绘制三维饼图
>> colormap hsv
3 y3 x* \$ _3 a$ D: B/ s执行结果如图7-38所示。

) b( u2 g5 ^: O+ x0 t7 w, E) a图7-38  三维饼图
3 E5 @! D# S6 f5 r4 K  r【实例分析】pie3的调用格式与pie类似。
, E( j! Y' a3 C, a) g. y# D: r5 i2 n

lili456

7.3.12  comet3——三维彗星图
【功能简介】绘制三维彗星图。
8 q% {) m0 V" v9 d4 l3 v- x3 G6 j【语法格式】
  K& a7 L& T- c- a+ H, N9 B7 H/ F& I1.comet3 (z)
6 j! ]3 V8 N. u. v5 [显示向量z的三维彗星图，X与Y轴范围由系统自动确定。
2.comet3(x,y,z)
* G2 s0 G5 I+ D1 E显示由x、y、z确定的三维曲线的彗星图。
3.comet3(x,y,z,p)
指定彗星体的长度为p*length(y)。
- [2 r1 q! i8 o: o9 V$ o& `5 k; [【实例7.38】绘制一个三维彗星图。
* \/ s" L7 }- S- Q8 ~>> t=-10*pi:pi/250:10*pi;
! h+ o2 }( ~, r8 q* K>> comet3((cos(2*t).^2).*sin(t),(sin(2*t).^2).*cos(t),t);                %绘制彗星图
执行结果如图7-39所示。

' t" N5 B# W& x7 j! ~图7-39  三维彗星图
【实例分析】彗星图的绘制是一个动态的过程，用户可按实例中的代码自行运行。
$ I4 K; {7 P( P6 m& c! G1 ?* H* I7.3.13  cylinder——生成圆柱图形
9 P2 }/ M% c2 S" o1 d. J+ Z3 Q【功能简介】绘制圆柱图形。
1.[X,Y,Z]=cylinder(r,n)
+ M8 H/ h: B0 ]$ g0 v4 l返回一个半径为r、高度为1的圆柱体的X、Y、Z轴坐标值，所绘制的圆柱体的圆周有指定的n个距离相同的点。用户可用surf或mesh函数和X、Y、Z参数画出圆柱。
! N* Q6 W; b, {, N5 E9 W2 P格式变体：
        [X,Y,Z]=cylinder(r)：所绘制的圆柱体的圆周有指定的20个距离相同的点。
        [X,Y,Z]= cylinder：半径采用默认值1。
2.cylinder(…)
/ R% S- _7 O& s# ]; u没有输出参量时，直接画出圆柱体。
【实例7.39】绘制一个两头粗、中间细的圆柱。
>> t=0:pi/10:2*pi;
>> [X,Y,Z]=cylinder(2+cos(t));                %用cylinder函数算出X、Y、Z参数
>> surf(X,Y,Z);                                        %画出圆柱体
>> axis square                                        %调整坐标轴
- K4 Y% |7 M# N2 f# f执行结果如图7-40所示。
# c# y7 ?* C8 C
2 ?4 i4 u" p8 v图7-40  绘制圆柱体
【实例分析】也可不带返回值，直接用cylinder(2+cos(t))画出圆柱。
" ^6 l" {( T- A9 Z7 F' ]/ V+ {

lili456

7.3.14  surfc——绘制阴影图及等高线
$ F5 Y9 B# e" L0 _+ [【功能简介】绘制阴影图及等高线。先画出三维曲面图，再在下方画出二维等高线。
【语法格式】
1 Z' Z. ?, N! A% K9 v1.surfc(X,Y,Z)
+ i0 M3 o3 \" c: U- `生成由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函数一致。
2.surfc(Z)
  T+ d0 K# V0 {, q% @3 lX与Y自动生成。[m,n]=size(Z)，则X=1:n，Y=1:m。颜色由高度决定。
0 ^6 H$ a* D7 h" ?* L- t2 |3.surfc(…,C)
2 W; N& |" U3 W/ ~图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
. A% e2 S* Q" p# ^1 v- U【实例7.40】绘制三维曲面及其二维等高线        。
# i0 E$ M$ |* ^$ o4 P  k1 o6 j>> [X,Y]=meshgrid(-8:.5:8);
>> R=sqrt(X.^2+Y.^2)+eps;
$ j+ H: S  v8 T' [! @' C( V1 m>> Z=sin(R)./R;
>> surfc (X,Y,Z);                        %绘制三维曲面及其二维等高线               
* v! h& J7 Y4 e. l8 {' h( [/ D执行结果如图7-41所示。
$ s  l( S/ s, ]; F( i
图7-41  绘制三维曲面及其二维等高线
【实例分析】suefc函数调用格式与surf相同，只是多画了一个二维等高线。
7.3.15  surfl——带光照模式的曲面图
【功能简介】绘制带光照模式的曲面图。
【语法格式】
1.surfl(X,Y,Z)
0 @' t9 \; C, Y4 V' c" j生成由X、Y和Z定义的曲面图，光照方向和系数采用默认值。
7 ^! G) i' E& A. w3 r* O) e) C( B+ [4 M, S格式变体：
0 x  I1 X5 z+ z# U* X. p4 E' [4 y        surfl (Z)：X、Y的值由系统自动生成，等于Z中元素的索引。
2.surfl(…,'light')
2 I0 _3 s, r% I用MATLAB光照对象生成一个带光照的曲面。
3.surfc(…,s)
0 Q0 P% P  s& f6 s+ S% n指定光源的方向s。s是长度为2或3的向量，即s=[azimuth,elevation]或[sx,sy,sz]。默认光源方向是从当前视角开始逆时针45度方向。
2 ]( K# L0 {& f% o6 y6 @4.surfc(…,s,k)
$ I* w4 u- Y7 t  r2 Q; y9 r$ _指定反射系数k，k是一个常量。k为长度为4的向量[ka,kd,ks,shine]，四个元素分别表示环境光系数、漫反射系数、镜面反射系数和镜面反射亮度，默认值为[0.55,0.6,0.4,10]。
1 {" f3 T9 b" S% e" [6 D【实例7.41】对peaks函数表示的曲面用surfl函数进行描绘。
; o% v  V2 s) }>> [x,y]=meshgrid(-3:1/8:3);
" I2 K0 }8 k. ]) K9 y>> z=peaks(x,y);                %peaks函数
" B- Z5 |4 D1 I& Z: o% b% z8 s  B>> surfl(x,y,z);                        %绘制带光照的曲面图
9 S5 M9 C3 p; J1 V>> shading interp;
>> colormap gray
- n$ O. W/ B# R- F7 n3 {执行结果如图7-42所示。
% |* g$ E6 g9 d# @/ |1 Q
! G- y! W$ w) B, |( y6 H2 }图7-42  绘制带光照的曲面图
9 I8 v; _1 Q9 s* O$ C【实例分析】命令colormap gray将色图指定为灰度色图。