2018年具有隔离和脉冲预防接种的SIQR模型的稳定性分析.docVIP

具有隔离和脉冲预防接种的SIQR模型的稳定性分析 　　第38卷第2期2012年4月曲阜师范大学学报JournalofQufuNormalUniversityVol．38No．2Apr．2012 　　具有隔离和脉冲预防接种的SIQR模型的稳定性分析 　　付 　　敏，宋 　　燕 　　* 　　(渤海大学数理学院，121000，辽宁省锦州市) 　　摘要:研究了具有指数出生，标准发生率的脉冲预防接种的SIQR模型．讨论了其无病周期解的存在性，给 　　出了无病周期解的全局稳定性分析．并利用数值仿真对结论进行了验证，最后获得了系统一致持续生存的条件． 　　关键词:隔离;接种;周期解;稳定性;一致持续生存 　　5337(2012)02-0023-06中图分类号:O193文献标识码:A文章编号:1001-为控制传染病的传播与蔓延，人们利用数学工具建立数学模型来研究疾病的传播规律．控制疾病最直接 　　的措施就是对染病人群进行隔离和对易感人群进行预防接种．对于具有隔离项的传染病模型已经有很多研 　　［4-6］ 　　1-3］，．董淑转研究了带有隔离项和连续预防接种的究成果，如文献［具有接种的传染病模型也备受关注 　　［7］ 　　传染病模型．但隔离和脉冲接种结合讨论的传染病模型较少．本文研究了一类具有指数出生，发生率为标 　　准发生率，具有隔离措施和脉冲接种的SIQR模型，利用频闪映射讨论了其无病周期解的存在性，利用Flo-quet定理证明了无病周期解的局部渐近稳定性，运用脉冲不等式证明了无病周期解的全局渐近稳定性，通过数值仿真验证了理论结果．最后，给出了系统一致持续生存的条件． 　　1模型建立与分析 　　b为出生率，d为自然死亡率，把总人口N分为易感者S，染病者I，隔离者Q和恢复者R，β为一个染病者 　　p为预防接种率．所具有的传染力，δ为隔离系数，γ，ε分别为染病者和隔离者的移出率， 　　具有指数出生，脉冲接种的SIQR模型的脉冲微分方程如下 　　?dS=bN－dS－βSI， 　　N?dt 　　?dISI 　　=β－(d+δ+γ)I，t≠tn，??dtN? 　　?dQ=δI－(d+ε)Q，?dt??dR=I+Q－dR． 　　γε?dt 　　+S(t)???I(t+)?+Q(t)?? 　　?R(t+) 　　(1) 　　=(1－p)S(t)，=I(t)，t=tn， 　　=Q(t)，n=0，1，2，…，=R(t)+pS(t)． 　　SIQR，i=，q=，r=， 　　NNNN 　　(2) 　　为处理方便，作下述归一化变换 　　s= 　　* 　　12-09收稿日期:2011-基金项目:辽宁省教育厅高等学校科学研究基金资助项目(2008009)． 　　1962-，作者简介:宋燕，女，教授;研究方向:非线性生物动力系统．Email:jzsongyan@163．com． 　　(2)变成则s+i+q+r≡1．从而方程组(1)， 　　?ds=b－bs－βsi，?dt?di?=βsi－(b+δ+γ)i，t≠tn，dt??dq 　　?dt=δi－(b+ε)q，? 　　?r=1－s－i－q． 　　+s(t)???i(t+)?+??q(t)?r(t+) 　　(3) 　　=(1－p)s(t)，=i(t)，t=tn，=q(t)，n=0，1，2，…，=r(t)+ps(t)． 　　(4) 　　3 　　i，q)∈［0，1］|0≤s+i+q≤1}．可行域为:D={(s， 　　2无病周期解的存在性和稳定性分析 　　(1)无病周期解的存在性 　　(4)变为当I(t)=Q(t)=0(t≥0)时，方程组(3)， 　　s=b－bs，t≠tn， 　　{ 　　s(t+)=(1－p)s(t)，t=tn， 　　(5) 　　方程组(5)在tn≤t≤tn+1上的解为 　　{ 　　+s(t)=1+(s(tn)－1)exp(－b(t－tn))，tn≤t＜tn+1， 　　s(t)=(1－p)s(t)， 　　+ 　　n 　　t=tn+1． 　　(6) 　　+ 　　令sn+1=s(tn+1)，由(6)式可得sn+1=(1－p)(1+(sn－1)exp(－bT))．设频闪映射F:sn→sn+1，满足sn+1=F(sn)=(1－p)(1+(sn－1)exp(－bT))，该式存在唯一的正不动点 　　(1－p)(exp(bT)－1) 　　s*(t)=． 　　exp(bT)+p－1 　　因为 　　dF(sn)ds 　　* 　　=(1－p)exp(－bT)＜1，(4)存在无病所以不动点是稳定的，即sn→s．故方程组(3)，sn=s* 　　～～ 　　T周期解(s(t)，0，0)，其中s(t)=1－ 　　pexp(bT) 　　exp(－b(t－tn))，tn＜t≤tn+1． 　　exp(bT)+p－1 　　(2)无病周期解的稳定性