Matlab100个实例程序-matlab简单程序实例.pdf

《Matlab100个实例程序-matlab简单程序实例.pdf》由会员分享，可在线阅读，更多相关《Matlab100个实例程序-matlab简单程序实例.pdf（167页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、Matlab100 个实例程序-matlab 简单程序实例程序代码：程序代码：(代码标记 code./codecode./code )1-32 是：图形应用篇33-66 是：界面设计篇67-84 是：图形处理篇85-100 是：数值分析篇实例 1：三角函数曲线（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);xlabel(自变量 X);ylabel(函数值

2、Y);title(SIN()函数曲线);grid on实例 2：三角函数曲线（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 onxlabel(自变量 X);ylabel(函数值 Y);title(三角函数);1Matlab100 个实例程序-matlab 简单程序实例实例 3：图形的叠加function shili03h0=figure(toolbar,none

3、,.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 onxlabel(自变量 X);ylabel(函数值 Y);title(三角函数);实例 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,hlin

4、e2=plotyy(x,y1,x,y2,semilogy,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;2Matlab100 个实例程序-matlab 简单程序实例x,y=meshgrid(t);subplot(2,2,1)plot(sin(t),cos(t)axis equals

5、ubplot(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 450 400,.name,实例 06);t=0:pi/10:2*pi;h=plot(t,sin(t);xlabel(t=0

6、 到 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(max(y)=y);text(x(imin),y(imin),.leftarrow 最小值=,num2str(y(imin),.fontsize,16)text(x(imax),y(imax),.3Matlab100 个实例程序-matlab 简单程序实例 leftarrow 最大值=,num2str(y(imax

7、),.fontsize,16)实例 7：条形图形function shili07h0=figure(toolbar,none,.position,200 150 450 350,.name,实例 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,yaxislocatio

8、n,right,color,none,xticklabel,)实例 8：区域图形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,profits1,facecolor,0.5 0.9 0.6,.edgecolor,b,.linewidth,3)hold onarea(x,p

9、rofits2,facecolor,0.9 0.85 0.7,.edgecolor,y,.linewidth,3)4Matlab100 个实例程序-matlab 简单程序实例hold onarea(x,profits3,facecolor,0.3 0.6 0.7,.edgecolor,r,.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 第一季度销量)gt

10、ext(leftarrow 第二季度销量)gtext(leftarrow 第三季度销量)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 68 45;68 54 69;x=sum(t);h=pie(x);textobjs=findobj(h,type,text);str1

11、=get(textobjs,string);val1=get(textobjs,extent);oldext=cat(1,val1:);names=商品一：;商品二：;商品三：;str2=strcat(names,str1);set(textobjs,string,str2)5Matlab100 个实例程序-matlab 简单程序实例val2=get(textobjs,extent);newext=cat(1,val2:);offset=sign(oldext(:,1).*(newext(:,3)-oldext(:,3)/2;pos=get(textobjs,position);textpos

12、=cat(1,pos:);textpos(:,1)=textpos(:,1)+offset;set(textobjs,position,num2cell(textpos,3,2)实例 10：阶梯图function shili10h0=figure(toolbar,none,.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=函数 e-(alpha*t)sinbeta*t 的阶梯图;gtext(gl

13、abel,fontsize,16)xlabel(t=0:10,fontsize,16)axis(0 10-1.2 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);6Matlab100 个实例程序-matlab 简单程序实例hold offh3=h1(1);h2;legend(h3,y1+y2,y1=

14、sin(x),y2=cos(x)xlabel(自变量 X);ylabel(函数值 Y);title(正弦函数与余弦函数的线性组合);实例 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;windpower=2 5 5 3 6 8 12 7 6 14 10 8;rdirection=winddirection*pi/180;x,y=pol2cart(rdirection

15、,windpower);compass(x,y);desc=风向和风力,北京气象台,10 月 1 日 0：00 到,10 月 1 日 12：00;gtext(desc)实例 13：轮廓图function shili13h0=figure(toolbar,none,.position,200 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)7Matlab100 个实例程序-matla

16、b 简单程序实例axis equalxlabel(实部,fontsize,16);ylabel(虚部,fontsize,16);h=polar(0 2*pi,0 1);delete(h)hold oncontour(x,y,abs(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

17、 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 off8Matlab100 个实例程序-matlab 简单程序实例实例 15：变换的傅立叶函数曲线function shili15h0=figure(toolbar,none,.position,200 150 450 250,.name,实例 15);axis equalm=movie

18、in(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,0 0 1 1);movie(m,30)实例 16：劳伦兹非线形方程的无序活动function shili15h0=figure(toolbar,none,.position,200 150 450 250,.name

19、,实例 15);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,0 0 1 1);movie(m,30)9Matlab100 个实例程序-matlab 简单程序实例实例 17：填充图function shili17h0=figure

20、(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,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=

21、exp(-x.*x);bar(x,y)title(2-D Bar Chart)subplot(2,2,2)x=-3:0.2:3;y=exp(-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)1 0Matlab100 个实例程序-matlab 简单程序实例subplot(2,2,4)x=-3:0.2:3;y=exp(-x.*x);barh(x,y)title(Horizontal Bar Chart)实例 19：三维曲线图function

22、 shili19h0=figure(toolbar,none,.position,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)subplot(2,1,2)x=linspace(0,2*pi

23、);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 轴);1 1Matlab100 个实例程序-matlab 简单程序实例ylabel(Y 轴);zlabel(Z 轴);title(Figure2:3-D Plot)实例 20：图形的隐藏属性function shili20h0=figure(toolbar,none,.position,200 150 450 300,.name,实例 20);subplot(1

24、,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实例 21：PEAKS 函数曲线function shili21h0=figure(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

25、=find(y0.8&y-0.6&x0.8&y-0.6&xn,.m=m+1;,.i=1;,.end,.end,.end);k2=uicontrol(parent,h0,.style,pushbutton,.position,180 100 50 30,.string,停止,.callback,.k=0;,.set(e1,string,m),.p=get(h,xdata);,.q=get(h,ydata);,.2 2Matlab100 个实例程序-matlab 简单程序实例 set(e2,string,p);,.set(e3,string,q);k3=uicontrol(parent,h0,.s

26、tyle,pushbutton,.position,280 100 50 30,.string,关闭,.callback,close);e1=uicontrol(parent,h0,.style,edit,.position,60 30 60 20);t1=uicontrol(parent,h0,.style,text,.string,循环次数,.position,60 50 60 20);e2=uicontrol(parent,h0,.style,edit,.position,180 30 50 20);t2=uicontrol(parent,h0,.style,text,.string,终

27、点的 X 坐标值,.position,155 50 100 20);e3=uicontrol(parent,h0,.style,edit,.position,300 30 50 20);t3=uicontrol(parent,h0,.style,text,.string,终点的 Y 坐标值,.position,275 50 100 20);实例 33：曲线转换按钮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 on2 3Matlab100 个

28、实例程序-matlab 简单程序实例huidiao=.if i=1,.i=0;,.y=cos(x);,.delete(h),.set(hm,string,正弦函数),.h=plot(x,y);,.grid on,.else if i=0,.i=1;,.y=sin(x);,.set(hm,string,余弦函数),.delete(h),.h=plot(x,y);,.grid on,.end,.end;hm=uicontrol(gcf,style,pushbutton,.string,余弦函数,.callback,huidiao);i=1;set(hm,position,250 20 60 20)

29、;set(gca,position,0.2 0.2 0.6 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,.;2 4Matlab100 个实例程序-matlab 简单程序实例huidiao2=.set(h_toggle1,value,0),.grid off,.;h_toggle1=uicontrol(g

30、cf,style,togglebutton,.string,grid on,.value,0,.position,20 45 50 20,.callback,huidiao1);h_toggle2=uicontrol(gcf,style,togglebutton,.string,grid off,.value,0,.position,20 20 50 20,.callback,huidiao2);set(gca,position,0.2 0.2 0.6 0.6)title(开关按钮的使用)实例 35：编辑框的使用h0=figure(toolbar,none,.position,200 150

31、350 250,.name,实例 35);f=Please input the letter;huidiao1=.g=upper(f);,.set(h2_edit,string,g),.;huidiao2=.g=lower(f);,.set(h2_edit,string,g),.;h1_edit=uicontrol(gcf,style,edit,.position,100 200 100 50,.HorizontalAlignment,left,.string,Please input the letter,.callback,f=get(h1_edit,string);,.backgroun

32、d,w,.2 5Matlab100 个实例程序-matlab 简单程序实例 max,5,.min,1);h2_edit=uicontrol(gcf,style,edit,.HorizontalAlignment,left,.position,100 100 100 50,.background,w,.max,5,.min,1);h1_button=uicontrol(gcf,style,pushbutton,.string,小写变大写,.position,100 45 100 20,.callback,huidiao1);h2_button=uicontrol(gcf,style,pushbu

33、tton,.string,大写变小写,.position,100 20 100 20,.callback,huidiao2);实例 36：弹出式菜单h0=figure(toolbar,none,.position,200 150 450 250,.name,实例 36);x=0:0.5:2*pi;y=sin(x);h=plot(x,y);grid onhm=uicontrol(gcf,style,popupmenu,.string,.sin(x)|cos(x)|sin(x)+cos(x)|exp(-sin(x),.position,250 20 50 20);set(hm,value,1)hu

34、idiao=.v=get(hm,value);,.switch v,.case 1,.delete(h),.y=sin(x);,.h=plot(x,y);,.2 6Matlab100 个实例程序-matlab 简单程序实例 grid on,.case 2,.delete(h),.y=cos(x);,.h=plot(x,y);,.grid on,.case 3,.delete(h),.y=sin(x)+cos(x);,.h=plot(x,y);,.grid on,.case 4,.delete(h),.y=exp(-sin(x);,.h=plot(x,y);,.grid on,.end;set(

35、hm,callback,huidiao)set(gca,position,0.2 0.2 0.6 0.6)title(弹出式菜单的使用)实例 37：滑标的使用h0=figure(toolbar,none,.position,200 150 450 250,.name,实例 37);x,y=meshgrid(-8:0.5:8);r=sqrt(x.2+y.2)+eps;z=sin(r)./r;h0=mesh(x,y,z);h1=axes(position,.0.2 0.2 0.5 0.5,.visible,off);htext=uicontrol(gcf,.units,points,.positi

36、on,20 30 45 15,.string,brightness,.style,text);2 7Matlab100 个实例程序-matlab 简单程序实例hslider=uicontrol(gcf,.units,points,.position,10 10 300 15,.min,-1,.max,1,.style,slider,.callback,.brighten(get(hslider,value);实例 38：多选菜单h0=figure(toolbar,none,.position,200 150 450 250,.name,实例 38);x,y=meshgrid(-8:0.5:8)

37、;r=sqrt(x.2+y.2)+eps;z=sin(r)./r;h0=mesh(x,y,z);hlist=uicontrol(gcf,style,listbox,.string,default|spring|summer|autumn|winter,.max,5,.min,1,.position,20 20 80 100,.callback,.k=get(hlist,value);,.switch k,.case 1,.colormap default,.case 2,.colormap spring,.case 3,.colormap summer,.case 4,.colormap au

38、tumn,.case 5,.colormap winter,.end);2 8Matlab100 个实例程序-matlab 简单程序实例实例 39：菜单控制的使用h0=figure(toolbar,none,.position,200 150 450 250,.name,实例 39);x=0:0.5:2*pi;y=cos(x);h=plot(x,y);grid onset(gcf,toolbar,none)hm=uimenu(label,example);huidiao1=.set(hm_gridon,checked,on),.set(hm_gridoff,checked,off),.grid

39、 on;huidiao2=.set(hm_gridoff,checked,on),.set(hm_gridon,checked,off),.grid off;hm_gridon=uimenu(hm,label,grid on,.checked,on,.callback,huidiao1);hm_gridoff=uimenu(hm,label,grid off,.checked,off,.callback,huidiao2);实例 40：UIMENU 菜单的应用h0=figure(toolbar,none,.position,200 150 450 250,.name,实例 40);h1=uim

40、enu(gcf,label,函数);h11=uimenu(h1,label,轮廓图,.callback,.set(h31,checked,on),.set(h32,checked,off),.x,y,z=peaks;,.contour3(x,y,z,30);h12=uimenu(h1,label,高斯分布,.2 9Matlab100 个实例程序-matlab 简单程序实例 callback,.set(h31,checked,on),.set(h32,checked,off),.mesh(peaks);,.axis tight);h13=uimenu(h1,label,Sinc 函数,.call

41、back,.set(h31,checked,on),.set(h32,checked,off),.x,y=meshgrid(-8:0.5:8);,.r=sqrt(x.2+y.2)+eps;,.z=sin(r)./r;,.mesh(x,y,z);h2=uimenu(gcf,label,色彩);hl2(1)=uimenu(h2,label,Default,.checked,on,.callback,.set(hl2,checked,off),.set(hl2(1),checked,on),.colormap(default);hl2(2)=uimenu(h2,label,spring,.callb

42、ack,.set(hl2,checked,off),.set(hl2(2),checked,on),.colormap(spring);hl2(3)=uimenu(h2,label,Summer,.callback,.set(hl2,checked,off),.set(hl2(3),checked,on),.colormap(summer);hl2(4)=uimenu(h2,label,Autumn,.callback,.set(hl2,checked,off),.3 0Matlab100 个实例程序-matlab 简单程序实例 set(hl2(4),checked,on),.colormap

43、(autumn);hl2(5)=uimenu(h2,label,Winter,.callback,.set(hl2,checked,off),.set(hl2(5),checked,on),.colormap(winter);h3=uimenu(gcf,label,坐标选项);h31=uimenu(h3,label,Axis on,.callback,.axis on,.set(h31,checked,on),.set(h32,checked,off);h32=uimenu(h3,label,Axis off,.callback,.axis off,.set(h32,checked,on),.

44、set(h31,checked,off);实例 41：除法计算器h=figure(toolbar,none,.position,200 150 450 250,.name,实例 41);h1=uicontrol(gcf,style,edit,.position,80 200 100 20,.HorizontalAlignment,right,.callback,m=get(h1,string);,.a=str2num(m););h2=uicontrol(gcf,style,edit,.HorizontalAlignment,right,.position,80 150 100 20,.call

45、back,n=get(h2,string);,.b=str2num(n););h3=uicontrol(gcf,style,text,.3 1Matlab100 个实例程序-matlab 简单程序实例 string,被除数,.position,80 230 100 20);h4=uicontrol(gcf,style,edit,.position,80 50 100 20);h5=uicontrol(gcf,style,pushbutton,.position,80 100 100 20,.string,=,.callback,.if b=0,.h7=errordlg(除数不能为 0！,err

46、or,on);,.else,.k=a/b;,.c=num2str(k);,.set(h4,string,c),.end);h8=uicontrol(gcf,style,text,.string,除数,.position,80 175 100 20);h9=uicontrol(gcf,style,text,.string,商,.position,80 75 100 20);实例 42：单选框的使用h0=figure(toolbar,none,.position,200 150 450 250,.name,实例 42);x=0:0.5:2*pi;y=sin(x);plot(x,y)grid ons

47、et(gcf,toolbar,none)g=set(gca,position,0.2 0.2 0.6 0.6);huidiao1=.grid on,.set(box_on,value,1),.set(box_off,value,0),;huidiao2=.3 2Matlab100 个实例程序-matlab 简单程序实例 grid off,.set(box_off,value,1),.set(box_on,value,0),;box_on=uicontrol(gcf,style,radio,.position,5 50 50 20,.string,grid on,.value,1,.callba

48、ck,huidiao1);box_off=uicontrol(gcf,style,radio,.position,5 20 50 20,.string,grid off,.value,0,.callback,huidiao2);title(无线按钮的使用)实例 43：添加环境效果h0=figure(toolbar,none,.position,198 56 450 468,.name,实例 43);h1=axes(parent,h0,.position,0.15 0.45 0.7 0.5,.visible,off);sphereh=findobj(type,surface);shading i

49、nterpaxis equall=light(position,0-2 1);k(1)=get(h,specularstrength);k(2)=get(h,diffusestrength);k(3)=get(h,specularexponent);k(4)=get(h,specularcolorreflectance);u1=uimenu(parent,h0,.label,灯光效果,.tag,u1,.backgroundcolor,0.75 0.75 0.75);u11=uimenu(parent,u1,.label,gouraud,.3 3Matlab100 个实例程序-matlab 简单

50、程序实例 tag,u11,.backgroundcolor,0.75 0.75 0.75,.callback,.set(u11,checked,on);,.set(u12,checked,off);,.set(h,facelighting,gouraud);u12=uimenu(parent,u1,.label,phong,.tag,u12,.backgroundcolor,0.75 0.75 0.75,.callback,.set(u11,checked,off);,.set(u12,checked,on);,.set(h,facelighting,gouraud);u2=uimenu(pa