Research Article  Open Access
Positive Periodic Solutions of an Epidemic Model with Seasonality
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 number is obtained. Moreover, only the basic reproduction number cannot ensure the existence of the positive equilibrium, which needs additional condition . 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 discretetime 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 businesscycle model with logistic population 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 , the population grows; that is, a critical mass of individuals for the disease to spread may be supported. If , 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 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 nonautonomous 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 where is the amplitude of seasonal variation in transmission (typically referred to as the “strength of seasonal forcing”) and is the period, which is a crude assumption for many infectious diseases [29–31]. When , there is no nonseasonal infections. Motivated by biological realism, some recent papers take the contact rate as , 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 [29]. Here, we take the form .
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 nonautonomous system in terms of global asymptotic stability of the diseasefree 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 incidence which is fit for the longterm 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, handfootmouth disease, tuberculosis, and AIDS [22, 32]. The model we employ is as follows: where all parameters are positive whose interpretations can be seen in Table 1.

Noticing the equations in model (1), we have
When there exists no disease, we have
Let , if , for , as ; that is, the population will grow and tend to a steady state . If , then which will cause the population to disappear. Thus, is the basic demographic reproduction number. From the above equation, the feasible region can be obtained: , 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: where is the origin, is the diseasefree equilibrium, and is the endemic equilibrium. Concretely, one can have where is the basic demographic reproduction number and is the basic reproduction number which can be obtained by the nextgeneration 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 and , where . So, we can obtain the following theorems.
Theorem 2. The system (1) has three equilibria: origin , diseasefree equilibrium , and the endemic equilibrium . always exists; if , exists; if and , exists.
Theorem 3. When and , is globally asymptotically stable.
Proof. By [33–35], we know that is locally asymptotically stable. Now we define a Lyapunov function
When and , the Lyapunov function satisfies
Moreover, only hold when . It is easy to verify that the diseasefree equilibrium point is the only fixed point of the system. Hence, applying the LyapunovLaSalle asymptotic stability theorem in [37, 38], the diseasefree equilibrium point 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 , is stable; see Figures 2(b) and 2(c).(2)When and , is stable; see Figure 2(d).(3)When , , and , is stable; see Figure 2(a).(4)When , , and , is stable; see Figure 2(e).
(a)
(b)
(c)
(d)
(e)
3. Nonautonomous Model and Analysis
3.1. The Basic Reproduction Number
Now, we consider the nonautonomous case of the model (1) when the transmission rate is periodic, which is given as follows: 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 . According to the basic reproduction number under nonautonomous system, we can refer to the method of [40, 41]. From the last section, we know that system (11) has one diseasefree equilibrium , where . By giving a new vector , we have
Taking the partial derivative of the above vectors about variables , and substituting the diseasefree equilibrium, we have
According to [41], denote to be the ordered Banach space of all periodic functions from to which is equipped with the maximum norm and the positive cone . Over the Banach space, we define a linear operator by 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
In order to give the expression of the basic reproduction number, we need to introduce the linear periodic system with parameter . Let , , be the evolution operator of system (16) on . In fact, , and , for all . By Theorems 2.1 and 2.2 in [41], the basic reproduction number also can be defined as such that , which can be straightforward to calculate.
3.2. Global Stability of the DiseaseFree Equilibrium
Theorem 4. The diseasefree equilibrium is globally asymptotically stable when and .
Proof. Theorem 2.2 in [41] implies that is locally asymptotically stable when and . So we only need to prove its global attractability. It is easy to know that . Thus,
The right comparison system can be written as
that is,
For (19), Lemma 2.1 in [42] shows that there is a positive periodic function such that is a solution of system (18), where . By Theorem 2.2 in [41], we know that when and , and , which implies 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 and , is globally attractive. Therefore, the proposition that 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 be the solution of system (11) with the initial value . By the fundamental existenceuniqueness theorem [45], is the unique solution of system (11) with .
Next, we need to introduce the Poincaré map associated with system (11); that is, where is the period. Theorem 1 implies that is positively invariant and is a dissipative point.
Now, we introduce two subsets of , and .
Lemma 5. (a) When and , there exists a such that when
for any , one has
where .
(b) When and , there exists a such that when
for any , one has
where .
Proof. (a) By Theorem 2.2 in [41], we know that when , . So there is a small enough positive number such that , where
If proposition (a) does not hold, there is some such that
We can assume that for all , . Applying the continuity of the solutions with respect to the initial values,
Let , where and . is the greatest integer which is not more than . Then, for any ,
So . Then, when ,
Thus, we can study the right linear system
For the system (30), there exists a positive periodic function such that is a solution of system (11), where . When , , which means that when , as . By the comparison principle [43], when , , , as . There appears a contradiction. Thus, the proposition (a) holds.
(b) When and , we have
So if , for as ; that is, .
Theorem 6. When and , there exists a such that any solution of system (11) with initial value satisfies and system (11) has at least one positive periodic solution.
Proof. From system (11),
for any , which shows that is positively invariant. Moreover, it is obvious to see that is relatively closed in . Denote
Next, we prove that
We only need to show that , which means that for any , , for all . If it does not hold, there exists a such that . Taking as the initial time and repeating the processes as in (33)–(35), we can have that , for all . Thus, , for all . There appears a contradiction, which means that the equality (37) holds. Therefore, is acyclic in . Obviously, when and , is acyclic in .
Furthermore, by Lemma 5, and are isolated invariant sets in , , and . By Theorem and Remark in [46], it can be obtained that is uniformly persistent with respect to ; that is, there exists a such that any solution of system (11) with the initial value satisfies
Applying Theorem in [46], has a fixed point
From (33), we know , for all . is also more than zero for all due to the periodicity. Similarly, for all , , . 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 , , . As described in the previous section, So , .
In this section, we adopt . Then, applying the numerical simulation to verify the above solution, we give the following conclusion:(1)when , is stable;(2)when and , is stable; see Figure 3(a);(3)when and , system (11) has at least one positive periodic solution; see Figure 3(b).
(a)
(b)
We can give more results about the conditions of existence of the positive periodic solution. When , , and , is stable; see Figures 4(a) and 4(b). When , , and , system (11) has at least one positive periodic solution, see Figure 5.
(a)
(b)
By numerical simulations, we can give that the conditions which ensure the existence of positive periodic solution are and or , , and . In fact, , , and are the sufficient conditions for and . As a result, the conditions and are broader.
4. Discussion
This paper considers a logistic growth system whose birth process incorporates densitydependent 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 thresholds and . Only cannot promise the existence of the endemic equilibrium which also needs . When and , 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 nonautonomous 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 diseasefree 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 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 and the average basic reproduction number which has been adopted by [27, 56–59]. In this paper, we can calculate the average reproduction number where . When , , and , we know that and . When , , , , and , , , then and . In that sense, it is confirmed that the basic reproduction number defined by [40] is more accurate than the average reproduction number 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. 20120210021, and the Opening Foundation of Institute of Information Economy, Hangzhou Normal University, Grant no. PD12001003002003.
References
 K. Dietz and J. A. P. Heesterbeek, “Bernoulli was ahead of modern epidemiology,” Nature, vol. 408, no. 6812, pp. 513–514, 2000. View at: Google Scholar
 K. Dietz and J. A. P. Heesterbeek, “Daniel Bernoulli's epidemiological model revisited,” Mathematical Biosciences, vol. 180, pp. 1–21, 2002. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 A. Tsoularis and J. Wallace, “Analysis of logistic growth models,” Mathematical Biosciences, vol. 179, no. 1, pp. 21–55, 2002. View at: Publisher Site  Google Scholar
 I. Nåsell, “On the quasistationary distribution of the stochastic logistic epidemic,” Mathematical Biosciences, vol. 156, no. 12, pp. 21–40, 1999. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 Y. Muroya, “Global attractivity for discrete models of nonautonomous logistic equations,” Computers and Mathematics with Applications, vol. 53, no. 7, pp. 1059–1073, 2007. View at: Publisher Site  Google Scholar
 S. Invernizzi and K. Terpin, “A generalized logistic model for photosynthetic growth,” Ecological Modelling, vol. 94, no. 23, pp. 231–242, 1997. View at: Publisher Site  Google Scholar
 Z. Min, W. BangJun, 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. View at: Publisher Site  Google Scholar
 L. I. Aniţa, S. Aniţa, and V. Arnǎutu, “Global behavior for an agedependent population model with logistic term and periodic vital rates,” Applied Mathematics and Computation, vol. 206, no. 1, pp. 368–379, 2008. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 S. Brianzoni, C. Mammana, and E. Michetti, “Nonlinear dynamics in a businesscycle model with logistic population growth,” Chaos, Solitons and Fractals, vol. 40, no. 2, pp. 717–730, 2009. View at: Publisher Site  Google Scholar
 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. View at: Google Scholar
 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. View at: Google Scholar
 O. N. Bjørnstad, B. F. Finkenstädt, 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. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 I. B. Schwartz, “Small amplitude, long period outbreaks in seasonally driven epidemics,” Journal of Mathematical Biology, vol. 30, no. 5, pp. 473–491, 1992. View at: Google Scholar
 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. View at: Google Scholar
 H. L. Smith, “Multiple stable subharmonics for a periodic epidemic model,” Journal of Mathematical Biology, vol. 17, no. 2, pp. 179–190, 1983. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 M. J. Keeling, P. Rohani, and B. T. Grenfell, “Seasonally forced disease dynamics explored as switching between attractors,” Physica D, vol. 148, no. 34, pp. 317–335, 2001. View at: Publisher Site  Google Scholar
 J. L. Aron and I. B. Schwartz, “Seasonality and perioddoubling bifurcations in an epidemic model,” Journal of Theoretical Biology, vol. 110, no. 4, pp. 665–679, 1984. View at: Google Scholar
 N. C. Grassly and C. Fraser, “Seasonal infectious disease epidemiology,” Proceedings of the Royal Society B, vol. 273, no. 1600, pp. 2541–2550, 2006. View at: Publisher Site  Google Scholar
 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. View at: Google Scholar
 O. Diekmann, J. A. P. Heesterbeek, and J. A. Metz, “On the definition and the computation of the basic reproduction ratio R_{0} in models for infectious diseases in heterogeneous populations,” Journal of mathematical biology, vol. 28, no. 4, pp. 365–382, 1990. View at: Google Scholar
 O. Diekmann, J. A. P. Heesterbeek, and M. G. Roberts, “The construction of nextgeneration matrices for compartmental epidemic models,” Journal of the Royal Society Interface, vol. 7, no. 47, pp. 873–885, 2010. View at: Publisher Site  Google Scholar
 P. van den Driessche and J. Watmough, “Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, pp. 29–48, 2002. View at: Publisher Site  Google Scholar
 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. View at: Google Scholar
 J. LaSalle and S. Lefschetz, Stability by Liapunov's Direct Method, Academic Press, New York, NY, USA, 1961.
 E. A. Barbashin, Introduction to the Theory of Stability, WoltersNoordhoff, Groningen, The Netherlands, 1970.
 D. Schenzle, “An agestructured model of pre and postvaccination measles transmission,” Mathematical Medicine and Biology, vol. 1, no. 2, pp. 169–191, 1984. View at: Publisher Site  Google Scholar
 N. Bacaër and S. Guernaoui, “The epidemic threshold of vectorborne diseases with seasonality: the case of cutaneous leishmaniasis in Chichaoua, Morocco,” Journal of Mathematical Biology, vol. 53, no. 3, pp. 421–436, 2006. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 H. L. Smith and P. Waltman, The Theory of the Chemostat, Cambridge University Press, New York, NY, USA, 1995.
 H. R. Thieme, “Convergence results and a PoincaréBendixson trichotomy for asymptotically automous differential equations,” Journal of Mathematical Biology, vol. 30, pp. 755–763, 1992. View at: Google Scholar
 L. Perko, Differential Equations and Dynamical Systems, Springer, New York, NY, USA, 2000.
 X. Q. Zhao, Dynamical Systems in Population Biology, Springer, New York, NY, USA, 2003.
 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. View at: Publisher Site  Google Scholar
 Z. G. Bai, Y. C. Zhou, and T. L. Zhang, “Existence of multiple periodic solutions for an SIR model with seasonality,” Nonlinear Analysis: Theory, Methods and Applications, vol. 74, no. 11, pp. 3548–3555, 2011. View at: Publisher Site  Google Scholar
 N. Bacaër, “Approximation of the basic reproduction number R_{0} for vectorborne diseases with a periodic vector population,” Bulletin of Mathematical Biology, vol. 69, no. 3, pp. 1067–1091, 2007. View at: Publisher Site  Google Scholar
 N. Bacaër and X. Abdurahman, “Resonance of the epidemic threshold in a periodic environment,” Journal of Mathematical Biology, vol. 57, no. 5, pp. 649–673, 2008. View at: Publisher Site  Google Scholar
 N. Bacaër, “Periodic matrix population models: growth rate, basic reproduction number, and entropy,” Bulletin of Mathematical Biology, vol. 71, no. 7, pp. 1781–1792, 2009. View at: Publisher Site  Google Scholar
 N. Bacaër 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 J. L. Liu, “Threshold dynamics for a HFMD epidemic model with periodic transmission rate,” Nonlinear Dynamics, vol. 64, no. 12, pp. 89–95, 2011. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 B. G. Williams and C. Dye, “Infectious disease persistence when transmission varies seasonally,” Mathematical Biosciences, vol. 145, no. 1, pp. 77–88, 1997. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. 23, pp. 116–129, 2009. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 GuiQuan 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.