一类随机SIR传染病模型的非标准数值离散化 (1).pdf

西 北 师 范 大 学 学 报(自然科学版)第6 0卷2 0 2 4年第1期 J o u r n a l o f N o r t h w e s t 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)V o l.6 0 2 0 2 4 N o.1 D O I:1 0.1 6 7 8 3/j.c n k i.n w n u z.2 0 2 4.0 1.0 0 6收稿日期:2 0 2 2 1 1 0 8;修改稿收到日期:2 0 2 3 0 5 2 2基金项目:国家自然科学基金资助项目(1 2 0 7 1 4 4 5,1 2 0 0 1 5 0 1)作者简介:谭伟(1 9 9 6),男,河北张家口人,硕士研究生.主要研究方向为传染病动力学.E m a i l:1 8 4 8 2 3 1 7 2 2q q.c o m*通信联系人,男,教授,博士.主要研究方向为传染病动力系统及复杂网络.E m a i l:l i u m a o x i n g 1 2 6.c o m一类随机S I R传染病模型的非标准数值离散化谭 伟1,刘茂省1,2*(1.中北大学 数学学院,山西 太原 0 3 0 0 5 1;2.北京建筑大学 理学院,北京 1 0 2 6 1 6)摘要:研究一个具有饱和发生率和疫苗接种率的随机离散S I R传染病模型的稳定性.首先,引入一个具有饱和发生率和疫苗接种率的确定性S I R模型,考虑到随机噪声对疾病传播有很大影响,应用非标准有限差分(N S F D)方法将模型离散化,最终得到一个随机离散的S I R传染病模型.这种离散方法是对系统的右侧进行局部离散,得出离散模型,然后系统左侧用广义前向差分法对一阶导数进行逼近,并且要选取恰当的分母函数.其次,应用L y a p u n o v函数和矩阵方法给出了系统平衡解稳定的充分条件,并且得到了非线性差分方程概率稳定的充分条件和线性差分方程渐近均方稳定的充分条件.最后,通过数值仿真对理论分析结论进行验证.关键词:饱和发生率;随机离散模型;非标准有限差分方法;L y a p u n o v函数;渐近均方稳定中图分类号:O 1 7 5.7 文献标志码:A 文章编号:1 0 0 1-9 8 8(2 0 2 4)0 1-0 0 3 0-0 9N o n s t a n d a r d n u m e r i c a l d i s c r e t i z a t i o n o fa s t o c h a s t i c S I R e p i d e m i c m o d e lT AN W e i1,L I U M a o-x i n g1,2(1.C o l l e g e o f M a t h e m a t i c s,N o r t h U n i v e r s i t y o f C h i n a,T a i y u a n 0 3 0 0 5 1,S h a n x i,C h i n a;(2.C o l l e g e o f S c i e n c e,B e i j i n g U n i v e r s i t y o f C i v i l E n g i n e e r i n g a n d A r c h i t e c t u r e,B e i j i n g 1 0 2 6 1 6,C h i n a)A b s t r a c t:T h i s p a p e r s t u d i e s t h e s t a b i l i t y o f a s t o c h a s t i c d i s c r e t e S I R 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 i n c i d e n c e a n d v a c c i n a t i o n r a t e.A d e t e r m i n i s t i c S I R m o d e l w i t h s a t u r a t e d i n c i d e n c e a n d v a c c i n a t i o n r a t e i s i n t r o d u c e d.C o n s i d e r i n g t h e g r e a t i n f l u e n c e o f s t o c h a s t i c n o i s e o n d i s e a s e t r a n s m i s s i o n,t h e m o d e l i s d i s c r e t i z e d b y n o n s t a n d a r d f i n i t e d i f f e r e n c e(N S F D)m e t h o d,a n d f i n a l l y a s t o c h a s t i c d i s c r e t e S I R e p i d e m i c m o d e l i s o b t a i n e d.T h i s d i s c r e t i z a t i o n m e t h o d i s t o l o c a l l y d i s c r e t i z e t h e r i g h t s i d e o f t h e s y s t e m t o o b t a i n t h e d i s c r e t e m o d e l,a n d t h e n u s e t h e g e n e r a l i z e d f o r w a r d d i f f e r e n c e m e t h o d t o a p p r o x i m a t e t h e f i r s t d e r i v a t i v e o n t h e l e f t s i d e o f t h e s y s t e m,a n d s e l e c t t h e a p p r o p r i a t e d e n o m i n a t o r f u n c t i o n.T h e s u f f i c i e n t c o n d i t i o n s f o r t h e s t a b i l i t y o f t h e e q u i l i b r i u m s o l u t i o n o f t h e s y s t e m a r e g i v e n b y u s i n g L y a p u n o v f u n c t i o n m e t h o d a n d m a t r i x m e t h o d,a n d t h e s u f f i c i e n t c o n d i t i o n s f o r t h e p r o b a b i l i s t i c s t a b i l i t y o f n o n l i n e a r d i f f e r e n c e e q u a t i o n s a n d t h e s u f f i c i e n t c o n d i t i o n s f o r t h e a s y m p t o t i c m e a n s q u a r e s t a b i l i t y o f l i n e a r d i f f e r e n c e e q u a t i o n s a r e a l s o p r o p o s e d.F i n a l l y,t h e c o n c l u s i o n i s v e r i f i e d b y n u m e r i c a l s i m u l a t i o n.K e y w o r d s:s a t u r a t i o n i n c i d e n c e;s t o c h a s t i c d i s c r e t e m o d e l;n o n s t a n d a r d f i n i t e d i f f e r e n c e m e t h o d;L y a p u n o v f u n c t i o n;a s y m p t o t i c m e a n s q u a r e s t a b i l i t y03 2 0 2 4年第1期 谭 伟等:一类随机S I R传染病模型的非标准数值离散化 2 0 2 4 N o.1N o n s t a n d a r d n u m e r i c a l d i s c r e t i z a t i o n o f a s t o c h a s t i c S I R e p i d e m i c m o d e l0 引言近几十年来,数学模型被认为是理解传染病传播机制和评估其控制策略的关键工具1-3.在这方面,Y u s u f等4引入并分析了一个具有疫苗接种和治疗的S I R流行模型,其中假设发病率是双线性的.但经过长期的验证发现,双线性发生率对于较大的种群规模来说并不符合实际,为了应对易感个体产生的心理或抑制效应,C a p a s s o等5提出了饱和发病率S I/(1+I),其中是疾病的传染率,是衡量由拥挤效应或社会行为变化引起的抑制效应的饱和因子.当前,饱和发病率已广泛应用于确定性传染病模型6-9和随机传染病模型1 0-1 2.本文研究具有饱和发生率和疫苗接种率的S I R确定性传染病模型1 3:dSdt=-S I1+I-(+)S,dIdt=S I1+I-(+)I,dRdt=S+I-R.(1)其中SS(t)为t时刻易感人群数量,II(t)为t时刻已感染人群数量,RR(t)为t时刻染病后免疫人群的数量,且S(0)0,I(0)0,R(0)0;为从外部进入种群的净迁移率,为人口的自然死亡率,为接种疫苗的人口比率,为患病种群恢复的概率,为因为感染疾病而死亡的人口概率,且所有参数都是正的.注意到系统(1)的第一和第二个方程中不包含R,系统(1)的第三个式子不影响其动力学行为,因此,我们只需研究以下系统:dSdt=-S I1+I-(+)S,dIdt=S I1+I-(+)I.(2)通过简单计算可知,系统(2)中存在一个无病平衡点E0和一个地方病平衡点E+,并且无论在什么条件下,无病平衡点E0=(S0,I0)=+,0 都存在.当R01时地方病平衡点E+=(S+,I+)存在,且S+=+(+)(+)+,I+=(+)(R0-1)(+)+,基本再生数R0=(+)(+).由于白噪声对疾病的传播有很大的影响,所以与确定性模型相比,随机模型更加符合实际.本文在模型(2)的基础上加入随机噪声,得到dSdt=-S I1+I-(+)S+1(S-S*)dB1(t),dIdt=S I1+I-(+)I+2(I-I*)dB2(t),(3)其中Bi(t)(i=1,2)为白噪声,也就是常说的布朗运动,2i(i=1,2)为白噪声Bi(t)的强度1 4.假设随机扰动强度与状态S(t)和I(t)偏离平衡点的偏差成正比,即当系统状态偏离平衡点的偏差增大时,随机扰动强度也随之增加1 5.传染病模型通常是一个非线性微分方程组,因此很难得到其解析解,通常需要数值近似来得到可靠的结果.然而,许多标准的数值方法,如E u l e r方法、R u n g e-K u t t a方法可能无法保持连续模型的动力学特性,例如收敛到错误的平衡点或数值不稳定1 6-1 8.为克服这些方法的缺点,M i c k e n s1 9提出了非标准有限差分(N S F D)方法.如果以下两个条件之一被满足,那么一阶微分方程组的这种数值格式即为N S F D格式:(i)用广义前向差分近似系统中的一阶导数dundtu(n+1)-u(n),其中un=u(tn),分母函数=(h)满足(h)=h+O(h2),h为时间步长.(i i)非线性项是非局部近似的,例如u2(tn)unun+1.然而,注意到N S F D格式对一个微分方程来说并不是唯一的,所以为了确保N S F D格式与原始连续模型具有相同的动力学性态,我们需要对N S F D格式进行动力学特性的分析.假定条件(i)成立,我们根据系统(3)提出以下N S F D格式:13西 北 师 范 大 学 学 报(自然科学版)第6 0卷 J o u r n a l o f N o r t h w e s t 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)V o l.6 0 S(n+1)-S(n)=-S(n)I(n)1+I(n)-(+)S(n)+11(S(n)-S*)1(n+1),I(n+1)-I(n)=S(n)I(n)1+I(n)-(+)I(n)+21(I(n)-I*)2(n+1),(4)其中分母函数是与h相关的函数,初始条件为S(0)=1(0),I(0)=2(0).系统(4)是对系统(3)进行广义前向差分和右侧局部近似得到的,通过一些简单的代数计算,我们可以 很 容 易 地 验 证,具 有N S F D格 式 的 模 型(4)与初始连续模型(2)具有完全相同的平衡点.数学期望我们通常由E来表示,i(n+1)(i=1,2)(nZ)是一个相互独立且服从标准正态分布N(0,1)的随机序列,并且对于i(n+1)(i=1,2,nZ),满足Ei(n)=0,E2i(n)=1,Ei(n)Ej(n),ij.(5)通过直接计算,容易验证分母函数2 0(h)=1-e-h.(6)对离散N S F D系统(4)可重新排列,得到其显式形式:S(n+1)=S(n)+-S(n)I(n)1+I(n)-(+)S(n)+1(S(n)-S*)1(n+1),I(n+1)=I(n)+S(n)I(n)1+I(n)-(+)I(n)+2(I(n)-I*)2(n+1).(7)近几十年来,具有噪声扰动的确定性模型的随机流行病模型2 1-2 3已经被广泛研究,并且越来越多的学者将差分方程模型2 4,2 5运用到传染病学的建模中.然而,多数学者都是将随机模型和差分格式分别建模来研究的2 1,但在传染病学中对差分方程所描述的随机离散模型的稳定性的研究很少,所以本文主要对随机离散的传染病模型进行研究.1 预备知识对系统(7)的地方病平衡点进行平移变换,令u(n)=S(n)-S+,m(n)=I(n)-I+可得u(n+1)=u(n)+-(u(n)+S+)(m(n)+I+)1+(m(n)+I+)-(u+v)(u(n)+S+)+1u(n)1(n+1),(8)m(n+1)=m(n)+(u(n)+S+)(m(n)+I+)1+(m(n)+I+)-(u+)(m(n)+I+)+2m(n)2(n+1).系统(8)的 零 解 等 价 于 系 统(7)的 正 解E+=(S+,I+).在平衡点E+=(S+,I+)处对系统(7)进行线性化,得到线性近似系统x(n+1)=1-I+1+I+-(+v)x(n)-S+(1+I+)2y(n)+1x(n)1(n+1),y(n+1)=I+1+I+x(n)+(9)1+S+(1+I+)2-(+)y(n)+2y(n)2(n+1).类似地进行平移变换.对无病平衡点,令u(n)=S(n)-/(+v),v(n)=I(n),则u(n+1)=u(n)+-u(n)+vm(n)1+m(n)-(u+v)u(n)+1u(n)1(n+1),(1 0)m(n+1)=m(n)+u(n)+v m(n)1+m(n)-(u+)m(n)+2m(n)2(n+1).23 2 0 2 4年第1期 谭 伟等:一类随机S I R传染病模型的非标准数值离散化 2 0 2 4 N o.1N o n s t a n d a r d n u m e r i c a l d i s c r e t i z a t i o n o f a s t o c h a s t i c S I R e p i d e m i c m o d e l在无病平衡点处对系统(7)线性化,得到相应的线性近似系统为x(n+1)=(1-(+v)x(n)-y(n)+1x(n)1(n+1),y(n+1)=1+-(+)y(n)+2y(n)2(n+1).(1 1)令(n)=(u(n),m(n)T,z(n)=(x(n),y(n)T,(n)=(1(n),2(n)T,T代表转置.下面引入了文献2 6 的一些定义和定理来研究系统的动态行为.定义1(文献2 6 定义7.1)对初始函数S(0)=1(0),I(0)=2(0),如果对10和20,存在一个0使得系统(8)或(1 0)的解(n)=(n,)满足不等式P s u pnZ(n)1 2和P0,存在()0,使得Ez(n)2,nZ,并且E(n)2Q.定理B(文献2 6 定理5.1)假设存在正定矩阵P,使得形为(1 3)的半正定矩阵D满足矩阵方程ATDA-D=-P,同时P要满足条件 Pd1 12100d2 222 ,(1 4)那么系统(1 2)的零解是渐近均方稳定的.证明 由(1 3)式,系统(1 2)可表示为z(n+1)=(A+(n+1)z(n),其中z(n)=x(n)y(n),(n+1)=11(n+1)0022(n+1).考虑L y a p u n o v函数V(n)=zT(n)D z(n),则V(n)=V(n+1)-V(n)=zT(n+1)D z(n+1)-zT(n)D z(n).然后计算V的期望,得到EV(n)=EzT(n+1)D z(n+1)-zT(n)D z(n)=EzT(n)(A+(n+1)T D(A+(n+1)-Dz(n)=EzT(n)ATDA-D+T(n+1)D(n+1)z(n)=E zT(n)-P+T(n+1)D(n+1)z(n).由关系式(5),可以求得33西 北 师 范 大 学 学 报(自然科学版)第6 0卷 J o u r n a l o f N o r t h w e s t 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)V o l.6 0 ET(+1)D)D(+1)=E 1w1(n+1)002w2(n+1)d1 1d1 2d1 2d2 2 1w1(n+1)002w2(n+1)=d1 12100d2 222 ,然后根据(1 4)式,可以得到EV(n)=EzT(n)-P+d1 12100d2 222 z(n)-c Ez(n)2,其中c是正数.依据定理A可知系统(1 2)的零解是渐近均方稳定的.】2 地方病平衡点的稳定性因为系统(9)是系统(1 2)的特殊形式,可以推导出系统(9)的系数矩阵A=A1 1A1 2A2 1A2 2 ,(1 5)其中A1 1=1-I+1+I+-(+v),A1 2=-S+(1+I+)2,A2 1=I+1+I+,A2 2=1+S+(1+I+)2-(+).从矩阵方程ATDA-D=-P可解出-p1 1=(a21 1-1)d1 1+2a1 1a2 1d1 2+a22 1d2 2,-p1 2=a1 1a1 2d1 1+(a2 1a1 2+a1 1a2 2-1)d1 2+a2 1a2 2d2 2,-p2 2=a21 2d1 1+2a1 2a2 2d1 2+(a22 2-1)d2 2.(1 6)将矩阵A的元素带入(1 6)式可得p1 1=-1-I+1+I+-(+v)2-1 d1 1-2I+h1-I+1+I+-(+v)1+I+d1 2-I+1+I+2d2 2,p1 2=1-I+1+I+-(+v)S+(1+I+)2 d1 1+I+1+I+S+(1+I+)2 d1 2-1-I+1+I+-(+v)1-S+(1+I+)2-(+)-1d1 2-I+1+I+1+S+(1+I+)2-(+)d2 2,p2 2=-S+(1+I+)2 2d1 1+2S+(1+I+)2 1+S+(1+I+)2-(+)d1 2-1+S+(1+I+)2-(+)2-1d2 2.(1 7)再由(1 4)式可导出 p1 1-d1 1210,p1 1p2 2-p21 20,(p1 1-d1 121)(p2 2-d2 222)-p21 20.(1 8)根据以上的分析我们得出以下推论.推论1 若存在一个正定矩阵P,使得条件(1 8)对于矩阵方程(1 7)的解(d1 1,d1 2,d2 2)被满足,并且矩阵D是半正定的,那么系统(9)的零解是渐近均方稳定的.因此,系统(8)的零解是依概率稳定的,这等同于系统(7)的正平衡点Ee是依概率稳定的.当R01时,系统(7)存在地方病平衡点E+=(S+,I+).对于数值仿真,假设初值S(0)=5,I(0)=2.(1 9)其它参数设置为=4,=0.0 1,=0.5,=0.5,=0.1,=0.2,=0.1,h=0.1.(2 0)通过这些参数能够求出分母函数为=0.0 9 7 6,基本再生数R0=4.0 8 1 61,从而地方病平衡点存在,且E+=(1.4 5 9 7,4.2 6 2 9).通过对(1 5)求解得到矩阵43 2 0 2 4年第1期 谭 伟等:一类随机S I R传染病模型的非标准数值离散化 2 0 2 4 N o.1N o n s t a n d a r d n u m e r i c a l d i s c r e t i z a t i o n o f a s t o c h a s t i c S I R e p i d e m i c m o d e lA=0.7 3 2 2-0.0 6 5 50.1 9 9 50.9 9 7 2 ;假设P为单位矩阵,即P=1 00 1 ,求解(1 7)式得到半正定矩阵D=7.8 6 9 07.3 8 4 47.3 8 4 4 1 2.3 9 8 1 .由推论1可知,若1,2满足条件11d1 1=1.1 4 1 0,21d2 2=0.9 0 9 1,(2 1)那么,系统(9)的零解是渐近均方稳定的.因此,系统(8)的零解是依概率稳定的,这等同于系统(7)的地方病平衡点E+是依概率稳定的.下面借助数值模拟证明之前的分析.(1 9)和(2 0)式给出了初值和一些参数值.图1显示了噪声扰动强度1=2=0的情况,系统(9)和线性系统(7)的解曲线分别由图1(a)和图1(b)表示.从中可以看出,线性系统(9)的解曲线最终收敛到0,系统(7)的解曲线最终收敛于地方病平衡点E+=图1 1=2=0时系统(9)和(7)的解曲线F i g 1 T h e s o l u t i o n o f s y s t e m(9)a n d(7)w h e n 1=2=0(1.4 5 9 7,4.2 6 2 9).因此,系统(9)的零解是渐近均方稳定的,系统(7)的地方病平衡点E+是依概率稳定的.对具有随机扰动的系统(7)的解进行模拟.(1 9)和(2 0)式已给出初值和其它的参数值.假设1,2的值满足条件(2 1),我们取1=2=0.8,得到具有噪声扰动的系统(7)的解轨迹,见图2.从图2可以看出,随机噪声并不影响曲线的最终走向,系统(7)的解轨迹最终还是收敛到地方病平衡点E+,这就表明系统(7)在正平衡点处是依概率稳定的.图2 1=2=0.8时系统(7)的解轨迹F i g 2 S o l u t i o n t r a j e c t o r i e s o f s y s t e m(7)w h e n 1=2=0.83 无病平衡点的稳定性由系统(1 1)可以写出它的系数矩阵A=1-(+v)-01+-(+).(2 2)同理,由方程(1 6)解得p1 1=-(1-(+)2-1)d1 1,p1 2=(1-(+)d1 1-(1-(+)1+-(+)-1d1 2,p2 2=-2d1 1+2 1+-(+)d1 2-1+-(+)2-1 d2 2.(2 3)53西 北 师 范 大 学 学 报(自然科学版)第6 0卷 J o u r n a l o f N o r t h w e s t 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)V o l.6 0 再由(1 8)式,我们同样能推导出以下结论.推论2 对于(2 2)式中的矩阵A,如果存在一个正定矩阵P,使条件(1 8)被方程(2 3)的解(d1 1,d1 2,d2 2)所满足,且矩阵D为半正定的,那么称系统(1 1)的零解是渐近均方稳定的.因此,系统(1 0)的零解是依概率稳定的,这等同于是依概率稳定的系统(7)的边界平衡点E0.当R01时,系统(7)不存在地方病平衡点,只有一个无病平衡点.取参数值=4,=0.0 1,=0.9,=0.0 1,=0.2,=0.1,h=0.1.(2 4)此时R0=0.1 6 71,所以只有一个无病平衡点存在,且E0=(2 0,0).借助(6)式求得分母函数=0.0 9 9 5,给定初值为S(0)=1 0,I(0)=5.由以上参数能够求出(2 2)式中的矩阵A=0.9 8 0 1-0.0 3 9 800.9 2 0 4 .仍取P为单位矩阵,借助(2 3)式求得半正定矩阵D=2 5.3 8 0 7-1 0.1 1 2 4-1 0.1 1 2 41 1.4 2 0 8 .并由推论1可知,若1,2满足条件:11d1 1=0.6 2 9 3,21时,考虑系统(7)中地方病平衡点存在时无病平衡点的稳定性.仍取(1 9)和(2 0)式中的初值和参数,此时求得(2 2)式中的矩阵A=0.9 3 1 7-0.3 9 0 401.3 2 2 1 .P仍为单位矩阵,推导出矩阵D=7.5 1 1 51 1.7 8 6 91 1.7 8 6 9 1 3.4 0 1 3 ,通过计算可知,矩阵D不是半正定的,所以根据推论1可知,系统(7)的无病平衡点不稳定.63 2 0 2 4年第1期 谭 伟等:一类随机S I R传染病模型的非标准数值离散化 2 0 2 4 N o.1N o n s t a n d a r d n u m e r i c a l d i s c r e t i z a t i o n o f a s t o c h a s t i c S I R e p i d e m i c m o d e l当随机噪声强度为0时,图5中代表染病者的曲线I趋向于无穷,所以此时系统(1 1)的零解是不稳定的.当加入一个较小的噪声强度后,系统的解轨迹是非常不规则的(图6).因此,当一个正平衡点存在时,无论是否存在随机扰动,系统(7)的边界平衡点都是不稳定的.图5 1=2=0时系统(1 1)的解曲线F i g 5 S o l u t i o n o f s y s t e m(1 1)w h e n 1=2=0图6 1=2=0.3时系统(7)的解轨迹F i g 6 S o l u t i o n t r a j e c t o r i e s o f s y s t e m(7)w h e n 1=2=0.34 结论本文讨论了一个具有饱和发生率和疫苗接种率的连续S I R传染病模型,鉴于疾病在传播过程中会受到各种外部因素的影响,例如随机白噪声的影响,我们在系统中加入随机扰动项.由于离散模型的动力学性态比连续模型表现的更加丰富和准确,我们还构造了保持连续模型基本性质的非标准有限差分格式,并选择了一个合适的分母函数.随后用L y a p u n o v函数方法探究了所构建的随机差分模型的动力学性质,建立了模型在平衡点处稳定的充分条件.数值模拟表明,当所选参数和噪声强度满足给定条件时,适当的噪声强度不会影响系统在平衡点处的稳定性;并且当正平衡点存在时,无论随机扰动强度是多少,边界平衡点的稳定性都不会被改变,在这种情况下,边界平衡点是不稳定的.参考文献:1 C A S T I L L O CHAV E Z C.M a t h e m a t i c a l M o d e l s i n P o p u l a t i o n B i o l o g y a n d E p i d e m i o l o g yM.N e w Y o r k:S p r i n g e r,2 0 1 2.2 B R AU E R F.M a t h e m a t i c a l e p i d e m i o l o g y:p a s t,p r e s e n t,a n d f u t u r eJ.I n f e c t i o u s D i s e a s e M o d e l l i n g,2 0 1 7(2):1 1 3.3 R E Z A P OUR S,MOHAMMA D I H,S AME I M E.S E I R e p i d e m i c m o d e l f o r C OV I D-1 9 t r a n s m i s s i o n b y C a p u t o d e r i v a t i v e o f f r a c t i o n a l o r d e rJ.A d v a n c e s i n D i f f e r e n c e E q u a t i o n s,2 0 2 0,2 0 2 0(1):4 9 0.4 YU S U F T T,B E NYAH F.O p t i m a l c o n t r o l o f v a c c i n a t i o n a n d t r e a t m e n t f o r a n S I R e p i d e m i o l o g i c a l m o d e lJ.W o r l d J o u r n a l o f M o d e l l i n g&S i m u l a t i o n,2 0 1 2,8(3):1 9 4.5 C A P A S S O V,S E R I O G.A g e n e r a l i z a t i o n o f t h e K e r m a c k-M c K e n d r i c k d e t e r m i n i s t i c e p i d e m i c m o d e lJ.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 7 8,4 2(1/2):4 3.6 XU R,MA Z,WANG Z.G l o b a l s t a b i l i t y o f a d e l a y e d S I R S e p i d e m i c m o d e l w i t h s a t u r a t i o n i n c i d e n c e a n d t e m p o r a r y i mm u n i t yJ.C o m p u t e r s&M a t h e m a t i c s w i t h A p p l i c a t i o n s,2 0 1 0,5 9(9):3 2 1 1.7 J ANA S,NAN D I S K,KA R T K.C o m p l e x d y n a m i c s o f a n S I R 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 i n c i d e n c e r a t e a n d t r e a t m e n tJ.A c t a B i o t h e o r e t i c a,2 0 1 6,6 4:6 5.8 MA THUR K S,NA R AYAN P.D y n a m i c s o f a n e p i d e m i c m o d e l w i t h v a c c i n a t i o n a n d s a t u r a t e d i n c i d e n c e r a t eJ.N a t u r e P u b l i c H e a l t h Em e r g e n c y C o l l e c t i o n,2 0 1 8,4(5):1.9 L I U J.B i f u r c a t i o n a n a l y s i s f o r a d e l a y e d S E I R 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 i n c i d e n c e a n d s a t u r a t e d t r e a t m e n t f u n c t i o nJ.J o u r n a l o f B i o l o g i c a l D y n a m i c s,2 0 1 9,1 3(1):4 6 1.1 0 R A J A S E KA R S P,P I T CHA I MAN I M,Z HU Q.D y n a m i c t h r e s h o l d p r o b e o f s t o c h a s t i c S I R m o d e l w i t h s a t u r a t e d i n c i d e n c e r a t e a n d s a t u r a t e d t r e a t m e n t f u n c t i o nJ.P h y s i c a A:S t a t i s t i c a l M e c h a n i c s a n d i t s A p p l i c a t i o n s,2 0 1 9,5 3 5:1 2 2 3 0 0.1 1 R A J A S E KA R S P,P I T CHA I MAN I M,Z HU Q.P r o g r e s s i v e d y n a m i c s o f a s t o c h a s t i c e p i d e m i c m o d e l w i t h l o g i s t i c g r o w t h a n d s a t u r a t e d t r e a t m e n tJ.P h y s i c a A:S t a t i s t i c a l M e c h a n i c s a n d i t s A p p l i c a t i o n s,2 0 1 9,5 3 6:1 2 2 6 4 9.73西 北 师 范 大 学 学 报(自然科学版)第6 0卷 J o u r n a l o f N o r t h w e s t 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)V o l.6 0 1 2 L I D,Z HAO Y,S ON G S.D y n a m i c a n a l y s i s o f s t o c h a s t i c v i r u s i n f e c t i o n m o d e l w i t h d e l a y e f f e c tJ.P h y s i c a A:S t a t i s t i c a l M e c h a n i c s a n d i t s A p p l i c a t i o n s,2 0 1 9,5 2 8:1 2 1 4 6 3.1 3 S UR YAN TO A,D A R T I I.O n t h e n o n s t a n d a r d n u m e r i c a l d i s c r e t i z