用Matlab 实现元胞自动机（网上收集、转载）技术 2010-03-13 13:04:28 阅读 121 评论 0 字号：大中小订阅 . file:life.m %% 初始化 m = 50; X = zeros(m,m); X(25,25) = 1; n = [m 1:m-1]; e = [2:m 1]; s = [2:m 1]; w = [m 1:m-1]; % 绘制初始图形 [i,j] = find(X); figure(gcf); plothandle = plot(i,j,'.', ... 'Color','blue', ... 'MarkerSize',12); axis([0 m+1 0 m+1]); %% 演化 for k = 1:50 %邻居数 N = X(n,:) + X(s,:) + X(:,e) + X(:,w) + ... X(n,e) + X(n,w) + X(s,e) + X(s,w); %概率阵 RAND = rand(m); %换代 X = X | (N.*RAND>0.99); %绘图 [i,j] = find(X); set(plothandle,'xdata',i,'ydata',j) drawnow pause(0.2) k end file 2: function sierpinski(n); % 使用元胞自动机生成 sierpinski 直角垫片 % Example: % sierpinski(256); % %算法见:孙博文,《分形算法与程序设计：用 Visual C++实现》 if nargin==0; n=256; end X=ones(n); X(1,n-1)=0; H=imshow(X,[]); set(gcf,'doublebuffer','on'); k=1; while k
