《MATLAB数据和函数的可视化.ppt》由会员分享，可在线阅读，更多相关《MATLAB数据和函数的可视化.ppt（106页珍藏版）》请在三一办公上搜索。
1、第6章 数据和函数的可视化,matlab语言的绘图功能,不仅能绘制几乎所有的标准图形，而且其表现形式也是丰富多样的。matlab语言不仅具有高层绘图能力，而且还具有底层绘图能力句柄绘图方法。在面向对象的图形设计基础上，使得用户可以用来开发各专业的专用图形。,6.1 引导,6.1.1 离散数据和离散函数的可视化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
2、plot命令基本调用格式,（1）plot 最基本的二维图形指令plot的功能：plot命令自动打开一个图形窗口Figure 用直线连接相邻两数据点来绘制图形根据图形坐标大小自动缩扩坐标轴，将数据标尺及单位标注自动加到两个坐标轴上，可自定坐标轴，可把x,y 轴用对数坐标表示,plot的调用格式,plot(x)缺省自变量绘图格式，x为向量,以x元素值为纵坐标，以相应元素下标为横坐标绘图 plot(x,y)基本格式，以y(x)的函数关系作出直角坐标图，如果y为nm的矩阵，则以x 为自变量，作出m条曲线plot(x1,y1,x2,y2)多条曲线绘图格式,plot(x,y,s)开关格式，开关量字符串s设
3、定曲线颜色和绘图方式，使用颜色字符串的前13个字母，如 yellowyel表示等。或plot(x1,y1,s1,x2,y2,s2，),6.2.2 曲线的色彩、线型和数据点型,1.色彩和线型色彩：符号：b,g,r,c,m,y,k,w 含义：蓝 绿 红 青 品 黄 黑 白线型：符号：-,:,-.,-含义：实线 虚线 点划 双划,6.2.2 曲线的色彩、线型和数据点型,2.数据点型,【例6.2.2.2-1】用图形演示平面上一个方块四个顶点在仿射投影（Affine Projection）下的位置、形状变化。,分别实现矩阵平移、旋转和刻度变换,%平面上四个点和它们构成的方块p1=-0.5,0,1;p2=
4、-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,dy;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
5、=TRSi*Sq;%进行仿射变换subplot(2,2,i)for k=1:4plot(W(1,k),W(2,k),ssk);%用不同点形标志四个顶点axis(-3,3,-1,5),axis equalhold onendplot(W(1,:),W(2,:)%连接四个顶点使得图形封闭grid on%title(tti)%hold off%end,6.2.3 坐标、刻度和分格线控制,1、坐标控制2、刻度、分格线和坐标框,坐标控制,【例6.2.3.1-1】观察各种轴控制指令的影响。演示采用长轴为3.25，短轴为1.15 的椭圆。注意：采用多子图表现时，图形形状不仅受“控制指令”影响，而且受整个图面
6、“宽高比”及“子图数目”的影响。本书这样处理，是出于篇幅考虑。读者欲想准确体会控制指令的影响，请在全图状态下进行观察。（图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)su
7、bplot(2,3,4),plot(x,y),axis image,box off,title(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)坐标框
8、boxbox on(off),【例6.2.3.2-1】通过绘制二阶系统阶跃响应，综合演示图形标识。本例比较综合，涉及的指令较广。请读者耐心读、实际做、再看例后说明，定会有匪浅收益。（图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,max(y),grid on%title(i
9、t 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);text(ts,0.85,cell_string
10、)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,xtick,-pi-pi/2 0 pi/2 pi
11、);set(gca,xticklabel,0);axis(-pi pi,-1 1);,6.2.4 图形标识,1、简捷指令形式titlexlabel,ylabel,zlabellegend(s1,s2,)text(xt,yt,s),2.精细指令形式(1)允许标识多行字符,(2)允许标识字体,风格及大小进行设置,(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
12、);y2=sin(t+0.5);plot(t,y)自动出现第一个窗口figure(2)plot(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 square,3.双纵坐标图,plotyy(X1,Y1,X2,Y2)plotyy(X1
13、,Y1,X2,Y2,FUN)plotyy(X1,Y1,X2,Y2,FUN1,FUN2),4.多子图,subplot 子图分割命令调用格式：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;zoom,6.3三维绘图的基本操作,6.3.1 三维线图指令plot3
14、6.3.2三维网线图和曲面图6.3.3透视、镂空和裁切,6.3.1 三维线图指令plot3,t=(0:0.02:2)*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(链,宝石),三维网线图和曲面图,用曲面图表现函数。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
15、与y 坐标为元素的下标 mesh(x,y,z)x,y,z分别为三维空 间的坐标位置,三维网线图作图要领,生成坐标 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 三维曲面绘图函数，与网格图 看起来一样与三维网
16、线图的区别：网线图：线条有颜色，空挡是黑色的（无颜色）曲面图：线条是黑色的，空挡有颜色（把线条之间的空挡填充颜色，沿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),透视、镂空和裁切,透视演示X0,Y0,Z0=sphere(30);X=2*X0;
17、Y=2*Y0;Z=2*Z0;surf(X0,Y0,Z0);shading interphold 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;
18、X,Y=meshgrid(x,y);ZZ=X.2-Y.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.1特殊图形指令例示,1.面域图area2.各种直方图bar,barh,bar3,bar3h3.饼图pie,pie34.填色
19、图fill,fill35.射线图compass和羽毛图feather6.Vornonoi图和三角剖分7.彩带图ribbon8.离散杆图stem,stem39.二维半图指令pcolor，contour，contourf10.散点图scatter，scatter3，plotmatrix11.不规则数据的网线图和曲面图trimesh，trisurf,1.面域图area,clf;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.各
20、种直方图bar,barh,bar3,bar3h,二维直方图有两种图型：垂直直方图和水平直方图。而每种图型又有两种表现模式：累计式：分组式。本例选其两种加以表现。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,pie3 用来表示各元素占总和的百分
21、数。该指令第二输入宗量为与第一宗量同长的0-1向量，1使对应扇块突出。a=1,1.6,1.2,0.8,2.1;subplot(1,2,1),pie(a,1 0 1 0 0);legend(1,2,3,4,5)subplot(1,2,2),pie3(a,a=min(a),colormap(cool),4.填色图fill,fill3,MATLAB画任意多边形的一种方法；保证绘图数据首尾重合，使勾画多边形封闭；使用图柄对图形的属性进行精细设置。clf;n=10;dt=2*pi/n;t=0:dt:2*pi;t=t,t(1);x=sin(t);y=cos(t);fill(x,y,c);axis offh
22、t=text(0,0,fontname隶书fontsize32十边形);set(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;rand(state,111)n=30;A
23、=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,contourf,伪彩图pcolor；等位线指令co
24、ntour、contourf；等位线标高指令 clabel 的配合使用和区别。（1）本例等位线指令中的第4输入宗量n设定高度的等级数，第5输入宗量设定等位线的线型、色彩。（2）左右两图的标高方法不同。左图的标识以“+”引导，水平放置。右图沿线布置。这是由clabel的调用格式不同产生的。（3）左右两图色彩的形成方法不同，色彩效果也不同。（4）在左图中，colorbar画出一根垂直色标尺，而caxis决定该色标尺的刻度,等高线修饰二维contour(Z,n)-绘制n条等高线C=contourc(Z,n)-计算n条等高线的坐标Clable(c)-给等高线加标注例、在二维平面上绘制peaks函数的1
25、0条等 高线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)colormap 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
26、.散点图scatter,scatter3,plotmatrix,x=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(del2(Z);meshz(X,Y,Z,C)；hold on,scatter3(X(:),Y(:),Z(:),filled)；hold off,colormap(hot)；,指令plotmatrix有两种基本调用方式：（1）对于数据矩阵维的X（p*n）和维的Y（p*m），调用格式plotmatrix(X,Y)将画出一个分割成（m*n）个子散点图。其中第（i，j）个子散点图是根据
27、Y第i列和X第j列数据画出的。（2）对于数据矩阵维的X，调用格式plotmatrix(X)将画出分割成（n*n）个子块的图。该图的对角块，画出的是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)subplot(1,3,3),plotmatrix(X,Y),用三角网线、曲面图表现函数。rand(seed,22),X=6*pi*(rand(20,10)-0.5);Y=6*pi*(r
28、and(20,10)-0.5);R=sqrt(X.2+Y.2)+eps;Z=sin(R)./R;tri=delaunay(X,Y);subplot(1,2,1),trimesh(tri,X,Y,Z)subplot(1,2,2),trisurf(tri,X,Y,Z)colormap(jet);brighten(0.5),6.4.2四维表现,1.用色彩表现函数特征2.切片图和切片等位线图3.,1.用色彩表现函数特征,surf（x,y,z）,surf(x,y,z,s)用色图阵表现函数的不同特征。演示：当三维网线图、曲面图的第四个输入宗量取一些特殊矩阵时，色彩就能表现或加强函数的某特征，如梯度、曲率、
29、方向导数等。,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,Y,
30、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.色图的变化3.影片动
31、画,彗星状轨迹图,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);plot3(x,y,
32、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.色图的变化,peaksspinmap,3.动画效果,动画生成的步骤创建帧矩阵 moviein对动画中的每一帧生成图形，并把它们放到帧矩阵中 getframe从帧矩阵中回放动画movie,clf
33、shg,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.1、视点控制和图形的旋动,1.视点控制viewView(az,
34、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.色图colormap两种基本着色技术RGB三元组色图矩阵和色图Colormap(CM)CM:autumn,bone,col
35、orcube,cool,.3.浓淡处理shadingShading optionsOptions-flat,interp,faceted,colormap(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 黄色 1 1 1 白色,调和色0.5 0.5 05 灰色0.5 0 0 暗红色1 0.62 0.4 铜色0
36、.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 optionsOptions:flat,gouraud,phong,none3.控制光效果的材质命
37、令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 interplight(position,2,-2,2,style,local)lightin
38、g phongmaterial(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 去掉黑色线条，根据小方块的值确定颜色,shading interp 颜色整体改变，根据小方块四角的值差补过度点的值确定颜
39、色,peaks(30);shading interp;colormap(hot),6.5.4 透明处理,6.8 函数绘图的简洁命令,ezplot 符号函数的简易绘图函数,ezplot的调用格式：ezplot(f)这里f为包含单个符号变量x的符号表达式，在x轴的默认范围-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】使用球坐标参量画
40、部分球壳。（图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(position,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);,
