第2个回答 2010-07-08
h=2*pi;
[r,theta]=meshgrid(linspace(0,1,50),linspace(0,4*pi,500));
x=r.*cos(theta);
y=r.*sin(theta);
z=h*theta/2/pi;
surf(x,y,z);
shading interp本回答被提问者采纳
第3个回答 2010-07-07
clc
clf
r=[0:0.001:1];
h=2*pi;
theta=[0:2*pi/1000:3*pi]';
x(length(theta),length(r))=0;
y(length(theta),length(r))=0;
z(length(theta),length(r))=0;
for k=1:length(r)
x(:,k)=r(k)*cos(theta);
y(:,k)=r(k)*sin(theta);
z(:,k)=h*theta/2/pi;
end
surf(x,y,z);
shading interp;
box off
grid off
colormap('default')
axis on
view(17,55)