元胞自动机的Matlab实现.pptVIP

元胞自动机的Matlab实现  MATLAB的编程考虑 1 图形用户界面 2 初始化（开辟空间初始值） 3 图形化（元胞空间） 4 核心代码（变换规则） 1 图形用户界面 %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]); 对象注释：字符串 2 初始化 矩阵和图像可以相互转化。 所以初始条件可以是矩阵，也可以是图形。 Cell：创建元胞数组的语法格式 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 3 图形化 矩阵和图像可以相互转化，所以矩阵的显示是可以直接实现的。如果矩阵cells的所有元素只包含两种状态且矩阵Z含有零，那么用image函数来显示cat命令建的RGB图像，并且能够返回句柄。 %build an image and display it imh = image(cat(3,cells,z,z)); set(imh, erasemode, none) axis equal axis tight 等长刻度坐标轴 设置当前数据为坐标范围 4 核心代码（变换规则） %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 drawno