元胞自动机——森林火灾-2

森之张卫东

元胞自动机——森林火灾-2

1. <p>
   
2. %元胞自动机——森林火灾
   
3. %时间：2015-8-29-六
   
4. %----------------------------------------------------
   
5. n=200;   %表示森林矩阵的尺寸：n x n
   
6. Pltg=5e-6;   %闪电的概率
   
7. Pgrw=1e-2;   %生长的概率
   
8. NW=[n,1:n-1];   %用于构造北邻居 veg(NW,:)和西邻居 veg(:,NW)
   
9. SE=[2:n,1];        %用于构造南邻居 veg(SE ,:)和西邻居 veg(:, SE)
   
10. veg=zeros(n);   %veg={0表示空，1表示火，2表示树}</p><p>imh=image ( cat( 3,(veg==1),(veg==2),zeros(n) ) );</p><p>for i=1:3000
    
11.     %sum为周围四个邻居中状态为火的数量
    
12.     num=             ( veg(NW,:)==1 ) +...
    
13.         ( veg(:,NW)==1 )    +    ( veg(:,SE)==1 )+...
    
14.                           ( veg(SE,:)==1 ) ;
    
15.     %规则体现
    
16.     veg=2*( (veg==2) | (veg==0 & rand(n)<Pgrw) )-...
    
17.         ( (veg==2) & (num>0 | rand(n)<Pltg) );
    
18.    
    
19.     set( imh,'cdata', cat( 3,(veg==1),(veg==2),zeros(n) ) );  % cdata 属性在:Graphics---images---image file operations中
    
20.     %下面这个为实验语句，简单说明：image---第一页数据为红色；第二页数据为绿色；第三页数据为蓝色(default情况下)
    
21.     %set( imh,'cdata', cat( 3,(veg==2),(veg==1),zeros(n) ) );
    
22.    
    
23.     drawnow
    
24.    
    
25.     pause(0.1)
    
26.    
    
27. end</p>

复制代码

小琳琳xiao

啦啦啦~~~我是默认签名(*^__^*)

小琳琳xiao

啦啦啦~~~我是默认签名(*^__^*)