如何用MATLAB繪制出sierpinski金字塔？

x=cosd([90?210?330]');

y=sind([90?210?330]');

tri=[1?2?3];

tnum=1;dnum=3;

for?ii=2:n

tri2=zeros(3*tnum,3);

nn=1;

for?jj=1:tnum

xnew=(x(tri(jj,[1?2?3]))+x(tri(jj,[2?3?1])))/2;

ynew=(y(tri(jj,[1?2?3]))+y(tri(jj,[2?3?1])))/2;

x=[x;xnew];

y=[y;ynew];

tri2(nn,:)=[tri(jj,1),dnum+1,dnum+3];

tri2(nn+1,:)=[dnum+1,tri(jj,2),dnum+2];

tri2(nn+2,:)=[dnum+3,dnum+2,tri(jj,3)];

nn=nn+3;

dnum=dnum+3;

end

tri=tri2;

tnum=3*tnum;

end

triplot(tri,x,y);

axis?image;

axis?off;