# Stability Results for an Age-Structured SIS Epidemic Model with Vector Population

Authors: He-Long Liu; Jing-Yuan Yu; Guang-Tian Zhu
Journal: Journal of Applied Mathematics
 (2015)

---

## Abstract

We formulate an age-structured SIS epidemic model with periodic parameters, which includes host population and vector population. The host population is described by two partial differential equations, and the vector population is described by a single ordinary differential equation. The existence problem for endemic periodic solutions is reduced to a fixed point problem of a nonlinear integral operator acting on locally integrable periodic functions. We obtain that if the spectral radius of the Fréchet derivative of the fixed point operator at zero is greater than one, there exists a unique endemic periodic solution, and we investigate the global attractiveness of disease-free steady state of the normalized system.


---

## Body

## 1. Introduction

Some infectious diseases such as malaria, Chagas disease, and dengue fever are transmitted via vector. The seasonality of vector population decides that the transmission of infectious diseases is periodic. So, in order to model the seasonal spread of infectious diseases many authors have studied differential equations systems with periodic parameters [1–6].In order to reflect the effect of demographic behavior of individuals, scholars have recognized that age-structured epidemic models are more realistic, since any disease prevention policy depends on the age structure of host population, and instantaneous death and infection rates depend on the age. Since the pioneer work of McKendrick [3], many authors have studied various age-structured epidemic models [2, 7–12]. However, due to their relatively complex form, the analyses of mathematical properties are difficult especially in the local and global stability of steady states.The model we study in this paper is an age-structured SIS epidemic model with vector population, host population is divided into two compartments, susceptibles and infectives, and we assume that recovered individuals cannot obtain immunity and directly go back to the susceptibles. The vector population is divided into two groups, susceptibles and infectives. Generally speaking, the number of susceptible population is much greater than the number of infective population in a vector population. So, we may assume that the change of infectives’ size does not affect the total number of susceptible vector population.The age-structured SIS epidemic models have been studied in [8, 9, 12–14]. The general form of these models to a periodic system was given in [6]; the paper analyzed age structured SIS models with seasonal periodicities and vertical transmission and studied the global stability of a nontrivial endemic periodic solution.The age-structured epidemic models with vector population have been studied in [15–18]. In [15], authors proved that the population dynamics of malaria, formulated vector-host model for malaria, and used the system of ordinary differential equations to describe the model. Paper [16] discussed a vector-host model for the spread of Chagas disease with infection-age. In [17], a deterministic model showed that the age-structured model underwent the phenomenon of backward bifurcation at R
0
=
1 under certain conditions and that the backward bifurcation feature was caused by malaria-induced mortality in humans. In [18], authors studied the existence and uniqueness of endemic periodic solution of an age-structured SIS epidemic model with periodic parameters.In our model, host population is described by two partial differential equations, and infective vector population is described by a single ordinary differential equation. Integrating the partial differential equations along the characteristic lines, we can normalize them to a partial differential equation of fraction of infected population and get an expression of the infective host population. Integrating ordinary differential equation, we obtain an expression of the infected vector population. Using these expressions, we can obtain an integral equation, which is a fixed point equation in locally integrable time-periodic functions. From the fixed point theory, we obtain that there exists a unique endemic periodic solution under certain conditions and investigate the global attractiveness of disease-free steady state of the normalized system.This paper is organized as follows: Section2 introduces an age-structured SIS epidemic model with vector population. In Section 3, we show the well-posedness of the time evolution problem. In Section 4, we prove the existence of endemic periodic solution of the system in case that threshold value is greater than one. In Section 5, we get that the nontrivial solution is unique if threshold value is greater than one. In Section 6, we study the global attractiveness of infection-free state E
0. Section 7 contains some discussions of the results.


## 2. The Model

In this section, we formulate an age-structured SIS model, which includes host population and vector population. The host population is divided into two classes: susceptible and infective. LetS
(
a
,
t
) and I
(
a
,
t
) be the age-densities of, respectively, the susceptible and infective host population at time t. In the vector population, we assume the number of susceptible individuals is far greater than the number of infective individuals, and let M
(
t
) (known T-periodic function) be the number of susceptible population and V
(
t
) the number of infective population. Let N
(
a
)
=
S
(
a
,
t
)
+
I
(
a
,
t
) be the density with respect to age of the total number of the host population and satisfy
(1)
N
a
=
μ
*
N
e
-
∫
0
a
μ
1
σ
d
σ
,
where the constant N is the total size of the host population, μ
* is the crude death rate of the host population, and μ
1
(
a
) denotes the instantaneous death rate at age a of the host population, is nonnegative, locally integrable on [
0
,
r
m
) (r
m
<
+
∞ denotes the maximum attainable age), and satisfies
(2)
∫
0
r
m
μ
1
σ
d
σ
=
+
∞
.
The crude death rate of the host population is determined such that
(3)
μ
*
∫
0
r
m
f
a
d
a
=
1
,
where f
(
a
)
=
e
-
∫
0
a
μ
1
(
σ
)
d
σ is the survival function. We have the relation
(4)
N
a
=
μ
*
N
f
a
.
Let δ
1 be the number of bites per vector per unit time and c be the proportion of infected bites that give rise to infection. Then the force of infection for the host population is defined by
(5)
δ
1
c
V
t
=
η
V
t
,
η
=
δ
1
c
.
Let δ
2 be the proportion of bites to infected hosts that give rise to infection in vector. Then the number of new infection of vectors per unit time from infected hosts is given by
(6)
δ
1
δ
2
∫
0
r
m
I
a
,
t
d
a
N
M
t
=
δ
∫
0
r
m
I
a
,
t
d
a
N
M
t
,
δ
=
δ
1
δ
2
.
Let γ
(
a
,
t
) be the age-specific recovery rate in the host population and μ
2
(
t
) the per capita death rate of vectors. Biologically speaking, γ
(
a
,
t
) and μ
2
(
t
) are periodic in time t, maybe different in period. But for theoretical analysis, we assume that their periods are the same. Moreover we assume that the death rate of the host population is not affected by the presence of the disease.With these assumptions, we obtain the following system of equations which describes the dynamics of the vector-host model:(7)
∂
S
a
,
t
∂
t
+
∂
S
a
,
t
∂
a
=
-
μ
1
a
+
η
V
t
S
a
,
t
+
γ
a
,
t
I
a
,
t
,
∂
I
a
,
t
∂
t
+
∂
I
a
,
t
∂
a
=
η
V
t
S
a
,
t
-
μ
1
a
+
γ
a
,
t
I
a
,
t
,
d
V
t
d
t
=
δ
∫
0
r
m
I
a
,
t
d
a
N
M
t
-
μ
2
t
V
t
,
with boundary and initial conditions:
(8)
S
0
,
t
=
μ
*
N
,
I
(
0
,
t
)
=
0
,
S
a
,
0
=
S
0
a
,
I
a
,
0
=
I
0
a
,
V
0
=
V
0
,
where
(9)
S
0
a
≥
0
,
I
0
(
a
)
≥
0
,
V
0
≥
0
,
S
0
a
+
I
0
a
=
N
a
.
To simplify the model, let
(10)
S
a
,
t
=
f
a
μ
*
N
s
a
,
t
,
I
a
,
t
=
f
a
μ
*
N
i
a
,
t
.
The system (7) and (8) can be rewritten as
(11)
∂
s
a
,
t
∂
t
+
∂
s
a
,
t
∂
a
=
-
η
V
t
s
a
,
t
+
γ
a
,
t
i
a
,
t
,
∂
i
a
,
t
∂
t
+
∂
i
a
,
t
∂
a
=
η
V
t
s
a
,
t
-
γ
a
,
t
i
a
,
t
,
d
V
t
d
t
=
δ
μ
*
∫
0
r
m
f
a
i
a
,
t
d
a
M
t
-
μ
2
(
t
)
V
(
t
)
,
with boundary and initial conditions:
(12)
s
0
,
t
=
1
,
i
0
,
t
=
0
,
s
a
,
0
=
s
0
a
,
i
a
,
0
=
i
0
a
,
V
0
=
V
0
.
The following results will be used in Section 3.From the third equation of the system (11), we have
(13)
d
V
t
d
t
≤
δ
M
t
-
μ
2
t
V
t
.
Integrating differential inequality (13), we have
(14)
V
t
≤
V
0
e
-
∫
0
t
μ
2
σ
d
σ
+
δ
e
-
∫
0
t
μ
2
σ
d
σ
∫
0
t
e
∫
0
τ
μ
2
σ
d
σ
M
τ
d
τ
.
Let t
=
m
T
+
t
′, 0
≤
t
′
<
T, m
∈
N
+. Since M
(
t
)
,
μ
2
(
t
) are T-periodic functions, we obtain
(15)
V
t
≤
V
0
e
-
∫
0
m
T
μ
2
(
σ
)
d
σ
+
δ
e
-
∫
0
m
T
μ
2
(
σ
)
d
σ
∫
0
(
m
+
1
)
T
e
∫
0
τ
μ
2
σ
d
σ
M
τ
d
τ
≤
V
0
e
-
m
∫
0
T
μ
2
(
σ
)
d
σ
+
δ
∫
0
T
e
∫
0
τ
μ
2
σ
d
σ
M
τ
d
τ
1
-
e
-
m
∫
0
T
μ
2
(
σ
)
d
σ
e
∫
0
T
μ
2
(
σ
)
d
σ
-
1
≤
V
0
+
∫
0
T
e
∫
0
τ
μ
2
σ
d
σ
M
τ
d
τ
δ
e
∫
0
T
μ
2
σ
d
σ
-
1
<
+
∞
.


## 3. Existence and Uniqueness of Solution

From the first equation and the second equation of the system (11) and (12), we obtain that s
(
a
,
t
)
+
i
(
a
,
t
)
=
1. So, the system (11) and (12) can be reduced to two equations for i
(
a
,
t
) and V
(
t
) as
(16)
∂
i
a
,
t
∂
t
+
∂
i
a
,
t
∂
a
=
η
V
t
1
-
i
a
,
t
-
γ
a
,
t
i
a
,
t
,
d
V
t
d
t
=
δ
μ
*
∫
0
r
m
f
a
i
a
,
t
d
a
M
t
-
μ
2
t
V
t
,
with boundary and initial conditions:
(17)
i
0
,
t
=
0
,
i
a
,
0
=
i
0
a
,
V
0
=
V
0
.We consider the initial-boundary value problem of the system (16) and (17) as an abstract Cauchy problem:
(18)
x
′
t
=
A
x
+
F
t
,
x
,
x
0
=
i
0
a
,
V
0
T
,
where
(19)
x
=
x
1
,
x
2
T
=
i
a
,
t
,
V
t
T
,
A
:
D
A
⟶
X
,
A
x
=
-
d
x
1
d
a
-
x
2
,
D
(
A
)
=
x
,
x
1
∈
W
1
0
,
+
∞
,
x
2
∈
R
,
x
1
0
=
0
,
X
=
L
1
0
,
r
m
×
R
,
endowed with the norm
(20)
x
=
x
1
+
x
2
,
x
1
=
∫
0
r
m
x
1
a
d
a
,
x
=
x
1
,
x
2
T
∈
X
.
Suppose γ
1
(
a
,
t
)
,
μ
1
(
a
)
∈
L
∞
[
0
,
+
∞
); we define
(21)
F
:
X
⟶
X
,
F
t
,
x
=
η
x
2
1
-
x
1
-
γ
a
,
t
x
1
δ
μ
*
∫
0
r
m
f
a
x
1
a
d
a
M
t
-
μ
2
t
-
1
x
2
.
The system (18) is a semilinear nonautonomous Cauchy problem, we easily obtain that the operator A is the infinitesimal generator of C
0-semigroup T
(
t
), t
≥
0, and F is continuously Fréchet differentiable on X. Then for each x
0
∈
D
(
A
), there exists a maximal interval of existence [
0
,
m
) and a unique mild solution t
→
x
(
t
,
x
0
) (see [8, 14, 19]), which satisfies (18), where either m
=
+
∞ or m
<
+
∞, in the case lim
t
→
m
⁡
x
(
t
,
x
0
)
=
lim
t
→
m
⁡
(
i
a
,
t
+
V
t
)
=
+
∞. Since i
(
a
,
t
)
≤
1, and from (15)  V
(
t
)
<
+
∞, we easily obtain m
=
+
∞. So, we have the following theorem.Theorem 1.
The initial-boundary value problem (18), that is, the system (16) and (17), has a unique nonnegative mild solution x
=
(
i
(
a
,
t
)
,
V
(
t
)
)
T.


## 4. Existence of Endemic Periodic Solution

In this section, we discuss existence of endemic periodic solution of the system (16) and (17). Let X
T
=
X
1
,
T
×
X
2
,
T, where X
1
,
T is the set of locally integrable T-periodic L
1
[
0
,
r
m
)-valued functions with norm
(22)
x
1
X
1
,
T
=
∫
0
T
x
1
(
t
)
d
t
=
∫
0
T
∫
0
r
m
x
1
a
,
t
d
a
d
t
,
x
1
∈
X
1
,
T
,
X
2
,
T is the set of locally integrable T-periodic functions with norm
(23)
x
2
X
2
,
T
=
∫
0
T
x
2
t
d
t
,
x
2
∈
X
2
,
T
,
and X
T
+
=
X
1
,
T
+
×
X
2
,
T
+; X
T
+, X
1
,
T
+, X
2
,
T
+ are their positive cone, respectively. We give the state space of the system (16) and (17) as follows:
(24)
X
T
*
=
x
∈
X
T
+
:
0
≤
x
1
a
,
t
≤
1
,
x
2
t
≥
0
.
If (
i
*
(
a
,
t
)
,
V
*
(
t
)
)
T
∈
X
T
*
-
{
0
} is an endemic T-periodic solution of the system (16) and (17), then it satisfies
(25)
∂
i
*
a
,
t
∂
t
+
∂
i
*
a
,
t
∂
a
=
η
V
*
t
1
-
i
*
a
,
t
-
γ
a
,
t
i
*
a
,
t
,
d
V
*
t
d
t
=
δ
μ
*
∫
0
r
m
f
a
i
*
a
,
t
d
a
M
t
-
μ
2
t
V
*
t
,
i
i
*
0
,
t
=
0
.
Integrating the first equation of (25) along the characteristic lines, we have
(26)
i
*
a
,
t
=
η
∫
0
a
V
*
t
-
a
+
σ
e
-
∫
σ
a
(
η
V
*
(
t
-
a
+
ρ
)
+
γ
(
t
-
a
+
ρ
,
ρ
)
)
d
ρ
d
σ
,
and 0
≤
i
*
(
a
,
t
)
≤
1.From the second equation of (25), we obtain(27)
V
*
t
=
μ
*
δ
e
-
∫
0
t
μ
2
τ
d
τ
·
∫
0
t
∫
0
r
m
f
a
i
*
a
,
τ
M
τ
e
∫
0
τ
‍
μ
2
ρ
d
ρ
d
a
d
τ
∫
0
r
m
+
V
0
*
,
where V
0
*
=
∫
0
T
∫
0
r
m
f
a
i
*
a
,
τ
M
τ
e
∫
0
τ
μ
2
ρ
d
ρ
d
a
d
τ
/
e
∫
0
T
μ
2
(
ρ
)
d
ρ
-
1.Substituting (26) into (27), we have
(28)
V
*
t
=
μ
*
δ
η
e
-
∫
0
t
μ
2
(
τ
)
d
τ
·
∫
0
t
M
(
τ
)
e
∫
0
τ
μ
2
(
ρ
)
d
ρ
∫
0
r
m
f
(
a
)
·
∫
0
a
V
*
τ
-
a
+
σ
·
e
-
∫
σ
a
η
V
*
τ
-
a
+
ρ
+
γ
τ
-
a
+
ρ
,
ρ
d
ρ
d
σ
d
a
d
τ
+
1
e
∫
0
T
μ
2
ρ
d
ρ
-
1
·
∫
0
T
M
(
τ
)
e
∫
0
τ
μ
2
(
ρ
)
d
ρ
∫
0
r
m
f
(
a
)
·
∫
0
a
V
*
τ
-
a
+
σ
·
e
-
∫
σ
a
η
V
*
τ
-
a
+
ρ
+
γ
τ
-
a
+
ρ
,
ρ
d
ρ
d
σ
d
a
d
τ
∫
0
a
=
μ
*
δ
η
e
-
∫
0
t
μ
2
(
τ
)
d
τ
·
∫
0
t
M
(
τ
)
e
∫
0
τ
μ
2
(
ρ
)
d
ρ
∫
0
r
m
f
(
a
)
·
∫
0
a
V
*
τ
-
σ
·
e
-
∫
a
-
σ
a
η
V
*
τ
-
a
+
ρ
+
γ
τ
-
a
+
ρ
,
ρ
d
ρ
d
σ
d
a
d
τ
+
1
e
∫
0
T
μ
2
ρ
d
ρ
-
1
·
∫
0
T
M
τ
e
∫
0
τ
μ
2
ρ
d
ρ
∫
0
r
m
f
a
·
∫
0
a
‍
V
*
(
τ
-
σ
)
·
e
-
∫
a
-
σ
a
‍
η
V
*
τ
-
a
+
ρ
+
γ
τ
-
a
+
ρ
,
ρ
d
ρ
d
σ
d
a
d
τ
∫
0
T
=
μ
*
δ
η
e
-
∫
0
t
μ
2
(
τ
)
d
τ
·
η
∫
0
t
M
(
τ
)
e
∫
0
τ
μ
2
(
ρ
)
d
ρ
∫
0
r
m
f
(
a
)
·
∫
0
a
V
*
τ
-
σ
·
e
-
∫
0
σ
η
V
*
τ
-
ρ
+
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
σ
d
a
d
τ
+
1
e
∫
0
T
μ
2
ρ
d
ρ
-
1
·
∫
0
T
M
(
τ
)
e
∫
0
τ
μ
2
(
ρ
)
d
ρ
∫
0
r
m
f
(
a
)
·
∫
0
a
V
*
(
τ
-
σ
)
·
e
-
∫
0
σ
η
V
*
τ
-
ρ
+
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
σ
d
a
d
τ
∫
0
T
.
According to expression of V
*
(
t
), we define a nonlinear positive operator
(29)
Φ
V
t
=
μ
*
δ
η
e
-
∫
0
t
μ
2
(
τ
)
d
τ
·
∫
0
t
M
(
τ
)
e
∫
0
τ
μ
2
(
ρ
)
d
ρ
∫
0
r
m
f
(
a
)
·
∫
0
a
V
(
τ
-
σ
)
·
e
-
∫
0
σ
η
V
τ
-
ρ
+
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
σ
d
a
d
τ
+
1
e
∫
0
T
μ
2
ρ
d
ρ
-
1
·
∫
0
T
M
τ
e
∫
0
τ
μ
2
ρ
d
ρ
∫
0
r
m
f
a
·
∫
0
a
V
(
τ
-
σ
)
∫
0
T
·
e
-
∫
0
σ
η
V
τ
-
ρ
+
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
σ
d
a
d
τ
on X
2
,
T
+.IfΦ has a nontrivial fixed point V
*
=
Φ
(
V
*
)
,
V
*
∈
X
2
,
T
+
-
{
0
}, then from (26) there is a T-periodic solution i
*
(
a
,
t
)
∈
X
1
,
T
+
-
{
0
}. So, the system (16) and (17) has an endemic T-periodic solution in a weak sense.In the following, we investigate such a fixed pointV
*
=
Φ
(
V
*
), V
*
∈
X
2
,
T
+
-
{
0
}. First, we define a positive bounded linear operator Ψ
:
X
2
,
T
+
→
X
2
,
T
+ as
(30)
Ψ
V
t
=
μ
*
δ
η
e
-
∫
0
t
μ
2
τ
d
τ
·
∫
0
t
M
(
τ
)
e
∫
0
τ
μ
2
(
ρ
)
d
ρ
∫
0
r
m
f
(
a
)
·
∫
0
a
V
τ
-
σ
e
-
∫
0
σ
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
σ
d
a
d
τ
+
1
e
∫
0
T
μ
2
ρ
d
ρ
-
1
·
∫
0
T
M
(
τ
)
e
∫
0
τ
μ
2
(
ρ
)
d
ρ
∫
0
r
m
f
(
a
)
·
∫
0
a
V
τ
-
σ
e
-
∫
0
σ
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
σ
d
a
d
τ
which is the Fréchet derivative of operator Φ at V
=
0, and it is a majorant of Φ; that is, Φ
≤
Ψ.Next we prove the following lemmas.Lemma 2.
The operatorΦ is monotone nondecreasing and uniformly bounded for V
∈
X
2
,
T
+.Proof.
ForV
∈
X
2
,
T
+, from (29) we obtain
(31)
Φ
V
t
=
μ
*
δ
e
-
∫
0
t
μ
2
(
τ
)
d
τ
·
∫
0
t
M
(
τ
)
e
∫
0
τ
μ
2
(
ρ
)
d
ρ
∫
0
r
m
f
(
a
)
·
1
-
e
-
∫
0
a
η
V
τ
-
ρ
+
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
a
d
τ
-
∫
0
t
‍
M
τ
e
∫
0
τ
‍
μ
2
ρ
d
ρ
∫
0
r
m
‍
f
a
·
∫
0
a
‍
γ
(
τ
-
σ
,
a
-
σ
)
·
e
-
∫
0
σ
‍
η
V
τ
-
ρ
+
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
σ
d
a
d
τ
∫
0
a
+
μ
*
δ
e
-
∫
0
t
‍
μ
2
τ
d
τ
e
∫
0
T
‍
μ
2
ρ
d
ρ
-
1
·
∫
0
T
M
(
τ
)
e
∫
0
τ
μ
2
(
ρ
)
d
ρ
∫
0
r
m
f
(
a
)
·
1
-
e
-
∫
0
a
η
V
τ
-
ρ
+
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
a
d
τ
-
∫
0
T
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
ρ
d
ρ
∫
0
r
m
‍
f
(
a
)
·
∫
0
a
‍
γ
(
τ
-
σ
,
a
-
σ
)
·
e
-
∫
0
σ
‍
η
V
τ
-
ρ
+
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
σ
d
a
d
τ
∫
0
T
.
From (31), we have that Φ is monotone nondecreasing and obtain
(32)
Φ
V
t
≤
μ
*
δ
e
-
∫
0
t
μ
2
τ
d
τ
·
∫
0
t
M
τ
e
∫
0
τ
‍
μ
2
ρ
d
ρ
∫
0
r
m
f
a
d
a
d
τ
+
μ
*
δ
e
-
∫
0
t
μ
2
τ
d
τ
e
∫
0
T
μ
2
ρ
d
ρ
-
1
·
∫
0
T
M
τ
e
∫
0
τ
μ
2
ρ
d
ρ
∫
0
r
m
f
a
d
a
d
τ
.
Since μ
*
∫
0
r
m
‍
f
(
a
)
d
a
=
1, similar to (15), we have
(33)
Φ
V
t
≤
2
δ
e
∫
0
T
μ
2
(
ρ
)
d
ρ
-
1
∫
0
T
M
(
τ
)
e
∫
0
τ
μ
2
(
ρ
)
d
ρ
d
τ
<
+
∞
.
From (33), we obtain that Φ is uniformly bounded.Lemma 3.
Letρ
(
Ψ
) be the spectral radius of operator Ψ. If ρ
(
Ψ
)
>
0, then it is a positive eigenvalue of Ψ associated with a positive eigenvector V
1
∈
X
2
,
T
+
-
{
0
}.Proof.
It is easy to get thatΨ is a linear map from X
2
,
T
+ into itself and leaves the cone invariant.
Let(
Ψ
V
)
(
t
)
=
(
Ψ
1
V
)
(
t
)
+
(
Ψ
2
V
)
(
t
), where
(34)
Ψ
1
V
t
=
μ
*
δ
η
e
-
∫
0
t
μ
2
τ
d
τ
·
∫
0
t
M
(
τ
)
e
∫
0
τ
μ
2
(
ρ
)
d
ρ
∫
0
r
m
f
(
a
)
·
∫
0
a
V
τ
-
σ
e
-
∫
0
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
σ
d
a
d
τ
,
Ψ
2
V
t
=
μ
*
δ
η
e
-
∫
0
t
μ
2
τ
d
τ
e
∫
0
T
μ
2
ρ
d
ρ
-
1
·
∫
0
T
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
f
(
a
)
·
∫
0
a
V
τ
-
σ
e
-
∫
0
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
σ
d
a
d
τ
.
(
Ψ
1
V
)
(
t
) can be rewritten as
(35)
Ψ
1
V
t
=
μ
*
δ
η
e
-
∫
0
t
‍
μ
2
τ
d
τ
·
∫
0
t
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
‍
∫
σ
r
m
‍
f
(
a
)
·
V
τ
-
σ
e
-
∫
0
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
a
d
σ
d
τ
=
μ
*
δ
η
e
-
∫
0
t
‍
μ
2
τ
d
τ
·
∫
0
t
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
τ
-
r
m
τ
‍
V
(
σ
)
∫
τ
-
σ
r
m
‍
f
(
a
)
·
e
-
∫
0
τ
-
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
a
d
σ
d
τ
.
We extend the domain of l
(
a
) and define l
(
a
)
=
0 for a
<
0 or a
≥
r
m; then we have
(36)
μ
*
δ
η
e
-
∫
0
t
‍
μ
2
(
τ
)
d
τ
∫
0
t
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
-
∞
τ
-
r
m
‍
V
(
σ
)
∫
τ
-
σ
r
m
‍
f
(
a
)
·
e
-
∫
0
τ
-
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
a
d
σ
d
τ
=
0
,
μ
*
δ
η
e
-
∫
0
t
‍
μ
2
(
τ
)
d
τ
∫
0
t
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
τ
+
∞
‍
V
(
σ
)
∫
τ
-
σ
r
m
‍
f
(
a
)
·
e
-
∫
0
τ
-
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
a
d
σ
d
τ
=
0
.
From (35) and (36), we have
(37)
∫
τ
-
r
m
τ
‍
V
σ
∫
τ
-
σ
r
m
‍
f
a
e
-
∫
0
τ
-
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
a
d
σ
=
∫
-
∞
+
∞
‍
V
σ
∫
τ
-
σ
r
m
‍
f
a
e
-
∫
0
τ
-
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
a
d
σ
=
∑
n
=
-
∞
+
∞
‍
∫
n
T
n
+
1
T
‍
V
σ
∫
τ
-
σ
r
m
‍
f
a
e
-
∫
0
τ
-
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
a
d
σ
=
∑
n
=
-
∞
+
∞
‍
∫
0
T
‍
V
σ
+
n
T
∫
τ
-
n
T
-
σ
r
m
‍
f
a
·
e
-
∫
0
τ
-
n
T
-
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
a
d
σ
=
∫
0
T
‍
V
(
σ
)
∑
n
=
-
∞
+
∞
‍
∫
τ
-
n
T
-
σ
r
m
‍
f
(
a
)
e
-
∫
0
τ
-
n
T
-
σ
‍
γ
(
τ
-
ρ
,
a
-
ρ
)
d
ρ
d
a
d
σ
.
Combining (35) and (37), we obtain
(38)
Ψ
1
V
t
=
μ
*
δ
η
e
-
∫
0
t
‍
μ
2
τ
d
τ
·
∫
0
t
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
·
∫
0
T
‍
V
(
σ
)
·
∑
n
=
-
∞
+
∞
‍
∫
τ
-
n
T
-
σ
r
m
‍
f
a
·
e
-
∫
0
τ
-
n
T
-
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
a
d
σ
d
τ
=
∫
0
T
‍
K
1
(
t
,
σ
)
V
(
σ
)
d
σ
,
where
(39)
K
1
t
,
σ
=
μ
*
δ
η
e
-
∫
0
t
‍
μ
2
τ
d
τ
·
∫
0
t
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
·
∑
n
=
-
∞
+
∞
‍
∫
τ
-
n
T
-
σ
r
m
‍
f
a
e
-
∫
0
τ
-
n
T
-
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
a
d
τ
.
K
1
(
t
,
σ
) is well defined and is uniformly bounded; the sum in the expression of K
1
(
t
,
σ
) is a finite sum due to f
(
a
)
=
0, a
≤
τ
-
n
T
-
σ
<
0 or a
≥
τ
-
n
T
-
σ
>
r
m.
Hence,Ψ
1 can be regarded as an operator on L
1
[
0
,
T
]. From the well-known compactness criteria in L
1 [20], we obtain that Ψ
1 is compact.
Similar toΨ
1, we get the expression of Ψ
2 as
(40)
Ψ
2
V
t
=
∫
0
T
‍
K
2
t
,
σ
V
σ
d
σ
,
where
(41)
K
2
t
,
σ
=
μ
*
δ
η
e
-
∫
0
t
‍
μ
2
τ
d
τ
e
∫
0
T
‍
μ
2
ρ
d
ρ
-
1
·
∫
0
T
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
·
∑
n
=
-
∞
+
∞
‍
∫
τ
-
n
T
-
σ
r
m
‍
f
a
e
-
∫
0
τ
-
n
T
-
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
a
d
τ
,
K
2
(
t
,
σ
) is well defined and is uniformly bounded, and Ψ
2 is compact in L
1
[
0
,
T
].
Since(
Ψ
V
)
(
t
)
=
(
Ψ
1
V
)
(
t
)
+
(
Ψ
2
V
)
(
t
), we obtain that Ψ is positive, linear, and compact. If ρ
(
Ψ
)
>
0, from the Krein-Rutman theorem [21], there is V
1
∈
X
2
,
T
+
-
{
0
} such that
(42)
Ψ
V
1
t
=
ρ
Ψ
V
1
t
holds.Using the above lemmas, we have the following theorem.Theorem 4.
Ifρ
(
Ψ
)
>
1 holds, then operator Φ has at least one nontrivial fixed point V
*:
(43)
V
*
=
Φ
V
*
,
V
*
∈
X
2
,
T
*
-
0
;
that is, the system (16) and (17) has at least an endemic T-periodic solution.Proof.
From Lemma3, we have
(44)
Ψ
V
1
t
=
ρ
Ψ
V
1
t
=
∫
0
T
‍
K
(
t
,
σ
)
V
1
(
σ
)
d
σ
,
V
1
(
t
)
∈
X
2
,
T
*
-
{
0
}
,
where K
(
t
,
σ
)
=
K
1
(
t
,
σ
)
+
K
2
(
t
,
σ
). Since K
1
(
t
,
σ
)
,
K
2
(
t
,
σ
) are uniformly bounded, let
(45)
K
¯
=
sup
⁡
0
≤
t
≤
T
,
0
≤
σ
≤
r
m
K
t
,
σ
<
+
∞
.
We have
(46)
ρ
Ψ
V
1
t
≤
K
¯
∫
0
T
‍
V
1
(
σ
)
d
σ
=
K
¯
V
1
(
t
)
X
2
,
T
.
Let
(47)
V
*
t
=
ρ
Ψ
ln
⁡
ρ
Ψ
η
K
¯
r
m
V
1
(
t
)
V
1
(
t
)
X
2
,
T
∈
X
2
,
T
*
-
{
0
}
,
whose positivity follows from ρ
(
Ψ
)
>
1. From (29), (44), and (46), we have
(48)
Φ
V
*
t
=
μ
*
δ
η
e
-
∫
0
t
‍
μ
2
(
τ
)
d
τ
·
∫
0
t
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
‍
f
(
a
)
·
∫
0
a
‍
V
*
(
τ
-
σ
)
·
e
-
∫
0
σ
‍
η
V
*
τ
-
ρ
+
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
σ
d
a
d
τ
+
1
e
∫
0
T
‍
μ
2
ρ
d
ρ
-
1
·
∫
0
T
‍
M
τ
e
∫
0
τ
‍
μ
2
ρ
d
ρ
∫
0
r
m
‍
f
a
·
∫
0
a
‍
V
*
(
τ
-
σ
)
·
e
-
∫
0
σ
‍
(
η
V
*
(
τ
-
ρ
)
+
γ
(
τ
-
ρ
,
a
-
ρ
)
)
d
ρ
d
σ
d
a
d
τ
∫
0
t
≥
μ
*
δ
η
e
-
∫
0
t
‍
μ
2
τ
d
τ
·
∫
0
t
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
‍
f
(
a
)
·
∫
0
a
‍
V
*
(
τ
-
σ
)
·
e
-
ln
⁡
ρ
Ψ
/
K
¯
r
m
V
1
t
X
2
,
T
∫
0
r
m
‍
ρ
Ψ
V
1
τ
-
ρ
d
ρ
·
e
-
∫
0
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
σ
d
a
d
τ
+
1
e
∫
0
T
‍
μ
2
ρ
d
ρ
-
1
·
∫
0
T
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
‍
f
(
a
)
·
∫
0
a
‍
V
*
(
τ
-
σ
)
·
e
-
ln
⁡
ρ
Ψ
/
K
¯
r
m
V
1
t
X
2
,
T
∫
0
r
m
‍
ρ
Ψ
V
1
τ
-
ρ
d
ρ
·
e
-
∫
0
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
d
σ
d
a
d
τ
∫
0
t
≥
e
-
ln
⁡
ρ
(
Ψ
)
(
Ψ
V
*
)
(
t
)
=
V
*
(
t
)
.
From Lemma 2, we obtain that operator Φ is monotone nondecreasing. So, we define a monotone sequence
(49)
V
n
+
1
=
Φ
(
V
n
)
,
V
1
≤
V
2
≤
⋯
≤
V
n
≤
⋯
.
From Lemma 2, we have that sequence {
V
n
} is bounded above. From B. Levi’s theorem, there exists V
*
(
t
)
∈
X
2
,
T
*
-
{
0
} such that lim
n
→
+
∞
⁡
V
n
(
t
)
=
V
*
(
t
) and V
*
(
t
)
=
(
Φ
V
*
)
(
t
).


## 5. Uniqueness of Endemic Periodic Solution

If(
i
1
*
(
a
,
t
)
,
V
1
*
(
t
)
)
,
(
i
2
*
(
a
,
t
)
,
V
2
*
(
t
)
) are the endemic T-periodic solutions of the system (16) and (17), respectively. Let
(50)
G
j
(
t
)
=
δ
μ
*
∫
0
r
m
‍
f
a
i
j
*
a
,
t
d
a
M
(
t
)
,
j
=
1,2
.
Since M
(
t
)
,
i
(
a
,
t
), j
=
1,2 are T-periodic functions, then G
j
(
t
), j
=
1,2, are periodic functions. From (16), we have
(51)
V
j
*
(
t
)
=
V
0
,
j
e
-
∫
0
t
‍
μ
2
(
σ
)
d
σ
+
e
-
∫
0
t
‍
μ
2
(
σ
)
d
σ
∫
0
t
‍
G
j
(
τ
)
e
∫
0
τ
‍
μ
2
(
σ
)
d
σ
d
τ
,
where V
0
,
j
=
∫
0
T
‍
∫
0
r
m
‍
f
a
i
*
a
,
τ
M
τ
e
∫
0
τ
‍
μ
2
ρ
d
ρ
d
a
d
τ
/
e
∫
0
T
‍
μ
2
(
ρ
)
d
ρ
-
1
=
∫
0
T
‍
G
j
(
τ
)
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
d
τ/e
∫
0
T
‍
μ
2
(
ρ
)
d
ρ
-
1, j
=
1,2.Considering thatG
j
(
t
)
,
μ
2
(
t
) are T-periodic functions, we obtain
(52)
0
<
a
1
j
=
e
-
∫
0
T
‍
μ
2
(
σ
)
d
σ
∫
0
T
‍
G
j
(
τ
)
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
d
τ
e
∫
0
T
‍
μ
2
(
ρ
)
d
ρ
-
1
≤
V
j
*
(
t
)
≤
∫
0
T
‍
G
j
τ
M
τ
e
∫
0
τ
‍
μ
2
ρ
d
ρ
d
τ
e
∫
0
T
‍
μ
2
ρ
d
ρ
-
1
+
∫
0
T
‍
G
j
(
τ
)
e
∫
0
τ
‍
μ
2
(
σ
)
d
σ
d
τ
=
a
2
j
<
+
∞
,
j
=
1,2
.
For a number ν
∈
(
0,1
) and the nonlinear positive operator Φ, we have
(53)
Φ
ν
V
j
*
t
-
ν
Φ
V
j
*
t
=
μ
*
δ
η
ν
e
-
∫
0
t
‍
μ
2
τ
d
τ
·
∫
0
t
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
‍
f
(
a
)
·
∫
0
a
‍
V
j
*
τ
-
σ
e
-
∫
0
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
·
e
-
η
ν
∫
0
σ
‍
V
j
*
τ
-
ρ
d
ρ
-
e
-
η
∫
0
σ
‍
V
J
*
τ
-
ρ
d
ρ
d
σ
d
a
d
τ
+
1
e
∫
0
T
‍
μ
2
ρ
d
ρ
-
1
·
∫
0
T
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
‍
f
(
a
)
·
∫
0
a
‍
V
j
*
τ
-
σ
e
-
∫
0
σ
‍
γ
τ
-
ρ
,
a
-
ρ
d
ρ
·
e
-
η
ν
∫
0
σ
‍
V
j
*
τ
-
ρ
d
ρ
-
e
-
η
∫
0
σ
‍
V
j
*
τ
-
ρ
d
ρ
d
σ
d
a
d
τ
∫
0
t
=
μ
*
δ
η
ν
-
∫
0
t
‍
μ
2
τ
d
τ
·
∫
0
t
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
‍
f
(
a
)
·
∫
0
a
‍
V
j
*
(
τ
-
σ
)
·
e
-
∫
0
σ
‍
η
V
j
*
τ
-
ρ
+
γ
τ
-
ρ
,
a
-
ρ
d
ρ
·
e
η
1
-
ν
∫
0
σ
‍
V
j
*
τ
-
ρ
d
ρ
-
1
d
σ
d
a
d
τ
+
1
e
∫
0
T
‍
μ
2
ρ
d
ρ
-
1
·
∫
0
T
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
‍
f
(
a
)
·
∫
0
a
‍
V
j
*
τ
-
σ
·
e
-
∫
0
σ
‍
η
V
j
*
τ
-
ρ
+
γ
τ
-
ρ
,
a
-
ρ
d
ρ
·
e
η
1
-
ν
∫
0
σ
‍
V
j
*
τ
-
ρ
d
ρ
-
1
d
σ
d
a
d
τ
∫
0
t
.
From (52), we have
(54)
e
η
(
1
-
ν
)
∫
0
σ
‍
V
j
*
(
τ
-
ρ
)
d
ρ
-
1
≥
e
η
(
1
-
ν
)
a
1
j
σ
-
1
.
Combine (53) and (54) to get
(55)
Φ
(
ν
V
j
*
)
(
t
)
-
ν
Φ
V
j
*
t
≥
c
j
>
0
,
j
=
1,2
,
where
(56)
c
j
=
μ
*
δ
η
ν
-
∫
0
t
‍
μ
2
τ
d
τ
·
∫
0
t
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
‍
f
(
a
)
·
∫
0
a
‍
V
j
*
(
τ
-
σ
)
·
e
-
∫
0
σ
‍
η
V
j
*
τ
-
ρ
+
γ
τ
-
ρ
,
a
-
ρ
d
ρ
e
η
1
-
ν
a
1
j
σ
-
1
d
σ
d
a
d
τ
+
1
e
∫
0
T
‍
μ
2
ρ
d
ρ
-
1
·
∫
0
T
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
‍
f
(
a
)
·
∫
0
a
‍
V
j
*
(
τ
-
σ
)
·
e
-
∫
0
σ
‍
η
V
j
*
τ
-
ρ
+
γ
τ
-
ρ
,
a
-
ρ
d
ρ
·
e
η
1
-
ν
a
1
j
σ
-
1
d
σ
d
a
d
τ
∫
0
t
.
Now we use (52) and (55) to investigate the uniqueness of endemic periodic solution.From (52) we have
(57)
V
1
*
(
t
)
≥
a
11
=
a
11
a
22
-
1
a
22
≥
a
11
a
22
-
1
V
2
*
(
t
)
.
If a
11
a
22
-
1
<
1, we take ν
1
=
a
11
a
22
-
1. From (55) we have
(58)
Φ
ν
1
V
j
*
t
≥
ν
1
Φ
V
j
*
t
+
c
j
(
1
)
,
c
j
(
1
)
>
0
,
j
=
1,2
.
Noting V
j
*
(
t
)
=
(
Φ
V
j
*
)
(
t
), j
=
1,2, from Lemma 2 and (52) we obtain
(59)
V
1
*
t
=
(
Φ
V
1
*
)
(
t
)
≥
(
Φ
(
ν
1
V
2
*
)
)
(
t
)
≥
ν
1
(
Φ
V
2
*
)
(
t
)
+
c
2
(
1
)
=
ν
1
V
2
*
(
t
)
+
c
2
(
1
)
a
22
-
1
a
22
≥
ν
1
V
2
*
(
t
)
+
c
2
(
1
)
a
22
-
1
V
2
*
(
t
)
=
(
ν
1
+
c
2
(
1
)
a
22
-
1
)
V
2
*
(
t
)
.
If ν
1
+
c
2
(
1
)
a
22
-
1
<
1, we take ν
2
=
ν
1
+
c
2
(
1
)
a
22
-
1 and repeat the above process to have
(60)
V
1
*
t
≥
ν
2
+
c
2
1
a
22
-
1
V
2
*
(
t
)
=
ν
1
+
c
2
1
a
22
-
1
+
c
2
1
a
22
-
1
V
2
*
t
.
Since lim
m
→
+
∞
⁡
(
ν
1
+
m
c
2
(
1
)
a
22
-
1
)
=
+
∞, then there exists a number m
0
∈
N
+ such that ν
1
+
m
0
c
2
(
1
)
a
22
-
1
≥
1, which implies that
(61)
V
1
*
(
t
)
≥
ν
1
+
m
0
c
2
1
a
22
-
1
V
2
*
t
.
So, if a
11
a
22
-
1
≥
1 from (57) and if a
11
a
22
-
1
<
1 from (61), we always have
(62)
V
1
*
t
≥
V
2
*
t
.
Exchanging the role of V
1
*
(
t
) and V
2
*
(
t
), similar to the above arguments, we can obtain
(63)
V
2
*
t
≥
a
12
=
a
12
a
21
-
1
a
21
≥
a
12
a
21
-
1
V
1
*
t
,
(64)
V
2
*
t
≤
V
1
*
t
.
Equations (62) and (64) imply that V
1
*
(
t
)
=
V
2
*
(
t
). From (26), we have i
1
*
(
a
,
t
)
=
i
2
*
(
a
,
t
).From Theorem4 and the above arguments, we obtain the theorem.Theorem 5.
Ifρ
(
Ψ
)
>
1 holds, the system (16) and (17) has a unique endemic T-periodic solution in a weak sense.


## 6. Stability of the Disease-Free Steady State

The system (16) always has the infection-free steady state:
(65)
E
0
=
i
*
,
V
*
=
0,0
.
In this section, we will prove that the disease-free steady state is global attractiveness. From the first equation of (16), we have
(66)
i
(
a
,
t
)
=
η
∫
0
a
‍
e
-
∫
σ
a
‍
η
V
t
-
a
+
τ
+
γ
τ
,
t
-
a
+
τ
d
τ
V
t
-
a
+
σ
d
σ
,
vvvvvvvvvvvvvvvvvvvvvvvvvvv
t
-
a
>
0
,
e
-
∫
0
t
‍
(
η
V
(
a
-
t
+
τ
)
+
γ
(
τ
,
a
-
t
+
τ
)
)
d
τ
i
0
(
a
-
t
)
+
η
∫
0
t
‍
e
-
∫
σ
t
‍
(
η
V
(
a
-
t
+
τ
)
+
γ
(
τ
,
a
-
t
+
τ
)
)
d
τ
V
·
(
a
-
t
+
σ
)
d
σ
,
a
-
t
≥
0
.
From the second equation of (16), we have
(67)
V
t
=
μ
*
δ
e
-
∫
0
t
‍
μ
2
τ
d
τ
·
∫
0
t
‍
∫
0
r
m
‍
f
a
i
a
,
τ
M
τ
e
∫
0
τ
‍
μ
2
ρ
d
ρ
d
a
d
τ
+
V
0
=
μ
*
δ
e
-
∫
0
t
‍
μ
2
τ
d
τ
·
∫
0
t
‍
∫
0
τ
‍
f
a
i
a
,
τ
M
τ
e
∫
0
τ
‍
μ
2
ρ
d
ρ
d
a
d
τ
+
∫
τ
r
m
‍
f
a
i
a
,
τ
M
τ
e
∫
0
τ
‍
μ
2
ρ
d
ρ
d
a
d
τ
+
V
0
∫
0
t
.Substituting (66) into (67), we have
(68)
V
t
=
μ
*
δ
e
-
∫
0
t
‍
μ
2
τ
d
τ
F
1
t
+
F
2
t
+
F
3
t
+
V
0
,
where
(69)
F
1
t
=
η
∫
0
t
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
ρ
d
ρ
∫
0
τ
‍
f
(
a
)
·
∫
0
a
‍
e
-
∫
σ
a
‍
(
η
V
(
τ
-
a
+
ξ
)
+
γ
(
ξ
,
τ
-
a
+
ξ
)
)
d
ξ
·
V
τ
-
a
+
σ
d
σ
d
a
d
τ
,
F
2
t
=
∫
0
t
‍
M
τ
e
∫
0
τ
‍
μ
2
ρ
d
ρ
∫
τ
r
m
‍
f
a
·
e
-
∫
0
τ
‍
η
V
a
-
τ
+
ξ
+
γ
ξ
,
a
-
τ
+
ξ
d
ξ
i
0
a
-
τ
d
a
d
τ
,
F
3
t
=
η
∫
0
t
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
ρ
d
ρ
∫
τ
r
m
‍
f
(
a
)
·
∫
0
τ
‍
e
-
∫
σ
τ
‍
(
η
V
(
a
-
τ
+
ξ
)
+
γ
(
ξ
,
a
-
τ
+
ξ
)
)
d
ξ
·
V
a
-
τ
+
σ
d
σ
d
a
d
τ
.
For t
>
r
m, by variable substitution, F
1
(
t
)
,
F
2
(
t
)
,
F
3
(
t
) can be rewritten as
(70)
F
1
t
=
η
∫
0
r
m
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
ρ
d
ρ
∫
0
τ
‍
f
(
a
)
·
∫
0
a
‍
e
-
∫
0
σ
‍
(
η
V
(
τ
-
ξ
)
+
γ
(
a
-
ξ
,
τ
-
ξ
)
)
d
ξ
·
V
τ
-
σ
d
σ
d
a
d
τ
+
H
(
t
)
,
where
(71)
H
t
=
η
∫
r
m
t
‍
M
τ
e
∫
0
τ
‍
μ
2
ρ
d
ρ
∫
0
τ
‍
f
a
·
∫
0
a
‍
e
-
∫
0
σ
‍
(
η
V
(
τ
-
ξ
)
+
γ
(
a
-
ξ
,
τ
-
ξ
)
)
d
ξ
·
V
τ
-
σ
d
σ
d
a
d
τ
,
F
2
t
=
∫
0
r
m
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
ρ
d
ρ
∫
τ
r
m
‍
f
(
a
)
·
e
-
∫
0
τ
‍
η
V
a
-
τ
+
ξ
+
γ
ξ
,
a
-
τ
+
ξ
d
ξ
i
0
a
-
τ
d
a
d
τ
+
∫
r
m
t
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
ρ
d
ρ
∫
τ
r
m
‍
f
(
a
)
·
e
-
∫
0
τ
‍
η
V
a
-
τ
+
ξ
+
γ
ξ
,
a
-
τ
+
ξ
d
ξ
i
0
a
-
τ
d
a
d
τ
,
F
3
t
=
η
∫
0
r
m
‍
M
τ
e
∫
0
τ
‍
μ
2
ρ
d
ρ
∫
τ
r
m
‍
f
a
·
∫
0
τ
‍
e
-
∫
0
σ
‍
η
V
a
-
ξ
+
γ
τ
-
ξ
,
a
-
ξ
d
ξ
·
V
a
-
σ
d
σ
d
a
d
τ
+
η
∫
r
m
t
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
τ
r
m
‍
f
(
a
)
·
∫
0
τ
e
-
∫
0
σ
‍
η
V
a
-
ξ
+
γ
τ
-
ξ
,
a
-
ξ
d
ξ
V
a
-
σ
d
σ
d
a
d
τ
.
Noting f
(
a
)
=
0, a
≥
r
m, for τ
≥
r
m we have
(72)
∫
0
τ
‍
f
a
∫
0
a
‍
e
-
∫
0
σ
‍
η
V
τ
-
ξ
+
γ
a
-
ξ
,
τ
-
ξ
d
ξ
V
τ
-
σ
d
σ
d
a
=
∫
0
r
m
‍
f
a
∫
0
a
‍
e
-
∫
0
σ
‍
η
V
τ
-
ξ
+
γ
a
-
ξ
,
τ
-
ξ
d
ξ
V
τ
-
σ
d
σ
d
a
,
∫
τ
r
m
‍
f
(
a
)
e
-
∫
0
τ
‍
(
η
V
(
a
-
τ
+
ξ
)
+
γ
(
ξ
,
a
-
τ
+
ξ
)
)
d
ξ
i
0
(
a
-
τ
)
d
a
=
0
,
∫
τ
r
m
‍
f
(
a
)
∫
0
τ
‍
e
-
∫
0
σ
‍
(
η
V
(
a
-
ξ
)
+
γ
(
τ
-
ξ
,
a
-
ξ
)
)
d
ξ
V
(
a
-
σ
)
d
σ
d
a
=
0
.
From (15) we obtain that V
(
t
) is bounded above. Moreover, from (72) we have
(73)
lim
t
→
+
∞
μ
*
δ
e
-
∫
0
t
‍
μ
2
τ
d
τ
F
1
t
-
H
t
=
0
,
lim
t
→
+
∞
μ
*
δ
e
-
∫
0
t
‍
μ
2
τ
d
τ
F
2
t
=
0
,
lim
t
→
+
∞
μ
*
δ
e
-
∫
0
t
‍
μ
2
τ
d
τ
F
3
t
=
0
.
Before proving the global attractiveness of the disease-free steady state, we first prove a lemma.Lemma 6.
Consider(74)
lim
sup
⁡
t
→
+
∞
μ
*
δ
H
t
≤
μ
*
δ
η
V
∞
e
∫
0
T
‍
μ
2
ρ
d
ρ
-
1
·
∫
0
T
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
‍
f
(
a
)
·
∫
0
a
‍
e
-
∫
0
σ
‍
γ
a
-
ξ
,
τ
-
ξ
d
ξ
d
σ
d
a
d
τ
,
where V
∞
=
lim
sup
⁡
t
→
+
∞
⁡
V
(
t
).Proof.
For any givenε
>
0, we choose n
0
∈
N
+ such that V
(
t
)
<
V
∞
+
ε, for t
>
n
0
T
-
r
m. Let t
=
n
T
+
t
′, 0
≤
t
′
<
T, n
>
n
0; we have
(75)
μ
*
δ
H
t
≤
μ
*
δ
η
e
-
∫
0
n
T
‍
μ
2
τ
d
τ
·
∫
r
m
(
n
+
1
)
T
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
‍
f
(
a
)
·
∫
0
a
‍
e
-
∫
0
σ
‍
γ
a
-
ξ
,
τ
-
ξ
d
ξ
V
τ
-
σ
d
σ
d
a
d
τ
≤
μ
*
δ
η
e
-
∫
0
n
T
‍
μ
2
τ
d
τ
·
∫
r
m
n
0
T
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
a
‍
f
(
a
)
e
-
∫
0
σ
‍
γ
(
a
-
ξ
,
τ
-
ξ
)
d
ξ
·
V
τ
-
σ
d
σ
d
a
d
τ
+
sup
⁡
n
0
T
-
r
m
≤
t
<
+
∞
V
(
t
)
·
∫
n
0
T
(
n
+
1
)
T
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
‍
f
(
a
)
·
∫
0
a
‍
e
-
∫
0
σ
‍
γ
a
-
ξ
,
τ
-
ξ
d
ξ
d
σ
d
a
d
τ
≤
μ
*
δ
η
e
-
∫
0
n
T
‍
μ
2
τ
d
τ
·
∫
r
m
n
0
T
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
a
‍
f
(
a
)
e
-
∫
0
σ
‍
γ
(
a
-
ξ
,
τ
-
ξ
)
d
ξ
·
V
τ
-
σ
d
σ
d
a
d
τ
+
V
∞
+
ε
·
∫
0
(
n
+
1
)
T
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
‍
f
(
a
)
·
∫
0
a
‍
e
-
∫
0
σ
‍
γ
a
-
ξ
,
τ
-
ξ
d
ξ
d
σ
d
a
d
τ
.
Since M
(
t
)
,
μ
2
(
t
)
,
γ
(
a
,
t
) are T-periodic functions, we obtain
(76)
e
-
∫
0
n
T
‍
μ
2
(
τ
)
d
τ
=
e
-
n
∫
0
T
‍
μ
2
(
τ
)
d
τ
,
(77)
∫
0
n
+
1
T
‍
M
τ
e
∫
0
τ
‍
μ
2
ρ
d
ρ
∫
0
r
m
‍
f
a
·
∫
0
a
‍
e
-
∫
0
σ
‍
γ
a
-
ξ
,
τ
-
ξ
d
ξ
d
σ
d
a
d
τ
=
e
n
∫
0
T
‍
μ
2
ρ
d
ρ
-
1
e
∫
0
T
‍
μ
2
ρ
d
ρ
-
1
·
∫
0
T
‍
M
τ
e
∫
0
τ
‍
μ
2
ρ
d
ρ
∫
0
r
m
‍
f
a
·
∫
0
a
‍
e
-
∫
0
σ
‍
γ
a
-
ξ
,
τ
-
ξ
d
ξ
d
σ
d
a
d
τ
.
From (75), (76), we take the limit supreme on both sides of (75) when t
→
+
∞  (
n
→
+
∞
) and have
(78)
lim
sup
⁡
t
→
+
∞
μ
*
δ
H
t
≤
μ
*
δ
η
V
∞
+
ε
e
∫
0
T
‍
μ
2
ρ
d
ρ
-
1
·
∫
0
T
‍
M
(
τ
)
e
∫
0
τ
‍
μ
2
(
ρ
)
d
ρ
∫
0
r
m
‍
f
(
a
)
·
∫
0
a
‍
e
-
∫
0
σ
‍
γ
a
-
ξ
,
τ
-
ξ
d
ξ
d
σ
d
a
d
τ
.
Since ε is arbitrary, the lemma follows.
Let(79)
R
0
T
=
μ
*
δ
η
e
∫
0
T
‍
μ
2
ρ
d
ρ
-
1
·
∫
0
T
‍
M
τ
e
∫
0
τ
‍
μ
2
ρ
d
ρ
∫
0
r
m
‍
f
a
·
∫
0
a
‍
e
-
∫
0
σ
‍
γ
a
-
ξ
,
τ
-
ξ
d
ξ
d
σ
d
a
d
τ
.
Using Lemma 6, we have the following theorem.Theorem 7.
IfR
0
(
T
)
<
1 holds, then the infection-free steady state E
0 is global attractiveness; that is,
(80)
lim
t
→
+
∞
(
i
(
a
,
t
)
,
V
(
t
)
)
=
(
0,0
)
.Proof.
From (73) and Lemma 6, we take the limit supreme on both sides of (68) when t
→
+
∞  (
n
→
+
∞
) and have
(81)
V
∞
=
lim
sup
⁡
t
→
+
∞
V
t
=
lim
sup
⁡
t
→
+
∞
μ
*
δ
e
-
∫
0
t
‍
μ
2
τ
d
τ
·
F
1
t
+
F
2
t
+
F
3
t
+
V
0

δ
e
-
∫
0
t
‍
μ
2
τ
d
τ
=
lim
sup
⁡
t
→
+
∞
μ
*
δ
e
-
∫
0
t
‍
μ
2
τ
d
τ
H
t
≤
V
∞
R
0
T
.
If R
0
(
T
)
<
1, then V
∞
=
0, which implies lim
t
→
+
∞
⁡
V
(
t
)
=
0. When t
>
r
m
≥
a, from (66) we have
(82)
i
(
a
,
t
)
=
η
∫
0
a
‍
e
-
∫
σ
a
‍
η
V
t
-
a
+
τ
+
γ
τ
,
t
-
a
+
τ
d
τ
V
t
-
a
+
σ
d
σ
.
Since V
∞
=
0, similar to the proof of Lemma 6, for any given ε
>
0, we choose n
0
∈
N
+ such that V
(
t
)
<
ε, for t
>
n
0
T
-
r
m. From (82) to get
(83)
i
(
a
,
t
)
≤
η
r
m
ε
.
Since ε is arbitrary, inequality (83) implies
(84)
lim
t
→
+
∞
i
(
a
,
t
)
=
0
.
This completes our proof.


## 7. Discussion

In this paper, in order to reflect the dependence of vector-borne diseases progress on seasonality of vector population and chronological age of the host population, our SIS epidemic model has periodic parameters, and the host population is structured by the chronological age. So, the SIS epidemic model has a relatively complex form. Due to the difficulty in the analysis, to our knowledge, we do not obtain threshold value for the existence and uniqueness of nontrivial endemic periodic solution of the age-structured SIS epidemic model and do not prove the infection-free state is globally stable and only get the condition under which the infection-free state is global attractiveness. In addition, for our age-dependent case, the stability of the endemic periodic solution(
i
*
(
a
,
t
)
,
V
*
(
t
)
) for ρ
(
Ψ
)
>
1 has been left as an open problem.



---

*Source: 838312-2015-02-24.xml*