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

lili456

7.3  三维图形
6 O1 t- c0 P& m* r) {. N3 N三维图形表现的是三个变量之间的依赖关系，使用户能够直观看出二元函数数值的变化趋势。下面的小节将介绍三维图形的绘制函数。
  D) C, A+ J! V7.3.1  plot3——绘制三维曲线
【功能简介】绘制三维曲线。
& U5 j: y, ?, |+ T, K% C【语法格式】
1.plot3(X1,Y1,Z1,…)
+ y- ^8 z: O* x( v. }+ V  z# u9 OX1、Y1、Z1是同型的向量或矩阵，函数在三维空间中绘制出一条或多条曲线，绘制多条曲线时，曲线的条数等于矩阵的列数。
7 ^+ r9 W$ I& L- i格式变体：
7 P8 o) ]$ o/ ?0 ?" J3 B0 j' z        plot3(X1,Y1,Z1,LineSpec)：用参数LineSpec指定画线的线型、标记符号和颜色。
2.plot3(…,'PropertyName','PropertyValue')
设置所绘图形指定属性的属性值。
6 |3 G9 f, \" ~% l【实例7.27】绘制三维螺旋图。
9 u2 |" ]0 s) e+ Q- |6 [>> t=0:pi/50:10*pi;
' [' t$ _5 N( V3 n- C! b>> plot3(sin(t),cos(t),t);        %三维螺旋图
>> grid on;
/ U! l: L4 I9 y1 h4 ~>> axis square;
8 O2 g2 ^3 d# d: s2 J# b1 X6 x执行结果如图7-28所示。

0 y1 ]& ?; D1 ~ 图7-28  三维螺旋图
【实例分析】三维螺旋图的参数方程形式为x=sin(t),y=cos(t),z=t。
* i& j6 ~2 r8 Y) F2 g/ B1 _

lili456

7.3.2  mesh——绘制三维网格图
【功能简介】绘制三维网格图。
【语法格式】
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。颜色由高度决定。
6 n4 E; t5 j4 P. [5 `3.mesh(…,C)
图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
0 S% u/ W: a" p5 H+ Y( A: e【实例7.28】绘制三维网格图。
& X' H1 t9 d# a  p0 e4 ?>> [X,Y]=meshgrid(-8:.5:8);        %构造X、Y矩阵
>> R=sqrt(X.^2+Y.^2)+eps;
# J3 |6 s7 n4 x9 h+ w. g. m>> Z=sin(R)./R;
  f# H2 I$ Z, @# e& f>> mesh(X,Y,Z);                        %绘制三维网格图
+ U8 B# r" J5 x5 g4 M% j+ @执行结果如图7-29所示。
( {' N- ~2 A/ X) o, R5 M6 N
7 \' j- U) n. T6 p( j: Z
图7-29  三维网格图
5 e( l; d1 {  E9 ~【实例分析】不指定颜色时，网格颜色由高度决定。
( z  }! W% S8 v9 U0 p9 t1 e

lili456

7.3.3  surf——三维曲面图
4 U+ m4 g; Z/ i4 B0 t% g* J$ h6 z【功能简介】绘制三维曲面图。
【语法格式】
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)
7 @7 A9 K- e  {% i# X) D. ^" K% l图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.29】绘制三维曲面图。
>> [X,Y]=meshgrid(-8:.5:8);        %生成X、Y矩阵
>> R=sqrt(X.^2+Y.^2)+eps;
+ G$ U1 G+ [5 d3 g" r" p1 @>> Z=sin(R)./R;
' N' \; `- e5 g, B) `>> surf(X,Y,Z);                                %绘制三维曲面图
* M8 u4 Z" p3 q8 A# Z9 ^) [6 k, v执行结果如图7-30所示。

图7-30  三维曲面图
【实例分析】surf函数与mesh函数调用格式是相同的。

lili456

7.3.4  contour3——三维等高线绘制
* C7 C3 a5 L  U  D9 |【功能简介】绘制三维等高线。
# U& ^: M- x+ F【语法格式】
' B6 q0 M+ X( c1.contour3(Z)
; l! S+ E- B: C画出矩阵Z的等高线图，矩阵Z至少为2×2大小，系统按矩阵Z中元素距离XY平面的高度，自动选择等高线的条数进行绘制。X和Y的值是自动确定的，[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
5 j6 V: x: M! R: `  k格式变体：
- g% [7 x3 G% U* \        contour3(Z,n)：画出矩阵Z的n条等高线。
1 h1 ^( T2 `0 u) W        contour3(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
2 [8 h2 _6 R8 w- B; h2.contour3(X,Y,Z)
画出矩阵Z的等高线图，X轴和Y轴的范围由参数X、Y指定。如果X为矩阵，则使用X(1,来定义X轴的范围，如果Y为矩阵，则使用Y(:,1)来定义Y轴的范围。如果X与Y均为矩阵，则两者必须同型。
5 E5 ~4 K7 J8 Y/ t$ m& M: r/ l格式变体：
+ m2 L  Q" _5 {' H        contour3(X,Y,Z,n)：画出矩阵Z的n条等高线。
        contour3(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
3.contour3(…,LineSpec)
参数LineSpec指定线型、标记符号和颜色。
【实例7.30】绘制等高线。
>> [X,Y]=meshgrid([-2:.25:2]);
; R, O1 P$ U7 c- `& r: _% c>> Z=X.*exp(-X.^2-Y.^2);
>> contour3(X,Y,Z,30);        %绘制等高线
>> colormap cool;
执行结果如图7-31所示。
6 W# @2 }0 |3 r, Q) q: ?5 [
5 A$ s& C7 x' v 图7-31  contour3函数绘制等高线
! f2 Y4 z5 n5 d【实例分析】等高线绘制有contour3和contour两个函数均可完成，contour函数绘制的是二维的等高线。
/ r; O" g- @# t7 ?" L

lili456

7.3.5  contour——曲面的等高线
, y+ a, M# @' o0 J2 f% ]【功能简介】绘制曲面的等高线。
【语法格式】
1.contour(Z)
画出矩阵Z的二维等高线图，矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
- D( I- n# x4 @" n' x6 ]* F格式变体：
        contour(Z,n)：画出矩阵Z的n条等高线。
        contour(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
2.contour(X,Y,Z)
画出矩阵Z的二维等高线图，X轴和Y轴的范围由参数X、Y指定。如果X与Y均为矩阵，则两者必须同型且单调递增。
格式变体：
        contour(X,Y,Z,n)：画出矩阵Z的n条等高线。
: w7 [3 W; y' [  @+ F. i7 T6 Q        contour(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
$ T4 t6 x, [9 c  S- }) p3.contour(…,LineSpec)
9 M9 [3 U7 d! F参数LineSpec指定线型、标记符号和颜色。
【实例7.31】画出曲面的等高线图。
>> [X,Y]=meshgrid([-2:.25:2]);
>> Z=X.*exp(-X.^2-Y.^2);
  }' N+ z- g* M% \1 w) i6 e! c>> contour(X,Y,Z,30);        %画出曲面的二维等高线
执行结果如图7-32所示。
6 V7 M5 z: W9 `
! v, ]% e8 P1 m4 x 图7-32  曲面的等高线
【实例分析】contour函数画出的是二维图形，contour3函数画出的是三维图形。

lili456

7.3.6  clabel——等高线高度标签
# p! R1 H8 G- W9 A【功能简介】在二维等高线图中绘制等高线高度标签。
. B# b- s, E/ V# }- U, {& k【语法格式】
1.clabel(C,h)
C为标签矩阵，h为等高线句柄。函数旋转标签到适当角度，在空间允许的情况下插入标签。
' @; o# C# @% h) |0 E格式变体：
        clabel(C,h,v)：在高度v处插入标签C。
        clabel(C,h,'manual')：手动添加标签到句柄h指定的等高线图中，用鼠标左键或空格键在最接近位置上放置标签，用回车键结束操作。
2.clabel(C)
添加标签到当前等高线图中，随机选择标签位置。
格式变体：
! D5 `+ g& V2 \9 T5 w        clabel(C,v)：在v指定的高度处添加标签。
        clabel(C,'manual')：手动添加标签到当前等高线图中。
( y; u/ E6 h7 D( o6 E【实例7.32】给等高线做标注。
>> [x,y,z]=peaks;
! z! J, Z: r' C& j+ k1 n* }8 i' R3 ?! N>> C=contour(x,y,z);        %绘制等高线
' R  P1 V# P. }  G# ~>> clabel(C);                %给等高线添加标注
" z$ l4 H) Q. a! j执行结果如图7-33所示。

7 p" m; B4 W6 {; s& `( b0 ^. }图7-33  给等高线添加标注
【实例分析】peaks是MATLAB自带函数。
7.3.7  contourc——等高线图形计算
【功能简介】计算等高线矩阵C，用函数contour、contour3和contourf来显示。
# v0 ~; k! i) v, @5 ]3 [【语法格式】
1.C=contourc(Z)
9 o- I0 P2 C: @! {/ f1 a从矩阵Z中计算等高矩阵C，Z可视为XY平面的高度矩阵。等高线的数量和高度值是系统自动确定的。Z至少为2×2矩阵，至少包含两个不同的值，X和Y的范围是1:n与1:m，其中[m,n]=size(Z)。
. D% B% }) `3 ^, B格式变体：
! F3 H6 o  w! [2 a# S6 T5 Z* g6 x        C=contour(Z,n)：确定等高线条数为n，返回等高矩阵。
        C=contour(Z,v)：在高度v处计算等高线。
/ h  V5 V0 x/ _: E" l2.C=contourc(X,Y,Z)
5 t5 x' w' o/ k8 i+ }在X和Y上计算Z的等高矩阵C，Z可视为XY平面的高度矩阵。等高线的数量和高度值是系统自动确定的。Z至少为2×2矩阵，至少包含两个不同的值。
格式变体：
2 i. h& E1 o& j6 V  F6 E: U        C=contour(X,Y,Z,n)：确定等高线条数为n，返回等高矩阵。
        C=contour(X,Y,Z,v)：在高度v处计算等高线。
7 m% x# k7 `; x3 u' A! Q) t【实例7.33】计算peaks函数的等高矩阵。
2 P% E& K4 w6 B4 _>> a=peaks;
+ _# b: h! c$ a>> c=contourc(a,10);                %计算10条等高线的等高矩阵
>> s=size(c)                        %矩阵维度
s =
     2   800
6 F) Q2 |2 {6 x( P6 [>> contour(a,c)                        %画出等高线
  T" v2 D! s( H3 \; R8 G; k6 w执行结果如图7-34所示。


图7-34  计算等高矩阵再绘制等高线
【实例分析】矩阵C是一个2×m矩阵，m的值取决于数据矩阵及等高线的条数。

lili456

7.3.8  fill3——填充三维图
【功能简介】填充三维多边形。
" F; E7 U% y+ v1 i2 Y【语法格式】
9 }0 e0 E$ ~! \& T1.fill3(X,Y,Z,C)
; P+ f+ s# ?/ b参数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)。
0 r* C( e: j* v格式变体：
        fill3(X1,Y1,Z1,C1,X2,Y2,Z2,C2,…)：绘制多个三维多边形。
& {$ ]) X% J( C( w& n6 l4 d  }7 N% A2.fill3(X,Y,Z,ColorSpec)
$ {/ A& E4 |$ L' t0 |用参数ColorSpec指定的颜色对三维多边形进行填充。
【实例7.34】用渐变的颜色来对三维多边形进行着色。
>> X=[0,1,1,2;1,1,2,2;0,0,1,1];
( }. _$ K/ a7 i( Y. P1 S9 ^>> 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];
% X) d0 `6 E+ E5 j! T% q>> C=[.5,1,1,.5;1,.5,.5,.1667;.3333,.3333,.5,.5]        %指定颜色
C =
    0.5000    1.0000    1.0000    0.5000
& ^0 M) j  h! v6 Q# j; @    1.0000    0.5000    0.5000    0.1667
, _& J0 s( X* @- e    0.3333    0.3333    0.5000    0.5000
( C  s& H6 p; C) {; c* F- z>> fill3(X,Y,Z,C)                                                %填充三维多边形
执行结果如图7-35所示。
5 {4 T$ V6 D0 t6 x" a
8 I9 E. H& P$ p% Z图7-35  填充三维多边形
1 o' J/ Z# B$ U, b2 x) {【实例分析】X、Y、Z决定三维多边形的形状，C确定颜色。
7.3.9  sphere——绘制球体
& E% D6 M6 p/ c8 U1 v) O% [1 m【功能简介】绘制球体。
【语法格式】
1.sphere(n)
# C2 j& `0 H8 E8 B3 e在三维直角坐标系中绘制由n×n个面组成的单位球体。
格式变体：
        sphere：默认n=20，绘制单位球体。
7 o; Z0 M& w3 m! L. [; q  i% g2.[X,Y,Z]=sphere(n)
返回三个阶数为(n+1)×(n+1)的坐标矩阵。该命令不画图，只是返回矩阵，要画图可以调用mesh(X,Y,Z)或surf(X,Y,Z)。
【实例7.35】绘制多个球体。
4 r* p" F2 V% k8 Z4 }# c9 }- H( F  |>> [x,y,z]=sphere;
* k& J  n, f1 R7 E# y6 Z1 L>> surf(x,y,z);                %画出第一个球体
>> hold on;
>> surf(x+3,y+2,z);        %画出第二个球体
>> surf(x,y-1,z+2);        %画出第三个球体
>> daspect([1 1 1])
' o, |. r# }5 {' Z  v5 [' ]# B执行结果如图7-36所示。
3 ]- o  [% ?6 x; Y5 T' @$ C9 P
图7-36  绘制多个球体
【实例分析】sphere绘制单位球体，半径为1。
! _4 Q) I% E6 _; Q3 `( @

lili456

7.3.10  contourf——填充二维等高线
【功能简介】填充二维等高线。
【语法格式】
1.contourf(Z)
( k* M* `- j4 {2 w. W+ a画出矩阵Z的二维等高线图，再对等高线之间的区域进行填充，填充的颜色使用当前窗口的颜色映射表。矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
6 j1 N; y- L8 C# u& Y6 |+ W格式变体：
        contourf(Z,n)：画出矩阵Z的n条等高线并对等高线间的区域进行填充。
        contourf(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
2.contourf(X,Y,Z)
画出矩阵Z的等高线图并进行填充，X轴和Y轴的范围由参数X、Y指定。如果X与Y均为矩阵，则两者必须同型且单调递增。
4 |) [' \* o; d0 ], K# v1 j格式变体：
        contourf(X,Y,Z,n)：画出矩阵Z的n条等高线并进行填充。
5 \" ?( a( {2 S2 x' ~4 s& a3 L        contourf(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
【实例7.36】画出等高线并进行填充。
>> a=peaks;
  O8 c9 `6 D2 L% v, c>> contourf(a);                %画出用颜色填充的等高线
6 S% O; Y( [$ H% h9 l& V0 J执行结果如图7-37所示。

图7-37  用颜色填充的等高线图
  A0 x5 g( g" T. H6 \' q' w【实例分析】用颜色对等高线间的区域进行填充，使图形更直观。
5 H3 Q& Z. L; }! R7.3.11  pie3——三维饼图
4 b6 `, e3 T/ _, j/ ?1 N5 p$ b9 W2 w/ b【功能简介】绘制三维饼图。
* H2 l7 H* _3 t! j【语法格式】
7 G' I, `+ b. G5 q1.pie3(X)
& M$ V& G! W5 k' L' i7 p使用X中的数据绘制一个三维饼图，X中的每一个元素都是饼图的一个部分。元素所占扇形面积的大小由X(i)/sum(X)决定，如果sum(X)<1，则只绘制不完整的三维饼形图。
$ x  Q! Y' ]8 x9 K- K& |+ p2.pie3(X,explode)
指定三维饼图中的每一部分是否分离出来。explode是与X同型的数组，用非零值表示相应位置的X中的元素分离出来显示。
【实例7.37】绘制三维饼图。
>> x=[1,3,.7,2.5,2];
& W' c, M2 C, l/ S" C! Y* s>> explode=[0,1,0,0,1];        %设置第二个元素和最后一个元素对应的区域分离出来独立显示
>> pie3(x,explode);                %绘制三维饼图
  z: m/ h- f6 y8 T# B- t>> colormap hsv
! r% b- H0 T" l8 g2 B( J执行结果如图7-38所示。
9 Q' y! C0 P( h5 ]$ Y& \; e
1 {0 ]' U; @0 s5 @图7-38  三维饼图
; k3 z: W3 O0 A5 q$ F( _6 e【实例分析】pie3的调用格式与pie类似。
$ N2 v* h" V6 ^7 |) o( _9 e4 ]

lili456

7.3.12  comet3——三维彗星图
. `+ Z* \; ~: E6 t) d9 \3 K1 _【功能简介】绘制三维彗星图。
/ \9 c5 n# W4 g! ~) [3 [4 ^  L【语法格式】
1 a9 m6 b( P- g: ]- p7 i1.comet3 (z)
显示向量z的三维彗星图，X与Y轴范围由系统自动确定。
2.comet3(x,y,z)
显示由x、y、z确定的三维曲线的彗星图。
  P/ [: Q' K" n4 h. ~3.comet3(x,y,z,p)
指定彗星体的长度为p*length(y)。
" t  f6 H% ^8 N3 ~# i) P5 r# b【实例7.38】绘制一个三维彗星图。
>> t=-10*pi:pi/250:10*pi;
>> comet3((cos(2*t).^2).*sin(t),(sin(2*t).^2).*cos(t),t);                %绘制彗星图
/ B) ?1 B4 a4 H8 {8 U/ u执行结果如图7-39所示。

4 I$ j. W4 h1 f图7-39  三维彗星图
【实例分析】彗星图的绘制是一个动态的过程，用户可按实例中的代码自行运行。
  B. B8 X# D7 t, {% T' G) L7.3.13  cylinder——生成圆柱图形
【功能简介】绘制圆柱图形。
% ]) r  {0 P% U* s. \- u* ~# U1.[X,Y,Z]=cylinder(r,n)
' n; A; O) {" G0 l4 `! T% g返回一个半径为r、高度为1的圆柱体的X、Y、Z轴坐标值，所绘制的圆柱体的圆周有指定的n个距离相同的点。用户可用surf或mesh函数和X、Y、Z参数画出圆柱。
格式变体：
$ @7 c5 T  i5 D1 O/ f" P' `        [X,Y,Z]=cylinder(r)：所绘制的圆柱体的圆周有指定的20个距离相同的点。
        [X,Y,Z]= cylinder：半径采用默认值1。
- w, ]4 v6 [7 h+ N& S4 v% i2.cylinder(…)
: O- C' n6 {6 ]没有输出参量时，直接画出圆柱体。
【实例7.39】绘制一个两头粗、中间细的圆柱。
6 n; }7 F: j  {5 j1 _- U- {; I>> t=0:pi/10:2*pi;
- O9 Y: ?, u8 y. r+ y>> [X,Y,Z]=cylinder(2+cos(t));                %用cylinder函数算出X、Y、Z参数
>> surf(X,Y,Z);                                        %画出圆柱体
>> axis square                                        %调整坐标轴
执行结果如图7-40所示。

图7-40  绘制圆柱体
【实例分析】也可不带返回值，直接用cylinder(2+cos(t))画出圆柱。
8 b7 R# H* I; B# ]# V3 x

lili456

7.3.14  surfc——绘制阴影图及等高线
; J: f5 z8 ^# A- p# J2 D7 @* y" T/ @【功能简介】绘制阴影图及等高线。先画出三维曲面图，再在下方画出二维等高线。
) ~7 {4 C0 i1 H+ k: Y' Z9 C【语法格式】
  C8 T  B, G& _9 r$ k1.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。颜色由高度决定。
* R, J% e6 _6 C1 O, n- _3.surfc(…,C)
图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
  d4 \  {1 U* s7 m0 i& K, ?' Z& @【实例7.40】绘制三维曲面及其二维等高线        。
4 {# P. j8 K! Z* D>> [X,Y]=meshgrid(-8:.5:8);
+ b6 b8 v  H. y( J' t& _( d>> R=sqrt(X.^2+Y.^2)+eps;
# u3 v  X4 S+ p. z>> Z=sin(R)./R;
3 v* Q6 V5 ~9 w4 G>> surfc (X,Y,Z);                        %绘制三维曲面及其二维等高线               
执行结果如图7-41所示。

$ u' x: O5 g0 o: d8 R; S图7-41  绘制三维曲面及其二维等高线
【实例分析】suefc函数调用格式与surf相同，只是多画了一个二维等高线。
' t9 Y7 p/ h5 @, w7.3.15  surfl——带光照模式的曲面图
5 k( ?0 g1 }# D# ^) g2 ?7 E% b6 x【功能简介】绘制带光照模式的曲面图。
【语法格式】
1.surfl(X,Y,Z)
# ~+ v% i) _4 g1 L生成由X、Y和Z定义的曲面图，光照方向和系数采用默认值。
  n0 j2 H# B) _4 w8 {$ v& u* T格式变体：
1 b1 e: R* l! u# r        surfl (Z)：X、Y的值由系统自动生成，等于Z中元素的索引。
$ A+ \! P' p9 [% q- X3 ^2.surfl(…,'light')
用MATLAB光照对象生成一个带光照的曲面。
3 ?! b; B' E* [. }# B! R3.surfc(…,s)
9 j" E! k1 h0 r指定光源的方向s。s是长度为2或3的向量，即s=[azimuth,elevation]或[sx,sy,sz]。默认光源方向是从当前视角开始逆时针45度方向。
. T8 T+ G2 U2 x0 k# I7 @" y4.surfc(…,s,k)
2 o( T% o  U4 X1 Z$ `. C* {指定反射系数k，k是一个常量。k为长度为4的向量[ka,kd,ks,shine]，四个元素分别表示环境光系数、漫反射系数、镜面反射系数和镜面反射亮度，默认值为[0.55,0.6,0.4,10]。
【实例7.41】对peaks函数表示的曲面用surfl函数进行描绘。
1 N5 _6 Y0 n" x$ }( Q' i% G>> [x,y]=meshgrid(-3:1/8:3);
; [- r) a3 W& M) J- U; A* E>> z=peaks(x,y);                %peaks函数
( a) R; G1 x. P" n>> surfl(x,y,z);                        %绘制带光照的曲面图
! j  Y6 Q% A6 @- R7 i  T>> shading interp;
>> colormap gray
& \( k0 s& v7 t执行结果如图7-42所示。

  k* S! o' c2 W' x图7-42  绘制带光照的曲面图
【实例分析】命令colormap gray将色图指定为灰度色图。
  ~7 ?0 W6 s5 J# r. {) Z