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('双曲