matlab绘图总结(个人总结)

1、 1，单矢量绘图x=1 2 3 4 5 6 7 8 9 10;y=0 0.6 2.3 5 8.3 11.7 15 17.7 19.4 20;plot(x,y)title('简单绘图举例');xlabel('单元下标');ylabel('给定的矢量');grid2，双矢量绘图x=0:0.05:4*pi;y=sin(x);plot(x,y)3，对数坐标绘图对数坐标分为x轴对数（semilogx），y轴对数（semilogy）及双对数（loglog）三个函数，用法和plot命令相似。4，极坐标绘图极坐标绘图的用法和基本绘图命令plot（）相似，其格式为

2、polar(theta,rho)t=0:.01:2*pi;polar(t,sin(2*t).*cos(2*t) x=0:pi/50:2*pi;y1=sin(x);y2=0.6*sin(x);y3=0.3*sin(x);plot(x,y1,x,y2,x,y3) x1=0:pi/50:2*pi;x2=pi/4:pi/50:2*pi+pi/4;x3=pi/2:pi/50:2*pi+pi/2; % 问题出在哪里？！ （已解决）y1=cos(x1);y2=cos(x2);y3=cos(x3);plot(x1,y1,x2,y2,x3,y3) x(1,:)=0:pi/50:2*pi;x(2,:)=pi/4:

3、pi/50:2*pi+pi/4;x(3,:)=pi/2:pi/50:2*pi+pi/2; % 问题出在哪呢？！（已解决）y=sin(x(1,:);plot(x,y) x=0:pi/50:2*pi;y(1,:)=sin(x);y(2,:)=0.6*sin(x);y(3,:)=0.3*sin(x);plot(x,y) x(:,1)=0:pi/50:2*pi'x(:,2)=pi/4:pi/50:2*pi+pi/4'x(:,3)=pi/2:pi/50:2*pi+pi/2' % 这个又有问题了！y(:,1)=sin(x(:,1);y(:,2)=0.6*sin(x(:,1);y(:

4、,2)=0.3*sin(x(:,1);plot(x,y) 4，多组变量绘图x=0:pi/50:2*pi;y1=sin(x);y2=0.6*sin(x);y3=0.3*sin(x);plot(x,y1,x,y2,x,y3)5，双y轴绘图x=0:0.3:12;y=exp(-0.3*x).*sin(x)+0.5;plotyy(x,y,x,y,'plot','stem') t=0:900;A=1000;a=0.005;b=0.005;z1=A*exp(-a*t);z2=sin(b*t);plotyy(t,z1,t,z2,'semilogy','p

5、lot') 5,图的形式和颜色matlab 提供了四种线形和十三种标记点类型MATLAB在多组变量绘图时,可将曲线以不同的颜色,不同的线型及标记点表示出来.这些选项如下表所示: 各种颜色属性选项 'r' 红色 'm' 粉红 'g' 绿色 'c' 青色 'b' 兰色 'w' 白色 'y' 黄色 'k' 黑色 各种线型属性选项 '-' 实线 '-' 虚线 ':' 点线 '-.' 点划线 各种标记点属性

6、选项 '.' 用点号绘制各数据点 '' 用上三角绘制各数据点 '+' 用'+'号绘制各数据点 'v' 用下三角绘制各数据点 '*' 用'*'号绘制各数据点 '>' 用右三角绘制各数据点 ' .' 用'.'号绘制各数据点 '<' 用左三角绘制各数据点 's'或squar 用正方形绘制各数据点'p' 用五角星绘制各数据点'd'或diamond用菱形绘制各数据点 &#

7、39;h' 用六角星绘制各数据点 这些选项可以连在一起用,如:'-.g'表示绘制绿色的点划线,'g+'表示用绿色的'+'号绘制曲线.t=0:pi/100:2*pi;y=sin(t);y2=sin(t-0.25);y3=sin(t-0.5);plot(t,y,'m-',t,y2,'k+',t,y3,'g:') t=0:pi/20:2*pi;y=sin(t);plot(t,y,'d') 6，图形的控制与表现axis 人工选择坐标轴尺寸clf 清图形窗口ginput 利用鼠标的十字

8、准线输入hold 保持图形shg 显示图形窗口subplot 将图形窗口分成N块子窗口 t=0:pi/20:4*pi;y=sin(t);plot(t,y)grid t=0:pi/20:2*pi;x,y=meshgrid(t);subplot(2,2,1)plot(sin(t),cos(t),axis equalsubplot(2,2,2)z=sin(x)+cos(y);plot(t,z)axis(0 2*pi -2 2)subplot(2,2,3)z=sin(x).*cos(y);plot(t,z)axis(0 2*pi -1 1)subplot(2,2,4)z=sin(x).2-cos(y)

9、.2;plot(t,z),axis(0 2*pi -1 1) x,y,z=peaks;contour(x,y,z,20,'k')hold onpcolor(x,y,z)shading interphold off 坐标轴控制指令axis(xmin xmax ymin ymax)x=0:.025:pi/2;plot(x,tan(x),'-ro')axis(0 pi/2 0 5) 图形的标注 title 标题xlabel x轴标注ylabel y轴标注text 任意定位的标注gtext 鼠标定位标注legend 标注图例加注坐标轴标题和图形标识，在图形中加注文本t=

10、0:pi/100:2*pi;y=sin(t);plot(t,y)axis(0 2*pi -1 1)xlabel('0leqittrmleqpi','FontSize',16)ylabel('sin(t)','FontSize',20)title('正弦函数图形','FontName','隶书','FontSize',20)text(3*pi/4,sin(3*pi/4),'leftarrowsin(t)=0.707','FontSize'

11、,16)text(pi,sin(pi),'leftarrowsin(t)=0','FontSize',16)text(5*pi/4,sin(5*pi/4),'sin(t)=-0.707rightarrow','FontSize',16,'HorizontalAlignment','right')指定 TeX字符t=0:pi/100:2*pi;alpha=-0.8;beta=15;y=sin(beta*t).*exp(alpha*t);plot(t,y)title('itAe-alphaitt

12、sinbetaittalpha<<beta')xlabel('时间 us.'),ylabel('幅值') 在图形中添加图例框x=0:pi/100:2*pi;y1=sin(x);y2=0.6*sin(x);y3=0.3*sin(x);plot(x,y1,'-ro',x,y2,'-g*',x,y3,'bd')legend('曲线1','曲线2','曲线3') 特殊图形1，条形图（1）二维条形图函数 意义bar 二维垂直条形图bath 二维水平条形图ba

13、r3 三维垂直条形图bar3h 三维水平条形图 y=27 38 55 37 98 45 67 43 54;bar(y)bar(y,'stack')barh(y)y=9 8 6;2 5 8;6 2 9;5 8 7;9 4 2;bar(y)bar(y,'stack')barh(y)bar3(y)bar3h(y)bar(y,'stack')barh(y,'stack') 指定x坐标的条形图x=1 2 4 7 10;y=9 8 6;2 5 8;6 2 9;5 8 7;9 4 2;bar(x,y)figure,barh(x,y) 可以把二维

14、曲线当做二维条形线的包络线 x=0:pi/10:2*pi;y=sin(x);bar(x,y)bar(x,y,'r') y=9 8 6;2 5 8;6 2 9;5 8 7;9 4 2;bar3(y,'group') % 分组形式的三维图figure,bar3(y) % 分列形式的三维图 有两组生物医学的实验数据，一组表示物质成分（TCE），一组表示温度（temp），数据是在35天里每隔5天的取样，将物质成分和温度与时间的关系画在同一张图中TCE=515 429 370 250 135 120 60 120;temp=29 23 27 25 20 23 23 27;

15、days=0:5:35;bar(days,temp,'c')xlabel('Day')ylabel('Temperature(oC)')h1=gca; % 获得当前轴对象句柄h2=axes('Position',get(h1,'position') ); % 建立新的轴对象句柄plot(days,TCE,'LineWidth',3)gridset(h2,'YAxisLocation','right','Color','none',

16、9;XTickLabel',)set(h2,'XLim',get(h1,'XLim'),'Layer','top')set(h2,'XLim',get(h1,'XLim'),'Layer','top') %完成一张图形的绘制text(11,380,'Concentration','Rotation',-57,'FontSize',16)ylabel('TCE Concentration (PPM)'

17、;)title('Bioremediation','FontSize',16) 饼图1，不分离的饼图x=5.5 74.7 44.5 33.2 46.6;pie(x)2,带分离的饼图x=5.5 74.7 44.5 33.2 46.6;pie(x,0 1 0 0 1)3，不完整的饼图当x的元素之和不为1时绘制的是不完整的饼图x=0.2 0.3 0.4;pie(x)4，三维饼图pie3(1 2 3 4 5, 0 0 0 0 1)legend('矿大','地质大学','石油大学','海洋大学','其

18、他') clear,clc,close all; x=0:0.01:2; y=-(x-1).2; h1=plot(x,y); hold on; h2=plot(1,0,'o','MarkerSize',15, 'MarkerEdgecolor','k', 'MarkerFacecolor','g'); legend(h1,h2,'-(x-1)2','max(fx)'); 其他图形 1,直方图yn=randn(10000,1);hist(yn)figure,hi

19、st(yn,20) y=randn(10000,3);hist(y)2,用杆状图表现离散数据x=0:0.2:10;y=exp(-0.3*x).*sin(x);stem(x,y)figure,stem(x,y,':sr')3,阶梯图的表现方法alpha=0.01;beta=0.5;t=0:10;f=exp(-alpha*t).*sin(beta*t);stairs(t,f)hold onplot(t,f,':*')hold offlabel='函数 e-(alpha*t)sinbeta*t的阶梯图'text(0.5,-0.2,label,'

20、FontSize',14)xlabel('t=0:10','FontSize',14)axis(0 10 -1.2 1.2) 4，彩色分点图t=0:pi/10:2*pi;y=sin(t);scatter(t,y)figure,scatter(t,y,'v')figure,scatter(t,y,(abs(y)+2).4,'filled')figure,scatter(t,y,30,y,'v','filled')figure,scatter(t,y,(t+1).3,t,'filled&

21、#39;) 三维图形1，三维曲线图t=0:pi/50:10*pi;x=exp(-t/15).*sin(2*t);y=exp(-t/15).*cos(2*t);z=t;plot3(x,y,t)axis square;grid on当x,y和z都是矩阵时的三维曲线X,Y=meshgrid(-2:0.1:2);Z=X.*exp(-X.2-Y.2);plot3(X,Y,Z)grid on t=0:pi/20:5*pi;plot3(sin(t),cos(t),t,':*r') x,y,z=peaks(10);surf(x,y,z)2,三维网格图mesh(peaks(20) meshc(p

22、eaks(20) Waterfall(peaks(30) p=peaks(30);meshz(p)3，着色表面图surf 绘制着色表面图的基本命令 surfc 绘制带等高线的着色表面图surfl 绘制表面图的光照效应 x=-1.5:0.3:1.5;y=-1:0.2:1;x,y=meshgrid(x,y);z=sqrt(4-x.2/9-y.2/4);surf(x,y,z)surfc(x,y,z)surfl(x,y,z)colormap(hot)surfl(x,y,z)shading facetedshading flat surfl(x,y,z)shading interp surfl(x,y,

23、z)shading flat t=0:pi/12:3*pi;r=abs(exp(-0.25*t).*sin(t);x,y,z=cylinder(r,30);surfl(x,y,z)colormap(pink)shading interp 表面图形的透明处理surfl(x,y,z)hidden off surfl(x,y,z)hidden on x,y=meshgrid(-8:.5:8);R=sqrt(x.2+y.2)+eps;z=sin(R)./R;mesh(z)colormap(0 0 1)hidden off 表面图形的颜色映射 x,y=meshgrid(-8:.5:8);R=sqrt(x

24、.2+y.2)+eps;z=sin(R)./R;mesh(z)colormap(0 0 1) % RGB三基色 x,y=meshgrid(-8:.5:8);R=sqrt(x.2+y.2)+eps;z=sin(R)./R;mesh(z)colormap(copper)colormap(bone)colormap(gray)colormap(hot)colormap(pink)colormap(cool)colormap(hsv)colormap(jet)colormap(spring)colormap(summer)colormap(autumn)colormap(winter)colormap

25、(flag)colormap(prism)colormap(copper)colormap(colorcube)colormap(lines) 切片图slice(X,Y,Z,v,xi,yi,zi) % 其中X,Y,Z为使用meshgrid生成的三维网格坐标矩阵，v为所绘制图形的函数，xi，yi，zi为切片位置x=-2:0.1:2;y=-2:0.25;z=-2:0.25:2;X,Y,Z=meshgrid(x,y,z);V=X.*exp(-X.2-Y.2-Z.2);xi=-0.7,0.7;yi=0.5;zi=-0.5;slice(X,Y,Z,V,xi,yi,zi)xlabel('x'

26、;);ylabel('y');zlabel('z');colorbar('horiz')view(-30,45) 图形的高级处理技术1，视图 控制视图和旋转的函数包括view，viewmtx，rotate，和rotate3d等1，视图方向view(az,el)peaks(50)view(90,0)view(30,45)view(2) peaks(50)rotate3d 2,光照（1），light 创建光源light函数为光源创建函数，其调用格式为：light('color',option1,'style',opti

27、on2,'poisition',option3)在执行该命令函数之前，图形采用的是强度相等的漫射光，一旦执行该命令，则激活图形上各子对象所有与光有关的属性，其参数的缺省设置为：无穷远处的白光光源穿过1 0 1射向坐标原点 （2），light 光照模式light options 设置照明模式 flat 入射光均匀照射在对象的每个小面上，该选项为缺省模式 phong 先对顶点颜色插值，在计算各像素点的反光，该选项表现效果最好，但占用机时较多 gouraud 对顶点处法线插值，再对由顶点勾画的面进行插值，选项用于曲面的表现（3），material 材料反射系数material opt

28、ions 使用预定的材料反射系数shiny 有光泽 明亮 镜反射份额较大dull 暗淡 柔和 漫反射份额较大 无镜面亮点metal 有金属光泽 镜反射份额大 背景光和漫反射份额小default 缺省方式（4）material(ka,kd,ks,n,sc)对五大反射要素进行特别设置ka 环境光kd 漫反射强度ks 镜射光强度n 镜面指数sc 镜面颜色的反射系数 clft=linspace(0,2*pi,100);r=1-exp(-t/2).*cos(4*t);%旋转母线x,y,z=cylinder(r,60); % 产生旋转柱面数据ii=find(x<0&y<0); % 确定

29、x-y平面上第四象限上的数据下标z(ii)=NaN; %剪切surf(x,y,z)colormap(pink)shading interplight('position',-3 -1 3,'style','local') % 设置光源material(0.5 0.4 0.3 10 0.3) % 设置表面反射 显示效果增强后的二元函数表面图形图及光照x,y=meshgrid(-8:.1:8);r=sqrt(x.2+y.2)+eps;z=2*sin(r)./r;surf(x,y,z)shading interp;daspect(5 5 11)axis

30、 tightview(-50,30)light('position',-10 -10 2) 绘制双光照射的球sphere(36);axis equallight('Position',1 3 2);light('Position',-3 -1 3);material shinyaxis offrotate3d x,map=imread('C:UsersuserDesktop新建文件夹单单.jpg'); % 读取图形文件image(x);colormap(map) % 显示图像axis image off % 关掉坐标轴 3,句柄图

31、形（1），自定义坐标轴的显示刻度t=0:pi/20:2*pi;y=sin(t);plot(t,y) t=-pi:pi/20:pi;y=sin(t);plot(t,y)set(gca,'xtick',-pi:pi/2:pi,'xticklabel','-pi','-pi/2','0','pi/2','pi') % 此句还可以写为：set(gca,'XtickLabeks','-pi|-pi/2|0|pi/2|pi'),如果是数字也可以直接写成set(gc

32、a,'XtickLabeks',1;5;10;20)axis (-pi pi -1 1)grid 2，在图形的上层显示网格建立句柄为h的三维曲面h=surf(sphere (30) figure,h=surf(sphere (30)set(h,'EdgeColor',0.8 0.8 0.8) rotate3d figure,h=surf(sphere (30)set(h,'EdgeColor','g') rotate3d 将坐标轴显示在图形的上层 pcolor(peaks);shading interpset(gca,'L

33、ayer','top')grid 用鼠标确定在曲线上画点的位置t=0:pi/20:2*pi;plot(t,cos(t) hold onset(gcf,'WindowButtonDownFcn','p=get(gca,''CurrentPoint'');','px=p(1,1);py=cos(px);','plot(t,cos(t);','plot(px,py,''*r'');')axis(axis) 绘制三个Y轴的曲线ax,hl

34、ines=plotyyy(x1,y1,x2,y2,x3,y3,ylabels)其中x1,y1;x2,y2和x3,y3分别为三组绘图数据，ylabels为3*1阶的单元数组，用于存储三个Y轴的标识，输出参数ax,hlines分别为三个轴和三条曲线的图形句柄。该函数是由加拿大人Denis Gilbert编写的function ax,hlines=plotyyy(x1,y1,x2,y2,x3,y3,ylabels)if nargin=6 ylabels1='' ylabels2='' ylabels3=''elseif nargin>7 erro

35、r('Too many input arguments')elseif nargin<6 error('Not enough input arguments')figure('units','normalized','DefaultAxesXMinorTicck','on','DefaultAxesYMinorTicck','on');ax,hlines(1),hlines(2)=plotyy(x1,y1,x2,y2);cfig=get(gcf,'col

36、or');pos=0.1 0.1 0.7 0.8;offset=pos(3)/5.5; % 使用plotyy（）函数绘制前两条曲线 pos(3)=pos(3)-offset/2;set(ax,'position',pos);% 改变已绘制曲线的比例pos3=pos(1) pos(2) pos(3)+offset pos(4);limx1=get(ax(1),'xlim');limx3=limx1(1) limx1(1)+1.2*(limx1(2)-limx1(1); % 确定第三个轴的位置ax(3)=axes('Position',pos

37、3,'box','off','Color','none','XColor','k','YColor','r','xtick','xlim',limx3,'yaxislocation','right');% 绘制第三个轴hlines(3)=line(x3,y3,'color','r','Parent',ax(3);limy3=get(ax(3),'YL

38、im');% 绘制第三条曲线line(limx1(2) limx3(2),limy3(1) limy3(1),'Color',cfig,'Parent',ax(3),'Clipping','off');axes(ax(2)% 在第二，第三轴之间消隐x轴set(get(ax(1),'ylabel'),'string',ylabels1)set(get(ax(2),'ylabel'),'string',ylabels2)set(get(ax(3),'yla

39、bel'),'string',ylabels3) x1=0:pi/100:10;y1=sin(x1);x2=0:200;y2=log(x2);x3=1:100;y3=5.*x3;ylabels=1 1 1'ax,hlines=plotyyy(x1,y1,x2,y2,x3,y3) 狐狸的存在对兔子自然增长率的影响homer=1;d=0.5;a=0.1;b=0.01;R(1)=25;F(1)=5;years=20;h=0.1;k=1; for years=0:h:years k=k+1; ii=k-1; RK1=r*R(ii)-a*R(ii)*F(ii); RK2=

40、r*(R(ii)+h/2)-a*(R(ii)+h/2)*(F(ii)+h/2*RK1); RK3=r*(R(ii)+h/2)-a*(R(ii)+h/2)*(F(ii)+h/2*RK2); RK4=r*(R(ii)+h)-a*(R(ii)+h)*(F(ii)+h*RK3); R(k)=R(ii)+h/6*(RK1+2*RK2+2*RK3+RK4);FK1=-d*F(ii)+b*R(ii)*F(ii);FK2=-d*(F(ii)+h/2*FK1)+b*(R(ii)+h/2)*(F(ii)+h/2*FK1);FK3=-d*(F(ii)+h/2*FK2)+b*(R(ii)+h/2)*(F(ii)+h/

41、2*FK2);FK4=-d*(F(ii)+h*FK3)+b*(R(ii)+h)*(F(ii)+h*FK1);F(k)=F(ii)+h/6*(FK1+2*FK2+2*FK3+FK4);endplot(R,F)xlabel('Rabbit');ylabel('Fox');zlabel('Years')grid onfigurehold onplot(R,'-r')plot(F,'-b')grid ontitle('Red line stands for the Rabbits and Bule lines st

42、ands for the Fox')xlabel('Years');ylabel('Animal Number');hold off 计算和绘制serpinski分形曲线x0=0 4 8;y0=0 4*sqrt(3) 0;g=1;fill(x0(1) x0(2) x0(3),y0(1) y0(2) y0(3),'g');hold onfor m=1:4 for I=1:g j=3*(I-1); k=9*(I-1); ax=(x0(j+1)+x0(j+2)/2; bx=(x0(j+2)+x0(j+3)/2; cx=(x0(j+1)+x0(j+3)/2; x(k+1)=x0(j+1);x(k+2)=ax;x(k+3)=cx;x(k+4)=x0(j+2);x(k+5)=ax;x(k+6)=bx;x(k+7)=x0(j+3);x(k+8)=bx;x(k+9)=cx; ay=(y0(j+1)+y0(j+2)/2; by=(y0(j+2)+y0(j+3)/2; cy=(y0(j+1)+y0(j+3)/2; y(k+1)=y0(j+1);y(k+2)=ay;y(k+3)