Matlab100个实例程序

程序代码：(代码标记[code]...[/code])1-32是：图形应用篇33-66是：界面设计篇67-84是：图形处理篇85-100是：数值分析篇实例1：三角函数曲线（1）functionshili01h0=figure('toolbar','none',...'position',[19856350300],...'name','实例01');h1=axes('parent',h0,...'visible','off');x=-pi:0.05:pi;y=sin(x);plot(x,y);xlabel('自变量X');ylabel('函数值Y');title('SIN()函数曲线');gridon实例2：三角函数曲线（2）functionshili02h0=figure('toolbar','none',...'position',[200150450350],...'name','实例02');x=-pi:0.05:pi;y=sin(x)+cos(x);plot(x,y,'-*r','linewidth',1);gridonxlabel('自变量X');ylabel('函数值Y');title('三角函数');实例3：图形的叠加functionshili03h0=figure('toolbar','none',...'position',[200150450350],...'name','实例03');x=-pi:0.05:pi;y1=sin(x);y2=cos(x);plot(x,y1,...'-*r',...x,y2,...'--og');gridonxlabel('自变量X');ylabel('函数值Y');title('三角函数');实例4：双y轴图形的绘制functionshili04h0=figure('toolbar','none',...'position',[200150450250],...'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,'semilogy','plot');axes(haxes(1))ylabel('semilogplot');axes(haxes(2))ylabel('linearplot');实例5：单个轴窗口显示多个图形functionshili05h0=figure('toolbar','none',...'position',[200150450250],...'name','实例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])实例6：图形标注functionshili06h0=figure('toolbar','none',...'position',[200150450400],...'name','实例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)实例7：条形图形functionshili07h0=figure('toolbar','none',...'position',[200150450350],...'name','实例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',[])实例8：区域图形functionshili08h0=figure('toolbar','none',...'position',[200150450250],...'name','实例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);实例9：饼图的绘制functionshili09h0=figure('toolbar','none',...'position',[200150450250],...'name','实例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]))实例10：阶梯图functionshili10h0=figure('toolbar','none',...'position',[200150450400],...'name','实例10');a=0.01;b=0.5;t=0:10;f=exp(-a*t).*sin(b*t);stairs(t,f)holdonplot(t,f,':*')holdoffglabel='函数e^{-(\alpha*t)}sin\beta*t的阶梯图';gtext(glabel,'fontsize',16)xlabel('t=0:10','fontsize',16)axis([010-1.21.2])实例11：枝干图functionshili11h0=figure('toolbar','none',...'position',[200150450350],...'name','实例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('函数值Y');title('正弦函数与余弦函数的线性组合');实例12：罗盘图functionshili12h0=figure('toolbar','none',...'position',[200150450250],...'name','实例12');winddirection=[54246584256122356212532434254];windpower=[255368127614108];rdirection=winddirection*pi/180;[x,y]=pol2cart(rdirection,windpower);compass(x,y);desc={'风向和风力','北京气象台','10月1日0：00到','10月1日12：00'};gtext(desc)实例13：轮廓图functionshili13h0=figure('toolbar','none',...'position',[200150450250],...'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)axisequalxlabel('实部','fontsize',16);ylabel('虚部','fontsize',16);h=polar([02*pi],[01]);delete(h)holdoncontour(x,y,abs(f),20)实例14：交互式图形functionshili14h0=figure('toolbar','none',...'position',[200150450250],...'name','实例14');axis([010010]);holdonx=[];y=[];n=0;disp('单击鼠标左键点取需要的点');disp('单击鼠标右键点取最后一个点');but=1;whilebut==1[xi,yi,but]=ginput(1);plot(xi,yi,'bo')n=n+1;disp('单击鼠标左键点取下一个点');x(n,1)=xi;y(n,1)=yi;e