【元胞自动机】基于元胞自动机模拟森林大火matlab源码

1 森林火灾元胞自动机原理

在元胞自动机模型中，空间被离散成网格，每一个网格被称为元胞。森林火灾元胞有三种状态：树，火（正在燃烧的树）和空（空地）状态。元胞下一时刻状态的更新规则如下： 树变火：一棵树，其上下左右若有一个状态为火，下一刻就会变成火。或者一棵树遇上闪电，下一刻就会变成火。由于遇上闪电着火的概率Plight很小。

火变空：火在下一时刻会变成空。

空变树：空地下一时刻会以一个很小的概率Pgrowth长出新树。

改进模型会考虑树的对角位置有没有着火。或者会考虑风向（比如吹西风（火从东吹向西），火的西边着火的机率会变大（顺风），火的东边着火的几率变小（逆风）），这里盗个图：

图a是基础元胞自动机，图b是考虑对角的元胞自动机，图c是吹西风的元胞自动机。 本人这里对基础元胞自动机，考虑对角情况的元胞自动机，考虑吹西风的元胞自动机三种模型进行仿真.

2 部分代码

P=[];for m=1:10clear D T fire_time_lightning fire_time_itself aspect tdata Index;%% O