sierpinski_zd4.asv

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

function Sierpinski_zd4(x,y,d,n); 
% Sierpinskidt地毯 
% 掉用格式 Sierpinski_zd4(x,y,d,n) 
% x 为正方形的第1个顶点的横坐标 
% x 为正方形的第1个顶点的纵坐标 
% d 为初始正方形边长 
% n为迭代次数 
% Examples: 
%       Sierpinski_zd4(0,0,1,4) 
if nargin==0;
    x=0;y=0;d=1;n=4;
end
hold on;
for p=1:n; 
   a1=[]; 
   b1=[]; 
   for q=1:length(x); 
       x1=x(q)+[0,d/3,2*d/3,0,2*d/3,0,d/3,2*d/3]; 
       y1=y(q)+[0,0,0,d/3,d/3,2*d/3,2*d/3,2*d/3]; 
       a1=[a1,x1]; 
       b1=[b1,y1]; 
   end 
   d=d/3; 
   x=a1; 
   y=b1; 
end 
for q=1:length(x); 
   fill(x(q)+[0,d,d,0,0],y(q)+[0,0,d,d,0],'b') 
   hold on 
end 
hold off 
axis off 
axis equal 
set(findobj(gcf,'type','patch'),'edgecolor','none')