# SVEIRS: A New Epidemic Disease Model with Time Delays and Impulsive Effects

Authors: Tongqian Zhang; Xinzhu Meng; Tonghua Zhang
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

We first propose a new epidemic disease model governed by system of impulsive delay differential equations. Then, based on theories for impulsive delay differential equations, we skillfully solve the difficulty in analyzing the global dynamical behavior of the model with pulse vaccination and impulsive population input effects at two different periodic moments. We prove the existence and global attractivity of the “infection-free” periodic solution and also the permanence of the model. We then carry out numerical simulations to illustrate our theoretical results, showing us that time delay, pulse vaccination, and pulse population input can exert a significant influence on the dynamics of the system which confirms the availability of pulse vaccination strategy for the practical epidemic prevention. Moreover, it is worth pointing out that we obtained an epidemic control strategy for controlling the number of population input.


---

## Body

## 1. Introduction

In epidemic modeling, susceptible-infectious-recovered type of models is well known [1–18] although such models very often ignore the incubation period in the development of mathematical models for some diseases. However, recent research shows for certain diseases, such as smallpox, rabies, BSE, and some skin diseases, the incubation period has significant effect on the epidemic dynamics so that it is nonnegligible. The incubation period varies greatly from a couple of days (e.g., H1N1 outbreaking worldwide has generally an incubation period of one to seven days) to several years (e.g., AIDS virus sometimes can be several years). When taking the incubation period into account in the development of models, we reach SEIR model, which is short for susceptible, exposed, infectious, and recovered [19–30]. And some researchers used time delay to describe the incubation period; for example, Cooke [31], Beretta and Takeuchi [4], Takeuchi et al. [32], and Ma et al. [5] studied a SIR model with time delay and nonlinear incidence rate β
S
(
t
)
I
(
t
-
τ
). Liu et al. [33, 34] used a nonlinear incidence rate β
S
p
(
t
)
I
q
(
t
), and Meng et al. [35] and Jiang et al. [30], respectively, studied an impulsively vaccinating SIR model with nonlinear incidences β
S
q
(
t
)
I
(
t
-
τ
) and β
S
q
(
t
-
τ
)
I
(
t
-
τ
), which are better to describe the spread process of diseases than linear one.In order to prevent infectious diseases, [36, 37] suggested that vaccination to the susceptible population is an important strategy. The traditional vaccinations are applied to each individual, while impulsive ones are to periodically vaccinate people within certain age groups [7–10, 38]. Some diseases may have a vaccination period after being cured but may cause losing immunity gradually. In this case, people might be infected again. So it is of great significance to investigate epidemic models with time delay and impulsive effects due to the incubation period and vaccination period [26–29]. For some certain regional systems, the immigrations can be periodic impulsive population input because the immigratory population might be susceptible. Certainly two different impulsive effects for periodic vaccination and population input do not usually happen simultaneously. Therefore, motivated by Jiang et al. [30] and Song et al. [19], we built a new mathematical model: susceptible, vaccinated, exposed, infectious, recovered, and susceptible epidemic model with two time delays and two nonlinear incidences with pulse vaccination and a constant periodic population input at two different moments as follows:
(1)
d
S
(
t
)
d
t
=
-
b
S
(
t
)
-
β
S
p
(
t
)
I
(
t
)
+
γ
I
(
t
-
ω
)
e
-
b
ω
,
d
V
(
t
)
d
t
=
-
δ
β
V
q
(
t
)
I
(
t
)
-
γ
1
V
(
t
)
-
b
V
(
t
)
,
d
E
(
t
)
d
t
=
-
b
E
+
β
S
p
(
t
)
I
(
t
)
+
δ
β
V
q
(
t
)
I
(
t
)
-
β
e
-
b
τ
S
p
(
t
)
I
(
t
-
τ
)
-
δ
β
e
-
b
τ
V
q
(
t
)
I
(
t
-
τ
)
,
d
I
(
t
)
d
t
=
β
e
-
b
τ
S
p
(
t
)
I
(
t
-
τ
)
+
δ
β
e
-
b
τ
V
q
(
t
)
I
(
t
-
τ
)
-
(
γ
+
b
+
α
)
I
(
t
)
,
d
R
(
t
)
d
t
=
γ
1
V
(
t
)
+
γ
I
(
t
)
-
b
R
(
t
)
-
γ
I
(
t
-
ω
)
e
-
b
ω
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
Δ
S
(
t
)
=
-
θ
S
(
t
)
,
Δ
V
(
t
)
=
θ
S
(
t
)
,
Δ
E
(
t
)
=
0
,
Δ
I
(
t
)
=
0
,
Δ
R
(
t
)
=
0
,
iiiiiiiiiiiiiiii
t
=
(
n
+
l
-
1
)
T
,
Δ
S
(
t
)
=
μ
,
Δ
V
(
t
)
=
0
,
Δ
E
(
t
)
=
0
,
Δ
I
(
t
)
=
0
,
Δ
R
(
t
)
=
0
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiii
t
=
n
T
.
Here all parameters of system (1) are nonnegative constants. For the significance of parameters in (1), please see literatures Jiang et al. [30] and Song et al. [19]. Terms β
S
p
I and V
q
I are the nonlinear incidence rates, and in our paper we only discuss the case
(2)
1
≤
q
≤
p
.


## 2. Preliminaries

LetN
(
t
)
=
S
(
t
)
+
V
(
t
)
+
E
(
t
)
+
I
(
t
)
+
R
(
t
), and then it is easy to see that N
(
t
) satisfies the following:
(3)
N
′
(
t
)
≤
b
(
1
-
N
(
t
)
)
,
lim
⁡
t
→
∞
sup
⁡
N
(
t
)
≤
1
.
Hence, for time t which is large, we obtain 0
≤
S
(
t
)
+
V
(
t
)
+
I
(
t
)
≤
1. Let ϖ
=
max
⁡
{
τ
,
ω
} and C
+
=
{
φ
=
(
φ
1
(
s
)
,
…
,
φ
5
(
s
)
)
∈
C
:
φ
i
(
0
)
>
0
}; here φ
i
(
s
)
>
0 is bounded function on interval [
-
ϖ
,
0
]. Since variable R
(
t
) only appears in the fifth equation, system (1) can be further reduced as
(4)
d
S
(
t
)
d
t
=
-
b
S
(
t
)
-
β
S
p
(
t
)
I
(
t
)
+
γ
I
(
t
-
ω
)
e
-
b
ω
,
d
V
(
t
)
d
t
=
-
δ
β
V
q
(
t
)
I
(
t
)
-
γ
1
V
(
t
)
-
b
V
(
t
)
,
d
I
(
t
)
d
t
=
β
e
-
b
τ
S
p
(
t
)
I
(
t
-
τ
)
+
δ
β
e
-
b
τ
V
q
(
t
)
I
(
t
-
τ
)
-
(
γ
+
b
+
α
)
I
(
t
)
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
Δ
S
(
t
)
=
-
θ
S
(
t
)
,
Δ
V
(
t
)
=
θ
S
(
t
)
,
Δ
I
(
t
)
=
0
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
t
=
(
n
+
l
-
1
)
T
,
Δ
S
(
t
)
=
μ
,
Δ
V
(
t
)
=
0
,
Δ
I
(
t
)
=
0
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
t
=
n
T
,
with the initial conditions
(5)
(
φ
1
(
s
)
,
φ
2
(
s
)
,
φ
4
(
s
)
)
∈
C
+
,
φ
i
(
0
)
>
0
,
i
=
1,2
,
4
.Lemma 1 (see [39, 40]).
For the following impulse differential inequalities(6)
s
′
(
t
)
≤
(
≥
)
q
(
t
)
s
(
t
)
+
r
(
t
)
,
t
≠
t
k
,
s
(
t
k
+
)
≤
(
≥
)
b
k
s
(
t
k
)
+
p
k
,
t
=
t
k
,
k
∈
N
,
where q
(
t
), r
(
t
)
∈
C
(
R
+
,
R
), b
k
≥
0, and p
k are constants.
Assume the following:(
A
0
)
the sequence{
t
k
} satisfies 0
≤
t
0
<
t
1
<
t
2
<
⋯, with lim
⁡
t
→
∞
t
k
=
∞;
(
A
1
)
w
∈
P
C
′
(
R
+
,
R
) and s
(
t
) is left-continuous at t
k
,
k
∈
N.Then
(7)
s
(
t
)
≤
(
≥
)
s
(
t
0
)
∏
t
0
<
t
k
<
t
‍
b
k
exp
⁡
(
∫
t
0
t
‍
q
(
u
)
d
u
)
+
∑
t
0
<
t
k
<
t
‍
(
∏
t
k
<
t
j
<
t
‍
b
j
exp
⁡
(
∫
t
k
t
‍
q
(
u
)
d
u
)
)
p
k
+
∫
t
0
t
‍
∏
u
<
t
k
<
t
‍
b
k
exp
⁡
(
∫
u
t
‍
q
(
θ
)
d
θ
)
r
(
u
)
d
u
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
t
≥
t
0
.Lemma 2 (see [41]).
For the following delay differential equation(8)
d
z
(
t
)
d
t
=
a
z
(
t
-
θ
)
-
b
z
(
t
)
,
where a, b, and θ are all positive constants and z
(
t
)
>
0 for t
∈
[
-
θ
,
0
], then we have
(9)
lim
⁡
t
→
∞
z
(
t
)
=
{
0
,
i
f
a
<
b
;
+
∞
,
i
f
a
>
b
.Lemma 3 (see [42]).
The following system,(10)
d
x
(
t
)
d
t
=
-
b
x
(
t
)
,
d
y
(
t
)
d
t
=
-
(
a
+
b
)
y
(
t
)
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiii
t
≠
n
T
,
t
≠
(
n
+
l
-
1
)
T
,
Δ
x
(
t
)
=
-
θ
x
(
t
)
,
Δ
y
(
t
)
=
θ
x
(
t
)
,
t
=
(
n
+
l
-
1
)
T
,
Δ
x
(
t
)
=
μ
,
Δ
y
(
t
)
=
0
,
t
=
n
T
,
has a unique positive T-periodic solution:
(11)
x
*
(
t
)
=
{
μ
exp
⁡
(
-
b
(
t
-
(
n
-
1
)
T
)
)
1
-
(
1
-
θ
)
exp
⁡
(
-
b
T
)
,
t
∈
(
(
n
-
1
)
T
,
(
n
+
l
-
1
)
T
]
,
μ
(
1
-
θ
)
exp
⁡
(
-
b
(
t
-
(
n
-
1
)
T
)
)
1
-
(
1
-
θ
)
exp
⁡
(
-
b
T
)
,
t
∈
(
(
n
+
l
-
1
)
T
,
≤
n
T
]
,
y
*
(
t
)
=
μ
θ
exp
⁡
(
-
b
l
T
)
exp
⁡
(
-
(
a
+
b
)
(
t
-
(
n
+
l
-
1
)
T
)
)
(
1
-
exp
⁡
(
-
(
a
+
b
)
T
)
)
(
1
-
(
1
-
θ
)
exp
⁡
(
-
b
T
)
)
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiii
t
∈
(
(
n
+
l
-
1
)
T
,
(
n
+
l
)
T
]
,
and we further have x
(
t
)
→
x
*
(
t
) and y
(
t
)
→
y
*
(
t
) as t
→
+
∞.


## 3. The Existence and Global Attractivity of “Infection-Free” Periodic Solution

### 3.1. Existence

In this section, we are committed to investigate the existence of “infection-free” periodic solution. In this case, we have(12)
I
(
t
)
=
0
,
t
≥
0
.
From systems (4) and (12), we obtain
(13)
d
S
(
t
)
d
t
=
-
b
S
(
t
)
,
d
V
(
t
)
d
t
=
-
(
γ
1
+
b
)
V
(
t
)
,
iiiiiiiiiiiiiii
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
n
∈
N
,
Δ
S
(
t
)
=
-
θ
S
(
t
)
,
Δ
V
(
t
)
=
θ
S
(
t
)
,
iiiiiiiiiiiiiiiii
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
S
(
t
)
=
μ
,
Δ
V
(
t
)
=
0
,
t
=
n
T
,
n
∈
N
.
By Lemma 3, system (13) has a unique positive T-periodic solution:
(14)
S
*
(
t
)
=
{
μ
exp
⁡
(
-
b
(
t
-
(
n
-
1
)
T
)
)
1
-
(
1
-
θ
)
exp
⁡
(
-
b
T
)
,
t
∈
(
(
n
-
1
)
T
,
(
n
+
l
-
1
)
T
]
,
μ
(
1
-
θ
)
exp
⁡
(
-
b
(
t
-
(
n
-
1
)
T
)
)
1
-
(
1
-
θ
)
exp
⁡
(
-
b
T
)
,
t
∈
(
(
n
+
l
-
1
)
T
,
≤
n
T
]
,
V
*
(
t
)
=
μ
θ
exp
⁡
(
-
b
l
T
)
exp
⁡
(
-
(
a
+
b
)
(
t
-
(
n
+
l
-
1
)
T
)
)
(
1
-
exp
⁡
(
-
(
a
+
b
)
T
)
)
(
1
-
(
1
-
θ
)
exp
⁡
(
-
b
T
)
)
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiii
t
∈
(
(
n
+
l
-
1
)
T
,
(
n
+
l
)
T
]
.
Furthermore, we can prove that it is the unique globally asymptotically stable positive periodic solution of system (4). We summarize this conclusion in the following lemma.Lemma 4.
The system (4) has an “infection-free” periodic solution (
S
*
(
t
)
,
V
*
(
t
)
,
0
), for t
∈
(
(
n
+
l
-
1
)
T
,
(
n
+
l
)
T
] and n
∈
N; for any solution (
S
(
t
)
,
V
(
t
)
,
I
(
t
)
) of it, the following holds true:
(15)
S
(
t
)
⟶
S
*
(
t
)
,
V
(
t
)
⟶
V
*
(
t
)
as t
→
∞.This lemma indicates that in between the vaccination the susceptible and vaccinated populations oscillate with periodT in synchronization with the periodic pulse vaccination. Next we prove the global attractivity of such solution.


### 3.2.  Global Attractivity

In this section, we will prove our main result on the global attractivity of the infection-free solution. It is stated in the following theorem.Theorem 5.
The system (4) has a unique infection-free periodic solution (
S
*
(
t
)
,
V
*
(
t
)
,
0
), and when it exists, it is globally attractive if
(16)
R
1
<
1
,
where
(17)
R
1
=
β
e
-
b
τ
(
A
1
p
+
δ
A
2
q
)
γ
+
b
+
α
,
with
(18)
A
1
=
γ
e
-
b
ω
b
+
μ
e
b
T
e
b
T
-
1
,
A
2
=
θ
e
-
b
l
T
(
1
-
θ
)
(
1
-
(
1
-
θ
)
e
-
b
T
)
(
1
-
e
-
(
γ
1
+
b
)
T
)
.Proof.
Let(
S
(
t
)
,
V
(
t
)
,
I
(
t
)
) be a solution of (4) satisfied initial condition (5). Since R
1
<
1, one can choose an ɛ
>
0 small enough such that
(19)
β
e
-
b
τ
(
(
Δ
1
)
p
+
δ
(
Δ
2
)
q
)
-
(
γ
+
b
+
α
)
<
0
,
where
(20)
Δ
1
=
γ
e
-
b
ω
b
+
μ
e
b
T
e
b
T
-
1
+
ɛ
,
Δ
2
=
θ
e
-
b
l
T
(
1
-
θ
)
(
1
-
(
1
-
θ
)
e
-
b
T
)
(
1
-
e
-
(
γ
1
+
b
)
T
)
+
ɛ
.
For n
>
n
1, we have
(21)
iiiiiiiiiiiii
d
S
(
t
)
d
t
≤
b
-
b
S
(
t
)
+
γ
e
-
b
ω
,
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
n
∈
N
,
Δ
S
(
t
)
=
-
θ
S
(
t
)
,
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
S
(
t
)
=
μ
,
t
=
n
T
,
n
∈
N
.
By impulsive differential inequality Lemma 1, we have
(22)
S
(
t
)
≤
S
(
n
1
T
+
)
∏
n
1
T
+
<
n
T
<
t
‍
exp
⁡
(
∫
n
1
T
t
‍
(
-
b
)
d
s
)
+
∑
n
1
T
<
n
T
<
t
‍
(
∏
n
T
<
t
j
<
t
‍
exp
⁡
(
∫
n
T
t
‍
(
-
b
)
d
s
)
)
μ
+
∫
n
1
T
t
‍
∏
s
<
n
T
<
t
‍
exp
⁡
(
∫
s
t
‍
(
-
b
)
d
θ
)
γ
e
-
b
ω
d
s
=
S
1
+
S
2
+
S
3
,
where
(23)
S
1
=
S
(
n
1
T
+
)
∏
n
1
T
+
<
n
T
<
t
‍
exp
⁡
(
∫
n
1
t
/
T
‍
(
-
b
)
d
T
ξ
)
=
S
(
n
1
T
+
)
e
-
b
(
t
-
n
1
T
)
,
S
2
=
∑
n
1
T
<
n
T
<
t
‍
(
∏
n
T
<
t
j
<
t
‍
exp
⁡
(
∫
n
T
t
‍
(
-
b
)
d
s
)
)
μ
=
∑
n
1
T
<
n
T
<
t
‍
(
μ
e
-
b
(
t
-
n
T
)
)
=
μ
e
-
b
(
t
-
n
1
)
T
-
e
-
b
(
t
-
(
n
+
n
1
)
)
T
1
-
e
b
T
,
S
3
=
∫
n
1
T
t
‍
∏
s
<
n
T
<
t
‍
exp
⁡
(
∫
s
t
‍
(
-
b
)
d
θ
)
γ
e
-
b
ω
d
s
=
γ
e
-
b
ω
∫
n
1
T
t
‍
∏
s
<
n
T
<
t
‍
exp
⁡
(
∫
s
t
‍
(
-
b
)
d
θ
)
d
s
=
γ
e
-
b
ω
e
-
b
t
b
∫
n
1
T
t
‍
∏
s
<
n
T
<
t
‍
e
b
s
d
(
b
s
)
=
γ
e
-
b
ω
e
-
b
t
b
(
e
b
t
-
e
b
n
1
T
)
.
Thus
(24)
S
(
t
)
≤
S
1
+
S
2
+
S
3
=
S
(
n
1
T
+
)
e
-
b
(
t
-
n
1
T
)
+
μ
e
-
b
(
t
-
n
1
)
T
-
e
-
b
(
t
-
(
n
+
n
1
)
)
T
1
-
e
b
T
+
γ
e
-
b
ω
e
-
b
t
b
(
e
b
t
-
e
b
n
1
T
)
≤
e
-
b
t
S
(
n
1
T
+
)
e
n
1
b
T
+
μ
e
-
b
(
t
-
n
1
T
)
1
-
e
b
T
+
γ
e
-
b
ω
b
e
-
b
(
t
-
n
1
T
)
+
γ
e
-
b
ω
b
+
μ
e
b
T
e
b
T
-
1
,
and then we have
(25)
lim
⁡
t
→
∞
sup
⁡
S
(
t
)
<
γ
e
-
b
ω
b
+
μ
e
b
T
e
b
T
-
1
.
Thus there exists a positive integer n
2
>
n
1 and constant ɛ
>
0 small enough such that, for all t
>
n
2
T,
(26)
S
(
t
)
≤
γ
e
-
b
ω
b
+
μ
e
b
T
e
b
T
-
1
+
ɛ
=
Δ
1
.
Forn
>
n
1, system (4) yields
(27)
d
V
(
t
)
d
t
≤
-
(
γ
1
+
b
)
V
(
t
)
,
t
≠
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
V
(
t
)
=
θ
S
(
t
)
,
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
.
We obtain the following comparison impulsive differential system:
(28)
d
x
(
t
)
d
t
=
-
(
γ
1
+
b
)
x
(
t
)
,
t
≠
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
x
(
t
)
=
θ
S
(
t
)
,
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
.
By Lemma 3, the system has a periodic solution given by
(29)
x
*
(
t
)
=
θ
e
-
b
l
T
e
-
(
γ
1
+
b
)
(
t
-
(
n
+
l
-
1
)
T
)
(
1
-
θ
)
(
1
-
(
1
-
θ
)
e
-
b
T
)
(
1
-
e
-
(
γ
1
+
b
)
T
)
,
iiiiiiiiiiiiiiiiiiii
t
∈
(
(
n
+
l
-
1
)
T
,
(
n
+
l
)
T
]
,
x
(
0
+
)
=
θ
e
-
b
l
T
(
1
-
θ
)
(
1
-
(
1
-
θ
)
e
-
b
T
)
(
1
-
e
-
(
γ
1
+
b
)
T
)
,
which is globally asymptotically stable.
Now, assume thatx
(
t
) is the solution of system (28) with initial value x
(
0
+
)
=
V
0. Then by Lemma 1, we know there exists a positive integer n such that
(30)
V
(
t
)
<
x
(
t
)
<
x
*
(
t
)
+
ɛ
,
t
∈
(
n
T
,
(
n
+
1
)
T
]
.
Hence,
(31)
V
(
t
)
<
x
(
t
)
<
x
*
(
t
)
+
ɛ
<
θ
e
-
b
l
T
(
1
-
θ
)
(
1
-
(
1
-
θ
)
e
-
b
T
)
(
1
-
e
-
(
γ
1
+
b
)
T
)
+
ɛ
0
=
Δ
2
.
From (27), (31), and the third equation in (4), for t
>
n
2
T
+
τ we have
(32)
d
I
(
t
)
d
t
≤
β
e
-
b
τ
(
Δ
1
p
+
δ
Δ
2
q
)
I
(
t
-
τ
)
-
(
γ
+
b
+
α
)
I
(
t
)
.
Consider the comparison equation:
(33)
d
y
(
t
)
d
t
≤
β
e
-
b
τ
(
Δ
1
p
+
δ
Δ
2
q
)
y
(
t
-
τ
)
-
(
γ
+
b
+
α
)
y
(
t
)
.
From (19), we have
(34)
β
e
-
b
τ
(
Δ
1
p
+
δ
Δ
2
q
)
-
(
γ
+
b
+
α
)
<
0
.
According to Lemma 2, we then obtain
(35)
lim
⁡
t
→
∞
y
(
t
)
=
0
.
Notice the fact that I
(
s
)
=
y
(
s
)
=
ϕ
3
(
s
)
>
0 for all s
∈
[
-
τ
,
0
] and I
(
t
)
≥
0, and the comparison theorem implies I
(
t
)
→
0 as t
→
∞. Without loss of generality, we may assume that 0
<
I
(
t
)
<
ɛ
1 for all t
≥
0. By using the first and second equations in (4), we reach
(3.2)
d
S
(
t
)
d
t
≥
-
b
S
(
t
)
-
β
ɛ
1
S
p
(
t
)
,
d
V
(
t
)
d
t
≥
-
δ
β
ɛ
1
V
q
(
t
)
-
γ
1
V
(
t
)
-
b
V
(
t
)
,
iiiiii
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
n
∈
N
,
Δ
S
(
t
)
=
-
θ
S
(
t
)
,
Δ
V
(
t
)
=
θ
S
(
t
)
,
iiiiiiiiiiiiiiiii
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
S
(
t
)
=
μ
,
Δ
V
(
t
)
=
0
,
t
=
n
T
,
n
∈
N
.
For 1
<
q
<
p, we have
(37)
d
S
(
t
)
d
t
≥
-
b
S
(
t
)
-
β
ɛ
1
S
(
t
)
,
d
V
(
t
)
d
t
≥
-
δ
β
ɛ
1
V
(
t
)
-
γ
1
V
(
t
)
-
b
V
(
t
)
,
iiiiii
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
n
∈
N
,
Δ
S
(
t
)
=
-
θ
S
(
t
)
,
Δ
V
(
t
)
=
θ
S
(
t
)
,
iiiiiiiiiiiiiiiii
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
S
(
t
)
=
μ
,
Δ
V
(
t
)
=
0
,
t
=
n
T
,
n
∈
N
considering the following system:
(38)
d
f
(
t
)
(
t
)
d
t
=
-
b
f
(
t
)
-
β
ɛ
1
f
(
t
)
,
d
g
(
t
)
d
t
=
-
δ
β
ɛ
1
g
(
t
)
-
γ
1
g
(
t
)
-
b
g
(
t
)
,
iiiii
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
n
∈
N
,
Δ
f
(
t
)
=
-
θ
f
(
t
)
,
Δ
g
(
t
)
=
θ
f
(
t
)
,
iiiiiiiiiiiiiiiiii
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
f
(
t
)
=
μ
,
Δ
g
(
t
)
=
0
,
t
=
n
T
,
n
∈
N
.
We obtain
(39)
f
~
(
t
)
=
{
(
b
+
β
ɛ
1
)
μ
(
b
+
β
ɛ
1
)
(
1
-
(
1
-
θ
)
e
-
(
b
+
β
ɛ
1
)
T
)
×
e
-
(
b
+
β
ɛ
1
)
(
t
-
(
n
-
1
)
T
)
,
t
∈
(
(
n
-
1
)
T
,
(
n
+
l
-
1
)
T
]
,
(
b
+
β
ɛ
1
)
μ
(
1
-
θ
)
e
-
(
b
+
β
ɛ
1
)
l
T
(
b
+
β
ɛ
1
)
(
1
-
(
1
-
θ
)
e
-
(
b
+
β
ɛ
1
)
T
)
×
e
-
(
b
+
β
ɛ
1
)
(
t
-
(
n
+
l
-
1
)
T
)
,
t
∈
(
(
n
+
l
-
1
)
T
,
n
T
]
,
g
~
(
t
)
=
(
(
θ
(
b
+
β
ɛ
1
)
μ
e
-
(
b
+
β
ɛ
1
)
l
T
-
β
ɛ
1
θ
(
1
-
e
-
(
b
+
β
ɛ
1
)
T
)
)
iii
(
A
)
×
e
-
(
γ
1
+
b
+
δ
β
ɛ
1
)
(
t
-
(
n
+
l
-
1
)
T
)
)
×
(
(
b
+
β
ɛ
1
)
(
1
-
(
1
-
θ
)
e
-
(
b
+
β
ɛ
1
)
T
)
iiiiii
×
(
1
-
e
-
(
γ
1
+
b
)
T
)
)
-
1
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiii
t
∈
(
(
n
+
l
-
1
)
T
,
(
n
+
l
)
T
]
.
Now by using comparison theorem of impulsive equations, for anyɛ
2
>
0 there exists a T
1
>
0 such that
(40)
S
(
t
)
>
f
~
(
t
)
-
ɛ
2
,
V
(
t
)
>
g
~
(
t
)
-
ɛ
2
,
for t
>
T
1. On the other side, from the first and second equations of (4), we have
(41)
iiiiiiiiiiiii
d
S
(
t
)
d
t
≤
-
b
S
(
t
)
+
γ
ɛ
1
e
-
b
ω
,
iiiiiiiiiiiii
d
V
(
t
)
d
t
≤
-
γ
1
V
(
t
)
-
b
V
(
t
)
,
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
n
∈
N
,
Δ
S
(
t
)
=
-
θ
S
(
t
)
,
Δ
V
(
t
)
=
θ
S
(
t
)
,
iiiiiiiiiiiiiiiii
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
S
(
t
)
=
b
,
Δ
V
(
t
)
=
0
,
t
=
n
T
,
n
∈
N
.
Then we have S
(
t
)
≤
h
~
(
t
), V
(
t
)
≤
g
~
(
t
) and h
~
(
t
)
→
S
*
(
t
), g
~
(
t
)
→
V
*
(
t
), as ɛ
1
→
0, where (
h
~
(
t
)
,
g
~
(
t
)
) is a unique positive periodic solution of
(42)
iiiiiiiiiiiii
d
h
(
t
)
d
t
=
-
b
h
(
t
)
+
γ
ɛ
1
e
-
b
ω
,
iiiiiiiiiiiiii
d
g
(
t
)
d
t
=
-
γ
1
g
(
t
)
-
b
g
(
t
)
,
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
n
∈
N
,
Δ
h
(
t
)
=
-
θ
h
(
t
)
,
Δ
g
(
t
)
=
θ
h
(
t
)
,
iiiiiiiiiiiiiii
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
h
(
t
)
=
μ
,
Δ
g
(
t
)
=
0
,
t
=
n
T
,
n
∈
N
,
from which we have that, for n
T
<
t
≤
(
n
+
1
)
T,
(43)
h
~
(
t
)
=
{
γ
ɛ
1
e
-
b
ω
θ
e
-
b
(
1
-
l
)
T
+
b
μ
b
(
1
-
(
1
-
θ
)
e
-
b
T
)
×
e
-
b
(
t
-
(
n
-
1
)
T
)
-
γ
ɛ
1
e
-
b
ω
b
,
(
n
-
1
)
T
<
t
≤
(
n
+
l
-
1
)
T
,
b
μ
(
1
-
θ
)
e
-
b
l
T
+
γ
ɛ
1
e
-
b
ω
θ
b
(
1
-
(
1
-
θ
)
e
-
b
T
)
×
e
-
b
(
t
-
(
n
+
l
-
1
)
T
)
-
γ
ɛ
1
e
-
b
ω
b
,
(
n
+
l
-
1
)
T
<
t
≤
n
T
,
g
~
(
t
)
=
(
θ
b
μ
e
-
b
l
T
-
γ
ɛ
1
e
-
b
ω
θ
(
1
-
e
-
b
T
)
)
e
-
(
γ
1
+
b
)
(
t
-
(
n
+
l
-
1
)
T
)
b
(
1
-
(
1
-
θ
)
e
-
b
T
)
(
1
-
e
-
(
γ
1
+
b
)
T
)
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
(
n
+
l
-
1
)
T
<
t
≤
(
n
+
l
)
T
.
Applying the comparison theorem again, for any ɛ
2
>
0, there exists a T
2
>
0 such that
(44)
S
(
t
)
<
h
~
(
t
)
-
ɛ
2
,
V
(
t
)
<
g
~
(
t
)
-
ɛ
2
,
for t
>
T
2. Let ɛ
1
→
0, and then from (40) and (44) we have
(45)
S
*
(
t
)
-
ɛ
2
<
S
(
t
)
<
S
*
(
t
)
-
ɛ
2
,
V
*
(
t
)
-
ɛ
2
<
V
(
t
)
<
V
*
(
t
)
-
ɛ
2
,
for t large enough, which implies S
(
t
)
→
S
*
(
t
), V
(
t
)
→
V
*
(
t
) as t
→
∞. This completes the proof.Corollary 6.
Ifτ
>
τ
* or μ
<
μ
*, then the infection-free periodic solution (
S
*
(
t
)
,
V
*
(
t
)
,
0
) is globally attractive, where the critical values are given below:
(46)
τ
*
=
1
b
ln
⁡
β
(
A
1
p
+
δ
A
2
q
)
γ
+
b
+
α
,
μ
*
=
(
1
-
e
-
b
T
)
×
(
γ
+
b
+
α
β
e
b
τ
-
δ
A
2
q
p
-
γ
e
-
b
ω
b
-
μ
e
b
T
e
b
T
-
1
)
.


## 3.1. Existence

In this section, we are committed to investigate the existence of “infection-free” periodic solution. In this case, we have(12)
I
(
t
)
=
0
,
t
≥
0
.
From systems (4) and (12), we obtain
(13)
d
S
(
t
)
d
t
=
-
b
S
(
t
)
,
d
V
(
t
)
d
t
=
-
(
γ
1
+
b
)
V
(
t
)
,
iiiiiiiiiiiiiii
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
n
∈
N
,
Δ
S
(
t
)
=
-
θ
S
(
t
)
,
Δ
V
(
t
)
=
θ
S
(
t
)
,
iiiiiiiiiiiiiiiii
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
S
(
t
)
=
μ
,
Δ
V
(
t
)
=
0
,
t
=
n
T
,
n
∈
N
.
By Lemma 3, system (13) has a unique positive T-periodic solution:
(14)
S
*
(
t
)
=
{
μ
exp
⁡
(
-
b
(
t
-
(
n
-
1
)
T
)
)
1
-
(
1
-
θ
)
exp
⁡
(
-
b
T
)
,
t
∈
(
(
n
-
1
)
T
,
(
n
+
l
-
1
)
T
]
,
μ
(
1
-
θ
)
exp
⁡
(
-
b
(
t
-
(
n
-
1
)
T
)
)
1
-
(
1
-
θ
)
exp
⁡
(
-
b
T
)
,
t
∈
(
(
n
+
l
-
1
)
T
,
≤
n
T
]
,
V
*
(
t
)
=
μ
θ
exp
⁡
(
-
b
l
T
)
exp
⁡
(
-
(
a
+
b
)
(
t
-
(
n
+
l
-
1
)
T
)
)
(
1
-
exp
⁡
(
-
(
a
+
b
)
T
)
)
(
1
-
(
1
-
θ
)
exp
⁡
(
-
b
T
)
)
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiii
t
∈
(
(
n
+
l
-
1
)
T
,
(
n
+
l
)
T
]
.
Furthermore, we can prove that it is the unique globally asymptotically stable positive periodic solution of system (4). We summarize this conclusion in the following lemma.Lemma 4.
The system (4) has an “infection-free” periodic solution (
S
*
(
t
)
,
V
*
(
t
)
,
0
), for t
∈
(
(
n
+
l
-
1
)
T
,
(
n
+
l
)
T
] and n
∈
N; for any solution (
S
(
t
)
,
V
(
t
)
,
I
(
t
)
) of it, the following holds true:
(15)
S
(
t
)
⟶
S
*
(
t
)
,
V
(
t
)
⟶
V
*
(
t
)
as t
→
∞.This lemma indicates that in between the vaccination the susceptible and vaccinated populations oscillate with periodT in synchronization with the periodic pulse vaccination. Next we prove the global attractivity of such solution.


## 3.2.  Global Attractivity

In this section, we will prove our main result on the global attractivity of the infection-free solution. It is stated in the following theorem.Theorem 5.
The system (4) has a unique infection-free periodic solution (
S
*
(
t
)
,
V
*
(
t
)
,
0
), and when it exists, it is globally attractive if
(16)
R
1
<
1
,
where
(17)
R
1
=
β
e
-
b
τ
(
A
1
p
+
δ
A
2
q
)
γ
+
b
+
α
,
with
(18)
A
1
=
γ
e
-
b
ω
b
+
μ
e
b
T
e
b
T
-
1
,
A
2
=
θ
e
-
b
l
T
(
1
-
θ
)
(
1
-
(
1
-
θ
)
e
-
b
T
)
(
1
-
e
-
(
γ
1
+
b
)
T
)
.Proof.
Let(
S
(
t
)
,
V
(
t
)
,
I
(
t
)
) be a solution of (4) satisfied initial condition (5). Since R
1
<
1, one can choose an ɛ
>
0 small enough such that
(19)
β
e
-
b
τ
(
(
Δ
1
)
p
+
δ
(
Δ
2
)
q
)
-
(
γ
+
b
+
α
)
<
0
,
where
(20)
Δ
1
=
γ
e
-
b
ω
b
+
μ
e
b
T
e
b
T
-
1
+
ɛ
,
Δ
2
=
θ
e
-
b
l
T
(
1
-
θ
)
(
1
-
(
1
-
θ
)
e
-
b
T
)
(
1
-
e
-
(
γ
1
+
b
)
T
)
+
ɛ
.
For n
>
n
1, we have
(21)
iiiiiiiiiiiii
d
S
(
t
)
d
t
≤
b
-
b
S
(
t
)
+
γ
e
-
b
ω
,
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
n
∈
N
,
Δ
S
(
t
)
=
-
θ
S
(
t
)
,
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
S
(
t
)
=
μ
,
t
=
n
T
,
n
∈
N
.
By impulsive differential inequality Lemma 1, we have
(22)
S
(
t
)
≤
S
(
n
1
T
+
)
∏
n
1
T
+
<
n
T
<
t
‍
exp
⁡
(
∫
n
1
T
t
‍
(
-
b
)
d
s
)
+
∑
n
1
T
<
n
T
<
t
‍
(
∏
n
T
<
t
j
<
t
‍
exp
⁡
(
∫
n
T
t
‍
(
-
b
)
d
s
)
)
μ
+
∫
n
1
T
t
‍
∏
s
<
n
T
<
t
‍
exp
⁡
(
∫
s
t
‍
(
-
b
)
d
θ
)
γ
e
-
b
ω
d
s
=
S
1
+
S
2
+
S
3
,
where
(23)
S
1
=
S
(
n
1
T
+
)
∏
n
1
T
+
<
n
T
<
t
‍
exp
⁡
(
∫
n
1
t
/
T
‍
(
-
b
)
d
T
ξ
)
=
S
(
n
1
T
+
)
e
-
b
(
t
-
n
1
T
)
,
S
2
=
∑
n
1
T
<
n
T
<
t
‍
(
∏
n
T
<
t
j
<
t
‍
exp
⁡
(
∫
n
T
t
‍
(
-
b
)
d
s
)
)
μ
=
∑
n
1
T
<
n
T
<
t
‍
(
μ
e
-
b
(
t
-
n
T
)
)
=
μ
e
-
b
(
t
-
n
1
)
T
-
e
-
b
(
t
-
(
n
+
n
1
)
)
T
1
-
e
b
T
,
S
3
=
∫
n
1
T
t
‍
∏
s
<
n
T
<
t
‍
exp
⁡
(
∫
s
t
‍
(
-
b
)
d
θ
)
γ
e
-
b
ω
d
s
=
γ
e
-
b
ω
∫
n
1
T
t
‍
∏
s
<
n
T
<
t
‍
exp
⁡
(
∫
s
t
‍
(
-
b
)
d
θ
)
d
s
=
γ
e
-
b
ω
e
-
b
t
b
∫
n
1
T
t
‍
∏
s
<
n
T
<
t
‍
e
b
s
d
(
b
s
)
=
γ
e
-
b
ω
e
-
b
t
b
(
e
b
t
-
e
b
n
1
T
)
.
Thus
(24)
S
(
t
)
≤
S
1
+
S
2
+
S
3
=
S
(
n
1
T
+
)
e
-
b
(
t
-
n
1
T
)
+
μ
e
-
b
(
t
-
n
1
)
T
-
e
-
b
(
t
-
(
n
+
n
1
)
)
T
1
-
e
b
T
+
γ
e
-
b
ω
e
-
b
t
b
(
e
b
t
-
e
b
n
1
T
)
≤
e
-
b
t
S
(
n
1
T
+
)
e
n
1
b
T
+
μ
e
-
b
(
t
-
n
1
T
)
1
-
e
b
T
+
γ
e
-
b
ω
b
e
-
b
(
t
-
n
1
T
)
+
γ
e
-
b
ω
b
+
μ
e
b
T
e
b
T
-
1
,
and then we have
(25)
lim
⁡
t
→
∞
sup
⁡
S
(
t
)
<
γ
e
-
b
ω
b
+
μ
e
b
T
e
b
T
-
1
.
Thus there exists a positive integer n
2
>
n
1 and constant ɛ
>
0 small enough such that, for all t
>
n
2
T,
(26)
S
(
t
)
≤
γ
e
-
b
ω
b
+
μ
e
b
T
e
b
T
-
1
+
ɛ
=
Δ
1
.
Forn
>
n
1, system (4) yields
(27)
d
V
(
t
)
d
t
≤
-
(
γ
1
+
b
)
V
(
t
)
,
t
≠
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
V
(
t
)
=
θ
S
(
t
)
,
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
.
We obtain the following comparison impulsive differential system:
(28)
d
x
(
t
)
d
t
=
-
(
γ
1
+
b
)
x
(
t
)
,
t
≠
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
x
(
t
)
=
θ
S
(
t
)
,
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
.
By Lemma 3, the system has a periodic solution given by
(29)
x
*
(
t
)
=
θ
e
-
b
l
T
e
-
(
γ
1
+
b
)
(
t
-
(
n
+
l
-
1
)
T
)
(
1
-
θ
)
(
1
-
(
1
-
θ
)
e
-
b
T
)
(
1
-
e
-
(
γ
1
+
b
)
T
)
,
iiiiiiiiiiiiiiiiiiii
t
∈
(
(
n
+
l
-
1
)
T
,
(
n
+
l
)
T
]
,
x
(
0
+
)
=
θ
e
-
b
l
T
(
1
-
θ
)
(
1
-
(
1
-
θ
)
e
-
b
T
)
(
1
-
e
-
(
γ
1
+
b
)
T
)
,
which is globally asymptotically stable.
Now, assume thatx
(
t
) is the solution of system (28) with initial value x
(
0
+
)
=
V
0. Then by Lemma 1, we know there exists a positive integer n such that
(30)
V
(
t
)
<
x
(
t
)
<
x
*
(
t
)
+
ɛ
,
t
∈
(
n
T
,
(
n
+
1
)
T
]
.
Hence,
(31)
V
(
t
)
<
x
(
t
)
<
x
*
(
t
)
+
ɛ
<
θ
e
-
b
l
T
(
1
-
θ
)
(
1
-
(
1
-
θ
)
e
-
b
T
)
(
1
-
e
-
(
γ
1
+
b
)
T
)
+
ɛ
0
=
Δ
2
.
From (27), (31), and the third equation in (4), for t
>
n
2
T
+
τ we have
(32)
d
I
(
t
)
d
t
≤
β
e
-
b
τ
(
Δ
1
p
+
δ
Δ
2
q
)
I
(
t
-
τ
)
-
(
γ
+
b
+
α
)
I
(
t
)
.
Consider the comparison equation:
(33)
d
y
(
t
)
d
t
≤
β
e
-
b
τ
(
Δ
1
p
+
δ
Δ
2
q
)
y
(
t
-
τ
)
-
(
γ
+
b
+
α
)
y
(
t
)
.
From (19), we have
(34)
β
e
-
b
τ
(
Δ
1
p
+
δ
Δ
2
q
)
-
(
γ
+
b
+
α
)
<
0
.
According to Lemma 2, we then obtain
(35)
lim
⁡
t
→
∞
y
(
t
)
=
0
.
Notice the fact that I
(
s
)
=
y
(
s
)
=
ϕ
3
(
s
)
>
0 for all s
∈
[
-
τ
,
0
] and I
(
t
)
≥
0, and the comparison theorem implies I
(
t
)
→
0 as t
→
∞. Without loss of generality, we may assume that 0
<
I
(
t
)
<
ɛ
1 for all t
≥
0. By using the first and second equations in (4), we reach
(3.2)
d
S
(
t
)
d
t
≥
-
b
S
(
t
)
-
β
ɛ
1
S
p
(
t
)
,
d
V
(
t
)
d
t
≥
-
δ
β
ɛ
1
V
q
(
t
)
-
γ
1
V
(
t
)
-
b
V
(
t
)
,
iiiiii
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
n
∈
N
,
Δ
S
(
t
)
=
-
θ
S
(
t
)
,
Δ
V
(
t
)
=
θ
S
(
t
)
,
iiiiiiiiiiiiiiiii
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
S
(
t
)
=
μ
,
Δ
V
(
t
)
=
0
,
t
=
n
T
,
n
∈
N
.
For 1
<
q
<
p, we have
(37)
d
S
(
t
)
d
t
≥
-
b
S
(
t
)
-
β
ɛ
1
S
(
t
)
,
d
V
(
t
)
d
t
≥
-
δ
β
ɛ
1
V
(
t
)
-
γ
1
V
(
t
)
-
b
V
(
t
)
,
iiiiii
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
n
∈
N
,
Δ
S
(
t
)
=
-
θ
S
(
t
)
,
Δ
V
(
t
)
=
θ
S
(
t
)
,
iiiiiiiiiiiiiiiii
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
S
(
t
)
=
μ
,
Δ
V
(
t
)
=
0
,
t
=
n
T
,
n
∈
N
considering the following system:
(38)
d
f
(
t
)
(
t
)
d
t
=
-
b
f
(
t
)
-
β
ɛ
1
f
(
t
)
,
d
g
(
t
)
d
t
=
-
δ
β
ɛ
1
g
(
t
)
-
γ
1
g
(
t
)
-
b
g
(
t
)
,
iiiii
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
n
∈
N
,
Δ
f
(
t
)
=
-
θ
f
(
t
)
,
Δ
g
(
t
)
=
θ
f
(
t
)
,
iiiiiiiiiiiiiiiiii
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
f
(
t
)
=
μ
,
Δ
g
(
t
)
=
0
,
t
=
n
T
,
n
∈
N
.
We obtain
(39)
f
~
(
t
)
=
{
(
b
+
β
ɛ
1
)
μ
(
b
+
β
ɛ
1
)
(
1
-
(
1
-
θ
)
e
-
(
b
+
β
ɛ
1
)
T
)
×
e
-
(
b
+
β
ɛ
1
)
(
t
-
(
n
-
1
)
T
)
,
t
∈
(
(
n
-
1
)
T
,
(
n
+
l
-
1
)
T
]
,
(
b
+
β
ɛ
1
)
μ
(
1
-
θ
)
e
-
(
b
+
β
ɛ
1
)
l
T
(
b
+
β
ɛ
1
)
(
1
-
(
1
-
θ
)
e
-
(
b
+
β
ɛ
1
)
T
)
×
e
-
(
b
+
β
ɛ
1
)
(
t
-
(
n
+
l
-
1
)
T
)
,
t
∈
(
(
n
+
l
-
1
)
T
,
n
T
]
,
g
~
(
t
)
=
(
(
θ
(
b
+
β
ɛ
1
)
μ
e
-
(
b
+
β
ɛ
1
)
l
T
-
β
ɛ
1
θ
(
1
-
e
-
(
b
+
β
ɛ
1
)
T
)
)
iii
(
A
)
×
e
-
(
γ
1
+
b
+
δ
β
ɛ
1
)
(
t
-
(
n
+
l
-
1
)
T
)
)
×
(
(
b
+
β
ɛ
1
)
(
1
-
(
1
-
θ
)
e
-
(
b
+
β
ɛ
1
)
T
)
iiiiii
×
(
1
-
e
-
(
γ
1
+
b
)
T
)
)
-
1
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiii
t
∈
(
(
n
+
l
-
1
)
T
,
(
n
+
l
)
T
]
.
Now by using comparison theorem of impulsive equations, for anyɛ
2
>
0 there exists a T
1
>
0 such that
(40)
S
(
t
)
>
f
~
(
t
)
-
ɛ
2
,
V
(
t
)
>
g
~
(
t
)
-
ɛ
2
,
for t
>
T
1. On the other side, from the first and second equations of (4), we have
(41)
iiiiiiiiiiiii
d
S
(
t
)
d
t
≤
-
b
S
(
t
)
+
γ
ɛ
1
e
-
b
ω
,
iiiiiiiiiiiii
d
V
(
t
)
d
t
≤
-
γ
1
V
(
t
)
-
b
V
(
t
)
,
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
n
∈
N
,
Δ
S
(
t
)
=
-
θ
S
(
t
)
,
Δ
V
(
t
)
=
θ
S
(
t
)
,
iiiiiiiiiiiiiiiii
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
S
(
t
)
=
b
,
Δ
V
(
t
)
=
0
,
t
=
n
T
,
n
∈
N
.
Then we have S
(
t
)
≤
h
~
(
t
), V
(
t
)
≤
g
~
(
t
) and h
~
(
t
)
→
S
*
(
t
), g
~
(
t
)
→
V
*
(
t
), as ɛ
1
→
0, where (
h
~
(
t
)
,
g
~
(
t
)
) is a unique positive periodic solution of
(42)
iiiiiiiiiiiii
d
h
(
t
)
d
t
=
-
b
h
(
t
)
+
γ
ɛ
1
e
-
b
ω
,
iiiiiiiiiiiiii
d
g
(
t
)
d
t
=
-
γ
1
g
(
t
)
-
b
g
(
t
)
,
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
n
∈
N
,
Δ
h
(
t
)
=
-
θ
h
(
t
)
,
Δ
g
(
t
)
=
θ
h
(
t
)
,
iiiiiiiiiiiiiii
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
h
(
t
)
=
μ
,
Δ
g
(
t
)
=
0
,
t
=
n
T
,
n
∈
N
,
from which we have that, for n
T
<
t
≤
(
n
+
1
)
T,
(43)
h
~
(
t
)
=
{
γ
ɛ
1
e
-
b
ω
θ
e
-
b
(
1
-
l
)
T
+
b
μ
b
(
1
-
(
1
-
θ
)
e
-
b
T
)
×
e
-
b
(
t
-
(
n
-
1
)
T
)
-
γ
ɛ
1
e
-
b
ω
b
,
(
n
-
1
)
T
<
t
≤
(
n
+
l
-
1
)
T
,
b
μ
(
1
-
θ
)
e
-
b
l
T
+
γ
ɛ
1
e
-
b
ω
θ
b
(
1
-
(
1
-
θ
)
e
-
b
T
)
×
e
-
b
(
t
-
(
n
+
l
-
1
)
T
)
-
γ
ɛ
1
e
-
b
ω
b
,
(
n
+
l
-
1
)
T
<
t
≤
n
T
,
g
~
(
t
)
=
(
θ
b
μ
e
-
b
l
T
-
γ
ɛ
1
e
-
b
ω
θ
(
1
-
e
-
b
T
)
)
e
-
(
γ
1
+
b
)
(
t
-
(
n
+
l
-
1
)
T
)
b
(
1
-
(
1
-
θ
)
e
-
b
T
)
(
1
-
e
-
(
γ
1
+
b
)
T
)
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
(
n
+
l
-
1
)
T
<
t
≤
(
n
+
l
)
T
.
Applying the comparison theorem again, for any ɛ
2
>
0, there exists a T
2
>
0 such that
(44)
S
(
t
)
<
h
~
(
t
)
-
ɛ
2
,
V
(
t
)
<
g
~
(
t
)
-
ɛ
2
,
for t
>
T
2. Let ɛ
1
→
0, and then from (40) and (44) we have
(45)
S
*
(
t
)
-
ɛ
2
<
S
(
t
)
<
S
*
(
t
)
-
ɛ
2
,
V
*
(
t
)
-
ɛ
2
<
V
(
t
)
<
V
*
(
t
)
-
ɛ
2
,
for t large enough, which implies S
(
t
)
→
S
*
(
t
), V
(
t
)
→
V
*
(
t
) as t
→
∞. This completes the proof.Corollary 6.
Ifτ
>
τ
* or μ
<
μ
*, then the infection-free periodic solution (
S
*
(
t
)
,
V
*
(
t
)
,
0
) is globally attractive, where the critical values are given below:
(46)
τ
*
=
1
b
ln
⁡
β
(
A
1
p
+
δ
A
2
q
)
γ
+
b
+
α
,
μ
*
=
(
1
-
e
-
b
T
)
×
(
γ
+
b
+
α
β
e
b
τ
-
δ
A
2
q
p
-
γ
e
-
b
ω
b
-
μ
e
b
T
e
b
T
-
1
)
.


## 4. Permanence

In this section, we discuss the permanence of the infectious population. First, we introduce the following definition.Definition 7.
System (4) is said to be permanent if there exist positive constants m
i, M
i, i
=
1,2
,
3 (independent of initial value), and a finite time T
0, which may depend on the initial condition, such that every positive solution (
S
(
t
)
,
V
(
t
)
,
I
(
t
)
) with initial condition (5) satisfies m
1
≤
S
(
t
)
≤
M
1, m
2
≤
I
(
t
)
≤
M
2, m
3
≤
V
(
t
)
≤
M
3 for all t
>
T
0.Let(47)
S
*
=
γ
+
b
+
α
β
e
-
b
τ
p
,
V
*
=
γ
1
+
b
δ
β
e
-
b
τ
q
,
m
*
=
(
1
/
T
)
ln
⁡
(
(
R
1
(
e
b
T
-
1
+
θ
)
+
1
-
θ
+
e
b
T
)
/
2
)
-
b
β
,
R
1
=
β
e
-
b
τ
γ
+
b
+
α
p
(
μ
(
1
-
θ
)
e
-
b
T
1
-
(
1
-
θ
)
e
-
b
T
)
,
R
2
=
β
e
-
b
τ
γ
1
+
b
q
μ
θ
e
-
(
γ
1
+
b
+
b
l
)
T
(
1
-
e
-
(
γ
1
+
b
)
T
)
(
1
-
(
1
-
θ
)
e
-
b
T
)
,
R
2
=
min
⁡
{
R
1
,
R
2
}
.
Then we have our main result of this section.Theorem 8.
Let1
≤
q
≤
p, if R
2
>
1, and then there exists a positive constant η small enough such that
(48)
I
(
t
)
≥
min
⁡
{
η
m
*
2
,
η
m
*
e
-
(
γ
+
b
+
α
)
ϖ
}
=
m
1
with t large enough.Proof.
As before, we suppose thatX
(
t
)
=
(
S
(
t
)
,
V
(
t
)
,
I
(
t
)
) is a positive solution of system (4) with initial condition (5). Then for t
≥
0, we construct a function as follows:
(49)
U
(
t
)
=
I
(
t
)
+
V
(
t
)
+
β
e
-
b
τ
(
S
*
)
p
×
∫
t
-
τ
t
‍
I
(
ϱ
)
d
ϱ
+
δ
β
e
-
b
τ
(
V
*
)
q
∫
t
-
τ
t
‍
V
(
ϱ
)
d
ϱ
.
And then differentiating U
(
t
) along the trajectory of (4) yields
(50)
U
˙
(
t
)
=
I
˙
(
t
)
+
V
˙
(
t
)
+
β
e
-
b
τ
(
S
*
)
p
I
(
t
)
-
β
e
-
b
τ
(
S
*
)
p
I
(
t
-
τ
)
+
δ
β
e
-
b
τ
(
V
*
)
q
I
(
t
)
-
δ
β
e
-
b
τ
(
V
*
)
q
I
(
t
-
τ
)
=
β
e
-
b
τ
(
S
p
(
t
)
-
(
S
*
)
p
)
I
(
t
-
τ
)
+
β
e
-
b
τ
(
V
q
(
t
)
-
(
V
*
)
q
)
I
(
t
-
τ
)
+
(
β
e
-
b
τ
(
S
*
)
p
-
(
γ
+
b
+
α
)
)
I
(
t
)
+
(
δ
β
e
-
b
τ
(
V
*
)
q
-
(
γ
1
+
b
)
)
I
(
t
)
=
β
e
-
b
τ
(
S
p
(
t
)
-
(
S
*
)
p
)
I
(
t
-
τ
)
+
β
e
-
b
τ
(
V
q
(
t
)
-
(
V
*
)
q
)
I
(
t
-
τ
)
=
β
e
-
b
τ
(
S
p
-
1
(
t
)
+
S
p
-
2
(
t
)
S
*
+
⋯
iiiiiiiiii
+
S
(
t
)
(
S
*
)
p
-
2
+
(
S
*
)
p
-
1
)
×
(
S
(
t
)
-
S
*
)
I
(
t
-
τ
)
+
δ
β
e
-
b
τ
(
V
q
-
1
(
t
)
+
V
q
-
2
(
t
)
V
*
+
⋯
iiiiiiiiiiiiiii
+
V
(
t
)
(
V
*
)
q
-
2
+
(
V
*
)
q
-
1
)
×
(
V
(
t
)
-
V
*
)
I
(
t
-
τ
)
for t
≥
0. Let
(51)
m
*
=
(
1
/
T
)
ln
⁡
(
(
R
1
(
e
b
T
-
1
+
θ
)
+
1
-
θ
+
e
b
T
)
/
2
)
-
b
β
,
S
*
=
γ
+
b
+
α
β
e
-
b
τ
p
.
Since R
2
>
1, we get R
1
>
1, R
2
>
1. Then we have m
*
>
0. And from R
1
>
1, we can get
(52)
β
e
-
b
τ
γ
+
b
+
α
p
(
μ
(
1
-
θ
)
e
-
b
T
1
-
(
1
-
θ
)
e
-
b
T
)
>
1
.
Thus, we have
(53)
μ
(
1
-
θ
)
e
-
b
T
1
-
(
1
-
θ
)
e
-
b
T
>
γ
+
b
+
α
β
e
-
b
τ
p
=
S
*
.
Form R
2
>
1, we have
(54)
β
e
-
b
τ
γ
1
+
b
q
μ
θ
e
-
(
γ
1
+
b
+
b
l
)
T
(
1
-
e
-
(
γ
1
+
b
)
T
)
(
1
-
(
1
-
θ
)
e
-
b
T
)
>
1
;
that is,
(55)
μ
θ
e
-
(
γ
1
+
b
+
b
l
)
T
(
1
-
e
-
(
γ
1
+
b
)
T
)
(
1
-
(
1
-
θ
)
e
-
b
T
)
>
γ
1
+
b
β
e
-
b
τ
q
=
V
*
.
We can take η small enough such that
(56)
μ
(
1
-
θ
)
e
-
(
β
η
m
*
+
b
)
T
1
-
(
1
-
θ
)
e
-
(
β
η
m
*
+
b
)
T
>
S
*
,
μ
θ
e
-
(
β
(
δ
+
l
)
η
m
*
+
γ
1
+
b
+
b
l
)
T
(
1
-
e
-
(
δ
β
η
m
*
+
γ
1
+
b
)
T
)
(
1
-
(
1
-
θ
)
e
-
(
β
η
m
*
+
b
)
T
)
>
V
*
.
Thus we can choose ɛ
1
,
ɛ
2
>
0 to be small enough such that
(57)
S
*
<
μ
(
1
-
θ
)
e
-
(
β
η
m
*
+
b
)
T
1
-
(
1
-
θ
)
e
-
(
β
η
m
*
+
b
)
T
-
ɛ
1
≡
S
Δ
,
V
*
<
μ
θ
e
-
(
β
(
δ
+
l
)
η
m
*
+
γ
1
+
b
+
b
l
)
T
(
1
-
e
-
(
δ
β
η
m
*
+
γ
1
+
b
)
T
)
(
1
-
(
1
-
θ
)
e
-
(
β
η
m
*
+
b
)
T
)
-
ɛ
2
≡
V
Δ
.
Then we claim that there exists anm
2
>
0 such that I
(
t
)
>
m
2 for t is large enough. We next prove this claim in two steps.
Step I. For any positive constant t
0, that I
(
t
)
≤
η
m
* for all t
≥
t
0 is not true.
Otherwise, there is a positive constantt
0, such that I
(
t
)
≤
η
m
* for all t
≥
t
0. First, if I
(
t
)
<
η
m
* for all t
≥
t
0, it follows from the first, fourth, and fifth equations of (4) that, for t
≥
t
0,
(58)
d
S
(
t
)
d
t
≥
-
(
β
η
m
*
+
b
)
S
(
t
)
,
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
Δ
S
(
t
)
=
-
θ
S
(
t
)
,
t
=
(
n
+
l
-
1
)
T
,
Δ
S
(
t
)
=
μ
,
t
=
n
T
.
By Lemma 1, there exists T
1
>
t
0
+
τ so that for t
>
T
1
(59)
S
(
t
)
>
μ
(
1
-
θ
)
e
-
(
β
η
m
*
+
b
)
T
1
-
(
1
-
θ
)
e
-
(
β
η
m
*
+
b
)
T
-
ɛ
≡
S
Δ
.
Similarly, from the second and the fourth equations of (4), we have
(60)
d
V
(
t
)
d
t
≥
-
(
δ
β
η
m
*
+
γ
1
+
b
)
V
(
t
)
,
t
≠
(
n
+
l
-
1
)
T
,
Δ
V
(
t
)
=
θ
S
(
t
)
,
t
=
(
n
+
l
-
1
)
T
,
and for t
>
T
1,
(61)
V
(
t
)
≥
μ
θ
e
-
(
β
(
δ
+
l
)
η
m
*
+
γ
1
+
b
+
b
l
)
T
(
1
-
e
-
(
δ
β
η
m
*
+
γ
1
+
b
)
T
)
(
1
-
(
1
-
θ
)
e
-
(
β
η
m
*
+
b
)
T
)
-
ɛ
≡
V
Δ
.
Then, by (50), for t
≥
T
1,
(62)
U
˙
(
t
)
=
β
e
-
b
τ
(
S
p
-
1
(
t
)
+
S
p
-
2
(
t
)
S
*
+
⋯
iiiiiiiiii
+
S
(
t
)
(
S
*
)
p
-
2
+
(
S
*
)
p
-
1
)
×
(
S
(
t
)
-
S
*
)
I
(
t
-
τ
)
+
δ
β
e
-
b
τ
(
V
q
-
1
(
t
)
+
V
q
-
2
(
t
)
V
*
+
⋯
iiiiiiiiiiiiiii
+
V
(
t
)
(
V
*
)
q
-
2
+
(
V
*
)
q
-
1
)
×
(
V
(
t
)
-
V
*
)
I
(
t
-
τ
)
>
p
β
e
-
b
τ
(
S
*
)
p
-
1
(
S
Δ
-
S
*
)
I
(
t
-
τ
)
+
q
δ
β
e
-
b
τ
(
V
*
)
q
-
1
(
V
Δ
-
V
*
)
I
(
t
-
τ
)
.
Let
(63)
I
L
=
min
⁡
t
∈
[
T
1
,
T
1
+
τ
]
I
(
t
)
.
We can prove that I
(
t
)
≥
I
L for all t
≥
T
1. Otherwise, there exists a nonnegative constant T
2 such that I
(
t
)
≥
I
L for t
∈
[
T
1
,
T
1
+
τ
+
T
2
]
,
I
(
T
1
+
τ
+
T
2
)
=
I
L, and I
˙
(
T
1
+
τ
+
T
2
)
≤
0. Then from the second equation of (4) and (37), we easily see that
(64)
I
˙
(
T
1
+
τ
+
T
2
)
≥
(
β
e
-
b
τ
S
p
(
t
)
+
δ
β
e
-
b
τ
V
q
(
t
)
-
(
γ
+
b
+
α
)
)
I
L
=
(
γ
+
b
+
α
)
(
β
e
-
b
τ
S
p
(
t
)
γ
+
b
+
α
+
δ
β
e
-
b
τ
V
q
(
t
)
γ
+
b
+
α
-
1
)
I
L
>
(
γ
+
b
+
α
)
(
(
S
Δ
S
*
)
p
+
γ
1
+
b
γ
+
b
+
α
(
V
Δ
V
*
)
q
-
1
)
I
L
>
0
,
which is a contradiction. Hence I
(
t
)
≥
I
L
>
0 for all t
≥
T
1. Equation (62) implies
(65)
d
U
(
t
)
d
t
>
p
β
e
-
b
τ
(
S
*
)
p
-
1
(
S
Δ
-
S
*
)
I
(
t
-
τ
)
+
q
δ
β
e
-
b
τ
(
V
*
)
q
-
1
(
V
Δ
-
V
*
)
I
(
t
-
τ
)
>
0
.
It then follows that U
(
t
)
→
+
∞ as t
→
+
∞. This is a contradiction to U
(
t
)
≤
(
α
+
γ
+
γ
1
+
2
b
)
τ
+
2. Therefore, for any positive constant t
0, the inequality I
(
t
)
<
η
m
* cannot hold for all t
≥
t
0.
Step II. From Step I, we only need to consider the followng: (i) I
(
t
)
>
η
m
* for all t large enough and (ii) I
(
t
) oscillates about η
m
* for all large t. However, Case (i) is obvious in the result of this theorem, so we only need to consider Case (ii), in which we will show that I
(
t
)
≥
m
1 for all large t where
(66)
m
1
=
min
⁡
{
η
m
*
2
,
η
m
*
e
-
(
γ
+
b
+
α
)
ϖ
}
.
First, we notice there exist two positive constants t
¯, φ such that
(67)
I
(
t
¯
)
=
I
(
t
¯
+
φ
)
=
I
*
,
I
(
t
)
<
η
m
*
,
for
t
¯
<
t
<
t
¯
+
φ
.
Second, because I
(
t
) is bounded continuous function and I
(
t
) has no pulse, we can get that I
(
t
) is uniformly continuous. Therefore there exists a constant T
3 (with 0
<
T
3
<
ϖ and T
3 is independent of the choice of t
¯) such that I
(
t
)
>
η
m
*
/
2 for all t
¯
≤
t
≤
t
¯
+
T
3.
Ifφ
≤
T
3, our aim is obtained.
IfT
3
<
φ
≤
ϖ, from the second equation of (4) we have that I
˙
(
t
)
≥
-
(
γ
+
b
+
α
)
I
(
t
) for t
¯
<
t
≤
t
¯
+
φ. Then we have I
(
t
)
≥
η
m
*
e
-
(
γ
+
b
+
α
)
ϖ for t
¯
<
t
≤
t
¯
+
φ
≤
t
¯
+
ϖ since I
(
t
¯
)
=
η
m
*. It is clear that I
(
t
)
≥
m
1 for t
¯
<
t
≤
t
¯
+
φ.
Ifφ
≥
ϖ, then we have I
(
t
)
≥
m
2 for t
¯
<
t
≤
t
¯
+
ϖ. We then can easily prove I
(
t
)
≥
m
1 for t
¯
+
ϖ
≤
t
≤
t
¯
+
φ. Since the interval [
t
¯
,
t
¯
+
φ
] is arbitrarily chosen, we know that I
(
t
)
≥
m
1 holds for t large enough. Finally, noticing the choice of m
1 is independent of the positive solution of (4), we completed our proof.Theorem 9.
Let1
≤
q
≤
p, if R
2
>
1, and then system (4) is permanent.Proof.
Suppose thatX
(
t
)
=
(
S
(
t
)
,
V
(
t
)
,
I
(
t
)
) is a positive solution of system (4) with initial conditions (5). Then from system (4), we have
(68)
iiiiiiiiii
d
S
(
t
)
d
t
≥
-
(
b
+
β
)
S
(
t
)
,
iiiiiiiiii
d
V
(
t
)
d
t
≥
-
(
δ
β
+
γ
1
+
b
)
V
(
t
)
,
t
≠
(
n
+
l
-
1
)
T
,
t
≠
n
T
,
n
∈
N
,
Δ
S
(
t
)
=
-
θ
S
(
t
)
,
Δ
V
(
t
)
=
θ
S
(
t
)
,
iiiiiiiiiiiiiiii
t
=
(
n
+
l
-
1
)
T
,
n
∈
N
,
Δ
S
(
t
)
=
μ
,
Δ
V
(
t
)
=
0
,
t
=
n
T
,
n
∈
N
.
As what we did in the proof of Theorem 5, we can prove that there exist t large enough and ɛ
>
0 small enough such that
(69)
S
(
t
)
≥
μ
(
1
-
θ
)
e
-
(
b
+
β
)
T
1
-
(
1
-
θ
)
e
-
(
b
+
β
)
T
-
ɛ
=
m
3
,
V
(
t
)
≥
μ
θ
e
-
(
b
+
γ
1
+
δ
β
+
(
b
+
β
)
l
)
T
(
1
-
(
1
-
θ
)
e
-
(
b
+
β
)
T
)
(
1
-
e
-
(
b
+
γ
1
+
δ
β
)
T
)
-
ɛ
2
=
m
4
.
Then for D
=
{
(
S
,
V
,
I
)
∈
R
+
3
∣
S
(
t
)
+
V
(
t
)
+
I
(
t
)
≤
1
}, by Theorem 8, we have
(70)
m
3
≤
S
(
t
)
≤
1
,
m
1
≤
I
(
t
)
≤
1
,
m
4
≤
V
(
t
)
≤
1
for t large enough. Thus the system (4) is uniformly permanent.


## 5. Numerical Simulations and Discussions

Next, we carry out numerical simulations to illustrate the theoretical results obtained in the previous sections. We first set the parameters as follows:b
=
0.2, β
=
0.5, α
=
0.05, γ
=
0.04, δ
=
0.02, γ
1
=
0.06, p
=
1.5, q
=
1.25, T
=
1.5, ω
=
1, τ
=
1, θ
=
0.4, and μ
=
1.4. Straightforward calculation shows R
2
=
2.6155
>
1. Then by Theorem 8, the disease will be permanent (please see Figures 1(a), 1(b), 1(c), and 1(d)). In order to show the effect of τ, we decrease τ to 4, and other parameters are the same with those in Figure 1, and the infection-free periodic solution of system (4) is globally attractive. This phenomenon is also seen from our theoretical analysis as in this case R
1
=
0.0339
<
1 and then according to Theorem 5, the disease will be eradicated; please see Figure 2(a).The results of numerical simulation on the threshold valuesR
2
=
2.6155
>
1, where p
=
1.5, q
=
1.25.
(a)
Time series ofS, V, and I. The disease will be permanent
(b)
Phase portrait ofS
(
t
), V
(
t
) of the system  (4)
(c)
Phase portrait ofS
(
t
), I
(
t
) of the system  (4)
(d)
Phase portrait ofS
(
t
), V
(
t
), and I
(
t
) of the system  (4)The results of numerical simulation on the threshold valuesR
1
=
0.0339
<
1, where p
=
1.5, q
=
1.25.
(a)
Time series ofS, V, and I. The disease dies out
(b)
Phase portrait ofS
(
t
), V
(
t
) of the system  (4)
(c)
Phase portrait ofS
(
t
), I
(
t
) of the system  (4)
(d)
Phase portrait ofS
(
t
), V
(
t
), and I
(
t
) of the system  (4)If we keepτ
=
ω
=
1 and μ
=
1, as the same with those in Figure 1, but increase vaccination proportion of susceptible persons θ to 0.9, then the disease will be eradicated; see Figure 3(a). If we keep τ
=
ω
=
1 and θ
=
0.4 and decrease μ to 0.2, then the disease also will be eradicated; see Figure 4(a).The results of numerical simulation on the threshold valuesR
1
=
0.0449
<
1, where p
=
1.5, q
=
1.25.
(a)
Time series ofS, V, and I. The disease dies out
(b)
Phase portrait ofS
(
t
), V
(
t
) of the system  (4)
(c)
Phase portrait ofS
(
t
), I
(
t
) of the system  (4)
(d)
Phase portrait ofS
(
t
), V
(
t
), and I
(
t
) of the system  (4)The results of numerical simulation on the threshold valuesR
1
=
0.0352
<
1, where p
=
1.5, q
=
1.25.
(a)
Time series ofS, V, and I. The disease dies out
(b)
Phase portrait ofS
(
t
), V
(
t
) of the system  (4)
(c)
Phase portrait ofS
(
t
), I
(
t
) of the system  (4)
(d)
Phase portrait ofS
(
t
), V
(
t
), and I
(
t
) of the system  (4)And if we keepτ
=
ω
=
4, μ
=
1 but decrease θ to 0.1, then the disease will be permanent; see Figure 5. If we keep τ
=
ω
=
4 and θ
=
0.4 and increase μ to 2, then the disease also will be permanent; see Figure 6. For details please see Table 1.Table 1
ω
τ
θ
μ
R
i
Status of the disease
1
1
0.4
1.4
R
2
=
2.6155
>
1
Permanence
4
4
0.4
1.4
R
1
=
0.0339
<
1
Eradication
1
1
0.9
1.4
R
1
=
0.0449
<
1
Eradication
1
1
0.4
0.2
R
1
=
0.0352
<
1
Eradication
4
4
0.1
1.4
R
2
=
1.1036
>
1
Permanence
4
4
0.4
2
R
2
=
2.3121
>
1
PermanenceThe results of numerical simulation on the threshold valuesR
2
=
1.1036
>
1, where p
=
1.5, q
=
1.25.
(a)
Time series ofS, V, and I. The disease will be permanent
(b)
Phase portrait ofS
(
t
), V
(
t
) of the system  (4)
(c)
Phase portrait ofS
(
t
), I
(
t
) of the system  (4)
(d)
Phase portrait ofS
(
t
), V
(
t
), and I
(
t
) of the system  (4)The results of numerical simulation on the threshold valuesR
2
=
2.3121
>
1, where p
=
1.5, q
=
1.25.
(a)
Time series ofS, V, and I. The disease will be permanent
(b)
Phase portrait ofS
(
t
), V
(
t
) of the system  (4)
(c)
Phase portrait ofS
(
t
), I
(
t
) of the system  (4)
(d)
Phase portrait ofS
(
t
), V
(
t
), and I
(
t
) of the system  (4)Lastly, we conclude our paper as follows. In this paper, we proposed an SVEIRS model, which is a new epidemic model with periodic pulse vaccination and pulse population input at two different fixed moments. Our primary result is to investigate the effect of impulsive vaccination, pulse population input, and time delays to the dynamics of population model. With the help of comparison theorems, we proved the existence of the “infection-free” periodic solution and obtained the conditions for global attractivity of the “infection-free” periodic solution and the conditions for the permanence of the system. All the theoretical results show that we believe it might be helpful in disease control: people can select appropriate vaccination rate and population input rate according to our theoretical results to control diseases.



---

*Source: 542154-2014-05-26.xml*