www.gusucode.com > matlab通信工程仿真源码（张德丰等编著）程序书籍 > matlab_code/matlab通信工程仿真源码（张德丰等编著）/第4章/li4_34.m

    x=linspace(0,2*pi,100);
R=1-exp(-x/2).*cos(4*x);
[X,Y,Z]=cylinder(R,60);
I=find(X<0 & Y<0);
Z(I)=NaN;
surf(X,Y,Z);
colormap(spring);
shading interp
light('position',[-3,-1,3],'style','local')
material([0.5,0.4,0.3,10,0.3])