matlab編程100例之圖形應(yīng)用篇

圖形應(yīng)用篇實(shí)例1：三角函數(shù)曲線（1）functionshili01h0=figure('toolbar','none',...'position',[19856350300],...'name','實(shí)例01');h1=axes('parent',h0,...'visible','off');x=-pi:0.05:pi;y=sin(x);plot(x,y);xlabel('自變量X');ylabel('函數(shù)值Y');title('SIN()函數(shù)曲線');gridon實(shí)例2：三角函數(shù)曲線（2）functionshili02h0=figure('toolbar','none',...'position',[200150450350],...'name','實(shí)例02');x=-pi:0.05:pi;y=sin(x)+cos(x);plot(x,y,'-*r','linewidth',1);gridonxlabel('自變量X');ylabel('函數(shù)值Y');title('三角函數(shù)');實(shí)例3：圖形的疊加functionshili03h0=figure('toolbar','none',...'position',[200150450350],...'name','實(shí)例03');x=-pi:0.05:pi;y1=sin(x);y2=cos(x);plot(x,y1,...'-*r',...x,y2,...'--og');gridonxlabel('自變量X');ylabel('函數(shù)值Y');title('三角函數(shù)');實(shí)例4：雙y軸圖形的繪制functionshili04h0=figure('toolbar','none',...'position',[200150450250],...'name','實(shí)例04');x=0:900;a=1000;b=0.005;y1=2*x;y2=cos(b*x);[haxes,hline1,hline2]=plotyy(x,y1,x,y2,'semilogy','plot');axes(haxes(1))ylabel('semilogplot');axes(haxes(2))ylabel('linearplot');實(shí)例5：?jiǎn)蝹€(gè)軸窗口顯示多個(gè)圖形functionshili05h0=figure('toolbar','none',...'position',[200150450250],...'name','實(shí)例05');t=0:pi/10:2*pi;[x,y]=meshgrid(t);subplot(2,2,1)plot(sin(t),cos(t))axisequalsubplot(2,2,2)z=sin(x)-cos(y);plot(t,z)axis([02*pi-22])subplot(2,2,3)h=sin(x)+cos(y);plot(t,h)axis([02*pi-22])subplot(2,2,4)g=(sin(x).^2)-(cos(y).^2);plot(t,g)axis([02*pi-11])實(shí)例6：圖形標(biāo)注functionshili06h0=figure('toolbar','none',...'position',[200150450400],...'name','實(shí)例06');t=0:pi/10:2*pi;h=plot(t,sin(t));xlabel('t=0到2\pi','fontsize',16);ylabel('sin(t)','fontsize',16);title('\it{從0to2\pi的正弦曲線}','fontsize',16)x=get(h,'xdata');y=get(h,'ydata');imin=find(min(y)==y);imax=find(max(y)==y);text(x(imin),y(imin),...['\leftarrow最小值=',num2str(y(imin))],...'fontsize',16)text(x(imax),y(imax),...['\leftarrow最大值=',num2str(y(imax))],...'fontsize',16)實(shí)例7：條形圖形functionshili07h0=figure('toolbar','none',...'position',[200150450350],...'name','實(shí)例07');tiao1=[56254822454541445745512];tiao2=[4748575854526548];t=0:7;bar(t,tiao1)xlabel('X軸');ylabel('TIAO1值');h1=gca;h2=axes('position',get(h1,'position'));plot(t,tiao2,'linewidth',3)set(h2,'yaxislocation','right','color','none','xticklabel',[])實(shí)例8：區(qū)域圖形functionshili08h0=figure('toolbar','none',...'position',[200150450250],...'name','實(shí)例08');x=91:95;profits1=[8875849377];profits2=[5164545668];profits3=[4254342524];profits4=[263818154];area(x,profits1,'facecolor',[0.50.90.6],...'edgecolor','b',...'linewidth',3)holdonarea(x,profits2,'facecolor',[0.90.850.7],...'edgecolor','y',...'linewidth',3)holdonarea(x,profits3,'facecolor',[0.30.60.7],...'edgecolor','r',...'linewidth',3)holdonarea(x,profits4,'facecolor',[0.60.50.9],...'edgecolor','m',...'linewidth',3)holdoffset(gca,'xtick',[91:95])set(gca,'layer','top')gtext('\leftarrow第一季度銷量')gtext('\leftarrow第二季度銷量')gtext('\leftarrow第三季度銷量')gtext('\leftarrow第四季度銷量')xlabel('年','fontsize',16);ylabel('銷售量','fontsize',16);實(shí)例9：餅圖的繪制functionshili09h0=figure('toolbar','none',...'position',[200150450250],...'name','實(shí)例09');t=[542135;685435;452512;486845;685469];x=sum(t);h=pie(x);textobjs=findobj(h,'type','text');str1=get(textobjs,{'string'});val1=get(textobjs,{'extent'});oldext=cat(1,val1{:});names={'商品一：';'商品二：';'商品三：'};str2=strcat(names,str1);set(textobjs,{'string'},str2)val2=get(textobjs,{'extent'});newext=cat(1,val2{:});offset=sign(oldext(:,1)).*(newext(:,3)-oldext(:,3))/2;pos=get(textobjs,{'position'});textpos=cat(1,pos{:});textpos(:,1)=textpos(:,1)+offset;set(textobjs,{'position'},num2cell(textpos,[3,2]))實(shí)例10：階梯圖functionshili10h0=figure('toolbar','none',...'position',[200150450400],...'name','實(shí)例10');a=0.01;b=0.5;t=0:10;f=exp(-a*t).*sin(b*t);stairs(t,f)holdonplot(t,f,':*')holdoffglabel='函數(shù)e^{-(\alpha*t)}sin\beta*t的階梯圖';gtext(glabel,'fontsize',16)xlabel('t=0:10','fontsize',16)axis([010-1.21.2])實(shí)例11：枝干圖functionshili11h0=figure('toolbar','none',...'position',[200150450350],...'name','實(shí)例11');x=0:pi/20:2*pi;y1=sin(x);y2=cos(x);h1=stem(x,y1+y2);holdonh2=plot(x,y1,'^r',x,y2,'*g');holdoffh3=[h1(1);h2];legend(h3,'y1+y2','y1=sin(x)','y2=cos(x)')xlabel('自變量X');ylabel('函數(shù)值Y');title('正弦函數(shù)與余弦函數(shù)的線性組合');實(shí)例12：羅盤圖functionshili12h0=figure('toolbar','none',...'position',[200150450250],...'name','實(shí)例12');winddirection=[54246584256122356212532434254];windpower=[255368127614108];rdirection=winddirection*pi/180;[x,y]=pol2cart(rdirection,windpower);compass(x,y);desc={'風(fēng)向和風(fēng)力','北京氣象臺(tái)','10月1日0：00到','10月1日12：00'};gtext(desc)實(shí)例13：輪廓圖functionshili13h0=figure('toolbar','none',...'position',[200150450250],...'name','實(shí)例13');[th,r]=meshgrid((0:10:360)*pi/180,0:0.05:1);[x,y]=pol2cart(th,r);z=x+i*y;f=(z.^4-1).^(0.25);contour(x,y,abs(f),20)axisequalxlabel('實(shí)部','fontsize',16);ylabel('虛部','fontsize',16);h=polar([02*pi],[01]);delete(h)holdoncontour(x,y,abs(f),20)實(shí)例14：交互式圖形functionshili14h0=figure('toolbar','none',...'position',[200150450250],...'name','實(shí)例14');axis([010010]);holdonx=[];y=[];n=0;disp('單擊鼠標(biāo)左鍵點(diǎn)取需要的點(diǎn)');disp('單擊鼠標(biāo)右鍵點(diǎn)取最后一個(gè)點(diǎn)');but=1;whilebut==1[xi,yi,but]=ginput(1);plot(xi,yi,'bo')n=n+1;disp('單擊鼠標(biāo)左鍵點(diǎn)取下一個(gè)點(diǎn)');x(n,1)=xi;y(n,1)=yi;endt=1:n;ts=1:0.1:n;xs=spline(t,x,ts);ys=spline(t,y,ts);plot(xs,ys,'r-');holdoff實(shí)例14：交互式圖形functionshili14h0=figure('toolbar','none',...'position',[200150450250],...'name','實(shí)例14');axis([010010]);holdonx=[];y=[];n=0;disp('單擊鼠標(biāo)左鍵點(diǎn)取需要的點(diǎn)');disp('單擊鼠標(biāo)右鍵點(diǎn)取最后一個(gè)點(diǎn)');but=1;whilebut==1[xi,yi,but]=ginput(1);plot(xi,yi,'bo')n=n+1;disp('單擊鼠標(biāo)左鍵點(diǎn)取下一個(gè)點(diǎn)');x(n,1)=xi;y(n,1)=yi;endt=1:n;ts=1:0.1:n;xs=spline(t,x,ts);ys=spline(t,y,ts);plot(xs,ys,'r-');holdoff實(shí)例15：變換的傅立葉函數(shù)曲線functionshili15h0=figure('toolbar','none',...'position',[200150450250],...'name','實(shí)例15');axisequalm=moviein(20,gcf);set(gca,'nextplot','replacechildren')h=uicontrol('style','slider','position',...[1001050020],'min',1,'max',20)forj=1:20plot(fft(eye(j+16)))set(h,'value',j)m(:,j)=getframe(gcf);endclf;axes('position',[0011]);movie(m,30)實(shí)例16：勞倫茲非線形方程的無(wú)序活動(dòng)functionshili15h0=figure('toolbar','none',...'position',[200150450250],...'name','實(shí)例15');axisequalm=moviein(20,gcf);set(gca,'nextplot','replacechildren')h=uicontrol('style','slider','position',...[1001050020],'min',1,'max',20)forj=1:20plot(fft(eye(j+16)))set(h,'value',j)m(:,j)=getframe(gcf);endclf;axes('position',[0011]);movie(m,30)實(shí)例17：填充圖functionshili17h0=figure('toolbar','none',...'position',[200150450250],...'name','實(shí)例17');t=(1:2:15)*pi/8;x=sin(t);y=cos(t);fill(x,y,'r')axissquareofftext(0,0,'STOP',...'color',[111],...'fontsize',50,...'horizontalalignment','center')例18：條形圖和階梯形圖functionshili18h0=figure('toolbar','none',...'position',[200150450250],...'name','實(shí)例18');subplot(2,2,1)x=-3:0.2:3;y=exp(-x.*x);bar(x,y)title('2-DBarChart')subplot(2,2,2)x=-3:0.2:3;y=exp(-x.*x);bar3(x,y,'r')title('3-DBarChart')subplot(2,2,3)x=-3:0.2:3;y=exp(-x.*x);stairs(x,y)title('StairChart')subplot(2,2,4)x=-3:0.2:3;y=exp(-x.*x);barh(x,y)title('HorizontalBarChart')實(shí)例19：三維曲線圖functionshili19h0=figure('toolbar','none',...'position',[200150450400],...'name','實(shí)例19');subplot(2,1,1)x=linspace(0,2*pi);y1=sin(x);y2=cos(x);y3=sin(x)+cos(x);z1=zeros(size(x));z2=0.5*z1;z3=z1;plot3(x,y1,z1,x,y2,z2,x,y3,z3)gridonxlabel('X軸');ylabel('Y軸');zlabel('Z軸');title('Figure1:3-DPlot')subplot(2,1,2)x=linspace(0,2*pi);y1=sin(x);y2=cos(x);y3=sin(x)+cos(x);z1=zeros(size(x));z2=0.5*z1;z3=z1;plot3(x,z1,y1,x,z2,y2,x,z3,y3)gridonxlabel('X軸');ylabel('Y軸');zlabel('Z軸');title('Figure2:3-DPlot')實(shí)例20：圖形的隱藏屬性functionshili20h0=figure('toolbar','none',...'position',[200150450300],...'name','實(shí)例20');subplot(1,2,1)[x,y,z]=sphere(10);mesh(x,y,z)axisofftitle('Figure1:Opaque')hiddenonsubplot(1,2,2)[x,y,z]=sphere(10);mesh(x,y,z)axisofftitle('Figure2:Transparent')hiddenoff實(shí)例21PEAKS函數(shù)曲線functionshili21h0=figure('toolbar','none',...'position',[200100450450],...'name','實(shí)例21');[x,y,z]=peaks(30);subplot(2,1,1)x=x(1,:);y=y(:,1);i=find(y>0.8&y<1.2);j=find(x>-0.6&x<0.5);z(i,j)=nan*z(i,j);surfc(x,y,z)xlabel('X軸');ylabel('Y軸');zlabel('Z軸');title('Figure1:surfc函數(shù)形成的曲面')subplot(2,1,2)x=x(1,:);y=y(:,1);i=find(y>0.8&y<1.2);j=find(x>-0.6&x<0.5);z(i,j)=nan*z(i,j);surfl(x,y,z)xlabel('X軸');ylabel('Y軸');zlabel('Z軸');title('Figure2:surfl函數(shù)形成的曲面')實(shí)例22：片狀圖functionshili22h0=figure('toolbar','none',...'position',[200150550350],...'name','實(shí)例22');subplot(1,2,1)x=rand(1,20);y=rand(1,20);z=peaks(x,y*pi);t=delaunay(x,y);trimesh(t,x,y,z)hiddenofftitle('Figure1:TriangularSurfacePlot');subplot(1,2,2)x=rand(1,20);y=rand(1,20);z=peaks(x,y*pi);t=delaunay(x,y);trisurf(t,x,y,z)title('Figure1:TriangularSurfacePlot');實(shí)例23：視角的調(diào)整functionshili23h0=figure('toolbar','none',...'position',[200150450350],...'name','實(shí)例23');x=-5:0.5:5;[x,y]=meshgrid(x);r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;subplot(2,2,1)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure1')view(-37.5,30)subplot(2,2,2)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure2')view(-37.5+90,30)subplot(2,2,3)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure3')view(-37.5,60)subplot(2,2,4)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure4')view(180,0)實(shí)例24：向量場(chǎng)的繪制functionshili24h0=figure('toolbar','none',...'position',[200150450350],...'name','實(shí)例24');subplot(2,2,1)z=peaks;ribbon(z)title('Figure1')subplot(2,2,2)[x,y,z]=peaks(15);[dx,dy]=gradient(z,0.5,0.5);contour(x,y,z,10)holdonquiver(x,y,dx,dy)holdofftitle('Figure2')subplot(2,2,3)[x,y,z]=peaks(15);[nx,ny,nz]=surfnorm(x,y,z);surf(x,y,z)holdonquiver3(x,y,z,nx,ny,nz)holdofftitle('Figure3')subplot(2,2,4)x=rand(3,5);y=rand(3,5);z=rand(3,5);c=rand(3,5);fill3(x,y,z,c)gridontitle('Figure4')實(shí)例25：燈光定位functionshili25h0=figure('toolbar','none',...'position',[200150450250],...'name','實(shí)例25');vert=[111;121;221;211;112;122;222;212];fac=[1234;2673;4378;1584;1265;5678];gridoffsphere(36)h=findobj('type','surface');set(h,'facelighting','phong',...'facecolor',...'interp',...'edgecolor',[0.40.40.4],...'backfacelighting',...'lit')holdonpatch('faces',fac,'vertices',vert,...'facecolor','y');light('position',[132]);light('position',[-3-13]);materialshinyaxisvis3doffholdoff實(shí)例26：柱狀圖functionshili26h0=figure('toolbar','none',...'position',[20050450450],...'name','實(shí)例26');subplot(2,1,1)x=[521873986555432];bar(x)xlabel('X軸');ylabel('Y軸');title('第一子圖');subplot(2,1,2)y=[521873986555432];barh(y)xlabel('X軸');ylabel('Y軸');title('第二子圖');實(shí)例27：設(shè)置照明方式functionshili27h0=figure('toolbar','none',...'position',[200150450350],...'name','實(shí)例27');subplot(2,2,1)sphereshadingflatcamlightleftcamlightrightlightingflatcolorbaraxisofftitle('Figure1')subplot(2,2,2)sphereshadingflatcamlightleftcamlightrightlightinggouraudcolorbaraxisofftitle('Figure2')subplot(2,2,3)sphereshadinginterpcamlightrightcamlightleftlightingphongcolorbaraxisofftitle('Figure3')subplot(2,2,4)sphereshadingflatcamlightleftcamlightrightlightingnonecolorbaraxisofftitle('Figure4')實(shí)例28：羽狀圖functionshili28h0=figure('toolbar','none',...'position',[200150450350],...'name','實(shí)例28');subplot(2,1,1)alpha=90:-10:0;r=ones(size(alpha));m=alpha*pi/180;n=r*10;[u,v]=pol2cart(m,n);feather(u,v)title('羽狀圖')axis([020010])subplot(2,1,2)t=0:0.5:10;x=0.05+i;y=exp(-x*t);feather(y)title('復(fù)數(shù)矩陣的羽狀圖')實(shí)例29：立體透視（1）functionshili29h0=figure('toolbar','none',...'position',[200150450250],...'name','實(shí)例29');[x,y,z]=meshgrid(-2:0.1:2,...-2:0.1:2,...-2:0.1:2);v=x.*exp(-x.^2-y.^2-z.^2);gridonfori=-2:0.5:2;h1=surf(linspace(-2,2,20),...linspace(-2,2,20),...zeros(20)+i);rotate(h1,[1-11],30)dx=get(h1,'xdata');dy=get(h1,'ydata');dz=get(h1,'zdata');delete(h1)slice(x,y,z,v,[-22],2,-2)holdonslice(x,y,z,v,dx,dy,dz)holdoffaxistightview(-5,10)drawnowend實(shí)例30：立體透視（2）functionshili30h0=figure('toolbar','none',...'position',[200150450250],...'name','實(shí)例30');[x,y,z]=meshgrid(-2:0.1:2,...-2:0.1:2,...-2:0.1:2);v=x.*exp(-x.^2-y.^2-z.^2);[dx,dy,dz]=cylinder;slice(x,y,z,v,[-22],2,-2)fori=-2:0.2:2h=surface(dx+i,dy,dz);rotate(h,[100],90)xp=get(h,'xdata');yp=get(h,'ydata');zp=get(h,'zdata');delete(h)holdonhs=slice(x,y,z,v,xp,yp,zp);axistightxlim([-33])view(-10,35)drawnowdelete(hs)holdoffend實(shí)例31：表面圖形functionshili31h0=figure('toolbar','none',...'position',[200150550250],...'name','實(shí)例31');subplot(1,2,1)x=rand(100,1)*16-8;y=rand(100,1)*16-8;r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;xlin=linspace(min(x),max(x),33);ylin=linspace(min(y),max(y),33);[X,Y]=meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y,'cubic');mesh(X,Y,Z)axistightholdonplot3(x,y,z,'.','Markersize',20)subplot(1,2,2)k=5;n=2^k-1;theta=pi*(-n:2:n)/n;phi=(pi/2)*(-n:2:n)'/n;X=cos(phi)*cos(theta);Y=cos(phi)*sin(theta);Z=sin(phi)*ones(size(theta));colormap([000;111])C=hadamard(2^k);surf(X,Y,Z,C)axissquare實(shí)例32：沿曲線移動(dòng)的小球h0=figure('toolbar','none',...'position',[19856408468],...'name','實(shí)例32');h1=axes('parent',h0,...'position',[0.150.450.70.5],...'visible','on');t=0:pi/24:4*pi;y=sin(t);plot(t,y,'b')n=length(t);h=line('color',[00.50.5],...'linestyle','.',...'markersize',25,...'erasemode','xor'