森林火灾模型代码（matlab）

%森林着火模型
n = 300; %总规模
Plight = 0.000005; Pgrowth = 0.01; %生长和雷电概率
UL = [n 1:n-1]; DR = [2:n 1];
veg = zeros(n, n); %初始化
img = image(cat(3, veg, veg, veg));
%循环模拟开始
for i = 1 : 1 : 3000
   sum = (veg(UL, :)==1) + (veg(DR, :)==1) + (veg(:, UL)==1) + (veg(:, DR)==1);
   veg = 2 * (veg == 2) -...
       ((veg == 2) & (sum > 0 | rand(n, n) < Plight)) +...
       2 * ((veg == 0) & rand(n, n) < Pgrowth); %更新树林
   
   set(img, 'cdata', cat(3, (veg == 1), (veg == 2), zeros(n)));
   drawnow 
end