clc,clear
% 森林火灾
% 演化规则
% 1. 如树4个邻居中有火，树下一时刻变火
% 2. 火下一时刻变空
% 3. 树有低概率变火;(闪电燃烧)
% 4. 空有低概率变树;(新树生长)
n=200; % 森林边长
pltg=0.0005; %闪电概率
pgrw=0.01; %生长概率
NW=[n,1:n-1];  %第一行的北面是第n行，第一列的西面是第n列 
SE=[2:n,1];    %第一行南面是第2行，第一列东面是第2列
veg=zeros(n); %0空，1火，2树
imh=image(cat(3,(veg==1),(veg==2),zeros(n)));  %cat将数组组合起来
for i=1:3000
    % num周围4个邻居为火状态的总量
    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)))  % 三维分别代表R\G\B，火红色，树绿色
    drawnow
end

效果如下：