2012年9月全国数学建模竞赛必备的MATLAB绘图程序集锦.docVIP

2012数学建模必备的MATLAB绘图程序集锦 Y=sin(x)在（0，2）上的图像：（ MATLAB语句：x=[0:0.01:2*pi];y=sin(x);plot(y) ） [x=0:pi/100:2*pi; y1=sin(x);y2=sin(x-0.25);y3=sin(x-0.5);plot(x,y1,b-.,x,y2,r*--,x,y3,g-.)] 【x=0:i/100:2*pi; y1=sin(x);y2=cos(x);y3=cos(x).^3+sin(x).^3;plot(x,y1,b*--,x,y2,rx-.,x,y3,ko-) 】 【 x=-2:0.01:2;y=exp(x);plot(x,y,b),... xlabel(自变量x), ylabel(因变量y),title(指数函数y=e^2的图像) 】 2、滤波器数学模型的三维曲面绘制：（见教材39页） 上图程序：[x,y]=meshgrid(0:31);n=2;D0=200;... D=sqrt((x-16).^2+(y-16).^2);z=1./(1+D.^(2*n)/D0);mesh(x,y,z),axis([0,31,0,31,0,1]) 程序2：[x,y]=meshgrid(0:31);n=2;D0=200;... D=sqrt((x-16).^2+(y-16).^2);z=1./(1+D.^(2*n)/D0);surf(x,y,z),axis([0,31,0,31,0,1]) Waterfall()绘制如上： contour3()绘制的三维等高线如下： 3.显函数方程y=sin(tan(x))-tan(sin(x))在x:[-pi,pi]上的图像：（P31） 【 x=-pi:0.05:pi;y=sin(tan(x))-tan(sin(x));plot(x,y) 】 4、P34页【ezplot(x^2*sin(x+y^2)+y^2*exp(x+y)+5*cos(x^2+y))】 【 ezplot(x^2*sin(x+y^2)+y^2*exp(x+y)+5*cos(x^2+y),[-10,10]) 】 5、绘制二元函数的图像：(P40) 【[x,y]=meshgrid(-2:0.1:2);z=1./(sqrt((1-x).^2+y.^2))+1./(sqrt((1+x).^2+y.^2));surf(x,y,z),shading flat】 6、syms x; fun=sin(x)*exp(-x)/(x+1)^3; % fun 应该为字符串 maple(with(numtheory):); f=maple([cfe:=cfrac( fun ,x,10)]) n=collect(maple(nthnumer,cfe,8),x); % 分子多项式合并同类项 d=collect(maple(nthdenom,cfe,8),x); [n,d]=numden(n/d); G=n/d; latex(G) n=collect(maple(nthnumer,cfe,10),x); d=collect(maple(nthdenom,cfe,10),x); [n,d]=numden(n/d); G1=n/d; latex(G1) % 显示从略 ezplot(fun,[0,2]), hold on; ezplot(G,[0,2]); ezplot(G1,[0,2]) figure; ezplot(fun,[0,5]), hold on; ezplot(G,[0,5]); ezplot(G1,[0,5]) 7、【 x=0:0.002:2; y=exp(-x).*sin(5*x); r=0.05*randn(size(x)); y1=y+r; plot(x,y1) 】 8、滤波效果图：【 b=1.2296e-6*conv([1 4 6 4 1],[1 3 3 1]); a=conv([1,-0.7265],... conv([1,-1.488,0.5644],conv([1,-1.595,0.6769],[1,-1.78,0.8713]))); x=0:0.002:2; y=exp(-x).*sin(5*x); r=0.05*randn(size(x)); y1=y+r; [h,w]=freqz(b,a,100); plot(w,abs(h)) % 放大倍数绘制 figure; y2=filter(b,a,y1); plot(x,y1,x,y2) % 滤波效果 】 9、【 x=[-5:.02:5]; y1=[]; y2=[]; mu1=[-1,0,0,0,1];