L=55;L2=15;
for r=0:25
for a=0:(pi*5/180):pi;
for s=0:(pi*5/180):2*pi;
x=L2.*sin(s).*sin(a)+L.*sin(s).*(1-cos(a))./a;
y=L2.*cos(s).*sin(a)+L.*cos(s).*(1-cos(a))./a;
z=L2.*cos(a)+L.*sin(a)./a+r;
plot3(x,y,z);
hold on;
end
end
end
L=55;L2=15;
for r=0:25
for a=0:(pi*5/180):pi;
for s=0:(pi*5/180):2*pi;
x=L2.*sin(s).*sin(a)+L.*sin(s).*(1-cos(a))./a;
y=L2.*cos(s).*sin(a)+L.*cos(s).*(1-cos(a))./a;
z=L2.*cos(a)+L.*sin(a)./a+r;
plot3(x,y,z);
hold on;
end
end
end
转载于:https://my.oschina.net/donngchao/blog/746564