matlab實用教程100例

1、matlab實用程序百例1-32是：圖形應(yīng)用篇33-66是：界面設(shè)計篇67-84是：圖形處理篇85-100是：數(shù)值分析篇實例1：三角函數(shù)曲線（1）function shili01h0=figure('toolbar','none',. 'position',198 56 350 300,. 'name','實例01');h1=axes('parent',h0,. 'visible','off');x=-pi:0.05:pi;y=sin(x);plot(x,y);xla

2、bel('自變量X');ylabel('函數(shù)值Y');title('SIN( )函數(shù)曲線');grid on實例2：三角函數(shù)曲線（2）function shili02h0=figure('toolbar','none',. 'position',200 150 450 350,. 'name','實例02');x=-pi:0.05:pi;y=sin(x)+cos(x);plot(x,y,'-*r','linewidth',1);grid

3、 onxlabel('自變量X');ylabel('函數(shù)值Y');title('三角函數(shù)');實例3：圖形的疊加function shili03h0=figure('toolbar','none',. 'position',200 150 450 350,. 'name','實例03');x=-pi:0.05:pi;y1=sin(x);y2=cos(x);plot(x,y1,. '-*r',. x,y2,. '-og');grid on

4、xlabel('自變量X');ylabel('函數(shù)值Y');title('三角函數(shù)');實例4：雙y軸圖形的繪制function shili04h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name','實例04');x=0:900;a=1000;b=0.005;y1=2*x;y2=cos(b*x);haxes,hline1,hline2=plotyy(x,y1,x,y2,'semilog

5、y','plot');axes(haxes(1)ylabel('semilog plot');axes(haxes(2)ylabel('linear plot');實例5：單個軸窗口顯示多個圖形function shili05h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name','實例05');t=0:pi/10:2*pi;x,y=meshgrid(t);subplot(2,2,1)p

6、lot(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)h=sin(x)+cos(y);plot(t,h)axis(0 2*pi -2 2)subplot(2,2,4)g=(sin(x).2)-(cos(y).2);plot(t,g)axis(0 2*pi -1 1)實例6：圖形標注function shili06h0=figure('toolbar','none',. 'position',200 150 4

7、50 400,. 'name','實例06');t=0:pi/10:2*pi;h=plot(t,sin(t);xlabel('t=0到2pi','fontsize',16);ylabel('sin(t)','fontsize',16);title('it從 0to2pi 的正弦曲線','fontsize',16)x=get(h,'xdata');y=get(h,'ydata');imin=find(min(y)=y);imax=find

8、(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) 實例7：條形圖形function shili07h0=figure('toolbar','none',. 'position',200 150 450 350,. 'nam

9、e','實例07');tiao1=562 548 224 545 41 445 745 512;tiao2=47 48 57 58 54 52 65 48;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',

10、9;color','none','xticklabel',)實例8：區(qū)域圖形function shili08h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name','實例08');x=91:95;profits1=88 75 84 93 77;profits2=51 64 54 56 68;profits3=42 54 34 25 24;profits4=26 38 18 15 4;area(x,pro

11、fits1,'facecolor',0.5 0.9 0.6,. 'edgecolor','b',. 'linewidth',3)hold onarea(x,profits2,'facecolor',0.9 0.85 0.7,. 'edgecolor','y',. 'linewidth',3)hold onarea(x,profits3,'facecolor',0.3 0.6 0.7,. 'edgecolor','r',.

12、 'linewidth',3)hold onarea(x,profits4,'facecolor',0.6 0.5 0.9,. 'edgecolor','m',. 'linewidth',3)hold offset(gca,'xtick',91:95)set(gca,'layer','top')gtext('leftarrow第一季度銷量')gtext('leftarrow第二季度銷量')gtext('leftarrow第三季度

13、銷量')gtext('leftarrow第四季度銷量')xlabel('年','fontsize',16);ylabel('銷售量','fontsize',16); 實例9：餅圖的繪制function shili09h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name','實例09');t=54 21 35; 68 54 35; 45 25 12; 48

14、68 45; 68 54 69;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

15、,'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)實例10：階梯圖function shili10h0=figure('toolbar','none'

16、,. 'position',200 150 450 400,. 'name','實例10');a=0.01;b=0.5;t=0:10;f=exp(-a*t).*sin(b*t);stairs(t,f)hold onplot(t,f,':*')hold offglabel='函數(shù)e-(alpha*t)sinbeta*t的階梯圖'gtext(glabel,'fontsize',16)xlabel('t=0:10','fontsize',16)axis(0 10 -1.2

17、1.2)實例11：枝干圖function shili11h0=figure('toolbar','none',. 'position',200 150 450 350,. 'name','實例11');x=0:pi/20:2*pi;y1=sin(x);y2=cos(x);h1=stem(x,y1+y2);hold onh2=plot(x,y1,'r',x,y2,'*g');hold offh3=h1(1);h2;legend(h3,'y1+y2','y1=si

18、n(x)','y2=cos(x)')xlabel('自變量X');ylabel('函數(shù)值Y');title('正弦函數(shù)與余弦函數(shù)的線性組合');實例12：羅盤圖function shili12h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name','實例12');winddirection=54 24 65 84 256 12 235 62 125 324 34 254;

19、windpower=2 5 5 3 6 8 12 7 6 14 10 8;rdirection=winddirection*pi/180;x,y=pol2cart(rdirection,windpower);compass(x,y);desc='風(fēng)向和風(fēng)力', '北京氣象臺', '10月1日0：00到', '10月1日12：00'gtext(desc)實例13：輪廓圖function shili13h0=figure('toolbar','none',. 'position',200

20、 150 450 250,. 'name','實例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)axis equalxlabel('實部','fontsize',16);ylabel('虛部','fontsize',16);h=polar(0 2*pi,0 1);delete(h)hold oncontour(x,y,abs

21、(f),20)實例14：交互式圖形function shili14h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name','實例14');axis(0 10 0 10);hold onx=;y=;n=0;disp('單擊鼠標左鍵點取需要的點');disp('單擊鼠標右鍵點取最后一個點');but=1;while but=1 xi,yi,but=ginput(1); plot(xi,yi,'bo'

22、) n=n+1; disp('單擊鼠標左鍵點取下一個點'); 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-');hold off實例14：交互式圖形function shili14h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name','實例14');axis(0 10 0

23、10);hold onx=;y=;n=0;disp('單擊鼠標左鍵點取需要的點');disp('單擊鼠標右鍵點取最后一個點');but=1;while but=1 xi,yi,but=ginput(1); plot(xi,yi,'bo') n=n+1; disp('單擊鼠標左鍵點取下一個點'); 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-');hold off實例15：變換的傅

24、立葉函數(shù)曲線function shili15h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name','實例15');axis equalm=moviein(20,gcf);set(gca,'nextplot','replacechildren')h=uicontrol('style','slider','position',. 100 10 500 20,'

25、;min',1,'max',20)for j=1:20 plot(fft(eye(j+16) set(h,'value',j) m(:,j)=getframe(gcf);endclf;axes('position',0 0 1 1);movie(m,30)實例16：勞倫茲非線形方程的無序活動function shili15h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name','實例15'

26、;);axis equalm=moviein(20,gcf);set(gca,'nextplot','replacechildren')h=uicontrol('style','slider','position',. 100 10 500 20,'min',1,'max',20)for j=1:20 plot(fft(eye(j+16) set(h,'value',j) m(:,j)=getframe(gcf);endclf;axes('position&#

27、39;,0 0 1 1);movie(m,30)實例17：填充圖function shili17h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name','實例17');t=(1:2:15)*pi/8;x=sin(t);y=cos(t);fill(x,y,'r')axis square offtext(0,0,'STOP',. 'color',1 1 1,. 'fontsize',

28、50,. 'horizontalalignment','center')例18：條形圖和階梯形圖function shili18h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name','實例18');subplot(2,2,1)x=-3:0.2:3;y=exp(-x.*x);bar(x,y)title('2-D Bar Chart')subplot(2,2,2)x=-3:0.2:3;y=exp(

29、-x.*x);bar3(x,y,'r')title('3-D Bar Chart')subplot(2,2,3)x=-3:0.2:3;y=exp(-x.*x);stairs(x,y)title('Stair Chart')subplot(2,2,4)x=-3:0.2:3;y=exp(-x.*x);barh(x,y)title('Horizontal Bar Chart') 實例19：三維曲線圖function shili19h0=figure('toolbar','none',. 'posi

30、tion',200 150 450 400,. 'name','實例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)grid onxlabel('X軸');ylabel('Y軸');zlabel('Z軸');title('Figure1:3-D Plot')subpl

31、ot(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)grid onxlabel('X軸');ylabel('Y軸');zlabel('Z軸');title('Figure2:3-D Plot')實例20：圖形的隱藏屬性function shili20h0=figure('toolbar','none'

32、;,. 'position',200 150 450 300,. 'name','實例20');subplot(1,2,1)x,y,z=sphere(10);mesh(x,y,z)axis offtitle('Figure1:Opaque')hidden onsubplot(1,2,2)x,y,z=sphere(10);mesh(x,y,z)axis offtitle('Figure2:Transparent')hidden off實例21PEAKS函數(shù)曲線function shili21h0=figure(

33、9;toolbar','none',. 'position',200 100 450 450,. 'name','實例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軸'

34、);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ù)形成的曲面') 實例22：片狀圖function shili22h0=figure(

35、9;toolbar','none',. 'position',200 150 550 350,. 'name','實例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)hidden offtitle('Figure1:Triangular Surface Plot');subplot(1,2,2)x=rand(1,20);y=rand(1,20);z=peaks(x,y*p

36、i);t=delaunay(x,y);trisurf(t,x,y,z)title('Figure1:Triangular Surface Plot');實例23：視角的調(diào)整function shili23h0=figure('toolbar','none',. 'position',200 150 450 350,. 'name','實例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

37、(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&#

38、39;)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)實例24：向量場的繪制function shili24h0=figure('toolbar','none',. 'p

39、osition',200 150 450 350,. 'name','實例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)hold onquiver(x,y,dx,dy)hold offtitle('Figure2')subplot(2,2,3)x,y,z=peaks(15);nx,ny,nz=surfnorm(x,y,z)

40、;surf(x,y,z)hold onquiver3(x,y,z,nx,ny,nz)hold offtitle('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)grid ontitle('Figure4')實例25：燈光定位function shili25h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name

41、9;,'實例25');vert=1 1 1;1 2 1; 2 2 1;2 1 1; 1 1 2;1 2 2; 2 2 2;2 1 2;fac=1 2 3 4;2 6 7 3; 4 3 7 8;1 5 8 4; 1 2 6 5;5 6 7 8;grid offsphere(36)h=findobj('type','surface');set(h,'facelighting','phong',. 'facecolor',. 'interp',. 'edgecolor',0

42、.4 0.4 0.4,. 'backfacelighting',. 'lit')hold onpatch('faces',fac,'vertices',vert,. 'facecolor','y');light('position',1 3 2);light('position',-3 -1 3);material shinyaxis vis3d offhold off實例26：柱狀圖function shili26h0=figure('toolbar'

43、;,'none',. 'position',200 50 450 450,. 'name','實例26');subplot(2,1,1)x=5 2 1 8 7 3 9 8 6 5 5 5 4 3 2;bar(x)xlabel('X軸');ylabel('Y軸');title('第一子圖');subplot(2,1,2)y=5 2 1 8 7 3 9 8 6 5 5 5 4 3 2;barh(y)xlabel('X軸');ylabel('Y軸');tit

44、le('第二子圖');實例27：設(shè)置照明方式function shili27h0=figure('toolbar','none',. 'position',200 150 450 350,. 'name','實例27');subplot(2,2,1)sphereshading flatcamlight leftcamlight rightlighting flatcolorbaraxis offtitle('Figure1')subplot(2,2,2)sphereshading f

45、latcamlight leftcamlight rightlighting gouraudcolorbaraxis offtitle('Figure2')subplot(2,2,3)sphereshading interpcamlight rightcamlight leftlighting phongcolorbaraxis offtitle('Figure3')subplot(2,2,4)sphereshading flatcamlight leftcamlight rightlighting nonecolorbaraxis offtitle('

46、Figure4')實例28：羽狀圖function shili28h0=figure('toolbar','none',. 'position',200 150 450 350,. 'name','實例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(0 20 0 10)subplot(2,1,2)t

47、=0:0.5:10;x=0.05+i;y=exp(-x*t);feather(y)title('復(fù)數(shù)矩陣的羽狀圖')實例29：立體透視（1）function shili29h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name','實例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);grid onfor i=-2:0.5:2;

48、h1=surf(linspace(-2,2,20),. linspace(-2,2,20),. zeros(20)+i); rotate(h1,1 -1 1,30) dx=get(h1,'xdata'); dy=get(h1,'ydata'); dz=get(h1,'zdata'); delete(h1) slice(x,y,z,v,-2 2,2,-2) hold on slice(x,y,z,v,dx,dy,dz) hold off axis tight view(-5,10) drawnowend實例30：立體透視（2）function sh

49、ili30h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name','實例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,-2 2,2,-2)for i=-2:0.2:2 h=surface(dx+i,dy,dz); rotate(h,1 0 0,90) xp=get(h,'

50、;xdata'); yp=get(h,'ydata'); zp=get(h,'zdata'); delete(h) hold on hs=slice(x,y,z,v,xp,yp,zp); axis tight xlim(-3 3) view(-10,35) drawnow delete(hs) hold offend實例31：表面圖形function shili31h0=figure('toolbar','none',. 'position',200 150 550 250,. 'name'

51、,'實例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)axis tighthold onplot3(x,y,z,'.','Markersize',

52、20)subplot(1,2,2)k=5;n=2k-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(0 0 0;1 1 1)C=hadamard(2k);surf(X,Y,Z,C)axis square實例32：沿曲線移動的小球h0=figure('toolbar','none',. 'position',198 56 408 468,.

53、 'name','實例32');h1=axes('parent',h0,. 'position',0.15 0.45 0.7 0.5,. 'visible','on');t=0:pi/24:4*pi;y=sin(t);plot(t,y,'b')n=length(t);h=line('color',0 0.5 0.5,. 'linestyle','.',. 'markersize',25,. 'erasemode&

54、#39;,'xor');k1=uicontrol('parent',h0,. 'style','pushbutton',. 'position',80 100 50 30,. 'string','開始',. 'callback',. 'i=1;',. 'k=1;,',. 'm=0;,',. 'while 1,',. 'if k=0,',. 'break,',. 'e

55、nd,',. 'if k=0,',. 'set(h,''xdata'',t(i),''ydata'',y(i),',. 'drawnow;,',. 'i=i+1;,',. 'if i>n,',. 'm=m+1;,',. 'i=1;,',. 'end,',. 'end,',. 'end');k2=uicontrol('parent',h0,. &

56、#39;style','pushbutton',. 'position',180 100 50 30,. 'string','停止',. 'callback',. 'k=0;,',. 'set(e1,''string'',m),',. 'p=get(h,''xdata'');,',. 'q=get(h,''ydata'');,',. 'set

57、(e2,''string'',p);,',. 'set(e3,''string'',q)');k3=uicontrol('parent',h0,. 'style','pushbutton',. 'position',280 100 50 30,. 'string','關(guān)閉',. 'callback','close');e1=uicontrol('parent',h0

58、,. 'style','edit',. 'position',60 30 60 20);t1=uicontrol('parent',h0,. 'style','text',. 'string','循環(huán)次數(shù)',. 'position',60 50 60 20);e2=uicontrol('parent',h0,. 'style','edit',. 'position',180 30 50 2

59、0);t2=uicontrol('parent',h0,. 'style','text',. 'string','終點的X坐標值',. 'position',155 50 100 20);e3=uicontrol('parent',h0,. 'style','edit',. 'position',300 30 50 20);t3=uicontrol('parent',h0,. 'style','t

60、ext',. 'string','終點的Y坐標值',. 'position',275 50 100 20);實例33：曲線轉(zhuǎn)換按鈕h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name','實例33');x=0:0.5:2*pi;y=sin(x);h=plot(x,y);grid onhuidiao=. 'if i=1,',. 'i=0;,',. '

61、y=cos(x);,',. 'delete(h),',. 'set(hm,''string'',''正弦函數(shù)''),',. 'h=plot(x,y);,',. 'grid on,',. 'else if i=0,',. 'i=1;,',. 'y=sin(x);,',. 'set(hm,''string'',''余弦函數(shù)''),',.

62、'delete(h),',. 'h=plot(x,y);,',. 'grid on,',. 'end,',. 'end'hm=uicontrol(gcf,'style','pushbutton',. 'string','余弦函數(shù)',. 'callback',huidiao);i=1;set(hm,'position',250 20 60 20);set(gca,'position',0.2 0.2 0.6

63、 0.6)title('按鈕的使用')hold on實例34：柵格控制按鈕h0=figure('toolbar','none',. 'position',200 150 450 250,. 'name','實例34');x=0:0.5:2*pi;y=sin(x);plot(x,y)huidiao1=. 'set(h_toggle2,''value'',0),',. 'grid on,',. ;huidiao2=. 'set(h_toggle1,''value'',0),',. 'grid off,',. ;h_toggle1=uicont