数学建模案例MATLAB实用程序百例.docx

1、实例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)

2、;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(三角函数)

3、;实例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

4、(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)ax

5、is(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到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);

6、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,ti

7、ao1) 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=8875849377;profits2=5164545668;profits3=425434

8、2524;profits4=263818154;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,. li

9、newidth,3) hold offset(gca,xtick,91:95)set(gca,layer,top) gtext(leftarrow第一季度销量) gtext(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;68543

10、5;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

11、(:,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,.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

12、(t,f,:*) hold offglabel=函数e-(alpha*t)sinbeta*t的阶梯图; gtext(glabel,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); ho

13、ld offh3=h1(1);h2;legend(h3,y1+y2,y1=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 84256 12 235 62125 324 34 254;windpower=2 5 5 36 8 12 76 14 10 8; rdirection=winddirection*pi/180; x,y=pol2cart(rdirection,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);