【摘 要】本文建立并研究了一类宿主具有垂直传染和预防接种的的媒介传染病模型，给出了疾病流行与否的阈值并讨论了平衡点的存在性。运用Hurwitz判据证明了当基本再生数R01时，无病平衡点是局部渐近稳定的;当R01时，地方病平衡点也是局部渐近稳定的，并对证明结果进行了数值模拟。
【关键词】媒介传染;垂直传染;动力学系统;基本再生数
Analysis of a Kind of Host―Vector Epidemic Model with Vertical Transmission
XIAO Ya-nan XUE Ya-kui
（School of Science， North University of China， Taiyuan Shanxi 030051， China）
【Abstract】The Host-Vector epidemic model with vertical transmission and vaccination is formulated and studied. The threshold is identified which determines the outcome of disease and the existence of the equilibrium is discussed. The model is shown that the disease free equilibrium is locally stability when R01 and the endemic equilibrium is locally asymptotically stable when R01. The numerical simulations are carried out to prove the results.
【Key words】Host-vector; Vertical transmission; Dynamic modeling; Basic reproduction number
0 引言
1 动力学模型
考虑宿主为具有垂直传染和预防接种的SIR模型，假设染病媒介叮咬易感人群使人患病，而易感媒介通过叮咬染病人群而被感染。将人口分为易感类、染病者类和康复类，用SH（t），IH（t），RH（t）分别表示t时刻的易感，染病，康复类的人群数量，NH（t）=SH（t）+IH（t）+RH（t）。假设不考虑因病死亡，用b和b′分别表示非染病者（S+R）和染病者I的出生率系数;d和d′是相应的死亡率系数，γ表示染病者的康复率，q是垂直传染率（p+q=1），m是对易感者和康复者的新生儿进行预防接种的比例。对于媒介，分为易感，感染两类，用SV（t），IV（t）分别表示t时刻的易感、染病媒介的数量，NV（t）=SV（t）+IV（t）。μ是媒介的出生和死亡率系数。β1表示染病媒介对易感人群的感染率系数，βV表示染病人群对易感媒介的感染率。假设康复者终身免疫该疾病，令b=b′，d=d′，根据仓室模型思想可以得到如下传染病动力学模型：
■
（1）
由于■，■，故人口总数和媒介总数均为常数，不妨设为NH=SH+IH+RH=1，NV=IV+RV=1，我们只需考虑如下系统：
■
（2）
2 再生数的表达式和平衡点的稳定性
系统
（2）总存在无病平衡点E0=（1-m，0，0）。由传染病动力学知识，令基本再生数R0=■。基本再生数表示在发病初期，当所有人均为易感者时，一个病人在患病期间内所传染的人数[3-4]。关于无病平衡点的稳定性我们有以下的结论：
定理一 对于系统
（2），当R01时，无病平衡点E0是局部渐近稳定的。
证明 该系统在无病平衡点E0的雅可比矩阵为
■，
定理二 对于系统
（2），当R01时，存在唯一的地方病平衡点E*。
证明 求模型
（2）的地方病平衡点E■（S■■，I■■，R■■，I■■）。将S■■，I■■，R■■，I■■代入模型
（2）中得
（1-m）b（1-I■■）+pb′I■■-β1S■■I■■-bS■■=0
（3）
β1S■■I■■-（pb′+γ）I■■=0
（4）
β2（1-I■■）I■■-μI■■=0
（5）
由
（5）可得■，由
（4）可得■，将S■■，I■■代入
（3）可以得到 -β1■I■■+[pb′-b（1-m）]■-■+b（1-m）=0
即■
定理三 对于系统
（2），当R01时，地方病平衡点E*是局部渐近稳定的。
证明 该系统在地方病平衡点E*的雅可比矩阵为
■
对应的特征方程为
3 数值模拟
图1
图2
图3
图4
【参考文献】
[2]B.Buonomo， Cr.Vargas-De-Len. Global Analysis of a Host-Vector Model with Infectious force in Latent and Infected Period[J]. Acta Analysis Functionalis Applicata， 2012，14（04）.
[4]H.W. Hethcote. The mathematics of infectious diseases[J]. SIAM Rev， 2000，42：599-653.
[5]Z.E. Ma， Y.C. Zhou， W.D. Wen. The mathematical modeling and research on dynamics of infectious diseases[M]. Beijing：Science Press， 2004.