loukong.m

来自「matlab宝典,电子工业出版社,包含该书的源码4-6章」· M 代码 · 共 16 行

% 脚本文件loukong.m
t=linspace(0,2*pi,100);
%旋转母线
r=1-exp(-t/2).*cos(4*t); 
%产生旋转柱面的数据
[x,y,z]=cylinder(r,60);
%确定第四象限的数据下标
ii=find(x<0&y<0);
%进行图形剪切
z(ii)=NaN;
surf(x,y,z);
colormap(hsv),shading interp
%设置光源
light('position',[-3,1,3],'style','local');
%设置表面反射
material([0.5,0.4,0.3,10,0.3])