【元胞自动机】基于元胞自动机实现传染病传播模拟附matlab代码

简介: 【元胞自动机】基于元胞自动机实现传染病传播模拟附matlab代码

✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信。

🍎个人主页:Matlab科研工作室

🍊个人信条:格物致知。

更多Matlab仿真内容点击👇

智能优化算法       神经网络预测       雷达通信      无线传感器        电力系统

信号处理              图像处理               路径规划       元胞自动机        无人机

⛄ 内容介绍

元胞自动机,亦被称为细胞自动机。CA的经典案例是定义一个网格,网格上的每个点代表一个有限数量的状态中的细胞。过渡规则同时应用到每一个细胞。典型的转换规则依赖于细胞和它的(4个或8个)近邻的状态,虽然临近的细胞也同样使用。CA的应用在并行计算研究、物理模拟和生物模拟等领域。在数学建模中,一般是借鉴元胞自动机的概念,应用于具体的适合于机理建模的问题中。这类问题的典型特征是,所研究的问题是一个系统问题,系统是由若干个一个或几个不同类的对象组成,经典的模型不适应。典型的问题如滴滴打车问题、开发小区问题。

 

MATLAB机理建模方法在数学建模中,如果遇到一个非典型的数学建模问题(非数据、优化、连续、评价),那么这种情况下,通常需要用到机理建模方法了。

机理建模概述机理建模就是根据对现实对象特性的认识,分析其因果关系,找出反映内部机理的规则,然后建立规则的数学模型。机理建模的经典案例有很多,比如万有引力公式的推导过程。机理建模常见的有两类,一类是推导法机理建模,类似于微分方程建模,常用于动力学的建模过程,比如化学中反应动力学,还有各种场的方程,比如压力场、热场方程等;一类是包含一个或几个类别对象的复杂系统问题,常通过元胞自动机一仿真法来进行机理建模。下面将介绍这两类机理建模的具体MATLB实现过程。

·需求分析

1 元胞自动机发展历程

最初的元胞自动机是由冯· 诺依曼在 1950 年代为模拟生物 细胞的自我复制而提出的. 但是并未受到学术界重视.

1970 年, 剑桥大学的约翰 · 何顿 · 康威设计了一个电脑游戏 “生命游戏” 后, 元胞自动机才吸引了科学家们的注意.

1983 年 S.Wolfram 发表了一系列论文. 对初等元胞机 256 种 规则所产生的模型进行了深入研究, 并用熵来描述其演化行 为, 将细胞自动机分为平稳型, 周期型, 混沌型和复杂型.

2 对元胞自动机的初步认识

元胞自动机(CA)是一种用来仿真局部规则和局部联系的方法。典型的元胞自动机是定义在网格上的,每一个点上的网格代表一个元胞与一种有限的状态。变化规则适用于每一个元胞并且同时进行。典型的变化规则,决定于元胞的状态,以及其( 4 或 8 )邻居的状态。

3 元胞的变化规则&元胞状态

典型的变化规则,决定于元胞的状态,以及其(4 或 8 )邻居的状态。

4 元胞自动机的应用

元胞自动机已被应用于物理模拟,生物模拟等领域。

5 元胞自动机的matlab编程

结合以上,我们可以理解元胞自动机仿真需要理解三点。一是元胞,在matlab中可以理解为矩阵中的一点或多点组成的方形块,一般我们用矩阵中的一点代表一个元胞。二是变化规则,元胞的变化规则决定元胞下一刻的状态。三是元胞的状态,元胞的状态是自定义的,通常是对立的状态,比如生物的存活状态或死亡状态,红灯或绿灯,该点有障碍物或者没有障碍物等等。

⛄ 部分代码

close;clear;clc%参数初始化n = 250;chushiqianfuzhe = 5;% 邻居采用 8邻居式% 0 :代表易感者% 1 :代表潜伏者% 2 :代表患病者% 3 :代表免疫者% 4 :代表死亡者% 5 :  代表虚拟值    end  figure(1);  hold on;  set(c,'cdata',cat(3,R,G,B));  %image(cat(3,R,G,B));  drawnow;  figure(2);  plot(yiganzhe);  hold on;  plot(qianfuzhe);  hold on;  plot(huanbingzhe);  hold on;  plot(mianyizhe);  hold on  plot(siwangzhe);  title(['时间: ',num2str(i)]);  legend(['无病者数量',num2str(yiganzhe(end))],['潜伏者数量',num2str(qianfuzhe(end))],['患病者数量',num2str(huanbingzhe(end))],['免疫者数量',num2str(mianyizhe(end))],['死亡者数量',num2str(siwangzhe(end))]);  %顺序和plot顺序一致  pause(0.0001);end

⛄ 运行结果

这次的主要内容是现在最最最热门的传染病问题,众所周知今年的新冠病毒来势凶猛,无论是在国内还是国外都造成了巨大的影响。用模型来预测或者说模拟人群的感染程度也是不久以前非常热门的工作。一般会先用推导法建模,找出事物之间本质的定量联系,然后再看看是否适合数值仿真,如果适合仿真,此时就可以考虑使用元胞自动机一仿真法了。这两种建模方法往往相辅相成,推导法为仿真提供理论基础,仿真法为推导法提供验证和改进依据,两种方法相得益彰,不断促进模型的提升。

⛄ 参考文献

[1] 李光亮,温利华,闫俊霞,等.基于元胞自动机的传染病传播模拟研究[J].长江大学学报(自科版), 2013(9):85-87.DOI:10.3969/j.issn.1673-1409.2013.09.029.

[2] 李眩,童百利,吴晓兵.基于元胞自动机的传染病传播和控制仿真模型[J].辽宁工业大学学报:自然科学版, 2020, 40(5):6.DOI:10.15916/j.issn1674-3261.2020.05.003.

⛳️ 代码获取关注我

❤️部分理论引用网络文献,若有侵权联系博主删除
❤️ 关注我领取海量matlab电子书和数学建模资料

🍅 仿真咨询

1.卷积神经网络(CNN)、LSTM、支持向量机(SVM)、最小二乘支持向量机(LSSVM)、极限学习机(ELM)、核极限学习机(KELM)、BP、RBF、宽度学习、DBN、RF、RBF、DELM实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断
2.图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知
3.旅行商问题(TSP)、车辆路径问题(VRP、MVRP、CVRP、VRPTW等)、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划
4.无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配
5.传感器部署优化、通信协议优化、路由优化、目标定位
6.信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号
7.生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化
8.微电网优化、无功优化、配电网重构、储能配置
9.元胞自动机交通流 人群疏散 病毒扩散 晶体生长


相关文章
|
4月前
|
安全
【2023高教社杯】D题 圈养湖羊的空间利用率 问题分析、数学模型及MATLAB代码
本文介绍了2023年高教社杯数学建模竞赛D题的圈养湖羊空间利用率问题,包括问题分析、数学模型建立和MATLAB代码实现,旨在优化养殖场的生产计划和空间利用效率。
232 6
【2023高教社杯】D题 圈养湖羊的空间利用率 问题分析、数学模型及MATLAB代码
|
4月前
|
存储 算法 搜索推荐
【2022年华为杯数学建模】B题 方形件组批优化问题 方案及MATLAB代码实现
本文提供了2022年华为杯数学建模竞赛B题的详细方案和MATLAB代码实现,包括方形件组批优化问题和排样优化问题,以及相关数学模型的建立和求解方法。
143 3
【2022年华为杯数学建模】B题 方形件组批优化问题 方案及MATLAB代码实现
|
4月前
|
数据采集 存储 移动开发
【2023五一杯数学建模】 B题 快递需求分析问题 建模方案及MATLAB实现代码
本文介绍了2023年五一杯数学建模竞赛B题的解题方法,详细阐述了如何通过数学建模和MATLAB编程来分析快递需求、预测运输数量、优化运输成本,并估计固定和非固定需求,提供了完整的建模方案和代码实现。
111 0
【2023五一杯数学建模】 B题 快递需求分析问题 建模方案及MATLAB实现代码
|
7月前
|
数据可视化
Matlab模拟森林火灾传播
Matlab模拟森林火灾传播
|
7月前
|
数据安全/隐私保护
耐震时程曲线,matlab代码,自定义反应谱与地震波,优化源代码,地震波耐震时程曲线
地震波格式转换、时程转换、峰值调整、规范反应谱、计算反应谱、计算持时、生成人工波、时频域转换、数据滤波、基线校正、Arias截波、傅里叶变换、耐震时程曲线、脉冲波合成与提取、三联反应谱、地震动参数、延性反应谱、地震波缩尺、功率谱密度
基于混合整数规划的微网储能电池容量规划(matlab代码)
基于混合整数规划的微网储能电池容量规划(matlab代码)
|
7月前
|
数据安全/隐私保护
地震波功率谱密度函数、功率谱密度曲线,反应谱转功率谱,matlab代码
地震波格式转换、时程转换、峰值调整、规范反应谱、计算反应谱、计算持时、生成人工波、时频域转换、数据滤波、基线校正、Arias截波、傅里叶变换、耐震时程曲线、脉冲波合成与提取、三联反应谱、地震动参数、延性反应谱、地震波缩尺、功率谱密度
|
7月前
|
算法 调度
面向配电网韧性提升的移动储能预布局与动态调度策略(matlab代码)
面向配电网韧性提升的移动储能预布局与动态调度策略(matlab代码)
|
7月前
|
算法 调度
含多微网租赁共享储能的配电网博弈优化调度(含matlab代码)
含多微网租赁共享储能的配电网博弈优化调度(含matlab代码)
|
7月前
|
Serverless
基于Logistic函数的负荷需求响应(matlab代码)
基于Logistic函数的负荷需求响应(matlab代码)

热门文章

最新文章