基于SEIR微分方程模型对疫情传播的简单预测

这篇具有很好参考价值的文章主要介绍了基于SEIR微分方程模型对疫情传播的简单预测。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

目录

一、模型的建立

传染病模型概念

模型假设

SEIR模型

模型中涉及的函数S(t)、E(t)、I(t)、R(t)

更改后的微分方程

二、模型的求解

三、模型的缺点

祝语


随着疫情的再次爆发,全国疫情防控再次进入紧张状态,疫情预测分析成为数学建模问题中的一个热点问题,本文基于微分方程的SEIR模型对疫情做出简单预测。


提示:以下是本篇文章正文内容,下面案例可供参考

一、模型的建立

传染病模型概念

   传染病基本数学模型主要是研究传染病的动力学原理、传播方式、空间范围、传播速度等问题。常见的传染病模型按传染病类型分为SI、SIR、SEIR模型等,对疫情的研究,本文选取SEIR模型,相比较SIR模型,能更好的预测疫情变化趋势。

模型假设

1.假设初始时刻易染者为总数N

2.假设病毒时间尺度远小于个体生命周期,即不考虑个体自然出生率和自然死亡率

3.建设各个个体间接触机会均等

SEIR模型

假设人群总体为N,定义四类人群,如下表所示。

                                                                表1  SEIR 模型的符号定义

符号 名称 解释
S 易染者 健康状态,可被感染的个体
E 暴露者  即与感染者接触过的个体(潜伏期)
I 感染者· 处于感染状态的个体还能够感染将康状态的个体
R 治愈者 恢复状态


在病毒最开始的时候S=N,然后S以每天α的速度变到E,E以每天σ的速度变到I,I又以每天β的速度变到R:

基于SEIR微分方程模型对疫情传播的简单预测        

SEIR模型对应的经典微分方程为:
                                                             

                                   

                                                              

                                                              

                                                基于SEIR微分方程模型对疫情传播的简单预测

模型中涉及的函数S(t)、E(t)、I(t)、R(t)

  S(t)的意思是第t天健康个体的数量,E(t)是第t天接触过感染者的个体数量,I(t)是第t天感染个体的数量,R(t)是第t天免疫个体的数量,N(t)是整个种群的数量,在假设情况下固定不变为N,取。

  该模型为经典传染病模型,由于疫情感染因素较普通传染病较为复杂,且具有潜伏期,因此引入三个额外的参数,,a.,,a分别为潜伏期自愈转化为康复者,即接触过感染者但未被感染,患病者治愈后转为康复者的概率,潜伏者对易感人群的传染率。

更改后的微分方程

                                                              

                                                     基于SEIR微分方程模型对疫情传播的简单预测              

                                                      

                                                      基于SEIR微分方程模型对疫情传播的简单预测

                                                      基于SEIR微分方程模型对疫情传播的简单预测

二、模型的求解

         对于模型的求解,我们需要拟合六个参数:即易感者初始值S(0),潜伏者初始值E(0),传染率β,潜伏率a,康复率。查找近七天的数据,运用MATLAB拟合工具箱对数据进行拟合。

     图1 全国4月10至16日确诊人数及拟合

基于SEIR微分方程模型对疫情传播的简单预测


得到I(t)的拟合表达式:

                                                         基于SEIR微分方程模型对疫情传播的简单预测

继而得出:

                                                                    

下面拟合近七天的治愈人数:

                                                                   图二 近七天治愈人数

基于SEIR微分方程模型对疫情传播的简单预测

 可以看到拟合效果不佳,可能由于科技,防控和疫情变化因素均很大,我们缩短时间尺度

基于SEIR微分方程模型对疫情传播的简单预测

 可以看到近六天的拟合效果明显好于近七天的,得出R(t)的函数表达式:

                                                                            基于SEIR微分方程模型对疫情传播的简单预测

                                                                             

    原方程组化为:

                                                                                       

                                                                基于SEIR微分方程模型对疫情传播的简单预测  

                                                                               

                                                                        基于SEIR微分方程模型对疫情传播的简单预测

    最后用MATLAB编程计算得出结果:         

                                                                                 

                                                                                  

                                                                                   

                                                                                   

                                                                                     a=0.05

                                                                       

        MATLAB得出最终预测结果:

基于SEIR微分方程模型对疫情传播的简单预测

 

 

三、模型的缺点

1.在SIR模型的基础上,SEIR模型在模型中加入了E, 具有潜伏期且潜伏期不具有传染性传染病适用 SEIR模型。新型冠状病毒在潜伏期也具有传染性,并不满足SEIR模型的适用条件,因此该模型存在较大局限性

2.疫情传播过程中人与人之间的接触机会并不均等

祝语

  新型肺炎看出人生百态,是对政府、民众、科研工作者最大的考题。第一时间公开数据信息、治疗方案、数理模型、预测结果等,是科研人最大的贡献。面对谣言和恐慌,科学分析、知识传递,也尤为重要。相信全国齐心协力,必能获得抗击新型肺炎的最终胜利。文章来源地址https://www.toymoban.com/news/detail-470492.html

到了这里,关于基于SEIR微分方程模型对疫情传播的简单预测的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 通过Matlab编程分析微分方程、SS模型、TF模型、ZPK模型的关系

    以最简单的单自由度振动模型为例: 以上表示u(t)线性组合输入系统(这里是3u(t))时求系统的响应(即输出函数y(t)) SS模型也可转成TF模型: tf(ss(A,B,C,D)) TF转零极点增益ZPK模型 [z p k]=tf2zp([3],[1 0 4]) z = Empty matrix: 0-by-1 p = 0 + 2.0000i 0 - 2.0000i k = 3 即 还可以用residue函数将传递函数

    2024年02月11日
    浏览(42)
  • 美赛BOOM数学建模4-2微分方程传染病预测模型

    注明:本文根据数学建模BOOM网课简单整理,自用 ❑从最简单的指数传播模型说起 • 不同类型传染病的发病机理和传播途径各有特点 • 有的传染病,在得过一次后可获得 免疫力 ,但有的则不会 • 有的传染病具有 潜伏期 ,有的则没有 • 需要对不同类型的传染病建立相应

    2024年02月08日
    浏览(47)
  • 【数学建模】常用微分方程模型 + 详细手写公式推导 + Matlab代码实现

    微分方程基本概念 微分方程在数学建模中的应用 微分方程常用模型(人口增长模型、传染病模型) 2022.06.19 微分方程,是指含有未知函数及其导数的关系式。解微分方程就是找出未知函数。 微分方程是伴随着微积分学一起发展起来的。微积分学的奠基人Newton和Leibniz的著作中

    2024年02月09日
    浏览(67)
  • 基于MATLAB的微分方程的解析解与欧拉算法的数值解(附完整代码)

    正常的求解微分方程的MATLAB格式如下: 如果需要指明自变量,则如下: 格式中的 fi 既可以描述微分方程,又可以描述 初始条件 或 边界条件 。 描述微分方程的MATLAB格式为: D4y=7 ; 描述条件的MATLAB格式为: D2y(2)=3 ; 输入信号u(t)如下: 求解如下微分方程的通解 解: 此题需

    2023年04月09日
    浏览(47)
  • 0702可分类变量的微分方程-微分方程

    本节至第四节我们学习的都是一阶微分方程 ​ y ′ = f ( x , y ) y^{\\\'}=f(x,y) y ′ = f ( x , y ) (2-1) 一阶微分方程对称形式 p ( x , y ) d x + Q ( x , y ) d y = 0 ( 2 − 2 ) p(x,y)dx+Q(x,y)dy=0qquad (2-2) p ( x , y ) d x + Q ( x , y ) d y = 0 ( 2 − 2 ) 若以x为自变量,y为因变量,则 d y d x = − P ( x , y ) Q (

    2024年02月04日
    浏览(52)
  • 常微分方程建模R包ecode(一)——构建常微分方程系统

    常微分方程在诸多研究领域中有着广泛应用,本文希望向大家介绍笔者于近期开发的R包 ecode ,该包 采用简洁易懂的语法帮助大家在R环境中构建常微分方程 ,并便利地调用R图形接口,研究常微分方程系统的相速矢量场、平衡点、稳定点等解析性质,或进行数值模拟,进行敏

    2024年02月16日
    浏览(44)
  • 【数学建模】常微分,偏微分方程

    普通边界   已知t0时刻的初值    ode45()  龙格-库塔法 一阶,高阶都一样 如下: s(1) = y , s(2)=y\\\'  s(3) = x , s(4)=x\\\'   分段边界 非匿名函数    手写改进的ode45()函数代码 复杂边界值(即已知初始值,也知道末尾值),用bvp4c()函数 1. pdepe()函数 椭圆-抛物线型 控制方程  左边界

    2024年02月09日
    浏览(45)
  • (矩阵)一阶微分方程和伯努利方程

    伯努利方程的标准形式: 伯努利方程解法: 方程两边同时除以y的n次, 做变量替换y-z: 转换为线性微分方程: 最后换回原来的变量即可得到伯努利方程。 一阶线性微分方程的标准形式: 当Q(x)=0,为齐次方程;当Q(x)≠0,为非齐次方程。 已知如下矩阵,求解一阶线性微分方

    2024年02月05日
    浏览(58)
  • 高等数学(微分方程)

    x y ′ ′ ′ + ( y ′ ) 3 + y 4 xy\\\'\\\'\\\'+(y\\\')^3+y^4 x y ′′′ + ( y ′ ) 3 + y 4 quad quad 三阶 y ′ = 2 x y\\\'=2x y ′ = 2 x quad quad quad quad quad quad 一阶 d y = 2 x d x dy=2xdx d y = 2 x d x quad quad quad quad 一阶 ( y ′ ′ ) 5 + 2 y ′ = 3 (y\\\'\\\')^5+2y\\\'=3 ( y ′′ ) 5 + 2 y ′ = 3 quad quad quad 二阶 quad 例1: 已知

    2024年02月10日
    浏览(48)
  • matlab解微分方程

    f=@(变量) 表达式; x1为2 3 4 5;x2为3 4 5 6的情况下求解函数f的值 用“dsolve” step1: 申明自变量和因变量 syms y(x) step2:编程 得到: step1: 申明自变量和因变量 syms y(x) step2:编程 得到 step1.写函数文件 step2.主函数 相当于定义了一个新向量y,然后列 匿名函数 ,方程的 左边都是一阶

    2024年02月13日
    浏览(81)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包