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

lili456

7.3  三维图形
三维图形表现的是三个变量之间的依赖关系，使用户能够直观看出二元函数数值的变化趋势。下面的小节将介绍三维图形的绘制函数。
7.3.1  plot3——绘制三维曲线
【功能简介】绘制三维曲线。
  x5 V3 E; i4 z2 _4 p5 v' m4 u【语法格式】
' g$ J9 R6 A+ z5 B; _/ u1.plot3(X1,Y1,Z1,…)
X1、Y1、Z1是同型的向量或矩阵，函数在三维空间中绘制出一条或多条曲线，绘制多条曲线时，曲线的条数等于矩阵的列数。
! y3 e3 Q3 V9 U! g" H# U格式变体：
        plot3(X1,Y1,Z1,LineSpec)：用参数LineSpec指定画线的线型、标记符号和颜色。
- X' e$ r- y$ n2.plot3(…,'PropertyName','PropertyValue')
& z& e% G3 S- `! n3 t  F' @设置所绘图形指定属性的属性值。
【实例7.27】绘制三维螺旋图。
>> t=0:pi/50:10*pi;
>> plot3(sin(t),cos(t),t);        %三维螺旋图
$ B' k; Y% c) `, n+ L) b0 a>> grid on;
>> axis square;
  ?7 F) n, W; w/ Z+ D执行结果如图7-28所示。

( s5 l. {5 V& S0 y: w- @, W* ^ 图7-28  三维螺旋图
# V/ ?3 ^. {# ]9 @6 N【实例分析】三维螺旋图的参数方程形式为x=sin(t),y=cos(t),z=t。

lili456

7.3.2  mesh——绘制三维网格图
8 ]) Z5 `$ {' Y【功能简介】绘制三维网格图。
【语法格式】
1.mesh(X,Y,Z)
1 h4 B* I- f; ~/ k生成由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! T  d* j; p* kX与Y自动生成。[m,n]=size(Z)，则X=1:n，Y=1:m。颜色由高度决定。
3.mesh(…,C)
图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
7 z0 w) b  V7 i3 [! D. }6 _( E【实例7.28】绘制三维网格图。
  ?/ G0 p! V0 z& q>> [X,Y]=meshgrid(-8:.5:8);        %构造X、Y矩阵
>> R=sqrt(X.^2+Y.^2)+eps;
% C  S$ }) E0 P6 r$ N>> Z=sin(R)./R;
>> mesh(X,Y,Z);                        %绘制三维网格图
+ K0 `1 H, U2 [. p3 U& x6 G- E执行结果如图7-29所示。
2 c& T, ^8 X+ I; e
0 W! P2 }8 J+ t: {7 Q  k
3 Y3 j$ O! Y) ?# E3 E1 n图7-29  三维网格图
【实例分析】不指定颜色时，网格颜色由高度决定。

lili456

7.3.3  surf——三维曲面图
* I% Q+ W* j! y9 t5 @' t【功能简介】绘制三维曲面图。
% g- b8 P4 k; A- C, \8 m2 G【语法格式】
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。颜色由高度决定。
3.surf (…,C)
0 a3 c4 K$ R( O4 @6 J, T. I2 N图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.29】绘制三维曲面图。
$ Y* ]7 M0 Z" g6 l8 C/ `2 M7 ^>> [X,Y]=meshgrid(-8:.5:8);        %生成X、Y矩阵
>> R=sqrt(X.^2+Y.^2)+eps;
8 R+ n  A% W, M& n>> Z=sin(R)./R;
>> surf(X,Y,Z);                                %绘制三维曲面图
5 q& V7 h/ ~! `, M执行结果如图7-30所示。

. W: h0 u2 o3 [6 A1 d. J- P图7-30  三维曲面图
【实例分析】surf函数与mesh函数调用格式是相同的。

lili456

7.3.4  contour3——三维等高线绘制
【功能简介】绘制三维等高线。
. {! ^5 L. e: p: {# l5 {$ d. U+ G【语法格式】
1.contour3(Z)
画出矩阵Z的等高线图，矩阵Z至少为2×2大小，系统按矩阵Z中元素距离XY平面的高度，自动选择等高线的条数进行绘制。X和Y的值是自动确定的，[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
: ]% u* T2 e7 c' [0 X! z1 ]格式变体：
: r9 l: k4 S# Y; m& I        contour3(Z,n)：画出矩阵Z的n条等高线。
: V' i9 {5 Y  G  |        contour3(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
2.contour3(X,Y,Z)
5 N! q. S3 Z4 E# W1 O9 q5 W画出矩阵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)
- D& s- O. K* I$ D参数LineSpec指定线型、标记符号和颜色。
- O2 {5 O+ }  x- `$ x! H【实例7.30】绘制等高线。
( Q7 U( a' y% r; O' R>> [X,Y]=meshgrid([-2:.25:2]);
4 S  Y9 j0 D  |1 K8 S( o. i2 y9 Q+ M>> Z=X.*exp(-X.^2-Y.^2);
>> contour3(X,Y,Z,30);        %绘制等高线
>> colormap cool;
执行结果如图7-31所示。
) w* i: }4 N, }9 }. S
图7-31  contour3函数绘制等高线
【实例分析】等高线绘制有contour3和contour两个函数均可完成，contour函数绘制的是二维的等高线。

lili456

7.3.5  contour——曲面的等高线
【功能简介】绘制曲面的等高线。
【语法格式】
1.contour(Z)
画出矩阵Z的二维等高线图，矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
% y1 n- B# t* L% m. w( E格式变体：
6 N1 l  ^/ D" `& f4 c! I3 k/ }' N        contour(Z,n)：画出矩阵Z的n条等高线。
        contour(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
2.contour(X,Y,Z)
5 T3 I: o7 ^* L) y' Z- a画出矩阵Z的二维等高线图，X轴和Y轴的范围由参数X、Y指定。如果X与Y均为矩阵，则两者必须同型且单调递增。
# h+ P  ~. h! o6 {. F9 ~格式变体：
        contour(X,Y,Z,n)：画出矩阵Z的n条等高线。
        contour(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
" ?, W: Y- P( z" F/ U3.contour(…,LineSpec)
: ~  O8 P! T2 I参数LineSpec指定线型、标记符号和颜色。
) v6 u. Z# Q4 Z; C) K, o; J【实例7.31】画出曲面的等高线图。
/ A7 D: S( ^9 x- b>> [X,Y]=meshgrid([-2:.25:2]);
>> Z=X.*exp(-X.^2-Y.^2);
- y' S, o& }/ [>> contour(X,Y,Z,30);        %画出曲面的二维等高线
: ]$ b) ^+ j, Q' e执行结果如图7-32所示。

图7-32  曲面的等高线
2 V0 n9 n% t. R8 I6 ?【实例分析】contour函数画出的是二维图形，contour3函数画出的是三维图形。
" C" M- I3 i: [- T8 x8 z

lili456

7.3.6  clabel——等高线高度标签
【功能简介】在二维等高线图中绘制等高线高度标签。
8 U2 B3 a* z) o. z【语法格式】
1.clabel(C,h)
; M3 R/ Q1 X$ _C为标签矩阵，h为等高线句柄。函数旋转标签到适当角度，在空间允许的情况下插入标签。
格式变体：
        clabel(C,h,v)：在高度v处插入标签C。
$ K7 T2 w1 A: h. z        clabel(C,h,'manual')：手动添加标签到句柄h指定的等高线图中，用鼠标左键或空格键在最接近位置上放置标签，用回车键结束操作。
2.clabel(C)
6 Q3 V4 A7 N/ z" M% }6 B/ \0 @* Q添加标签到当前等高线图中，随机选择标签位置。
+ J. S% q0 }1 _" ^; ]; r) [格式变体：
' ?) R% ?: ?+ V. G3 @& r        clabel(C,v)：在v指定的高度处添加标签。
        clabel(C,'manual')：手动添加标签到当前等高线图中。
【实例7.32】给等高线做标注。
>> [x,y,z]=peaks;
>> C=contour(x,y,z);        %绘制等高线
6 p  o' Y4 O9 ~>> clabel(C);                %给等高线添加标注
执行结果如图7-33所示。
6 P8 y, b# J. E2 u" K0 k
0 z( y- K! e$ a) t. N7 o0 B* d图7-33  给等高线添加标注
【实例分析】peaks是MATLAB自带函数。
7.3.7  contourc——等高线图形计算
& M1 e0 O- W2 l' x【功能简介】计算等高线矩阵C，用函数contour、contour3和contourf来显示。
【语法格式】
) w5 l, w2 n, i2 S  r! ?1.C=contourc(Z)
8 Z& V9 u& L  V( s$ g* F- K从矩阵Z中计算等高矩阵C，Z可视为XY平面的高度矩阵。等高线的数量和高度值是系统自动确定的。Z至少为2×2矩阵，至少包含两个不同的值，X和Y的范围是1:n与1:m，其中[m,n]=size(Z)。
格式变体：
4 b7 ]' i9 z3 j1 G; }        C=contour(Z,n)：确定等高线条数为n，返回等高矩阵。
        C=contour(Z,v)：在高度v处计算等高线。
2.C=contourc(X,Y,Z)
在X和Y上计算Z的等高矩阵C，Z可视为XY平面的高度矩阵。等高线的数量和高度值是系统自动确定的。Z至少为2×2矩阵，至少包含两个不同的值。
3 d2 q. r6 |) V, c; P$ q格式变体：
        C=contour(X,Y,Z,n)：确定等高线条数为n，返回等高矩阵。
        C=contour(X,Y,Z,v)：在高度v处计算等高线。
【实例7.33】计算peaks函数的等高矩阵。
* v- b& J3 O; R% N5 |>> a=peaks;
>> c=contourc(a,10);                %计算10条等高线的等高矩阵
! T3 ~- L+ o1 o! i# b$ q+ D  j. G>> s=size(c)                        %矩阵维度
: {, R" X: b$ W7 z! {, cs =
/ ~% u0 z# h! O7 q     2   800
9 Q: a  ^( @# O0 E2 f8 r8 P# }. Q% }9 ?>> contour(a,c)                        %画出等高线
6 K& u; G2 M- I% c) v, e  ^1 P; R执行结果如图7-34所示。


图7-34  计算等高矩阵再绘制等高线
( ^& l% W  D7 p4 {) k' c8 M, ?$ b【实例分析】矩阵C是一个2×m矩阵，m的值取决于数据矩阵及等高线的条数。
+ J4 P; u# e, N/ [8 H5 {

lili456

7.3.8  fill3——填充三维图
【功能简介】填充三维多边形。
6 T5 h0 Z2 X6 d0 N% z6 I【语法格式】
1.fill3(X,Y,Z,C)
; _5 R" B( M! e7 L参数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)。
8 e$ T# }* L0 _$ m6 W$ a格式变体：
        fill3(X1,Y1,Z1,C1,X2,Y2,Z2,C2,…)：绘制多个三维多边形。
+ v5 ]. h: R8 J1 f8 G! a3 v* h3 j2.fill3(X,Y,Z,ColorSpec)
6 }" p5 B+ g( ]$ _; n2 o" a用参数ColorSpec指定的颜色对三维多边形进行填充。
5 W" F! b; Z2 N【实例7.34】用渐变的颜色来对三维多边形进行着色。
: t+ }1 C! w9 h>> 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];
. b; ^* G7 t, p* S0 j$ c8 A>> Z=[1,1,1,1;1,0,1,0;0,0,0,0];
/ `& f7 s; f5 v1 [! K# S6 B>> C=[.5,1,1,.5;1,.5,.5,.1667;.3333,.3333,.5,.5]        %指定颜色
C =
    0.5000    1.0000    1.0000    0.5000
    1.0000    0.5000    0.5000    0.1667
    0.3333    0.3333    0.5000    0.5000
>> fill3(X,Y,Z,C)                                                %填充三维多边形
执行结果如图7-35所示。

图7-35  填充三维多边形
1 q) {) C/ @2 ]$ j: P! M# ]【实例分析】X、Y、Z决定三维多边形的形状，C确定颜色。
5 b/ O$ |3 l4 V  y& w$ u% G$ d7.3.9  sphere——绘制球体
# l% c4 z- b) R: v【功能简介】绘制球体。
5 i5 R. t/ ~  X  n8 m1 e- H【语法格式】
1.sphere(n)
3 m3 F4 Q0 K/ M1 X在三维直角坐标系中绘制由n×n个面组成的单位球体。
$ m3 Z5 h" Y  @  F0 _格式变体：
        sphere：默认n=20，绘制单位球体。
  L. }& y2 A5 q0 C/ F* a- {8 d, h& P2.[X,Y,Z]=sphere(n)
返回三个阶数为(n+1)×(n+1)的坐标矩阵。该命令不画图，只是返回矩阵，要画图可以调用mesh(X,Y,Z)或surf(X,Y,Z)。
, E" ]. T9 B. a2 n【实例7.35】绘制多个球体。
>> [x,y,z]=sphere;
>> surf(x,y,z);                %画出第一个球体
) ^7 H3 y7 G& Z* y3 l>> hold on;
>> surf(x+3,y+2,z);        %画出第二个球体
6 E& G$ E! o4 a$ B>> surf(x,y-1,z+2);        %画出第三个球体
>> daspect([1 1 1])
执行结果如图7-36所示。

图7-36  绘制多个球体
【实例分析】sphere绘制单位球体，半径为1。

lili456

7.3.10  contourf——填充二维等高线
7 q, V% K9 L! j4 d: O5 F【功能简介】填充二维等高线。
【语法格式】
1.contourf(Z)
6 A% ^0 w$ H! D5 ]0 y: Q( m画出矩阵Z的二维等高线图，再对等高线之间的区域进行填充，填充的颜色使用当前窗口的颜色映射表。矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
8 Z8 Q8 d- G8 |" _3 M% i: a格式变体：
7 Y* g9 G5 E8 }5 [        contourf(Z,n)：画出矩阵Z的n条等高线并对等高线间的区域进行填充。
        contourf(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
  C8 N+ M* G4 U6 p7 c2 J2.contourf(X,Y,Z)
& D! L4 A- C5 H! i, O) s' t, q  d画出矩阵Z的等高线图并进行填充，X轴和Y轴的范围由参数X、Y指定。如果X与Y均为矩阵，则两者必须同型且单调递增。
格式变体：
        contourf(X,Y,Z,n)：画出矩阵Z的n条等高线并进行填充。
2 V  U. {) r$ x* U        contourf(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
! L/ p' V  h! m, Q. @【实例7.36】画出等高线并进行填充。
>> a=peaks;
9 [  z& n- }% u1 i- z$ R* I>> contourf(a);                %画出用颜色填充的等高线
7 U; I4 g6 n( f$ O执行结果如图7-37所示。
" C) }0 Q, t% T1 C! J0 j
+ R1 p% C! f4 T' Y6 h图7-37  用颜色填充的等高线图
【实例分析】用颜色对等高线间的区域进行填充，使图形更直观。
( D* z3 z! b% e0 }, Y7.3.11  pie3——三维饼图
【功能简介】绘制三维饼图。
2 Y' k& b% \5 q* r! s( b【语法格式】
6 o2 o6 }0 e. Y; [1.pie3(X)
* c# S( ~& I1 n! P- ^使用X中的数据绘制一个三维饼图，X中的每一个元素都是饼图的一个部分。元素所占扇形面积的大小由X(i)/sum(X)决定，如果sum(X)<1，则只绘制不完整的三维饼形图。
; _/ s, K. {/ u2 ~8 R0 e6 h2.pie3(X,explode)
( z$ u( M) }4 `" p* A" {指定三维饼图中的每一部分是否分离出来。explode是与X同型的数组，用非零值表示相应位置的X中的元素分离出来显示。
【实例7.37】绘制三维饼图。
>> x=[1,3,.7,2.5,2];
( Z: M$ n+ o  Z6 s( [9 U>> explode=[0,1,0,0,1];        %设置第二个元素和最后一个元素对应的区域分离出来独立显示
>> pie3(x,explode);                %绘制三维饼图
- {4 F4 j. Y6 p1 M9 j>> colormap hsv
执行结果如图7-38所示。
  V, u' _* E8 [. m: q/ A
4 j* ]4 J; u' B! l% E图7-38  三维饼图
【实例分析】pie3的调用格式与pie类似。

lili456

7.3.12  comet3——三维彗星图
【功能简介】绘制三维彗星图。
7 Q9 Y5 L: U$ b7 }【语法格式】
6 C; L" [0 L, W& A  J, @1.comet3 (z)
显示向量z的三维彗星图，X与Y轴范围由系统自动确定。
2.comet3(x,y,z)
; w- \) U/ m3 _9 O* E0 t. b) G显示由x、y、z确定的三维曲线的彗星图。
% p/ L! N( x3 {) Z1 L' f! m3.comet3(x,y,z,p)
指定彗星体的长度为p*length(y)。
" R" I1 p$ H) u4 L! X1 \【实例7.38】绘制一个三维彗星图。
>> t=-10*pi:pi/250:10*pi;
: L: d1 q* Y* k* }' @>> comet3((cos(2*t).^2).*sin(t),(sin(2*t).^2).*cos(t),t);                %绘制彗星图
执行结果如图7-39所示。
* x' c  U" q9 I  W( y1 E! g
图7-39  三维彗星图
【实例分析】彗星图的绘制是一个动态的过程，用户可按实例中的代码自行运行。
8 \$ M& w8 {; z# }9 J7 V7.3.13  cylinder——生成圆柱图形
; b2 V: Z# ]( a! |【功能简介】绘制圆柱图形。
2 B" M% w* w0 x% V6 j; c+ l; M1.[X,Y,Z]=cylinder(r,n)
返回一个半径为r、高度为1的圆柱体的X、Y、Z轴坐标值，所绘制的圆柱体的圆周有指定的n个距离相同的点。用户可用surf或mesh函数和X、Y、Z参数画出圆柱。
格式变体：
        [X,Y,Z]=cylinder(r)：所绘制的圆柱体的圆周有指定的20个距离相同的点。
% p4 p' I3 x. C. l& `% k& A5 |        [X,Y,Z]= cylinder：半径采用默认值1。
2.cylinder(…)
8 r" c0 ?0 v" R% a/ r4 e没有输出参量时，直接画出圆柱体。
& i1 s* K+ F* V# G  o+ `【实例7.39】绘制一个两头粗、中间细的圆柱。
>> t=0:pi/10:2*pi;
! y* B# }7 V6 x' ?3 o, j3 N7 D1 @>> [X,Y,Z]=cylinder(2+cos(t));                %用cylinder函数算出X、Y、Z参数
>> surf(X,Y,Z);                                        %画出圆柱体
>> axis square                                        %调整坐标轴
0 P( ?9 q# i9 Z( `" F$ ]  f执行结果如图7-40所示。

* l9 }% |% v$ \图7-40  绘制圆柱体
【实例分析】也可不带返回值，直接用cylinder(2+cos(t))画出圆柱。
; F/ N& P: m* [2 m! }

lili456

7.3.14  surfc——绘制阴影图及等高线
# p0 p2 _8 n5 r3 G/ H2 }' N【功能简介】绘制阴影图及等高线。先画出三维曲面图，再在下方画出二维等高线。
* U/ C1 s' J. `: _8 Q- c【语法格式】
1.surfc(X,Y,Z)
* `* q9 N9 i( M- ^2 t生成由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函数一致。
- U# t% q+ s9 X2.surfc(Z)
3 e( R9 u  i" e+ V+ S" J( `, N2 l: }X与Y自动生成。[m,n]=size(Z)，则X=1:n，Y=1:m。颜色由高度决定。
3 k9 |( W& ?  m; U# \2 a2 J- |5 ~3.surfc(…,C)
图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.40】绘制三维曲面及其二维等高线        。
6 O2 M. c, e$ e) S>> [X,Y]=meshgrid(-8:.5:8);
- |. y0 U( y, Z0 W/ X>> R=sqrt(X.^2+Y.^2)+eps;
- C4 Q6 I" P- ^# m! Q# H9 z: Y>> Z=sin(R)./R;
>> surfc (X,Y,Z);                        %绘制三维曲面及其二维等高线               
/ q+ O$ s) J) R$ w' a9 n) J7 X- b8 P执行结果如图7-41所示。

# O  L4 m  S& M: J0 L图7-41  绘制三维曲面及其二维等高线
【实例分析】suefc函数调用格式与surf相同，只是多画了一个二维等高线。
7.3.15  surfl——带光照模式的曲面图
【功能简介】绘制带光照模式的曲面图。
9 S& ^, z* y- B1 h; Y# a【语法格式】
' f* H, z* O3 V' X$ b1.surfl(X,Y,Z)
/ t! z. i: F9 B% a3 s/ @生成由X、Y和Z定义的曲面图，光照方向和系数采用默认值。
格式变体：
8 Y- G  ]$ W) h4 @        surfl (Z)：X、Y的值由系统自动生成，等于Z中元素的索引。
2.surfl(…,'light')
5 ^$ r- t& ?0 a, k0 d用MATLAB光照对象生成一个带光照的曲面。
8 I. N% I; r. |' a3 x6 f; N. f3.surfc(…,s)
指定光源的方向s。s是长度为2或3的向量，即s=[azimuth,elevation]或[sx,sy,sz]。默认光源方向是从当前视角开始逆时针45度方向。
4.surfc(…,s,k)
指定反射系数k，k是一个常量。k为长度为4的向量[ka,kd,ks,shine]，四个元素分别表示环境光系数、漫反射系数、镜面反射系数和镜面反射亮度，默认值为[0.55,0.6,0.4,10]。
【实例7.41】对peaks函数表示的曲面用surfl函数进行描绘。
* R, S& W! X& Z( j7 y. T6 f>> [x,y]=meshgrid(-3:1/8:3);
>> z=peaks(x,y);                %peaks函数
>> surfl(x,y,z);                        %绘制带光照的曲面图
>> shading interp;
>> colormap gray
执行结果如图7-42所示。
' F' x. [2 e+ g) ~
图7-42  绘制带光照的曲面图
【实例分析】命令colormap gray将色图指定为灰度色图。
# F- n  S3 p+ j( ]