元胞自动机序.docVIP

元胞自动机程序 用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.*RAND0.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 kn; X(k+1,1:end-1)=xor(X(k,1:end-1),X(k,2:end)); X(k+1,n)=1; set(H,CData,X); pause(0.1); k=k+1; end file 3: function CA_sim_cloud; % 使用元胞自动机模拟地球卫星的云图 % % reference: % Piazza, E.; Cuccoli, F.; % Cellular Automata Simulation of Clouds in Satellite Images, % Geoscience and Remote Sensing Symposium, 2001. IGARSS 01. % IEEE 2001 International Volume 4, 9-13 July 2001 Page(s): % 1722 - 1724 vol.4 Digital Object Identifier 10.1109/IGARSS. % 2001.977050 time=888; % 程序执行步数 M=240; N=320; S=round(rand(M,N)*15); p=[1,2,1,6,6,1,2,1]; p=sum(tril(meshgrid(p)),2)/20; rand(state,0); SS=S; R=rand(M,N); G=R; B=R; C=cat(3,R,G,B); fig=figure; set(fig,DoubleBuffer,on); mov = avifile(example2.avi); cc=imshow(C,[]); set(gcf,Position,[13 355 157 194]) x1=(1:3)+round(M/2);y1=(1:3)+round(N/3); x2=(1:3)+round(M/3);y2=(1:3)+round(N/2); x3=(1:3)+round(M/1.5);y3=(1:3)+round(N/2); q=0; qq=15/4; while qtime; SS=zeros(M,N); for k=1:15; r=rand(M,N); % 生成几率r K=zeros(M+2,N+2); T=(S-k=0); % 粒子数矩阵 K(2:end-1,2:en