matlab實用程序百例1

實例 1：三角函數(shù)曲線（1） function shili01 h0=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('函數(shù)值 Y'); title('SIN( )函數(shù)曲線'); grid on
實例 2：三角函數(shù)曲線（2） function shili02 h0=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 on xlabel('自變量 X'); ylabel('函數(shù)值 Y'); title('三角函數(shù)');
實例 3：圖形的疊加 function shili03 h0=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 xlabel('自變量 X'); ylabel('函數(shù)值 Y'); title('三角函數(shù)');
實例 4：雙 y 軸圖形的繪制 function shili04 h0=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,'semilogy','plot'); axes(haxes(1)) ylabel('semilog plot'); axes(haxes(2)) ylabel('linear plot');
實例 5：單個軸窗口顯示多個圖形 function shili05 h0=figure('toolbar','none',... 'position',[200 150 450 250],... 'name','實例 05'); t=0:pi/10:2*pi; [x,y]=meshgrid(t); subplot(2,2,1) plot(sin(t),cos(t)) axis equal subplot(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：圖形標(biāo)注 function shili06 h0=figure('toolbar','none',... 'position',[200 150 450 400],... '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：條形圖形 function shili07 h0=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,'yaxislocation','right','color','none','xticklabel',[])
實例 8：區(qū)域圖形 function shili08 h0=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 on area(x,profits2,'facecolor',[0.9 0.85 0.7],... 'edgecolor','y',... 'linewidth',3) hold on area(x,profits3,'facecolor',[0.3 0.6 0.7],... 'edgecolor','r',... 'linewidth',3) hold on area(x,profits4,'facecolor',[0.6 0.5 0.9],... 'edgecolor','m',... 'linewidth',3) hold off set(gca,'xtick',[91:95]) set(gca,'layer','top') gtext('\leftarrow 第一季度銷量') gtext('\leftarrow 第二季度銷量') gtext('\leftarrow 第三季度銷量') gtext('\leftarrow 第四季度銷量') xlabel('年','fontsize',16); ylabel('銷售量','fontsize',16);
實例 9：餅圖的繪制 function shili09 h0=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=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：階梯圖 function shili10 h0=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 on plot(t,f,':*') hold off glabel='函數(shù) e^{-(\alpha*t)}sin\beta*t 的階梯圖'; gtext(glabel,'fontsize',16) xlabel('t=0:10','fontsize',16) axis([0 10 -1.2 1.2]) 實例 11：枝干圖 function shili11 h0=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 on h2=plot(x,y1,'^r',x,y2,'*g'); hold off h3=[h1(1);h2]; legend(h3,'y1+y2','y1=sin(x)','y2=cos(x)') xlabel('自變量 X'); ylabel('函數(shù)值 Y'); title('正弦函數(shù)與余弦函數(shù)的線性組合');