# Positive Periodic Solutions of an Epidemic Model with Seasonality

Authors: Gui-Quan Sun; Zhenguo Bai; Zi-Ke Zhang; Tao Zhou; Zhen Jin
Journal: The Scientific World Journal
 (2013)

---

## Abstract

An SEI autonomous model with logistic growth rate and its corresponding nonautonomous model are investigated. For the autonomous case, we give the attractive regions of equilibria and perform some numerical simulations. Basic demographic reproduction numberRd is obtained. Moreover, only the basic reproduction number R0 cannot ensure the existence of the positive equilibrium, which needs additional condition Rd>R1. For the nonautonomous case, by introducing the basic reproduction number defined by the spectral radius, we study the uniform persistence and extinction of the disease. The results show that for the periodic system the basic reproduction number is more accurate than the average reproduction number.


---

## Body

## 1. Introduction

Bernoulli was the first person to use mathematical method to evaluate the effectiveness of inoculation for smallpox [1–6]. Then in 1906, Hawer studied the regular occurrence of measles by a discrete-time model. Moreover, Ross [3, 4] adopted the continuous model to study the dynamics of malaria between mosquitoes and humans in 1916 and 1917. In 1927, Kermack and McKendrick [5, 6] extended the above works and established the threshold theory. So far, mathematical models have gotten great development and have been used to study population dynamics, ecology, and epidemic, which can be classified in terms of different aspects. From the aspect of the incidence of infectious diseases, there are bilinear incidence, standard incidence, saturating incidence, and so on. According to the type of demographic import, the constant import, the exponential import, and the logistic growth import are the most common forms. The simple exponential growth models can provide an adequate approximation to population growth for the initial period. If no predation or intraspecific competition for populations is included, the population can continue to increase. However, it is impossible to grow immoderately due to the intraspecific competition for environmental resources such as food and habitat. So, for this case, logistic model is more reasonable and realistic which has been adopted and studied [7–18]. Moreover, due to its rich dynamics, the logistic models have been applied to many fields. Fujikawa et al. [9] applied the logistic model to show Escherichia coli growth. Invernizzi and Terpin [14] used a generalized logistic model to describe photosynthetic growth and predict biomass production. Min et al. [15] used logistic dynamics model to describe coalmining cities' economic growth mechanism and sustainable development. There is a good fit in simulating the coalmining cities’ growth and development track based on resource development cycle. Banaszak et al. [17] investigated logistic models in flexible manufacturing, and Brianzoni et al. [18] studied a business-cycle model with logistic population growth. Muroya [13] investigated discrete models of non-autonomous logistic equations. As a result, this paper builds an SEI ordinary differential model with the logistic growth rate and the standard incidence.For general epidemic models, we mainly study their threshold dynamics, that is, the basic reproduction number which determines whether the disease can invade the susceptible population successfully. However, for the system with logistic growth rate, besides the basic reproduction number, the qualitative dynamics are controlled by a demographic thresholdRd which has a similar meaning and is called as the basic demographic reproduction number. If Rd>1, the population grows; that is, a critical mass of individuals for the disease to spread may be supported. If Rd<1, the population will not survive; that is, not enough mass of individuals may be supported for the disease to spread. For this case, the dynamical behavior of disease will be decided by two thresholds R0 and Rd.It is well known that many diseases exhibit seasonal fluctuations, such as whooping cough, measles, influenza, polio, chickenpox, mumps, and rabies [19–22]. Seasonally effective contact rate [22–26], periodic changing in the birth rate [27], and vaccination program [28] are often regarded as sources of periodicity. Seasonally effective contact rate is related to the behavior of people and animals, the temperature, and the economy. Due to the existence of different seasons, people have different activities which may lead to a different contact rate. Because of various factors, the economy in a different season has a very big difference. Therefore, this paper studies the corresponding non-autonomous system which is obtained by changing the constant transmission rate into the periodic transmission rate. Seasonal transmission is often assumed to be sinusoidal (cosine function has the same meaning), such that λ(t)=λ(1+ηsin(πt/b)) where η is the amplitude of seasonal variation in transmission (typically referred to as the “strength of seasonal forcing”) and 2b is the period, which is a crude assumption for many infectious diseases [29–31]. When η=0, there is no nonseasonal infections. Motivated by biological realism, some recent papers take the contact rate as λ(t)=λ(1+ηterm(t)), where term is a periodic function which is +1 during a period of time and −1 during other time. More natural term can be written as λ(t)=λ(1+η)term(t)  [29]. Here, we take the form β(t)=a[1+bsin(πt/10)].The paper is organized as follows. In Section2, we introduce an autonomous model and analyze the equilibria and their respective attractive region. In Section 3, we study the non-autonomous system in terms of global asymptotic stability of the disease-free equilibrium and the existence of positive periodic solutions. Moreover, numerical simulations are also performed. In Section 4, we give a brief discussion.


## 2. Autonomous Model and Analysis

### 2.1. Model Formulation

The model is a system of SEI ordinary differential equations, whereS is the susceptible, E is the exposed, I is the infected, and N=S+E+I. This system considers the logistic growth rate and the standard incidence which is fit for the long-term growth of many large populations. The incubation period is considered for many diseases which do not develop symptoms immediately and need a period of time to accumulate a pathogen quantity for clinical outbreak, such as rabies, hand-foot-mouth disease, tuberculosis, and AIDS [22, 32]. The model we employ is as follows:
(1)dSdt=rN(1-Nk)-βSIN-mS,dEdt=βSIN-σE-mE,dIdt=σE-mI-μI,
where all parameters are positive whose interpretations can be seen in Table 1.Table 1
Descriptions and values of parameters in model (1).
Parameter
Interpretation
Value
r
The intrinsic growth rate
k
The carrying capacity
100000
β
The transmission rate
m
The natural mortality rate
0.1
σ
Clinical outcome rate
0.2
μ
The disease-induced mortality rate
0.1
a
The baseline contact rate
b
The magnitude of forcingNoticing the equations in model (1), we have
(2)dNdt=rN(1-Nk)-mN-μI.When there exists no disease, we have(3)dNdt=rN(1-Nk)-mN=[r(1-Nk)-m]N.LetRd=r/m, if Rd>1, N→N0=(1-m/r)k for N(0)>0, as t→+∞; that is, the population will grow and tend to a steady state N0. If Rd<1, then dN/dt<0 which will cause the population to disappear. Thus, Rd is the basic demographic reproduction number. From the above equation, the feasible region can be obtained: X={(S,E,I)∣S,E,I≥0,0≤S+E+I≤(1-m/r)k}, where r>m.Theorem 1.
The regionX is positively invariant with respect to system (1).


### 2.2. Dynamical Analysis

Let the right hand of system (1) to be zero; it is easy to see that system (1) has three equilibria:
(4)O=(0,0,0),E0=((1-mr)k,0,0),E*=(S*,E*,I*),
where O is the origin, E0 is the disease-free equilibrium, and E* is the endemic equilibrium. Concretely, one can have
(5)S*=[mβ(m+σ+μ)+μ(m+μ)(m+σ)(R0-1)]N*β[(m+μ)(m+σ)(R0-1)+m(m+σ+μ)],(6)E*=[r(1-N*/k)-m](m+μ)N*σμ,(7)I*=(m+μ)(m+σ)(R0-1)N*β(m+σ+μ),(8)N*=(k[βm(m+σ+μ)(Rd-1)+μ(m+σ)(m+μ)(1-R0)])×(rβ(m+σ+μ))-1,
where Rd=r/m is the basic demographic reproduction number and R0=βσ/(m+σ)(m+μ) is the basic reproduction number which can be obtained by the next-generation matrix method [33–35]. The introduction of the basic demographic reproduction number can be found in [36].Moreover, from (6) and (7), the conditions of the endemic equilibrium to exist are R0>1 and Rd>R1, where R1=1+μ(m+μ)(m+σ)(R0-1)/mβ(m+σ+μ). So, we can obtain the following theorems.Theorem 2.
The system (1) has three equilibria: origin O, disease-free equilibrium E0, and the endemic equilibrium E*. O always exists; if Rd>1, E0 exists; if R0>1 and Rd>R1, E* exists.Theorem 3.
WhenRd>1 and R0<1, E0 is globally asymptotically stable.Proof.
By [33–35], we know that E0 is locally asymptotically stable. Now we define a Lyapunov function
(9)V=E+m+σσI≥0.
WhenRd>1 and R0<1, the Lyapunov function satisfies
(10)V˙=E˙+m+σσI˙=βSIN-(m+σ)(m+μ)σI≤[β-(m+σ)(m+μ)σ]I=(m+σ)(m+μ)σ[R0-1]I≤0.
Moreover,V˙=0 only hold when I=0. It is easy to verify that the disease-free equilibrium point E0 is the only fixed point of the system. Hence, applying the Lyapunov-LaSalle asymptotic stability theorem in [37, 38], the disease-free equilibrium point E0 is globally asymptotically stable.Since the proof of the stability of equilibriaO and E* is more difficult, we only give some numerical results.In sum, we can show the respective basins of attraction of the three equilibria which can be seen in Figure1 and confirmed in Figure 2. (1)
WhenRd<1, O is stable; see Figures 2(b) and 2(c).(2)
WhenRd>1 and R0<1, E0 is stable; see Figure 2(d).(3)
WhenRd>1, R0>1, and Rd<R1, O is stable; see Figure 2(a).(4)
WhenRd>1, R0>1, and Rd>R1, E* is stable; see Figure 2(e).Figure 1
R
d in terms of R0. The region of O1, O2, and O3 is the basin of attraction of equilibrium O; the region of E0 is the basin of attraction of equilibrium E0; the region of E* is the basin of attraction of equilibrium E*; the line L is Rd=R1=1+μ(m+μ)(m+σ)(R0-1)/mβ(m+σ+μ).The phase curves of the system under different initial conditions. (a) In the regionO3 with r=0.101 and β=0.5; (b) in the region O2 with r=0.05 and β=0.35; (c) in the region O1 with r=0.08 and β=0.2; (d) in the region E0 with r=0.13 and β=0.2; (e) in the region E* with r=0.13 and β=0.7. The value of other parameters can be seen in Table 1.
(a)
(b)
(c)
(d)
(e)


## 2.1. Model Formulation

The model is a system of SEI ordinary differential equations, whereS is the susceptible, E is the exposed, I is the infected, and N=S+E+I. This system considers the logistic growth rate and the standard incidence which is fit for the long-term growth of many large populations. The incubation period is considered for many diseases which do not develop symptoms immediately and need a period of time to accumulate a pathogen quantity for clinical outbreak, such as rabies, hand-foot-mouth disease, tuberculosis, and AIDS [22, 32]. The model we employ is as follows:
(1)dSdt=rN(1-Nk)-βSIN-mS,dEdt=βSIN-σE-mE,dIdt=σE-mI-μI,
where all parameters are positive whose interpretations can be seen in Table 1.Table 1
Descriptions and values of parameters in model (1).
Parameter
Interpretation
Value
r
The intrinsic growth rate
k
The carrying capacity
100000
β
The transmission rate
m
The natural mortality rate
0.1
σ
Clinical outcome rate
0.2
μ
The disease-induced mortality rate
0.1
a
The baseline contact rate
b
The magnitude of forcingNoticing the equations in model (1), we have
(2)dNdt=rN(1-Nk)-mN-μI.When there exists no disease, we have(3)dNdt=rN(1-Nk)-mN=[r(1-Nk)-m]N.LetRd=r/m, if Rd>1, N→N0=(1-m/r)k for N(0)>0, as t→+∞; that is, the population will grow and tend to a steady state N0. If Rd<1, then dN/dt<0 which will cause the population to disappear. Thus, Rd is the basic demographic reproduction number. From the above equation, the feasible region can be obtained: X={(S,E,I)∣S,E,I≥0,0≤S+E+I≤(1-m/r)k}, where r>m.Theorem 1.
The regionX is positively invariant with respect to system (1).


## 2.2. Dynamical Analysis

Let the right hand of system (1) to be zero; it is easy to see that system (1) has three equilibria:
(4)O=(0,0,0),E0=((1-mr)k,0,0),E*=(S*,E*,I*),
where O is the origin, E0 is the disease-free equilibrium, and E* is the endemic equilibrium. Concretely, one can have
(5)S*=[mβ(m+σ+μ)+μ(m+μ)(m+σ)(R0-1)]N*β[(m+μ)(m+σ)(R0-1)+m(m+σ+μ)],(6)E*=[r(1-N*/k)-m](m+μ)N*σμ,(7)I*=(m+μ)(m+σ)(R0-1)N*β(m+σ+μ),(8)N*=(k[βm(m+σ+μ)(Rd-1)+μ(m+σ)(m+μ)(1-R0)])×(rβ(m+σ+μ))-1,
where Rd=r/m is the basic demographic reproduction number and R0=βσ/(m+σ)(m+μ) is the basic reproduction number which can be obtained by the next-generation matrix method [33–35]. The introduction of the basic demographic reproduction number can be found in [36].Moreover, from (6) and (7), the conditions of the endemic equilibrium to exist are R0>1 and Rd>R1, where R1=1+μ(m+μ)(m+σ)(R0-1)/mβ(m+σ+μ). So, we can obtain the following theorems.Theorem 2.
The system (1) has three equilibria: origin O, disease-free equilibrium E0, and the endemic equilibrium E*. O always exists; if Rd>1, E0 exists; if R0>1 and Rd>R1, E* exists.Theorem 3.
WhenRd>1 and R0<1, E0 is globally asymptotically stable.Proof.
By [33–35], we know that E0 is locally asymptotically stable. Now we define a Lyapunov function
(9)V=E+m+σσI≥0.
WhenRd>1 and R0<1, the Lyapunov function satisfies
(10)V˙=E˙+m+σσI˙=βSIN-(m+σ)(m+μ)σI≤[β-(m+σ)(m+μ)σ]I=(m+σ)(m+μ)σ[R0-1]I≤0.
Moreover,V˙=0 only hold when I=0. It is easy to verify that the disease-free equilibrium point E0 is the only fixed point of the system. Hence, applying the Lyapunov-LaSalle asymptotic stability theorem in [37, 38], the disease-free equilibrium point E0 is globally asymptotically stable.Since the proof of the stability of equilibriaO and E* is more difficult, we only give some numerical results.In sum, we can show the respective basins of attraction of the three equilibria which can be seen in Figure1 and confirmed in Figure 2. (1)
WhenRd<1, O is stable; see Figures 2(b) and 2(c).(2)
WhenRd>1 and R0<1, E0 is stable; see Figure 2(d).(3)
WhenRd>1, R0>1, and Rd<R1, O is stable; see Figure 2(a).(4)
WhenRd>1, R0>1, and Rd>R1, E* is stable; see Figure 2(e).Figure 1
R
d in terms of R0. The region of O1, O2, and O3 is the basin of attraction of equilibrium O; the region of E0 is the basin of attraction of equilibrium E0; the region of E* is the basin of attraction of equilibrium E*; the line L is Rd=R1=1+μ(m+μ)(m+σ)(R0-1)/mβ(m+σ+μ).The phase curves of the system under different initial conditions. (a) In the regionO3 with r=0.101 and β=0.5; (b) in the region O2 with r=0.05 and β=0.35; (c) in the region O1 with r=0.08 and β=0.2; (d) in the region E0 with r=0.13 and β=0.2; (e) in the region E* with r=0.13 and β=0.7. The value of other parameters can be seen in Table 1.
(a)
(b)
(c)
(d)
(e)


## 3. Nonautonomous Model and Analysis

### 3.1. The Basic Reproduction Number

Now, we consider the non-autonomous case of the model (1) when the transmission rate is periodic, which is given as follows:
(11)dSdt=rN(1-Nk)-β(t)SIN-mS,dEdt=β(t)SIN-σE-mE,dIdt=σE-mI-μI,
where β(t) is a periodic function which is proposed by [39]. In the subsequent section, we will discuss the dynamical behavior of the system (11).For system (11), firstly we can give the basic reproduction number R0. According to the basic reproduction number under non-autonomous system, we can refer to the method of  [40, 41]. From the last section, we know that system (11) has one disease-free equilibrium E0=(N0,0,0), where N0=(1-m/r)k. By giving a new vector x=(E,I), we have
(12)F=(β(t)SIN0),V=(mE+σEmI+μI-σE),V-=(mE+σEmI+μI),V+=(0σE).Taking the partial derivative of the above vectors about variablesE, I and substituting the disease-free equilibrium, we have
(13)F(t)=(0β(t)00),V(t)=(m+σ0-σm+μ).According to [41], denote Cω to be the ordered Banach space of all ω-periodic functions from ℝ to ℝ4 which is equipped with the maximum norm ∥·∥ and the positive cone Cω+:={ϕ∈Cω:ϕ(t)≥0,∀t∈ℝ+}. Over the Banach space, we define a linear operator L:Cω→Cω by
(14)(Lϕ)(t)=∫0∞‍Y(t,t-a)F(t-a)ϕ(t-a)da,∀t∈ℝ+,ϕ∈Cω,
where L is called the next infection operator and the interpretation of Y(t,t-a), ϕ(t-a) can be seen in [41]. Then the spectral radius of L is defined as the basic reproduction number
(15)R0:=ρ(L).In order to give the expression of the basic reproduction number, we need to introduce the linearω-periodic system
(16)dwdt=[-V(t)+F(t)λ]w,t∈ℝ+,
with parameter λ∈ℝ. Let W(t,s,λ), t≥s, be the evolution operator of system (16) on ℝ2. In fact, W(t,s,λ)=Φ(F/λ)-V(t), and ΦF-V(t)=W(t,0,1), for all t≥0. By Theorems 2.1 and 2.2 in [41], the basic reproduction number also can be defined as λ0 such that ρ(Φ(F/λ0)-V(ω))=1, which can be straightforward to calculate.


### 3.2. Global Stability of the Disease-Free Equilibrium

Theorem 4.
The disease-free equilibriumE0 is globally asymptotically stable when R0<1 and Rd>1.Proof.
Theorem 2.2 in [41] implies that E0 is locally asymptotically stable when R0<1 and Rd>1. So we only need to prove its global attractability. It is easy to know that S(t)≤N0=(1-(m/r))k. Thus,
(17)dEdt≤β(t)I-(m+σ)E,dIdt=σE-mI-μI.
The right comparison system can be written as(18)dEdt=β(t)I-(m+σ)E,dIdt=σE-mI-μI;
that is,
(19)dhdt=(F(t)-V(t))h(t),h(t)=(E(t),I(t)).
For (19), Lemma 2.1 in [42] shows that there is a positive ω-periodic function h^(t)=(E(t),I(t))T such that h(t)=epth^(t) is a solution of system (18), where p=(1/ω)lnρ(ΦF-V(ω)). By Theorem 2.2 in  [41], we know that when R0<1 and Rd>1, ρ(ΦF-V(ω))<1 and p<0, which implies h(t)→0 as t→∞. Therefore, the zero solution of system (18) is globally asymptotically stable. By the comparison principle [43] and the theory of asymptotic autonomous systems [44], when R0<1 and Rd>1, E0 is globally attractive. Therefore, the proposition that E0 is globally asymptotically stable holds.


### 3.3. Existence of Positive Periodic Solutions

Before the proof of the existence of positive periodic solutions, we firstly introduce some denotations. Letu(t,x0) be the solution of system (11) with the initial value x0=(S(0),E(0),I(0)). By the fundamental existence-uniqueness theorem [45], u(t,x0) is the unique solution of system (11) with u(0,x0)=x0.Next, we need to introduce the Poincaré mapP:X→X associated with system (11); that is,
(20)P(x0)=u(ω,x0),∀x0∈X,
where ω is the period. Theorem 1 implies that X is positively invariant and P is a dissipative point.Now, we introduce two subsets ofX, X0:={(S,E,I)∈X:E>0,I>0} and ∂X0=X∖X0.Lemma 5.
(a) WhenR0>1 and r>m+μ, there exists a δ>0 such that when
(21)∥(S(0),E(0),I(0))-E0∥≤δ
for any (S(0),E(0),I(0))∈X0, one has
(22)limsupm→∞d[Pm(S(0),E(0),I(0)),E0]≥δ,
where E0=(N0,0,0).
(b) WhenR0>1 and r>m+μ, there exists a δ>0 such that when
(23)∥(S(0),E(0),I(0))-O∥≤δ
for any (S(0),E(0),I(0))∈X0, one has
(24)limsupm→∞d[Pm(S(0),E(0),I(0)),O]≥δ,
where O=(0,0,0).Proof.
(a) By Theorem 2.2 in [41], we know that when R0>1, ρ(ΦF-V(ω))>1. So there is a small enough positive number ϵ such that ρ(ΦF-V-Mϵ(ω))>1, where
(25)Mϵ=(0β(t)ϵIN000).
If proposition (a) does not hold, there is some(S(0),E(0),I(0))∈X0 such that
(26)limsupm→∞d(Pm(S(0),E(0),I(0)),E0)<δ.
We can assume that for allm≥0, d(Pm(S(0),E(0),I(0)),E0)<δ. Applying the continuity of the solutions with respect to the initial values,
(27)∥u(t,Pm(S(0),E(0),I(0)))-u(t,E0)∥≤ϵ,11111111111111111111∀m≥0,∀t1∈[0,ω].
Lett=mω+t1, where t1∈[0,ω] and m=[t/ω]. m=[t/ω] is the greatest integer which is not more than t/ω. Then, for any t≥0,
(28)∥u(t,(S(0),E(0),I(0)))-u(t,E0)∥=∥u(t1,Pm(S(0),E(0),I(0)))-u(t1,E0)∥≤ϵ.
SoN0-ϵ≤S(t)≤N0+ϵ. Then, when limsupm→∞d(Pm(S(0),E(0),I(0)),E0)<δ,
(29)dEdt≥β(t)I-β(t)ϵIN0-(m+σ)E,dIdt=σE-mI-μI.
Thus, we can study the right linear system(30)dEdt=β(t)I-β(t)ϵIN0-(m+σ)E,dIdt=σE-mI-μI.
For the system (30), there exists a positive ω-periodic function g^(t)=(E(t),I(t))T such that g(t)=eptg^(t) is a solution of system (11), where p=(1/ω)lnρ(ΦF-V-Mϵ(ω)). When R0>1, ρ(ΦF-V-Mϵ(ω))>1, which means that when g(0)>0, g(t)→∞ as t→∞. By the comparison principle [43], when E(0)>0, I(0)>0, E(t)→∞, I(t)→∞ as t→∞. There appears a contradiction. Thus, the proposition (a) holds.
(b) WhenR0>1 and r>m+μ, we have
(31)dNdt=rN(1-Nk)-mN-μI≥[r(1-Nk)-m-μ]N>0.
So ifRd>1, N→N0 for N(0)>0 as t→+∞; that is, WS(O)∩X0=∅.Theorem 6.
WhenR0>1 and r>m+μ, there exists a δ>0 such that any solution (S(t),E(t),I(t)) of system (11) with initial value (S(0),E(0),I(0))∈{(S,E,I)∈X:E>0,I>0} satisfies
(32)liminft→∞E(t)≥δ,liminft→∞I(t)≥δ,
and system (11) has at least one positive periodic solution.Proof.
From system (11),
(33)S(t)=e-∫0t‍(β(s)I(s)/N+m)ds×[S(0)+∫0t‍rN(s)(1-N(s)k)×e∫0s‍(β(c)I(c)/N(c)+m)dcds∫0t‍rN(s)(1-N(s)k)]>0,∀t>0,(34)E(t)=e-(m+σ)t[E(0)+∫0t‍β(s)S(s)I(s)N(s)e(m+σ)sds]>0,∀t>0,(35)I(t)=e-(m+μ)t[I(0)+∫0t‍σE(s)e(m+μ)sds]>0,∀t>0,
for any (S(0),E(0),I(0))∈X0, which shows that X0 is positively invariant. Moreover, it is obvious to see that ∂X0 is relatively closed in X. Denote
(36)M∂={(S(0),E(0),I(0))∈∂X0:Pm(S(0),E(0),I(0))∈∂X0,∀m≥0}.
Next, we prove that(37)M∂={(S,0,0)∈X:S≥0}.
We only need to show thatM∂⊆{(S,0,0)∈X:S≥0}, which means that for any (S(0),E(0),I(0))∈∂X0, E(mω)=I(mω)=0, for all m≥0. If it does not hold, there exists a m1≥0 such that (E(m1ω),I(m1ω))T>0.  Taking m1ω as the initial time and repeating the processes as in (33)–(35), we can have that (S(t),E(t),I(t))T>0, for all t>m1ω. Thus, (S(t),E(t),I(t))∈X0, for all t>m1ω. There appears a contradiction, which means that the equality (37) holds. Therefore, E0 is acyclic in ∂X0. Obviously, when R0>1 and r>m+μ, O is acyclic in ∂X0.
Furthermore, by Lemma5, E0=(N0,0,0) and O=(0,0,0) are isolated invariant sets in X, WS(E0)∩X0=∅, and WS(O)∩X0=∅. By Theorem 1.3.1 and Remark 1.3.1 in [46], it can be obtained that P is uniformly persistent with respect to (X0,∂X0); that is, there exists a δ>0 such that any solution (S(t),E(t),I(t)) of system (11) with the initial value (S(0),E(0),I(0))∈{(S,E,I)∈X:E>0,I>0} satisfies
(38)liminft→∞E(t)≥δ,liminft→∞I(t)≥δ.
Applying Theorem1.3.6 in [46], P has a fixed point
(39)(S*(0),E*(0),I*(0))∈X0.
From (33), we know S*>0,  for all t∈[0,ω]. S*(t) is also more than zero for all t>0 due to the periodicity. Similarly, for all t≥0, E*(t)>0, I*(t)>0. Therefore, it can be obtained that one of the positive ω-periodic solutions of system (11) is (S*(t),E*(t),I*(t)).


### 3.4. Numerical Simulations

Firstly, we give some notations. Ifg(t) is a periodic function with period ω, we define g-=(1/ω)∫0ω‍g(t)dt, gl=mint∈[0,ω]g(t), gu=maxt∈[0,ω]g(t). As described in the previous section,
(40)R1(t)=1+μ(m+μ)(m+σ)(R0-1)mβ(t)(m+σ+μ).
So R1l=1+μ(m+μ)(m+σ)(R0l-1)/mβu(m+σ+μ), R1u=1+μ(m+μ)(m+σ)(R0u-1)/mβl(m+σ+μ).In this section, we adoptβ(t)=a[1+bsin(πt/10)]. Then, applying the numerical simulation to verify the above solution, we give the following conclusion:(1)
whenRd<1, O is stable;(2)
whenRd>1 and R0<1, E0 is stable; see Figure 3(a);(3)
whenR0>1 and r>m+μ, system (11) has at least one positive periodic solution; see Figure 3(b).Phase plane ofS(t) and I(t). (a) When the parameter values are r=0.13, a=0.3, and b=0.2, R0=0.9029<1, Rd=1.3>1. (b) When the parameter values are r=0.3, a=0.7, and b=0.2, R0=2.1067 and r>m+μ. The value of other parameters can be seen in Table 1.
(a)
(b)We can give more results about the conditions of existence of the positive periodic solution.(1′)
WhenRd>1, R0>1, and Rd<R1u, O is stable; see Figures 4(a) and 4(b).(2′)
WhenRd>1, R0>1, and Rd>R1u, system (11) has at least one positive periodic solution, see Figure 5.Phase plane ofS(t) and I(t). (a) When the parameter values are r=0.11, a=0.9, b=0.2 and μ=0.3, R0=1.4991, Rd=1.1>1, Rd<R1u=1.4159, and Rd<R1l=1.2773. (b) When the parameter values are r=0.126, a=0.7, b=0.2, and μ=0.3, R0=2.1067, Rd=1.26>1, Rd<R1u=1.2964, and Rd>R1l=1.1976. The value of other parameters can be seen in Table 1.
(a)
(b)Figure 5
Phase plane ofS(t) and I(t). When the parameter values are r=0.13, a=0.36, and b=0.2, R0=1.0834, Rd=1.3>1, Rd>R1u=1.0435, and Rd>R1l=1.029. The value of other parameters can be seen in Table 1.By numerical simulations, we can give that the conditions which ensure the existence of positive periodic solution areR0>1 and r>m+μ or Rd>1, R0>1, and Rd>R1u. In fact, Rd>1, R0>1, and Rd>R1u are the sufficient conditions for R0>1 and r>m+μ. As a result, the conditions R0>1 and r>m+μ are broader.


## 3.1. The Basic Reproduction Number

Now, we consider the non-autonomous case of the model (1) when the transmission rate is periodic, which is given as follows:
(11)dSdt=rN(1-Nk)-β(t)SIN-mS,dEdt=β(t)SIN-σE-mE,dIdt=σE-mI-μI,
where β(t) is a periodic function which is proposed by [39]. In the subsequent section, we will discuss the dynamical behavior of the system (11).For system (11), firstly we can give the basic reproduction number R0. According to the basic reproduction number under non-autonomous system, we can refer to the method of  [40, 41]. From the last section, we know that system (11) has one disease-free equilibrium E0=(N0,0,0), where N0=(1-m/r)k. By giving a new vector x=(E,I), we have
(12)F=(β(t)SIN0),V=(mE+σEmI+μI-σE),V-=(mE+σEmI+μI),V+=(0σE).Taking the partial derivative of the above vectors about variablesE, I and substituting the disease-free equilibrium, we have
(13)F(t)=(0β(t)00),V(t)=(m+σ0-σm+μ).According to [41], denote Cω to be the ordered Banach space of all ω-periodic functions from ℝ to ℝ4 which is equipped with the maximum norm ∥·∥ and the positive cone Cω+:={ϕ∈Cω:ϕ(t)≥0,∀t∈ℝ+}. Over the Banach space, we define a linear operator L:Cω→Cω by
(14)(Lϕ)(t)=∫0∞‍Y(t,t-a)F(t-a)ϕ(t-a)da,∀t∈ℝ+,ϕ∈Cω,
where L is called the next infection operator and the interpretation of Y(t,t-a), ϕ(t-a) can be seen in [41]. Then the spectral radius of L is defined as the basic reproduction number
(15)R0:=ρ(L).In order to give the expression of the basic reproduction number, we need to introduce the linearω-periodic system
(16)dwdt=[-V(t)+F(t)λ]w,t∈ℝ+,
with parameter λ∈ℝ. Let W(t,s,λ), t≥s, be the evolution operator of system (16) on ℝ2. In fact, W(t,s,λ)=Φ(F/λ)-V(t), and ΦF-V(t)=W(t,0,1), for all t≥0. By Theorems 2.1 and 2.2 in [41], the basic reproduction number also can be defined as λ0 such that ρ(Φ(F/λ0)-V(ω))=1, which can be straightforward to calculate.


## 3.2. Global Stability of the Disease-Free Equilibrium

Theorem 4.
The disease-free equilibriumE0 is globally asymptotically stable when R0<1 and Rd>1.Proof.
Theorem 2.2 in [41] implies that E0 is locally asymptotically stable when R0<1 and Rd>1. So we only need to prove its global attractability. It is easy to know that S(t)≤N0=(1-(m/r))k. Thus,
(17)dEdt≤β(t)I-(m+σ)E,dIdt=σE-mI-μI.
The right comparison system can be written as(18)dEdt=β(t)I-(m+σ)E,dIdt=σE-mI-μI;
that is,
(19)dhdt=(F(t)-V(t))h(t),h(t)=(E(t),I(t)).
For (19), Lemma 2.1 in [42] shows that there is a positive ω-periodic function h^(t)=(E(t),I(t))T such that h(t)=epth^(t) is a solution of system (18), where p=(1/ω)lnρ(ΦF-V(ω)). By Theorem 2.2 in  [41], we know that when R0<1 and Rd>1, ρ(ΦF-V(ω))<1 and p<0, which implies h(t)→0 as t→∞. Therefore, the zero solution of system (18) is globally asymptotically stable. By the comparison principle [43] and the theory of asymptotic autonomous systems [44], when R0<1 and Rd>1, E0 is globally attractive. Therefore, the proposition that E0 is globally asymptotically stable holds.


## 3.3. Existence of Positive Periodic Solutions

Before the proof of the existence of positive periodic solutions, we firstly introduce some denotations. Letu(t,x0) be the solution of system (11) with the initial value x0=(S(0),E(0),I(0)). By the fundamental existence-uniqueness theorem [45], u(t,x0) is the unique solution of system (11) with u(0,x0)=x0.Next, we need to introduce the Poincaré mapP:X→X associated with system (11); that is,
(20)P(x0)=u(ω,x0),∀x0∈X,
where ω is the period. Theorem 1 implies that X is positively invariant and P is a dissipative point.Now, we introduce two subsets ofX, X0:={(S,E,I)∈X:E>0,I>0} and ∂X0=X∖X0.Lemma 5.
(a) WhenR0>1 and r>m+μ, there exists a δ>0 such that when
(21)∥(S(0),E(0),I(0))-E0∥≤δ
for any (S(0),E(0),I(0))∈X0, one has
(22)limsupm→∞d[Pm(S(0),E(0),I(0)),E0]≥δ,
where E0=(N0,0,0).
(b) WhenR0>1 and r>m+μ, there exists a δ>0 such that when
(23)∥(S(0),E(0),I(0))-O∥≤δ
for any (S(0),E(0),I(0))∈X0, one has
(24)limsupm→∞d[Pm(S(0),E(0),I(0)),O]≥δ,
where O=(0,0,0).Proof.
(a) By Theorem 2.2 in [41], we know that when R0>1, ρ(ΦF-V(ω))>1. So there is a small enough positive number ϵ such that ρ(ΦF-V-Mϵ(ω))>1, where
(25)Mϵ=(0β(t)ϵIN000).
If proposition (a) does not hold, there is some(S(0),E(0),I(0))∈X0 such that
(26)limsupm→∞d(Pm(S(0),E(0),I(0)),E0)<δ.
We can assume that for allm≥0, d(Pm(S(0),E(0),I(0)),E0)<δ. Applying the continuity of the solutions with respect to the initial values,
(27)∥u(t,Pm(S(0),E(0),I(0)))-u(t,E0)∥≤ϵ,11111111111111111111∀m≥0,∀t1∈[0,ω].
Lett=mω+t1, where t1∈[0,ω] and m=[t/ω]. m=[t/ω] is the greatest integer which is not more than t/ω. Then, for any t≥0,
(28)∥u(t,(S(0),E(0),I(0)))-u(t,E0)∥=∥u(t1,Pm(S(0),E(0),I(0)))-u(t1,E0)∥≤ϵ.
SoN0-ϵ≤S(t)≤N0+ϵ. Then, when limsupm→∞d(Pm(S(0),E(0),I(0)),E0)<δ,
(29)dEdt≥β(t)I-β(t)ϵIN0-(m+σ)E,dIdt=σE-mI-μI.
Thus, we can study the right linear system(30)dEdt=β(t)I-β(t)ϵIN0-(m+σ)E,dIdt=σE-mI-μI.
For the system (30), there exists a positive ω-periodic function g^(t)=(E(t),I(t))T such that g(t)=eptg^(t) is a solution of system (11), where p=(1/ω)lnρ(ΦF-V-Mϵ(ω)). When R0>1, ρ(ΦF-V-Mϵ(ω))>1, which means that when g(0)>0, g(t)→∞ as t→∞. By the comparison principle [43], when E(0)>0, I(0)>0, E(t)→∞, I(t)→∞ as t→∞. There appears a contradiction. Thus, the proposition (a) holds.
(b) WhenR0>1 and r>m+μ, we have
(31)dNdt=rN(1-Nk)-mN-μI≥[r(1-Nk)-m-μ]N>0.
So ifRd>1, N→N0 for N(0)>0 as t→+∞; that is, WS(O)∩X0=∅.Theorem 6.
WhenR0>1 and r>m+μ, there exists a δ>0 such that any solution (S(t),E(t),I(t)) of system (11) with initial value (S(0),E(0),I(0))∈{(S,E,I)∈X:E>0,I>0} satisfies
(32)liminft→∞E(t)≥δ,liminft→∞I(t)≥δ,
and system (11) has at least one positive periodic solution.Proof.
From system (11),
(33)S(t)=e-∫0t‍(β(s)I(s)/N+m)ds×[S(0)+∫0t‍rN(s)(1-N(s)k)×e∫0s‍(β(c)I(c)/N(c)+m)dcds∫0t‍rN(s)(1-N(s)k)]>0,∀t>0,(34)E(t)=e-(m+σ)t[E(0)+∫0t‍β(s)S(s)I(s)N(s)e(m+σ)sds]>0,∀t>0,(35)I(t)=e-(m+μ)t[I(0)+∫0t‍σE(s)e(m+μ)sds]>0,∀t>0,
for any (S(0),E(0),I(0))∈X0, which shows that X0 is positively invariant. Moreover, it is obvious to see that ∂X0 is relatively closed in X. Denote
(36)M∂={(S(0),E(0),I(0))∈∂X0:Pm(S(0),E(0),I(0))∈∂X0,∀m≥0}.
Next, we prove that(37)M∂={(S,0,0)∈X:S≥0}.
We only need to show thatM∂⊆{(S,0,0)∈X:S≥0}, which means that for any (S(0),E(0),I(0))∈∂X0, E(mω)=I(mω)=0, for all m≥0. If it does not hold, there exists a m1≥0 such that (E(m1ω),I(m1ω))T>0.  Taking m1ω as the initial time and repeating the processes as in (33)–(35), we can have that (S(t),E(t),I(t))T>0, for all t>m1ω. Thus, (S(t),E(t),I(t))∈X0, for all t>m1ω. There appears a contradiction, which means that the equality (37) holds. Therefore, E0 is acyclic in ∂X0. Obviously, when R0>1 and r>m+μ, O is acyclic in ∂X0.
Furthermore, by Lemma5, E0=(N0,0,0) and O=(0,0,0) are isolated invariant sets in X, WS(E0)∩X0=∅, and WS(O)∩X0=∅. By Theorem 1.3.1 and Remark 1.3.1 in [46], it can be obtained that P is uniformly persistent with respect to (X0,∂X0); that is, there exists a δ>0 such that any solution (S(t),E(t),I(t)) of system (11) with the initial value (S(0),E(0),I(0))∈{(S,E,I)∈X:E>0,I>0} satisfies
(38)liminft→∞E(t)≥δ,liminft→∞I(t)≥δ.
Applying Theorem1.3.6 in [46], P has a fixed point
(39)(S*(0),E*(0),I*(0))∈X0.
From (33), we know S*>0,  for all t∈[0,ω]. S*(t) is also more than zero for all t>0 due to the periodicity. Similarly, for all t≥0, E*(t)>0, I*(t)>0. Therefore, it can be obtained that one of the positive ω-periodic solutions of system (11) is (S*(t),E*(t),I*(t)).


## 3.4. Numerical Simulations

Firstly, we give some notations. Ifg(t) is a periodic function with period ω, we define g-=(1/ω)∫0ω‍g(t)dt, gl=mint∈[0,ω]g(t), gu=maxt∈[0,ω]g(t). As described in the previous section,
(40)R1(t)=1+μ(m+μ)(m+σ)(R0-1)mβ(t)(m+σ+μ).
So R1l=1+μ(m+μ)(m+σ)(R0l-1)/mβu(m+σ+μ), R1u=1+μ(m+μ)(m+σ)(R0u-1)/mβl(m+σ+μ).In this section, we adoptβ(t)=a[1+bsin(πt/10)]. Then, applying the numerical simulation to verify the above solution, we give the following conclusion:(1)
whenRd<1, O is stable;(2)
whenRd>1 and R0<1, E0 is stable; see Figure 3(a);(3)
whenR0>1 and r>m+μ, system (11) has at least one positive periodic solution; see Figure 3(b).Phase plane ofS(t) and I(t). (a) When the parameter values are r=0.13, a=0.3, and b=0.2, R0=0.9029<1, Rd=1.3>1. (b) When the parameter values are r=0.3, a=0.7, and b=0.2, R0=2.1067 and r>m+μ. The value of other parameters can be seen in Table 1.
(a)
(b)We can give more results about the conditions of existence of the positive periodic solution.(1′)
WhenRd>1, R0>1, and Rd<R1u, O is stable; see Figures 4(a) and 4(b).(2′)
WhenRd>1, R0>1, and Rd>R1u, system (11) has at least one positive periodic solution, see Figure 5.Phase plane ofS(t) and I(t). (a) When the parameter values are r=0.11, a=0.9, b=0.2 and μ=0.3, R0=1.4991, Rd=1.1>1, Rd<R1u=1.4159, and Rd<R1l=1.2773. (b) When the parameter values are r=0.126, a=0.7, b=0.2, and μ=0.3, R0=2.1067, Rd=1.26>1, Rd<R1u=1.2964, and Rd>R1l=1.1976. The value of other parameters can be seen in Table 1.
(a)
(b)Figure 5
Phase plane ofS(t) and I(t). When the parameter values are r=0.13, a=0.36, and b=0.2, R0=1.0834, Rd=1.3>1, Rd>R1u=1.0435, and Rd>R1l=1.029. The value of other parameters can be seen in Table 1.By numerical simulations, we can give that the conditions which ensure the existence of positive periodic solution areR0>1 and r>m+μ or Rd>1, R0>1, and Rd>R1u. In fact, Rd>1, R0>1, and Rd>R1u are the sufficient conditions for R0>1 and r>m+μ. As a result, the conditions R0>1 and r>m+μ are broader.


## 4. Discussion

This paper considers a logistic growth system whose birth process incorporates density-dependent effects. This type of model has a rich dynamical behavior and practical significance. By analyzing its equilibria and respective attractive region, we find that the dynamical behavior of a disease will be determined by two thresholdsR0 and Rd. Only R0>1 cannot promise the existence of the endemic equilibrium which also needs Rd>R1. When R0>1 and Rd<R1, the solutions of the system (1) will tend to the origin O. It is caused by the phenomenon that the death number due to disease cannot be supplemented by the birth number promptly. Finally, all people are infected and die out. The fact interpreted by this model is more reasonable. Theoretically, we prove the global asymptotic stability of the disease-free equilibrium and give respective attractive regions of equilibria.Seasonally effective contact rate is the most common form which may be related to various factors, and thus this paper studies the corresponding non-autonomous system which is obtained by changing the constant transmission rate of the above system into the periodic transmission rate. For the periodic systems, their dynamical behaviors, especially the basic reproduction number, have been investigated in depth by [41, 47–55] which provide many methods that we can utilize. For the obtained periodic model, by analyzing the global asymptotic stability of the disease-free equilibrium and the existence of positive periodic solution, we have the similar results as the autonomous system. The dynamic behavior of disease will be decided by two conditions R0>1 and r>m+μ that show that when the disease is prevalent, the birth rate should be larger than the death rate to guarantee the sustainable growth of population. Otherwise, the population will disappear. In addition, we will evaluate and compare the basic reproduction number R0 and the average basic reproduction number R-0 which has been adopted by [27, 56–59]. In this paper, we can calculate the average reproduction number
(41)R-0=β-σ(m+σ)(m+μ),
where β-=(1/20)∫020‍β(t)dt. When r=0.13, a=0.3, and b=0.2, we know that R0=0.9029 and R-0=1. When r=0.13, k=100000, a=0.36, b=0.2, and m=0.1, σ=0.2, μ=0.1, then R0=1.0834 and R-0=1.2. In that sense, it is confirmed that the basic reproduction number R0 defined by [40] is more accurate than the average reproduction number R-0 which overestimates the risk of disease.It should be noted that we live in a spatial world and it is a natural phenomenon that a substance goes from high density regions to low density regions. As a result, epidemic models should include spatial effects. In a further study, we need to investigate spatial epidemic models with seasonal factors.



---

*Source: 470646-2013-11-10.xml*