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

lili456

7.3  三维图形
, h! m% x0 w, d三维图形表现的是三个变量之间的依赖关系，使用户能够直观看出二元函数数值的变化趋势。下面的小节将介绍三维图形的绘制函数。
* n8 m# a) G# N7.3.1  plot3——绘制三维曲线
【功能简介】绘制三维曲线。
【语法格式】
: \1 r0 ]- @6 C2 H  w! ^1.plot3(X1,Y1,Z1,…)
+ ^4 Z: O  r6 n$ [1 w; ?0 W* QX1、Y1、Z1是同型的向量或矩阵，函数在三维空间中绘制出一条或多条曲线，绘制多条曲线时，曲线的条数等于矩阵的列数。
) L1 W/ }8 ^2 U% \格式变体：
        plot3(X1,Y1,Z1,LineSpec)：用参数LineSpec指定画线的线型、标记符号和颜色。
) w( _6 J! g% m' R2.plot3(…,'PropertyName','PropertyValue')
设置所绘图形指定属性的属性值。
) Z3 q: w0 |1 J$ D4 Y1 E! D【实例7.27】绘制三维螺旋图。
0 p) s! y1 `" E. l6 b3 V7 ^>> t=0:pi/50:10*pi;
>> plot3(sin(t),cos(t),t);        %三维螺旋图
8 E7 m9 F( R" o4 J) h' v+ P! Z3 y>> grid on;
>> axis square;
执行结果如图7-28所示。
( H7 o0 i" d: ?; \' F7 c
图7-28  三维螺旋图
/ M3 ?: e; h+ F2 z* K  Q【实例分析】三维螺旋图的参数方程形式为x=sin(t),y=cos(t),z=t。

lili456

7.3.2  mesh——绘制三维网格图
【功能简介】绘制三维网格图。
【语法格式】
1.mesh(X,Y,Z)
' L1 M  C/ u& a, k6 S; F生成由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定义， 颜色与高度成比例。
# t' D( q3 @6 M: C2.mesh(Z)
X与Y自动生成。[m,n]=size(Z)，则X=1:n，Y=1:m。颜色由高度决定。
3.mesh(…,C)
图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.28】绘制三维网格图。
. g( Z0 q  G6 V/ q* @$ y>> [X,Y]=meshgrid(-8:.5:8);        %构造X、Y矩阵
4 `  w8 s3 g' q, y" G>> R=sqrt(X.^2+Y.^2)+eps;
>> Z=sin(R)./R;
>> mesh(X,Y,Z);                        %绘制三维网格图
' @: H" G9 m5 {, n* V执行结果如图7-29所示。
$ N' n2 v2 P. J  o( l7 E

图7-29  三维网格图
【实例分析】不指定颜色时，网格颜色由高度决定。

lili456

7.3.3  surf——三维曲面图
+ f% i. y0 l% u9 H. l9 i- u【功能简介】绘制三维曲面图。
【语法格式】
- G# t: ^1 F& S. `; D1.surf(X,Y,Z)
+ A/ T- k7 S* [生成由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)
6 L( W6 c& u5 J4 S7 A8 t0 O3 S1 k( e" I5 D图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.29】绘制三维曲面图。
>> [X,Y]=meshgrid(-8:.5:8);        %生成X、Y矩阵
>> R=sqrt(X.^2+Y.^2)+eps;
) Q/ i9 J3 k/ Y+ l$ O2 b# t>> Z=sin(R)./R;
; a) L4 ]- i$ q' Y4 n>> surf(X,Y,Z);                                %绘制三维曲面图
执行结果如图7-30所示。

& S/ a" w+ g+ B6 }图7-30  三维曲面图
【实例分析】surf函数与mesh函数调用格式是相同的。

lili456

7.3.4  contour3——三维等高线绘制
& _% _$ j) p* E4 H5 P【功能简介】绘制三维等高线。
【语法格式】
1.contour3(Z)
画出矩阵Z的等高线图，矩阵Z至少为2×2大小，系统按矩阵Z中元素距离XY平面的高度，自动选择等高线的条数进行绘制。X和Y的值是自动确定的，[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
- e, P% k. j: o* s: n格式变体：
        contour3(Z,n)：画出矩阵Z的n条等高线。
9 s  V8 k" a0 c, f  C9 u        contour3(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
2.contour3(X,Y,Z)
9 R, G0 r% T  x) V8 r$ i) G画出矩阵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)。
" T7 E+ z- n" s) v% h: J3.contour3(…,LineSpec)
参数LineSpec指定线型、标记符号和颜色。
, {4 C! J! V8 n" j9 Y, B& Y/ K( q【实例7.30】绘制等高线。
6 B  _5 T8 v0 L>> [X,Y]=meshgrid([-2:.25:2]);
+ k8 [& J& x1 g8 a" N- q>> Z=X.*exp(-X.^2-Y.^2);
0 I* D) M! d3 [1 G% l: W>> contour3(X,Y,Z,30);        %绘制等高线
$ e& q; J7 k+ M4 ]4 z>> colormap cool;
; `& W) {" U6 S" u+ v执行结果如图7-31所示。

图7-31  contour3函数绘制等高线
【实例分析】等高线绘制有contour3和contour两个函数均可完成，contour函数绘制的是二维的等高线。
7 q. ^5 E, s) M4 P

lili456

7.3.5  contour——曲面的等高线
【功能简介】绘制曲面的等高线。
【语法格式】
: D( N* P$ b$ A- |1.contour(Z)
/ \; R* W8 t% o4 T! Z! G' n' c% R  d画出矩阵Z的二维等高线图，矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
格式变体：
3 V; i- h4 U8 V8 N        contour(Z,n)：画出矩阵Z的n条等高线。
- ]1 M. f# I. ?" ~! X1 {% ]$ F/ s        contour(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
5 q4 J: S; A, v5 j* n8 j" C. N2.contour(X,Y,Z)
. D4 F( g0 \" j" o0 O" v: Z! X画出矩阵Z的二维等高线图，X轴和Y轴的范围由参数X、Y指定。如果X与Y均为矩阵，则两者必须同型且单调递增。
8 c& G' d& @: \0 J1 z格式变体：
8 J9 T6 ^" u3 K0 b8 {9 d1 J        contour(X,Y,Z,n)：画出矩阵Z的n条等高线。
        contour(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
7 v9 u' u; U8 I6 e) S" L3.contour(…,LineSpec)
3 P# D6 b' }3 m7 e参数LineSpec指定线型、标记符号和颜色。
1 A  V0 O1 C7 s6 v/ b! q9 k% y【实例7.31】画出曲面的等高线图。
>> [X,Y]=meshgrid([-2:.25:2]);
* V# v7 z3 o$ _! L# F>> Z=X.*exp(-X.^2-Y.^2);
( C8 Q" [- w* I) `% P1 @; X  O>> contour(X,Y,Z,30);        %画出曲面的二维等高线
, R  J( s. G0 x* Z* W执行结果如图7-32所示。
' L7 ~2 L# I. s, I/ I4 @
7 V" l- H* [  n0 P( R$ x: Y 图7-32  曲面的等高线
5 k# [4 C0 N0 \; r5 u/ i* `; M【实例分析】contour函数画出的是二维图形，contour3函数画出的是三维图形。

lili456

7.3.6  clabel——等高线高度标签
【功能简介】在二维等高线图中绘制等高线高度标签。
【语法格式】
1.clabel(C,h)
( `, \3 ?7 Z# I7 x, ]0 d; {C为标签矩阵，h为等高线句柄。函数旋转标签到适当角度，在空间允许的情况下插入标签。
- c7 e' [- P* X/ F格式变体：
        clabel(C,h,v)：在高度v处插入标签C。
        clabel(C,h,'manual')：手动添加标签到句柄h指定的等高线图中，用鼠标左键或空格键在最接近位置上放置标签，用回车键结束操作。
2.clabel(C)
添加标签到当前等高线图中，随机选择标签位置。
2 Z( F$ H5 x; T0 a& k& F2 `格式变体：
2 J' ^  Z; ?# L: X( D+ u        clabel(C,v)：在v指定的高度处添加标签。
' o% J- h8 V9 d& s        clabel(C,'manual')：手动添加标签到当前等高线图中。
【实例7.32】给等高线做标注。
+ B" y$ v" i. j. l6 ?1 d6 \% x>> [x,y,z]=peaks;
7 c3 N9 w) w3 d) D5 P: [>> C=contour(x,y,z);        %绘制等高线
+ f& t! p# l* c$ B0 [6 _4 b) ?; H>> clabel(C);                %给等高线添加标注
执行结果如图7-33所示。
" \* p/ C5 `0 g) @/ S3 W$ L- O! b
图7-33  给等高线添加标注
7 A# D: i8 N+ q1 x& T9 C" h【实例分析】peaks是MATLAB自带函数。
7.3.7  contourc——等高线图形计算
$ a* u5 e1 O% |' Q2 f6 I【功能简介】计算等高线矩阵C，用函数contour、contour3和contourf来显示。
【语法格式】
1.C=contourc(Z)
+ g( I) a% g/ i1 T6 \! ~+ R$ E$ S  Q从矩阵Z中计算等高矩阵C，Z可视为XY平面的高度矩阵。等高线的数量和高度值是系统自动确定的。Z至少为2×2矩阵，至少包含两个不同的值，X和Y的范围是1:n与1:m，其中[m,n]=size(Z)。
$ h& f0 z: d8 @# M5 W7 u格式变体：
3 }" w7 m' D- @1 N0 T% I        C=contour(Z,n)：确定等高线条数为n，返回等高矩阵。
& ?# ?; R9 |5 z        C=contour(Z,v)：在高度v处计算等高线。
2.C=contourc(X,Y,Z)
; Y+ Z1 P, E; ^, _$ e在X和Y上计算Z的等高矩阵C，Z可视为XY平面的高度矩阵。等高线的数量和高度值是系统自动确定的。Z至少为2×2矩阵，至少包含两个不同的值。
! P8 [- f: Y4 T& E" L格式变体：
$ M5 W- W# e" |2 Z* }7 @. r5 V4 I        C=contour(X,Y,Z,n)：确定等高线条数为n，返回等高矩阵。
, M6 |: X) i) n3 K! `0 [" i        C=contour(X,Y,Z,v)：在高度v处计算等高线。
【实例7.33】计算peaks函数的等高矩阵。
9 u' N0 l& |2 [4 V>> a=peaks;
>> c=contourc(a,10);                %计算10条等高线的等高矩阵
6 T, {$ `( s2 s& D# \# V>> s=size(c)                        %矩阵维度
s =
     2   800
" u4 r5 _  l# }4 {>> contour(a,c)                        %画出等高线
执行结果如图7-34所示。
; w" w0 D# W, h* I3 g
3 _5 x! W& \, x$ w
% A7 z7 q: c/ Z4 e6 N1 p3 S图7-34  计算等高矩阵再绘制等高线
【实例分析】矩阵C是一个2×m矩阵，m的值取决于数据矩阵及等高线的条数。

lili456

7.3.8  fill3——填充三维图
# {0 B5 l' R" F+ \9 \【功能简介】填充三维多边形。
【语法格式】
; c9 `  h! o0 X. \4 N1.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)。
格式变体：
6 L! x# }( g- @1 s9 r8 q9 T1 p4 r3 T        fill3(X1,Y1,Z1,C1,X2,Y2,Z2,C2,…)：绘制多个三维多边形。
2.fill3(X,Y,Z,ColorSpec)
用参数ColorSpec指定的颜色对三维多边形进行填充。
' l" q% K: J9 Y【实例7.34】用渐变的颜色来对三维多边形进行着色。
>> X=[0,1,1,2;1,1,2,2;0,0,1,1];
6 s/ z6 Z, n. t7 @% o3 c& x>> 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];
: p, o5 o- }5 _8 P# q>> 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
) r3 u! X& X3 u6 K  _4 K# v    0.3333    0.3333    0.5000    0.5000
( [3 E" G6 u, o+ v>> fill3(X,Y,Z,C)                                                %填充三维多边形
" C+ w* b2 A! {执行结果如图7-35所示。

/ c# h( U. J9 |# V* u/ e图7-35  填充三维多边形
【实例分析】X、Y、Z决定三维多边形的形状，C确定颜色。
7.3.9  sphere——绘制球体
7 l+ p: C; H* `7 O% g- k4 ^【功能简介】绘制球体。
+ ^1 N+ D# `5 M3 x+ }. ^- Z  i【语法格式】
$ j9 }6 I% W9 b2 T1 W1.sphere(n)
在三维直角坐标系中绘制由n×n个面组成的单位球体。
格式变体：
        sphere：默认n=20，绘制单位球体。
: u/ M7 h, o! g. h  e2.[X,Y,Z]=sphere(n)
返回三个阶数为(n+1)×(n+1)的坐标矩阵。该命令不画图，只是返回矩阵，要画图可以调用mesh(X,Y,Z)或surf(X,Y,Z)。
" p: q/ k+ K9 e( G$ R* k* w. ~【实例7.35】绘制多个球体。
>> [x,y,z]=sphere;
>> surf(x,y,z);                %画出第一个球体
>> hold on;
>> surf(x+3,y+2,z);        %画出第二个球体
: \4 O+ T' p7 W- f7 i( w5 l>> surf(x,y-1,z+2);        %画出第三个球体
/ o* @) w" w4 u% V, q- m>> daspect([1 1 1])
执行结果如图7-36所示。
: u4 I0 b$ }6 U/ j$ I
- H  k  t7 d# v图7-36  绘制多个球体
【实例分析】sphere绘制单位球体，半径为1。

lili456

7.3.10  contourf——填充二维等高线
【功能简介】填充二维等高线。
【语法格式】
/ M  x0 j/ ^# g& _1.contourf(Z)
! o$ x, Z3 H* G& W5 C* w画出矩阵Z的二维等高线图，再对等高线之间的区域进行填充，填充的颜色使用当前窗口的颜色映射表。矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
格式变体：
, d% s( W7 Y( o        contourf(Z,n)：画出矩阵Z的n条等高线并对等高线间的区域进行填充。
        contourf(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
5 c+ E4 d4 t1 t3 V% ?2.contourf(X,Y,Z)
画出矩阵Z的等高线图并进行填充，X轴和Y轴的范围由参数X、Y指定。如果X与Y均为矩阵，则两者必须同型且单调递增。
  {7 y. d. q, e$ C格式变体：
        contourf(X,Y,Z,n)：画出矩阵Z的n条等高线并进行填充。
6 o+ P6 H! G. r  D, Z, X        contourf(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
【实例7.36】画出等高线并进行填充。
4 m* s8 f/ L. q  p# L# v>> a=peaks;
>> contourf(a);                %画出用颜色填充的等高线
3 |/ [( Q; @! H! s5 T- x5 L- K执行结果如图7-37所示。

图7-37  用颜色填充的等高线图
5 ^; y" M% N9 o0 [, q【实例分析】用颜色对等高线间的区域进行填充，使图形更直观。
  @% t. \, [1 H2 }0 c5 q: A; f. C7.3.11  pie3——三维饼图
1 t- \: f1 J  p, |0 n7 D$ D【功能简介】绘制三维饼图。
# s; ?6 M+ o* [  g* r7 }【语法格式】
1.pie3(X)
使用X中的数据绘制一个三维饼图，X中的每一个元素都是饼图的一个部分。元素所占扇形面积的大小由X(i)/sum(X)决定，如果sum(X)<1，则只绘制不完整的三维饼形图。
3 f: x1 S3 f. t' ?) Y2.pie3(X,explode)
指定三维饼图中的每一部分是否分离出来。explode是与X同型的数组，用非零值表示相应位置的X中的元素分离出来显示。
/ ^- o9 k0 }' W2 u- h3 |【实例7.37】绘制三维饼图。
>> x=[1,3,.7,2.5,2];
>> explode=[0,1,0,0,1];        %设置第二个元素和最后一个元素对应的区域分离出来独立显示
4 g( w* b2 O" C>> pie3(x,explode);                %绘制三维饼图
>> colormap hsv
执行结果如图7-38所示。
& \/ i- ^$ L+ D" c
图7-38  三维饼图
【实例分析】pie3的调用格式与pie类似。

lili456

7.3.12  comet3——三维彗星图
【功能简介】绘制三维彗星图。
, R& p) [: V" p) A( w& ?* T1 n【语法格式】
1.comet3 (z)
显示向量z的三维彗星图，X与Y轴范围由系统自动确定。
9 w8 T1 r2 e; z0 ~9 a/ }2.comet3(x,y,z)
" R. A# w2 O0 s" X% |1 b2 {显示由x、y、z确定的三维曲线的彗星图。
7 ], J# M0 R( E! c. a) F) g  J3.comet3(x,y,z,p)
* A( o7 w2 v5 Q. R4 s/ c0 ^" W) @指定彗星体的长度为p*length(y)。
- X1 u6 N$ u8 A. `【实例7.38】绘制一个三维彗星图。
+ p" O+ [  t6 g3 B>> t=-10*pi:pi/250:10*pi;
$ r" Z5 O, j  h4 J- x/ m>> comet3((cos(2*t).^2).*sin(t),(sin(2*t).^2).*cos(t),t);                %绘制彗星图
执行结果如图7-39所示。

图7-39  三维彗星图
. b, b0 V+ E  Z2 P' k" z【实例分析】彗星图的绘制是一个动态的过程，用户可按实例中的代码自行运行。
7.3.13  cylinder——生成圆柱图形
5 t, G1 Y" ?; X) Q8 J: W( h: m【功能简介】绘制圆柱图形。
- }, [/ Q8 R8 z0 ^) x) m1.[X,Y,Z]=cylinder(r,n)
% R. s" ]5 Z- O7 u" k7 m返回一个半径为r、高度为1的圆柱体的X、Y、Z轴坐标值，所绘制的圆柱体的圆周有指定的n个距离相同的点。用户可用surf或mesh函数和X、Y、Z参数画出圆柱。
格式变体：
        [X,Y,Z]=cylinder(r)：所绘制的圆柱体的圆周有指定的20个距离相同的点。
% b  e# o3 V1 h& I: g        [X,Y,Z]= cylinder：半径采用默认值1。
2.cylinder(…)
没有输出参量时，直接画出圆柱体。
【实例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                                        %调整坐标轴
执行结果如图7-40所示。

8 k8 p! U# ?( s$ S图7-40  绘制圆柱体
/ N6 u: t/ t: ^/ j1 `1 V$ v# [【实例分析】也可不带返回值，直接用cylinder(2+cos(t))画出圆柱。
+ x' a, C1 F) L2 S

lili456

7.3.14  surfc——绘制阴影图及等高线
【功能简介】绘制阴影图及等高线。先画出三维曲面图，再在下方画出二维等高线。
, w2 q# E: x  c! O【语法格式】
1.surfc(X,Y,Z)
  M% ]4 x6 }) Q& }& J  [+ B; Y6 X生成由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函数一致。
) d8 o; ?  N; X4 G) z4 q2.surfc(Z)
9 i( c" u( s5 i( hX与Y自动生成。[m,n]=size(Z)，则X=1:n，Y=1:m。颜色由高度决定。
9 B  A" [6 o* N" S$ E6 F5 h3.surfc(…,C)
图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.40】绘制三维曲面及其二维等高线        。
3 j8 d6 F  k7 [# m* l>> [X,Y]=meshgrid(-8:.5:8);
! h5 S( f9 }; L6 i>> R=sqrt(X.^2+Y.^2)+eps;
, N7 p/ I2 t$ a; J; \3 B6 t# y>> Z=sin(R)./R;
>> surfc (X,Y,Z);                        %绘制三维曲面及其二维等高线               
执行结果如图7-41所示。
; l. c' O( I% ^6 c/ o
图7-41  绘制三维曲面及其二维等高线
【实例分析】suefc函数调用格式与surf相同，只是多画了一个二维等高线。
7.3.15  surfl——带光照模式的曲面图
5 j3 Z- K2 M- l- E# h* D【功能简介】绘制带光照模式的曲面图。
- d' S1 Z: R5 h  K3 r3 y( |2 v* X【语法格式】
  H0 d2 `) o  P+ I* @$ q2 w: |6 q1.surfl(X,Y,Z)
生成由X、Y和Z定义的曲面图，光照方向和系数采用默认值。
3 j8 i' c% I. f$ T格式变体：
, r& r9 K( f5 n. i& j, k        surfl (Z)：X、Y的值由系统自动生成，等于Z中元素的索引。
4 G+ S3 q8 i5 G& a$ H( P1 n4 l2.surfl(…,'light')
/ a2 r0 _) p( W7 A1 z) {& r用MATLAB光照对象生成一个带光照的曲面。
3.surfc(…,s)
指定光源的方向s。s是长度为2或3的向量，即s=[azimuth,elevation]或[sx,sy,sz]。默认光源方向是从当前视角开始逆时针45度方向。
, a9 p$ W" B! W8 c- T9 ~3 a$ u1 }, a4.surfc(…,s,k)
指定反射系数k，k是一个常量。k为长度为4的向量[ka,kd,ks,shine]，四个元素分别表示环境光系数、漫反射系数、镜面反射系数和镜面反射亮度，默认值为[0.55,0.6,0.4,10]。
4 U( T/ ^, K6 {& x1 N【实例7.41】对peaks函数表示的曲面用surfl函数进行描绘。
% V6 K& d* s5 f% ^>> [x,y]=meshgrid(-3:1/8:3);
>> z=peaks(x,y);                %peaks函数
) c* ?* E3 q6 V. }6 s/ i) v>> surfl(x,y,z);                        %绘制带光照的曲面图
>> shading interp;
% u1 ~0 M- @* K1 `7 }7 a>> colormap gray
0 ]% U. F( [6 R4 {1 H) C6 d  t" s执行结果如图7-42所示。
/ I/ o  i/ ^" T& ], p' l- `7 U
( U( i" u' i3 F% W图7-42  绘制带光照的曲面图
【实例分析】命令colormap gray将色图指定为灰度色图。