《数学建模案例MATLAB实用程序百例》由会员分享，可在线阅读，更多相关《数学建模案例MATLAB实用程序百例（206页珍藏版）》请在金锄头文库上搜索。
1、file:/E|/Document/发展篇/Mh1=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 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(自
2、变量X);ylabel(函数值Y);title(三角函数);file:/E|/Document/发展篇/Mx=-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,hli
3、ne1,hline2=plotyy(x,y1,x,y2,semilogy,plot);axes(haxes(1)ylabel(semilog plot);file:/E|/Document/发展篇/M实例5：单个轴窗口显示多个图形function shili05h0=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 equalsubplot(2,2,2)z=sin(x)-cos(y);plot(
4、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,.file:/E|/Document/发展篇/Mt=0:pi/10:2*pi;h=plot(t,sin(t);xlabel(t=0到2pi,fontsize,16);ylabel(sin(t),fontsize,16);t
5、itle(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),.leftarrow最大值=,num2str(y(imax),.fontsize,16) 实例7：条形图形function shili07h0=figure(toolbar,none,.position,200 150
6、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值);file:/E|/Document/发展篇/Mh2=axes(position,get(h1,position);plot(t,tiao2,linewidth,3)set(h2,yaxislocation,right,color,none,xticklabel,)实例8：区域图形function shili08h0=figure(too
7、lbar,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
8、,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,.file:/E|/Document/发展篇/Mylabel(销售量,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=s
9、um(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);file:/E|/Document/发展篇/Mnewext=cat(1,val2:);offset=sign(oldext(:,1).*(newext(:,3)-oldext(:,3)/2;pos=get(textobjs,position);textpos=cat(1,po
10、s:);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(glabel,fontsize
11、,16)xlabel(t=0:10,fontsize,16)axis(0 10 -1.2 1.2)file:/E|/Document/发展篇/Mx=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);hold offh3=h1(1);h2;legend(h3,y1+y2,y1=sin(x),y2=cos(x)xlabel(自变量X);ylabel(函数值Y);title(正弦函数与余弦函数的线性组合); 实例12：罗盘图function shili12h0=figure(toolbar
12、,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;file:/E|/Document/发展篇/Mx,y=pol2cart(rdirection,windpower);compass(x,y);desc=风向和风力,北京气象台,10月1日0：00到,10月1日12：00;gtext(desc)实例13：轮廓图function shili13h0=figure(toolbar,n
13、one,.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)axis equalxlabel(实部,fontsize,16);ylabel(虚部,fontsize,16);h=polar(0 2*pi,0 1);delete(h)hold oncontour(x,y,abs(f),20)file:/E|/Document/发展篇/Maxis(0 10 0 10
14、);hold onx=;y=;n=0;disp(单击鼠标左键点取需要的点);disp(单击鼠标右键点取最后一个点);but=1;while but=1xi,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 off 实例14：交互式图形file:/E|/Document/发展篇/Maxis(0 10 0 10);hold onx=;y=;n=0;disp(单击鼠标左键点取需要的点);disp(单击鼠标右键点取最后一个点);but=1;while but=1xi,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 off
