【SEIRS流行病模型】通过经典确定性模型模拟流行病:SI,SIS,SIR,SIRS,SEIR,SEIR,SEIRS(Matlab代码实现)

 👨‍🎓个人主页:研学社的博客    

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势:🌞🌞🌞博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。

⛳️座右铭:行百里者,半于九十。

📋📋📋本文目录如下:🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、数据、文档


💥1 概述

过去几十年生态学在连续动力系统的研究结果已相当丰富.在连续动力学中总是假设被描述的生命现象是连续发生的,比如动物的生育、死亡,或是两个相对独立的环境中的迁移过程,还有人类对可更新自然资源的管理和利用方面,如对一些人工养殖的动物的放养和收获,以及在传染病控制方面实
施的人工接种免疫等行为都在数学模型里被描述为连续发生的.一般来说,这些动物的生理行为以及人类的干扰行为本身并不是连续的,而是离散的·疫苗接种是控制疾病传播的很有效的途径,数学模型对接种策略的设计和估计起到帮助作用.如果成功接种人口的比例高于某一个临界值时传统的接种策略会导致疾病根除.例如麻疹的接种比例约等于95%,而在实践中,对人口密集地区实施接种是又困难,花费又高.

通过该类流行病,人们可以通过几个经典的确定性流行病模型模拟对流行病的评估,而无需生命动力学:
SI -- 易感>感染性
SIS -- 易感性>感染性>易感性 SIR -- 易感>感染性去除
SIRS -- 易感>感染性>移除>易感 SEIR -- 易感

>暴露>感染->移除
SEIRS -- 易感->暴露->感染->去除->易感

基本术语是:混合良好的人口,被划分为几个不连贯的阶层;流行病,即控制个人在阶层之间过渡的过程。经典模型包括以下类别:•S类易感个体,即那些可能变得具有传染性的人•E类暴露个体,即那些被感染但尚未具有传染性的人•I类感染个体,即那些可以传播疾病的人•R类移除个体,即那些不能再传播疾病的人。类别的数量取决于模型,但所有模型都包括S类和I类。

所有实现的模型都是确定性的,即流行病动力学是由微分方程给出的。

因此,人口是封闭的,被视为一个连续体,而流行病的进程在时间上是连续的。每个类都是同质的,具有时间无关的性质,即一个类中的所有个体都是平等的;因此,类之间的转换速率是常数。

请注意,经典的流行病模型描述的是个体在阶级之间的动态过渡,而不是阶级内部的过程。

 

📚2 运行结果

 

 

 

 

 

部分代码:

Amin = 5;  % data less than Amin will be deleted from file begining

% get table 
T = readtable(fileXLSX);

% get data
A =  T{:,2:end};
[nrow,ncol] = size(A);

% table column names
name = string(T.Properties.VariableNames);
name = name(2:end)';

% start date
date0 = datenum(T{1,1}); %datenum(txt{2,1},'dd.mm.yyyy'); %datenum('2019/12/31');

% end date
%date1 = date0 + nrow - 1;

%functions names
ffname = strings(ncol,1);
nn = 0;
for n = 1:ncol
    nname = name{n};
    if strcmp("",nname) 
        continue
    end
    nn = nn + 1;
    nname = strrep(nname,' ','_');
    nname = strrep(nname,'-','_');  
    nname = strrep(nname,'''','_');      
    nname = strrep(nname,'(','_');       
    nname = strrep(nname,')','_');      
    fname = sprintf('getData%s.m',nname);
    fid = fopen(fullfile(path,fname),'w');
    if fid < 0
        fprintf('***Fail to open %s\n',fname);
        continue
    end

🎉3 参考文献

部分理论来源于网络,如有侵权请联系删除。

[1]宋宇恒,仲崇阳,韩七星.非线性SEIR流行病模型的平稳分布[J].东北师大学报(自然科学版),2022,54(01):14-18.DOI:10.16163/j.cnki.22-1123/n.2022.01.004.

[2]罗芬,向中义.一类SEIRS流行病模型的研究[J].湖北民族学院学报(自然科学版),2008(02):141-145.

[3]徐文雄,张太雷.一类非线性SEIRS流行病传播数学模型[J].西北大学学报(自然科学版),2004(06):627-630.

🌈4 Matlab代码、数据、文档

  • 1
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

荔枝科研社

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值