元胞自动机：森林火灾模拟（Matlab：heatmap、colormap）

最新推荐文章于 2023-07-28 16:11:51 发布
最新推荐文章于 2023-07-28 16:11:51 发布
阅读量2.5k 收藏 30
点赞数 6
分类专栏： 北师大 数学模型 文章标签： 算法 matlab
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_55818063/article/details/123757223
版权

目录

1.元胞自动机（Cellular Automation,CA）

2.森林火灾

3.伪代码

4.元胞自动机模拟森林火灾（Matlab实现）

（1）定义森林火灾函数

①空位生长树木 （0 --> 1）

②周围树木燃烧 （1-->-1）

③燃烧树木清除（-1-->0）

④雷电击中正常树木 （1-->-1）

（2）定义主函数

①遍历部分

②趋势曲线绘制部分

（3）执行命令

5.森林火灾动态仿真

1.元胞自动机（Cellular Automation,CA）

该部分详细内容请参考如下链接：

元胞自动机：森林火灾模拟（Python：numpy、seaborn）_Destiny坠明的博客-CSDN博客

2.森林火灾

该部分详细内容请参考如下链接：

元胞自动机：森林火灾模拟（Python：numpy、seaborn）_Destiny坠明的博客-CSDN博客

3.伪代码

该部分详细内容请参考如下链接：

元胞自动机：森林火灾模拟（Python：numpy、seaborn）_Destiny坠明的博客-CSDN博客

4.元胞自动机模拟森林火灾（Matlab实现）

（1）定义森林火灾函数

<Forest_Fire.m>

Forest_Fire函数如下，传入参数为：

（1）上一时刻的燃烧矩阵；

（2）空格位 树的生长概率p；

（3）正常位 树的燃烧概率f；

假定：空格位 --> 0 正常位 --> 1 燃烧位 --> -1

总量：空格位 --> C 正常位 --> G 燃烧位 --> R

输出参数为：规则演化后的燃烧矩阵；

该函数可分为4个操作部分：

function matrix=Forest_Fire(current_matrix,p,f)

    global area;
    
    matrix=current_matrix;
    
    % (1)空位生长树木 （0 --> 1）   储存位置
    i_C_indexes = [];
    j_C_indexes = [];

    for i = 1:1:area   % 行循环
        for j = 1:1:area   % 列循环
            if (matrix(i, j) == 0 & rand()<p)
                i_C_indexes(end+1)=i;
                j_C_indexes(end+1)=j;

            end
        end
    end

    
    % (2)周围树木燃烧   (1 --> -1)   储存位置，并存储上一时刻的燃烧数据
    fire_memory=find(matrix==-1);
    i_indexes=[];
    j_indexes=[];

    for i = 1:1:area   % 行循环
        for j = 1:1:area   % 列循环
            if (matrix(i,j)==-1)
                
                if((1<=i-1) & (i-1<=area) & (1<=j) & (j<= area))
                    if (matrix(i-1,j)==1)
                        i_indexes(end+1)=i-1;
                        j_indexes(end+1)=j;
                    end
                end
                
                if((1<=i+1) & (i+1<=area) & (1<=j) & (j<= area))
                    if (matrix(i+1,j)==1)
                        i_indexes(end+1)=i+1;
                        j_indexes(end+1)=j;
                    end
                end
                
                if((1<=i) & (i<=area) & (1<=j-1) & (j-1<= area))
                    if (matrix(i,j-1)==1)
                        i_indexes(end+1)=i;
                        j_indexes(end+1)=j-1;
                    end
                end

                if((1<=i) & (i<=area) & (1<=j+1) & (j+1<= area))
                    if (matrix(i,j+1)==1)
                        i_indexes(end+1)=i;
                        j_indexes(end+1)=j+1;
                    end
                end
            end
        end
    end
    
    for k = 1:1: length(i_indexes)
        matrix(i_indexes(k),j_indexes(k))=-1;
    end
    

    % (3)燃烧树木清除 （-1 --> 0）
    matrix(fire_memory)=0;

    
    % (4)雷电击中正常树木 （1 --> -1）    储存位置
    i_indexes = [];
    j_indexes = [];
    
    for i = 1:1:area   % 行循环
        for j = 1:1:area   % 列循环
            if matrix(i, j) == 1
                if(1<=i-1 & i-1<=area & 1<=j & j<= area)
                    if(matrix(i-1, j) == -1)
                        continue
                    end
                end
                
                if(1<=i+1 & i+1<=area & 1<=j & j<= area)
                    if(matrix(i+1, j) == -1)
                        continue
                    end
                end
                
                if(1<=i & i<=area & 1<=j-1 & j-1<= area)
                    if(matrix(i, j-1) == -1)
                        continue
                    end
                end
                
                if(1<=i & i<=area & 1<=j+1 & j+1<= area)
                    if(matrix(i, j+1) == -1)
                        continue
                    end
                end
            
                if rand()<f
                    i_indexes(end+1)=i;
                    j_indexes(end+1)=j;
                end
                
            end
        end
    end
    
    for k = 1:1: length(i_indexes)
        matrix(i_indexes(k),j_indexes(k))=-1;
    end

    % 完成空位生长树木的操作
    for k = 1:1: length(i_C_indexes)
        matrix(i_C_indexes(k),j_C_indexes(k))=1;
    end

 

area为全局变量。（森林的边长）

function matrix=Forest_Fire(current_matrix,p,f)

    global area;
    
    matrix=current_matrix;

全局变量声明规则：

如在主函数中，你要设置n为全局变量，则需声明 global n;
然后在子函数里面你用到了n这个全局变量，则需在子函数中再次声明 global n;
这样在子函数中，就可以使用n这个全局变量了。

 

①空位生长树木 （0 --> 1）

    % (1)空位生长树木 （0 --> 1）   储存位置
    i_C_indexes = [];
    j_C_indexes = [];

    for i = 1:1:area   % 行循环
        for j = 1:1:area   % 列循环
            if (matrix(i, j) == 0 & rand()<p)
                i_C_indexes(end+1)=i;
                j_C_indexes(end+1)=j;

            end
        end
    end

        每当一个新的演化时间点到来时，首先要执行 空位生长树木 （0 --> 1） 的操作。若这一操作滞后，则会导致前一时刻燃烧的树木在这一时刻烧尽后重新生长为树木，故这一演化的优先级较高。并且需要将符合条件的坐标存储（而不是直接将矩阵元素更改），若不如此处理，则会导致中心灰烬点重焕生机后再次被周围的余火吞噬，无法达到 空位生长树木 （0 --> 1） 的效果。

        综上，该演化优先级最高并且需要储存符合要求的坐标位置。
 

②周围树木燃烧 （1-->-1）

    % (2)周围树木燃烧   (1 --> -1)   储存位置，并存储上一时刻的燃烧数据
    fire_memory=find(matrix==-1);
    i_indexes=[];
    j_indexes=[];

    for i = 1:1:area   % 行循环
        for j = 1:1:area   % 列循环
            if (matrix(i,j)==-1)
                
                if((1<=i-1) & (i-1<=area) & (1<=j) & (j<= area))
                    if (matrix(i-1,j)==1)
                        i_indexes(end+1)=i-1;
                        j_indexes(end+1)=j;
                    end
                end
                
                if((1<=i+1) & (i+1<=area) & (1<=j) & (j<= area))
                    if (matrix(i+1,j)==1)
                        i_indexes(end+1)=i+1;
                        j_indexes(end+1)=j;
                    end
                end
                
                if((1<=i) & (i<=area) & (1<=j-1) & (j-1<= area))
                    if (matrix(i,j-1)==1)
                        i_indexes(end+1)=i;
                        j_indexes(end+1)=j-1;
                    end
                end

                if((1<=i) & (i<=area) & (1<=j+1) & (j+1<= area))
                    if (matrix(i,j+1)==1)
                        i_indexes(end+1)=i;
                        j_indexes(end+1)=j+1;
                    end
                end
            end
        end
    end
    
    for k = 1:1: length(i_indexes)
        matrix(i_indexes(k),j_indexes(k))=-1;
    end

        在这一部分需要注意边界特殊情况导致的矩阵索引越界报错。

        存储上一时刻的燃烧数据 fire_memory 的目的是为了在之后的树木烧尽操作中区分出上一时刻燃烧的树木，避免与当前时刻新燃烧的树木产生混淆。

        储存符合条件坐标的原因是避免当前时刻新燃烧的树木引燃树木（这是一个很严重的bug），所以需要先存储所有坐标信息，再统一赋新值处理。
 

注意：Matlab在判断一个标量是否在一个区间内，要将两个边界的条件拆开写！

③燃烧树木清除（-1-->0）

    % (3)燃烧树木清除 （-1 --> 0）
    matrix(fire_memory)=0;

        直接将上一操作存储的燃烧数据坐标 fire_memory 对应位置的元素归零，实现燃烧树木清除 （-1 --> 0）效果。

注：fire_memory存储的是一个向量（而非矩阵），

每一个数值标量对应着一个索引（Matlab的存储结构特性）。

④雷电击中正常树木 （1-->-1）

    % (4)雷电击中正常树木 （1 --> -1）    储存位置
    i_indexes = [];
    j_indexes = [];
    
    for i = 1:1:area   % 行循环
        for j = 1:1:area   % 列循环
            if matrix(i, j) == 1
                if(1<=i-1 & i-1<=area & 1<=j & j<= area)
                    if(matrix(i-1, j) == -1)
                        continue
                    end
                end
                
                if(1<=i+1 & i+1<=area & 1<=j & j<= area)
                    if(matrix(i+1, j) == -1)
                        continue
                    end
                end
                
                if(1<=i & i<=area & 1<=j-1 & j-1<= area)
                    if(matrix(i, j-1) == -1)
                        continue
                    end
                end
                
                if(1<=i & i<=area & 1<=j+1 & j+1<= area)
                    if(matrix(i, j+1) == -1)
                        continue
                    end
                end
            
                if rand()<f
                    i_indexes(end+1)=i;
                    j_indexes(end+1)=j;
                end
                
            end
        end
    end
    
    for k = 1:1: length(i_indexes)
        matrix(i_indexes(k),j_indexes(k))=-1;
    end

        这一部分的算法实现过程与操作2类似，不再赘述。

    % 完成空位生长树木的操作
    for k = 1:1: length(i_C_indexes)
        matrix(i_C_indexes(k),j_C_indexes(k))=1;
    end

        函数的最后完成第一步 空位生长树木 （0 --> 1） 的操作，最终返回一个矩阵（二维数组）。

（2）定义主函数

<main.m>

main函数如下，传入参数为：

（1）森林边长area；

（2）演替次数N；

（3）空格位 树的生长概率p；

（4）正常位 树的燃烧概率f；

输出参数为：最终的燃烧矩阵；

该函数可分为2个操作部分：

function matrix=main(area,N,p,f)
    C_count=[];
    G_count=[];
    R_count=[];
    
    map = [239/255,29/255,31/255;
            165/255,165/255,165/255;
            28/255,172/255,76/255];
    
    matrix=ones(area,area);
    
    for time=1:1:N
        matrix=Forest_Fire(matrix,p,f);
        
        C_count(end+1)=length(find(matrix==0));
        G_count(end+1)=length(find(matrix==1));
        R_count(end+1)=length(find(matrix==-1));
    
   
        hfig=figure(1);
        set(hfig, 'position', get(0,'ScreenSize'));
        h=heatmap(matrix);
        colormap(map);
        %h.CellLabelFormat = '%d';
        %h.Title = 'Forest Fire';
        %h.XLabel = 'X';
        %h.YLabel = 'Y';
        h.CellLabelColor = 'none';

    end
    
        
    T=1:1:N;
    hfig=figure(2);
    set(hfig, 'position', get(0,'ScreenSize'));
    hold on;
    plot(T,C_count,'LineWidth',2,'Color',[165/255,165/255,165/255]);
    plot(T,G_count,'LineWidth',2,'Color',[28/255,172/255,76/255]);
    plot(T,R_count,'LineWidth',2,'Color',[239/255,29/255,31/255]);
    
    grid on;
    
    legend('Vacancy','Trees','Burning');
    
    hold off;

 

热力图自定义颜色：

map= [239/255,29/255,31/255;
           165/255,165/255,165/255;
           28/255,172/255,76/255];

colormap(map);

function matrix=main(area,N,p,f)
    C_count=[];
    G_count=[];
    R_count=[];
    
    map = [239/255,29/255,31/255;
            165/255,165/255,165/255;
            28/255,172/255,76/255];

    matrix=ones(area,area);

①遍历部分

    for time=1:1:N
        matrix=Forest_Fire(matrix,p,f);
        
        C_count(end+1)=length(find(matrix==0));
        G_count(end+1)=length(find(matrix==1));
        R_count(end+1)=length(find(matrix==-1));
    
        hfig=figure(1);
        set(hfig, 'position', get(0,'ScreenSize'));
        h=heatmap(matrix);
        colormap(map);
        %h.CellLabelFormat = '%d';
        %h.Title = 'Forest Fire';
        %h.XLabel = 'X';
        %h.YLabel = 'Y';
        h.CellLabelColor = 'none';

    end

调整figure尺寸（铺满屏幕）：

hfig=figure(1);
set(hfig, 'position', get(0,'ScreenSize'));

②趋势曲线绘制部分

    T=1:1:N;
    hfig=figure(2);
    set(hfig, 'position', get(0,'ScreenSize'));
    hold on;
    plot(T,C_count,'LineWidth',2,'Color',[165/255,165/255,165/255]);
    plot(T,G_count,'LineWidth',2,'Color',[28/255,172/255,76/255]);
    plot(T,R_count,'LineWidth',2,'Color',[239/255,29/255,31/255]);
    
    grid on;
    
    legend('Vacancy','Trees','Burning');
    
    hold off;

        这一段代码负责绘制出树木三种状态的数量变化趋势图。

 

（3）执行命令

<Cellular_Automation>

global area;
area=100;

final_matrix=main(area,100,0.25,0.01);


% figure(2);
% 
% h=heatmap(final_matrix);
% 
% h.CellLabelFormat = '%d';
% h.Title = 'Forest Fire';
% h.XLabel = 'X';
% h.YLabel = 'Y';
% 
% %h.CellLabelColor = 'none';
% 
% map = [239/255,29/255,31/255;
%     165/255,165/255,165/255;
%     28/255,172/255,76/255];
% 
% colormap(map)

 

5.森林火灾动态仿真

Matlab森林火灾仿真

Destiny坠明
关注
  • 6
    点赞
  • 踩
  • 30
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
专栏目录
【数学建模】基于matlab GUI森林求火问题【含Matlab源码 4001期】
订阅付费专栏（初级版），可赠送该付费专栏任意代码1份；
02-23 74
森林求火问题 完整的代码，方可运行；可提供运行操作视频！适合小白！
matlab元胞自动机森林火灾高仿真代码
08-05
对元胞自动机森林火灾模型进行了深度优化修复了较多bug，权衡了树的生长和火灾，达到一个较好的仿真效果，适合初学！
(完整版)元胞自动机-森林火灾模型MATLAB代码.pdf
07-13
(完整版)元胞自动机-森林火灾模型MATLAB代码.pdf(完整版)元胞自动机-森林火灾模型MATLAB代码.pdf(完整版)元胞自动机-森林火灾模型MATLAB代码.pdf(完整版)元胞自动机-森林火灾模型MATLAB代码.pdf(完整版)元胞自动机-森林火灾模型MATLAB代码.pdf(完整版)元胞自动机-森林火灾模型MATLAB代码.pdf(完整版)元胞自动机-森林火灾模型MATLAB代码.pdf(完整版)元胞自动机-森林火灾模型MATLAB代码.pdf
【元胞自动机】基于元胞自动机模拟森林救火问题附matlab代码
最新发布
m0_57702748的博客
07-28 335
卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断。图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知。
森林火灾模型代码（matlab）
pp2373886592的博客
08-03 4448
%森林着火模型 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(:,
元胞自动机森林火灾_hurtn3k_森林火灾_火灾_森林火灾程序_元胞自动机模拟森林火灾_
09-30
在本主题中，我们聚焦于一个特定的应用：使用元胞自动机模拟森林火灾。森林火灾是一个典型的复杂系统，涉及多种变量，如风向、湿度、温度、植被类型等。通过元胞自动机模拟，我们可以理解和预测火灾可能的行为，这...
元胞自动机.zip_fire_matlab_ship8bl_元胞自动机_森林火灾程序
07-14
采用MATLAB语言编写的森林火灾元胞自动机程序，使用过程中，可以根据具体的情况，进行相应参数的修改。
基于元胞自动机模拟森林火灾附MATLAB代码
06-29
在这个特定的案例中，我们关注的是使用元胞自动机来模拟森林火灾的过程。元胞自动机通过定义简单的规则在网格上的每个单元格（元胞）进行迭代更新，从而能够展现复杂的集体行为，如火势蔓延。 MATLAB是一种强大的...
基于元胞自动机模拟晶体生成附matlab代码
07-05
本资源包含基于MATLAB 2019a的元胞自动机模拟晶体生成的代码，对于本科及硕士阶段的学生来说，是进行教研学习的理想材料。MATLAB是一种强大的数学计算软件，特别适合于数值分析、矩阵运算、图形绘制等任务，同时也...
基于元胞自动机森林火灾模型及应用（Matlab代码实现）
weixin_46039719的博客
12-10 836
元胞自动机是一种具有时空设计能力的动力学模型,按照一定的局部规则,在离散的时间维上演化的动力学系统”﹐其特点是:散布于规格网格每一元胞均取有限的离散状态,各元胞遵循相同的演化规则进行同步演化且仅和它邻近的元胞发生相互作用。元胞自动机最初由冯﹒诺依曼于20世纪50年代提出,用以模拟生命系统所具有的自复制功能。随着人工生命和复杂科学的发展,元胞自动机的应用领域迅速扩展,已经成功运用于许多自然现象的模拟,如，Is ing模型,交通流的仿真模拟﹐量子元胞自动机[.等。部分理论来源于网络，如有侵权请联系删除。
元胞自动机：森林火灾模拟（Python：numpy、seaborn）
Destiny坠明的博客
03-25 8730
目录 1.元胞自动机（Cellular Automation,CA） 2.森林火灾 3.伪代码 4.元胞自动机模拟森林火灾（Python实现） （1）导入相关的库 （2）定义森林火灾函数（Forest_Fire函数） ①空位生长树木 （0 --> 1） ②周围树木燃烧 （1 --> -1） ③烧树木清除 （-1 --> 0） ④雷电击中正常树木 （1 --> -1） （3）定义主函数（main函数） （4）执行命令 5.森林火灾可视化呈现 （1）森林火灾各
森林救火模型matlab,元胞自动机模拟森林火灾--matlab实现
weixin_35122299的博客
03-27 1516
因为建模需要，学习了元胞自动机，并尝试了用其模拟森林火灾，相对于退火和遗传，元胞自动机更容易理解，并且其广泛应用于交通问题，在以往的建模比赛中占了很大的比重。一.元胞自动机定义和适用范围不同于一般的动力学模型，元胞自动机不是由严格定义的物理方程或函数确定，而是用一系列模型构造的规则构成。凡是满足这些规则的模型都可以算作是元胞自动机模型。因此，元胞自动机是一类模型的总称，或者说是一个方法框架。其特点...
森林救火模型matlab,数学建模1-森林救火模型
weixin_33093323的博客
03-16 2299
本篇文章是我简书上的第一篇创作，由于今年大三一年即将度过，感觉时光匆忙，以前的有许多的不努力和懊悔，可是现在没有办法进行补偿，所以想通过这几篇文章帮助正在上学的小伙伴。希望看到这篇文章的人儿，停下你的脚步，珍惜时间，砥砺前行。与此同时，我也希望自己将这件有意义的分享知识的事情坚持下去。同时也希望大家能够不吝赐教，帮助我改进文章质量。什么是数学建模，如果用正规一点的文本来讲的话，应该这样叙述：数学建...
数学建模森林着火问题Matlab,数学建模森林救火问题.doc
weixin_32894769的博客
03-21 1997
数学建模森林救火问题.doc森林救火问题的研究【摘要】：森林救火问题是一个优化问题，经过分析我们决定采用极值法和定积分的方法来求森林烧毁的面积，从而解决该问题，通过对问题的剖析，得出表达式：救火的总费用=单位森林面积损失费×损失面积每个队员的单位时间灭火费用人数灭火时间+单位人数一次性支出×参加救火的消防员人数.对各个量进行分析，得知森林损失面积较为难求，于是我们将其单独考虑。在有风的情况下，火势...
元胞自动机-森林火灾模拟
席木木的博客
03-25 6400
引入：     元胞自动机，英文名及缩写：cellular automata，CA。最初是由冯诺依曼在二十世纪五十年代为模拟生物自保的自我复制而提出的，但是当时并未受到重视。后来才逐渐发展起来，著名的“生命游戏”就是其一个应用案例。但今天想和大家分享的案例并不是生命游戏，而是来自近期课程实验中的内容：根据元胞自动机原理编程实现“森林火灾模拟”。 正文: ...
基于元胞自动机模拟森林火灾附matlab代码
m0_60703264的博客
06-17 54
1 算法原理元胞自动机（Cellular Automaton）是一种离散化的数学模型，由一系列相互作用的自动单元（称为元胞）组成。每个元胞处于离散状态，并通过一组规则和相邻元胞的交互进行更新。元胞自动机的原理包括以下几个关键要素：元胞：元胞是模型中最小的单位，通常表示为一个离散位置。每个元胞可以具有多种离散状态。状态：每个元胞具有一定的状态，例如二进制（0或1），有限的离散值等。状态可以表示元胞的属性或特征。规则：规则定义了元胞如何根据其自身状态和相邻元胞的状态来更新。
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值