MATLAB数据和函数的可视化.ppt

《MATLAB数据和函数的可视化.ppt》由会员分享，可在线阅读，更多相关《MATLAB数据和函数的可视化.ppt（106页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第6章 数据和函数的可视化matlab语言丰富的图形表现方法，使得数学计算结果可以方便地、多样性地实现了可视化，这是其它语言所不能比拟的。matlab语言的绘图功能 不仅能绘制几乎所有的标准图形，而且不仅能绘制几乎所有的标准图形，而且其表现形式也是丰富多样的。其表现形式也是丰富多样的。matlabmatlab语言不仅具有高层绘图能力，而语言不仅具有高层绘图能力，而且还具有底层绘图能力且还具有底层绘图能力句柄绘图方句柄绘图方法。法。在面向对象的图形设计基础上，使得用在面向对象的图形设计基础上，使得用户可以用来开发各专业的专用图形。户可以用来开发各专业的专用图形。6.1 引导6.1.1 离散数据和

2、离散函数的可视化6.1.2 连续函数的可视化6.1.3 可视化的一般步骤6.2 二维曲线绘图的基本操作6.2.1 plot命令基本调用格式6.2.2 曲线的色彩、线型和数据点型6.2.3 坐标、刻度和分格线6.2.4 图形标志6.2.5 多次叠绘、双纵坐标和多子图6.2.6 交互图形指令6.2.1 plot命令基本调用格式（1）plot 最基本的二维图形指令plot的功能：plot命令自动打开一个图形窗口Figure 用直线连接相邻两数据点来绘制图形根据图形坐标大小自动缩扩坐标轴，将数据标尺及单位标注自动加到两个坐标轴上，可自定坐标轴，可把x,y 轴用对数坐标表示plot的调用格式 plot(

3、x)缺省自变量绘图格式，x为向量,以x元素值为纵坐标，以相应元素下标为横坐标绘图 plot(x,y)基本格式，以y(x)的函数关系作出直角坐标图，如果y为nm的矩阵，则以x 为自变量，作出m条曲线plot(x1,y1,x2,y2)多条曲线绘图格式 plot(x,y,s)开关格式，开关量字符串s设定曲线颜色和绘图方式，使用颜色字符串的前13个字母，如 yellowyel表示等。或plot(x1,y1,s1,x2,y2,s2，)6.2.2 曲线的色彩、线型和数据点型1.色彩和线型色彩：符号：b,g,r,c,m,y,k,w 含义：蓝 绿 红 青 品 黄 黑 白线型：符号：-,:,-.,-含义：实线

4、虚线 点划 双划 6.2.2 曲线的色彩、线型和数据点型2.数据点型符号含义符号含义.实心黑点d菱形+十字符h六角心符*八线符o空心圆圈朝上三角符p五角星符朝右三角符x叉字符v朝下三角符【例6.2.2.2-1】用图形演示平面上一个方块四个顶点在仿射投影（Affine Projection）下的位置、形状变化。分别实现矩阵平移、旋转和刻度变换%平面上四个点和它们构成的方块平面上四个点和它们构成的方块p1=-0.5,0,1;p2=-0.5,1,1;p3=0.5,1,1;p4=0.5,0,1;Sq=p1,p2,p3,p4,p1;%平移投影平移投影：dx=0.5;dy=1;T=1,0,dx;0,1,d

5、y;0,0,1;%旋转投影：逆时针旋转旋转投影：逆时针旋转30th=pi/6;R=cos(th),-sin(th),0;sin(th),cos(th),0;0,0,1;%刻度投影：刻度投影：alpha=2;beta=3;S=alpha,0,0;0,beta,0;0,0,1;E=eye(3,3);%TRS=E,T,R,S;ss=r,rd,rp,rh;tt=Original Square,Translation,Rotation,Scaling;for i=1:4W=TRSi*Sq;%进行仿射变换进行仿射变换subplot(2,2,i)for k=1:4plot(W(1,k),W(2,k),ssk

6、);%用不同点形标志四个顶点用不同点形标志四个顶点axis(-3,3,-1,5),axis equalhold onendplot(W(1,:),W(2,:)%连接四个顶点使得图形封闭连接四个顶点使得图形封闭grid on%title(tti)%hold off%end6.2.3 坐标、刻度和分格线控制1、坐标控制2、刻度、分格线和坐标框坐标控制【例6.2.3.1-1】观察各种轴控制指令的影响。演示采用长轴为3.25，短轴为1.15 的椭圆。注意：采用多子图表现时，图形形状不仅受“控制指令”影响，而且受整个图面“宽高比”及“子图数目”的影响。本书这样处理，是出于篇幅考虑。读者欲想准确体会控制指

7、令的影响，请在全图状态下进行观察。（图6.2-6）t=0:2*pi/99:2*pi;x=1.15*cos(t);y=3.25*sin(t);%subplot(2,3,1),plot(x,y),axis normal,grid on,title(Normal and Grid on)subplot(2,3,2),plot(x,y),axis equal,grid on,title(Equal)subplot(2,3,3),plot(x,y),axis square,grid on,title(Square)subplot(2,3,4),plot(x,y),axis image,box off,t

8、itle(Image and Box off)subplot(2,3,5),plot(x,y),axis image fill,box offtitle(Image and Fill)subplot(2,3,6),plot(x,y),axis tight,box off,title(Tight)2、刻度、分格线和坐标框（1）刻度设置set(gca,xtick,xs,ytick,ys,ztick,zs);set(gca,xticklabel,sx,);(2)分格线gridgrid on(off)(3)坐标框boxbox on(off)【例例6.2.3.2-1】通过绘制二阶系统阶跃响应，综合演通过

9、绘制二阶系统阶跃响应，综合演示图形标识。本例比较综合，涉及的指令较广。请示图形标识。本例比较综合，涉及的指令较广。请读者耐心读、实际做、再看例后说明，定会有匪浅读者耐心读、实际做、再看例后说明，定会有匪浅收益。（图收益。（图6.2-7）clf;t=6*pi*(0:100)/100;y=1-exp(-0.3*t).*cos(0.7*t);tt=t(find(abs(y-1)0.05);ts=max(tt);plot(t,y,r-,LineWidth,3)axis(-inf,6*pi,0.6,inf)set(gca,Xtick,2*pi,4*pi,6*pi,Ytick,0.95,1,1.05,ma

10、x(y)grid on%title(it y=1-e-alphatcosomegat)%text(13.5,1.2,fontsize12alpha=0.3)%text(13.5,1.1,fontsize12omega=0.7)%hold on;plot(ts,0.95,bo,MarkerSize,10);hold offcell_string1=fontsize12uparrow;cell_string2=fontsize16 fontname隶书隶书镇定时间镇定时间;cell_string3=fontsize6 ;cell_string4=fontsize14rmt_s=num2str(ts

11、);text(ts,0.85,cell_string)xlabel(fontsize14 bft rightarrow)ylabel(fontsize14 bfy rightarrow)latextext=(s,p)(text(Interpreter,latex,String,s,Position,p,fontsize,14);latextext($-pi$,-pi-.25,-1.1);latextext($-fracpi2$,-pi/2-.25,-1.1);latextext($fracpi2$,pi/2-.1,-1.1);latextext($pi$,pi-.1,-1.1);set(gca

12、,xtick,-pi -pi/2 0 pi/2 pi);set(gca,xticklabel,0 );axis(-pi pi,-1 1);6.2.4 图形标识图形标识1、简捷指令形式titlexlabel,ylabel,zlabellegend(s1,s2,)text(xt,yt,s)2.精细指令形式精细指令形式(1)允许标识多行字符允许标识多行字符指令指令arg取值取值 举举 例例示例指令示例指令效果效果单单行行arg任何任何合法合法字符字符Single lineSingle line多多行行arg1,arg2元胞数组元胞数组,标志标志,Multiline元胞数组元胞数组标志标志Multil

13、inearg1;arg2元元 胞胞 数数 组组;标标 志志;Multiline元元 胞胞 数数 组组标标 志志Multiline(2)允许标识字体允许标识字体,风格及大小进行设风格及大小进行设置置字字体体指令指令arg取值取值 举举 例例示例指令示例指令效果效果名名称称fontnameargarial;courier;roman;宋体宋体;隶书隶书;黑体黑体fontnamecourierExample 1风风格格argbf-黑体黑体lt-斜体斜体sl斜体斜体rm正体正体bfExample 3大大小小fontsizearg正整数正整数缺省值为缺省值为0fontsize14 Example 5(3

14、)允许使用上下标允许使用上下标指令指令arg取值取值 举举 例例示例指令示例指令效效果果上标上标arg任何合法字符任何合法字符ite-tsint下标下标_arg任何合法字符任何合法字符xchi_alpha2(3)(4)允许标识希腊字符和其他特殊字符允许标识希腊字符和其他特殊字符6.2.5 多次叠绘、双纵坐标和多子图1、多次叠绘hold hold on(off)2.多窗口绘图 figure(n)创建窗口函数，n为窗 口顺序号。t=0:pi/100:2*pi;y=sin(t);y1=sin(t+0.25);y2=sin(t+0.5);plot(t,y)自动出现第一个窗口figure(2)plot(

15、t,y1)在第二窗口绘图figure(3)plot(t,y2)在第三窗口绘图例：t=0:0.1:10y1=sin(t);y2=cos(t);plot(t,y1,r,t,y2,b-);x=1.7*pi;1.6*pi;y=-0.3;0.8;s=sin(t);cos(t);text(x,y,s);title(正弦和余弦曲线);legend(正弦,余弦)xlabel(时间t),ylabel(正弦、余弦)gridaxis square3.双纵坐标图plotyy(X1,Y1,X2,Y2)plotyy(X1,Y1,X2,Y2,FUN)plotyy(X1,Y1,X2,Y2,FUN1,FUN2)4.多子图sub

16、plot 子图分割命令调用格式：subplot(m,n,p)按从左至右，从上至下排列subplot(position,left bottom width height)6.2.6 交互式图形指令ginputx,y=ginput(n)用鼠标从二维图形上获取n个点的数据坐标(x,y).gtextgtext(arg)用鼠标把字符串元胞数组arg放置到图形上.legendlegend(arg,pos);legend off;zoom6.3三维绘图的基本操作6.3.1 三维线图指令plot36.3.2三维网线图和曲面图6.3.3透视、镂空和裁切6.3.1 三维线图指令plot3t=(0:0.02:2)*

17、pi;x=sin(t);y=cos(t);z=cos(2*t);plot3(x,y,z,b-,x,y,z,bd)view(-82,58)box onlegend(链,宝石)6.3.2三维网线图和曲面图三维网线图和曲面图用曲面图表现函数。clf,x=-4:4;y=x;X,Y=meshgrid(x,y);Z=X.2+Y.2;surf(X,Y,Z);hold on,colormap(hot)stem3(X,Y,Z,bo)（1）三维网线图mesh 三维网线绘图函数调用格式：mesh(z)z为nm的矩阵，x与y 坐标为元素的下标 mesh(x,y,z)x,y,z分别为三维空 间的坐标位置三维网线图作图要

18、领生成坐标 X,Y=meshgrid(x,y)表达式点运算 Z=X.2+Y.2X，Y是nm的矩阵，维数可任定X nY mZ默认方位角：37。5，俯角30。nm维 meshgrid 网线坐标值计算函数 z=f(x,y)根据x,y坐标找出z的高度例：绘制z=x2+y2的三维网线图形 x=-5:5;y=x;X,Y=meshgrid(x,y)Z=X.2+Y.2 mesh(X,Y,Z)坐标矩阵坐标向量纵坐标矩阵绘图函数（2）三维曲面图surf 三维曲面绘图函数，与网格图 看起来一样与三维网线图的区别：网线图：线条有颜色，空挡是黑色的（无颜色）曲面图：线条是黑色的，空挡有颜色（把线条之间的空挡填充颜色，沿

19、z轴按每一网格变化调用格式：surf(x,y,z)绘制三维曲面图，x,y,z 为图形坐标向量例：X,Y,Z=peaks(30)peaks为matlab自动生成的三维测试图形surf(X,Y,Z)surfc(X,Y,Z)带等高线的曲面图 X,Y,Z=peaks(30);surfc(X,Y,Z)surfl(X,Y,Z)被光照射带阴影 的曲面图X,Y,Z=peaks(30);surfl(X,Y,Z)6.3.3透视、镂空和裁切透视、镂空和裁切透视演示透视演示X0,Y0,Z0=sphere(30);X=2*X0;Y=2*Y0;Z=2*Z0;surf(X0,Y0,Z0);shading interphol

20、d on,mesh(X,Y,Z),colormap(hot),hold offhidden offaxis equal,axis off 演示：如何利用演示：如何利用“非数非数”NaN，对图形进，对图形进行行镂空镂空处理。处理。P=peaks(30);P(18:20,9:15)=NaN;surfc(P);colormap(summer)；light(position,50,-10,5)；lighting flat；material(0.9,0.9,0.6,15,0.4)；表现切面（将切面部分强制为零）clf,x=-8:0.05:8;y=x;X,Y=meshgrid(x,y);ZZ=X.2-Y.

21、2;ii=find(abs(X)6|abs(Y)6);ZZ(ii)=zeros(size(ii);surf(X,Y,ZZ);shading interp;colormap(copper);light(position,0,-15,1);lighting phong;material(0.8,0.8,0.5,10,0.5)6.4特殊图形和高维可视化特殊图形和高维可视化6.4.1特殊图形指令例示特殊图形指令例示6.4.2四维表现四维表现6.4.3动态图形动态图形6.4.16.4.1特殊图形指令例示特殊图形指令例示1.1.面域图面域图areaarea2.2.各种直方图各种直方图bar,barh,ba

22、r3,bar3hbar,barh,bar3,bar3h3.3.饼图饼图pie,pie3pie,pie34.4.填色图填色图fill,fill3fill,fill35.5.射线图射线图compasscompass和羽毛图和羽毛图featherfeather6.Vornonoi6.Vornonoi图和三角剖分图和三角剖分7.7.彩带图彩带图ribbonribbon8.8.离散杆图离散杆图stem,stem3stem,stem39.9.二维半图指令二维半图指令pcolorpcolor，contourcontour，contourfcontourf10.10.散点图散点图scatterscatter，

23、scatter3scatter3，plotmatrixplotmatrix11.11.不规则数据的网线图和曲面图不规则数据的网线图和曲面图trimeshtrimesh，trisurftrisurf1.面域图areaclf;x=-2:2Y=3,5,2,4,1;3,4,5,2,1;5,4,3,2,5Cum_Sum=cumsum(Y)area(x,Y,0)legend(因素因素A,因素因素B,因素因素C),grid on,colormap(spring)2.各种直方图bar,barh,bar3,bar3h二维直方图有两种图型：二维直方图有两种图型：垂直直方图垂直直方图和和水平直水平直方图方图。而每种

24、图型又有两种表现模式：累计式：。而每种图型又有两种表现模式：累计式：分组式。本例选其两种加以表现。分组式。本例选其两种加以表现。x=-2:2;Y=3,5,2,4,1;3,4,5,2,1;5,4,3,2,5;subplot(1,2,1),bar(x,Y,stacked)xlabel(x),ylabel(Sigma y),colormap(cool)legend(因素因素A,因素因素B,因素因素C)subplot(1,2,2),barh(x,Y,grouped)xlabel(y),ylabel(x)3.饼图pie,pie3饼图指令pie,pie3pie,pie3 用来表示各元素占总和的百分数。该指

25、令第二输入宗量为与第一宗量同长的0-1向量，1使对应扇块突出。a=1,1.6,1.2,0.8,2.1;a=1,1.6,1.2,0.8,2.1;subplot(1,2,1),subplot(1,2,1),pie(a,1 0 1 0 0);pie(a,1 0 1 0 0);legend(1,2,3,4,5)legend(1,2,3,4,5)subplot(1,2,2),subplot(1,2,2),pie3(a,a=pie3(a,a=min(amin(a),),colormap(coolcolormap(cool)4.填色图填色图fill,fill3MATLAB画任意多边形的一种方法；保证绘图画任

26、意多边形的一种方法；保证绘图数据首尾重合，使勾画多边形封闭；使用图柄数据首尾重合，使勾画多边形封闭；使用图柄对图形的属性进行精细设置。对图形的属性进行精细设置。clf;nclf;n=10;=10;dt dt=2*pi/n;=2*pi/n;t=0:dt:2*pi;t=0:dt:2*pi;t=t,t(1);t=t,t(1);x=x=sin(t);ysin(t);y=cos(tcos(t););fill(x,y,c);axisfill(x,y,c);axis off offht=text(0,0,fontnameht=text(0,0,fontname隶书隶书fontsize32fontsize32

27、十边形十边形););set(ht,Color,k,HorizontalAlignment,Centerset(ht,Color,k,HorizontalAlignment,Center)射线图compass和羽毛图feather compass和feather指令的区别。t=-pi/2:pi/12:pi/2;r=ones(size(t);x,y=pol2cart(t,r);subplot(1,2,1),compass(x,y),title(Compass)subplot(1,2,2),feather(x,y),title(Feather)6.Voronoi图和三角剖分图和三角剖分clf;ran

28、d(state,111)n=30;A=rand(n,1)-0.5;B=rand(n,1)-0.5;T=delaunay(A,B);T=T T(:,1);voronoi(A,B)hold on;axis squarefill(A(T(10,:),B(T(10,:),y);voronoi(A,B)7.彩带图ribbon x,y=meshgrid(-3:.5:3,-3:.1:3);z=peaks(x,y);ribbon(y,z)xlabel(X)ylabel(Y)zlabel(Z)colormap hsv 离散杆图stem,stem3 8.二维半图指令二维半图指令pcolor,contour,con

29、tourf 伪彩图伪彩图pcolor；等位线指令；等位线指令contour、contourf；等位；等位线标高指令线标高指令 clabel 的配合使用和区别。的配合使用和区别。（1）本例等位线指令中的第）本例等位线指令中的第4输入宗量输入宗量n设定高度的等级设定高度的等级数，第数，第5输入宗量设定等位线的线型、色彩。输入宗量设定等位线的线型、色彩。（2）左右两图的标高方法不同。左图的标识以）左右两图的标高方法不同。左图的标识以“+”引引导，水平放置。右图沿线布置。这是由导，水平放置。右图沿线布置。这是由clabel的调用的调用格式不同产生的。格式不同产生的。（3）左右两图色彩的形成方法不同，色

30、彩效果也不同。）左右两图色彩的形成方法不同，色彩效果也不同。（4）在左图中，）在左图中，colorbar画出一根垂直色标尺，而画出一根垂直色标尺，而caxis决定该色标尺的刻度决定该色标尺的刻度 等高线修饰二维contour(Z,n)-绘制n条等高线C=contourc(Z,n)-计算n条等高线的坐标Clable(c)-给等高线加标注例、在二维平面上绘制peaks函数的10条等 高线contour(peaks,10);C=contourc(peaks,10);clabel(C)clf;clear;X,Y,Z=peaks(40);n=4;subplot(1,2,1),pcolor(X,Y,Z)c

31、olormap jet,shading interphold on,C=contour(X,Y,Z,n,k:);clabel(C)zmax=max(max(Z);zmin=min(min(Z);caxis(zmin,zmax)colorbarhold off,subplot(1,2,2)C,h,CF=contourf(X,Y,Z,n,k:);clabel(C,h)10.散点图散点图scatter,scatter3,plotmatrixx=3*pi*(-1:0.2:1);y=x;X,Y=meshgrid(x,y);R=sqrt(X.2+Y.2)+eps;Z=sin(R)./R;C=abs(del

32、2(Z);meshz(X,Y,Z,C)；hold on,scatter3(X(:),Y(:),Z(:),filled)；hold off,colormap(hot)；指令指令plotmatrixplotmatrix有两种基本调用方式：有两种基本调用方式：（1 1）对于数据矩阵维的）对于数据矩阵维的X X（p*np*n）和维的）和维的Y Y（p*mp*m），调用格式），调用格式plotmatrix(X,Yplotmatrix(X,Y)将画出一将画出一个分割成（个分割成（m*nm*n）个子散点图。其中第（）个子散点图。其中第（i i，j j）个子散点图是根据个子散点图是根据Y Y第第i i列和列和

33、X X第第j j列数据画出的。列数据画出的。（2 2）对于数据矩阵维的）对于数据矩阵维的X X，调用格式，调用格式plotmatrix(Xplotmatrix(X)将画出分割成（将画出分割成（n*nn*n）个子块的图。）个子块的图。该图的对角块，画出的是该图的对角块，画出的是X X每列的数据的频数每列的数据的频数直方图；而其他子块是相应列构成的散点图。直方图；而其他子块是相应列构成的散点图。randn(seed,1111)X=randn(100,2);Y=randn(100,2);subplot(1,3,1),plotmatrix(X)subplot(1,3,2),plotmatrix(X,X

34、)subplot(1,3,3),plotmatrix(X,Y)用三角网线、曲面图表现函数。rand(seed,22),rand(seed,22),X=6*pi*(rand(20,10)-0.5);X=6*pi*(rand(20,10)-0.5);Y=6*pi*(rand(20,10)-0.5);Y=6*pi*(rand(20,10)-0.5);R=sqrt(X.2+Y.2)+eps;Z=R=sqrt(X.2+Y.2)+eps;Z=sin(Rsin(R)./R;)./R;tri=tri=delaunay(X,Ydelaunay(X,Y););subplot(1,2,1),trimesh(tri,

35、X,Y,Z)subplot(1,2,1),trimesh(tri,X,Y,Z)subplot(1,2,2),trisurf(tri,X,Y,Z)subplot(1,2,2),trisurf(tri,X,Y,Z)colormap(jet);brighten(0.5)colormap(jet);brighten(0.5)6.4.2四维表现1.用色彩表现函数特征2.切片图和切片等位线图3.1.用色彩表现函数特征surf（x,y,z）,surf(x,y,z,s)用色图阵表现函数的不同特征。演示：当三维网线图、曲面图的第四个输入宗量取一些特殊矩阵时，色彩就能表现或加强函数的某特征，如梯度、曲率、方向导数

36、等。2.切片图和切片等位线图切片图和切片等位线图利用slice和contourslice 表现MATLAB提供的无限大水体中水下射流速度数据flow。flow是一组定义在三维空间上的函数数据。本例将动用切片、视角、色彩和透明度等综合手段观察定义在三维空间上的函数。%clf;X,Y,Z,V=flow;x1=min(min(min(X);x2=max(max(max(X);y1=min(min(min(Y);y2=max(max(max(Y);z1=min(min(min(Z);z2=max(max(max(Z);sx=linspace(x1+1.2,x2,5);sy=0;sz=0;slice(X

37、,Y,Z,V,sx,sy,sz);view(-33,36);shading interp;colormap jet;%以下为透明度设置以下为透明度设置alpha(color)alphamap(rampdown)alphamap(increase,.1)colorbaraxis off clf;v1=min(min(min(V);v2=max(max(max(V);cv=linspace(v1,v2,15);contourslice(X,Y,Z,V,sx,sy,sz,cv);view(-12,30)colormap jet;colorbar;box on 6.4.3 动态图形1.彗星状轨迹图2.

38、色图的变化3.影片动画彗星状轨迹图 shg;n=10;t=n*pi*(0:0.0005:1);x=sin(t);y=cos(t);plot(x,y,g);axis square;hold oncomet(x,y,0.01);hold off卫星返回地球的运动轨线示意。shg;R0=1;a=12*R0;b=9*R0;T0=2*pi;T=5*T0;dt=pi/100;t=0:dt:T;f=sqrt(a2-b2);th=12.5*pi/180;E=exp(-t/20);x=E.*(a*cos(t)-f);y=E.*(b*cos(th)*sin(t);z=E.*(b*sin(th)*sin(t);pl

39、ot3(x,y,z,g)X,Y,Z=sphere(30);X=R0*X;Y=R0*Y;Z=R0*Z;grid on,hold on,surf(X,Y,Z),shading interpx1=-18*R0;x2=6*R0;y1=-12*R0;y2=12*R0;z1=-6*R0;z2=6*R0;axis(x1 x2 y1 y2 z1 z2)view(117 37),comet3(x,y,z,0.02),hold off 2.色图的变化peaksspinmap3.动画效果动画生成的步骤1.创建帧矩阵 moviein2.对动画中的每一帧生成图形，并把它们放到帧矩阵中 getframe3.从帧矩阵中回放

40、动画movieclfshg,x=3*pi*(-1:0.05:1);y=x;X,Y=meshgrid(x,y);R=sqrt(X.2+Y.2)+eps;Z=sin(R)./R;h=surf(X,Y,Z);colormap(jet);axis off;n=12;mmm=moviein(n);for i=1:n rotate(h,0 0 1,25);mmm(:,i)=getframe;endmovie(mmm,5,10)6.5三维图形的精细控制6.5.1、视点控制和图形的旋动视点控制和图形的旋动6.5.2、色彩控制色彩控制 6.5.3、照明和材质处理照明和材质处理 6.5.4、透明处理透明处理6.5

41、.1、视点控制和图形的旋动视点控制和图形的旋动1.视点控制viewView(az,el)通过方位角，俯视角设置azimuth，方位；elevation，俯视角View(vx,vy,vz)通过直角坐标设置6.5.1、视点控制和图形的旋动视点控制和图形的旋动2.图形旋动rotateRotate(h,direction,alpha,origin)Direction设置theta,phi,alpha设置旋转角度6.5.2、色彩控制色彩控制1.用色风格colordefColordef optionsWhitebgWhitebg(fig)Whitebg(C)6.5.2、色彩控制色彩控制2.色图colorm

42、ap两种基本着色技术RGB三元组色图矩阵和色图Colormap(CM)CM:autumn,bone,colorcube,cool,.3.浓淡处理shadingShading optionsOptions-flat,interp,facetedcolormap(R,G,B)色图设定函数matlab 的颜色数据集合为红、绿、兰三颜色矩阵R，G，B，维数m3r,g,b在0 1区间连续取值，理论上颜色种类可达无穷多种matlab使用三维向量表示一种颜色，常用颜色数据见下表 饱和色 0 0 0 黑色 0 0 1 兰色 0 1 0 绿色 0 1 1 浅兰 1 0 0 红色 1 0 1 粉红 1 1 0 黄

43、色 1 1 1 白色 调和色0.5 0.5 05 灰色0.5 0 0 暗红色1 0.62 0.4 铜色0.49 1 0.8 浅绿0.49 1 0.83 宝石兰matlab 的色图函数：hsv 饱和值色图 gray 线性灰度色图 hot 暖色色图 cool 冷色色图 bone 兰色调灰色图 copper 铜色色图 pink 粉红色图 prism 光谱色图 jet 饱和值色图II flag 红、白、蓝交替色图6.5.3 照明和材质处理1.灯光lightLight(color,option1,style,option2,position,option3)2.照明模式lightingLighting

44、optionsOptions:flat,gouraud,phong,none3.控制光效果的材质命令materialMaterial optionsshiny,dull,metal,defaultMaterial(ka kd ks n sc)-ka-背景光强度，kd-漫反射强度ks-有向反射强度,n-镜面指数，sc-反射系数4.surfl指令的精细调用格式Surfl(X,Y,Z,S,K,method)【例6.5.3.4-1】灯光、照明、材质指令所表现的图形。clf;X,Y,Z=sphere(40);colormap(jet)subplot(1,2,1);surf(X,Y,Z);shading

45、interplight(position,2,-2,2,style,local)lighting phong material(0.5,0.3,0.5,10,0.5)subplot(1,2,2);surf(X,Y,Z,-Z);shading flatlight;lighting flatlight(position,-1,-1,-2,color,y)light(position,-1,0.5,1,style,local,color,w)material(0.4,0.5,0.3,10,0.3)shading faceted 网格修饰，缺省方式shading flat 去掉黑色线条，根据小方块的值

46、确定颜色shading interp 颜色整体改变，根据小方块四角的值差补过度点的值确定颜色peaks(30);shading interp;colormap(hot)6.5.4 透明处理6.8 函数绘图的简洁命令指令名含义可执行示例ezcontour 等位线ezcontouf 画填色等位线ezmesh画网线图ezmeshc画带等位线的曲线图ezplot画二维曲线ezplot3画三维曲线ezpolar画极坐标曲线ezsurf画曲面图ezsurfc画带等位线的曲面图ezplot 符号函数的简易绘图函数ezplot的调用格式：ezplot(f)这里f为包含单个符号变量x的符号表达式，在x轴的默认范

47、围 -2*pi 2*pi内绘制f(x)的函数图ezplot(f,xmin,xmax)给定区间ezplot(f,xmin,xmax,figure(n)指定绘图窗口绘图。ezplot(sin(x)ezplot(sin(x),cos(y),-4*pi 4*pi,figure(2)Ezsurf函数应用【例6.8.2-2】使用球坐标参量画部分球壳。（图6.8-3）x=cos(s)*cos(t);y=cos(s)*sin(t);z=sin(s);ezsurf(x,y,z,0,pi/2,0,3*pi/2)view(17,40);shading interp;colormap(spring)light(pos

48、ition,0,0,-10,style,local)light(position,-1,-0.5,2,style,local)material(0.5,0.5,0.5,10,0.3)一个例子load woman;A=imread(cameraman.tif);A=double(A);B=imread(rice.png)B=double(B)Z=zeros(size(X);Y=ones(size(X);xx,yy=meshgrid(linspace(0,1,256);surface(Y,xx,yy,flipud(X),FaceColor,texturemap,edgecolor,none,cdatamapping,direct);colormap(map);view(3);xlim(0,1);surface(xx,Y,yy,flipud(A),FaceColor,texturemap,edgecolor,none,cdatamapping,direct);