3D_game of life_3d元胞自动机_三维生命游戏_gameoflife_生命游戏

function life3D(action) play= 1; stop=-1; if nargin<1, action='initialize'; end; if strcmp(action,'initialize'), figNumber=figure( ... 'Name','Life3D: Conway''s Game of Life in 3-Dimesions', ... 'NumberTitle','off', ... 'DoubleBuffer','on', ... 'Visible','off', ... 'Color','white', ... 'BackingStore','off'); axes( ... 'Units','normalized', ... 'Position',[0.05 0.05 0.75 0.90], ... 'Visible','off', ... 'DrawMode','fast', ... 'Color','none',... 'NextPlot','add'); % 'NextPlot','replace' ); text(0,0,{'Press the "Start" button to see the Game of Life demo' 'Use the slider to change the number of initial cells.'}, ... 'HorizontalAlignment','center'); axis([-1 1 -1 1]); %=================================== % Information for all buttons labelColor=[0.8 0.8 0.8]; yInitPos=0.90; xPos=0.85; btnLen=0.10; btnWid=0.10; % Spacing between the button and the next command's label spacing=0.05; %==================================== % The CONSOLE frame frmBorder=0.02; yPos=0.05-frmBorder; frmPos=[xPos-frmBorder yPos btnLen+2*frmBorder 0.9+2*frmBorder]; h=uicontrol( ... 'Style','frame', ... 'Units','normalized', ... 'Position',frmPos, ... 'BackgroundColor',[0.50 0.50 0.50]); %==================================== % The START button btnNumber=1; yPos=0.90-(btnNumber-1)*(btnWid+spacing); labelStr='Start'; cmdStr='start'; callbackStr='life3D(''start'');'; % Generic button information btnPos=[xPos yPos-spacing btnLen btnWid]; startHndl=uicontrol( ... 'Style','pushbutton', ... 'Units','normalized', ... 'Position',btnPos, ... 'String',labelStr, ... 'Interruptible','on', ... 'Callback',callbackStr); %==================================== % The STOP button btnNumber=2; yPos=0.90-(btnNumber-1)*(btnWid+spacing); labelStr='Stop'; % Setting userdata to -1 (=stop) will stop the demo. callbackStr='set(gca,''Userdata'',-1)'; % Generic button information btnPos=[xPos yPos-spacing btnLen btnWid]; stopHndl=uicontrol( ... 'Style','pushbutton', ... 'Units','normalized', ... 'Position',btnPos, ... 'Enable','off', ... 'String',labelStr, ... 'Callback',callbackStr); %==================================== % The NumberOfCell Slider labelStr='# Cells'; callbackStr='h=findobj(gcf,''Tag'',''StartCells'');set(h,''Tooltip'',sprintf(''Initial Cells: %4d'',floor(get(h,''Value''))));'; infoHndl=uicontrol( ... 'Style','slider', ... 'Units','normalized', ... 'Position',[xPos btnWid+btnWid*2+0.05 btnLen/4 0.10*3], ... 'String',labelStr, ... 'SliderStep',[.01 .1], ... 'Max',100, ... 'Min',1, ... 'Value',20, ... 'Tooltip','Initial Cells: 20',... 'Tag','StartCells', ... 'Callback',callbackStr); %==================================== % The INFO button labelStr='Info'; callbackStr='life3D(''info'')'; infoHndl=uicontrol( ... 'Style','push', ... 'Units','normalized', ... 'Position',[xPos 0.20 btnLen 0.10], ... 'String',labelStr, ... 'Callback',callbackStr); %==================================== % The CLOSE button labelStr='Close'; callbackStr='close(gcf)'; closeHndl=uicontrol( ... 'Style','push', ... 'Units','normalized', ... 'Position',[xPos 0.05 btnLen 0.10], ... 'String',labelStr, ... 'Callback',callbackStr); % Uncover the figure hndlList=[startHndl stopHndl infoHndl closeHndl]; set(figNumber,'Visible','on', ... 'UserData',hndlList); view(3); StartCells = 20; m = 19; colormap(jet(m)); colorbar; elseif strcmp(action,'start'), objh=findobj(gcf,'Tag','StartCells'); StartCells = get(objh,'Value'); m = 19; cla; axHndl=gca; figNumber=gcf; hndlList=get(figNumber,'Userdata'); startHndl=hndlList(1); stopHndl=hndlList(2); infoHndl=hndlList(3); closeHndl=hndlList(4); set([startHndl closeHndl infoHndl],'Enable','off'); set(stopHndl,'Enable','on'); % ====== Start of Demo set(axHndl, ... 'UserData',play, ... 'DrawMode','fast', ... 'Visible','on'); box off X = zeros(m,m,m); p = -1:1; for count=1:StartCells, kx=floor(rand*(m-4))+2; ky=floor(rand*(m-4))+2; kz=floor(rand*(m-4))+2; X(kx+p,ky+p,kz+p)=(rand(3,3,3)>0.5); end; % The following statements plot the initial configuration. % The "find" function returns the indices of the nonzero elements. plothandle=[]; figure(gcf); hold on colour=jet(m); for k=1:m [i,j] = find(X(:,:,k)); % kd = sub2ind(size(X(i,j,k)),i,j) kd = k(ones(size(i))); if isempty(i) i = 1; j = 1; kd = NaN; end plothandle(k) = line(i,j,kd); set(plothandle(k),'linestyle','none',... 'Marker','.',... 'MarkerSize',20*2,... 'MarkerFaceColor',colour(k,:),... 'MarkerEdgeColor',colour(k,:),... 'EraseMode','normal'); end drawnow axis([0 m+1 0 m+1 0 m+1]); hold off % Whether cells stay alive, die, or generate new cells depends % upon how many of their eight possible neighbors are alive. % Here we generate index vectors for four of the eight neighbors. % We use periodic (torus) boundary conditions at the edges of the universe. n = [m 1:m-1]; e = [2:m 1]; s = [2:m 1]; w = [m 1:m-1]; u = [m 1:m-1]; d = [2:m 1]; rotate3d on while get(axHndl,'UserData')==play, % How many of eight+5+5 neighbors are alive. (only the ones that share at least one border % N = X(n,:,:) + X(s,:,:) + X(:,e,:) + X(:,w,:) + ... % X(n,e,:) + X(n,w,:) + X(s,e,:) + X(s,w,:); % 3D Version % Use Euclidean distance d1 = 1/1; % Adjacent cells (Share 4 vertix) d2 = 1/sqrt(2); % Diagonal cells (Share 2 vertix) d3 = 1/sqrt(3); % Double Diagonal cells (Share 1 vertix) % How many of 9+9+8 neighbors are alive. (Only the ones that share at least one vertix) N = X(n,:,:)*d1 + X(s,:,:)*d1 + X(:,e,:)*d1 + X(:,w,:)*d1 + ... X(n,e,:)*d2 + X(n,w,:)*d2 + X(s,e,:)*d2 + X(s,w,:)*d2 + ... X(n,:,u)*d2 + X(s,:,u)*d2 + X(:,e,u)*d2 + X(:,w,u)*d2 + ... X(n,:,d)*d2 + X(s,:,d)*d2 + X(:,e,d)*d2 + X(:,w,d)*d2 + ... X(n,e,u)*d3 + X(s,e,u)*d3 + X(s,w,u)*d3 + X(n,w,u)*d3 + ... X(n,e,d)*d3 + X(s,e,d)*d3 + X(s,w,d)*d3 + X(n,w,d)*d3 + ... X(:,:,u)*d1 + X(:,:,d)*d1; % A live cell with two live neighbors, or any cell with three % neigbhors, is alive at the next time step. Xold = X; CL = 3.0; % Minimun # cells to create a new one CH = 4.0; % Maximum # cells to create a new one KL = 4.0; % Minimun # cells to kill one X = ( ( ((N >= CL) & (N <= CH)) ) &... % Create cells ~( X & ((N > KL) )) ) ; % Kill cells by overpopulation Xsum = sum(X(:)); if (Xsum>prod(size(X))/8); X = ( X & (rand(size(X))>0.1) ); % Expontaneous cell annihilation end if Xold==X % if no change the exit break; end % Update plot. for k=1:m [i,j] = find(X(:,:,k)); set(plothandle(k),'xdata',i,'ydata',j,'zdata',k(ones(size(i)))); end set(plothand