t=0:pi/20:2*pi;
x=sin(t)+1;
y=cos(t);
z=linspace(-4,4,length(t));
X=meshgrid(x);Y=meshgrid(y);Z=[meshgrid(z)]';
mesh(X,Y,Z)%第一个圆柱面
xlabel('x')
ylabel('y')
zlabel('z')
axis equal
t=0:pi/20:2*pi;
x=sin(t)+1;
y=cos(t);
z=linspace(-4,4,length(t));
X=meshgrid(x);Y=meshgrid(y);Z=[meshgrid(z)]';
mesh(X,Y,Z)%第一个圆柱面
xlabel('x')
ylabel('y')
zlabel('z')
axis equal