数学建模案例matlab实用程序百例

实例 1：三角函数曲线(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( 函数值 Y );title( SIN( )函数曲 线 ); grid on实例 2：三角函数曲线(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 onxlabel( 自变量 X );ylabel( 函数值 Y );title( 三角函数 );实例 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 xlabel( 自变量 X );ylabel( 函数值 Y );title( 三角函数 );实例 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：图形标注 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 到 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 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, yaxislocation , 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,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 ,. 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 第三季度销量 ) gtext( \leftarrow 第四季度销量 ) xlabel( 年 , fontsiz