数学建模模型11——元胞自动机

模拟森林火灾

代码

clc;
clear;
close;

%0为空地 1为燃烧的树 2为正常的树
%初始情况是森林数目有一定的覆盖率 某些树木着火
length=300; % 森林长
width=500; % 森林宽
coverage_rate=0.8;% 初始森林覆盖率
growth_rate=0.0001;% 新树生长的概率
spread_rate=0.7;% 火情扩散率


tes=rand(length,width);
tt=(tes<growth_rate);
u=sum(tt(:));


SUM=zeros(length,width);%计算每个点有多少个邻居着火
a=randperm(length*width);
forest=reshape(2*(a<=coverage_rate*length*width),length,width);

% initial_fire_tree_number=10;
% 随机生成着火的树的位置【不太准确 应该使得着火的树应该比较集中？】
% b=