(CA学习资料)美国康奈尔大学BioNB441在Matlab中的元胞自动机.doc

... WORD格式可编辑版 网站翻译内容 美国康奈尔大学BioNB441在Matlab中的元胞自动机 介绍 元胞自动机（CA）是用于计算计划利用当地的规则和本地通信。普遍CA定义一个网格，网格上的每个点代表一个有限数量的状态中的细胞。过渡规则同时应用到每一个细胞。?典型的转换规则依赖于细胞和它的（4个或8个）近邻的状态，虽然临近的细胞也同样使用。 CA的应用在并行计算研究，物理模拟和生物模拟。这个页面将考虑如何写出高效的MATLAB代码的CA的实施和看一些有趣的规则。 Matlab代码注意事项 以下注意事项将说明使用Matlab程序计算康威的生命。部分的代码的解释如下。 矩阵和图像可以被转换为一个另一个，所以显示是为straighforward的的。如果阵列细胞包含二进制的每一个电池单体的状态和数组Z含有零，那么cat命令建立一个RGB图像，显示图像命令。图像命令也返回一个句柄。 imh = image(cat(3,cells,z,z)); set(imh, erasemode, none) axis equal axis tight 矩阵和图像可以被转换为一个另一个，那么初始条件，可以计算矩阵或图形命令。下面的代码生成一个元素为零的数组，细胞状态初始化为零，然后使细胞的交叉状态=1的数组的中心。使用（渗流群集）下面的例子中的一个图形命令来初始化CA阵列。 z = zeros(n,n); cells = z; cells(n/2,.25*n:.75*n) = 1; cells(.25*n:.75*n,n/2) = 1; Matlab代码进行矢量化，以减少开销。下面的两个语句计算的总和近邻的计算CA规则。代码利用Matlab的非常灵活的索引指定的近邻。 x = 2:n-1; y = 2:n-1; 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(x+1,y-1) + cells(x+1,y+1); cells = (sum==3) | (sum==2 cells); 添加一个简单的GUI是很容易的。在这个例子中，三个按钮和一个文本字段implmented。三个按钮允许用户运行，停止，或选择退出。文本字段中显示模拟执行的步骤数。 %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）初始化，程序进入一个循环的测试中设置的每个按钮的回调函数的变量的状态下降。就目前而言，只是看在while循环，if语句的嵌套结构。程序循环，直到按下退出按钮。另外两个按钮导致当推一个if语句来执行。 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) + c