数理学院 | 添加收藏
您当前的位置:首页 > 建模提高

matlab经典算法的程序——绘图

 

 

程序代码:

clear
clc
subplot(1,2,1);
t = 0:0.01:2*pi;
x = cos(2*t).*(cos(t).^2);
y = sin(2*t).*(sin(t).^2);
comet(x,y)
subplot(1,2,2);
t = -10*pi:pi/250:10*pi;
comet3((cos(2*t).^2).*sin(t),(sin(2*t).^2).*cos(t),t)

clear
clc
subplot(2,2,1)
contour3(peaks,50); %画出曲面在三度空间中的等高线
axis([-inf inf -inf inf -inf inf]);
subplot(2,2,2)
contour(peaks, 50); %画出曲面等高线在XY平面的投影
subplot(2,2,3)
t=linspace(0,20*pi, 501);
plot3(t.*sin(t), t.*cos(t), t);% 画出三度空间中的曲线
subplot(2,2,4)
plot3(t.*sin(t), t.*cos(t), t, t.*sin(t), t.*cos(t), -t);% 同时画出两条三度空间中的曲线


clear
clc
t=(0:0.02:2)*pi;
x=sin(t);
y=cos(t);
z=cos(2*t);
plot3(x,y,z,'b-',x,y,z,'bd')
view([-82,58]);
box on
legend('链','宝石');

clear
clc
subplot(3,3,1)
ezplot('cos(x)')
subplot(3,3,2)
ezplot('cos(x)', [0, pi])
subplot(3,3,3)
ezplot('1/y-log(y)+log(-1+y)+x - 1')
subplot(3,3,4)
ezplot('x^2 - y^2 - 1')
subplot(3,3,5)
ezplot('x^2 + y^2 - 1',[-1.25,1.25]);
axis equal
subplot(3,3,6)
ezplot('x^3 + y^3 - 5*x*y + 1/5',[-3,3])
subplot(3,3,7)
ezplot('x^3 + 2*x^2 - 3*x + 5 - y^2')
subplot(3,3,8)
ezplot('sin(t)','cos(t)')
subplot(3,3,9)
ezplot('sin(3*t)*cos(t)','sin(3*t)*sin(t)',[0,pi])


clear
clc
subplot(2,2,1), fplot('humps',[0 1])
subplot(2,2,2), fplot('abs(exp(-j*x*(0:9))*ones(10,1))',[0 2*pi])
subplot(2,2,3), fplot('[tan(x),sin(x),cos(x)]',2*pi*[-1 1 -1 1])
subplot(2,2,4), fplot('sin(1 ./ x)', [0.01 0.1],1e-3)


clear
clc
[X0,Y0,Z0]=sphere(30); %产生单位球面的三维坐标
X=2*X0;Y=2*Y0;Z=2*Z0; ? %产生半径为2的球面的三维坐标
clf
subplot(1,2,1);
surf(X0,Y0,Z0); ? %画单位球面
shading interp %采用插补明暗处理
hold on,mesh(X,Y,Z),colormap(hot),hold off %采用hot色图
hidden off ? %产生透视效果
axis equal,axis off ? %不显示坐标轴
title('透视图')
subplot(1,2,2);
surf(X0,Y0,Z0); ? %画单位球面
shading interp %采用插补明暗处理
hold on,mesh(X,Y,Z),colormap(hot),hold off %采用hot色图
hidden on ? %产生消隐效果
axis equal,axis off ? %不显示坐标轴
title('消隐图')


clear
clc
[x,y,z] =peaks;
subplot(2,3,1);
meshz(x,y,z); %曲面加上围裙,即给出曲面和零平面
axis([-inf inf -inf inf -inf inf]);
subplot(2,3,2);
waterfall(x,y,z); %在x方向产生水流效果
axis([-inf inf -inf inf -inf inf]);
subplot(2,3,3);
meshc(x,y,z); %同时画出网状图与等高线
axis([-inf inf -inf inf -inf inf]);
subplot(2,3,4);
surfc(x,y,z); %同时画出曲面图与等高线
axis([-inf inf -inf inf -inf inf]);
subplot(2,3,5)
surfl(x,y,z); %给出带光照效果的彩色表面图
axis([-inf inf -inf inf -inf inf]);
subplot(2,3,6)
contourf(x,y,z);
axis([-inf inf -inf inf -inf inf]);


clear
clc
x=-2:0.01:2;
[x,y]=meshgrid(x,x); %x和y都是401x401的矩阵
r=sqrt(x.^2+x.^2)+eps;
z=sinc(r);
subplot(2,1,1);
mesh(z);
subplot(2,1,2);
surf(x,y,z);


clear
clc
x=[11.4 23.5 35.4 15.6];
explode=zeros(size(x));
[c,offset]=min(x);
explode(offset)=c;
pie(x,explode)


clear
clc
t=0:0.001:10;
y=sin(t);
% plot(t,y);
Y=sin(10*t);
c=y.*Y;
plot(t,y,'r:',t,c,'b')


t=-10:1:10;
subplot(2,2,1);
bar(t,cos(t));
subplot(2,2,2);
compass(t,cos(t));
subplot(2,2,3);
rose(t,cos(t));
subplot(2,2,4);
fill(t,cos(t),'b');

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 


 


  友情链接
网站首页 | 关于我们 | 我要投稿 | 联系我们 | 数模论坛
Copyright@2010 苏州科技学院数理学院 All Rights Reserved.