Hindawi Publishing Corporation The Scientific World Journal Volume 2013, Article ID 470646, 10 pages http://dx.doi.org/10.1155/2013/470646
Research Article Positive Periodic Solutions of an Epidemic Model with Seasonality Gui-Quan Sun,1,2,3 Zhenguo Bai,4 Zi-Ke Zhang,5,6,7 Tao Zhou,6 and Zhen Jin1,2 1
Complex Sciences Center, Shanxi University, Taiyuan 030006, China School of Mathematical Sciences, Shanxi University, Taiyuan, Shanxi 030006, China 3 Department of Mathematics, North University of China, Taiyuan, Shanxi 030051, China 4 Department of Applied Mathematics, Xidian University, Xiβan, Shaanxi 710071, China 5 Institute of Information Economy, Hangzhou Normal University, Hangzhou 310036, China 6 Web Sciences Center, University of Electronic Science and Technology of China, Chengdu, Sichuan 610054, China 7 Department of Physics, University of Fribourg, Chemin du MusΒ΄ee 3, 1700 Fribourg, Switzerland 2
Correspondence should be addressed to Gui-Quan Sun;
[email protected] Received 17 August 2013; Accepted 12 September 2013 Academic Editors: P. Leach, O. D. Makinde, and Y. Xia Copyright Β© 2013 Gui-Quan Sun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 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 number π
π is obtained. Moreover, only the basic reproduction number π
0 cannot ensure the existence of the positive equilibrium, which needs additional condition π
π > π
1 . 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.
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
2 growth. Muroya [13] investigated discrete models of nonautonomous 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 threshold π
π which has a similar meaning and is called as the basic demographic reproduction number. If π
π > 1, the population grows; that is, a critical mass of individuals for the disease to spread may be supported. If π
π < 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 π
0 and π
π . 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 π(π‘) = π(1 + π sin(ππ‘/π)) where π is the amplitude of seasonal variation in transmission (typically referred to as the βstrength of seasonal forcingβ) and 2π 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 π(π‘) = π(1 + π term(π‘)), 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 π(π‘) = π(1 + π)term(π‘) [29]. Here, we take the form π½(π‘) = π[1 + π sin(ππ‘/10)]. The paper is organized as follows. In Section 2, 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, where π is the susceptible, πΈ is the exposed, πΌ is the infected, and π = π + πΈ + πΌ. This system considers the logistic growth rate and the standard
The Scientific World Journal Table 1: Descriptions and values of parameters in model (1). Parameter π π π½ π π π π π
Interpretation The intrinsic growth rate The carrying capacity The transmission rate The natural mortality rate Clinical outcome rate The disease-induced mortality rate The baseline contact rate The magnitude of forcing
Value 100000 0.1 0.2 0.1
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 β ) β β ππ, ππ‘ π π ππΈ π½ππΌ = β ππΈ β ππΈ, ππ‘ π
(1)
ππΌ = ππΈ β ππΌ β ππΌ, ππ‘ where all parameters are positive whose interpretations can be seen in Table 1. Noticing the equations in model (1), we have π ππ = ππ (1 β ) β ππ β ππΌ. ππ‘ π
(2)
When there exists no disease, we have π π ππ = ππ (1 β ) β ππ = [π (1 β ) β π] π. ππ‘ π π
(3)
Let π
π = π/π, if π
π > 1, π β π0 = (1 β π/π)π for π(0) > 0, as π‘ β +β; that is, the population will grow and tend to a steady state π0 . If π
π < 1, then ππ/ππ‘ < 0 which will cause the population to disappear. Thus, π
π is the basic demographic reproduction number. From the above equation, the feasible region can be obtained: π = {(π, πΈ, πΌ) | π, πΈ, πΌ β₯ 0, 0 β€ π + πΈ + πΌ β€ (1 β π/π)π}, where π > π. Theorem 1. The region π 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: π = (0, 0, 0) , πΈ0 = ((1 β
π ) π, 0, 0) , π
πΈβ = (πβ , πΈβ , πΌβ ) ,
(4)
The Scientific World Journal
3
where π is the origin, πΈ0 is the disease-free equilibrium, and πΈβ is the endemic equilibrium. Concretely, one can have [ππ½ (π + π + π) + π (π + π) (π + π) (π
0 β 1)] πβ , π½ [(π + π) (π + π) (π
0 β 1) + π (π + π + π)] (5) πΈβ =
[π (1 β πβ /π) β π] (π + π) πβ , ππ
(6)
(π + π) (π + π) (π
0 β 1) πβ , π½ (π + π + π)
(7)
πΌβ =
πβ = (π [π½π (π + π + π) (π
π β 1) +π (π + π) (π + π) (1 β π
0 )])
(8)
β1
Γ (ππ½ (π + π + π)) , where π
π = π/π is the basic demographic reproduction number and π
0 = π½π/(π+π)(π+π) 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 π
0 > 1 and π
π > π
1 , where π
1 = 1 + π(π + π)(π + π)(π
0 β 1)/ππ½(π + π + π). So, we can obtain the following theorems. Theorem 2. The system (1) has three equilibria: origin π, disease-free equilibrium πΈ0 , and the endemic equilibrium πΈβ . π always exists; if π
π > 1, πΈ0 exists; if π
0 > 1 and π
π > π
1 , πΈβ exists. Theorem 3. When π
π > 1 and π
0 < 1, πΈ0 is globally asymptotically stable.
O3 Rd
πβ =
L
Eβ
E0 1
O1
0
O2
1 R0
Figure 1: π
π in terms of π
0 . The region of π1 , π2 , and π3 is the basin of attraction of equilibrium π; the region of πΈ0 is the basin of attraction of equilibrium πΈ0 ; the region of πΈβ is the basin of attraction of equilibrium πΈβ ; the line πΏ is π
π = π
1 = 1 + π(π + π)(π + π)(π
0 β 1)/ππ½(π + π + π).
asymptotic stability theorem in [37, 38], the disease-free equilibrium point πΈ0 is globally asymptotically stable. Since the proof of the stability of equilibria π and πΈβ 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 Figure 1 and confirmed in Figure 2. (1) When π
π < 1, π is stable; see Figures 2(b) and 2(c). (2) When π
π > 1 and π
0 < 1, πΈ0 is stable; see Figure 2(d). (3) When π
π > 1, π
0 > 1, and π
π < π
1 , π is stable; see Figure 2(a). (4) When π
π > 1, π
0 > 1, and π
π > π
1 , πΈβ is stable; see Figure 2(e).
Proof. By [33β35], we know that πΈ0 is locally asymptotically stable. Now we define a Lyapunov function π=πΈ+
π+π πΌ β₯ 0. π
3. Nonautonomous Model and Analysis (9)
When π
π > 1 and π
0 < 1, the Lyapunov function satisfies π+π Μ πΌ πΜ = πΈΜ + π =
=
π½ (π‘) ππΌ π ππ = ππ (1 β ) β β ππ, ππ‘ π π
π½ππΌ (π + π) (π + π) β πΌ π π
β€ [π½ β
(π + π) (π + π) ]πΌ π
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:
(10)
(π + π) (π + π) [π
0 β 1] πΌ π
β€ 0. Moreover, πΜ = 0 only hold when πΌ = 0. It is easy to verify that the disease-free equilibrium point πΈ0 is the only fixed point of the system. Hence, applying the Lyapunov-LaSalle
ππΈ π½ (π‘) ππΌ = β ππΈ β ππΈ, ππ‘ π
(11)
ππΌ = ππΈ β ππΌ β ππΌ, ππ‘ where π½(π‘) 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 π
0 . 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)
The Scientific World Journal 80
80
70
70
60
60
50
50 I(t)
I(t)
4
40
40
30
30
20
20
10
10
0
0
50
100
R0 = 1.6667 Rd = 1.01
150 S(t)
200
250
0
300
R1 = 1.2
0
50
100
R0 = 1.667 Rd = 0.5
150 S(t)
200
250
300
R1 = 1.0714
(a)
(b)
80
600
70
500
60 400 I(t)
I(t)
50 40 30
300 200
20 100
10 0
0
50
100
150
200
250
0
300
0
0.5
1
1.5
2 S(t)
S(t) R0 = 0.6667 Rd = 0.8
R0 = 0.6667 Rd = 1.3
R1 = 0.75
2.5
3
3.5
4 Γ104
R1 = 0.75
(c)
(d)
600 500
I(t)
400 300 200 100 0 100
200
R0 = 2.333 Rd = 1.3
300
400
500 S(t)
600
700
800
900
R1 = 1.2857 (e)
Figure 2: The phase curves of the system under different initial conditions. (a) In the region π3 with π = 0.101 and π½ = 0.5; (b) in the region π2 with π = 0.05 and π½ = 0.35; (c) in the region π1 with π = 0.08 and π½ = 0.2; (d) in the region πΈ0 with π = 0.13 and π½ = 0.2; (e) in the region πΈβ with π = 0.13 and π½ = 0.7. The value of other parameters can be seen in Table 1.
The Scientific World Journal
5
has one disease-free equilibrium πΈ0 = (π0 , 0, 0), where π0 = (1 β π/π)π. By giving a new vector π₯ = (πΈ, πΌ), we have π½ (π‘) ππΌ πΉ = ( π ), 0
ππΈ + ππΈ π=( ), ππΌ + ππΌ β ππΈ
ππΈ + ππΈ ), πβ = ( ππΌ + ππΌ
to prove its global attractability. It is easy to know that π(π‘) β€ π0 = (1 β (π/π))π. Thus, ππΈ β€ π½ (π‘) πΌ β (π + π) πΈ, ππ‘ ππΌ = ππΈ β ππΌ β ππΌ. ππ‘
(12)
0 ). ππΈ
π+ = (
The right comparison system can be written as ππΈ = π½ (π‘) πΌ β (π + π) πΈ, ππ‘
Taking the partial derivative of the above vectors about variables πΈ, πΌ and substituting the disease-free equilibrium, we have 0 π½ (π‘) πΉ (π‘) = ( ), 0 0 π+π 0 π (π‘) = ( ). βπ π + π
(13)
According to [41], denote πΆπ to be the ordered Banach space of all π-periodic functions from R to R4 which is equipped with the maximum norm β β
β and the positive cone πΆπ+ := {π β πΆπ : π(π‘) β₯ 0, βπ‘ β R+ }. Over the Banach space, we define a linear operator πΏ : πΆπ β πΆπ by (πΏπ) (π‘) β
= β« π (π‘, π‘ β π) πΉ (π‘ β π) π (π‘ β π) ππ, 0
βπ‘ β R+ , π β πΆπ , (14)
where πΏ is called the next infection operator and the interpretation of π(π‘, π‘ β π), π(π‘ β π) can be seen in [41]. Then the spectral radius of πΏ is defined as the basic reproduction number π
0 := π (πΏ) .
(15)
In order to give the expression of the basic reproduction number, we need to introduce the linear π-periodic system ππ€ πΉ (π‘) = [βπ (π‘) + ] π€, ππ‘ π
π‘ β R+ ,
(17)
(16)
with parameter π β R. Let π(π‘, π , π), π‘ β₯ π , be the evolution operator of system (16) on R2 . In fact, π(π‘, π , π) = Ξ¦(πΉ/π)βπ (π‘), and Ξ¦πΉβπ(π‘) = π(π‘, 0, 1), for all π‘ β₯ 0. By Theorems 2.1 and 2.2 in [41], the basic reproduction number also can be defined as π 0 such that π(Ξ¦(πΉ/π 0 )βπ (π)) = 1, which can be straightforward to calculate. 3.2. Global Stability of the Disease-Free Equilibrium Theorem 4. The disease-free equilibrium πΈ0 is globally asymptotically stable when π
0 < 1 and π
π > 1. Proof. Theorem 2.2 in [41] implies that πΈ0 is locally asymptotically stable when π
0 < 1 and π
π > 1. So we only need
ππΌ = ππΈ β ππΌ β ππΌ; ππ‘
(18)
that is, πβ = (πΉ (π‘) β π (π‘)) β (π‘) , ππ‘
β (π‘) = (πΈ (π‘) , πΌ (π‘)) .
(19)
For (19), Lemma 2.1 in [42] shows that there is a positive π-periodic function Μβ(π‘) = (πΈ(π‘), πΌ(π‘))π such that β(π‘) = πππ‘ Μβ(π‘) is a solution of system (18), where π = (1/π) ln π(Ξ¦πΉβπ(π)). By Theorem 2.2 in [41], we know that when π
0 < 1 and π
π > 1, π(Ξ¦πΉβπ(π)) < 1 and π < 0, which implies β(π‘) β 0 as π‘ β β. 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 π
0 < 1 and π
π > 1, πΈ0 is globally attractive. Therefore, the proposition that πΈ0 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. Let π’(π‘, π₯0 ) be the solution of system (11) with the initial value π₯0 = (π(0), πΈ(0), πΌ(0)). By the fundamental existence-uniqueness theorem [45], π’(π‘, π₯0 ) is the unique solution of system (11) with π’(0, π₯0 ) = π₯0 . Next, we need to introduce the PoincarΒ΄e map π : π β π associated with system (11); that is, π (π₯0 ) = π’ (π, π₯0 ) ,
βπ₯0 β π,
(20)
where π is the period. Theorem 1 implies that π is positively invariant and π is a dissipative point. Now, we introduce two subsets of π, π0 := {(π, πΈ, πΌ) β π : πΈ > 0, πΌ > 0} and ππ0 = π \ π0 . Lemma 5. (a) When π
0 > 1 and π > π + π, there exists a πΏ > 0 such that when σ΅©σ΅© σ΅© σ΅©σ΅©(π (0) , πΈ (0) , πΌ (0)) β πΈ0 σ΅©σ΅©σ΅© β€ πΏ
(21)
for any (π(0), πΈ(0), πΌ(0)) β π0 , one has lim sup π [ππ (π (0) , πΈ (0) , πΌ (0)) , πΈ0 ] β₯ πΏ, πββ
where πΈ0 = (π0 , 0, 0).
(22)
6
The Scientific World Journal
(b) When π
0 > 1 and π > π + π, there exists a πΏ > 0 such that when β(π (0) , πΈ (0) , πΌ (0)) β πβ β€ πΏ
(23)
for any (π(0), πΈ(0), πΌ(0)) β π0 , one has lim sup π [ππ (π (0) , πΈ (0) , πΌ (0)) , π] β₯ πΏ, πββ
π(0) > 0, π(π‘) β β as π‘ β β. By the comparison principle [43], when πΈ(0) > 0, πΌ(0) > 0, πΈ(π‘) β β, πΌ(π‘) β β as π‘ β β. There appears a contradiction. Thus, the proposition (a) holds. (b) When π
0 > 1 and π > π + π, we have ππ π = ππ (1 β ) β ππ β ππΌ ππ‘ π
(24)
π β₯ [π (1 β ) β π β π] π > 0. π
where π = (0, 0, 0).
(31)
Proof. (a) By Theorem 2.2 in [41], we know that when π
0 > 1, π(Ξ¦πΉβπ(π)) > 1. So there is a small enough positive number π such that π(Ξ¦πΉβπβππ (π)) > 1, where
So if π
π > 1, π β π0 for π(0) > 0 as π‘ β +β; that is, π (π) β© π0 = 0.
π½ (π‘) ππΌ π0 ) . ππ = ( 0 0
Theorem 6. When π
0 > 1 and π > π + π, there exists a πΏ > 0 such that any solution (π(π‘), πΈ(π‘), πΌ(π‘)) of system (11) with initial value (π(0), πΈ(0), πΌ(0)) β {(π, πΈ, πΌ) β π : πΈ > 0, πΌ > 0} satisfies
0
(25)
If proposition (a) does not hold, there is some (π(0), πΈ(0), πΌ(0)) β π0 such that lim sup π (ππ (π (0) , πΈ (0) , πΌ (0)) , πΈ0 ) < πΏ. πββ
π
lim inf πΈ (π‘) β₯ πΏ,
We can assume that for all π β₯ 0, π(π (π(0), πΈ(0), πΌ(0)), πΈ0 ) < πΏ. Applying the continuity of the solutions with respect to the initial values,
βπ β₯ 0, βπ‘1 β [0, π] .
0
Γπ > 0,
π (π ) ) π
π
(29)
π‘
0
> 0,
(33) ππ ]
βπ‘ > 0,
πΈ (π‘) = πβ(π+π)π‘ [πΈ (0) + β«
π½ (π ) π (π ) πΌ (π ) (π+π)π π ππ ] π (π ) (34)
βπ‘ > 0, π‘
πΌ (π‘) = πβ(π+π)π‘ [πΌ (0) + β« ππΈ (π ) π(π+π)π ππ ] 0
> 0,
for any (π(0), πΈ(0), πΌ(0)) β π0 , which shows that π0 is positively invariant. Moreover, it is obvious to see that ππ0 is relatively closed in π. Denote
ππ (π (0) , πΈ (0) , πΌ (0)) β ππ0 , βπ β₯ 0} . (30)
For the system (30), there exists a positive π-periodic Μ Μ is a function π(π‘) = (πΈ(π‘), πΌ(π‘))π such that π(π‘) = πππ‘ π(π‘) solution of system (11), where π = (1/π) ln π(Ξ¦πΉβπβππ (π)). When π
0 > 1, π(Ξ¦πΉβπβππ (π)) > 1, which means that when
(35)
βπ‘ > 0,
ππ = {(π (0) , πΈ (0) , πΌ (0)) β ππ0 :
Thus, we can study the right linear system
ππΌ = ππΈ β ππΌ β ππΌ. ππ‘
π‘
β«0 (π½(π)πΌ(π)/π(π)+π)ππ
So π0 β π β€ π(π‘) β€ π0 + π. Then, when lim supπ β β π(ππ (π(0), πΈ(0), πΌ(0)), πΈ0 ) < πΏ,
π½ (π‘) ππΌ ππΈ β (π + π) πΈ, = π½ (π‘) πΌ β ππ‘ π0
π‘
π (π‘) = πβ β«0 (π½(π )πΌ(π )/π+π)ππ
(27)
σ΅© σ΅©σ΅© σ΅©σ΅©π’ (π‘, (π (0) , πΈ (0) , πΌ (0))) β π’ (π‘, πΈ0 )σ΅©σ΅©σ΅© (28) σ΅© σ΅© = σ΅©σ΅©σ΅©π’ (π‘1 , ππ (π (0) , πΈ (0) , πΌ (0))) β π’ (π‘1 , πΈ0 )σ΅©σ΅©σ΅© β€ π.
ππΌ = ππΈ β ππΌ β ππΌ. ππ‘
(32)
Proof. From system (11),
Γ [π (0) + β« ππ (π ) (1 β
Let π‘ = ππ + π‘1 , where π‘1 β [0, π] and π = [π‘/π]. π = [π‘/π] is the greatest integer which is not more than π‘/π. Then, for any π‘ β₯ 0,
π½ (π‘) ππΌ ππΈ β (π + π) πΈ, β₯ π½ (π‘) πΌ β ππ‘ π0
π‘ββ
and system (11) has at least one positive periodic solution. (26)
π
σ΅©σ΅© σ΅© π σ΅©σ΅©π’ (π‘, π (π (0) , πΈ (0) , πΌ (0))) β π’ (π‘, πΈ0 )σ΅©σ΅©σ΅© β€ π,
lim inf πΌ (π‘) β₯ πΏ,
π‘ββ
(36)
Next, we prove that ππ = {(π, 0, 0) β π : π β₯ 0} .
(37)
We only need to show that ππ β {(π, 0, 0) β π : π β₯ 0}, which means that for any (π(0), πΈ(0), πΌ(0)) β ππ0 , πΈ(ππ) = πΌ(ππ) = 0, for all π β₯ 0. If it does not hold, there exists a π1 β₯ 0 such that (πΈ(π1 π), πΌ(π1 π))π > 0. Taking π1 π as
The Scientific World Journal
7 Γ104 2.4
Γ104 2.5
2.3 2.2
2
2.1 2
I(t)
I(t)
1.5
1.9 1.8
1
1.7 1.6
0.5
1.5 0
1
1.5
2
2.5
3
S(t)
3.5
1.4 1.5
2
2.5
Γ104
3
S(t)
3.5 Γ104
S(0) = 25000, E(0) = 10000, I(0) = 24000 S(0) = 35000, E(0) = 20000, I(0) = 15000 S(0) = 15000, E(0) = 10000, I(0) = 20000
S(0) = 25000, E(0) = 10000, I(0) = 24000 S(0) = 35000, E(0) = 20000, I(0) = 15000 S(0) = 15000, E(0) = 10000, I(0) = 20000 (a)
(b)
Figure 3: Phase plane of π(π‘) and πΌ(π‘). (a) When the parameter values are π = 0.13, π = 0.3, and π = 0.2, π
0 = 0.9029 < 1, π
π = 1.3 > 1. (b) When the parameter values are π = 0.3, π = 0.7, and π = 0.2, π
0 = 2.1067 and π > π + π. The value of other parameters can be seen in Table 1.
the initial time and repeating the processes as in (33)β(35), we can have that (π(π‘), πΈ(π‘), πΌ(π‘))π > 0, for all π‘ > π1 π. Thus, (π(π‘), πΈ(π‘), πΌ(π‘)) β π0 , for all π‘ > π1 π. There appears a contradiction, which means that the equality (37) holds. Therefore, πΈ0 is acyclic in ππ0 . Obviously, when π
0 > 1 and π > π + π, π is acyclic in ππ0 . Furthermore, by Lemma 5, πΈ0 = (π0 , 0, 0) and π = (0, 0, 0) are isolated invariant sets in π, ππ (πΈ0 ) β© π0 = 0, and ππ (π) β© π0 = 0. By Theorem 1.3.1 and Remark 1.3.1 in [46], it can be obtained that π is uniformly persistent with respect to (π0 , ππ0 ); that is, there exists a πΏ > 0 such that any solution (π(π‘), πΈ(π‘), πΌ(π‘)) of system (11) with the initial value (π(0), πΈ(0), πΌ(0)) β {(π, πΈ, πΌ) β π : πΈ > 0, πΌ > 0} satisfies lim inf πΈ (π‘) β₯ πΏ, π‘ββ
lim inf πΌ (π‘) β₯ πΏ. π‘ββ
(38)
Applying Theorem 1.3.6 in [46], π has a fixed point (πβ (0) , πΈβ (0) , πΌβ (0)) β π0 .
(39)
From (33), we know πβ > 0, for all π‘ β [0, π]. πβ (π‘) is also more than zero for all π‘ > 0 due to the periodicity. Similarly, for all π‘ β₯ 0, πΈβ (π‘) > 0, πΌβ (π‘) > 0. Therefore, it can be obtained that one of the positive π-periodic solutions of system (11) is (πβ (π‘), πΈβ (π‘), πΌβ (π‘)). 3.4. Numerical Simulations. Firstly, we give some notations. If π(π‘) is a periodic function with period π, we define π = π (1/π) β«0 π(π‘)ππ‘, ππ = minπ‘β[0,π] π(π‘), ππ’ = maxπ‘β[0,π] π(π‘). As described in the previous section, π
1 (π‘) = 1 +
π (π + π) (π + π) (π
0 β 1) . ππ½ (π‘) (π + π + π)
(40)
So π
1π = 1 + π(π + π)(π + π)(π
0π β 1)/ππ½π’ (π + π + π), π
1π’ = 1 + π(π + π)(π + π)(π
0π’ β 1)/ππ½π (π + π + π). In this section, we adopt π½(π‘) = π[1 + π sin(ππ‘/10)]. Then, applying the numerical simulation to verify the above solution, we give the following conclusion: (1) when π
π < 1, π is stable; (2) when π
π > 1 and π
0 < 1, πΈ0 is stable; see Figure 3(a); (3) when π
0 > 1 and π > π + π, system (11) has at least one positive periodic solution; see Figure 3(b). We can give more results about the conditions of existence of the positive periodic solution. (1σΈ ) When π
π > 1, π
0 > 1, and π
π < π
1π’ , π is stable; see Figures 4(a) and 4(b). (2σΈ ) When π
π > 1, π
0 > 1, and π
π > π
1π’ , system (11) has at least one positive periodic solution, see Figure 5. By numerical simulations, we can give that the conditions which ensure the existence of positive periodic solution are π
0 > 1 and π > π + π or π
π > 1, π
0 > 1, and π
π > π
1π’ . In fact, π
π > 1, π
0 > 1, and π
π > π
1π’ are the sufficient conditions for π
0 > 1 and π > π + π. As a result, the conditions π
0 > 1 and π > π + π 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
8
The Scientific World Journal Γ104 2.5
2
2
1.5
1.5 I(t)
I(t)
Γ104 2.5
1
1
0.5
0.5
0
0
0.5
1
1.5
2
2.5
3
S(t)
3.5 Γ104
0
0
0.5
1
1.5
2
2.5
S(t)
S(0) = 25000, E(0) = 10000, I(0) = 24000 S(0) = 35000, E(0) = 20000, I(0) = 15000 S(0) = 15000, E(0) = 10000, I(0) = 20000
S(0) = 25000, E(0) = 10000, I(0) = 24000 S(0) = 35000, E(0) = 20000, I(0) = 15000 S(0) = 15000, E(0) = 10000, I(0) = 20000
(a)
(b)
3
3.5 Γ104
Figure 4: Phase plane of π(π‘) and πΌ(π‘). (a) When the parameter values are π = 0.11, π = 0.9, π = 0.2 and π = 0.3, π
0 = 1.4991, π
π = 1.1 > 1, π
π < π
1π’ = 1.4159, and π
π < π
1π = 1.2773. (b) When the parameter values are π = 0.126, π = 0.7, π = 0.2, and π = 0.3, π
0 = 2.1067, π
π = 1.26 > 1, π
π < π
1π’ = 1.2964, and π
π > π
1π = 1.1976. The value of other parameters can be seen in Table 1.
Γ104 7000 6000
I(t)
5000 4000 3000 2000 1000 0
0.5
1
1.5 S(t)
2 Γ104
S(0) = 22000, E(0) = 5000, I(0) = 1000 S(0) = 4000, E(0) = 2000, I(0) = 1000 S(0) = 10000, E(0) = 10000, I(0) = 6000
Figure 5: Phase plane of π(π‘) and πΌ(π‘). When the parameter values are π = 0.13, π = 0.36, and π = 0.2, π
0 = 1.0834, π
π = 1.3 > 1, π
π > π
1π’ = 1.0435, and π
π > π
1π = 1.029. The value of other parameters can be seen in Table 1.
be determined by two thresholds π
0 and π
π . Only π
0 > 1 cannot promise the existence of the endemic equilibrium which also needs π
π > π
1 . When π
0 > 1 and π
π < π
1 , the solutions of the system (1) will tend to the origin π. 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 diseasefree 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 π
0 > 1 and π > π + π 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 π
0 and the average basic reproduction number π
0 which has been adopted by [27, 56β 59]. In this paper, we can calculate the average reproduction number π
0 =
π½π , (π + π) (π + π)
(41)
The Scientific World Journal 20
where π½ = (1/20) β«0 π½(π‘)ππ‘. When π = 0.13, π = 0.3, and π = 0.2, we know that π
0 = 0.9029 and π
0 = 1. When π = 0.13, π = 100000, π = 0.36, π = 0.2, and π = 0.1, π = 0.2, π = 0.1, then π
0 = 1.0834 and π
0 = 1.2. In that sense, it is confirmed that the basic reproduction number π
0 defined by [40] is more accurate than the average reproduction number π
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.
Acknowledgments The research was partially supported by the National Natural Science Foundation of China under Grant nos. 11301490, 11301491, 11331009, 11147015, 11171314, 11101251, and 11105024, Natural Science Foundation of ShanXi Province Grant no. 2012021002-1, and the Opening Foundation of Institute of Information Economy, Hangzhou Normal University, Grant no. PD12001003002003.
References [1] K. Dietz and J. A. P. Heesterbeek, βBernoulli was ahead of modern epidemiology,β Nature, vol. 408, no. 6812, pp. 513β514, 2000. [2] K. Dietz and J. A. P. Heesterbeek, βDaniel Bernoulliβs epidemiological model revisited,β Mathematical Biosciences, vol. 180, pp. 1β21, 2002. [3] R. Ross, βAn application of the theory of probabilities to the study of a priori pathometry (part I),β Proceedings of the Royal Society A, vol. 92, no. 638, pp. 204β230, 1916. [4] R. Ross and H. P. Hudson, βAn application of the theory of probabilities to the study of a priori pathometry (part II),β Proceedings of the Royal Society A, vol. 93, no. 650, pp. 212β225, 1917. [5] W. O. Kermack and A. G. McKendrick, βA contribution to the mathematical theory of epidemics (part I),β Proceedings of the Royal Society A, vol. 115, no. 772, pp. 700β721, 1927. [6] W. O. Kermack and A. G. McKendrick, βA contribution to the mathematical theory of epidemics (part II),β Proceedings of the Royal Society A, vol. 138, no. 834, pp. 55β83, 1932. [7] A. Tsoularis and J. Wallace, βAnalysis of logistic growth models,β Mathematical Biosciences, vol. 179, no. 1, pp. 21β55, 2002. [8] I. NΛasell, βOn the quasi-stationary distribution of the stochastic logistic epidemic,β Mathematical Biosciences, vol. 156, no. 1-2, pp. 21β40, 1999. [9] H. Fujikawa, A. Kai, and S. Morozumi, βA new logistic model for Escherichia coli growth at constant and dynamic temperatures,β Food Microbiology, vol. 21, no. 5, pp. 501β509, 2004. [10] L. Berezansky and E. Braverman, βOscillation properties of a logistic equation with several delays,β Journal of Mathematical Analysis and Applications, vol. 247, no. 1, pp. 110β125, 2000. [11] M. Tabata, N. Eshima, and I. Takagi, βThe nonlinear integropartial differential equation describing the logistic growth of human population with migration,β Applied Mathematics and Computation, vol. 98, pp. 169β183, 1999.
9 [12] L. Korobenko and E. Braverman, βOn logistic models with a carrying capacity dependent diffusion: stability of equilibria and coexistence with a regularly diffusing population,β Nonlinear Analysis: Real World Applications, vol. 13, no. 6, pp. 2648β2658, 2012. [13] Y. Muroya, βGlobal attractivity for discrete models of nonautonomous logistic equations,β Computers and Mathematics with Applications, vol. 53, no. 7, pp. 1059β1073, 2007. [14] S. Invernizzi and K. Terpin, βA generalized logistic model for photosynthetic growth,β Ecological Modelling, vol. 94, no. 2-3, pp. 231β242, 1997. [15] Z. Min, W. Bang-Jun, and J. Feng, βCoalmining citiesβ economic growth mechanism and sustainable development analysis based on logistic dynamics model,β Procedia Earth and Planetary Science, vol. 1, no. 1, pp. 1737β1743, 2009. [16] L. I. AnitΒΈa, S. AnitΒΈa, and V. ArnΛautu, βGlobal behavior for an age-dependent population model with logistic term and periodic vital rates,β Applied Mathematics and Computation, vol. 206, no. 1, pp. 368β379, 2008. [17] Z. A. Banaszak, X. Q. Tang, S. C. Wang, and M. B. Zaremba, βLogistics models in flexible manufacturing,β Computers in Industry, vol. 43, no. 3, pp. 237β248, 2000. [18] S. Brianzoni, C. Mammana, and E. Michetti, βNonlinear dynamics in a business-cycle model with logistic population growth,β Chaos, Solitons and Fractals, vol. 40, no. 2, pp. 717β730, 2009. [19] W. P. London and J. A. Yorke, βRecurrent outbreaks of measles, chickenpox and mumps. I. Seasonal variation in contact rates,β The American Journal of Epidemiology, vol. 98, no. 6, pp. 453β 468, 1973. [20] S. F. Dowell, βSeasonal variation in host susceptibility and cycles of certain infectious diseases,β Emerging Infectious Diseases, vol. 7, no. 3, pp. 369β374, 2001. [21] O. N. BjΓΈrnstad, B. F. FinkenstΒ¨adt, and B. T. Grenfell, βDynamics of measles epidemics: estimating scaling of transmission rates using a Time series SIR model,β Ecological Monographs, vol. 72, no. 2, pp. 169β184, 2002. [22] J. Zhang, Z. Jin, G. Q. Sun, X. Sun, and S. Ruan, βModeling seasonal rabies epidemics in China,β Bulletin of Mathematical Biology, vol. 74, no. 5, pp. 1226β1251, 2012. [23] I. B. Schwartz, βSmall amplitude, long period outbreaks in seasonally driven epidemics,β Journal of Mathematical Biology, vol. 30, no. 5, pp. 473β491, 1992. [24] I. B. Schwartz and H. L. Smith, βInfinite subharmonic bifurcation in an SEIR epidemic model,β Journal of Mathematical Biology, vol. 18, no. 3, pp. 233β253, 1983. [25] H. L. Smith, βMultiple stable subharmonics for a periodic epidemic model,β Journal of Mathematical Biology, vol. 17, no. 2, pp. 179β190, 1983. [26] J. Dushoff, J. B. Plotkin, S. A. Levin, and D. J. D. Earn, βDynamical resonance can account for seasonality of influenza epidemics,β Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 48, pp. 16915β16916, 2004. [27] J. L. Ma and Z. Ma, βEpidemic threshold conditions for seasonally forced SEIR models,β Mathematical Biosciences and Engineering, vol. 3, no. 1, pp. 161β172, 2006. [28] D. J. D. Earn, P. Rohani, B. M. Bolker, and B. T. Grenfell, βA simple model for complex dynamical transitions in epidemics,β Science, vol. 287, no. 5453, pp. 667β670, 2000.
10 [29] M. J. Keeling, P. Rohani, and B. T. Grenfell, βSeasonally forced disease dynamics explored as switching between attractors,β Physica D, vol. 148, no. 3-4, pp. 317β335, 2001. [30] J. L. Aron and I. B. Schwartz, βSeasonality and period-doubling bifurcations in an epidemic model,β Journal of Theoretical Biology, vol. 110, no. 4, pp. 665β679, 1984. [31] N. C. Grassly and C. Fraser, βSeasonal infectious disease epidemiology,β Proceedings of the Royal Society B, vol. 273, no. 1600, pp. 2541β2550, 2006. [32] J. Zhang, Z. Jin, G. Q. Sun, T. Zhou, and S. Ruan, βAnalysis of rabies in China: transmission dynamics and control,β PLoS ONE, vol. 6, no. 7, Article ID e20891, 2011. [33] O. Diekmann, J. A. P. Heesterbeek, and J. A. Metz, βOn the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations,β Journal of mathematical biology, vol. 28, no. 4, pp. 365β382, 1990. [34] O. Diekmann, J. A. P. Heesterbeek, and M. G. Roberts, βThe construction of next-generation matrices for compartmental epidemic models,β Journal of the Royal Society Interface, vol. 7, no. 47, pp. 873β885, 2010. [35] P. van den Driessche and J. Watmough, βReproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,β Mathematical Biosciences, vol. 180, pp. 29β48, 2002. [36] F. Berezovsky, G. Karev, B. J. Song, and C. C. Chavez, βA simple epidemic model with surprising dynamics,β Mathematical Biosciences and Engineering, vol. 2, pp. 133β152, 2005. [37] J. LaSalle and S. Lefschetz, Stability by Liapunovβs Direct Method, Academic Press, New York, NY, USA, 1961. [38] E. A. Barbashin, Introduction to the Theory of Stability, WoltersNoordhoff, Groningen, The Netherlands, 1970. [39] D. Schenzle, βAn age-structured model of pre- and postvaccination measles transmission,β Mathematical Medicine and Biology, vol. 1, no. 2, pp. 169β191, 1984. [40] N. BacaΒ¨er and S. Guernaoui, βThe epidemic threshold of vector-borne diseases with seasonality: the case of cutaneous leishmaniasis in Chichaoua, Morocco,β Journal of Mathematical Biology, vol. 53, no. 3, pp. 421β436, 2006. [41] W. D. Wang and X. Q. Zhao, βThreshold dynamics for compartmental epidemic models in periodic environments,β Journal of Dynamics and Differential Equations, vol. 20, no. 3, pp. 699β717, 2008. [42] F. Zhang and X. Q. Zhao, βA periodic epidemic model in a patchy environment,β Journal of Mathematical Analysis and Applications, vol. 325, no. 1, pp. 496β516, 2007. [43] H. L. Smith and P. Waltman, The Theory of the Chemostat, Cambridge University Press, New York, NY, USA, 1995. [44] H. R. Thieme, βConvergence results and a PoincarΒ΄e-Bendixson trichotomy for asymptotically automous differential equations,β Journal of Mathematical Biology, vol. 30, pp. 755β763, 1992. [45] L. Perko, Differential Equations and Dynamical Systems, Springer, New York, NY, USA, 2000. [46] X. Q. Zhao, Dynamical Systems in Population Biology, Springer, New York, NY, USA, 2003. [47] Z. G. Bai and Y. C. Zhou, βThreshold dynamics of a bacillary dysentery model with seasonal fluctuation,β Discrete and Continuous Dynamical Systems B, vol. 15, no. 1, pp. 1β14, 2011. [48] Z. G. Bai, Y. C. Zhou, and T. L. Zhang, βExistence of multiple periodic solutions for an SIR model with seasonality,β Nonlinear
The Scientific World Journal
[49]
[50]
[51]
[52]
[53]
[54]
[55]
[56]
[57]
[58]
[59]
Analysis: Theory, Methods and Applications, vol. 74, no. 11, pp. 3548β3555, 2011. N. BacaΒ¨er, βApproximation of the basic reproduction number R0 for vector-borne diseases with a periodic vector population,β Bulletin of Mathematical Biology, vol. 69, no. 3, pp. 1067β1091, 2007. N. BacaΒ¨er and X. Abdurahman, βResonance of the epidemic threshold in a periodic environment,β Journal of Mathematical Biology, vol. 57, no. 5, pp. 649β673, 2008. N. BacaΒ¨er, βPeriodic matrix population models: growth rate, basic reproduction number, and entropy,β Bulletin of Mathematical Biology, vol. 71, no. 7, pp. 1781β1792, 2009. N. BacaΒ¨er and E. H. A. Dads, βGenealogy with seasonality, the basic reproduction number, and the influenza pandemic,β Journal of Mathematical Biology, vol. 62, no. 5, pp. 741β762, 2011. L. J. Liu, X. Q. Zhao, and Y. C. Zhou, βA tuberculosis model with seasonality,β Bulletin of Mathematical Biology, vol. 72, no. 4, pp. 931β952, 2010. J. L. Liu, βThreshold dynamics for a HFMD epidemic model with periodic transmission rate,β Nonlinear Dynamics, vol. 64, no. 1-2, pp. 89β95, 2011. Y. Nakata and T. Kuniya, βGlobal dynamics of a class of SEIRS epidemic models in a periodic environment,β Journal of Mathematical Analysis and Applications, vol. 363, no. 1, pp. 230β 237, 2010. B. G. Williams and C. Dye, βInfectious disease persistence when transmission varies seasonally,β Mathematical Biosciences, vol. 145, no. 1, pp. 77β88, 1997. I. A. Moneim, βThe effect of using different types of periodic contact rate on the behaviour of infectious diseases: a simulation study,β Computers in Biology and Medicine, vol. 37, no. 11, pp. 1582β1590, 2007. C. L. Wesley and L. J. S. Allen, βThe basic reproduction number in epidemic models with periodic demographics,β Journal of Biological Dynamics, vol. 3, no. 2-3, pp. 116β129, 2009. D. Greenhalgh and I. A. Moneim, βSIRS epidemic model and simulations using different types of seasonal contact rate,β Systems Analysis Modelling Simulation, vol. 43, no. 5, pp. 573β 600, 2003.