元胞自动机——森林火灾-2 - 数学建模社区-数学中国

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

复制代码