sierpinski_zd.m

来自「包括Koch曲线、Levy 曲线、分形树、Sierpinski三角形」· M 代码 · 共 24 行

function Sierpinski_zd(n); 
% 中点法生成Sierpinski地毯 
% Examples: 
%      Sierpinski_zd(200) 
a=0; 
b=100; 
c=50+i*sqrt(3)*50; 
z=15+15i; 
zp=[]; 
for k=1:n; 
   N=randperm(100); 
   for p=N; 
       if p<33; 
           z=(z+b)/2; 
       elseif p<66 
           z=(z+c)/2; 
       else 
           z=(z+a)/2; 
       end 
       zp=[zp,z]; 
   end 
end 
plot(zp,'.','markersize',2)