24小时热门版块排行榜

返回列表

wangyongjian

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 2103
红花: 2
帖子: 16
在线: 69.7小时
虫号: 1896594
注册: 2012-07-17
性别: GG
专业: 金属结构材料

[求助] MATLAB生命的游戏代码求助

各位大侠，我刚刚接触MATLAB，自己看了一下生命的游戏的代码。有几个困惑还请大侠们给解答一下，谢过了。
首先，
%set a few cells to one
cells(n/2,.25*n:.75*n) = 1;
cells(.25*n:.75*n,n/2) = 1;

cells = (rand(n,n))<.5 ;
%how long for each case to stability or simple oscillators
根据上面的代码可以知道，cells(n/2,.25*n:.75*n) = 1;cells(.25*n:.75*n,n/2) = 1是给cells矩阵赋值，但是我觉得cells = (rand(n,n))<.5 就已经可以完成任务了，而且加上cells = (rand(n,n))<.5 这一句后，cells(n/2,.25*n:.75*n) = 1;cells(.25*n:.75*n,n/2) = 1的效果也体现不出来啊？另外，cells = (rand(n,n))<.5 我觉得就只有赋值功能，但是看注释好像还有作为评判标准的意思（%how long for each case to stability or simple oscillators）？
第二，%nearest neighbor sum
      sum(x,y) = cells(x,y-1) + cells(x,y+1) +cells(x-1, y) + cells(x+1,y) + cells(x-1,y-1) + cells(x-1,y+1) +cells(3:n,y-1) + cells(x+1,y+1);
这个是计算一个方格周围近邻的个数， cells(3:n,y-1)是什么意思，看到很多计算近邻的代码都是这样写的，我觉得应该是 cells(x+1,y-1)啊。
第三， % The CA rule
cells = (sum==3) | (sum==2 & cells)。这句就不明白了，我觉得应该是cells（x，y） = (sum（x，y）==3) | (sum（x，y）==2 & cells)。。？

完整代码如下

%Conway's life with GUI
clf
clear all
%=============================================
%build the GUI
%define the plot button
plotbutton=uicontrol('style','pushbutton',...
'string','Run', ...
'fontsize',12, ...
'position',[100,400,50,20], ...
'callback', 'run=1;');
%define the stop button
erasebutton=uicontrol('style','pushbutton',...
'string','Stop', ...
'fontsize',12, ...
'position',[200,400,50,20], ...
'callback','freeze=1;');
%define the Quit button
quitbutton=uicontrol('style','pushbutton',...
'string','Quit', ...
'fontsize',12, ...
'position',[300,400,50,20], ...
'callback','stop=1;close;');
number = uicontrol('style','text', ...
'string','1', ...
'fontsize',12, ...
'position',[20,400,50,20]);

%=============================================
%CA setup
n=128;
%initialize the arrays
z = zeros(n,n);
cells = z;
sum = z;
%set a few cells to one
cells(n/2,.25*n:.75*n) = 1;
cells(.25*n:.75*n,n/2) = 1;
%cells(.5*n-1,.5*n-1)=1;
%cells(.5*n-2,.5*n-2)=1;
%cells(.5*n-3,.5*n-3)=1;
cells = (rand(n,n))<.5 ;
%how long for each case to stability or simple oscillators
%build an image and display it
imh = image(cat(3,cells,z,z));
set(imh, 'erasemode', 'none')
axis equal
axis tight
%index definition for cell update
x = 2:n-1;
y = 2:n-1;
%Main event loop
stop= 0; %wait for a quit button push
run = 0; %wait for a draw
freeze = 0; %wait for a freeze
while (stop==0)
if (run==1)
      %nearest neighbor sum
      sum(x,y) = cells(x,y-1) + cells(x,y+1) + ...
         cells(x-1, y) + cells(x+1,y) + ...
         cells(x-1,y-1) + cells(x-1,y+1) + ...
         cells(3:n,y-1) + cells(x+1,y+1);
      % The CA rule
      cells = (sum==3) | (sum==2 & cells);
      %draw the new image
      set(imh, 'cdata', cat(3,cells,z,z) )
      %update the step number diaplay
      stepnumber = 1 + str2num(get(number,'string'));
      set(number,'string',num2str(stepnumber))
end
if (freeze==1)
      run = 0;
      freeze = 0;
end
      drawnow  %need this in the loop for controls to work

end

回复此楼