The aim of this paper is to develop two delayed SEIR epidemic models with nonlinear incidence rate, continuous treatment, and impulsive vaccination for a class of epidemic with latent period and vertical transition. For continuous treatment, we obtain a basic reproductive number ℜ0 and prove the global stability by using the Lyapunov functional method. We obtain two thresholds ℜ* and ℜ* for impulsive vaccination and prove that if ℜ*<1, then the disease-free periodic solution is globally attractive and if ℜ*>1, then the disease is permanent by using the comparison theorem of impulsive differential equation. Numerical simulations indicate that pulse vaccination strategy or a longer latent period will make the population size infected by a disease decrease.
1. Introduction
Mathematical models describing the population dynamics of infectious diseases have been playing an important role in understanding epidemiological patterns and disease control. Researchers have studied the epidemic models by ordinary differential equations [1–3] and the references cited therein. A customarily epidemic model is susceptible, infectious, and recovered model (SIR for short) [4–7]. But in real life, many diseases have a period of incubation time inside the hosts before the hosts become infectious; if we include incubation period of the hosts, the model is described as SEIR model. As tuberculosis (TB), measles and so on, a susceptible individual becomes exposed (infected but not infective) by adequate contact with an infectious individual. SEIR infections disease model has been studied by many authors for its important biological meaning [8–13]. In [13], the authors considered the following delayed SEIR epidemic model:
(1)S′(t)=Λ-μS(t)-βS(t)I(t)1+αI(t)+δR(t),E′(t)=βS(t)I(t)1+αI(t)-βe-μτS(t-τ)I(t-τ)1+αI(t-τ)-μE,I′(t)=βe-μτS(t-τ)I(t-τ)1+αI(t-τ)-(μ+γ)I,R′(t)=γI-(μ+δ)R(t),
where S(t), E(t), I(t), and R(t) represent the number of individuals who are susceptible, exposed, infected, and removed, respectively. The parameters Λ, β, γ, and μ are positive constants, and here Λ is the constant recruitment rate into the population, β is the contact rate, μ is the birth and death rate, γ is the removal rate. τ>0 represents a time delay describing the latent period of the disease and the term (βe-μτS(t-τ)I(t-τ))/(1+αI(t-τ)) represents the individuals surviving in the latent period τ and becoming infective at time t. The sufficient conditions are obtained for the global asymptotic stability of the endemic equilibrium.
In the study of epidemic model, the spread of an infectious disease is a crucial issue, which depends on both the population behavior and the infectivity of the disease. These two aspects are captured in the incidence rate of a disease. In many epidemiological models, the incidence rate is described as mass action incidence with bilinear interactions given by βSI, where β is the probability of transmission per contact and S and I represent the susceptible and infected populations, respectively. This contact law is more appropriate for a few of infected individuals; when the size of infected individuals is increasing, the underlying assumption of homogeneous mixing may not be valid. In fact, with the increase of infected populations, the susceptible individual will take measure to prevent unbounded contact rates. In [14], Anderson and May proposed a saturated incidence rate of the form βSI/(1+αS) in which βSI measures the infection force of the disease and 1/(1+αS) measures the inhibition effect from the behavioral change of the susceptible individuals. The same as the nonlinear incidence rates of the form kIpSq were investigated by Liu et al. [15, 16].
In real life, some diseases may be transferred through horizontal transmission and vertical transmission (disease is the passing of an infection to offspring of infected parents). The offspring of infected parents may already be infected with the disease at birth, so many infections in nature transmit through both horizontal and vertical modes, such as tuberculosis (TB), rubella, hepatitis B, and AIDS [17–21].
Vaccination and treatment are important strategy for the elimination of infectious diseases. Recently, pulse vaccination has been confirmed as an effective method to prevent the spread of the disease [22–24]. Theoretical results show that the pulse vaccination strategy can be distinguished from the conventional strategies in leading to disease eradication at relatively low values of vaccination [25]. The study of vaccination, treatment, and associated behavioral changes related to disease transmission has been the subject of intense theoretical analysis.
The literature on SEIR model with nonlinear incidence, constant infectious period, impulsive vaccination, dealing with the analysis of disease that is vertically and horizontally transmitted is not extensive [17, 19]. But, in fact, under the situation of disease with vertical transmission, the continuous treatment should be considered for the infected, and impulsive vaccination to the susceptible, newborns of the susceptible, exposed and the removed, and newborns of infected which not be vertical infected.
Motivated by the literature above, we introduce delay epidemic models with nonlinear incidence rates of the form βSpI, and we also considered the constant latency period and vertically and horizontally in (2). The purpose of this paper is to study the nonlinear dynamics of system, and we consider two different strategies to the model which are constant treatment and pulse vaccination to the newborns and susceptible,
(2)S′(t)=bm(S+R+E)-βSpI-bS+q′δI,E′(t)=βSpI-βe-bτSp(t-τ)I(t-τ)-bE,I′(t)=βe-bτSp(t-τ)I(t-τ)+qδI-δI-γI-dTI,R′(t)=γI-bR+bm′(S+R+E)+dTI.
The basic assumptions are as follows.
The total population size at time t (day) is denoted by N=S+E+I+R. For N˙=0, this shows that the total population has a constant size. Without loss of generality, we assume in this paper N=1. The newborns of S,E, and R are susceptible individuals, and the newborns of I who are not vertically infected are also susceptible individuals.
The positive constant b (per day) denotes the death rate and birth rate of susceptible, exposed, and recovered individuals. The positive constant δ (per day) denotes the death rate and birth rate of infective individuals. The positive constant γ (per day) is the natural recovery rate of infective individuals. The positive constant q (q≤1) (per day) is the vertical transmission rate, and note q′=1-q, q′<q (per day), and then 0<q′<1. Fraction m′ of all newborns with mothers in the susceptible, exposed, and recovered classes are vaccinated and appeared in the recovered class, while the remaining fraction, m=1-m′, appears in the susceptible class; suppose bm>q′δ.(d/T) is the proportion of those cured successfully.
The incidence rate is described by a nonlinear function βSpI where β (per day) is a positive constant describing the infection rate. τ>0 is the length of the latent period and the term βe-bτSp(t-τ)I(t-τ) reflects the fact that an individual is surviving in the latent period τ and becoming infective at time t.
The remaining part of this paper is organized as follows. In Section 2, we investigate the global stability of the endemic equilibrium of (2) by using Rouches theorem and Lyapunov-LaSalle type theorem. The global asymptotic stability of disease-free periodic solution and the conditions for the permanence of the disease by comparison techniques are described in Section 3. Numerical simulations are presented in Section 4. In Section 5, we conclude this paper with some remarks.
2. Continuous Treatment Strategy of the SEIR Model
In this section, we consider a continuous treatment of SEIR model with constant latent period and nonlinear incidence rate. By using S+R+E=1-I, notice that first and third equations of system (2) do not contain the variables E and R; therefore, system (2) is equivalent to the following 2-dimensional system:
(3)S′(t)=bm(1-I)-βSpI-bS+q′δI,I′(t)=βe-bτSp(t-τ)I(t-τ)+qδI-δI-γI-dTI.
2.1. Disease-Free Equilibrium and Its Stability
First, we define
(4)ℜ0=βmpe-bτδ-qδ+γ+(d/T).
Theorem 1.
If ℜ0<1, the disease-free equilibrium E0(m,0) of system (3) is locally asymptotically stable for all τ≥0; if ℜ0>1, the disease-free equilibrium E0(m,0) is unstable.
Proof.
Steady states of system satisfy the following system of equations:
(5)bm(1-I)-βSpI-bS+q′δI=0,βe-bτSpI+qδI-δI-γI-dTI=0.
Obviously, E0(m,0) is the disease-free equilibrium of (3). In order to analyze the behavior of the system (3) near E0, we linearize the system about the equilibrium point; let S(t)=X(t)+m, I(t)=Y(t)(6)X′(t)=-bX(t)-(βmp+bm-q′δ)Y(t),Y′(t)=-(δ-qδ+γ+dT)Y(t)+βe-bτmpY(t-τ).λ1=-b<0 is one of the eigenvalues of the linearization of system (6) near the steady state E0, and the other eigenvalue λ2 is determined by equation:
(7)λ-βmpe-(b+λ)τ+δ-qδ+γ+dT=0.
Let
(8)f(λ)=λ-βmpe-(b+λ)τ+δ-qδ+γ+dT;
if ℜ0>1, it is easy to show that, for λ real,
(9)f(0)=(δ-qδ+γ+dT)(1-ℜ0)<0,limλ→+∞f(λ)=+∞;
hence, f(λ)=0 has a positive real root. Therefore, if ℜ0>1, the disease-free equilibrium E0(m,0) is unstable.
If ℜ0<1, we prove that the disease-free equilibrium E0(m,0) is locally stable. Otherwise, Reλ≥0. We note that
(10)Reλ=(δ-qδ+γ+dT)(ℜ0e-Reλτcos(lmλτ)-1)≤(δ-qδ+γ+dT)(ℜ0-1),
a contradiction. Hence, the disease-free equilibrium E0 is locally asymptotically stable if ℜ0<1.
Theorem 2.
If ℜ0<1, the disease-free equilibrium E0(m,0) of system (3) is globally asymptotically stable for all τ≥0.
To proof the global stability of the disease-free equilibrium E0(m,0), we choose Lyapunov function
(11)V(t)=X(t)+Y(t)+βe-bτmp∫t-τtY(ξ)dξ,
and it is easy to prove V′(t)<0, limt→∞V(t)=0; it follows that limt→∞X(t)=0, limt→∞Y(t)=0.
2.2. Endemic Equilibrium and Its Stability
If ℜ0>1, then system (3) has a unique positive equilibrium E1(S*,I*), where
(12)S*=((δ-qδ+γ+(d/T))ebτβ)1/p,I*=b(bm-q′δ)+(δ-qδ+γ+(d/T))ebτ×(m-((δ-qδ+γ+(d/T))ebτβ)1/p).
Theorem 3.
If ℜ0>1, conditions (17) and (22) are satisfied, then for τ≥0 the endemic equilibrium E1(S*,I*) of system (3) is locally asymptotically stable.
Proof.
Let S(t)=X(t)+S*, I(t)=Y(t)+I*; the linearized system is obtained
(13)X′(t)=-(b+βpI*S*(p-1))X(t)-(βS*p+bm-q′δ)Y(t),Y′(t)=βpe-bτI*S*(p-1)X(t-τ)+βe-bτS*pY(t-τ)-(δ-qδ+γ+dT)Y(t).
From the linearized system we obtain the characteristic equation
(14)λ2+pλ+r+(gλ+q)e-λτ=0,
where
(15)p=δ-qδ+γ+dT+βpS*(p-1)I*+b,r=(b+βpI*S*(p-1))(δ-qδ+γ+dT),g=-βS*pe-bτ,q=-(b+βpI*S*(p-1))βe-bτS*p+(βS*p+bm-q′δ)βpe-bτI*S*(p-1).
For τ=0 the characteristic equation becomes
(16)λ2+(p+g)λ+(r+q)=0,
and we can see that both roots are negative and real if and only if
(17)p+g>0,r+q>0.
Now for τ≠0, if λ=ωi is a root of (16), we have
(18)-ω2+qe-ωτi+pωi+r+gωe-ωτi=0.
Separating the real and imaginary parts, we have
(19)r-ω2+gωsin(ωτ)+qcosωτ=0,pω+gωcosωτ-qsinωτ=0.
Adding both equations and regrouping by powers of ω, we obtain the following fourth degree polynomial
(20)ω4+(p2-g2-2r)ω2+r2-q2=0,
from which we have
(21)ω2=g2-p2+2r±(g2-p2+2r)2-4(r2-q2)2.
It follows that if
(22)p2-g2-2r>0,r2-q2>0,
are satisfied, (20) does not have positive solutions, and the characteristic equation (14) does not have purely imaginary roots. Inequalities in (17) and (22) guarantee that all roots of (14) have no positive roots. According to Rouche’s theorem, Theorem 3 is proved.
Subsequently we discuss the sufficient conditions under which the endemic equilibrium is globally asymptotically stable for the system (3). For S+E+I+R=1, hence, the dynamics of system (3) in the first octant of R4 is equivalent to that of the following system:
(23)S′(t)=bm(1-I)-βSpI-bS+q′δI,E′(t)=βSpI-βe-bτSp(t-τ)I(t-τ)-bE,I′(t)=βe-bτSp(t-τ)I(t-τ)+qδI-δI-γI-dTI.
The initial conditions for system (23) take the form
(24)S(ξ)=φ1(ξ),E(ξ)=φ2(ξ),I(ξ)=φ3(ξ),φ1(0)>0,φ2(0)>0,φ3(0)>0,
where (φ1(ξ),φ2(ξ),φ3(ξ))∈C([-τ,0],R+03), the space of continuous functions mapping the interval [-τ,0] into R+03, where R+03={(x1,x2,x3):xi≥0,i=1,2,3}.
For continuity of the initial conditions, we require
(25)E(0)=∫-τ0βebξS(ξ)pI(ξ)dξ.
It is well known by the fundamental theory of functional differential equations [26], the system (23) has a unique solution (S(t),E(t),I(t)) satisfying the initial conditions. It is easy to show that all solutions of system (23) with initial conditions are defined on [0,+∞) and remain positive for all t≥0.
Lemma 4 (see [25]).
Let the initial condition be S(ξ)=S(0)>0, E(ξ)=E(0)>0 and I(0)>0, for all ξ∈[-τ,0). Then S(t)≤max{1,S(0)+E(0)+I(0)}=M.
Theorem 5.
Let the initial condition be S(ξ)=S(0)>0, E(ξ)=E(0)>0 and I(0)>0, for all ξ∈[-τ,0). Further suppose ℜ0>1; then for any infectious period τ satisfying
(26)τ>max{βMp-1(I*p+2Mp+3M)2(δ-qδ+γ+dT)+ρ(bm-q′δ)1blnβMp-1(ρM+2pI*+3ρpI*)βpI*Mp-1-4bρ-ρ(bm-q′δ),1blnβMp-1(ρ+1)(I*p+M)4ρb+2b-βpI*Mp-1,1blnβMp-1(I*p+2Mp+3M)2(δ-qδ+γ+dT)+ρ(bm-q′δ)},
where M=max{1,S(0)+E(0)+I(0)}, the endemic equilibrium is globally asymptotically stable.
Proof.
Let S(t)=X(t)+S*, E(t)=Y(t)+E*, I(t)=Z(t)+I*; the linearized system is:
(27)X′(t)=-(b+βpI*S*(p-1))X(t)+(q′δ-βS*p-bm)Z(t),Y′(t)=βpI*S*(p-1)X(t)+βS*pZ(t)-βpe-bτI*S*(p-1)X(t-τ)-βe-bτS*pZ(t-τ)-bY(t),Z′(t)=βpe-bτI*S*(p-1)X(t-τ)+βe-bτS*pZ(t-τ)-(δ-qδ+γ+dT)Z(t).
The trivial solution of system (27) is globally asymptotically stable and is equivalent to the fact that the endemic equilibrium (S*,E*,I*) of system (23) is globally asymptotically stable. We will employ Lyapunov functional technique to prove it.
Now let us introduce the following functions:
(28)V1(t)=12ρ(X(t)+Y(t))2+12(Y2(t)+Z2(t)),V2(t)=(ρ+1)βpe-bτI*M(p-1)∫t-τtX2(ξ)dξ+(ρ+1)βe-bτMp∫t-τtZ2(ξ)dξ,
where ρ>0 is an arbitrary real constant. Choosing ρ=βpS*p/(bm-q′δ), the derivative of V1(t) is
(29)V1′(t)=ρ(X(t)+Y(t))(X′(t)+Y′(t))+Y(t)Y′(t)+Z(t)Z′(t)=ρ(X(t)+Y(t))×[-bX(t)+(q′δ-bm)Z(t)-bY(t)-βpe-bτI*S*(p-1)X(t-τ)-βe-bτS*pZ(t-τ)(q′δ-bm)Z(t)]+Y(t)[βpI*S*(p-1)X(t)+βS*pZ(t)ccccccc-βpe-bτI*S*(p-1)X(t-τ)ccccccc-βe-bτS*pZ(t-τ)-bY(t)]+Z(t)[dTβpe-bτI*S*(p-1)X(t-τ)ccccccc+βe-bτS*pZ(t-τ)ccccccc-(δ-qδ+γ+dT)Z(t)]=-bρX2(t)-(bρ+b)Y2(t)-[δ-qδ+γ+dT]Z2(t)+[βpI*S*(p-1)-2ρb]X(t)Y(t)+ρ(q′δ-bm)Z(t)X(t)-ρβpe-bτI*S*(p-1)X(t)X(t-τ)-βρe-bτS*pX(t)Z(t-τ)-βρe-bτS*pY(t)Z(t-τ)-ρβpe-bτI*S*(p-1)Y(t)X(t-τ)-βpe-bτI*S*(p-1)Y(t)X(t-τ)-βe-bτS*pY(t)Z(t-τ)+βpe-bτI*S*(p-1)X(t-τ)Z(t)+βe-bτS*pZ(t-τ)Z(t),
and applying Cauchy-Chwartz inequality to all product terms, we obtain the following expression:
(30)V1′(t)≤-[2bρ-12βpI*S*(p-1)+12ρ(bm-q′δ)-12βρe-bτMp-1(I*p+M)]X2(t)-[2ρb+b-12pβI*Mp-1-12(ρ+1)βe-bτMp-1(I*p+M)]Y2(t)-[δ-qδ+γ+dT+12ρ(bm-q′δ)-12βe-bτMp-1(I*p+M)]Z2(t)+(ρ+1)βpe-bτI*M(p-1)X2(t-τ)+(ρ+1)βe-bτMpZ2(t-τ).
We choose Lyapunov function to be the form V(t)=V1(t)+V2(t), and we get
(31)V′(t)=V1′(t)+(ρ+1)βpe-bτI*M(p-1)×(X2(t)-X2(t-τ))+(ρ+1)βe-bτMp×(Z2(t)-Z2(t-τ)).
Substituting this in the inequality for V1(t), we get
(32)V′(t)≤-[2bρ+12ρ(bm-q′δ)-12βpI*M*(p-1)-12βe-bτMp-1(ρM+2pI*+3ρpI*)]X2(t)-[2ρb+b-12pβI*Mp-1-12(ρ+1)βe-bτMp-1(I*p+M)]Y2(t)-[δ-qδ+γ+dT+12ρ(bm-q′δ)-12βe-bτMp-1(I*p+2Mρ+3M)]Z2(t).
The right-hand expression of the above inequality is always negative provided that (26) holds. A direct application of the Lyapunov-LaSalle type theorem shows that limt→∞X(t)→0, limt→∞Y(t)→0, limt→∞Z(t)→0. The proof is complete.
3. Continuous Treatment and Pulse Vaccination Strategies
When continuous treatment and pulse vaccination strategies are included in the SEIR epidemic model with the nonlinear infectious force and vertical transmission, it can be written as follows:
(33)S′(t)=bm(S+R+E)-βSpI-bS+q′δI,E′(t)=βSpI-βe-bτSp(t-τ)I(t-τ)-bE,I′(t)=βe-bτSp(t-τ)I(t-τ)+qδI-δI-γI-dTI,R′(t)=γI-bR+bm′(S+R+E)+dTI.t≠kT,k∈Z+,S(t+)=(1-θ)S(t),E(t+)=E(t),I(t+)=I(t),R(t+)=R(t)+θS(t),t=kT,k∈Z+,
and S(t), E(t), I(t), and R(t) are the number of susceptible, exposed, infectious, and recovered at time t, respectively. θ is the proportion of those vaccinated successfully at kT, which is called pulse vaccination rate.
We also consider the following reduced systems:
(34)S′(t)=bm(1-I)-βSpI-bS+q′δI,I′(t)=βe-bτSp(t-τ)I(t-τ)+qδI-δI-γI-dTI,t≠kT,k∈Z+,S(t+)=(1-θ)S(t),I(t+)=I(t),t=kT,k∈Z+.
Let Ω be the following subset of R+2, Ω={(S,I)∈R+2|S≥0,I≥0,S+I≤1}. From biological considerations, we discuss system (34) in the closed set Ω. It can be verified that Ω is positively invariant with respect to system (34).
We first state two important lemmas which are useful in our following discussions.
Lemma 6 (see [12]).
Consider the following impulsive differential equation
(35)u′(t)=a-bu(t),t≠kT,u(t+)=(1-θ)u(t),t=kT,
where a>0, b>0, 0<θ<1. Then the above system has a unique positive periodic solution given by
(36)u*(t)=ab+(u¯-ab)e-b(t-kT),kT<t≤(k+1)T,
which is globally asymptotically stable, where
(37)u¯=ab(1-θ)(1-e-bT)1-(1-θ)e-bT.
Lemma 7 (see [27, 28]).
Consider the following equation:
(38)u′(t)=a1u(t-τ)-a2u(t),
where a1,a2,τ>0, u(t)>0 for -τ≤t≤0. One has
if a1<a2, then limt→∞u(t)=0;
if a1>a2, then limt→∞u(t)=+∞.
3.1. Global Stability of the Disease-Free Periodic Solution
Now we will prove the disease-free periodic solution (S(t),0) is global attractively. We first demonstrate the existence of the disease-free periodic solution, in which infectious individuals are entirely absent from the population permanently, that is, I(t)≡0 for all t>0. Under this condition, the growth of susceptible individuals must satisfy
(39)S′(t)=bm-bS,t≠kT,k∈Z+,S(t+)=(1-θ)S(t),t=kT,k∈Z+.
By Lemma 6, we obtain the periodic solution of system (39),
(40)S*(t)=m-mθe-b(t-nT)1-(1-θ)e-bT,kT<t≤(k+1)T,
and this solution is globally asymptotically stable. Hence, the system (34) has a disease-free periodic solution (S*(t),0).
Theorem 8.
Let (S(t),I(t)) be any solution of (34); then the disease-free periodic solution (S*(t),0) is globally asymptotically stable provided that
(41)ℜ*=mpe-bτβδ-qδ+γ+(d/T)(1-e-bT1-(1-θ)e-bT)p<1.
Proof.
Since ℜ*<1, we can choose ɛ0>0 sufficiently small such that
(42)e-bτβ(m(1-e-bT)1-(1-θ)e-bT+ɛ0)p<δ-qδ+γ+dT.
From the first equation of system (34), we have S′(t)<bm-bS, and then we consider the following comparison system with pulse:
(43)w′(t)=bm-bw(t),t≠kT,w(t+)=(1-θ)w(t),t=kT,w(0+)=S(0+).
In view of Lemma 6, we obtain S(t)≤w(t) and
(44)limt→∞w(t)=m-mθe-b(t-nT)1-(1-θ)e-bT=S*(t),kT<t≤(k+1)T.
There exists an integer k1, such that
(45)S(t)≤w(t)<S*(t)+ɛ0.
That is,
(46)S(t)<S*(t)+ɛ0≤m(1-e-bT)1-(1-θ)e-bT+ɛ0:=S¯,kT<t≤(k+1)T,k>k1.
Furthermore, from the second equation, we have
(47)I′(t)≤βe-bτS¯pI(t-τ)-(δ-qδ+γ+dT)I(t),ccccccccccccccccccccccccccct≠kT,k>k1,
and then we consider the following comparison equation:
(48)w1′(t)=βe-bτS¯pI(t-τ)-(δ-qδ+γ+dT)I(t),cccccccccccccccccccccccccccct≠kT,k>k1;
then, from (41) we have βe-bτS¯p<δ-pδ+γ+(d/T). In view of Lemma 7, we have w1′(t)<0, limk→∞I(t)=0. So there must exist an integer k2>k1, such that I(t)<ɛ1 for all t>k2T.
When t>k2T, from the first equation of system (34), we have
(49)S′(t)>(bm-(bm-q′δ)ɛ1)-(βS¯p-1ɛ1+b)S(t).
Consider the following comparison impulsive differential equation for all t>k2T(50)w2′(t)=(bm-(bm-q′δ)ɛ1)-(βS¯p-1ɛ1+b)w2(t),t≠kT,w2(t+)=(1-θ)w2(t),t=kT,w2(0+)=S(0+).
By Lemma 6, we have the unique periodic solution of system (50) given by
(51)w2*=bm-(bm-q′δ)ɛ1βS¯p-1ɛ1+b(1-θe-(ɛ1βS¯p-1+b)(t-nT)1-(1-θ)e-(ɛ1βS¯p-1+b)T),ccccccccccccccccccccccccccccccccccckT<t≤(k+1)T.
By the comparison theorem, there exists an integer k3>k2 such that
(52)S(t)>w2(t)>w2*-ɛ1,kT<t≤(k+1)T.
Because ɛ0 and ɛ1 are sufficiently small, it follows from (46) and (52) that limt→∞S(t)=S*(t). Therefore, the disease-free solution (S*(t),0) of system (34) is globally attractive. The proof is completed.
(i) If βmp≤ebτ(δ-qδ+γ+(d/T)), then the infection-free periodic solution is globally attractive.
(ii) If βmp>ebτ(δ-qδ+γ+(d/T)), then the infection-free periodic solution is globally attractive provided that θ>θ*.
Corollary 10.
(i) If βmp((1-e-bT)/(1-(1-θ)e-bT))p≤ebτ(δ-qδ+γ+(d/T)), then the infection-free periodic solution is globally attractive.
(ii) If βmp((1-e-bT)/(1-(1-θ)e-bT))p>ebτ(δ-qδ+γ+(d/T)), then the disease will be endemic and system (34) is permanent provided that τ>τ*.
Theorem 8 determines the global attractivity of (34) in Ω for the case ℜ*<1. Its epidemiological implication is that the infectious population vanishes in time so the disease dies out. Corollaries 9 and 10 imply that the disease will disappear if the vaccination rate or the length of latent period of the disease is large enough.
3.2. Persistent
In this section we say the disease is endemic if the infectious population persists above a certain positive level for sufficiently large time. The endemicity of the disease can be well captured and studied through the notion of uniform persistence.
Definition 11.
System (34) is said to be uniformly persistent if there exist positive constants Mi≥mi, i=1,2, (both are independent of the initial values), such that every solution (S(t),I(t)) with positive initial conditions of system (34) satisfies
(54)m1≤S(t)≤M1,m2≤I(t)≤M2.
If ℜ*>1, then there is a positive constant mI such that each positive solution (S(t),I(t)) of system (34) satisfies I(t)≥mI for all t sufficiently large.
Proof.
Let (S(t),I(t)) be any solution with initial values of system (34), and then it is obvious that S(t)≤1, I(t)≤1 for all t>0. We are left to prove there exist positive constants mS, mI and t0 (t0 is sufficiently large) such that S(t)≥mS, I(t)≥mI for all t>t0.
Firstly, from the first equation of system (34), we have
(56)S′(t)>q′δ-(β+b)S.
Consider the following comparison equations:
(57)X′(t)=q′δ-(β+b)X(t),t≠kT,X(t+)=(1-θ)X(t),t=kT.
By Lemma 6 and the comparison theorem [29], we know that for any sufficiently small ɛ>0, there exists a t0 (t0 is sufficiently large) such that
(58)S(t)≥X(t)>X*(t)-ɛ≥q′δβ+b((1-θ)(1-e-(β+b)T)1-(1-θ)e-(β+b)T)-ɛ=mS>0.
Now, we will prove that there exist mI>0 and a sufficiently large t0 such that I(t)≥mI for all t>t0. Since the proof is rather long, it will be convenient to divide it into two steps.
Step 1. Since ℜ*>1, there exist mI*>0, ɛ¯>0 sufficiently small such that
(59)βηpe-bτ-(δ+γ+dT-qδ)>0,
where η=(q′δ/(βmI*+b))((1-θ)(1-e-bT)/(1-(1-θ)e-bT))-ɛ¯.
We claim that for any t0>0, it is impossible that I(t)<mI* for all t≥t0. Suppose that the claim is not valid. There exists a t0>0 such that I(t)<mI* for all t≥t0, and then follows from the first equation of system (34) that for t≥t0,
(60)S′(t)>q′δ-(βmI*+b)S(t).
Consider the comparison impulsive system for t≥t0,
(61)X1′(t)=q′δ-(βmI*+b)X(t),t≠kT,X1(t+)=(1-θ)X1(t),t=kT.
According to Lemma 6, there exists T1>t0 such that
(62)S(t)>X1*(t)-ɛ¯≥q′δβmI*+b((1-θ)(1-e-β(mI*+b)T)1-(1-θ)e-β(mI*+b)T)-ɛ¯=η,
for all T1>t0. The second equation of system (34) can be translated into the following form:
(63)I′(t)=(βSp(t)e-bτ-δ+qδ-γ-dT)I(t)-βe-bτddt∫t-τtSp(ξ)I(ξ)dξ.
Define a function V(t) such that
(64)V(t)=I(t)+βe-bτ∫t-τtSp(ξ)I(ξ)dξ;
then the derivative of V(t) along the solution of system (34) is
(65)V′(t)=(βSp(t)e-bτ-δ+qδ-γ-dT)I(t)=(δ-qδ+γ+dT)(βSp(t)e-bτδ-qδ+γ+(d/T)-1)I(t),ccccccccccccccccccccccccccccccccccccccccct>T1.
From (55), we obtain V′(t)>0, t>T1, which implies that V(t)→∞, t→∞. This is contrary to V(t)<1+βτe-bτ. Hence, there exists a t1>0 such that I(t1)≥mI*.
Step 2. According to Step 1, for any positive solution (S(t),I(t)) of system (34), we are left to consider two cases. First, if I(t)>mI* for all t>t1, then our aim is obtained. Second I(t) oscillates about mI* for all large t. In this case, setting t*=inft>t1I(t)≤mI*, there are two possible cases for t*.
Define
(66)mI=min{mI*2,q1},q1=mI*e-(δ-qδ+γ+(d/T))τ.
We hope to show that I(t)≥mI for all large t. The conclusion is evident in the first case. For the second case, let t*>0 and ρ>0 satisfy I(t*)=I(t*+ρ)=mI*, and I(t)<I*, S(t)>η for t*<t<t*+ρ. Therefore, it is certain that there exists a g (0<g<τ) such that
(67)I(t)≥mI*2fort*<t<t*+g.
In this case, we will discuss three possible cases in terms of the sizes of g, ρ, and τ.
Case 1. If ρ≤g<τ, then I(t)≥(mI*/2) for t*<t<t*+ρ.
Case 2. If g≤ρ≤τ, then from the second equation of system (34), we can deduce I′(t)>-(δ-qδ+γ+(d/T))I(t) for t∈[t*,t*+τ] and I(t*)=mI*, and it is obvious that I(t)≥q1 for t*<t<t*+g.
Case 3. If g≤T≤ρ, we will consider the following two cases, respectively.
Subcase 3.1. For t*<t<t*+τ, it is easy to obtain I(t)>q1.
Subcase 3.2. For t*+τ<t<t*+ρ, it is easy to obtain I(t)>q1. Then, proceeding exactly as the proof for the above claim, we see that I(t)≥mI for t*+τ<t<t*+ρ. Since this kind of interval [t*,t*+ρ] is chosen in an arbitrary way (we only need t* to be large), we conclude that I(t)≥mI for all large t in the second case. In view of our above discussions, the choices of mI are independent of the positive solution, and we have proved that any positive solution of (34) satisfies I(t)≥mI for all large t. The proof is completed.
Set
(68)θ*=(ebT-1)(1-((ebτ(δ-qδ+γ+(d/T)))/β)1/p)((ebτ(δ-qδ+γ+(d/T)))/β)1/p+(ebT-1),τ*=1blnβebτ(δ-qδ+γ+(d/T))((1-θ)(1-e-bT)1-(1-θ)e-bT)p.
From Theorems 12, we also easily obtain the following results.
Corollary 13.
(i) If βe-bτ>δ-qδ+γ+(d/T), then the disease will be endemic and system (34) is permanent provided that θ<θ*.
(ii) If β((1-θ)(1-e-bT)/(1-(1-θ)e-bT))p>δ-qδ+γ+(d/T), then the disease will be endemic and system (34) is permanent provided that τ<τ*.
4. Numerical Simulations
In this section, we present some numerical simulations to demonstrate our theoretical results established in this paper.
Example 1. Letting b=0.5, m=0.85, m′=0.15, τ=0.5, β=0.7, q′=0.5, q=0.5, γ=0.05, δ=0.1, T=2, p=3, we consider the pulse vaccination strategy (see Figure 1), the result shows that the disease fades away when the proportion of those vaccinated successfully θ=0.5; but the disease will exist everlasting when the proportion of those vaccinated successfully θ=0.1. So, this verifies the results in Corollaries 9 and 13, for the epidemic disease with vertical transition periodical vaccination is an effective method to prevent the disease. Also, it can be seen that with the increase of θ which typically causes oscillation bigger of the susceptible.
Dynamical behavior of the system (34) with b=0.5, m=0.85, m′=0.15, τ=0.5, β=0.7, q′=0.5, q=0.5, γ=0.05, δ=0.1, T=2, and p=3. (a) Time-series of the susceptible population with θ=0.1,0.5 respectively. (b) Time-series of the infective population with θ=0.1,0.5.
Example 2. We use the same parameters as in Example 1 except choosing θ=0.1 and τ=0.1,0.8, respectively. Compare τ=0.1 with τ=0.8 is obviously that the longer of the latent period the lower of the infective number (see Figure 2), and this shows disease with long latent period disadvantage to the spread of the disease.
Dynamical behavior of the system (34) with b=0.5, m=0.85, m′=0.15, θ=0.1, β=0.7, q′=0.5, q=0.5, γ=0.05, δ=0.1, T=2, and p=3. (a) Time-series of the susceptible population with τ=0.1,0.8, respectively. (b) Time-series of the infective population with τ=0.1,0.8.
Example 3. For the same parameters as in Example 1 and take θ∈[0,0.1] and τ∈[0,0.8], we obtain the number of infected individual (see Figure 3), when t=100. With the increase of θ and τ the number of infected individual is decreasing, so these results indicate that it will be helpful to control the disease with vertical transition for bigger θ and τ.
The value of infectious individual with t=100, b=0.5, m=0.85, m′=0.15, β=0.7, q′=0.5, q=0.5, γ=0.05, δ=0.1, T=2, p=3, θ∈[0,0.1], τ∈[0,0.8].
5. Conclusion
In this paper, for a class of epidemic disease with latent period and vertical transition, we present two delayed SEIR epidemic models with nonlinear incidence rate based on the spread characters of the disease (such as tuberculosis). Our model is more approach to the realistic problem which is different from [17, 19]. Moreover, the methods in our model are different from the existing results because more factors are considered. When only considering constant treatment, we obtain basic reproductive number ℜ0 and prove the global stability by using the Lyapunov functional method. For the SEIR model with pulse vaccination we also get the theoretical result, if ℜ*<1, the disease-free periodic solution is globally attractive; and if ℜ*>1, the disease is permanent by using the comparison theorem of impulsive differential equation. By some simulation experiments, it clearly shows that the larger of the proportion of those vaccinated successfully the lower of the infective individuals, and the longer of the latent period the lower of the infective individuals. So these results demonstrate that it will be helpful to control the disease with vertical transition for bigger θ and τ.
Acknowledgment
The work is supported by the National Natural Science Foundation of China (no. 1124319).
AlexanderM. E.MoghadasS. M.Periodicity in an epidemic model with a generalized non-linear incidence2004189175962-s2.0-164234163410.1016/j.mbs.2004.01.003AndersonR. M.MayR. M.1991New York, NY, USAOxford University PressJinY.WangW.XiaoS.An SIRS model with a nonlinear incidence rate2007345148214972-s2.0-3425022063010.1016/j.chaos.2006.04.022KorobeinikovA.Lyapunov functions and global stability for SIR and SIRS epidemiological models with non-linear transmission20066836156262-s2.0-3374659460810.1007/s11538-005-9037-9NieL.TengZ.TorresA.Dynamic analysis of an SIR epidemic model with state dependent pulse vaccination2012134162116292-s2.0-8485655855710.1016/j.nonrwa.2011.11.019MengX.ChenL. S.The dynamics of a new SIR epidemic model concerning pulse vaccination strategy200819725825972-s2.0-3944911983810.1016/j.amc.2007.07.083XiaoD.RuanS.Global analysis of an epidemic model with nonmonotone incidence rate200720824194292-s2.0-3454713426110.1016/j.mbs.2006.09.025d'OnofrioA.Stability properties of pulse vaccination strategy in SEIR epidemic model2002179157722-s2.0-003629119410.1016/S0025-5564(02)00095-0SafiM. A.GarbaS. M.Global stability analysis of SEIR model with holling type II incidence function20122012882605210.1155/2012/826052GaoS.TengZ.XieD.The effects of pulse vaccination on SEIR model with two time delays20082011-22822922-s2.0-4464912222510.1016/j.amc.2007.12.019LiG.JinZ.Global stability of a SEIR epidemic model with infectious force in latent, infected and immune period2005255117711842-s2.0-1844437676710.1016/j.chaos.2004.11.062GaoS.ChenL. S.TengZ.Impulsive vaccination of an SEIRS model with time delay and varying total population size20076927317452-s2.0-3384725838310.1007/s11538-006-9149-xXuR.MaZ.Global stability of a delayed SEIRS epidemic model with saturation incidence rate2010611-22292392-s2.0-7795434925210.1007/s11071-009-9644-3AndersonR. M.MayR. M.Regulation and stability of host-parasite population interactions. I. Regulatory processes197847121924710.2307/3933LiuW. M.HethcoteH. W.LevinS. A.Dynamical behavior of epidemiological models with nonlinear incidence rates19872543593802-s2.0-002307021010.1007/BF00277162LiuW. M.LevinS. A.IwasaY.Influence of nonlinear incidence rates upon the behavior of SIRS epidemiological models19852321872042-s2.0-0022298258HerreraM.Modeling the spread of tuberculosis in semiclosed communities201320131964829110.1155/2013/648291d'OnofrioA.On pulse vaccination strategy in the SIR epidemic model with vertical transmission20051877297322-s2.0-1814441358510.1016/j.aml.2004.05.012LuZ.ChiX.ChenL. S.The effect of constant and pulse vaccination on SIR epidemic model with horizontal and vertical transmission2002369-10103910572-s2.0-0037195889KgosimoreM.LunguE. M.The effects of vaccination and treatment on the spread of HIV/AIDS20041243994172-s2.0-1214426108010.1142/S0218339004001294BoilyM. C.AndersonR. M.Sexual contact patterns between men and women and the spread of HIV-1 in urban centres in Africa1991842212472-s2.0-0026275496SabinA. B.Measles, killer of millions in developing countries: strategy for rapid elimination and continuing control1991711222-s2.0-0026035347de QuadrosC. A.AndrusJ. K.OliveJ.da SilveiraC. M.EikhofR. M.CarrascoP.FitzsimmonsJ. W.PinheiroF. P.Eradication of poliomyelitis: progress in the Americas19911032222292-s2.0-0025905391RamsayM.GayN.MillerE.RushM.WhiteJ.Morgan-CapnerP.BrownD.The epidemiology of measles in England and Wales: rationale for the 1994 national vaccination campaign1994412R141R1462-s2.0-0028771291YongzhenP.ShupingL.ChangguoL.ChenS.The effect of constant and pulse vaccination on an SIR epidemic model with infectious period2011358386638782-s2.0-7995509965810.1016/j.apm.2011.02.012HaleJ.1977Heidelberg, GermanySpringerKuangY.1993NewYork, NY, USAAcademic PressXiaoY. N.ChenL. S.Modeling and analysis of a predator-prey model with disease in the prey2001171159822-s2.0-003501695310.1016/S0025-5564(01)00049-9LakshmikanthamV.BainovD.SimeonovP. S.1989SingaporeWorld Scientific