元胞自动机——渗流集群

森之张卫东

1. <p>%渗流集群
   
2. %Percolation Cluster
   
3. clc
   
4. clear
   
5. close all
   
6. threshold=0.63;
   
7. %最初的图形是由图形操作决定的，以下程序设定坐
   
8. %标系为一个固定的尺寸，在坐标系里写入文本，然后
   
9. %获得并返回坐标系内容，并用getframe函数把它们
   
10. %写入一个矩阵。
    
11. ax=axes( 'units','pixels','position',[1,1,1050,400],'color','k' );%1050,400画布大小！
    
12. text( 'units','pixels', 'position',[50 270 0],...
    
13.     'string','ZhangWeiDong','color','w','fontname','helvetica',...
    
14.     'fontsize',50 )
    
15. text( 'units','pixels', 'position',[120 150 0],...
    
16.     'string','1302034220','color','w','fontname','helvetica',...
    
17.     'fontsize',50 )
    
18. initial=getframe( gca );</p><p>[a,b,c]=size(initial.cdata);
    
19. z=zeros(a,b);
    
20. cells=double( initial.cdata(:,:,1)==255 );
    
21. visit=z;
    
22. sum=z;
    
23. imh=image( cat(3,z,cells,z) );
    
24. set( imh,'erasemode','none' )</p><p>for i=1:100
    
25.     sum(2:a-1,2:b-1)=cells(2:a-1,1:b-2)+cells(2:a-1,3:b)+...
    
26.                              cells(1:a-2,2:b-1)+cells(3:a,2:b-1)+...
    
27.                              cells(1:a-2,1:b-2)+cells(1:a-2,3:b)+...
    
28.                              cells(3:a,1:b-2)+cells(3:a,3:b);
    
29.     pick=rand(a,b);
    
30.     %edge only
    
31.     cells=cells | ( (sum>=1)&(pick>=threshold)&(visit==0) );
    
32.     visit=(sum>=1);
    
33.    
    
34.     set( imh,'cdata',cat(3,z,cells,z) )
    
35.     drawnow
    
36.     %pause(0.5)
    
37. end</p><p>return
    
38. figure(2)
    
39. image( cat(3,z,cells,z) )</p><p><p><p><p><p><p>
    
40. </p>

复制代码