cell1.m

元胞自动机

% Percolation Cluster
clf
clear all
threshold = .63;
%
ax = axes('units','pixels','position',[1 1 500 400],'color','k');
text('units', 'pixels', 'position', [50,255,0],...
'string','GalaXy','color','w','fontname','helvetica','fontsize',100)
text('units', 'pixels', 'position', [120,120,0],...
'string','Cellular Automata','color','w','fontname','helvetica','fontsize',30)
initial = getframe(gca);
[a,b,c]=size(initial.cdata);
z=zeros(a,b);
cells = double(initial.cdata(:,:,1)==255);
visit = z ;
sum = z;
imh = image(cat(3,z,cells,z));
set(imh, 'erasemode', 'none')
%return
for i=1:100
sum(2:a-1,2:b-1) = cells(2:a-1,1:b-2) + cells(2:a-1,3:b) + ...
cells(1:a-2, 2:b-1) + cells(3:a,2:b-1) + ...
cells(1:a-2,1:b-2) + cells(1:a-2,3:b) + ...
cells(3:a,1:b-2) + cells(3:a,3:b);
pick = rand(a,b);
%edges only
%cells = (cells & (sum<8)) | ((sum>=1) & (pick>=threshold) & (visit==0)) ;
cells = cells | ((sum>=1) & (pick>=threshold) & (visit==0)) ;
visit = (sum>=1) ;%& (pick<threshold) ;
set(imh, 'cdata', cat(3,z,cells,z) )
drawnow
end
return
figure(2)
image(cat(3,z,cells,z))