Discrete Dynamics in Nature and Society

Volume 2013 (2013), Article ID 705601, 8 pages

http://dx.doi.org/10.1155/2013/705601

## Bifurcation and Chaotic Behavior of a Discrete-Time SIS Model

Department of Mathematics and Sciences, Hebei Institute of Architecture and Civil Engineering, Zhangjiakou, Hebei 075000, China

Received 5 January 2013; Accepted 8 April 2013

Academic Editor: Qingdu Li

Copyright © 2013 Junhong Li and Ning Cui. 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.

#### Abstract

The discrete-time epidemic model is investigated, which is obtained using the Euler method. It is verified that there exist some dynamical behaviors in this model, such as transcritical bifurcation, flip bifurcation, Hopf bifurcation, and chaos. The numerical simulations, including bifurcation diagrams and computation of Lyapunov exponents, not only show the consistence with the theoretical analysis but also exhibit the rich and complex dynamical behaviors.

#### 1. Introduction

Epidemic models have been widely used in different forms for studying epidemiological processes such as the spread of HIV [1], SARS [2], and influenza [3]. It is well known that dynamical systems with simple dynamical behavior in the constant parameter case display very complex behaviors including chaos when they are periodically perturbed [4, 5]. The continuous-time epidemic models have been widely investigated in many articles (e.g., [6–10] and the references cited therein). In recent years, we have found that more attention is paid to the discrete-time epidemic models (see [11–15] and the references cited therein). The reasons are as follows: first, difference models are more realistic than continuous differential ones because the epidemic statistics are compiled from given time intervals and are discontinuous. Second, the discrete-time models can provide natural simulators for the continuous cases. One can thus not only study the behaviors of the continuous-time model with good accuracy, but also assess the effect of lager time steps. At last, the use of discrete -time models makes it possible to use the entire arsenal of methods recently developed for the study of mappings and lattice equations, either from the integrability and/or chaos points of view.

On the other hand, daily treatments are frequently done for some infections, such as the group of those being responsible for the common cold, which do not confer any long lasting immunity. Such infections do not have a recovered state and individuals become susceptible again after infection. For such reasons, according to [16], we firstly consider the SIS epidemic model with nonlinear incidence rate: where denotes the susceptible population, is the infected population, and is the intrinsic birth rate constant. , are the carrying capacity and the infection rate, respectively. is the recovery rate constant ( is the average infective time).

Let We obtain the following system analogous to (1): where Applying Euler scheme to (3), we obtain the following equation:

This paper is organized as follows. In Section 2, we give sufficient conditions of existence for transcritical bifurcation, flip bifurcation, and Hopf bifurcation. In Section 3, a series of numerical simulations show that there are bifurcation and chaos in the discrete-time epidemic model. Finally, we give remarks to conclude this paper in Section 4.

#### 2. Bifurcations

It is easy to visualize that system (5) has three fixed points , , and when , are fixed. We can see that the fixed point is a saddle. In the following paper, we focus on investigating the bifurcations of , .

Theorem 1. *If and , (5) undergo a transcritical bifurcation at . Furthermore, the system has three fixed points when and has two fixed points when .*

*Proof. *The Jacobian matrix of (5) at takes the form
has eigenvalues , . And implies that .

Let
equation (5) becomes
where
By the following transformation:
equation (8) becomes
where
Then, we can consider
which must satisfy
Thus, we have
And (8) is restricted to the center manifold, which is given by
Since
system (5) undergoes a transcritical bifurcation at . This proves the theorem.

Theorem 2. *Equation (5) undergoes a flip bifurcation at when , . Furthermore, the stable periodic-2 point bifurcates from this fixed point.*

*Proof. *Let
system (5) becomes
where
By the following transformation:
equation (19) becomes
where
Then, we can consider , which must satisfy
Thus, we have , . We obtain the center manifold as follows:
And (19) is restricted to the center manifold, which is given by
Direct calculations show that
Hence, system (5) undergoes a flip bifurcation at . This completes the proof.

The positive fixed point is so important to the biological system that people usually are very interested in it. We will next pay attention to the only positive fixed point of (5).

Theorem 3. *Equation (5) undergoes a transcritical bifurcation at when , , and , where
*

Theorem 4. *If , , and , (5) undergoes a flip bifurcation at .*

Since the analysis is similar to the case at , the previously mentioned proofs are omitted.

We next give the condition of existence of Hopf bifurcation by using the hopf bifurcation theorem in [17].

The characteristic equation of the Jacobian matrix can be written as , where

The eigenvalues of the Jacobian matrix of (5) at are . The eigenvalues are complex conjugates for . We translate the fixed point to the origin by , , and the system (5) becomes The eigenvalues of the matrix associated with the linearized map (30) at fixed point are complex conjugates which are written as , where Now assume that , and let . Then we have Let The system (30) becomes where Notice that (34) is exactly on the center manifold in the form, in which the coefficient [18] is given by , where From the previous analysis, we have Theorem 5.

Theorem 5. *System (5) undergoes a Hopf bifurcation at fixed point , if , and
*

#### 3. Numerical Simulations

With development of scientific computation, computer becomes a powerful tool to study nonlinear systems, especially for a system without explicit solution. It not only is able to explore new complex dynamical behaviors, for example, periodic orbits and chaos in different regions [16, 19], but also can do rigorous analysis by combining modern dynamical systems theory and reliable computation, for example, topological horseshoes [20–23]. In this section, we use the bifurcation diagrams, Lyapunov exponents, and phase portraits to illustrate the previous analytic results and find new dynamics of the model (5) as the parameters vary. The attractors of (5) are also given by using the method in [20]. The bifurcation parameters are considered in the following three cases:(I)varying in range and fixing , ;(II)varying in range and fixing , ;(III)varying in range and fixing , .

*For case (I).* The bifurcation diagram of (5) in - and - space for , is given in Figures 1(a) and 1(b) to show the dynamical changes of susceptible and infective, respectively, as varies. There is obvious phenomenon of bifurcation when we select the stepsize to be . The spectrum of Lyapunov exponents of the system (5) with respect to parameter is given in Figure 1(c).

Moreover, we can see that the orbit with initial values approaches to the stable fixed point for approximately, and Hopf bifurcation occurs at . When increased to , (5) becomes stable. In Figures 1(b) and 1(c), we observe the period -4,8 windows within the chaotic regions and boundary crisis at . For the maximum Lyapunov exponents are positive which correspond to chaotic region. To well see the dynamics, the attractor in the system (5) and time series of and are given in Figures 1(d) and 1(e), respectively.

*For case (II).* Figure 2(a) is the bifurcation diagram of (5) for and and depicts that there are period -4,5 windows within the chaotic regions and boundary crisis at . Figure 2(b) shows the spectrum of Lyapunov exponents of the system (5) with respect to parameter . For the maximum Lyapunov exponents are positive which correspond to chaotic region.

*For case (III).* Figures 3(a) and 3(b) are the bifurcation diagrams of (5) for and , respectively, and depict that there are period -2,4 windows within the chaotic regions and boundary crisis at . Figure 3(c) shows the spectrum of Lyapunov exponents of (5) with respect to parameter . Figure 3(d) is the attractor of (5) with .

#### 4. Conclusion

In this paper, we investigate the behaviors of a discrete-time SIS epidemic model with nonlinear incidence rate, and find many complex and new interesting dynamical phenomena. Without the recovery rate of infectious hosts, (1) becomes the SI model (see [16]). Our theoretical analysis and numerical simulations have demonstrated that the model exhibits the variety of dynamical behaviors, which includes that the discrete epidemic model undergoes transcritical bifurcation, flip bifurcation, Hopf bifurcation, and chaos. Furthermore, chaos can cause the population to run a higher risk of extinction due to the unpredictability [24, 25]. Thus, how to control chaos in the epidemic model is very important, which needs further consideration.

#### Acknowledgments

The authors are very grateful to the reviewers for their valuable comments and suggestions. This work was supported by the Youth Science Foundations of Education Department of Hebei Province (nos. 2010233, 2011236).

#### References

- L. Cai, X. Li, M. Ghosh, and B. Guo, “Stability analysis of an HIV/AIDS epidemic model with treatment,”
*Journal of Computational and Applied Mathematics*, vol. 229, no. 1, pp. 313–323, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y.-H. Hsieh, C.-C. King, C. W. S. Chen, M.-S. Ho, S.-B. Hsu, and Y.-C. Wu, “Impact of quarantine on the 2003 SARS outbreak: a retrospective modeling study,”
*Journal of Theoretical Biology*, vol. 244, no. 4, pp. 729–736, 2007. View at Publisher · View at Google Scholar · View at MathSciNet - S. Iwami, Y. Takeuchi, X. Liu, and S. Nakaoka, “A geographical spread of vaccine-resistance in avian influenza epidemics,”
*Journal of Theoretical Biology*, vol. 259, no. 2, pp. 219–228, 2009. View at Publisher · View at Google Scholar · View at MathSciNet - D. Summers, J. G. Cranford, and B. P. Healey, “Chaos in periodically forced discrete-time ecosystem models,”
*Chaos, Solitons and Fractals*, vol. 11, no. 14, pp. 2331–2342, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Vandermeer, L. Stone, and B. Blasius, “Categories of chaos and fractal basin boundaries in forced predator-prey models,”
*Chaos, Solitons and Fractals*, vol. 12, no. 2, pp. 265–276, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Gray, D. Greenhalgh, X. Mao, and J. Pan, “The SIS epidemic model with Markovian switching,”
*Journal of Mathematical Analysis and Applications*, vol. 394, no. 2, pp. 496–516, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - A. Gray, D. Greenhalgh, L. Hu, X. Mao, and J. Pan, “A stochastic differential equation SIS epidemic model,”
*SIAM Journal on Applied Mathematics*, vol. 71, no. 3, pp. 876–902, 2011. View at Publisher · View at Google Scholar · View at MathSciNet - J. Liu and T. Zhang, “Bifurcation analysis of an SIS epidemic model with nonlinear birth rate,”
*Chaos, Solitons & Fractals*, vol. 40, no. 3, pp. 1091–1099, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Gao and S. Ruan, “An SIS patch model with variable transmission coefficients,”
*Mathematical Biosciences*, vol. 232, no. 2, pp. 110–115, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. Li and J. Cui, “The effect of constant and pulse vaccination on SIS epidemic models incorporating media coverage,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 14, no. 5, pp. 2353–2365, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Q. Wu and X. Fu, “Modelling of discrete-time SIS models with awareness interactions on degree-uncorrelated networks,”
*Physica A*, vol. 390, no. 3, pp. 463–470, 2011. View at Publisher · View at Google Scholar · View at Scopus - P. L. Salceanu and H. L. Smith, “Persistence in a discrete-time stage-structured fungal disease model,”
*Journal of Biological Dynamics*, vol. 3, no. 2-3, pp. 271–285, 2009. View at Publisher · View at Google Scholar · View at MathSciNet - J. E. Franke and A.-A. Yakubu, “Disease-induced mortality in density-dependent discrete-time S-I-S epidemic models,”
*Journal of Mathematical Biology*, vol. 57, no. 6, pp. 755–790, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Li, Z. Ma, and F. Brauer, “Global analysis of discrete-time SI and SIS epidemic models,”
*Mathematical Biosciences and Engineering*, vol. 4, no. 4, pp. 699–710, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. J. S. Allen, Y. Lou, and A. L. Nevai, “Spatial patterns in a discrete-time SIS patch model,”
*Journal of Mathematical Biology*, vol. 58, no. 3, pp. 339–375, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. Li, G.-Q. Sun, and Z. Jin, “Bifurcation and chaos in an epidemic model with nonlinear incidence rates,”
*Applied Mathematics and Computation*, vol. 216, no. 4, pp. 1226–1234, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Guckenheimer and P. Holmes,
*Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields*, vol. 42, Springer, New York, NY, USA, 1983. View at MathSciNet - S. Wiggins,
*Introduction to Applied Nonlinear Dynamical Systems and Chaos*, vol. 2, Springer, New York, NY, USA, 1990. View at MathSciNet - T. Chen, J. He, and Q. Yin, “Dynamics evolution of credit risk contagion in the CRT market,”
*Dynamics in Nature and Society*, vol. 2013, Article ID 206201, 9 pages, 2013. View at Publisher · View at Google Scholar - Q. Li, L. Zhang, and F. Yang, “An algorithm to automatically detect the Smale horseshoes,”
*Discrete Dynamics in Nature and Society*, Article ID 283179, 9 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Q. Li and X.-S. Yang, “A simple method for finding topological horseshoes,”
*International Journal of Bifurcation and Chaos in Applied Sciences and Engineering*, vol. 20, no. 2, pp. 467–478, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X.-S. Yang, “Topological horseshoes and computer assisted verification of chaotic dynamics,”
*International Journal of Bifurcation and Chaos in Applied Sciences and Engineering*, vol. 19, no. 4, pp. 1127–1145, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. Qingdu and T. Song, “Algorithm for finding horseshoes in three-dimensinal hyperchaltic maps and its application,”
*Acta Physica Sinica*, vol. 62, no. 2, Article ID 020510, 2013. View at Google Scholar - M. P. Hassell, H. N. Comins, and R. M. May, “Spatial structure and chaos in insect population dynamics,”
*Nature*, vol. 353, pp. 255–258, 1991. View at Google Scholar - A. A. Berryman and J. A. Millstein, “Are ecological systems chaotic and if not, why not?”
*Trends in Ecology & Evolution*, vol. 4, no. 1, pp. 26–28, 1989. View at Google Scholar