具有疫苗接种的COVID-19传播模型的动力学分析.pdf
《具有疫苗接种的COVID-19传播模型的动力学分析.pdf》由会员分享,可在线阅读,更多相关《具有疫苗接种的COVID-19传播模型的动力学分析.pdf(10页珍藏版)》请在咨信网上搜索。
1、Advances in Applied Mathematics 应用数学进展应用数学进展,2024,13(4),1187-1196 Published Online April 2024 in Hans.https:/www.hanspub.org/journal/aam https:/doi.org/10.12677/aam.2024.134109 文章引用文章引用:赵勇盛,张蒙.具有疫苗接种的COVID-19传播模型的动力学分析J.应用数学进展,2024,13(4):1187-1196.DOI:10.12677/aam.2024.134109 具有疫苗接种的具有疫苗接种的COVID-19传
2、播模型的动力学传播模型的动力学分析分析 赵勇盛赵勇盛*,张,张 蒙蒙 北京建筑大学理学院,北京 收稿日期:2024年3月9日;录用日期:2024年4月1日;发布日期:2024年4月10日 摘摘 要要 自自COVID-19出现以来,毒株不断变异,毒性变小、传播性增强,多数专家认为我们将与出现以来,毒株不断变异,毒性变小、传播性增强,多数专家认为我们将与COVID-19长期长期共存。由于是否出现重症患者很大程度取决于个体免疫力,因此处于高危人群的老人和儿童长期面临威共存。由于是否出现重症患者很大程度取决于个体免疫力,因此处于高危人群的老人和儿童长期面临威胁。此外由于没有治疗胁。此外由于没有治疗CO
3、VID-19的特效药物,加之多次感染新冠会对人体带来诸多后遗症,因此许多国的特效药物,加之多次感染新冠会对人体带来诸多后遗症,因此许多国家都在研发和接种针对最新变异毒株可以有效减少患重症和死亡的新型疫苗。并且由于毒株变异和免疫家都在研发和接种针对最新变异毒株可以有效减少患重症和死亡的新型疫苗。并且由于毒株变异和免疫逃逸,我们将持续研发有效且快速生效的新型疫苗。因此本文提出了一个具有疫苗接种时滞的逃逸,我们将持续研发有效且快速生效的新型疫苗。因此本文提出了一个具有疫苗接种时滞的SVIAR传传染病模型,计算了模型的控制再生数,研究了模型的动力学特征。文末的数值模拟验证了理论结果,并染病模型,计算了
4、模型的控制再生数,研究了模型的动力学特征。文末的数值模拟验证了理论结果,并且分析了疫苗接种率和疫苗生效时间对疫情防控的影且分析了疫苗接种率和疫苗生效时间对疫情防控的影响,给出了对应的疫苗接种策略。响,给出了对应的疫苗接种策略。关键词关键词 COVID-19,疫苗接种疫苗接种,时滞微分系统时滞微分系统,全局稳定性全局稳定性,数值模拟数值模拟 Dynamical Analysis of a COVID-19 Transmission Model with Vaccination Yongsheng Zhao*,Meng Zhang School of Science,Beijing Univers
5、ity of Civil Engineering and Architecture,Beijing Received:Mar.9th,2024;accepted:Apr.1st,2024;published:Apr.10th,2024 Abstract Since the emergence of COVID-19,the viral strains have continuously mutated,with reduced viru-lence and enhanced transmissibility.Most experts believe that we will coexist w
6、ith COVID-19 in *通讯作者。赵勇盛,张蒙 DOI:10.12677/aam.2024.134109 1188 应用数学进展 the long term.As the severity of infection largely depends on individual immunity,older adults and children in high-risk groups face prolonged threats.Additionally,due to the lack of specific drugs for treating COVID-19,and consid
7、ering that multiple infections with the novel coronavirus can lead to various lingering effects,many countries are developing and administering novel vac-cines targeting the latest mutant strains to effectively reduce severe cases and fatalities.Further-more,due to ongoing viral mutations and immune
8、 evasion,we will continue to research and de-velop efficient and rapidly effective new vaccines.Therefore,this article proposes an SVIAR infec-tious disease model with vaccination delays,calculates the models basic reproduction number,investigates its dynamic characteristics,and validates theoretica
9、l results through numerical si-mulations.The analysis also explores the impact of vaccination rates and vaccine efficacy duration on epidemic control,providing corresponding vaccination strategies.Keywords COVID-19,Vaccination,Delay Differential System,Global Stability,Numerical Simulation Copyright
10、 2024 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 自 2019 年 COVID-19 首次出现以来,截止到目前已经出现了众多的变异毒株,其中 JN.1 变异株作为全球占比迅速增加的流行毒株,截止到 2023 年 12 月其占比已经高达 30%,被全球超过至少 59 个国家和地区检
11、测到,主要分布在欧美地区。目前我国主流毒株仍为 EG.5 及其亚分支,但感染最新的 JN.1 毒株的人数也在增加。JN.1 毒株作为新冠病毒奥密克戎 BA.2.86 变异株的第二代亚分支,其毒性减弱、传播性更强,这也是一种自然选择的现象。因此有不少专家表示我们应当习惯于与新冠病毒长期共存。由于感染者症状的类型和严重程度通常更多地取决于个体免疫力,并且一项研究表明既往新冠感染与多种症状的患病率和严重程度增加有关1。因此作为免疫力较为低下的高风险人群,老人和小孩仍有患重症以及诸多后遗症的风险。鉴于目前没有可以治疗 COVID-19 的特效药物。疫苗接种成为了我们疫情防控的重要手段。因此全球各个国家
12、和地区都在致力于疫苗的研发和接种,并不断扩大疫苗接种覆盖率以此建立群体免疫2。根据最新声明,国际药品监管当局联盟(ICMRA)成员和世卫组织 COVID-19 疫苗组成技术咨询小组(TAG-CO-VAC)都认为 SARS-CoV-2 病毒似乎正在进化与之前病毒不同。因此有必要更新COVID-19 疫苗的组成。此外欧洲药品管理局去年建议批准了针对奥密克戎 XBB.1.5 的改良 COVID-19疫苗3。综上所述,由于病毒的变异和免疫逃逸,疫苗的效力也受到一定影响4-9。因此,研发和接种新型疫苗是应对 COVID-19 大流行的紧迫任务。也可以有效保障老人儿童等高危人群的生命健康。在10 11 1
13、2中的学者构建不同的舱室模型来研究疫苗接种对传染病动力学的影响,这些模型都能反映出疫苗接种的对传染病防控的作用和意义。在13 14中的学者在认为 COVID-19 传播过程中无症状感染者对疫情传播起到了重要作用,因此模型中至少应该加入无症状感染者这一显著特征。考虑到COVID-19在疫情传播过程中的延迟现象,在15 16 17 18 19中的学者建立了具有时滞的COVID-19传播模型,研究了时滞对疫情传播的影响,使结果更加直接和准确。在20中的学者建立了具有疫苗接种时滞的 COVID-19 传播模型,表明不同的疫苗生效时间对疾病能否有效控制起到了重要作用。通过上述学者们的研究以及结合实际案例
14、,本文创建了一个具有无症状感染者以及疫苗接种的Open AccessOpen Access赵勇盛,张蒙 DOI:10.12677/aam.2024.134109 1189 应用数学进展 SVIAR 传染病模型,用于研究 COVID-19 的传播,并对有症状感染者和无症状感染者的传播率进行区分。为使模型更符合实际情况,在模型中引入了时滞,考虑了疫苗从接种到产生抗体的时间间隔对疾病传播的影响。本文行文结构如下,第二节中建立了 COVID-19 传播模型并列出每个参数所代表的含义;第三节中求出系统的无病平衡点和地方病平衡点;第四节利用 Lyapunov 方法证明了模型平衡点的全局动力学性质;第五节对
15、模型进行数值模拟以检验理论分析的正确性和可靠性;第六节归纳结论给出疫苗接种策略。2.模型建立模型建立 根据 COVID-19 在人群中的传播特性,我们建立了一个 SVIAR 的数学模型用于研究 COVID-19 疫苗接种,其中()S t表示易感者人群,()V t表示疫苗接种人群,()I t表示有症状感染者人群,()A t表示无症状感染者人群,()R t表示恢复者人群,本文构造 COVID-19 传播模型为:()()()()()()()()()()()()()()()()()()()()()()()()()()()()()d,dd,dd1,dd,dd.dS tS tS tI tA tS ttV
16、tS tV ttI tS tI tA tI tI ttA tS tI tA tA tA ttT tI tA tR tt=+=+=+=+(1)我们假设模型(1)中所有的参数都是正的,各个参数的定义参见表 1。Table 1.The definition of parameters in a model 表表 1.模型中参数的定义 参数 定义 人类自然出生率 疫苗接种率 感染率 无症状感染者传播概率的调节因子 自然死亡率 有症状感染者恢复率 无症状感染者恢复率 无症状感染者比率 设()5,0,CC+=是 Banach 空间,所有连续函数的上确界范数从区间,0映射到5+,定义:()T12345,:.
17、CC +=0 对于模型(1)的适定性和耗散性我们有如下定理:定理定理 2.1 系统(1)过任意初值C+的解()()()()()()()T12345,u tututututut=唯一存在,并且是非负赵勇盛,张蒙 DOI:10.12677/aam.2024.134109 1190 应用数学进展 的,在5+最终有界。证明证明 参考20 21中时滞微分方程的基本理论,系统(1)经过初值5+的解在其最大存在区间)0,T上是唯一的。下面我们证明经过初值5+的解在其最大存在区间)0,T是非负的。对于任意的()0,T以及任意小的0,可知模型:()()()()()()()()()()()()()()()()()
18、()()()()()()()()()()()()d,dd,dd1,dd,dd.dS tS tS tI tA tS ttV tS tV ttI tS tI tA tI tI ttA tS tI tA tA tA ttR tI tA tR tt=+=+=+=+=+(2)通过C+的解为()()()()()()()T12345,u tututututut=,基于时滞微分方程解对参数的连续依赖性20 21,(),u t在)0,上是一致存在的,并且关于参数连续。当()0,0iu=,()1,2,3,4,5i=,有()0,0iu,于是我们可以得到(),u t 0,()0,t。从而对于)0,t,有()u t 0
19、。由的任意性,我们可以得到()u t 0,()0,tT。接下来证明模型(1)的有界性,定义函数()()()()()()H tS tV tI tA tR t=+,对函数()H t进行求导可得:()()()()()()()()()()()()dddddddddddd.H tS tV tI tA tR tttttttS tV tI tA tR tH t+=+=取极限可以得到:()limtH t=所以有()u t最终是有界的。由时滞微分方程解的延拓定理22 23 24知T=。于是定理2.1得证。其中模型的可行域为:()5,:;0S V I A RRSVIAR+=+3.阈值动力学阈值动力学 很容易得到系
20、统(1)有一个无病平衡点()T0,0,0,0E =+,利用下一代矩阵的方法求出控制再生数:()()()()()1.cR=+赵勇盛,张蒙 DOI:10.12677/aam.2024.134109 1191 应用数学进展 接下来令(1)式右侧为零,并通过化简可以得到系统的地方病平衡点,因此当且仅当1cR 时,系统存在唯一的地方病平衡点()T,*ES VIA R=。其中:()()()()()()()()()()*,11SV +=+()()()()()()*11111,1cIkkR+=+()()()()()()*,.11AIRI +=+()()()()()()()11.1k+=+4.全局稳定全局稳定
21、在本节中我们将利用 Lyapunov 函数法来讨论模型(1)的无病平衡点和地方病平衡点的全局稳定性,即全局渐近稳定。设()u t是模型(1)经过任何一个初始值5+的解,定义()是关于模型(1)的-极限集。记()1ln,0g xxx x=。于是我们对于无病平衡点0E有如下的结果。定理定理 4.1 对任意的0,当1cR 时,无病平衡点0E在5+中是 GAS,当1cR=时,是全局吸引的。证明证明:首先,我们证明当1cR 时,0E是全局吸引的。由定理 2.1 可以得知,()5+是紧的,记:()51:00+=我们在上定义如下 V 函数:()()()()()10340VS gS=+(3)当1t 时,V 沿
22、着tu的导数为:()()()()()()()00210.tcSSV uS tI tA tRSS=+(4)结合(3)式和(4)式,我们可以得到()。因此对1cR,V 在:1tut 是 Lyapunov 函数。由25中推 3.3,可得()有()0V=。假设()()()()()()T,tuS t V tI tA tR t=是系统(1)过任意()的解。由于()的不变性,我们可以得出结论,对于所有的t,()tu。所以得到,对任意的t,()0tV u=。因此,从(6)可知,所有t的()0S tS=。因此,它认为()()()0I tA tR t=,0VV=当所有t。因此,()0E=。接下来,我们证明0E是稳
23、定的。由于系统(1)的第二和第五个方程是独立的,因此系统(1)的其余三个方程可以构成一个独立的子系统。通过(3)、(4)和(26,引理 4.1),我们可以证明该子系统的无病平衡 点()T0,0,0,0E =+是稳定的。利用(26,定理 4.3)中类似的论证,我们可以得到0E是稳定 的。因此,当1cR 时,无病平衡点0E在5R+中是全局渐近稳定的。下面我们论证地方病平衡点的*E的全局稳定性,以下是我们的结果:赵勇盛,张蒙 DOI:10.12677/aam.2024.134109 1192 应用数学进展 定理定理 4.2 对任意的0,当1cR 时,地方病平衡点*E在53:0MR+=中是 GAS。证
24、明证明:由定理 3.1 可得,当1cR。接下来我们证明当1cR 时,地方病平衡点*E在 M 中是全局渐进稳定的。设()u t是模型(1)经过任何一个初始值M的解。显然 M 是系统(1)的正向不变集,且(),0u tt 0。记52:+=0。显然2M。定义如下 V 函数:()()()*314*1IAVS gI gA gSIAIAIA=+(5)让模型(1)右侧为 0 并经过化简可以得到:()*IAS+=+()*SIAA+=()()*1SIAI+=当1t 时,V 沿着tu的导数为:()()()()()()()()()2*0.A t IA I tS tSVggS tS tSIAA t I t=+(6)当
25、且仅当()()()*,S tSI tIA tA=时成立。因此0V。当且仅当()()()()*,S tSI tIA tAR tR=时等号成立。通过(5)、(6)和系统(1)的第二、第五个方程,我们可以得到()。因此对1cR,V 在:1tut 是 Lyapunov 函数。由25中推论 3.3,可得()有()0V=。假设()()()()()()T,tuS t V tI tA tR t=是系统(1)过任意()的解。由于()的不变性,我们可以得出结论,对于所有的t,()tu。所以得到,对任意的t,()0tV u=。因此,从(9)可知,()*S tS=()()*A t IA I t=,对于所有的t。因此:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 具有 疫苗 接种 COVID 19 传播 模型 动力学 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。