万圣节快到了，一起用MATLAB绘制一个可爱的南瓜灯叭

效果及原理

效果如下：

调一下数据还能改成三角眼： 原理 南瓜主体函数从知友 [九章算法] 的一张图而来，大体是瓜身瓜柄分段函数，然后绕着z轴旋转一周得到曲面，我对数值做了微调，原图及原始数据：

这里我故意保留了网格让南瓜看起来有一点布娃娃的感觉，(大家也可以根据自己喜好改写，例如将’EdgeColor’设置为’none’并打个光啥的)

% 构造网格
[t,p]=meshgrid(linspace(0,2*pi,200),linspace(0,pi-.05,200));
% 二元分段函数
foutline=@(t,p)(p.14).*(cos(20.*t)+70).*(p.*(pi-p+.1)).^.2;

% 球面坐标转化为X,Y,Z
h=cos(p).*foutline(t,p);
R=sin(p).*foutline(t,p);
R=R;
X=cos(t).*R;
Y=sin(t).*R;

% 颜色矩阵构造
CMap=ones([size(t),3]);
tMap=ones(size(t));
c1=[253,158,3]./255;
c2=[76,103,86]./255;
for i=1:3
    tMap(:,:)=c1(i);
    tMap(p50&maskX130&maskY45&maskX90&maskY85&maskX90&maskY45&maskX90&maskY80&maskX90&maskY45&maskX90&maskY-35))=1;
% mask(maskX>45&maskX90&maskY80&maskX90&maskY45&maskX90&maskY80&maskX90&maskY