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

lili456

7.3  三维图形
三维图形表现的是三个变量之间的依赖关系，使用户能够直观看出二元函数数值的变化趋势。下面的小节将介绍三维图形的绘制函数。
7.3.1  plot3——绘制三维曲线
0 b. q. s2 d3 d2 M【功能简介】绘制三维曲线。
【语法格式】
( j. r9 i: d0 T0 v1 x1.plot3(X1,Y1,Z1,…)
. n; x% A: A: Q4 L- N3 }' TX1、Y1、Z1是同型的向量或矩阵，函数在三维空间中绘制出一条或多条曲线，绘制多条曲线时，曲线的条数等于矩阵的列数。
格式变体：
) t) Q- Q& W6 A( s) ~5 ?        plot3(X1,Y1,Z1,LineSpec)：用参数LineSpec指定画线的线型、标记符号和颜色。
2.plot3(…,'PropertyName','PropertyValue')
: \' x7 G: v3 L0 H) H3 i设置所绘图形指定属性的属性值。
【实例7.27】绘制三维螺旋图。
: A0 M7 j4 S1 S1 [% S/ X>> t=0:pi/50:10*pi;
! V6 Z% d% K. ~3 `' j>> plot3(sin(t),cos(t),t);        %三维螺旋图
4 Q* |: h+ z5 X( h5 F>> grid on;
>> axis square;
: [0 c6 i4 I. P5 k" [" ~4 K. u4 n执行结果如图7-28所示。

图7-28  三维螺旋图
【实例分析】三维螺旋图的参数方程形式为x=sin(t),y=cos(t),z=t。

lili456

7.3.2  mesh——绘制三维网格图
% j' h1 X% Y8 q+ g+ J【功能简介】绘制三维网格图。
【语法格式】
2 y1 E1 d; p7 X3 A( I0 N+ Y1.mesh(X,Y,Z)
9 k/ l( ]3 R% O+ o: J* I生成由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。颜色由高度决定。
3 ^; [6 [3 ^5 A/ s, u3.mesh(…,C)
' U: }2 x: b) P图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.28】绘制三维网格图。
5 \9 W$ L- _8 @1 s>> [X,Y]=meshgrid(-8:.5:8);        %构造X、Y矩阵
>> R=sqrt(X.^2+Y.^2)+eps;
>> Z=sin(R)./R;
7 I4 g6 f" l1 ]0 B4 E>> mesh(X,Y,Z);                        %绘制三维网格图
执行结果如图7-29所示。

2 v4 z1 Q3 J3 b$ l4 O9 r
5 {/ b. A4 @: a: A: O7 r5 y图7-29  三维网格图
【实例分析】不指定颜色时，网格颜色由高度决定。
2 M" X( X# D  o) G1 g+ H

lili456

7.3.3  surf——三维曲面图
【功能简介】绘制三维曲面图。
) Z. B) x7 J8 w: P" J9 ^3 F- \4 s【语法格式】
1.surf(X,Y,Z)
+ z) |  z, F8 A1 _/ ?生成由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函数一致。
& u4 l: Y! y7 d3 l2.surf (Z)
X与Y自动生成。[m,n]=size(Z)，则X=1:n，Y=1:m。颜色由高度决定。
3.surf (…,C)
图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.29】绘制三维曲面图。
>> [X,Y]=meshgrid(-8:.5:8);        %生成X、Y矩阵
>> R=sqrt(X.^2+Y.^2)+eps;
  Z. r6 A' ]3 g6 W3 P>> Z=sin(R)./R;
>> surf(X,Y,Z);                                %绘制三维曲面图
执行结果如图7-30所示。
, z2 R" Y" V, U+ y" v  s& y
; b& y4 S3 S- C& Q图7-30  三维曲面图
【实例分析】surf函数与mesh函数调用格式是相同的。
) S$ X/ p: ~5 x. W

lili456

7.3.4  contour3——三维等高线绘制
【功能简介】绘制三维等高线。
; P% H$ E( O3 q【语法格式】
+ N5 w$ N& M8 l  L1.contour3(Z)
画出矩阵Z的等高线图，矩阵Z至少为2×2大小，系统按矩阵Z中元素距离XY平面的高度，自动选择等高线的条数进行绘制。X和Y的值是自动确定的，[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
格式变体：
        contour3(Z,n)：画出矩阵Z的n条等高线。
8 b3 \$ n" [2 A        contour3(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
5 d2 y0 s: Q; S+ A. i# C( J3 o2.contour3(X,Y,Z)
画出矩阵Z的等高线图，X轴和Y轴的范围由参数X、Y指定。如果X为矩阵，则使用X(1,来定义X轴的范围，如果Y为矩阵，则使用Y(:,1)来定义Y轴的范围。如果X与Y均为矩阵，则两者必须同型。
! T7 L1 ^6 y' I! b格式变体：
+ F! [3 r7 E. E1 @* J! C        contour3(X,Y,Z,n)：画出矩阵Z的n条等高线。
        contour3(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
( @( e% m4 p. M7 H5 W3.contour3(…,LineSpec)
( _+ O* D. p9 i# s  p/ f: K4 u& G6 p参数LineSpec指定线型、标记符号和颜色。
8 X" _: p- d6 v0 x. y! K! z8 V* f【实例7.30】绘制等高线。
>> [X,Y]=meshgrid([-2:.25:2]);
>> Z=X.*exp(-X.^2-Y.^2);
- N* b3 @- f+ ~2 |; _3 K7 U+ l>> contour3(X,Y,Z,30);        %绘制等高线
# S- F# Z+ w, ?& p>> colormap cool;
执行结果如图7-31所示。
( E( w8 ~1 y! S8 J$ G
% D' [3 \4 A& ]6 u+ n2 w/ b 图7-31  contour3函数绘制等高线
【实例分析】等高线绘制有contour3和contour两个函数均可完成，contour函数绘制的是二维的等高线。
) _" M: F* B+ `: z) x

lili456

7.3.5  contour——曲面的等高线
! A% R+ j/ j) T5 J. r6 z1 `【功能简介】绘制曲面的等高线。
' C% T# u, k- ^! U【语法格式】
& S. k8 d3 X% N8 P4 |% c1.contour(Z)
+ \$ F9 h: r2 }2 A! F$ e; C4 Y画出矩阵Z的二维等高线图，矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
格式变体：
1 a5 P9 h: W9 J        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条等高线。
: S4 p. A% D5 O2 I" n- x6 Y0 f        contour(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
3.contour(…,LineSpec)
7 E) A5 C( H8 c2 q: K& s参数LineSpec指定线型、标记符号和颜色。
: G+ Q- K" L: |# Y( J8 J) j8 N; x. k【实例7.31】画出曲面的等高线图。
/ J. Y& T4 o/ a$ b2 N>> [X,Y]=meshgrid([-2:.25:2]);
9 M, F) i8 E5 t% Y: `>> Z=X.*exp(-X.^2-Y.^2);
>> contour(X,Y,Z,30);        %画出曲面的二维等高线
执行结果如图7-32所示。
5 D( \* A+ T$ Q$ G6 F* B7 \7 n6 E
1 D7 {: k$ J4 t; d; f9 I2 K% B 图7-32  曲面的等高线
【实例分析】contour函数画出的是二维图形，contour3函数画出的是三维图形。

lili456

7.3.6  clabel——等高线高度标签
( A" k2 Q+ o& p2 j5 R【功能简介】在二维等高线图中绘制等高线高度标签。
【语法格式】
3 `  O' }' L, ^& u2 H9 `7 `3 ~4 Y1.clabel(C,h)
C为标签矩阵，h为等高线句柄。函数旋转标签到适当角度，在空间允许的情况下插入标签。
格式变体：
        clabel(C,h,v)：在高度v处插入标签C。
        clabel(C,h,'manual')：手动添加标签到句柄h指定的等高线图中，用鼠标左键或空格键在最接近位置上放置标签，用回车键结束操作。
7 F6 B7 A, d2 d; j1 s% e) U2.clabel(C)
添加标签到当前等高线图中，随机选择标签位置。
格式变体：
$ H- _  [7 W- G* v! ^, x9 s. p        clabel(C,v)：在v指定的高度处添加标签。
3 ~) H2 F) f9 Y$ K' \        clabel(C,'manual')：手动添加标签到当前等高线图中。
【实例7.32】给等高线做标注。
. n0 [$ j% D  K7 Q>> [x,y,z]=peaks;
$ q' U1 v! p: d  ?1 S1 ^9 o1 D, Q" @>> C=contour(x,y,z);        %绘制等高线
3 B' _1 B# ^1 E>> clabel(C);                %给等高线添加标注
执行结果如图7-33所示。
- E( z- u' T& y0 C" i! H
图7-33  给等高线添加标注
【实例分析】peaks是MATLAB自带函数。
" |, |7 b8 i! ]! Q& X: [! K7.3.7  contourc——等高线图形计算
【功能简介】计算等高线矩阵C，用函数contour、contour3和contourf来显示。
# j% x" m9 ~# H+ l# B% E【语法格式】
2 N' ?$ g- Q) _7 u1.C=contourc(Z)
从矩阵Z中计算等高矩阵C，Z可视为XY平面的高度矩阵。等高线的数量和高度值是系统自动确定的。Z至少为2×2矩阵，至少包含两个不同的值，X和Y的范围是1:n与1:m，其中[m,n]=size(Z)。
; g, X" v  U0 b. n$ Q格式变体：
        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矩阵，至少包含两个不同的值。
格式变体：
        C=contour(X,Y,Z,n)：确定等高线条数为n，返回等高矩阵。
) w& I& M2 l, b0 `8 [        C=contour(X,Y,Z,v)：在高度v处计算等高线。
【实例7.33】计算peaks函数的等高矩阵。
8 e/ E7 H0 |# C' j! ~8 j: k. a" N$ H>> a=peaks;
>> c=contourc(a,10);                %计算10条等高线的等高矩阵
>> s=size(c)                        %矩阵维度
s =
* {) v6 d3 d3 ~7 G     2   800
& ]: h  O4 t6 U% H! N: V5 l$ ?>> contour(a,c)                        %画出等高线
执行结果如图7-34所示。


图7-34  计算等高矩阵再绘制等高线
) X% t0 p# r! X1 r2 V( W【实例分析】矩阵C是一个2×m矩阵，m的值取决于数据矩阵及等高线的条数。
  j7 l  z" N! {3 R8 X, ]2 \+ v

lili456

7.3.8  fill3——填充三维图
【功能简介】填充三维多边形。
【语法格式】
, y. R$ w4 I* h; p0 ?; ?# W  m1 T1.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,…)：绘制多个三维多边形。
2.fill3(X,Y,Z,ColorSpec)
用参数ColorSpec指定的颜色对三维多边形进行填充。
【实例7.34】用渐变的颜色来对三维多边形进行着色。
>> X=[0,1,1,2;1,1,2,2;0,0,1,1];
8 N1 ?  z' x3 H3 y$ N>> Y=[1,1,1,1;1,0,1,0;0,0,0,0];
5 f# V% S  _" d8 J>> Z=[1,1,1,1;1,0,1,0;0,0,0,0];
% K' N- _, U, d8 }$ p>> C=[.5,1,1,.5;1,.5,.5,.1667;.3333,.3333,.5,.5]        %指定颜色
+ k7 c3 s! M" {( C) o8 Q7 U' kC =
- u, P/ {. Y6 J% v) l& [    0.5000    1.0000    1.0000    0.5000
) E# P' X+ R7 L( ~. z. o0 |    1.0000    0.5000    0.5000    0.1667
    0.3333    0.3333    0.5000    0.5000
! V/ z5 W6 i1 k! R2 @, v>> fill3(X,Y,Z,C)                                                %填充三维多边形
3 ~2 _. m5 b" w3 p执行结果如图7-35所示。
4 h/ G/ S6 c. d& H& n# w
图7-35  填充三维多边形
【实例分析】X、Y、Z决定三维多边形的形状，C确定颜色。
6 o. ^* v: z) b6 \8 m; n/ \" @' c7.3.9  sphere——绘制球体
3 Y* R, Z$ Q; a: z; d" O& j" s( z9 \3 Y【功能简介】绘制球体。
# a/ Y% M6 N% _7 ^6 Z2 q+ D【语法格式】
. y( e& ~- |, b3 F1.sphere(n)
在三维直角坐标系中绘制由n×n个面组成的单位球体。
! N9 h9 s% J: B% y5 m4 V/ J格式变体：
" `3 W; ]/ P* W        sphere：默认n=20，绘制单位球体。
3 b$ y, \/ a- {( R2.[X,Y,Z]=sphere(n)
返回三个阶数为(n+1)×(n+1)的坐标矩阵。该命令不画图，只是返回矩阵，要画图可以调用mesh(X,Y,Z)或surf(X,Y,Z)。
1 k1 Y' l4 ]4 l【实例7.35】绘制多个球体。
. W( U( p4 Y$ a>> [x,y,z]=sphere;
4 Q  z3 w: ]8 e2 {$ ~>> surf(x,y,z);                %画出第一个球体
( l. X" S, R, @& m4 m>> hold on;
>> surf(x+3,y+2,z);        %画出第二个球体
/ T7 ^1 R1 Q; E2 ^$ Q>> surf(x,y-1,z+2);        %画出第三个球体
8 k, ~5 P" y6 w( o7 U4 d( C>> daspect([1 1 1])
执行结果如图7-36所示。

. H0 s: @4 D( G6 h0 J6 n图7-36  绘制多个球体
【实例分析】sphere绘制单位球体，半径为1。

lili456

7.3.10  contourf——填充二维等高线
, A: G) n; v! X【功能简介】填充二维等高线。
: J# K: Q, l1 q" p' {; T【语法格式】
1.contourf(Z)
. Z* \6 P# h  s1 v( g8 J# ]. l) M  V画出矩阵Z的二维等高线图，再对等高线之间的区域进行填充，填充的颜色使用当前窗口的颜色映射表。矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
% l4 k. \4 g( [$ F格式变体：
1 ~7 [" f4 v8 a2 v5 r/ n        contourf(Z,n)：画出矩阵Z的n条等高线并对等高线间的区域进行填充。
        contourf(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
4 H) h% [+ }6 Q7 Z; Q( n8 D& [2.contourf(X,Y,Z)
画出矩阵Z的等高线图并进行填充，X轴和Y轴的范围由参数X、Y指定。如果X与Y均为矩阵，则两者必须同型且单调递增。
) ]) P) P* y7 g$ O/ E. v( u* w, Y格式变体：
! I" c* b* V( G. M* G; K. Z        contourf(X,Y,Z,n)：画出矩阵Z的n条等高线并进行填充。
        contourf(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
# d. S3 w# i8 o3 I【实例7.36】画出等高线并进行填充。
>> a=peaks;
>> contourf(a);                %画出用颜色填充的等高线
执行结果如图7-37所示。

; h) H- I% @: w% Q/ p图7-37  用颜色填充的等高线图
0 k! x  m% X! A4 g& e8 p! z【实例分析】用颜色对等高线间的区域进行填充，使图形更直观。
7.3.11  pie3——三维饼图
( {4 f- J8 J- w8 V) i6 |" q; q【功能简介】绘制三维饼图。
9 [( x- ?& G( g# {0 e【语法格式】
" T+ }3 Y5 i7 R. L, K1.pie3(X)
使用X中的数据绘制一个三维饼图，X中的每一个元素都是饼图的一个部分。元素所占扇形面积的大小由X(i)/sum(X)决定，如果sum(X)<1，则只绘制不完整的三维饼形图。
2.pie3(X,explode)
指定三维饼图中的每一部分是否分离出来。explode是与X同型的数组，用非零值表示相应位置的X中的元素分离出来显示。
【实例7.37】绘制三维饼图。
>> x=[1,3,.7,2.5,2];
. j1 Z+ ^5 D! ^) l  k8 T>> explode=[0,1,0,0,1];        %设置第二个元素和最后一个元素对应的区域分离出来独立显示
( Z5 b, {9 p* D2 V>> pie3(x,explode);                %绘制三维饼图
>> colormap hsv
执行结果如图7-38所示。
" e# Q. ~5 `1 ~1 e' U7 E% G  O
3 a6 @2 W1 \9 Y9 v2 J$ q" ^  j图7-38  三维饼图
8 d% b+ a/ M3 ^8 R% q' v【实例分析】pie3的调用格式与pie类似。

lili456

7.3.12  comet3——三维彗星图
【功能简介】绘制三维彗星图。
【语法格式】
0 w/ V/ {: J9 A6 `9 Z1.comet3 (z)
显示向量z的三维彗星图，X与Y轴范围由系统自动确定。
( i4 _/ _. W" M; _3 k* \1 \2.comet3(x,y,z)
: L7 s- ^4 M) y) f/ [4 O显示由x、y、z确定的三维曲线的彗星图。
4 e4 w, K" h* i3.comet3(x,y,z,p)
& k  U3 u+ e( j6 l$ Y指定彗星体的长度为p*length(y)。
6 H/ y4 w; k7 O8 E+ j0 C& S  b【实例7.38】绘制一个三维彗星图。
# O9 a2 o6 i* ^, a# X>> t=-10*pi:pi/250:10*pi;
>> comet3((cos(2*t).^2).*sin(t),(sin(2*t).^2).*cos(t),t);                %绘制彗星图
执行结果如图7-39所示。
4 o4 F. p, g2 z& W$ P
. ?7 k% p8 f, c" D# Y8 H图7-39  三维彗星图
3 X% k: @' `9 z7 I. Z/ ~' l【实例分析】彗星图的绘制是一个动态的过程，用户可按实例中的代码自行运行。
) W. V9 D# h: n( ?7.3.13  cylinder——生成圆柱图形
【功能简介】绘制圆柱图形。
1.[X,Y,Z]=cylinder(r,n)
返回一个半径为r、高度为1的圆柱体的X、Y、Z轴坐标值，所绘制的圆柱体的圆周有指定的n个距离相同的点。用户可用surf或mesh函数和X、Y、Z参数画出圆柱。
6 Z6 P. \' Q  @7 D/ x7 u" R) P, p) O格式变体：
        [X,Y,Z]=cylinder(r)：所绘制的圆柱体的圆周有指定的20个距离相同的点。
        [X,Y,Z]= cylinder：半径采用默认值1。
2.cylinder(…)
& V( {! h2 X& D) e5 |# `) f没有输出参量时，直接画出圆柱体。
【实例7.39】绘制一个两头粗、中间细的圆柱。
>> t=0:pi/10:2*pi;
>> [X,Y,Z]=cylinder(2+cos(t));                %用cylinder函数算出X、Y、Z参数
. D1 T* A: }% X  W9 h: d# [" V5 O6 _>> surf(X,Y,Z);                                        %画出圆柱体
! U5 F$ q& e6 f3 X6 R( J>> axis square                                        %调整坐标轴
执行结果如图7-40所示。
1 h9 D: I& X4 V
图7-40  绘制圆柱体
1 }; e* I' `9 q8 o! o* N  p4 O【实例分析】也可不带返回值，直接用cylinder(2+cos(t))画出圆柱。
9 q8 N- [9 ^" i% K, L

lili456

7.3.14  surfc——绘制阴影图及等高线
5 a7 x) n+ v1 B+ A【功能简介】绘制阴影图及等高线。先画出三维曲面图，再在下方画出二维等高线。
【语法格式】
% K! }6 h$ ]: ]+ c" D& d1.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。颜色由高度决定。
3.surfc(…,C)
图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.40】绘制三维曲面及其二维等高线        。
>> [X,Y]=meshgrid(-8:.5:8);
1 ~; V1 h' Y4 R4 x' n* |>> R=sqrt(X.^2+Y.^2)+eps;
+ H+ V0 E( ?+ a; Z& `1 s* ^>> Z=sin(R)./R;
" H% I2 q4 F" x* p>> surfc (X,Y,Z);                        %绘制三维曲面及其二维等高线               
# ~4 B- N/ v6 B' |执行结果如图7-41所示。

图7-41  绘制三维曲面及其二维等高线
9 ]; n3 }+ h  m& _& U6 {, [, `【实例分析】suefc函数调用格式与surf相同，只是多画了一个二维等高线。
7.3.15  surfl——带光照模式的曲面图
【功能简介】绘制带光照模式的曲面图。
4 t7 p0 F& V; j( A7 n1 L" R【语法格式】
1.surfl(X,Y,Z)
, r$ L1 s5 S5 P& I: T生成由X、Y和Z定义的曲面图，光照方向和系数采用默认值。
格式变体：
* x7 z: H1 W( F2 R2 J1 g        surfl (Z)：X、Y的值由系统自动生成，等于Z中元素的索引。
6 X+ b2 M/ c1 \2.surfl(…,'light')
& Y+ A- p$ c% v) @( S1 F, V用MATLAB光照对象生成一个带光照的曲面。
% v# f8 @& T: I+ @3.surfc(…,s)
指定光源的方向s。s是长度为2或3的向量，即s=[azimuth,elevation]或[sx,sy,sz]。默认光源方向是从当前视角开始逆时针45度方向。
4.surfc(…,s,k)
  ?$ M& M1 q) t( c% @' ?3 a6 B指定反射系数k，k是一个常量。k为长度为4的向量[ka,kd,ks,shine]，四个元素分别表示环境光系数、漫反射系数、镜面反射系数和镜面反射亮度，默认值为[0.55,0.6,0.4,10]。
【实例7.41】对peaks函数表示的曲面用surfl函数进行描绘。
9 V/ D6 E. m( Q. o>> [x,y]=meshgrid(-3:1/8:3);
3 u6 T5 D4 J* {>> z=peaks(x,y);                %peaks函数
/ i6 I4 V  r% J2 B>> surfl(x,y,z);                        %绘制带光照的曲面图
>> shading interp;
>> colormap gray
执行结果如图7-42所示。

; q2 w4 _0 z8 E; P/ u图7-42  绘制带光照的曲面图
【实例分析】命令colormap gray将色图指定为灰度色图。
$ P1 d! a, \' [: ^