MATLAB元胞自动机算法学习之森林火灾代码演示

转自：MATLAB演示元胞自动机算法
作者：Lindark
部分改动，侵删。

一、元胞自动机理论

元胞自动机与格子理论是一个非常好的模型，许多复杂的问题都可以通过它来建立模型，下面就简要介绍一下。

元胞自动机
实质上是定义在一个具有离散、有限状态的元胞组成的元胞空间上，并按照一定的局部规则，在离散的时间维度上演化的动力学系统。

元胞
元胞又可称为单元、细胞，是元胞自动机的最基本的组成部分。

元胞具有以下特点：

元胞自动机最基本的单元。
元胞有记忆贮存状态的功能。
所有元胞状态都按照元胞规则不断更新。

演化规则
中心元胞的下一个状态由中心元胞的当前状态和其邻居的当前状态按照一定的规则确定。

二、森林火灾的演示

下面就用MATLAB来演示森林火灾，以便更好地理解元胞自动机理论。

森林火灾的元胞自动机模型有三种状态：空位，燃烧着的树木及树木。则某元胞下一时刻状态由该时刻本身的状态和周围四个邻居的状态以一定的规则确定，规则如下：

如果某树木元胞的4个邻居有燃烧着的，那么该元胞下一时刻的状态是燃烧着的。
一个燃烧着的元胞在下一时刻变成空位。
所有树木元胞以一个低概率开始燃烧（模拟闪电引起的火灾）
所有空元胞以一个低概率变成树木（以模拟新的树木的生长）

clc
n = 300;     %元胞矩阵大小
Plight = .000005; Pgrowth = .01;   %Plight:模拟闪电引起树木燃烧的概率，Pgrowth：模拟树木生长的概率
UL = [n 1:n-1];DR = [2:n 1];    %为下面计算sum做铺垫
veg = zeros(n,n);        %初始化
% The value of veg:
% empty == 0
% burning == 1
% green == 2
imh = image(cat(3,veg,veg,veg));    %将矩阵用图像显示
for i = 1:300%计算每个元胞4个邻居的值(逻辑值)的和sum = (veg(UL,:) == 1) + (veg(:,UL) == 1) + (veg(DR,:) == 1) + (veg(:,DR) == 1);%根据规则更新森林矩阵：树 = 树 - 着火的树 + 新生的树veg =  2*(veg == 2) - ( (veg == 2) & (sum > 0 | (rand(n,n) < Plight)) ) + 2*( (veg == 0) & rand(n,n) < Pgrowth);%设置imh图像数据的颜色，veg=1为红色，veg=2为绿色,veg=0默认为黑色set(imh, 'cdata', cat(3, (veg == 1), (veg == 2), zeros(n)) )drawnowpause(0.01)
end

运行的结果如下，可以通过调整闪电和树木生长的概率来观察不同的结果。

MATLAB元胞自动机算法学习之森林火灾代码演示相关推荐

雨水冲刷沙堡模型——matlab元胞自动机
雨水冲刷沙堡模型,新手练习matlab元胞自动机写的. 这是代码 clear all; % 设置GUI按键 plotbutton=uicontrol('style','pushbutton','str ...
【机器学习】通俗的元胞自动机算法解析和应用
[机器学习]通俗的元胞自动机算法解析和应用文章目录 1 元胞自动机的定义 2 元胞自动机的组成 3 元胞自动机的特征 4 Python实现元胞自动机(生命游戏) 5 总结 6 Github(华盛顿州 ...
【Matlab元胞自动机】元胞自动机大型商场人流疏散【含源码 665期】
一.代码运行视频(哔哩哔哩) [Matlab元胞自动机]元胞自动机大型商场人流疏散[含源码 665期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1]贺琳,聂化 ...
MATLAB演示元胞自动机算法
一.元胞自动机理论元胞自动机与格子理论是一个非常好的模型,许多复杂的问题都可以通过它来建立模型,下面就简要介绍一下. 元胞自动机实质上是定义在一个具有离散.有限状态的元胞组成的元胞空间上,并按照一 ...
【元胞自动机】基于matlab元胞自动机3D森林火灾模型【含Matlab源码 656期】
⛄一.元胞自动机简介 1 元胞自动机发展历程最初的元胞自动机是由冯 · 诺依曼在 1950 年代为模拟生物细胞的自我复制而提出的. 但是并未受到学术界重视. 1970 年, 剑桥大学的约翰 · 何 ...
【元胞自动机】基于matlab元胞自动机多车道信号交叉口仿真【含Matlab源码 818期】
⛄一.元胞自动机简介 1 元胞自动机发展历程最初的元胞自动机是由冯 · 诺依曼在 1950 年代为模拟生物细胞的自我复制而提出的. 但是并未受到学术界重视. 1970 年, 剑桥大学的约翰 · 何 ...
matlab 元胞自动机
用MATLAB实现元胞自动机的事例 file 1:life.m %% 初始化 m = 50; X = zeros(m,m); X(25,25) = 1;n = [m 1:m-1];e = [2:m 1 ...
【元胞自动机】基于matlab元胞自动机模拟交通路况（含超车）【含Matlab源码 2389期】
⛄一.元胞自动机简介 1 元胞自动机发展历程最初的元胞自动机是由冯 · 诺依曼在 1950 年代为模拟生物细胞的自我复制而提出的. 但是并未受到学术界重视. 1970 年, 剑桥大学的约翰 · 何 ...
【元胞自动机】基于matlab元胞自动机模拟SEIR传播模型和采取隔离措施的SEIR模型【含Matlab源码 2181期】
⛄一.元胞自动机简介 1 元胞自动机发展历程最初的元胞自动机是由冯 · 诺依曼在 1950 年代为模拟生物细胞的自我复制而提出的. 但是并未受到学术界重视.1970 年, 剑桥大学的约翰 · 何顿 ...
【元胞自动机】基于matlab元胞自动机双车道交通流模型含靠右行驶【含Matlab源码 231期】
⛄一.元胞自动机简介 1 元胞自动机发展历程最初的元胞自动机是由冯 · 诺依曼在 1950 年代为模拟生物细胞的自我复制而提出的. 但是并未受到学术界重视. 1970 年, 剑桥大学的约翰 · 何 ...

MATLAB元胞自动机算法学习之森林火灾代码演示

一、元胞自动机理论

二、森林火灾的演示

MATLAB元胞自动机算法学习之森林火灾代码演示相关推荐

最新文章

热门文章