Discrete Dynamics in Nature and Society

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

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

Research Article

Bifurcation and Chaotic Behavior of a Discrete-Time SIS Model

Junhong Li and Ning Cui

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).

(a)

(b)

(c)

(d)

(e)

Figure 1: (a) Bifurcation diagram in space. (b) The bifurcation diagram in - plane and in - plane. (c) Spectrum of Lyapunov exponents corresponding to (b). (d) The attractor of (5) with (e) Time series of and of the same parameters.

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.

(a)

(b)

Figure 2

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 .

(a)

(b)

(c)

(d)

Figure 3

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.

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.

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.

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.

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.

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.

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.

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.

D. Gao and S. Ruan, “An SIS patch model with variable transmission coefficients,” Mathematical Biosciences, vol. 232, no. 2, pp. 110–115, 2011.

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.

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.

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.

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.

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.

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.

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.

J. Guckenheimer and P. Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields, vol. 42, Springer, New York, NY, USA, 1983.

S. Wiggins, Introduction to Applied Nonlinear Dynamical Systems and Chaos, vol. 2, Springer, New York, NY, USA, 1990.

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.

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.

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.

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.

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.

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.

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.