Question

matlab中圆周元胞自动机的程序

Answer 1

%CA setup
n=200;
%initialize the arrays
z=zeros(n,n);
cells=z;
sum=z;
%set a few cells to one
cells(100,100)=1;
%index definition for cell update
x=2:1:199;
y=2:1:199;
((x.^2)+(y.^2))<=100^2;
%nearest neighbor sum
sum(x,y)=cells(x,y+1)+cells(x,y-1)+cells(x+1,y)+cells(x-1,y);
% The CA rule
cells=(sum==1);
plot(x,y);