元胞自动机对传染病模型的模拟
人群种类分为:病人,潜伏期患者,移出者,群体移出者,死亡者和健康者。
病毒传播模型的动态模拟
在治愈率相对较低与防控措施不齐全的前提下,通过Matlab对病情的传播进行了模拟。
图 五 7病毒传播动态模拟图(1)
病情初始,病毒在人群中爆炸式扩散,大部分人群被病毒感染。
% 在Matlab中模拟病毒的传播模型
% 数值的建立
M=100; % 尺寸大小
beta=0.05; % 感染率
gama=0.01; % 免疫能力
% 建立网格
y = zeros(M, M); % 网格x元素的意思是: 0是健康人,1使感染者,2免疫者
% 设置初始网格x,在网格的中心有一圈感染者,半径为10个细胞
for i=1:Mfor j=1:Mdxx = i-M/2;dyy = j-M/2;d = sqrt(dxx*dxx+dyy*dyy);if ( d<10 )y(i,j)=1;endend
end
% 定义社区,也就是最近的8个邻居
lingju = [-1 -1; 0 -1; 1 -1; 1 0; 1 1; 0 1; -1 1; -1 0];
% 创造一个新的窗口
figure
hold on
% 主循环,迭代时间变量t
for t=1:100000% 遍历网格 x 中的所有单元格,对于索引 i 1从n 和 j 从1到 nfor i=1:Mfor j=1:M% 在邻居之间来回走动,传播疾病for k=1:8i2 = i+lingju(k, 1);j2 = j+lingju(k, 2);% 检查细胞是否在网格边界内if ( i2>=1 && j2>=1 && i2<=M && j2<=M )%如果细胞处于易感状态和邻近细胞% 被感染的传播感染的概率是betaif ( y(i,j)==0 && y(i2, j2)==1 )if ( rand<beta )y(i,j) = 1;endendendend% 如果被感染的人能够以伽马的概率从疾病中康复if ( y(i,j)==1 && rand<gama )y(i,j) = 2;endendend% 动态模拟clfimagesc(y, [0 2]) % 网格展示pause(0.01) % 暂停0.01scolormap([1 0 1; 1 1 0; 1 1 1]); % 定义0、1、2分别对应的颜色% 如果没有更多的感染者,就停止模拟if ( sum(y==1)==0 )break;end
end
%
元胞自动机对传染病模型的模拟相关推荐
- 学习常用模型及算法:2.元胞自动机和交通模型
元胞自动机在数模中的应用十分广泛. 元胞自动机原理 一维元胞自动机.给出任意一个状态,都能知道下一时刻的状态.规则已给出,总共有2^3 = 8种可能. 二维元胞自动机--生命游戏 元胞的状态是有限的- ...
- 元胞自动机交通流模型c++_MATLAB——含出入匝道的交织区快速路元胞自动机模型...
输出效果 本模型为包含一条变速车道和出入匝道的快速路交织区仿真模型. 车道编号 换道规则和不含出入匝道的路段的区别在于换道分为自由换道和强制换道.车辆在不同区域有不同换道状态,而车辆本身又有出匝道和直 ...
- Matlab基于元胞自动机实现NaSch模型
一.NaSch模型 二.部分代码 % 主程序: NaSch_3.m程序代码 % 单车道最大速度 3个元胞开口边界条件加速减速随机慢化 clf clear all %build the GUI %def ...
- 基于元胞自动机的森林火灾matlab模拟
先建立一个300*300的矩阵,模拟森林的各种状态,其中矩阵的每个元素有3个状态,分别是0.1.2,其中0代表此处没有树木生长,1代表此处的树木正在燃烧,2代表此处有树木. 这个模型的规则是, A.如 ...
- 元胞自动机模型_【ABM仿真模拟】第三章 元胞自动机 B
是新朋友吗?记得先点蓝字关注我哦- 第三章 元胞自动机 B 2020/01/23 3.3 投票模型(Voting) 课程导读 无论是国内国外,投票结果的预测都异常受到关注.在众多唐斯模型.中间选民模型 ...
- 元胞自动机:更接近人类思考的智能模型
前段时间,一则"计算机首次通过图灵测试"的新闻被众多媒体报道,似乎人工智能已然离我们不远了. 在未来,计算机真的会拥有人类一样的自主思考.学习能力吗,甚至出现未来的"机器 ...
- python元胞自动机模拟交通_结构专栏 | 解析DEFORM软件中的元胞自动机法
点击上方蓝色字体,关注我们 导语 金属材料的性能取决于内部的微观组织结构,而好的材料性能和价格是产品最大的优势.随着现代物理冶金.热成形技术.热处理技术和计算机技术的兴起与发展,使预测和控制金属材料热 ...
- 【元胞自动机】基于元胞自动机模拟大型商场人流疏散含Matlab源码
1 简介 基于元胞自动机的场强模型在二维平面行人流疏散问题的研究中已得到了广泛应用.已有模型主要描述行人基于出口位置并跟随其余行人进行疏散的行为特征,未充分考虑火灾蔓延和局部拥堵对行为选择的影响,难以 ...
- 元胞自动机模拟森林火灾--matlab实现
因为建模需要,学习了元胞自动机,并尝试了用其模拟森林火灾,相对于退火和遗传,元胞自动机更容易理解,并且其广泛应用于交通问题,在以往的建模比赛中占了很大的比重. 一.元胞自动机定义和适用范围 不同于一般 ...
最新文章
- 【原创】多台电脑如何公用一个键盘鼠标进行操作 - Mouse with Borders 软件
- cdn加速的原理是怎么样的
- 尼康d7200拍照_尼康D7500适合入门吗
- 国际农民丰收节贸易会-万祥军: 全球TOP10农企携手阿里
- Idea左侧显示目录结构和.java文件的解决办法
- Flink进行Kafka事实表与Mysql维度表Join(纯DDL/DML方式)
- 【加密算法】Base64
- 北美互联网公司“隐藏福利”大盘点:看完我彻底酸了...
- js中避免函数名和变量名跟别人冲突
- mysql+installer+community+5.7.9_win10系统,mysql-installer-community-5.7.19.0.msi安装
- md5后得到的32位字符串存储到mysql中太占空间了_面试官:你对MySQL高性能优化有什么规范建议?...
- 首都师范 博弈论 5 4 2 Shapley值应用案例
- 超火动态表白H5源码
- 如何交叉编译fio并移植到ARM、IOT上
- 一元稀疏多项式计算器
- 【模型检测学习笔记】8:无限字上ω正则LT性质的验证
- 基于改进YOLOv7和CRNN的管道裂缝检测系统(源码&教程)
- 【HUSTOJ】1022: 输出小写字母表
- 2018中南大学 计算机考研分数,2018中南大学考研分数线多少分?
- 基于安卓的校园订餐系统开发设计
热门文章
- matlab中数组的声明,如何使用Matlab声明矩阵数组
- Retrofit2.0初始化
- 小白终是踏上了这条不归路----小文的mysql学习笔记目录
- 一个漂亮的证明与作图:高斯的正十七边形
- 近3000款长久未更新App或将被苹果下架
- 计算机二级加拼音,计算机二级C语言上机题库超强拼音检索快速找到(多个版本,综合整理,十分强大)最好配合第二套一起使用...
- 工程技术TOP期刊《Engineering》观点评述:基于脑电图的脑-机接口系统在实用化进程中面临的挑战...
- [实训题目EmoProfo]基于深度学习的表情识别服务搭建(一)
- omniplan导出html,OmniPlan 3.14.4 最NB的项目管理流程软件
- 通过U盘的HDD启动安装原版Win7系统教程-Win7安装版