具有随机扰动的SIV传染病模型的动力学行为分析
- 格式:pdf
- 大小:1.14 MB
- 文档页数:5
媒体报道下随机SIQS流行病模型的动态行为研究
近期,一份最新的研究成果表明,SIQS流行病模型对传染病疫情的预测性和准确性更高。
此模型模拟了人们在患病前后实施的不同干预措施,包括隔离、检测、随机筛查等。
这份研究是由美国哈佛大学流行病学家Sarah Cobey等人发布在《科学》杂志上的。
研究团队通过对不同干预措施的模拟,分析不同干预措施对于传染病疫情的影响,以及各种策略所需要的资源和成本。
其中,随机筛查是该模型的最大亮点之一。
随机筛查意味着对未经过筛查的个体进行检测,无论其是否表现出患病症状。
这一策略常常被用于疾病控制的实践中,以遏制疫情进一步扩散。
该模型结果表明,在随机筛查与其他干预措施联合使用的情况下,病毒复制率将迅速下降。
此外,研究团队还发现,隔离与随机筛查相结合的策略,将完全控制疾病疫情,并几乎没有任何经济成本。
相较于传统流行病模型,SIQS模型更加接近真实的疫情传播情况,可以充分利用现代计算机大数据处理的优势,模拟更加复杂和多样化的情况,帮助政府和社会机构制定更加科学有效的疫情应对方案。
然而研究团队也强调,不同的疾病、不同的地区和不同的人群特征可能需要不同的控制策略,我们需要继续推进流行病学研究,探索出更加精准有效的控制策略。
总之,SIQS流行病模型的出现,使得我们在抗击疫情的过程中多了一种新的策略和参考,也为我们提供了更加准确、细致和科学的数据支撑。
㊀㊀㊀㊀㊀102㊀传染病SIR模型在多层网络上动力学行为分析传染病SIR模型在多层网络上动力学行为分析Һ曹㊀蓉1㊀唐㊀甜2㊀童细心∗1㊀(1.汕头职业技术学院自然科学系,广东㊀汕头㊀515041㊀㊀㊀㊀㊀㊀㊀㊀2.桂林电子科技大学电子电路国家级实验教学示范中心,广西㊀桂林㊀541004)㊀㊀ʌ摘要ɔ许多现实世界中的网络都是相互关联和依赖的,如水陆空组成的交通网络和不同群体构成的社会网络等,由此形成具有复杂的拓扑结构和动力学性态的多层次网络,但此类网络上的传染病动力学还缺乏相关研究结果.根据平均场理论提出了一个多层耦合网络上的传染病SIR模型,以刻画传染病在多个种群或社区之间的传播过程.首先给出模型的基本再生数和全局动力学行为,接着分析耦合结构和感染力对传播阈值和疫情的影响.结果表明发病率的相变取决于有网络拓扑结构合传播参数构成的临界值,而相互关联结构更容易造成疫情扩散,且内部接触比交叉接触更容易引起疾病暴发.ʌ关键词ɔ多层网络;传染病平均场模型;动力学行为;基本再生数ʌ基金项目ɔ广西自然科学基金(2017A030313699),汕头职业技术学院2016年院级科研课题(SZK2016Y15)引㊀言随着小世界和无标度性质在众多现实系统和结构中的发现,复杂网络近十多年逐渐变成了一个强大的科学工具,用以描述各种社会㊁经济㊁生物等系统的拓扑结构,能更好地拟合异质结构,并更客观地刻画动力学和结构的关联特点.因现实网络并不是单独存在,而是相互依赖,且一个大系统往往蕴含多个层带不同拓扑和功能的网络,由此形成了一类新的复杂系统,即多层耦合网络.这类网络广泛存在于大型系统,如物流网,交通网和社团等,有着广泛的应用背景和研究价值,比如分析路面的交通载量和电网的停电事件.探索多层耦合网络的结构特点和动力学性质是近十年的研究热点,此类网络正在不断地被用于刻画传染病在多种种群如人和动物之间的传播,并获得了一些新的进展,Allard等人利用键渗流理论模拟了多层网络上的传染病扩散过程;Funk等人根据物理机制建立了一个叠加网络上的疾病传播的建模方案,Dickison等人通过耦合网络上的传染病SIR模型,发现强耦合会使疾病扩散至整个网络,而弱耦合会出现一个混合相图,使疾病只能在一个子网络上传播.Saumell⁃Mendiola等人通过多层网络上的平均场模型,子网耦合能使疾病更容易暴发.传染病动力学模型也越来越多地被用于研究网络传播性质,但多层结构和接触模式对传播动力学的影响还缺乏系统分析.本文利用平均场近似在两层网络相互作用的网络上建立一个传染病SIR模型,利用微分方程理论和数值模拟分析模型的动力学性态,建立传染病动力学和层次网络结构的关系.主要回答2个问题:什么条件会引起传染病在层次网络上的暴发?子网内部和之间的耦合方式和节点的接触模型如何影响着传染病的传播和扩散?上述问题的答案有助于更好地认识传染病在多群体之间的传播规律.一㊁数学模型首先根据子系统的耦合关系,给出一个由两个相互依赖的子网络A和B构成的复杂网络,每个子网由一个群体(节点)以及它们的连接(边)构成.不同群体的接触模式的差异性对应两种不同的连接方式,一种表示同一子网内的个体连接(内部连边),另一种表示不同子网之间的个体连接(交叉连边),故每个节点对应两个度.因为群体内部和交叉的连接的异质性和多样性,所以整个网络具有异常复杂的拓扑结构.此网络可以表示多种现实系统,如性接触网络,其中A,B分别表示男性和女性的性接触网络,子网内部连接表示同性接触,交叉连接表示异性接触,如果节点是双性恋,则同时存在内部和交叉连接;也可表示媒介和宿主的接触模式,A表示宿主如人,B表示媒介如动物,交叉连接表示宿主和媒介之间的接触.下面用参数表示具体的网络结构.度(i,j)表示有i条边与子网A连,有j条边与子网B连.度(i,㊃)表示有i条边与A连,度(㊃,j)表示有j条边与B连.NXi,j表示子网X上度为(i,j)的节点数.SXi,j,IXi,j和RXi,j分别表示网X上度为(i,j)的易感㊁染病和康复的节点数量.PX(i,j)表示任取一个X子网的节点其度为(i,j)的概率,即度分布. k⓪11和 k⓪12分别表示子网A连接子网A和B内节点的平均度. k⓪21和 k⓪22分别表示子网B连接子网A和B内节点的平均度.根据以上定义,子网X上所有的易感㊁染病和康复的节点数分别是SX=ðiðjSXi,j,IX=ðiðjIXi,j,RX=ðiðjRXi,j,子网X的所有节点数是NX=SX+IX+RX.假设网络是度不相关的,即任何一个节点的连接度与邻居的连接度无关,则对于任意的i和j子网X的联合度分布分别为PX(i,j)=NXi,jNX,子网X边界度分布为PX(i,㊃)=ðjPX(i,j)和PX(㊃,j)=ðiPX(i,j),以上X可取为子网A和B.另外,两个子网的平均度(φ=1)和度的二阶矩(φ=2)为kφ⓪11=ðiφPA(i,㊃), kφ⓪12=ðjφPA(㊃,j),kφ⓪21=ðiφPB(i,㊃), kφ⓪22=ðjφPB(㊃,j),忽略个体的出生和死亡,则子网络的节点总数NA和NB保持不变.在子网间的交叉连接中,A连接B的边数等于B连接A的边数,由此可得NA k⓪12=NB k⓪21.这表明如果A连接B的平均度比B连接A的平均度大,则A的节点数大于B㊀㊀㊀103㊀㊀的节点数,这意味着异质网络上交叉连接不均匀的子网络拥有更多的节点数.下面利用仓室SIR模型拟合疾病在此网络上的传播,把群体按状态划分成三个仓室:染病者(S)㊁易感者(I)和康复者(R).在每一个时间步,各个节点处于这三种状态之一.传播过程如下:子网A(B)的一个易感节点通过连边与子网A和B的染病节点连接,分别以概率λ11和λ21(λ12和λ22)被感染变成染病者,经过平均染病期1μ11μ2()后康复.此动力学的时间演化过程可用以下高维的微分方程系统模拟:dIAg,h(t)dt=λ11gSAg,h(t)ðii(ðjIAi,j(t))ðii(ðjNAi,j)+λ21hSAg,h(t)ðii(ðjIBi,j(t))ðii(ðjNBi,j)-μ1IAg,h(t),dRAg,h(t)dt=μ1IAg,h(t),dIBk,l(t)dt=λ12kSBk,l(t)ðjj(ðiIAi,j(t))ðjj(ðiNAi,j)+λ22lSBk,l(t)ðjj(ðiIBi,j(t))ðjj(ðiNBi,j)-μ2IBk,l(t),dRAg,h(t)dt=μ2IAg,h(t),ìîíïïïïïïïïïïïïï(1)以上表示的是各个状态的人群数量变化的动力学模型,设sXg,h=SXg,hNXg,h,ρXg,h=RXg,hNXg,h和rXg,h=IXg,hNXg,h,表示对应的密度,则模型可简化成等价的密度模型:dρAg,h(t)dt=(λ11g11(t)+λ21h21(t))(1-ρAg,h(t)-rAg,h)-μ1ρAg,h(t),dρAg,h(t)dt=μ1ρAg,h(t),dρBk,l(t)dt=(λ12k12(t)+λ22l22(t))(1-ρBk,l(t)-rBk,l)-μ2ρBk,l(t),dρBg,h(t)dt=μ2ρBg,h(t),ìîíïïïïïïïïïï(2)模型(1)和(2)刻画了两层耦合网络上的传染病SIR传播模型.若只存在一个子网络A,那么模型(1)就变成经典网络上的SIR模型;另外,若不存在内部连接,则网络Π变成了一个二分图;如果λ22=0,即子网B内部的接触不传染疾病,那么模型也可以表示某些媒介宿主传染病如登革热,其中媒介(如蚊子)之间不能相互感染.二、数学分析首先我们分析传染病模型中的重要的参数,即基本再生数,其表示一个病人在其染病期平均能够感染的人数.基本再生数作为传播阈值,能够量化传染病潜在的传播能力.同质网络上的基本再生数可表示为R0=cβD,其中c是接触率,β是感染概率,D是平均感染期,而在异质网络上c表示网络的异构性c= k2⓪k⓪.大规模且有无标度结构的网络会使得c很大,由此导致R0>1,此时疾病将不可控.本文建立的模型忽略了种群的出生和死亡,根据SIR模型的特点,最后疾病会灭绝,只剩下易感节点和康复节点,但在疾病消亡之前,不同的网络结构和模型参数会决定是否有过大暴发,其条件即为基本再生数.下面利用下一代生成矩阵,求解模型的基本再生数,经过一系列的等价变换,下一代生成矩阵可写成Γ=λ11ði,ji2PA(i,j)μ1 k⓪11λ11ði,jijPA(i,j)μ1 k⓪110000λ21ði,ji2PB(i,j)μ2 k⓪21λ21ði,jijPB(i,j)μ2 k⓪21λ12ði,jijPA(i,j)μ1 k⓪12λ12ði,jj2PA(i,j)μ1 k⓪120000λ22ði,jijPB(i,j)μ2 k⓪22λ22ði,jj2PB(i,j)μ2 k⓪22æèççççççççççöø÷÷÷÷÷÷÷÷÷÷㊀㊀㊀㊀㊀104㊀㊀㊀当网络是度不相关时,矩阵可进一步简化为Γ=λ11 k2⓪11μ1 k⓪11λ11μ1 k⓪120000λ21 k2⓪21μ2 k⓪21λ21μ2k⓪22λ12μ1 k⓪11λ12 k2⓪12μ1 k⓪120000λ22μ2 k⓪21λ22 k2⓪22μ2 k⓪22æèççççççççççöø÷÷÷÷÷÷÷÷÷÷设ρ(Γ)为矩阵Γ的特征值,则模型的基本再生数为R0=ρ(Γ).由此得主要的理论结果:假设模型初始时刻的暴发规模不大,则当R0<1时,发病规模将会持续降低而至于灭绝;而若R0>1,则传染病会经历一个发病高峰期后慢慢消亡.基本再生数R0得不到其具体表达式,但由Perron⁃Frobenius定理可知,R0小于矩阵Γ每一行之和的最大值,并大于每一行之和的最小值.这意味着,基本再生数受到每个子网络的异质性参数c及平均度的控制,因此异质性强和平均度大的结构都会使得多层网络容易传播疾病.三㊁数值分析下面用数值模拟进一步探索传播动力学.对于动物传染病,宿主的接触模式通常具有异质性,而动物的交互作用具有同质性,对应了两种接触网络,随机网络和无标度网络.随机网络类似于均匀网络,其度分布满足泊松特性,而无标度网络是非常不均匀的,其度分布满足无标度特性.设A(B)是子网A(B)的内部接触模式,AB(BA)是子网A连接B(B连接A)的交叉接触模式.所有的图中子网A和B具有相同的节点数.基本再生数R0和总的感染人数ρA和ρB能反映流行病学的重要性质,图1和图2表明了感染概率和接触模式对它们的影响.由图1和图2可知,在任何的网络结构下感染率的增大都会使得R0和感染规模增大,而无标度结构会使其增加更快.若所有的子网具有相同的结构,则内部感染和交叉感染对R0的作用相等,但若增加内部感染率(λ11或λ22),会导致更大的R0和暴发规模,而交叉感染病(λ12和λ21)则作用不明显,主要原因是内部感染率有双重反应,感染同一子网的其他节点也可以被其他节点感染,而交叉感染率只有单个方向的作用.图1㊀感染率对基本再生数的影响,其中λ11=λ22(内部感染),λ12=λ21(交叉感染)图2㊀感染率对染病规模的影响ρ=ρA+ρB2四㊁结㊀论本文建立了一个两层耦合网络上的传播模型,并求解了模型的基本再生数,给出了疫情暴发的条件,并揭示了基本再生数与网络结构合传播参数的关系.结果表明传播阈值受网络的异质性和平均度控制,这与单个网络和二分图网络只受异质性影响不同,说明层次结构对传播动力学的特殊性.另外多层网络具有多个子结构,其中最异质的那个子结构网络对传播起决定作用.特别地,我们发现内部接触比交叉接触更容易导致疾病传播,这也说明了同性恋在性病传播中的关键作用.ʌ参考文献ɔ[1]汪小帆,李翔,陈关荣.网络科学导论[M].北京:高等教育出版社,2012.[2]郑国庆,唐清干,祝光湖.带接种免疫的网络传染病的有效度模型[J].数学的实践与认识,2015(15):315-322.[3]马知恩,周义仓,王稳地,等.传染病动力学的数学建模与研究[M].北京:科学出版社,2004.。
一类随机SEIQR传染病模型的动力学行为分析李雪;胡良剑【摘要】考虑疾病传播过程中随机因素的影响,研究一类具有潜伏期和隔离仓室的SIR传染病模型(SEIQR).通过构造Lyapunov函数并利用It^o公式,证明了随机SEIQR传染病模型存在唯一的全局正解,得到其关于相应的确定性模型的无病平衡点和地方病平衡点渐近稳定的充分性条件.【期刊名称】《纺织高校基础科学学报》【年(卷),期】2019(032)001【总页数】7页(P37-43)【关键词】随机传染病模型;It(^o)公式;Lyapunov函数;平衡点;渐近稳定【作者】李雪;胡良剑【作者单位】东华大学理学院,上海 201620;东华大学理学院,上海 201620【正文语种】中文【中图分类】O211.630 引言传染病一直以来是医学界普遍关注的重要问题,为了更好地预防和治疗传染病,学者们通过建立数学模型对传染病模型进行定性和定量研究。
最经典的传染病模型就是由Kermack和Mckendrick于1927年创立的SIR(susceptible-infective-removed)仓室模型[1]。
自此,很多传染病模型在此基础之上建立起来[2-3]。
然而,传染病的传播会受到外部环境的随机干扰。
研究在随机因素影响下的传染病模型所具有的性质对疾病的控制有重要意义。
近些年来,国内外学者对于随机传染病模型在生物,医药等领域的应用和研究取得了很多进展[4-6]。
考虑疾病潜伏期的SEIR模型加入随机扰动,已有大量研究[7-10]。
YANG等[11]考虑了随机SIR及SEIR传染病模型解的全局正性,并得到了求解关于相应的确定性系统的无病平衡点及地方病平衡点的渐近行为。
魏凤英等[12]进一步研究了一类具有非线性发病率的随机SEIR传染病模型的灭绝性及平稳分布等问题。
另一方面,也有部分学者考虑将具有隔离项的传染病模型加入随机扰动,ZHANG等[13]对一类非线性发病率的确定性和随机SIQS模型进行研究,给出了模型解是否稳定的阀值R0。
一类离散SIS传染病模型的动力学性态分析一类离散SIS传染病模型的动力学性态分析摘要:传染病是当前全球范围内严重威胁人类健康的重大公共卫生问题之一。
为了掌握传染病的传播规律,许多学者对传染病的传播进行了建模和分析。
本文基于离散SIS传染病模型,对该模型的动力学性态进行深入探讨。
研究结果表明,传染病的传播受到许多因素的影响,其中包括感染率、康复率和出生率。
关键词:离散SIS模型;传染病;动力学性态;感染率;康复率;出生率一、引言传染病是一种由病原体引起的在个体之间传播的疾病。
通过对传染病的传播规律进行建模和分析,可以更好地了解疾病的传播过程,以便制定相应的防控策略。
目前,传染病模型主要分为连续模型和离散模型两种。
离散模型对传染病的传播进行了离散化处理,适用于描述个体之间的接触和传播。
二、模型建立本文基于离散SIS传染病模型,假设人口总数为N,感染者的数目为I,易感者的数目为S。
模型的主要参数有感染率β、康复率γ和出生率μ。
感染率β表示一个感染者每天可以感染多少个易感者。
康复率γ表示一个感染者每天以一定的几率康复。
出生率μ表示一个单位时间内新生的人口增长率。
三、动力学分析1. 平衡点分析通过求解模型方程的平衡状态,可以得到模型的平衡点。
当感染者和易感者的数量不再发生变化时,称为一个平衡点。
平衡点对应着传染病的稳态,有利于分析传染病的传播趋势。
2. 稳定性分析通过线性化稳定性理论,可以判断平衡点的稳定性。
当平衡点是稳定的时,传染病的传播将会趋于平稳;当平衡点是不稳定的时,传染病的传播将会呈现周期性或者混沌行为。
3. 条件分析在模型分析中,除了考虑感染率、康复率和出生率等参数外,还需要对其他条件进行分析。
比如,感染者康复后是否具有免疫力、易感者的接触频率等。
这些条件对传染病的传播规律有重要影响。
四、数值模拟通过数值计算,我们可以更加详细地研究传染病模型的动力学行为。
通过改变感染率、康复率和出生率等参数,我们可以观察到传染病模型的不同状态,进一步验证模型的合理性。
两类传染病模型的动力学分析的开题报告
题目:两类传染病模型的动力学分析
背景和意义:
传染病是全球公共卫生问题中的重要因素,对人类健康和社会经济
持续稳定发展造成威胁。
为了更好地理解传染病的传播规律和控制传染
病的策略,建立和分析传染病模型已成为传染病研究的重要手段之一。
传染病模型可以分为两类:基于微观个体的代谢模型和基于宏观总体的
群体模型。
宏观群体模型常见的有SIR模型、SIS模型和SI模型等,这
些模型是基于群体总体的认知和随机性来分析传染病流行趋势的。
研究方法:
本研究将主要基于建立SIR模型和SIS模型,对两种不同的传染病
进行动力学分析。
针对不同的传染病,将确定它们在人群中的传染性、
潜伏期、传播途径等参数,搭建相应的数理模型,进而进行模拟和优化
分析。
预期结果:
通过对SIR模型和SIS模型的分析,本研究将获得本类模型的动力
学特征,如平衡状态稳定性等基本性质,从而为疫情防控提供科学、有
效的策略。
同时,针对上述两种模型,通过设定不同的参数,分析传染
病的基本再生数(R0)对疫情流行趋势的影响,从而为医学决策提供参
考依据。
研究意义:
本研究将为科学控制传染病的传播和防控提供探索,对传染病的预
防和控制具有重要意义。
此外,本研究将为公共卫生领域提供科学数据,有助于加强对传染病的认识和管理。
第28卷 第3期2023年6月 哈尔滨理工大学学报JOURNALOFHARBINUNIVERSITYOFSCIENCEANDTECHNOLOGY Vol.28No.3Jun.2023
一类随机SICR丙肝模型的动力学分析
康玉娇, 张太雷, 马怡婷(长安大学理学院,西安710064)
摘 要:为了研究现实中不确定因素对丙肝的影响,考虑了一类具有急性与慢性感染的随机SICR丙肝模型。首先,依据停时理论证明了全局正解的存在唯一性;其次,通过构造Lyapunov函数并结合伊藤公式,讨论了随机模型的解在确定模型的无病平衡点和正平衡点附近的振荡行为,并给出了随机模型解的平均持续和疾病灭绝的充分条件。最后,数值模拟考虑了噪声对模型的影响,结果表明:随机模型在确定模型的平衡点附近扰动,扰动强度与噪声强度成正相关,并且足够大的噪声使疾病灭绝,进一步证实了理论结果。关键词:随机模型;平衡点;振荡行为;伊藤公式;疾病持续与灭绝DOI:10.15938/j.jhust.2023.03.016
中图分类号:O175.1文献标志码:A文章编号:1007-2683(2023)03-0129-11
DynamicAnalysisofanSICRHepatitisCInfectionModelwithStochasticEffects
KANGYujiao, ZHANGTailei, MAYiting(SchoolofScience,Chang′anUniversity,Xi′an710064,China)
Abstract:InordertostudytheinfluenceofuncertainfactorsonhepatitisCinreality,arandomizedSICRhepatitisCmodelwithacuteandchronicinfectionwasconsidered.Firstly,theexistenceanduniquenessofglobalpositivesolutionsareprovedaccordingtothestoppingtimetheory;secondly,byconstructingLyapunovfunctionandcombiningItoformula,theoscillationbehaviorofthesolutionofthestochasticmodelnearthedisease-freeequilibriumpointandthepositiveequilibriumpointofthemodelisdiscussed,andthenthesufficientconditionsfortheaveragepersistenceanddiseaseextinctionofthesolutionofthestochasticmodelaregiven.Finally,thenumericalsimulationconsiderstheinfluenceofnoiseonthemodel.Theresultsshowthatthestochasticmodelisdisturbedneartheequilibriumpointofthedeterminedmodel,andthedisturbanceintensityispositivelycorrelatedwiththenoiseintensity.Moreover,thenoiseislargeenoughtomakethediseaseextinct,whichfurtherconfirmsthetheoreticalresults.Keywords:stochasticmodel;equilibriumpoint;oscillatorybehavior;Itoformula;diseasepersistenceandextinction
一个带有隔离和扰动的SIQ模型的动力学分析丁金梅【摘要】对包含隔离策略的传染病模型,通过修改传染率形式、加入随机扰动项等方面的改进建立了更加合理、更广泛适用的随机传染病模型.首先,利用It?公式、大数定律等结论证明了该随机模型正解的几乎必然全局存在性;其次利用It?公式、局部鞅不等式、Hasminskii定理等通过构造二阶连续可微函数得到了疾病指数灭绝的充分条件;最后利用Matlab软件给出仿真结果验证了所得结论.【期刊名称】《宁夏工程技术》【年(卷),期】2017(016)004【总页数】4页(P327-330)【关键词】传染病模型;隔离;随机扰动;灭绝性【作者】丁金梅【作者单位】宁夏大学人事处,宁夏银川 750021【正文语种】中文【中图分类】O175.13自Kermack和Mckendrick提出著名的传染病模型[1],利用微分方程相关理论研究传染病的传播规律、预测其传播趋势成为一个重要的研究热点。
许多学者对多种形式的传染病模型进行了深入的研究[2—6],特别是考虑年龄结构、时滞、扰动等因素的传染病模型更加符合实际情况,越来越得到大家的关注[7—10]。
隔离和免疫策略是抑制疾病扩散和爆发的两种重要手段。
Herbert和马知恩等[11]考虑了隔离策略对传染病传播的影响,就传染率具有的不同形式,讨论了6种传染病模型,得到了每个模型平衡点具有全局渐近稳定的条件。
其中一类具有双线性传染率的模型为式中:S(t),I(t),Q(t)分别代表种群中易感者、感染者和隔离者的数量。
模型中所有系数均为正常数,式中,Λ为再生率,指种群的移入者和新出生者的速率,假设为常数;β为传染率系数;μ为种群的自然死亡率;γ,ε分别为感染者和隔离者恢复为易感染者的比例;δ为感染者被隔离的比例。
考虑到某些传染病,如艾滋病、乙型肝炎等传播方式包含垂直传播,即可通过母婴直接传播;另对于SARS、MERS、登革热等重大疫情爆发时,各国(或地区)均会对移入者进行严格审查,对疑似者直接隔离。
For personal use only in study and research; not for commercial use传染病模型详解2.2.2 /,SI SIS SIR 经典模型经典的传播模型大致将人群分为传播态S ,易感染态I 和免疫态R 。
S 态表示该个体带有病毒或谣言的传播能力,一旦接触到易感染个体就会以一定概率导致对方成为传播态。
I 表示该个体没有接触过病毒或谣言,容易被传播态个体感染。
R 表示当经过一个或多个感染周期后,该个体永远不再被感染。
SI 模型考虑了最简单的情况,即一个个体被感染,就永远成为感染态,向周围邻居不断传 播病毒或谣言等。
假设个体接触感染的概率为β,总人数为 N ,在各状态均匀混合网络中建立传播模型如下:dS SI dt N I SId tN ββ⎧=-⎪⎪⎨⎪=⎪⎩ 从而得到(1)di i i dtβ=- 对此方程进行求解可得:0000(),01tti e i t i i i i e ββ==-+() 可见,起初绝大部分的个体为I 态,任何一个S 态个体都会遇到I 态个体并且传染给对方,网络中的S 态个数随时间成指数增长。
与此同时,随着I 态个体的减少,网络中S 态个 数达到饱和,逐渐网络中个体全部成为S 态。
然而在现实世界中,个体不可能一直都处于传播态。
有些节点会因为传播的能力和意愿 的下降,从而自动转变为永不传播的R 态。
而有些节点可能会从S 态转变I 态,因此简单的SI 模型就不能满足节点具有自愈能力的现实需求,因而出现SIS 模型和SIR 模型。
SIR 是研究复杂网络谣言传播的经典的模型。
采用与病毒传播相似的过程中的S ,I ,R 态 代表传播过程中的三种状态。
Zanetee ,Moreno 先后研究了小世界传播过程中的谣言传播。
Moreno 等人将人群分为S (传播谣言)、I (没有听到谣言),R (对谣言不再相信也不传播)。
假设没有听到谣言I 个体与S 个体接触,以概率()k λ变为S 个体,S 个体遇到S 个体 或R 个体以概率()k α变为R ,如图 2.9 所示。
一类具有传播媒介的SIS传染病模型的动力学性态分析
代林烽;陈龙伟
【期刊名称】《理论数学》
【年(卷),期】2023(13)2
【摘要】随着传媒技术的进步和成熟,媒体对疾病疫情的报道产生的影响已经成为人们预防和控制疾病传播不可忽略的因素之一。
文中考虑了疾病信息意识影响,建立了一个新的具有媒介影响的SIS传染病模型。
首先,对建立的具有媒介影响的新的SIS传染病模型讨论了无病平衡点和地方病平衡点的存在性和基本再生数。
其次,应用Routh-Hurwitz判别准则分别证明了无病平衡点和地方病平衡点的局部渐近稳定性的满足条件。
紧接着构造了恰当Lyapunov函数并结合LaSalle不变原理,分别讨论了无病平衡点和地方病平衡点在不同阈值条件下的全局渐近稳定性。
最后进行了数值模拟,验证在不同阈值条件下系统的稳定性,比较参数在不同取值时平衡点的变化,分析传染病的发展规律。
数值模拟的结果与理论分析的结果是吻合的。
【总页数】13页(P294-306)
【作者】代林烽;陈龙伟
【作者单位】云南财经大学昆明
【正文语种】中文
【中图分类】O17
【相关文献】
1.一类具有免疫反应和抗逆转录病毒治疗的HIV病毒传染病模型的动力学性态分析
2.一类具有非线性发生率的随机SIS传染病模型阈值动力学行为研究
3.一类具有分布时滞和饱和发生率的海洛因传染病模型的全局动力学性态
4.一类具有非线性传染率和生育脉冲的随机SIS传染病模型的动力学分析
5.具有饱和治愈率的离散SIS传染病模型的动力学性态
因版权原因,仅展示原文概要,查看原文内容请购买。
传染性疾病的传播动力学与模型研究随着全球化进程的加速和人口流动的增加,传染性疾病的爆发和蔓延成为了一项全球性的挑战。
了解传染病传播的动力学以及建立相应的数学模型是预测、防控和治疗传染病的重要手段。
本文将探讨传染性疾病的传播动力学和模型研究。
一、传染性疾病的传播动力学传染性疾病的传播动力学研究主要关注疾病的传播途径、传播速度以及影响传播的因素。
首先,传染病的传播途径可以是经空气、直接接触、水源和食物等途径。
其次,传播速度取决于传染病的传染性和人群的接触频率。
最后,影响传播的因素包括人群的免疫力、卫生条件、医疗水平等。
二、传染性疾病的数学模型数学模型在研究传染性疾病传播动力学中起着至关重要的作用。
常用的数学模型有SIR模型、SEIR模型、SI模型等等。
SIR模型将人群分为易感者(Susceptible)、感染者(Infected)和康复者(Recovered),并通过一组微分方程来描述这三个群体之间的转变。
SEIR模型在SIR模型的基础上加入了潜伏期(Exposed)这一状态,更准确地描述了疾病传播过程。
SI模型则忽略了康复者,并假设感染者永远不康复。
三、数学模型的应用数学模型广泛应用于预测传染病的传播趋势、制定防控策略以及评估防控效果。
通过模拟不同因素的变化,可以得出疾病的传播规律和影响传播的关键因素。
在实际应用中,政府和卫生部门可以利用这些模型预测疫情走势,制定针对性的防控策略,减少疫情蔓延。
此外,通过模型的结果还可以评估不同的防控措施对疫情的影响,为防控工作提供科学依据。
四、模型研究的挑战和展望尽管数学模型在传染病研究中取得了一定的成果,但仍然面临着一些挑战。
首先,人群行为的不确定性使得模型的参数估计和预测存在一定的误差。
其次,缺乏准确的疾病数据和人口数据也限制了模型的精确性。
未来,我们需要加强数据收集和分析能力,提高模型的准确性和预测能力。
总之,传染性疾病的传播动力学与模型研究是预测、防控和治疗传染病的重要手段。
暋暋文章编号:1673灢5196(2018)05灢0150灢05具有随机扰动的SIV传染病模型的动力学行为分析黄灿云,郝一新,孟新友(兰州理工大学理学院,甘肃兰州暋730050)
摘要:研究具有随机扰动的SIV传染病模型的动力学行为.通过构造恰当的李雅普诺夫函数并运用伊藤公式,证
明了系统全局正解依概率1的存在性,给出了系统随机最终有界和随机持久的充分条件.利用Higham等提出的Milstein方法对所给的系统进行了数值模拟.
关键词:SIV传染病模型;随机扰动;全局正解;随机持久
中图分类号:O175暋暋文献标志码:A
DynamicbehavioranalysisofstochasticSIVepidemicmodelHUANGCan灢yun,HAOYi灢xin,MENGXin灢you(SchoolofScience,LanzhouUniv.ofTech.,Lanzhou暋730050,China)Abstract:DynamicbehaviorofSIVepidemicmodelwithstochasticperturbationswasinvestigated.By
constructingsuitableLyapunovfunctionsandapplyingIt^oformula,theexistenceofglobalpositivesolu灢
tionswithprobabilitytothesystemwasprovedandthesufficientconditionofstochasticultimatebounded灢
nessandstochasticpermanenceofthesystemwasgiven.EmployingMilsteinmethodpresentedby
Highametalii,thenumericalsimulationofthegivensystemwasperformed.Keywords:SIVepidemicmodel;stochasticperturbations;globalpositivesolution;stochasticperma灢
nence
暋暋传染病动力学的研究主要起源于20世纪初期,受经典的SIR模型研究的影响,大部分研究结果主要分析传染病的传播规律[1灢3].随着研究水平的提高,近年来,大部分学者开始研究疫苗接种对传染病预防和控制的影响,希望为传染病的防治决策提供理论基础和数量依据[4灢7].Kribs灢zaleta等[7]研究了具有疫苗接种的SIV模型,其研究结果表明疫苗的接种成功率对传染病的控制能起到决定性的作用.目前,为了刻画传染病的流行规律,大多数学者基于确定性动力学模型进行研究,主要依据是假设生物种群的个体数量足够大.根据大数定理,系统的状态会呈现出比较平稳的统计规律,从而可以把系统近似地看成确定性的,这就使得研究工作更加简便.事实上,这些确定性模型的研究成果的确对传染病的防治工作起到了积极作用.然而,在现实生活中,任何一个生物系统都不可暋暋收稿日期:2017灢03灢14暋暋基金项目:国家自然科学基金(11661050)暋暋作者简介:黄灿云(1965灢),男,河南淮阳人,博士,教授.能孤立地存在,而是处于一个复杂的相互联系的环
境中,各种形式的随机干扰无处不在.环境噪声会不同程度地影响到增长率、环境容量和系统的其它参数,而且实际资料也表明系统状态的随机波动非常明显,因此在某些情况下使用随机数学模型对传染病进行描述尤其必要[8灢10].受上述文献的启发,本文将随机扰动引入到Kribs灢zaleta等[7]所研究的SIV
确定性系统中,分析环境噪声对SIV模型动力学行为的影响.假设随机扰动是白噪声类型,且与S(t),
I(t),V(t)到稳态值S*,I*,V*
的距离以及对
dS(t)/dt,dI(t)/dt,dV(t)/dt的影响成正比,那么
具有随机扰动的SIV传染病模型为dS=(d-毬SI-(d+p)S+毭I+eV)dt+
暋暋氁1(S-S*)dB(t)
dI=(毬SI-(毭+d)I+氁毬VI)dt+
暋暋氁2(I-I*)dB(t)
dV=(pS-氁毬VI-(d+e)V)dt+
暋暋氁3(V-V*)dB(tìîíïïïïïïïï)(1)
第44卷第5期2018年10月兰暋州暋理暋工暋大暋学暋学暋报JournalofLanzhouUniversityofTechnologyVol.44No.5Oct.2018式中:B(t)是标准布朗运动;S(t)代表t时刻易感者数量;I(t)代表t时刻染病者数量;V(t)代表t时刻接种者数量;d是人口自然出生率和死亡率;毬是传染率;毭是染病者的恢复率;p代表接种率;e是疫苗消失的比率;氁是疫苗对染病者减少率的影响,氁=0
代表疫苗对阻止感染完全有效,氁=1表示疫苗完全无效,通常令0)是实的常数,氁2
i
表示随机干扰的强度;所有参数均为正值.当氁1
=
氁2=氁3
=0时,系统(1)是确定性的模型,Kribs灢zale灢
ta等[7]对其动力学行为进行了大量的研究.本文主
要分析系统(1)依概率1全局正解的存在性,给出系统随机最终有界和随机持久的充分条件.
1暋动力学行为
定理1暋对任意给定的初值(S(0),I(0),V
(0))暿
R3
+
,系统(1)存在唯一解(S(t),I(t),V(t)),且此解
依概率1停留在R3+中.证明暋由于系统(1)的系数满足局部李普希兹
条件,则系统(1)在[0,氂e)上存在唯一的局部解,其
中氂e是爆炸时间[10].由It^o公式知该局部解是正的,接下来只需证明此解是全局的,也就是说,氂e=
曓a.s..令n0>0充分大,使得S(0),I(0),V(0)都
在区间[1/n0,n0]内.对每一个整数n曒n0
,定义停时
氂n=inf{t暿[0,氂e]:S(t)烖1/n,()n
或
I(t)烖1/n,()n暋或暋V(t)烖1/n,()n}
(2)
令inf熈=曓(熈是空集),由于n曻曓时氂n是非减的,则存在极限氂曓=limn曻曓氂n
(3)
显然,氂曓曑氂ea.s..因此只需证明氂曓=曓a.s..否则,存在T>0和毰暿(0,1),使得P{氂曓曑T}>毰.因此,存在整数n1曒n0使得对任意的n曒n1
,
P{氂n曑T}曒毰(4)
构造函数V(S,I,V)=(S-1-lnS)+(I-1-lnI)+
暋暋暋暋(V-1-lnV),暋(S,I,V)暿R3+(5)
显然V(S,I,V)是非负函数.利用It^o公式可得dV(S,I,V)=[(1-1/S)(d-毬SI-(d+p)S+
毭I+eV)+(1-1/I)(毬SI-(毭+d)I+氁毬VI)+
(1-1/V)(pS-氁毬VI-(d+e)V)+
氁21(1-S*/S)2/2+氁22(1-I*/I)2/2+
氁23(1-V*/V)2/2]dt+[氁1(1-1/S)(S-S*
)
+
氁2(1-1/I)(I-I*)+氁3(1-1/V)(V-V*)]dB(t)=
[(毬I+3d+p+e+毭-毬S-氁毬V-d/S-毭I/S-eV/S-pS/V)+氁21(1-S*/S)2/2+氁22(1-I*I)2/2+氁23(1-V*/V)2/2]dt+[氁1(1-1/S)(S-S*)+氁2(1-1/I)(I-I*)+氁3(1-1/V)(V-V*)]dB(t)曑[(1+氁)毬I+3d+p+e+毭+氁21(1-S*/S)2/2+氁22(1-I*/I)2/2+氁23(1-V*/V)2/2]dt+[氁1(1-1/S)(S-S*)+氁2(1-1/I)(I-I*)+氁3(1-1/V)(V-V*)]dB(t)曑Mdt+[氁1(1-1/S)(S-S*)+氁2(1-1/I)(I-I*)+氁3(1-1/V)(V-V*)]dB(t)(6)式中:M是一个正常数.对上述不等式(6)两端从0到氂n暷T积分可得曇氂n暷T0dV(S(s),I(s),V(s))曑曇氂n暷T0Mds+曇氂n暷T0[氁1(1-1/S)(S-S*)+氁2(1-1/I)(I-I*)+氁3(1-1/V)(V-V*)]dB(s)炤V(S(氂n暷T),I(氂n暷T),V(氂n暷T))曑V(S(0),I(0),V(0))+M(氂n暷T)+曇氂n暷T0[氁1(1-1/S)(S-S*)+氁2(1-1/I)(I-I*)+氁3(1-1/V)(V-V*)]dB(s)(7)对式(7)取均值得E[V(S(氂n暷T),I(氂n暷T),V(氂n暷T))]曑V(S(0),I(0),V(0))+MT令毟n={氂n曑T},n曒n1,从式(4)有P(毟n)曒毰.对任何毻暿毟n,S(氂n,毻),I(氂n,毻),V(氂n,毻)或者等于n,或者等于1/n.因此,V(S(氂n,毻),I(氂n,毻),V(氂n,毻))不小于min{n-1-lnn,1/n-1-ln1/n},则有V(S(0),I(0),V(0))+MT曒E[1毟n(毻)V(S(氂n),I(氂n),V(氂n))]曒毰min{n-1-lnn,1/n-1-ln1/n}其中:1毟n是毟n的指标函数.令n曻曓,则有曓=V(S(0),I(0),V(0))+MT毮}暿
·151·第
5期暋暋暋暋暋暋暋暋暋暋黄灿云等:具有随机扰动的SIV传染病模型的动力学行为分析暋暋暋暋暋暋暋暋暋暋暋暋暋