sierpinski3_by_ca.m

来自「元胞自动机的matlab程序。 元胞自动机（CA）是一种用来仿真局部规则和局部」· M 代码 · 共 26 行

function sierpinski3_by_CA(n); 
% 使用元胞自动机生成sierpinski直角垫片 
if nargin==0; 
   n=256; 
end 
X=zeros(n); 
X(1,round(n/2))=1; 
H=imshow(X,[]); 
set(gcf,'doublebuffer','on'); 
k=1; 
while k<round(n/2); 
   X(k+1,2:end-1)=and(xor(X(k,1:end-2),X(k,3:end)),... 
       ~X(k,2:end-1)); 
   set(H,'CData',1-X); 
   pause(0.05); 
   k=k+1; 
end 
nm=round(n/2); 
k=1; 
while k<nm; 
   X(nm+k,1:end)=X(nm-k,1:end); 
   set(H,'CData',1-X); 
   pause(0.05); 
   k=k+1; 
end