元胞自动机程序.txt资料文档下载

1、end）; X（k+1,n）=1; set（H,CData,X）; pause（0.1）; k=k+1;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

2、（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,DoubleBuffermov = avifile（example2.avicc=imshow（C,）;Po

3、sition,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:q=0;qq=15/4;while q=0）; % K（2:end-1,2:end-1）=T; SS=K（1:end-2,1:end-2）.*（rp（1）+. K（1:end-2,2:end-1）.*（r=p（1）+. end-2,3:end）.*（r=p（7）+SS; end S=SS; %SS S（S15）=15; S（x1,y1）

4、=15; S（x2,y2）=15; S（x3,y3）=15; % G=（Sqq）; R=（Sqq & S0.99 | imag（pz）0.99 | real（pz） A=-A-pi;set（po,XData,real（pz）,YData,imag（pz）;pause（0.2）;po1=uicontrol（gcf,stylepushunitnormalizedposition,0.74,0.87,0.1,0.08,.stringstart,18,callback,）;set（po1,k=k;if k=1; stopelse set（po1,ses）;file 5:% DLA% clc;clear

5、;close all;S=ones（40,100）;% state matrixS（end,:）=0; % initial sttaeSs=zeros（size（S）+1,0）; % top line is origin of particleSs（2:end,:）=S; % showing matrixN=size（S,2）;II=imagesc（Ss）;axis equal;colormap（gray）while sum（1-S（1,:）0.5; y=1; x=round（rand*N-1+1）; % random position D=0; while D % random travel

6、 r=rand; if abs（x-1）0.1; SL=1; else SL=S（y,x-1）; end if abs（x-N） SR=1; SR=S（y,x+1）; if SL+SR+S（y+1,x）2.5; % check the neighbor: left, right, under D=1; S（y,x）=0; % stop in the position if r=1/3; % travel randomly x=x-1; elseif r=2/3; x=x+1; y=y+1; Ss（2: if xN+0.5; % out of the range else Ss（y,x）=0;

7、% to show the moving particle set（II,Ss）; % to show pause（0.1）; endfile 6:function sands（N）;% % % % :（） close allDoublebufferD=ones（N）;D1=D;X,Y=meshgrid（1:N）;Z=2*X-Y;p=fix（9.5*N/21）;D（Zp & ZB（fix（end/2）+1:B（:,fix（end/2）+1:end）=1;B=min（B,fliplr（B）;B（1:3,:% figure;% imshow（B,）% mov = avifile（example29

8、00.avikk=3800;for k=1:kk; B,Nu=duisha（D,B）; Bk=min（D,B）; imshow（Bk,） Ha（k）=Nu;% F = getframe（gca）;% mov = addframe（mov,F）;% mov = close（mov）;plot（Ha）function Y,Nu=duisha（D,B）;Dq=10*（1-D）;Bg=1-B+Dq;% S=zeros（size（B）;S1=S;S2=S;S（2:）=Bg（2:）-Bg（1:end-1,:S1（S=-1）=1;S2（1:）=S1（2:Bg（S1=1）=Bg（S1=1）;Bg（S2=1）=

9、Bg（S2=1）;% clear Sclear S1clear S2S（1:end-1）=Bg（1:end-1）+Bg（2:end-1）-Bg（2:end-2）;S1（S=2）=1;S2（2:end-2）=S1（1:end-1）;Bg（S1=1）=0;Bg（S2=1）=1;end）;end）=S1（1:Y=（1-Bg）.*D;Nu=prod（size（find（Y=0）;file 7:time=500;S=zeros（M,N）;R=1-S;G=S;B=S;x=round（M/2）;y=（1: SS=zeros（M,N）; for k=1: r=rand（M,N）; % r K=zeros（M+

10、2,N+2）; T=（S-k K（2: SS=K（1: K（1: K（2: K（3: end S=SS; S（S S（x,y）=15; G=（S5）; R=（S5 & C=double（cat（3,R,G,B）; set（cc, q=q+1; pause（0.2）; title（ Nu（q）=sum（S（1:file 8: （Came of Life）J. H. Conway2060（GarclnerM.19701971） 01;򻮷񣬺򻮷（）乹򵥡:（1）򻮷 （2）010l （3）8Moore （4） （）:- 

11、27283;򵥡򼸸򱣳󡭡 （Glider） matlab % % function Game_of_Life（n） % % Game_of_Life（100）; n=100;B=round（rand（n+2）;Z=B（2:H=imshow（Z,）;,241 132 560 420） xlabel（Please press space key and stop this program!,. ,12,colortitle（Game of lifeb s=get（gcf,currentkey if strcmp（s,space clc;k

12、=0; A=sumfun（B）; X=zeros（n）; X（Z=1 & （A=2 | A=3）=1; X（Z=0 & A=3）=1; B（2:end-1）=X; Z=X; set（H,1-X）; pause（0.5）;figure（gcf）;function S=sumfun（B）;% 8 S=B（1:end-1）+. B（3: B（2:end-2）+. end）+. B（1:file 9:x（i-1,i,i+1,t）0 0 11 0 0x（i,t+1）1x（i,t+1）0 : 0 0 1 1 0 0 1 1 0 function sierpinski3_by_CA（n）;% sierpin

13、ski3_by_CA（256）;X=zeros（n）;X（1,round（n/2）=1;round（n/2）; X（k+1,2:end-1）=and（xor（X（k,1:end-2）,X（k,3:end）,. X（k,2:end-1）; pause（0.05）;nm=round（n/2）;nm; X（nm+k,1:end）=X（nm-k,1:file 10:file 11:function Clumsychild（c1,m）;% 4 % Clumsychild（0.1,0.1）;zjliu2001 N=100;A=randperm（N2）;S=zeros（N）;S（A（1:3000）=3; %

14、 30%3 S（A（3001:5000）=1; % 20%1 S（A（5001:6000）=2; % 10%2 clear A;,159 42 567 427）;imagesc（1:4） YAxisLocationright,0.8,0.1,0.1,0.8） ,0.84,0.12,0.1,0.8） xtickytick,1:4）;yticklabel,num2str（0:3）;axes（,0.06,0.12,0.7,0.8）;H=imagesc（S）;,159 42 485 427）;title（c1=,num2str（c1）, m=,num2str（m）;while k S=evolveme

15、nt（S,c1,m）;,S）;function S=evolvement（S,c1,m）;P=zeros（size（S）;Da=rand（size（S）;Da（Da1-c1）=1;Da（Da1-c1）=0;P（S=1 | S=2）=1;R=round（rand（size（S）+1）;P=P.*R.*Da;V=round（rand（size（S）*3）+1;V=V.*P; % V: % 1 - up % 2 - down % 3 - left % 4 - right V（1,V（1,1:end）=1）=0;V（end,V（end,1:end）=2）=0;V（V（1:end,1）=3,1）=0;e

16、nd,end）=4,end）=0;% x,y=find（V=1）;DD=zeros（size（S）;DD（x-1,y）=P（x-1,y）;S（S=0 | S=2 & DD=1）=1;S（S=0 & DD=2）=2;x,y=find（V=2）;DD（x+1,y）=P（x+1,y）;x,y=find（V=3）;DD（x,y-1）=P（x,y-1）;x,y=find（V=4）;DD（x,y+1）=P（x,y+1）;Dr=rand（size（S）;S（S3 & Drm）=0;file 12:4 zz 4 zz 2005-06-16 08:36100100401230 11 223 1 12c1c2 2

17、 i1 202 2 3 3 12m c1 122c2 0 3 2嶼 D0.330 3310.220112 0.110220.440 0file 13: 2006-02-07 13:32% ,0.1,0.15,0.56,0.7）;N=50;A=ones（N,N,3）;xn=round（N+1）/2）;yn=xn;A（xn,yn,:A（10:13,10:13,2:3）=0;A（40:43,10:13,40:43,2:43,40:H=imshow（A）;p=1,-1,i,-i;ss=while k=1;t=round（3*rand）+1;x=real（p（t）;y=imag（p（t）;xn=xn+x;yn=yn+y;if xnN;xn=mod（xn,N）+1;if ynyn=mod（xn,N）+1;if A（xn,yn,2）=0;xn=xn-x;yn=yn-y;else A（xn,yn,:set（H,A）;,0.74,0.6,0.2,0.08,.,ss）;