2018年Matlab画函数图像.doc

Matlab画函数图像

用 Matlab 画函数图像

一、螺旋线

1. 静态螺旋线

a=0:0.1:20*pi;

h=plot3(a.*cos(a),a.*sin(a),2.*a,'b','linewidth',2);

axis([-50,50,-50,50,0,150]);

grid on

set(h,'erasemode','none','markersize',22);

xlabel('x轴 ');ylabel('y轴 ');zlabel('z轴 ');

title('静态螺旋线 ');

2. 动态螺旋线

t=0:0.1:10*pi;

i=1;

h=plot3(sin(t(i)),cos(t(i)),t(i),'*','erasemode','none');

grid on

axis([-2 2 -2 2 0 35])

for i=2:length(t)

set(h,'xdata',sin(t(i)),'ydata',cos(t(i)),'zdata',t(i)); drawnow

pause(0.01)

end

title('动态螺旋线 ');

(图略 )

3. 圆柱螺旋线

t=0:0.1:10*pi;

x=r.*cos(t);

y=r.*sin(t);

z=t;

plot3(x,y,z,'h','linewidth',2);

grid on

axis('square')

xlabel('x轴 ');ylabel('y轴 ');zlabel('z轴 ');

title('圆柱螺旋线

')

二、旋转抛物面

b=0:0.2:2*pi;

[X,Y]=meshgrid(-6:0.1:6);

Z=(X.^2+Y.^2)./4;

meshc(X,Y,Z);

axis('square')

xlabel('x轴 ');ylabel('y轴 ');zlabel('z轴 ');

title('旋转抛物面 ')

或直接用:ezsurfc('(X.^2+Y.^2)./4')

三、椭圆柱面

load clown

ezsurf('(2*cos(u))','4*sin(u)','v',[0,2*pi,0,2*pi])

view(-105,40) %视角处理

shading interp %灯光处理

colormap(map) %颜色处理

grid on %添加网格线

axis equal %使 x,y 轴比例一致

xlabel('x轴 ');ylabel('y轴 ');zlabel('z轴 '); %添加坐标轴说明 title('椭圆柱面 ') %添加标题

四、椭圆抛物面

b=0:0.2:2*pi;

[X,Y]=meshgrid(-6:0.1:6);

Z=X.^2./9+Y.^2./4;

meshc(X,Y,Z);

axis('square')

xlabel('x轴 ');ylabel('y轴 ');zlabel('z轴 '); title('椭圆抛物面 ')

或直接用:ezsurfc('X.^2./9+Y.^2./4')

五、 ' 双叶双曲面

ezsurf('8*tan(u)*cos(v)','8.*tan(u)*sin(v)','2.*sec(u)',[-pi./2,3*pi./2,0,2*pi]) axis equal

grid on

axis square

xlabel('x轴 ');ylabel('y轴 ');zlabel('z轴 ');

title('双叶双曲面 ')

六、双曲柱面

load clown

ezsurf('2*sec(u)','2*tan(u)','v',[-pi/2,pi/2,-3*pi,3*pi]) hold on %在原来的图上继续作图

ezsurf('2*sec(u)','2*tan(u)','v',[pi/2,3*pi/2,-3*pi,3*pi]) colormap(map)

shading interp

view(-15,30)

axis equal

grid on

axis equal

xlabel('x轴 ');ylabel('y轴 ');zlabel('z轴 ');

title('双曲