【元胞自动机】基于matlab元胞自动机模拟SEIR传播模型和采取隔离措施的SEIR模型【含Matlab源码 2181期】
⛄一、元胞自动机简介
1 元胞自动机发展历程
最初的元胞自动机是由冯 · 诺依曼在 1950 年代为模拟生物 细胞的自我复制而提出的. 但是并未受到学术界重视.1970 年, 剑桥大学的约翰 · 何顿 · 康威设计了一个电脑游戏 “生命游戏” 后, 元胞自动机才吸引了科学家们的注意.1983 年 S.Wolfram 发表了一系列论文. 对初等元胞机 256 种 规则所产生的模型进行了深入研究, 并用熵来描述其演化行 为, 将细胞自动机分为平稳型, 周期型, 混沌型和复杂型.
2 对元胞自动机的初步认识
元胞自动机(CA)是一种用来仿真局部规则和局部联系的方法。典型的元胞自动机是定义在网格上的,每一个点上的网格代表一个元胞与一种有限的状态。变化规则适用于每一个元胞并且同时进行。典型的变化规则,决定于元胞的状态,以及其( 4 或 8 )邻居的状态。
3 元胞的变化规则&元胞状态
典型的变化规则,决定于元胞的状态,以及其( 4 或 8 )邻居的状态。
4 元胞自动机的应用
元胞自动机已被应用于物理模拟,生物模拟等领域。
5 元胞自动机的matlab编程
结合以上,我们可以理解元胞自动机仿真需要理解三点。一是元胞,在matlab中可以理解为矩阵中的一点或多点组成的方形块,一般我们用矩阵中的一点代表一个元胞。二是变化规则,元胞的变化规则决定元胞下一刻的状态。三是元胞的状态,元胞的状态是自定义的,通常是对立的状态,比如生物的存活状态或死亡状态,红灯或绿灯,该点有障碍物或者没有障碍物等等。
6 一维元胞自动机——交通规则
定义:
6.1 元胞分布于一维线性网格上.
6.2 元胞仅具有车和空两种状态.
7 二维元胞自动机——生命游戏
定义:
7.1 元胞分布于二维方型网格上.
7.2 元胞仅具有生和死两种状态.
元胞状态由周围八邻居决定.
规则:
骷髅:死亡;笑脸:生存
周围有三个笑脸,则中间变为笑脸
少于两个笑脸或者多于三个,中间则变死亡。
8 什么是元胞自动机
离散的系统: 元胞是定义在有限的时间和空间上的, 并且元 胞的状态是有限.
动力学系统: 元胞自动机的举止行为具有动力学特征.
简单与复杂: 元胞自动机用简单规则控制相互作用的元胞 模拟复杂世界.
9 构成要素
(1)元胞 (Cell)
元胞是元胞自动机基本单元:
状态: 每一个元胞都有记忆贮存状态的功能.
离散: 简单情况下, 元胞只有两种可能状态; 较复杂情况下, 元胞具有多种状态.
更新: 元胞的状态都安照动力规则不断更新.
(2)网格 (Lattice)
不同维网格
常用二维网格
(3)邻居 (Neighborhood)
(4)边界 (Boundary)
反射型:以自己作为边界的状态
吸收型:不管边界(车开到边界就消失)
(5)规则(状态转移函数)
定义:根据元胞当前状态及其邻居状况确定下一时刻该元胞状态的动力学函数, 简单讲, 就是一个状态转移函数.
分类 :
总和型: 某元胞下时刻的状态取决于且仅取决于它所有邻居 的当前状态以及自身的当前状态.
合法型: 总和型规则属于合法型规则. 但如果把元胞自动机 的规则限制为总和型, 会使元胞自动机具有局限性.
(6)森林火灾
绿色:树木;红色:火;黑色:空地。
三种状态循环转化:
树:周围有火或者被闪电击中就变成火。
空地:以概率p变为树木
理性分析:红为火;灰为空地;绿是树
元胞三种状态的密度和为1
火转化为空地的密度等于空地转换为树的密度(新长出来的树等于烧没的树)
f是闪电的概率:远远小于树生成的概率;T s m a x T_{smax}T smax
是一大群树被火烧的时间尺度
程序实现
周期性边界条件
购进啊
其中的数字为编号
构建邻居矩阵
上面矩阵中的数字编号,对应原矩阵相同位置编号的上邻居编号,一 一对应
同样道理:
(7)交通概念
车距和密度
流量方程
守恒方程
时空轨迹(横轴是空间纵轴为时间)
红线横线与蓝色交点表示每个时间车的位置。
如果是竖线则表示车子在该位置对应的时间
宏观连续模型:
最常用的规则:
红色条表示速度是满的。
1 加速规则:不能超过v m a x ( 2 格 / s ) v_{max}(2格/s)v
max(2格/s)
2 防止碰撞:不能超过车距
理论分析:
结果分析: 密度与流量
第一个图:横坐标是归一化后的密度,纵坐标是车流量。第二个图:理论值与CA的结果
结果分析: 时空轨迹
中间的深色区域是交通堵塞的区域。
⛄二、部分源代码
%SEIR模型修正
%假设第5天开始采取隔离措施1月19日第一例确诊;1月23号武汉开始封城,此时其他省市也相应作出响应隔离措施
N=13026800;%人口数
E=0;%潜伏者
D=0;%死亡患者人数
I=1;%感染人数
S=N-I;%易感人数
R=0;%康复者人数
r=20;%感染者接触数量
B=0.6;%传染概率
a=0.17;%潜伏者转化为感染者的概率
r2=8;%潜伏者接触人数
B2=0.03;%潜伏者传染正常人的概率
y=0.05;%康复概率
k=0.0001;%日致死率
B3=0.1;%转阴率
T=1:200;
for idx=1:length(T)-1
if idx>=7
r=0.20;
r2=1;
y=0.1
end
if idx<7
S(idx+1)=S(idx)-rBS(idx)I(idx)/N-r2B2S(idx)E(idx)/N;%易感人群迭代
E(idx+1)=E(idx)+rBS(idx)I(idx)/N-aE(idx)+r2B2S(idx)E(idx)/N;%潜伏者迭代
I(idx+1)=I(idx)+aE(idx)-(k+y)I(idx);%感染人数迭代
R(idx+1)=R(idx)+yI(idx);%康复人数迭代
D(idx+1)=R(idx)+kI(idx);%死亡患者人数迭代
else
S(idx+1)=S(idx)-rBS(idx)I(idx)/N-r2B2S(idx)E(idx)/N+B3E(idx-6);%易感人群迭代
E(idx+1)=E(idx)+rBS(idx)I(idx)/N-aE(idx)+r2B2S(idx)E(idx)/N-B3E(idx-6);%潜伏者迭代
I(idx+1)=I(idx)+aE(idx)-(k+y)I(idx);%感染人数迭代
R(idx+1)=R(idx)+yI(idx);%康复人数迭代
D(idx+1)=R(idx)+kI(idx);%死亡患者人数迭代
end
end
⛄三、运行结果
⛄四、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1]党珊,蒋太刚,巫承军.基于元胞自动机方法的消防疏散仿真研究[J].现代电子技术. 2022,45(21)
[2]帅斌,秦梦瑶,许旻昊.基于元胞自动机的高速铁路列车运行仿真研究[J].计算机仿真. 2022,39(08)
[3]张睿洋.元胞自动机在两类传染病模型中的应用[J].现代信息科技. 2022,6(10)
3 备注
简介此部分摘自互联网,仅供参考,若侵权,联系删除
【元胞自动机】基于matlab元胞自动机模拟SEIR传播模型和采取隔离措施的SEIR模型【含Matlab源码 2181期】相关推荐
- 【Matlab水果识别】苹果质量检测及分级系统(带面板)【含GUI源码 1613期】
一.代码运行视频(哔哩哔哩) [Matlab水果识别]苹果质量检测及分级系统(带面板)[含GUI源码 1613期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1 ...
- 【Matlab生物电信号】生物电信号仿真【含GUI源码 684期】
一.代码运行视频(哔哩哔哩) [Matlab生物电信号]生物电信号仿真[含GUI源码 684期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1]董兵,超于毅,李 ...
- 【Matlab语音分析】语音信号分析【含GUI源码 1718期】
一.代码运行视频(哔哩哔哩) [Matlab语音分析]语音信号分析[含GUI源码 1718期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1]韩纪庆,张磊,郑铁 ...
- 【Matlab验证码识别】遗传算法和最大熵优化+大津法(OTSU)+自定义阈值数字验证码识别【含GUI源码 1694期】
一.代码运行视频(哔哩哔哩) [Matlab验证码识别]遗传算法和最大熵优化+大津法(OTSU)+自定义阈值数字验证码识别[含GUI源码 1694期] 二.matlab版本及参考文献 1 matlab ...
- 【Matlab人脸识别】BP神经网络人脸识别(含识别率)【含GUI源码 891期】
一.代码运行视频(哔哩哔哩) [Matlab人脸识别]BP神经网络人脸识别(含识别率)[含GUI源码 891期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1] ...
- 【Matlab人脸识别】形态学教室人数统计(带面板)【含GUI源码 1703期】
一.代码运行视频(哔哩哔哩) [Matlab人脸识别]形态学教室人数统计(带面板)[含GUI源码 1703期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1]孟 ...
- 【Matlab人脸识别】人脸实时检测与跟踪【含GUI源码 673期】
一.代码运行视频(哔哩哔哩) [Matlab人脸识别]人脸实时检测与跟踪[含GUI源码 673期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1]孟逸凡,柳益君 ...
- 【Matlab图像融合】小波变换遥感图像融合【含GUI源码 744期】
一.代码运行视频(哔哩哔哩) [Matlab图像融合]小波变换遥感图像融合[含GUI源码 744期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1] 包子阳,余 ...
- 【Matlab语音加密】语音信号加密解密(带面板)【含GUI源码 181期】
一.代码运行视频(哔哩哔哩) [Matlab语音加密]语音信号加密解密(带面板)[含GUI源码 181期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1]韩纪庆 ...
- 【Matlab身份证识别】身份证号码识别【含GUI源码 014期】
一.代码运行视频(哔哩哔哩) [Matlab身份证识别]身份证号码识别[含GUI源码 014期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1] 蔡利梅.MAT ...
最新文章
- Swift学习:自动引用计数
- android编译的tool版本有多少,macOS Mojave(10.14.2)系统上编译LunarG/VulkanTools工程的Android版本...
- linux 指令tftp传输文件_Linux tftp命令使用详解:在本机与tftp服务器之间使用TFTP协议传输文件...
- 学计算机怎么收费呢,学计算机怎么收费服务为先_贵阳经济技术学校
- 炫酷!用Python制作漂亮的流动桑基图
- eclipse4.3.1标准版安装freemarker插件
- 计算机学业水平测试字处理多少分,【计算机应用论文】计算机应用基础学业水平的测试问题(共3624字)...
- Python 进程池 multiprocessing.Pool - Python零基础入门教程
- Vue 脚手架配置代理
- python3字符串常用操作
- Facebook将用实际行动引导青少年远离有害内容
- 如何在ASP.NET Core中建立有效的分页
- HDU 4539 郑厂长系列故事——排兵布阵 —— 状压DP
- 盘点国内外私募基金业绩报酬计提方式
- java读CSV文件入数组
- 图片转字符画,GIF动画转GIF字符动画
- Uva 11137 Ingenuous Cubrency(整数划分方案 背包)
- 科学的工作方法—PDCA循环工作法
- python爬虫——正则表达式
- 掌握这些Revit常见术语,从此BIM建模再无压力