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

lili456

7.3  三维图形
& D2 @* Z  J, j/ X) z2 l三维图形表现的是三个变量之间的依赖关系，使用户能够直观看出二元函数数值的变化趋势。下面的小节将介绍三维图形的绘制函数。
* w2 s# b2 {2 D, y0 V) ?7.3.1  plot3——绘制三维曲线
【功能简介】绘制三维曲线。
【语法格式】
" G& |0 [; K% E& H6 W1.plot3(X1,Y1,Z1,…)
- A7 D# X* ^: {! N, fX1、Y1、Z1是同型的向量或矩阵，函数在三维空间中绘制出一条或多条曲线，绘制多条曲线时，曲线的条数等于矩阵的列数。
0 q/ W9 B0 }1 h1 z( i格式变体：
        plot3(X1,Y1,Z1,LineSpec)：用参数LineSpec指定画线的线型、标记符号和颜色。
8 c: [. E3 Y$ h( g" k2.plot3(…,'PropertyName','PropertyValue')
设置所绘图形指定属性的属性值。
【实例7.27】绘制三维螺旋图。
2 I& C& H) t5 D* t/ r* y9 E& P0 @>> t=0:pi/50:10*pi;
>> plot3(sin(t),cos(t),t);        %三维螺旋图
>> grid on;
2 _5 n" ?/ J3 L4 u5 H, s>> axis square;
执行结果如图7-28所示。
: ]4 T! X) J& w9 m! H
+ [1 H, a( r9 z 图7-28  三维螺旋图
+ d, }3 g1 J0 F' V; C【实例分析】三维螺旋图的参数方程形式为x=sin(t),y=cos(t),z=t。

lili456

7.3.2  mesh——绘制三维网格图
【功能简介】绘制三维网格图。
, B6 J, `4 Z. f  U# g1 M【语法格式】
1.mesh(X,Y,Z)
$ m: S& P* s  h2 p) i9 m, 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)]。网格线的颜色由Z定义， 颜色与高度成比例。
2.mesh(Z)
, P$ C# ?! G( ]8 ?X与Y自动生成。[m,n]=size(Z)，则X=1:n，Y=1:m。颜色由高度决定。
3.mesh(…,C)
. v5 h- Q1 a% |' h" O' H& K# ]图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
+ @0 q$ B1 J& Q# Q【实例7.28】绘制三维网格图。
>> [X,Y]=meshgrid(-8:.5:8);        %构造X、Y矩阵
>> R=sqrt(X.^2+Y.^2)+eps;
/ c+ C# ~5 Q$ N4 U2 J>> Z=sin(R)./R;
>> mesh(X,Y,Z);                        %绘制三维网格图
) R) H! [% J- A5 T1 G0 K: m执行结果如图7-29所示。


2 J0 L3 \- u# A5 O. G0 Q图7-29  三维网格图
【实例分析】不指定颜色时，网格颜色由高度决定。
) p1 x' k( N( a2 d% J) `5 A7 D6 T8 ?

lili456

7.3.3  surf——三维曲面图
) b* s" `) Q; ]7 ~& Y【功能简介】绘制三维曲面图。
【语法格式】
" h2 r2 _7 G' z" ~1.surf(X,Y,Z)
: |+ Y) ~" n2 N- E3 F5 |  z+ h) c生成由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)
% g+ r( ~$ c8 Z图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
& f. G. B! S# v& N7 u6 R; v【实例7.29】绘制三维曲面图。
' O) |! O. Q, d# M" i7 O+ k/ t>> [X,Y]=meshgrid(-8:.5:8);        %生成X、Y矩阵
>> R=sqrt(X.^2+Y.^2)+eps;
  _4 ~$ V% M. r1 l>> Z=sin(R)./R;
>> surf(X,Y,Z);                                %绘制三维曲面图
执行结果如图7-30所示。
0 {4 F$ F. R  L% J$ q
图7-30  三维曲面图
- s5 S5 R, r' a, K1 Z【实例分析】surf函数与mesh函数调用格式是相同的。

lili456

7.3.4  contour3——三维等高线绘制
( ]" D1 O% O# Z: u# b1 w8 v【功能简介】绘制三维等高线。
【语法格式】
1.contour3(Z)
& l5 a3 L, f- V画出矩阵Z的等高线图，矩阵Z至少为2×2大小，系统按矩阵Z中元素距离XY平面的高度，自动选择等高线的条数进行绘制。X和Y的值是自动确定的，[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
' i  G6 j. @* e格式变体：
. @  q2 p5 Z% B; t2 x; T- Q7 N$ W        contour3(Z,n)：画出矩阵Z的n条等高线。
        contour3(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
# N& j. q6 @4 m: O: D8 n2.contour3(X,Y,Z)
画出矩阵Z的等高线图，X轴和Y轴的范围由参数X、Y指定。如果X为矩阵，则使用X(1,来定义X轴的范围，如果Y为矩阵，则使用Y(:,1)来定义Y轴的范围。如果X与Y均为矩阵，则两者必须同型。
格式变体：
! ?8 T, Z' v1 |, x- l1 x4 h* [        contour3(X,Y,Z,n)：画出矩阵Z的n条等高线。
+ S. O6 \! q7 ^9 e9 f  I        contour3(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
3.contour3(…,LineSpec)
# p  v8 A. `/ U( ~" A9 ~参数LineSpec指定线型、标记符号和颜色。
【实例7.30】绘制等高线。
) @- ^5 Q0 n; J: L# z3 s2 Z>> [X,Y]=meshgrid([-2:.25:2]);
>> Z=X.*exp(-X.^2-Y.^2);
# D* W( a, y( Y. N>> contour3(X,Y,Z,30);        %绘制等高线
>> colormap cool;
/ ]1 j/ I  M, Q执行结果如图7-31所示。
; ]6 l( v' N! S( t4 ]
图7-31  contour3函数绘制等高线
, a  K) z* Z7 O' O【实例分析】等高线绘制有contour3和contour两个函数均可完成，contour函数绘制的是二维的等高线。
1 s( C* R4 o$ U) Z* F* ^6 y" r: Z" t

lili456

7.3.5  contour——曲面的等高线
【功能简介】绘制曲面的等高线。
7 U, z  v$ z& K7 V6 i) U9 e【语法格式】
1.contour(Z)
+ ?& f- l  S6 J: u2 @! k6 |画出矩阵Z的二维等高线图，矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
0 {5 l* T3 o  R% s格式变体：
        contour(Z,n)：画出矩阵Z的n条等高线。
& Y8 ], R; C& {5 y+ [        contour(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
! P$ _" J# [% t1 i2.contour(X,Y,Z)
画出矩阵Z的二维等高线图，X轴和Y轴的范围由参数X、Y指定。如果X与Y均为矩阵，则两者必须同型且单调递增。
格式变体：
        contour(X,Y,Z,n)：画出矩阵Z的n条等高线。
        contour(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
/ t5 t2 T. t9 W5 b0 v7 V3 V3.contour(…,LineSpec)
参数LineSpec指定线型、标记符号和颜色。
0 D  k' r; q+ q7 ?+ V. \【实例7.31】画出曲面的等高线图。
0 y3 v; |0 c6 f* |3 ]! L- h# w>> [X,Y]=meshgrid([-2:.25:2]);
9 c/ p1 V, l! E* h# @' x- k& l8 [>> Z=X.*exp(-X.^2-Y.^2);
' G* y7 q; t6 h) _' |>> contour(X,Y,Z,30);        %画出曲面的二维等高线
执行结果如图7-32所示。

# j5 g2 x9 U+ s. N# |  ~+ W% K: A& e 图7-32  曲面的等高线
【实例分析】contour函数画出的是二维图形，contour3函数画出的是三维图形。

lili456

7.3.6  clabel——等高线高度标签
; p) w" `% q# L7 k) W【功能简介】在二维等高线图中绘制等高线高度标签。
( J4 t2 q# [$ g1 F* B7 m6 P4 {% \【语法格式】
' X) \1 X* R9 |  _; ~& X6 L1.clabel(C,h)
7 J2 {! [  D5 O7 g3 a8 BC为标签矩阵，h为等高线句柄。函数旋转标签到适当角度，在空间允许的情况下插入标签。
6 C* A4 f( X9 c- L1 H( I( B+ i( Q格式变体：
, L4 ~: }. k. M4 H; o* x        clabel(C,h,v)：在高度v处插入标签C。
, H: x: o6 O$ v7 v; i        clabel(C,h,'manual')：手动添加标签到句柄h指定的等高线图中，用鼠标左键或空格键在最接近位置上放置标签，用回车键结束操作。
0 N! |- x' ]7 Y2.clabel(C)
* `* P. V1 v( z* P, M- L  x' D2 T+ p添加标签到当前等高线图中，随机选择标签位置。
格式变体：
        clabel(C,v)：在v指定的高度处添加标签。
        clabel(C,'manual')：手动添加标签到当前等高线图中。
【实例7.32】给等高线做标注。
" t# m# P2 u! @) b7 W>> [x,y,z]=peaks;
2 L! k3 h7 o5 i; z4 q6 ?% A" ?8 P>> C=contour(x,y,z);        %绘制等高线
>> clabel(C);                %给等高线添加标注
+ X2 S  M% a. f6 f) g: @8 u执行结果如图7-33所示。
( F( y7 {' w' O0 K
图7-33  给等高线添加标注
【实例分析】peaks是MATLAB自带函数。
- Z3 P8 H9 W6 c7 T7.3.7  contourc——等高线图形计算
7 c# a* T1 @3 ]: I2 Y【功能简介】计算等高线矩阵C，用函数contour、contour3和contourf来显示。
【语法格式】
1.C=contourc(Z)
从矩阵Z中计算等高矩阵C，Z可视为XY平面的高度矩阵。等高线的数量和高度值是系统自动确定的。Z至少为2×2矩阵，至少包含两个不同的值，X和Y的范围是1:n与1:m，其中[m,n]=size(Z)。
, ^1 G* x+ e7 W; b* _  ?0 W格式变体：
# V" n2 o3 ~5 n: N/ W        C=contour(Z,n)：确定等高线条数为n，返回等高矩阵。
+ D& U! F; I; f1 Z        C=contour(Z,v)：在高度v处计算等高线。
2.C=contourc(X,Y,Z)
在X和Y上计算Z的等高矩阵C，Z可视为XY平面的高度矩阵。等高线的数量和高度值是系统自动确定的。Z至少为2×2矩阵，至少包含两个不同的值。
* {7 K4 S* O, z* G9 l格式变体：
2 Q& Y6 a- W( H& ^6 Z; Y* F        C=contour(X,Y,Z,n)：确定等高线条数为n，返回等高矩阵。
        C=contour(X,Y,Z,v)：在高度v处计算等高线。
6 ^+ {5 Y* Y8 [2 Q9 x, U# K【实例7.33】计算peaks函数的等高矩阵。
>> a=peaks;
>> c=contourc(a,10);                %计算10条等高线的等高矩阵
>> s=size(c)                        %矩阵维度
4 \$ |% P1 ^7 C: p1 e: |( Ts =
     2   800
>> contour(a,c)                        %画出等高线
执行结果如图7-34所示。

" Q3 J, G) Z/ a  O; A$ g
图7-34  计算等高矩阵再绘制等高线
【实例分析】矩阵C是一个2×m矩阵，m的值取决于数据矩阵及等高线的条数。

lili456

7.3.8  fill3——填充三维图
$ @; B# v, w6 B$ K+ g! d【功能简介】填充三维多边形。
8 S- m. P% m2 H* f$ L【语法格式】
$ Q4 q: z, g3 ]) |4 y1.fill3(X,Y,Z,C)
+ H% z7 }" ?: T5 D( e& m* O参数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)。
# j) G% x; e2 \2 r3 {; _格式变体：
& e* B) J0 t+ e2 `6 ?* C0 e        fill3(X1,Y1,Z1,C1,X2,Y2,Z2,C2,…)：绘制多个三维多边形。
5 Y! I; P6 D$ t7 n2.fill3(X,Y,Z,ColorSpec)
) y3 K/ U9 g$ t) ~用参数ColorSpec指定的颜色对三维多边形进行填充。
3 i% ~9 g' Y0 l' h& o. g0 Y【实例7.34】用渐变的颜色来对三维多边形进行着色。
>> 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];
>> C=[.5,1,1,.5;1,.5,.5,.1667;.3333,.3333,.5,.5]        %指定颜色
C =
- r- S* k5 \, k$ X& T7 X    0.5000    1.0000    1.0000    0.5000
    1.0000    0.5000    0.5000    0.1667
( z# \( ~. m8 r0 b0 E  W( U0 h2 S    0.3333    0.3333    0.5000    0.5000
>> fill3(X,Y,Z,C)                                                %填充三维多边形
执行结果如图7-35所示。
/ o$ L8 O) ~5 E" z3 e
图7-35  填充三维多边形
  b! f8 Q. v% p  _4 C9 L  J【实例分析】X、Y、Z决定三维多边形的形状，C确定颜色。
: x. v' D; ]; C2 v; C7 k( _" U7.3.9  sphere——绘制球体
6 p' b1 q4 j4 P7 l【功能简介】绘制球体。
【语法格式】
( |) \9 N5 M) W- G3 [+ G1.sphere(n)
* o# W3 R6 b) Y, a在三维直角坐标系中绘制由n×n个面组成的单位球体。
格式变体：
        sphere：默认n=20，绘制单位球体。
+ L2 ]' V; e2 w. {% H: H) C2.[X,Y,Z]=sphere(n)
: g3 ?  ]: ]+ G* [2 [) @返回三个阶数为(n+1)×(n+1)的坐标矩阵。该命令不画图，只是返回矩阵，要画图可以调用mesh(X,Y,Z)或surf(X,Y,Z)。
【实例7.35】绘制多个球体。
>> [x,y,z]=sphere;
>> surf(x,y,z);                %画出第一个球体
>> hold on;
' j6 c1 w  P: y( [$ t>> surf(x+3,y+2,z);        %画出第二个球体
( z& S9 G$ m) ~. R>> surf(x,y-1,z+2);        %画出第三个球体
>> daspect([1 1 1])
执行结果如图7-36所示。
7 Y  `; R) m" B/ j" e" v
* ~* W2 \$ Y( I% M图7-36  绘制多个球体
【实例分析】sphere绘制单位球体，半径为1。
! A. S" `# ?9 `* D( l. m4 M  O

lili456

7.3.10  contourf——填充二维等高线
【功能简介】填充二维等高线。
% A& }; `+ A+ D. J$ o【语法格式】
/ t0 Z, i* c( U; I1.contourf(Z)
画出矩阵Z的二维等高线图，再对等高线之间的区域进行填充，填充的颜色使用当前窗口的颜色映射表。矩阵Z至少为2×2大小，可视为XY平面的高度矩阵。等高线的个数和值是基于Z的最大值和最小值自动选取的。[m,n]=size(Z)，X轴的范围是1:n，Y轴的范围是1:m。
# a+ |; `" N! V: u# L% u. |1 o格式变体：
' x% }. e6 H+ i/ t9 k2 Y% P8 [        contourf(Z,n)：画出矩阵Z的n条等高线并对等高线间的区域进行填充。
        contourf(Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
2.contourf(X,Y,Z)
画出矩阵Z的等高线图并进行填充，X轴和Y轴的范围由参数X、Y指定。如果X与Y均为矩阵，则两者必须同型且单调递增。
9 n, Y3 C% _7 K/ ]$ \. C: x& W6 |. U格式变体：
4 M4 g* F) _8 p6 C8 @0 R0 P( i        contourf(X,Y,Z,n)：画出矩阵Z的n条等高线并进行填充。
( q) ~' A+ S0 ^8 h2 H        contourf(X,Y,Z,v)：参数v确定绘制等高线的高度，等高线的条数等于length(v)。
4 @  O& y; q( E【实例7.36】画出等高线并进行填充。
>> a=peaks;
>> contourf(a);                %画出用颜色填充的等高线
执行结果如图7-37所示。

9 {8 r2 y( ]! H- K! m. G* Y+ r2 H图7-37  用颜色填充的等高线图
$ z% i1 j# @7 d7 C【实例分析】用颜色对等高线间的区域进行填充，使图形更直观。
7.3.11  pie3——三维饼图
【功能简介】绘制三维饼图。
【语法格式】
1.pie3(X)
7 r+ l( |! Q3 h# k/ _  v1 ^使用X中的数据绘制一个三维饼图，X中的每一个元素都是饼图的一个部分。元素所占扇形面积的大小由X(i)/sum(X)决定，如果sum(X)<1，则只绘制不完整的三维饼形图。
2.pie3(X,explode)
指定三维饼图中的每一部分是否分离出来。explode是与X同型的数组，用非零值表示相应位置的X中的元素分离出来显示。
% L3 r) Q' [( N【实例7.37】绘制三维饼图。
9 w+ I+ w! y( d' V: P0 Q$ I8 j>> x=[1,3,.7,2.5,2];
>> explode=[0,1,0,0,1];        %设置第二个元素和最后一个元素对应的区域分离出来独立显示
>> pie3(x,explode);                %绘制三维饼图
9 L* J; _; m, u>> colormap hsv
执行结果如图7-38所示。

图7-38  三维饼图
【实例分析】pie3的调用格式与pie类似。

lili456

7.3.12  comet3——三维彗星图
3 e: r1 ], p/ H, y& `0 N- w【功能简介】绘制三维彗星图。
6 X# v7 J1 T6 q2 u/ S1 m【语法格式】
+ w& ^) t4 F1 x1.comet3 (z)
显示向量z的三维彗星图，X与Y轴范围由系统自动确定。
2.comet3(x,y,z)
8 g+ v7 g5 V& N0 A; B0 }" I: a" y显示由x、y、z确定的三维曲线的彗星图。
3.comet3(x,y,z,p)
0 ~, u! U! }+ x6 S* Y( A指定彗星体的长度为p*length(y)。
【实例7.38】绘制一个三维彗星图。
>> t=-10*pi:pi/250:10*pi;
  V/ ^1 y' `( G0 `8 k>> comet3((cos(2*t).^2).*sin(t),(sin(2*t).^2).*cos(t),t);                %绘制彗星图
7 N( `# Z# M9 K- }. h0 E7 L1 e  d执行结果如图7-39所示。
, q- U9 I3 j5 ?/ @, h
图7-39  三维彗星图
【实例分析】彗星图的绘制是一个动态的过程，用户可按实例中的代码自行运行。
7.3.13  cylinder——生成圆柱图形
5 W& r3 i( l4 S. V; R( q6 ~0 N2 {【功能简介】绘制圆柱图形。
1.[X,Y,Z]=cylinder(r,n)
返回一个半径为r、高度为1的圆柱体的X、Y、Z轴坐标值，所绘制的圆柱体的圆周有指定的n个距离相同的点。用户可用surf或mesh函数和X、Y、Z参数画出圆柱。
格式变体：
        [X,Y,Z]=cylinder(r)：所绘制的圆柱体的圆周有指定的20个距离相同的点。
! X$ ?7 }6 U" l( `2 ?/ e9 t        [X,Y,Z]= cylinder：半径采用默认值1。
! R" _* s% c, O* Y" M2.cylinder(…)
没有输出参量时，直接画出圆柱体。
【实例7.39】绘制一个两头粗、中间细的圆柱。
# Q/ I+ U! [* E& z* U>> t=0:pi/10:2*pi;
: v  m; z% e7 q4 k- b: M>> [X,Y,Z]=cylinder(2+cos(t));                %用cylinder函数算出X、Y、Z参数
>> surf(X,Y,Z);                                        %画出圆柱体
2 J$ }' j8 _- {: _>> axis square                                        %调整坐标轴
" g. f% R( t4 g, X$ Y7 Q& A执行结果如图7-40所示。

  G3 j3 ]0 R4 F3 Y/ D图7-40  绘制圆柱体
6 X* L( x( S# a0 @1 M5 M7 @% c/ b5 r【实例分析】也可不带返回值，直接用cylinder(2+cos(t))画出圆柱。

lili456

7.3.14  surfc——绘制阴影图及等高线
1 [; a: n0 {; P【功能简介】绘制阴影图及等高线。先画出三维曲面图，再在下方画出二维等高线。
【语法格式】
7 E% `) v* r) ?: B8 Z* u, m1.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)
! z" q4 K# X8 @3 vX与Y自动生成。[m,n]=size(Z)，则X=1:n，Y=1:m。颜色由高度决定。
" w" C8 d( Q+ _9 x3.surfc(…,C)
. j" A4 _9 a+ k* G图形颜色由矩阵C决定。如果X、Y和Z也是矩阵，则四个矩阵必须同型。
【实例7.40】绘制三维曲面及其二维等高线        。
, T# \3 `8 Y0 O* c>> [X,Y]=meshgrid(-8:.5:8);
>> R=sqrt(X.^2+Y.^2)+eps;
>> Z=sin(R)./R;
2 k0 D; [; f! }0 t2 R! M7 H7 y>> surfc (X,Y,Z);                        %绘制三维曲面及其二维等高线               
执行结果如图7-41所示。
7 q3 @* @6 L! i
图7-41  绘制三维曲面及其二维等高线
  P: C& ?) E  }8 u【实例分析】suefc函数调用格式与surf相同，只是多画了一个二维等高线。
3 S( [$ _% L3 q/ r7.3.15  surfl——带光照模式的曲面图
) r* y2 Z( |; }% q. B【功能简介】绘制带光照模式的曲面图。
【语法格式】
1.surfl(X,Y,Z)
$ Z$ Z& B$ k; J/ L' L" y. G) C% P( d生成由X、Y和Z定义的曲面图，光照方向和系数采用默认值。
' _4 w4 H! J7 u1 n- Y' K- L) M格式变体：
; J* A# i) m* K' u( f        surfl (Z)：X、Y的值由系统自动生成，等于Z中元素的索引。
3 {6 [' p" r* k, ?) a* h1 l, z2.surfl(…,'light')
用MATLAB光照对象生成一个带光照的曲面。
3.surfc(…,s)
0 ^1 f$ V. J; e, x指定光源的方向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]。
, P# A7 c& J. q7 P- [7 f6 W【实例7.41】对peaks函数表示的曲面用surfl函数进行描绘。
>> [x,y]=meshgrid(-3:1/8:3);
>> z=peaks(x,y);                %peaks函数
! c( Q* ?( r; H  g! u6 T>> surfl(x,y,z);                        %绘制带光照的曲面图
7 m5 F, y/ r6 N>> shading interp;
+ q& K" v: N8 X5 M6 }9 l2 q1 P! S>> colormap gray
# r8 m, `" w8 R0 |执行结果如图7-42所示。

图7-42  绘制带光照的曲面图
【实例分析】命令colormap gray将色图指定为灰度色图。