媒体报道诱导的离散传染病切换系统建模与分析.pdf

D O I:1 0.3 9 6 9/j.i s s n.1 0 0 3-0 9 7 2.2 0 2 3.0 3.0 0 9 文章编号:1 0 0 3-0 9 7 2(2 0 2 3)0 3-0 3 9 6-0 7媒体报道诱导的离散传染病切换系统建模与分析覃文杰1,2*,张嘉敏1,向中义3(1.三峡大学 理学院,湖北 宜昌 4 4 3 0 0 2;2.云南民族大学 数学与计算机科学学院,云南 昆明 6 5 0 5 0 4;3.湖北民族大学 数学与统计学院,湖北 恩施 4 4 5 0 0 0)摘 要:针对疾病数据收集单位不连续的情况,运用E u l e r离散法建立了一类由媒体报道诱导的离散S I R传染病切换模型。从理论上研究了切换系统子系统地方病平衡点的存在性、稳定性等动力学行为。从数值上研究了切换系统的平衡点分支、单参数分支、初值敏感性。研究发现:关键参数的小扰动直接影响着疾病的暴发次数与频率,同时媒体报道可以降低疾病传染率,有效遏制疾病的流行与暴发。关键词:媒体报道;阈值策略;离散模型;切换系统;分支分析中图分类号:O 1 7 5 文献标识码:A开放科学(资源服务)标识码(O S I D):M o d e l i n g a n d A n a l y s i s o f a D i s c r e t e S w i t c h e d E p i d e m i c S y s t e m I n d u c e d b y M e d i a R e p o r t sQ I N W e n j i e1,2*,Z H A N G J i a m i n1,X I A N G Z h o n g y i3(1.C o l l e g e o f S c i e n c e,C h i n a T h r e e G o r g e s U n i v e r s i t y,Y i c h a n g 4 4 3 0 0 2,C h i n a;2.S c h o o l o f M a t h e m a t i c s a n d C o m p u t e r S c i e n c e,Y u n n a n M i n z u U n i v e r s i t y,K u n m i n g 6 5 0 5 0 4,C h i n a;3.S c h o o l o f M a t h e m a t i c s a n d S t a t i s t i c s,H u b e i M i n z u U n i v e r s i t y,E n s h i 4 4 5 0 0 0,C h i n a)A b s t r a c t:I n v i e w o f t h e d i s c o n t i n u i t y o f d i s e a s e d a t a c o l l e c t i o n u n i t s,E u l e r d i s c r e t e m e t h o d i s u s e d t o e s t a b l i s h a k i n d o f d i s c r e t e S I R i n f e c t i o u s d i s e a s e s w i t c h i n g m o d e l i n d u c e d b y m e d i a r e p o r t s.T h e e x i s t e n c e a n d s t a b i l i t y o f t h e e n d e m i c e q u i l i b r i a o f t w o s u b s y s t e m s a r e s t u d i e d t h e o r e t i c a l l y.B y n u m e r i c a l s i m u l a t i o n,t h e e q u i l i b r i a b i f u r c a t i o n,s e n s i t i v i t y a n a l y s i s o f t h e k e y p a r a m e t e r a n d i n i t i a l v a l u e a r e i n v e s t i g a t e d.T h e r e s u l t s s h o w t h a t t h e m i n o r p e r t u r b a t i o n o f s o m e p a r a m e t e r s c a n i n f l u e n c e t h e f r e q u e n c y a n d t i m e s o f t h e o u t b r e a k o f t h e i n f e c t i o u s d i s e a s e.M e a n w h i l e,t h e d i s e a s e t r a n s m i s s i o n r a t e c a n b e r e d u c e d t h r o u g h m e d i a r e p o r t s,a n d t h e o u t b r e a k o f t h e d i s e a s e c a n b e p r e v e n t e d e f f e c t i v e l y.K e y w o r d s:m e d i a c o v e r a g e;t h r e s h o l d s t r a t e g y;d i s c r e t e m o d e l;s w i t c h e d s y s t e m s;b i f u r c a t i o n a n a l y s i s0 引言传染病危害着人类生命健康、影响着社会经济发展,因此人们一直以来都在与传染病作斗争。通过发展传染病动力学模型来研究疾病的传播机制1,可为传染病的预防和控制提供理论和数据支持,具有重要的现实意义。影响传染病传播的因素有很多,近年来人们越来越关注媒体报道对人们行为方式的改变。在疾病暴发初期,媒体会通过网络、电视、报纸等不同的渠道实时报道相关情况,人们便会自觉采取相应的预防措施来阻止和减少疾病的传播,如出门戴口罩、少聚集、常消毒等2。因此国内外学者也持续关注媒体报道对疾病传播的影响。2 0 0 3年,RUAN等3用f(I)=k/(1+a I2)刻画媒体影响,得出基本再生数R0对疾病的暴发或消除起决定性作用。2 0 0 7年,L I U等4通过构造传染率e x p(a1E-a2I-a3H)来研究媒体或心理效应对传染病的影响。2 0 0 8年,C U I等5考虑含有媒体影响因子m的传染率e x p(-m I),建立S E I模型来讨论媒体报道对疾病传播的影响。2 0 1 3年,刘玉英等6构造具有饱和性的媒体影响因子函数来刻画媒体对传染病感染率的影响。其他有关媒体效应的传染病模型可见文献7-1 3。收稿日期:2 0 2 2-0 9-1 6;修订日期:2 0 2 2-1 1-2 5;*.通信联系人,E-m a i l:w e n j i e q i n h o t m a i l.c o m 基金项目:国家自然科学基金项目(1 2 2 6 1 1 0 4,1 2 0 2 6 2 3 3,1 2 0 2 6 2 2 1);云南省“兴滇英才支持计划”青年人才项目(X D Y C-QN R C-2 0 2 2-0 7 0 8)作者简介:覃文杰(1 9 8 3),男,湖北五峰人,副教授,博士,硕士研究生导师,主要从事生物数学研究。693信阳师范学院学报(自然科学版)J o u r n a l o f X i n y a n g N o r m a l U n i v e r s i t y第3 6卷 第3期 2 0 2 3年7月 N a t u r a l S c i e n c e E d i t i o n V o l.3 6 N o.3 J u l.2 0 2 3 然而,目前基于媒体影响的传染病模型主要集中在连续模型,对于媒体影响的离散传染病模型还并不多见。事实上,对于某些传播速度较慢或疾病数据收集单位不连续的情形,用离散模型来刻画传播规律将更为合理。本文正是基于此而建立媒体报道诱导的离散传染病切换模型,研究媒体报道因子对传染病传播的影响。1 模型的建立在K E RMA C K和MC K E N D R I C K1 4的S I R仓室模型基础上,考虑具有媒体报道影响的S I R连续模型:S=b S1-SK -I S I-d S,I=I S I-(d+)I,R=S-d R,(1)式中:S=S(t)、I=I(t)、R=R(t)分别表示t时间内易感者、感染者、恢复者的数量,b是内禀增长率,K是环境容纳量,d是自然死亡率,是因病死亡率,是感染者的自然恢复率,(I)表示疾病在媒体报道影响下的传染率。基于文献3-6,考虑媒体影响因子m,以及疾病达到暴发点时感染者数量的临界值Ic,将其构造为非线性饱和函数:(I)=1+m I,=0,IIc,1,IIc,这里所有参数都是非负数。借助E u l e r离散1 5-1 7,将模型(1)离散化,得到Sn+1=Sn+b Sn1-SnK -In SnIn-d Sn,In+1=In+In SnIn-(d+)In,Rn+1=Rn+In-d Rn,(2)由于模型(2)的前两个方程不含Rn,故简化模型(2)可得Sn+1=Sn+b Sn1-SnK -In SnIn-d Sn,In+1=In+In SnIn-(d+)In。(3)通常情形下,当感染人数较少(IIc)时,媒体不会关心此类传染病,(In)=,模型(3)为Sn+1=Sn+b Sn1-SnK -SnIn-d Sn,In+1=In+SnIn-(d+)In。(4)然而随着感染人数的不断增加并达到阈值Ic时,相关媒体会持续关注此类传染病的传播与流行,全媒体的宣传和报道会使公众关注疾病的动态,从而采取相应的防护措施,改变个人行为,此时模型(4)为Sn+1=Sn+b Sn1-SnK -1+m InSnIn-d Sn,In+1=In+1+m InSnIn-(d+)In。(5)综合模型(4)和(5),得到如下切换模型:Sn+1=Sn+b Sn1-SnK -1+m InSnIn-d Sn,In+1=In+1+m InSnIn-(d+)In,(6)式中:=0,IIc,1,IIc。离散模型(6)是一个基于阈值控制策略的动态切换系统,接下来将从理论和数值两方面来研究模型(6)的动力学行为。2 子系统的动力学行为为了叙述方便,将模型(6)改写成向量形式Z(t)=FG1(Z),ZG1,FG2(Z),ZG2,(7)式中:FG1(Z)=Sn+b Sn1-SnK -SnIn-d SnIn+SnIn-(d+)In ,FG2(Z)=Sn+b Sn1-SnK -1+m InSnIn-d SnIn+1+m InSnIn-(d+)In ,G1=ZR2+|I(Z)0,I0,G2=ZR2+|I(Z)0,S0,I0,I(Z)=In-Ic,Z=Sn,InT。特别地,切换系统定义在区域G1和G2的子系统分别记为FG1和FG2,下面将分析研究两个子系统正平衡点的存在性和稳定性。2.1 子系统FG1的动力学行为子系统FG1的平衡点由如下方程决定:b Sn1-SnK -SnIn-d Sn=0,SnIn-(d+)In=0,793覃文杰,张嘉敏,向中义.媒体报道诱导的离散传染病切换系统建模与分析显然,系统FG1具有零平衡点(0,0)、平衡点E0G1和唯一的地方病平衡点E*G1,其中E0G1=b-d Kb,0 ,E*G1=(S*G1,I*G1)=d+,bd+K2R0-1 ,R0=b-d Kbd+,R0为子系统FG1的基本再生数,且有b-d0。记D=b-d,M=d+。定理 1 若子系统FG1满足下列条件之一:1)1R01+1M,0D2,(8)2)D(M+2)DM+4R01+1M,2D41+1M ,(9)则子系统FG1的地方病平衡点E*G1是局部渐近稳定的。证明 子系统FG1在平衡点E*G1的雅可比矩阵为J1=1+D-b MKR0+1 -Mb MKR0-1 1 ,(1 0)由式(1 0)可得t r(J1)=2-DR0,d e t(J1)=1+DM1-1R0 。根据J u r y判据1 8,当|t r(J1)|d e t(J1)+1,且满足d e t(J1)0的基础上对其进行分析:(1)当=B2-4A C0,即R00-2mD2+m2D2时,考虑以下3种情况:若C=1-R00时方程(1 2)存在唯一的正实根,当-B/(2A)0时不存在正实根;若C=1-R00,则当-B/(2A)0时方程(1 2)存在两个正实根,当-B/(2A)1(或00)时,子系统FG2存在地方病平衡点E*G2。定理 2 若子系统FG2满足如下条件:2b M S*G2b(S*G2)2+KDM,2b M2S*G2b(M-1)(S*G2)2+K(D+1)M2-MK S*G2,(1 3)则子系统FG2的平衡点E*G2是局部渐近稳定。证明 由于E*G2=(S*G2,I*G2)满足S*G2=M1+m I*G2 ,I*G21+m I*G2=D-bKS*G2,(1 4)子系统FG2的雅可比矩阵为J2=1-bKS*G2-M2S*G2D-bKS*G21-M+M2S*G2 ,(1 5)893第3 6卷 第3期信阳师范学院学报(自然科学版)h t t p:/j o u r n a l.x y n u.e d u.c n2 0 2 3年7月显然,其特征方程为P()=2-t r(J2)+d e t(J2),(1 6)式中:t r(J2)=2-M+M2S*G2-bKS*G2,d e t(J2)=1-M 1-bKS*G2 +M2D+1S*G2-2bK 。由条件可知,P(1)0,P(-1)0,d e t(J2)1,即特征方程所有根的模都小于1,结合J u r y判据即可得子系统FG2的地方病平衡点E*G2是局部渐近稳定的。证毕。3 切换系统的复杂动力学行为由于离散切换系统各类平衡态存在的类型比较复杂,从理论上研究系统的动力学行为存在很多困难,因此本节将通过数值仿真方法研究传染病切换系统的复杂动力学行为,其中包括平衡点分支、单参数分支、初值敏感性和多吸引子共存等。3.1 切换系统的平衡点分支分析为了研究方便,给出切换系统真假平衡态的定义1 9-2 2:定义 1 若Z*是FG1的平衡点且I(Z*)0(或Z*是FG2的平衡点且I(Z*)0),则点Z*称为切换系统的真平衡态,记为ErG1(或ErG2);若Z*是FG1的平衡点且I(Z*)0(或Z*是FG2的平衡点且I(Z*)0),则点Z*称为切换系统的假平衡态,记为EvG1(或EvG2)。为了研究系统以内禀增长率b和阈值Ic为分支参数的平衡点分支,固定参数K=0.5,=1.8,m=0.2 8,d=0.0 2,=0.3 5,=0.1,如图1所示,当b0.6 5,2.5 0,Ic0.1,1.5 时,离散切换系统存在6种情形的平衡态:区域I-1(绿色为ErG1)和I-2(黄色为EvG1)为子系统FG1的平衡态区域;区域I I-1(红色为ErG1和ErG2共存)为子系统FG1的真平衡态ErG1和子系统FG2的真平衡态ErG2的共存区域;区域I I-2(蓝色为ErG1和EvG2共存)为子系统FG1的真平衡态ErG1和子系统FG2的假平衡态EvG2共存区域;区域I I-3(青蓝色为EvG1和ErG2共存)为子系统FG1的假平衡态EvG1和子系统FG2的真平衡态ErG2共存区域。区域I I-4(洋红色为EvG1和EvG2共存)为两个子系统的假平衡态共存区域;研究结果表明:通过媒体报道策略影响传染病的传播,要将离散切换系统的最终状态控制进入区域I I-4,即两个子系统都不存在真平衡态,换言之,突发传染疾病不会成为地方病,不会影响到人们的正常生产、生活。一定要避免区域I I-1及区域I I-3所示的情形发生,因为传染疾病一旦形成地方病,将直接威胁着人们的身体健康、经济的发展、社会的稳定。图1 系统以b和Ic为分支参数的平衡态分支图F i g.1 E q u i l i b r i a b i f u r c a t i o n d i a g r a m o f s y s t e mf o r p a r a m e t e r s b a n d Ic3.2 单参数分支分析为了进一步研究切换系统的敏感参数、关键参数对系统动力学行为的影响,本小节将分别以自然恢复率、传染率以及因病死亡率为分支参数进行系统的单参数分支分析。首先,固定参数b=2.8,K=1.8,=0.5 2,m=0.8,d=0.0 8,=0.2 5,Ic=1,以自然恢复率作为分支参数来研究切换系统的动力学行为,如图2所示。随着从0.0 5到0.4 0的不断变化,切换系统经历了周期解分支、倍周期解分支、混沌、周期减半、周期窗口、拟周期解等复杂的动力学行为,这也表明了突发性传染的自然恢复率在控制传染病传播方面起着相当重要的作用,在制定相关防控策略时一定要重点关注疾病的自然恢复率。其次,固定参数b=1.8,K=3,m=0.0 2,d=0.3,=0.3,=0.0 2,Ic=5,以疾病的传染率作为分支参数来研究切换系统的动力学行为,如图3所示。当0.6 0,0.8 8 不断变化时,切换系统呈现出稳定解、混沌解、周期倍增解、多重稳定解等复杂的动力学行为。这说明了传染率对疾病的传播有着重要影响,为了控制疾病传播稳定在可控的993覃文杰,张嘉敏,向中义.媒体报道诱导的离散传染病切换系统建模与分析范围内,需要采取措施降低传染率,如增大媒体宣传报道的力度就是间接降低疾病传染率的有效途径之一,这也是本文重点研究的内容之一。图2 系统以为分支参数的分支图F i g.2 B i f u r c a t i o n d i a g r a m o f s y s t e m f o r p a r a m e t e r 图3 系统以为分支参数的分支图F i g.3 B i f u r c a t i o n d i a g r a m o f s y s t e m f o r p a r a m e t e r 接下来,研究因病死亡率对系统动力学行为的影响,以作为分支参数来研究切换系统的动力学行 为。固 定b=1.8,K=2.8,=0.5 6,d=0.0 8,m=0.0 2,=0.0 0 1,Ic=2,如图4所示,随着因病死亡率的不断变化,切换系统本身的动力学行为发生了不同的变化状态。图4 系统以为分支参数的分支图F i g.4 B i f u r c a t i o n d i a g r a m o f s y s t e m f o r p a r a m e t e r 在图4(a)(=0.5)中,系统具有稳定解;当从0.5减小到0.2时,系统出现一个周期解,如图4(b)所示;当继续变小达到0.1或0.0 6时,系统会出现混沌解,如图4(c)和图4(d)所示。3.3 初值敏感性分析除了重要参数、敏感参数对疾病传播有着重要影响外,一定区域内易感者和感染者的初始密度对疾病的暴发状态也会产生不同的影响,本小节将系统研究切换系统在不同初值下的动态行为。固定参数b=1.6,K=2,=1.8 5,m=0.2,d=0.5,=0.8 5,=0.0 9,Ic=0.3 6来研究初值对离散切换系统动力学行为的影响,如图5所示。图5 系统的初始密度对疾病控制的影响F i g.5 T h e i m p a c t o f i n i t i a l d e n s i t i e s o f s y s t e m o n d i s e a s e c o n t r o l图5中所示情形的初值都在切换线以下,即系统的感染者并未达到疾病暴发的阈值。具体而言,图5(a)表示在初值(0.8 2 0 8,0.2 4 6 0)下,切换系统感染者密度未达到疾病暴发的阈值;图5(b)表示在初值(0.9 4 2 3,0.2 0 8 9)下,切换系统感染者密度仅仅一次达到了疾病暴发的阈值,即突发性传染疾病只暴发一次;图5(c)表示在初值(1.0 7 5 0,0.2 3 6 9)下,切换系统感染者密度两次达到了疾病暴发的阈值,即突发性传染疾病暴发两次;图5(d)表示在初值(1.3 0,0.1 8)下,切换系统感染者的密度多次达到了疾病暴发的阈值,即突发性传染疾病多次暴发。研究表明,为了抑制疾病的暴发流行,应该适当控制易感者和感染者的初始数量,让系统稳定在子系统FG1中,即图5(a)、图5(b)的情形。同时,为了控制疾病暴发,需要反复多次地对公众通过媒体宣传疾病传播的相关知识,使得人们对疾病防护引起相当的重视,从而改变人们的个人行为方式,达到控制疾病传播的目的。为了更进一步研究离散切换系统的初始密度对突发性传染疾病传播的影响,下面给出了系统的盆吸引区域图,即选取参数b=0.8,K=3,=1.8 5,m=0.5 7,d=0.5,=0.8 5,=0.2,Ic=1.8,给出系统在区域S00,2.5,I00,2.5的平面分支图,将5个区域分别记为I、I I、I I I、I V004第3 6卷 第3期信阳师范学院学报(自然科学版)h t t p:/j o u r n a l.x y n u.e d u.c n2 0 2 3年7月及V,如图6所示。其中区域I表示疾病不会暴发为地方病;区域I I表示疾病暴发一次后就不再流行;区域I I I表示疾病暴发二次就不再流行;区域I V表示疾病暴发三次就不再流行;区域V表示疾病会多次暴发流行。图6 系统的盆吸引域F i g.6 B a s i n o f a t t r a c t i o n s o f s y s t e m研究表明,当感染者数量超过疾病暴发的阈值时,要加大媒体的宣传力度、报道范围以及报道频率,使得公众提高警惕并养成良好的预防习惯,从而使得疾病即使在暴发多次后也要稳定于安全可控的范围之内。这也说明感染者和易感者的初始密度对后续疾病的流行趋势和程度会产生重要的影响。因此,为了更有效地控制疾病传播流行,需要对易感者和感染者的密度进行实时监控和上报。3.4 多吸引子共存分析由3.3节相关内容可知,不同的初值会导致系统发生复杂多变的动力学行为以及致使系统稳定在不同的吸引子上,这就是吸引子共存现象,本节将关注系统多吸引子共存的研究。切换系统的吸引子同样具有很强的初值敏感性,为了研究不同初始密度对共存吸引子的影响,在不同的初值状态下选取参数b=2.8,K=1.8,=1.2,m=0.5,d=0.0 8,=0.2 5,=0.4 6,Ic=3,如图7所示。图7(a)为切换系统周期解情形,图7(b)为系统拟周期解情形。很显然,随着时间序列的增大,图7(b)所示情形必将导致系统呈混沌状态,这必然会给疾病防控带来一系列挑战。进一步,继续研究重要参数疾病传染率对系统多吸引子的影响,在初值(S0,I0)=(2,0.1)的情形下,固定参数b=3,K=2,m=0.1,d=0.5,a=0.8 5,=0.0 4,Ic=1.8,通过参数的不断变化,切换系统呈现出不同的振幅和频率的吸引子。随着参加不断增大,切换系统吸引子波动范围也不断变大,如图8(a)(f)所示。研究表明,疾病传染率在传染病传播中起着重要作用,防控部门可以通过媒体的宣传和报道降低疾病传染率,以达到有效控制传染病传播流行的目的。图7 系统不同初值情形下的多吸引子共存F i g.7 C o e x i s t i n g a t t r a c t o r s o f s y s t e m w i t h d i f f e r e n t i n i t i a l v a l u e s图8 不同参数b下系统的吸引子F i g.8 T h e a t t r a c t o r s o f s y s t e m w i t h d i f f e r e n t p a r a m e t e r s 4 结语考虑媒体报道对人们行为改变的切换效应,建立了一类由媒体效应诱导的离散传染病切换模型。一方面,从理论上研究了模型两个子系统的动力学行为,包括无病平衡点、地方病平衡点的存在性以及地方病平衡点的局部稳定性。另一方面,通过数值方法研究了敏感参数对系统动力学行为的影响。研究发现:关键参数的小扰动直接影响着疾病的暴发次数与频率,同时媒体报道可以降低疾病传染率,有效遏制疾病的流行与暴发。事实上,决定系统切换的阈值相当关键,就目前掌握的文献,还没有对如何确定这个阈值进行相关的研究,这也是一个非常具有挑战的问题。同时结合具体的突发性传染疾病实例(如新型冠状病毒感染),在媒体效应的基础上考虑阈值对疫情防控的影响,将具有重要的理论和现实意义,这将是后续需要研究的重要课题。104覃文杰,张嘉敏,向中义.媒体报道诱导的离散传染病切换系统建模与分析参考文献:1 马知恩.传染病动力学的数学建模与研究M.北京:科学出版社,2 0 0 4.MA Z h i e n.M a t h e m a t i c a l m o d e l i n g a n d r e s e a r c h o n t h e d y n a m i c s o f i n f e c t i o u s d i s e a s e sM.B e i j i n g:S c i e n c e P r e s s,2 0 0 4.2 TON I A T.S o c i a l m e d i a i n p u b l i c h e a l t h:I s i t u s e d a n d i s i t u s e f u l?J.I n t e r n a t i o n a l J o u r n a l o f P u b l i c H e a l t h,2 0 1 4,5 9(6):8 8 9-8 9 1.3 R UAN S h i g u i,WANG W e n d i.D y n a m i c a l b e h a v i o r o f a n e p i d e m i c m o d e l w i t h a n o n l i n e a r i n c i d e n c e r a t eJ.J o u r n a l o f D i f f e r e n t i a l E q u a t i o n s,2 0 0 3,1 8 8(1):1 3 5-1 6 3.4 L I U R o n g s o n g,WU J i a n h o n g,Z HU H u a i p i n g.M e d i a/p s y c h o l o g i c a l i m p a c t o n m u l t i p l e o u t b r e a k s o f e m e r g i n g i n f e c t i o u s d i s e a s e sJ.C o m p u t a t i o n a l a n d M a t h e m a t i c a l M e t h o d s i n M e d i c i n e,2 0 0 7,8(3):1 5 3-1 6 4.5 C U I J i n g a n,S UN Y o n g h o n g,Z HU H u a i p i n g.T h e i m p a c t o f m e d i a o n t h e c o n t r o l o f i n f e c t i o u s d i s e a s e sJ.J o u r n a l o f D y n a m i c s a n d D i f f e r e n t i a l E q u a t i o n s,2 0 0 8,2 0(1):3 1-5 3.6 刘玉英,肖燕妮.一类受媒体影响的传染病模型的研究J.应用数学和力学,2 0 1 3,3 4(4):3 9 9-4 0 7.L I U Y u y i n g,X I AO Y a n n i.A n e p i d e m i c m o d e l w i t h s a t u r a t e d m e d i a/p s y c h o l o g i c a l i m p a c tJ.A p p l i e d M a t h e m a t i c s a n d M e c h a n i c s,2 0 1 3,3 4(4):3 9 9-4 0 7.7 王俊荣,窦霁虹,孙梦皎.受媒体报道影响的具有L o g i s t i c人口变化的S I R S传染病模型分析J.陕西师范大学学报(自然科学版),2 0 1 8,4 6(6):4 6-5 3.WANG J u n r o n g,D OU J i h o n g,S UN M e n g j i a o.A n a l y s i s o f a n S I R S e p i d e m i c m o d e l w i t h m e d i a c o v e r a g e a n d L o g i s t i c v a r i a b l e p o p u l a t i o nJ.J o u r n a l o f S h a a n x i N o r m a l U n i v e r s i t y(N a t u r a l S c i e n c e E d i t i o n),2 0 1 8,4 6(6):4 6-5 3.8 S AHU G P,DHA R J.D y n a m i c s o f a n S E Q I HR S e p i d e m i c m o d e l w i t h m e d i a c o v e r a g e,q u a r a n t i n e a n d i s o l a t i o n i n a c o mm u n i t y w i t h p r e-e x i s t i n g i mm u n i t yJ.J o u r n a l o f M a t h e m a t i c a l A n a l y s i s a n d A p p l i c a t i o n s,2 0 1 5,4 2 1(2):1 6 5 1-1 6 7 2.9 S ONG P e n g f e i,X I AO Y a n n i.G l o b a l H o p f b i f u r c a t i o n o f a d e l a y e d e q u a t i o n d e s c r i b i n g t h e l a g e f f e c t o f m e d i a i m p a c t o n t h e s p r e a d o f i n f e c t i o u s d i s e a s eJ.J o u r n a l o f M a t h e m a t i c a l B i o l o g y,2 0 1 8,7 6(5):1 2 4 9-1 2 6 7.1 0 李芳,刘茂省.具有媒体饱和效应影响的时滞S I S模型研究J.数学的实践与认识,2 0 1 7,4 7(1 0):2 1 5-2 2 1.L I F a n g,L I U M a o x i n g.A c l a s s o f d e l a y S I S m o d e l w i t h m e d i a s a t u r a t i o nJ.M a t h e m a t i c s i n P r a c t i c e a n d T h e o r y,2 0 1 7,4 7(1 0):2 1 5-2 2 1.1 1 B A S I R F A,R AY S,V E N TUR I NO E.R o l e o f m e d i a c o v e r a g e a n d d e l a y i n c o n t r o l l i n g i n f e c t i o u s d i s e a s e s:A m a t h e m a t i c a l m o d e lJ.A p p l i e d M a t h e m a t i c s a n d C o m p u t a t i o n,2 0 1 8,3 3 7:3 7 2-3 8 5.1 2 S ON G P e n g f e i,X I AO Y a n n i.A n a l y s i s o f a n e p i d e m i c s y s t e m w i t h t w o r e s p o n s e d e l a y s i n m e d i a i m p a c t f u n c t i o nJ.B u l l e t i n o f M a t h e m a t i c a l B i o l o g y,2 0 1 9,8 1(5):1 5 8 2-1 6 1 2.1 3 D E NG J i a w e i,TANG S a n y i,S HU H o n g y i n g.J o i n t i m p a c t s o f m e d i a,v a c c i n a t i o n a n d t r e a t m e n t o n a n e p i d e m i c F i l i p p o v m o d e l w i t h a p p l i c a t i o n t o C OV I D-1 9J.J o u r n a l o f T h e o r e t i c a l B i o l o g y,2 0 2 1,5 2 3:1 1 0 6 9 8.1 4 K E RMA C K W O,MC K E N D R I C K A G.C o n t r i b u t i o n s t o t h e m a t h e m a t i c a l t h e o r y o f e p i d e m i c s-IJ.B u l l e t i n o f M a t h e m a t i c a l B i o l o g y,1 9 9 1,5 3(1/2):3 3-5 5.1 5 A L L E N L J.S o m e d i s c r e t e-t i m e S I,S I R,a n d S I S e p i d e m i c m o d e l sJ.M a t h e m a t i c a l B i o s c i e n c e s,1 9 9 4,1 2 4(1):8 3-1 0 5.1 6 C A S T I L L O-CHAV E Z C,YAKU B U A.D i s p e r s a l,d i s e a s e a n d l i f e-h i s t o r y e v o l u t i o nJ.M a t h e m a t i c a l B i o s c i e n c e s,2 0 0 1,1 7 3(1):3 5-5 3.1 7 Z HOU Y i c a n g,F E R GO L A P.D y n a m i c s o f a d i s c r e t e a g e-s t r u c t u r e