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

lili456

7.3  三维图形
9 w% {, X# U% O& Q三维图形表现的是三个变量之间的依赖关系，使用户能够直观看出二元函数数值的变化趋势。下面的小节将介绍三维图形的绘制函数。
8 A% [6 X9 R) \# U9 T& G& a3 g7.3.1  plot3——绘制三维曲线
4 C  q* X( O  H4 N【功能简介】绘制三维曲线。
【语法格式】
1.plot3(X1,Y1,Z1,…)
X1、Y1、Z1是同型的向量或矩阵，函数在三维空间中绘制出一条或多条曲线，绘制多条曲线时，曲线的条数等于矩阵的列数。
2 b6 e3 ]4 L% s格式变体：
        plot3(X1,Y1,Z1,LineSpec)：用参数LineSpec指定画线的线型、标记符号和颜色。
* Y; S0 D* A( ]8 L8 ^2.plot3(…,'PropertyName','PropertyValue')
6 j' M: q; E# p- O2 n9 b6 J设置所绘图形指定属性的属性值。
" Y% U. V  G& X: r% R& {$ e【实例7.27】绘制三维螺旋图。
>> t=0:pi/50:10*pi;
" R' `: _: z* P8 A>> plot3(sin(t),cos(t),t);        %三维螺旋图
>> grid on;
>> axis square;
执行结果如图7-28所示。
  `, m3 D& Q9 p6 Y5 }- Z
# J( ?* K: Y- a2 A( _ 图7-28  三维螺旋图
% `) x/ n$ K) S! j2 |! ~0 E【实例分析】三维螺旋图的参数方程形式为x=sin(t),y=cos(t),z=t。
* a) o7 e9 U0 a6 Z

lili456

7.3.2  mesh——绘制三维网格图
【功能简介】绘制三维网格图。
【语法格式】
1.mesh(X,Y,Z)
* a. d' H2 M7 w7 E/ ^生成由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。颜色由高度决定。
& c! q* v/ z+ |3 k3.mesh(…,C)
图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.28】绘制三维网格图。
>> [X,Y]=meshgrid(-8:.5:8);        %构造X、Y矩阵
5 v! b8 J* J5 b0 N( v>> R=sqrt(X.^2+Y.^2)+eps;
  \9 X9 ~# e7 l% x- y$ p>> Z=sin(R)./R;
9 ]+ U: g% ?& v' r' D" Z>> mesh(X,Y,Z);                        %绘制三维网格图
执行结果如图7-29所示。
' }" A' a8 n+ q; N
6 u  Y. o& w9 r/ H4 d
/ ]; N2 P( R" R. Y8 A5 x% a' r图7-29  三维网格图
【实例分析】不指定颜色时，网格颜色由高度决定。

lili456

7.3.3  surf——三维曲面图
) C2 L( {) V& E' |; _4 X【功能简介】绘制三维曲面图。
【语法格式】
1.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函数一致。
2.surf (Z)
X与Y自动生成。[m,n]=size(Z)，则X=1:n，Y=1:m。颜色由高度决定。
; o9 b% ^8 m$ o+ k3 m  q3.surf (…,C)
* d7 A) d% W0 v7 D1 S0 M图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.29】绘制三维曲面图。
( }9 W. ?" `# F' H: T( n, ?, Y& B>> [X,Y]=meshgrid(-8:.5:8);        %生成X、Y矩阵
>> R=sqrt(X.^2+Y.^2)+eps;
8 T, e& j4 ?% Q2 l. E1 v$ D>> Z=sin(R)./R;
>> surf(X,Y,Z);                                %绘制三维曲面图
执行结果如图7-30所示。
+ f6 L' s# U+ P, z% N
8 X5 o! n6 E, m  s. g3 y. A图7-30  三维曲面图
【实例分析】surf函数与mesh函数调用格式是相同的。
' g2 P; m7 h, a! Z8 }8 H

lili456

7.3.4  contour3——三维等高线绘制
2 {* F" u+ O+ d1 @8 ?$ D【功能简介】绘制三维等高线。
2 I. Y5 m4 ]' Y) a) b* o- H2 }9 Q【语法格式】
+ b( I0 B; y( z( o1.contour3(Z)
6 [3 a5 f, }8 {- l  E, Q画出矩阵Z的等高线图，矩阵Z至少为2×2大小，系统按矩阵Z中元素距离XY平面的高度，自动选择等高线的条数进行绘制。X和Y的值是自动确定的，[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
- t2 n" h% f; M# Z' H: ^1 s" I* P格式变体：
; U( v: }" r! N- W+ {6 f        contour3(Z,n)：画出矩阵Z的n条等高线。
0 w6 Z) E* Q% W7 F' u        contour3(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
; g% G. l) X, r; b( W( x2.contour3(X,Y,Z)
画出矩阵Z的等高线图，X轴和Y轴的范围由参数X、Y指定。如果X为矩阵，则使用X(1,来定义X轴的范围，如果Y为矩阵，则使用Y(:,1)来定义Y轴的范围。如果X与Y均为矩阵，则两者必须同型。
$ d1 K5 A, ?! H) ^  N! u9 X格式变体：
        contour3(X,Y,Z,n)：画出矩阵Z的n条等高线。
        contour3(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
3.contour3(…,LineSpec)
参数LineSpec指定线型、标记符号和颜色。
8 m( x' k' e. A. s# k【实例7.30】绘制等高线。
% a- |+ z! ~+ e- \! h  n>> [X,Y]=meshgrid([-2:.25:2]);
' F6 j  C, H; d* B5 s>> Z=X.*exp(-X.^2-Y.^2);
( T# P0 d0 w; t7 r& |>> contour3(X,Y,Z,30);        %绘制等高线
( u0 X4 p. j+ U2 c# A+ A+ \>> colormap cool;
" y1 a9 p. t0 ?5 {5 O. L: b; ^1 F执行结果如图7-31所示。

. w7 J5 W) }/ D6 r 图7-31  contour3函数绘制等高线
; F0 U: s+ \8 o- p# Y  p【实例分析】等高线绘制有contour3和contour两个函数均可完成，contour函数绘制的是二维的等高线。

lili456

7.3.5  contour——曲面的等高线
【功能简介】绘制曲面的等高线。
$ ]: W  H5 @! P$ o( B2 d【语法格式】
4 k# ?$ ]- n/ z4 ]5 Z0 g' P2 P0 Y1.contour(Z)
画出矩阵Z的二维等高线图，矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
格式变体：
- N0 ~4 s6 Z3 b- Q5 s        contour(Z,n)：画出矩阵Z的n条等高线。
        contour(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
2.contour(X,Y,Z)
/ v6 W+ @) |. I9 R$ v画出矩阵Z的二维等高线图，X轴和Y轴的范围由参数X、Y指定。如果X与Y均为矩阵，则两者必须同型且单调递增。
格式变体：
        contour(X,Y,Z,n)：画出矩阵Z的n条等高线。
        contour(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
# z7 L! L5 R  d  R3.contour(…,LineSpec)
参数LineSpec指定线型、标记符号和颜色。
5 i$ I8 O6 @, v/ U【实例7.31】画出曲面的等高线图。
>> [X,Y]=meshgrid([-2:.25:2]);
>> Z=X.*exp(-X.^2-Y.^2);
>> contour(X,Y,Z,30);        %画出曲面的二维等高线
执行结果如图7-32所示。

图7-32  曲面的等高线
1 ~+ k5 s. s  B! N: l4 P【实例分析】contour函数画出的是二维图形，contour3函数画出的是三维图形。

lili456

7.3.6  clabel——等高线高度标签
【功能简介】在二维等高线图中绘制等高线高度标签。
6 a1 \9 m! q7 P1 o( I5 a【语法格式】
1.clabel(C,h)
/ ]6 N1 z$ n9 b8 M/ }' F* AC为标签矩阵，h为等高线句柄。函数旋转标签到适当角度，在空间允许的情况下插入标签。
  U1 b. Q, U, k4 f  L: `$ ~: R7 {格式变体：
        clabel(C,h,v)：在高度v处插入标签C。
        clabel(C,h,'manual')：手动添加标签到句柄h指定的等高线图中，用鼠标左键或空格键在最接近位置上放置标签，用回车键结束操作。
2.clabel(C)
添加标签到当前等高线图中，随机选择标签位置。
格式变体：
        clabel(C,v)：在v指定的高度处添加标签。
0 o4 G( N4 `" k        clabel(C,'manual')：手动添加标签到当前等高线图中。
. E& d8 Q& n% F  e【实例7.32】给等高线做标注。
>> [x,y,z]=peaks;
6 M) D! X! B" m>> C=contour(x,y,z);        %绘制等高线
% }. `8 y' M" b9 \! X! M; z>> clabel(C);                %给等高线添加标注
' k  B1 K' {. {2 [/ G$ O4 s! z执行结果如图7-33所示。
: K: }% c, X! Q# [- |9 c
8 N( D8 U5 F) Q- {* n图7-33  给等高线添加标注
【实例分析】peaks是MATLAB自带函数。
7 M# L! n4 ~) n; a7.3.7  contourc——等高线图形计算
【功能简介】计算等高线矩阵C，用函数contour、contour3和contourf来显示。
; F, U7 D) \6 G$ s3 v; s% f【语法格式】
1.C=contourc(Z)
从矩阵Z中计算等高矩阵C，Z可视为XY平面的高度矩阵。等高线的数量和高度值是系统自动确定的。Z至少为2×2矩阵，至少包含两个不同的值，X和Y的范围是1:n与1:m，其中[m,n]=size(Z)。
格式变体：
        C=contour(Z,n)：确定等高线条数为n，返回等高矩阵。
& @8 e5 @  S+ e4 x9 x        C=contour(Z,v)：在高度v处计算等高线。
1 W2 U' O/ L! \$ A2 g* L* n2.C=contourc(X,Y,Z)
在X和Y上计算Z的等高矩阵C，Z可视为XY平面的高度矩阵。等高线的数量和高度值是系统自动确定的。Z至少为2×2矩阵，至少包含两个不同的值。
( }: J6 o5 r& _0 I9 ?+ `+ M5 f! x格式变体：
1 l( h! n* I2 Q: b0 _) b  r        C=contour(X,Y,Z,n)：确定等高线条数为n，返回等高矩阵。
  {. K2 j- Y. e* L, M5 o        C=contour(X,Y,Z,v)：在高度v处计算等高线。
; e9 L% M. X1 t5 V7 l) V【实例7.33】计算peaks函数的等高矩阵。
>> a=peaks;
# v$ E$ ^; b0 _/ G>> c=contourc(a,10);                %计算10条等高线的等高矩阵
>> s=size(c)                        %矩阵维度
& [" {) y2 A8 E8 B- o, js =
     2   800
! h/ h1 n2 a7 k* A; ^>> contour(a,c)                        %画出等高线
# V. W' }& U+ e- W7 T执行结果如图7-34所示。

5 J2 g: I) q9 @0 U; f- a  p) D: Q
图7-34  计算等高矩阵再绘制等高线
【实例分析】矩阵C是一个2×m矩阵，m的值取决于数据矩阵及等高线的条数。

lili456

7.3.8  fill3——填充三维图
【功能简介】填充三维多边形。
' E) n& k- f: x( I  E【语法格式】
0 w8 g# Y( ]  `% f+ L- K' g& ?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)。
格式变体：
        fill3(X1,Y1,Z1,C1,X2,Y2,Z2,C2,…)：绘制多个三维多边形。
. ^! Q  D3 S1 T6 U  b, u6 V$ \' f2.fill3(X,Y,Z,ColorSpec)
! o& |  o# N7 j# r$ R- B用参数ColorSpec指定的颜色对三维多边形进行填充。
【实例7.34】用渐变的颜色来对三维多边形进行着色。
/ `0 N! P( a4 ~: O: `0 Z>> 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];
: L0 e4 |( {4 k1 i) Z>> C=[.5,1,1,.5;1,.5,.5,.1667;.3333,.3333,.5,.5]        %指定颜色
: F; q4 j2 f9 K7 f" s" w2 S% _) KC =
    0.5000    1.0000    1.0000    0.5000
. a5 }1 o4 A9 J+ ?9 H    1.0000    0.5000    0.5000    0.1667
8 [7 x3 S5 P. O3 x/ V" O    0.3333    0.3333    0.5000    0.5000
0 m, b* R9 a0 y8 J>> fill3(X,Y,Z,C)                                                %填充三维多边形
执行结果如图7-35所示。

图7-35  填充三维多边形
& @5 {, Y, o! {% H7 s9 [【实例分析】X、Y、Z决定三维多边形的形状，C确定颜色。
7.3.9  sphere——绘制球体
" [: z" O( {6 ~1 H; h" B2 d# p" Z【功能简介】绘制球体。
4 K! [* e7 Q5 z% S& t6 o3 \【语法格式】
1.sphere(n)
在三维直角坐标系中绘制由n×n个面组成的单位球体。
格式变体：
* U* p* E( N1 ]! W  N; Z8 s- n        sphere：默认n=20，绘制单位球体。
2.[X,Y,Z]=sphere(n)
. W' Q) J, e, e+ L5 r返回三个阶数为(n+1)×(n+1)的坐标矩阵。该命令不画图，只是返回矩阵，要画图可以调用mesh(X,Y,Z)或surf(X,Y,Z)。
【实例7.35】绘制多个球体。
>> [x,y,z]=sphere;
>> surf(x,y,z);                %画出第一个球体
>> hold on;
1 y. N+ z  A1 y+ l>> surf(x+3,y+2,z);        %画出第二个球体
5 t5 R# G0 a# x# F1 \7 k+ \>> surf(x,y-1,z+2);        %画出第三个球体
>> daspect([1 1 1])
执行结果如图7-36所示。
3 @7 K( _3 _$ M9 t  w3 V
" l, D+ s, s( a+ M; j  q7 F图7-36  绘制多个球体
【实例分析】sphere绘制单位球体，半径为1。
1 Q# Q' I0 |4 h5 I( H

lili456

7.3.10  contourf——填充二维等高线
% J1 L" {% d8 W8 C+ t& V, H, i4 y【功能简介】填充二维等高线。
8 e7 i* ^. C6 i$ W! r2 c$ \" g7 ]【语法格式】
9 t) l7 d, f3 x) c$ W8 ]1.contourf(Z)
画出矩阵Z的二维等高线图，再对等高线之间的区域进行填充，填充的颜色使用当前窗口的颜色映射表。矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
, V* c3 N0 g! N  k+ a格式变体：
* J* {% Q% t' s  a5 W' u" V3 n: R        contourf(Z,n)：画出矩阵Z的n条等高线并对等高线间的区域进行填充。
6 F, d. Y+ ?2 J7 F, s" Q        contourf(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
2.contourf(X,Y,Z)
! l7 a" A8 O; B1 S. E$ _; d0 L$ z" T画出矩阵Z的等高线图并进行填充，X轴和Y轴的范围由参数X、Y指定。如果X与Y均为矩阵，则两者必须同型且单调递增。
* g3 M6 O3 P% X- d* l格式变体：
        contourf(X,Y,Z,n)：画出矩阵Z的n条等高线并进行填充。
        contourf(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
  B# s  |/ V5 a+ t7 b【实例7.36】画出等高线并进行填充。
>> a=peaks;
4 ?5 f: }9 A2 ^* c6 W4 W: R  T>> contourf(a);                %画出用颜色填充的等高线
! r6 c; H4 _* A7 r7 Q% X  g执行结果如图7-37所示。

图7-37  用颜色填充的等高线图
【实例分析】用颜色对等高线间的区域进行填充，使图形更直观。
7.3.11  pie3——三维饼图
【功能简介】绘制三维饼图。
4 e& b$ l0 q7 O& ?【语法格式】
1.pie3(X)
使用X中的数据绘制一个三维饼图，X中的每一个元素都是饼图的一个部分。元素所占扇形面积的大小由X(i)/sum(X)决定，如果sum(X)<1，则只绘制不完整的三维饼形图。
2.pie3(X,explode)
指定三维饼图中的每一部分是否分离出来。explode是与X同型的数组，用非零值表示相应位置的X中的元素分离出来显示。
, g9 i: W7 H5 |【实例7.37】绘制三维饼图。
* G  {3 t8 J* W7 @& J2 S>> x=[1,3,.7,2.5,2];
>> explode=[0,1,0,0,1];        %设置第二个元素和最后一个元素对应的区域分离出来独立显示
$ w. Q: `# E* Q8 O! Z  r7 u! H! k( g( e>> pie3(x,explode);                %绘制三维饼图
( D; p7 ^- [. _2 Q0 C' V* t>> colormap hsv
6 ?7 q( V7 }2 ~6 t- t- `执行结果如图7-38所示。
  `6 |0 n. D% \7 q+ s
图7-38  三维饼图
【实例分析】pie3的调用格式与pie类似。

lili456

7.3.12  comet3——三维彗星图
【功能简介】绘制三维彗星图。
【语法格式】
$ O1 ~3 Z: J) |# B2 Q1.comet3 (z)
显示向量z的三维彗星图，X与Y轴范围由系统自动确定。
2.comet3(x,y,z)
显示由x、y、z确定的三维曲线的彗星图。
3.comet3(x,y,z,p)
指定彗星体的长度为p*length(y)。
' r6 s& p% T4 h0 l! B) }1 I( T【实例7.38】绘制一个三维彗星图。
/ _6 [# y6 [, B' s) ~1 E>> t=-10*pi:pi/250:10*pi;
5 W- ~( {6 [$ U: P  r, d  a) |( D0 C>> comet3((cos(2*t).^2).*sin(t),(sin(2*t).^2).*cos(t),t);                %绘制彗星图
执行结果如图7-39所示。
* h# o. H4 q. ]7 |2 s: A
6 g+ F5 t" A! h4 p7 ~图7-39  三维彗星图
4 p- k* T" G6 t0 B  F【实例分析】彗星图的绘制是一个动态的过程，用户可按实例中的代码自行运行。
) `/ s! A4 Z. f0 p$ U3 T0 V7.3.13  cylinder——生成圆柱图形
6 i' ~+ u# v/ s* U7 r【功能简介】绘制圆柱图形。
' s. O0 j% w- W% e# }; c  c1.[X,Y,Z]=cylinder(r,n)
8 h1 j3 m7 `3 C返回一个半径为r、高度为1的圆柱体的X、Y、Z轴坐标值，所绘制的圆柱体的圆周有指定的n个距离相同的点。用户可用surf或mesh函数和X、Y、Z参数画出圆柱。
格式变体：
        [X,Y,Z]=cylinder(r)：所绘制的圆柱体的圆周有指定的20个距离相同的点。
+ I9 v( J3 z" R4 P" A" h( ~        [X,Y,Z]= cylinder：半径采用默认值1。
2.cylinder(…)
没有输出参量时，直接画出圆柱体。
3 }( V' Z+ s1 }) H0 \【实例7.39】绘制一个两头粗、中间细的圆柱。
; U& w( P" d1 r2 B; O, d>> t=0:pi/10:2*pi;
. S1 K1 t/ p5 s* a  X>> [X,Y,Z]=cylinder(2+cos(t));                %用cylinder函数算出X、Y、Z参数
>> surf(X,Y,Z);                                        %画出圆柱体
>> axis square                                        %调整坐标轴
: n5 G; C6 `9 A执行结果如图7-40所示。
; U1 V2 S. q: x$ j+ K) j- j. W  D
% u, o# j/ [/ @9 ^4 s, f8 w6 e图7-40  绘制圆柱体
【实例分析】也可不带返回值，直接用cylinder(2+cos(t))画出圆柱。

lili456

7.3.14  surfc——绘制阴影图及等高线
! @6 a/ d. M2 y【功能简介】绘制阴影图及等高线。先画出三维曲面图，再在下方画出二维等高线。
! x( T. ?: r- b. e. A% `" L【语法格式】
3 I" i* e4 w, ]; k* d8 g4 m4 r1.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函数一致。
2.surfc(Z)
X与Y自动生成。[m,n]=size(Z)，则X=1:n，Y=1:m。颜色由高度决定。
7 H4 j( f- ~$ W* g8 K) D3.surfc(…,C)
$ A6 u1 _) P( X/ _% l图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.40】绘制三维曲面及其二维等高线        。
>> [X,Y]=meshgrid(-8:.5:8);
>> R=sqrt(X.^2+Y.^2)+eps;
>> Z=sin(R)./R;
>> surfc (X,Y,Z);                        %绘制三维曲面及其二维等高线               
$ c  N5 U+ T+ X% `* I( \7 m. Q8 K执行结果如图7-41所示。
9 r5 k" g7 S0 V0 v; s" \1 P
. B6 n2 S+ o4 S& }# {; E图7-41  绘制三维曲面及其二维等高线
; J8 M) z6 Z6 q, s, R【实例分析】suefc函数调用格式与surf相同，只是多画了一个二维等高线。
3 J. {1 P) V4 d1 G; N7 l% v7.3.15  surfl——带光照模式的曲面图
- f# C  R3 s# ~. W: K【功能简介】绘制带光照模式的曲面图。
【语法格式】
) {/ s6 }# a# B' e! g1.surfl(X,Y,Z)
& R& q6 m1 s, |4 ?5 D生成由X、Y和Z定义的曲面图，光照方向和系数采用默认值。
格式变体：
; F! y1 @1 n; q' N7 X        surfl (Z)：X、Y的值由系统自动生成，等于Z中元素的索引。
0 L; N! O6 e8 G2 M& k# B+ p2 W2.surfl(…,'light')
用MATLAB光照对象生成一个带光照的曲面。
/ N& W2 _  b# L6 x3.surfc(…,s)
指定光源的方向s。s是长度为2或3的向量，即s=[azimuth,elevation]或[sx,sy,sz]。默认光源方向是从当前视角开始逆时针45度方向。
9 s9 D% ^3 t( g3 U4.surfc(…,s,k)
  k  r8 _! _* R指定反射系数k，k是一个常量。k为长度为4的向量[ka,kd,ks,shine]，四个元素分别表示环境光系数、漫反射系数、镜面反射系数和镜面反射亮度，默认值为[0.55,0.6,0.4,10]。
【实例7.41】对peaks函数表示的曲面用surfl函数进行描绘。
>> [x,y]=meshgrid(-3:1/8:3);
7 i: |+ X) R( a7 F/ m) ]>> z=peaks(x,y);                %peaks函数
+ C1 F; H8 Q  A. U$ f2 ?>> surfl(x,y,z);                        %绘制带光照的曲面图
- ]0 p7 G7 E9 P6 A5 A2 B$ `8 w>> shading interp;
>> colormap gray
执行结果如图7-42所示。

图7-42  绘制带光照的曲面图
8 j& [+ O& Z* a: r【实例分析】命令colormap gray将色图指定为灰度色图。