基于元胞自动机模拟晶体生成附matlab代码

clear clc close all iterations = 100; beta = 0.95; %value for background level gamma = 0.0035; %constant added to receptive sites crystal_size = 101; %number of hex on top of center (note only use odd numbers) tic %just time stuff A = beta * ones(crystal_size*2 +1, crystal_size*2 + 1 ); %Array to store states A(crystal_size+1, crystal_size+1) = 1; %setting the center to ice x = -crystal_size:1:crystal_size; %array to store x positions of hexagons y = crystal_size:-1:-crystal_size; %array to store y positions of hexagons [X, Y] = meshgrid(x,y); %add offset to alternate rows and columns for i = 1:2:crystal_size*2 +1 X(i,:) = X(i,:) + 0.5; end Dis = sqrt(X.^2 + Y.^2); % drawing(A,X,Y,Dis,crystal_size); % drawnow for i = 1 : iterations [odd_row, even_row] = neighbor_array(A, beta); receptive_odd = sum((odd_row >= 1) , 3)>0; receptive_even = sum((even_row >= 1) , 3)>0; receptive_logic = zeros(size(A)); for k = 1:size(A,1) if mod(k,2) == 0 receptive_logic(k,:) = receptive_even(k,:); else receptive_logic(k,:) = receptive_odd(k,:); end end unreceptive_logic = receptive_logic==0; receptive = receptive_logic .* A; unreceptive = unreceptive_logic .* A; receptive_update = (receptive_logic * gamma) + receptive; [odd_row, even_row] = neighbor_array(unreceptive, beta); unreceptive_ne = zeros(size(odd_row)); for k = 1:size(A,1) if mod(k,2) == 0 unreceptive_ne(k,:,:) = even_row(k,:,:); else unreceptive_ne(k,:,:) = odd_row(k,:,:); end end averaged_neighbors = (unreceptive_ne(:,:,1)*0.5) + (unreceptive_ne(:,:,2)*(1/12)) + (unreceptive_ne(:,:,3)*(1/12)) ... + (unreceptive_ne(:,:,4)*(1/12)) + (unreceptive_ne(:,:,5)*(1/12)) + (unreceptive_ne(:,:,6)*(1/12)) + (unreceptive_ne(:,:,7)*(1/12)); final_update = averaged_neighbors + receptive_update; % if mod(i,10)==0 % clf % drawing(final_update,beta) % drawnow % end clf drawing(final_update,beta) drawnow A = final_update; end toc tic clf drawing(A,beta) drawnow toc