Matlab 绘制螺壳

 

 

t=linspace(0,2*pi,512);
[u,v]=meshgrid(t);
a=-0.4;b=0.5;c=0.1;
n=2;
x=(a*(1-v/(2*pi)).*(1+cos(u))+c).*cos(n*v);
y=(a*(1-v/(2*pi)).*(1+cos(u))+c).*sin(n*v);
z=b*v/(2*pi)+a*(1-v/(2*pi)).*sin(u);
surf(x,y,z,flip(y)) %第四个参数为颜色
axis off
axis equal
shading interp  %插值
material shiny;
colormap(hsv(1000)) %色条数量
lighting phong; %光照
camlight('left','infinite');view([-160 25])

　　将代码里的flip去掉可以得到更好看的颜色