Discrete Dynamics in Nature and Society

Volume 2015, Article ID 974868, 13 pages

http://dx.doi.org/10.1155/2015/974868

## Bifurcation Analysis and Chaos Control in a Discrete Epidemic System

^{1}Department of Mathematics and Information Science, Beifang University of Nationalities, Yinchuan, Ningxia 750021, China^{2}Department of Mathematics and Statistics, Zhoukou Normal University, Zhoukou, Henan 466001, China

Received 14 June 2015; Accepted 9 August 2015

Academic Editor: Garyfalos Papashinopoulos

Copyright © 2015 Wei Tan 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.

#### Abstract

The dynamics of discrete *SI* epidemic model, which has been obtained by the forward Euler scheme, is investigated in detail. By using the center manifold theorem and bifurcation theorem in the interior , the specific conditions for the existence of flip bifurcation and Neimark-Sacker bifurcation have been derived. Numerical simulation not only presents our theoretical analysis but also exhibits rich and complex dynamical behavior existing in the case of the windows of period-1, period-3, period-5, period-6, period-7, period-9, period-11, period-15, period-19, period-23, period-34, period-42, and period-53 orbits. Meanwhile, there appears the cascade of period-doubling 2, 4, 8 bifurcation and chaos sets from the fixed point. These results show the discrete model has more richer dynamics compared with the continuous model. The computations of the largest Lyapunov exponents more than 0 confirm the chaotic behaviors of the system , . Specifically, the chaotic orbits at an unstable fixed point are stabilized by using the feedback control method.

#### 1. Introduction

Within the field of epidemic theory, both the continuous-time model and the discrete-time model were described by the different equations. The continuous-time epidemic models have been widely investigated by many researchers (e.g., [1–10]). A Leslie-Gower predator-prey model with disease in prey incorporating a prey refuge was investigated in [7]. A ratio-dependent predator-prey model with Allee effect and disease in prey was investigated in [8]. Meanwhile, plenty of scholars used some discretization scheme to the continuous-time systems to study the consistency, convergence, permanence, and stability of the discrete system in [11–16]. In the recent years, there appear to be a number of articles about the flip bifurcation, fold bifurcation, and pitchfork bifurcation and in the sense of Marottos chaos of the discrete system which were presented in [17–28]. Both a discrete-time predator-prey system without Allee effect and a discrete-time predator-prey system with Allee effect were studied in [22, 23]. Dynamics of a system of three interacting populations with Allee effects and stocking was researched in [26]. Stability analysis and rich oscillation patterns in discrete-time FitzHugh-Nagumo excitable system with delayed coupling were studied in [24].

In this paper by using the center manifold theorem and bifurcation theorem, the local stability and bifurcation behaviors of a simple epidemic model with surprising dynamic is investigated in detail. The paper is organized as follows. Section 2 gives the specific conditions for existence of flip bifurcation and Neimark-Sacker bifurcation. Section 3 analyzes the two bifurcations by choosing the time step size as the parameter of the bifurcation. Section 4 verifies the theoretical analysis and displays the complex and surprising dynamics. In Section 5, chaos approaching to the unable fixed point is controlled by using the feedback control method.

A simple epidemic model with surprising dynamic was showed in [1] aswhere is the total population size; the birth process incorporates density dependent effects through a logistic equation with intrinsic growth rate and the carrying capacity , and other parameters are positive constant.

It is obvious that system (1) is a two-dimensional nonlinear ODE, where it is an impossible task to obtain its analytical solutions to study its qualitative properties including the global stability, so it is necessary to solve numerical solutions or approximate solutions of system (1) according to different discrete difference schemes. Therefore, in order to study the dynamics behavior of system (1), we will focus on the complex dynamical behaviors of a simple discrete epidemic model approached by the forward Euler scheme.

In this paper, we apply the forward Euler scheme to system (1) to get the following form:where is the step size.

#### 2. Stability of Fixed Points

The fixed points of system (2) are and , where satisfyThe Jacobian matrix of (2) at any point is written asThe characteristic equation of Jacobian matrix (4) is given bywhereLet us give a lemma as similar as in [21].

Lemma 1. *Suppose and are two roots of ; the local stability of the positive point can be determined by the modules of and *(i)*If and , then not only is a sink, but also it is locally asymptotically stable.*(ii)*If and , then not only is a source, but also it is unstable.*(iii)*If and and , then is a saddle.*(iv)*If or , then is nonhyperbolic.*

*In order to discuss the stability of the fixed point of system (2), the following lemma is also needed which can be easily proved by the relation between roots and coefficient of a quadratic equation [21].*

*Lemma 2. Let . Suppose that ; and are root of . Then(i) and if and ;(ii) and and if and only if ;(iii) and if and only if and ;(iv) and if and only if and ;(v) or and and are complex conjugate if and only if and .*

*Proposition 3. The eigenvalues of the fixed point are , .(i) is sink if and .(ii) is source if and .(iii) is nonhyperbolic if and or .(iv) is a saddle, except for that values of parameters which (i) to (iii).LetIf term (iii) of Proposition 3 holds, then one of the eigenvalues of the fixed point is −1 and the magnitude of the other is not equal to 1. The point undergoes flip bifurcation when the parameter changes in small neighborhood of .*

*The characteristic equation of the Jacobian matrix of system (2) at the positive is written aswhereLetTherefore,According to Lemma 2, we get the following proposition.*

*Proposition 4. There exists the fixed point for every possible parameter:(i) is a sink if satisfying(i.1) and ;(i.2) and ;(ii) is a source if satisfying(ii.1) and ;(ii.2) and ;(iii) is nonhyperbolic if satisfying(iii.1) and ;(iii.2) and ;(iv) is a saddle, except for that values of parameters which (i) to (iii). If (iii.1) holds, then and .LetIf (iii.2) holds, then and are a pair of complex conjugate numbers and and .*

Let

*3. Bifurcation Analysis on the Positive *

*3. Bifurcation Analysis on the Positive*

*In this section, we will investigate the flip bifurcation and Neimark-Sacker bifurcation corresponding to the positive fixed point.*

*3.1. Flip Bifurcation*

*3.1. Flip Bifurcation*

*Consider system (2) with arbitrary parameter , thenwhere*

*Now, give a perturbation of model (14) as follows:where is a limited perturbation parameter.*

*Let , . After the transformation of the fixed to , as follows:whereand .*

*Let*

*Consider the following translation:*

*Taking on both sides of (17), obtainwhere*

*Apply the center manifold theorem to (17) at the origin in the limited neighborhood of . The center manifold can be approximately presented asfor , sufficiently small.*

*We assume a center manifold of the formwhich must satisfy*

*By the simple calculations, obtain*

*Consider the map restricted to the center manifold :where*

*According to flip bifurcation, we get*

*After simple calculations, we obtain*

*Theorem 5. If , system (2) will undergo flip bifurcation at the positive point . Also the period-2 orbits that bifurcate from fixed are stable (resp., unstable) if (resp., ).*

*3.2. Neimark-Sacker Bifurcation*

*3.2. Neimark-Sacker Bifurcation*

*Consider model (2) with arbitrary parameter , which is described by is the fixed point of model (2), which is given by (3) and .*

*We consider the perturbation of model (31) as follows:where is a limited perturbation parameter.*

*Let , . After the transformation of the fixed to , as follows:whereand .*

*The characteristic equation of model (31) at is as follows:*

*Since the parameter , the characteristic values of are a pair of complex conjugate numbers and with modulus 1:*

*Now, we have*

*When varies in limited neighborhood of , let , . Neimark-Sacker bifurcation requires that when , then (), which is equivalent to . Since and , therefore ; then , which follows that*

*Let*

*Consider the translation below:*

*Taking on both sides of (19), obtainwhere*

*Now*

*According to Neimark-Sacker bifurcation, the discriminatory quantities are given by where*

*Then we can conclude the following theorem.*

*Theorem 6. If condition (39) holds, , when the parameter alters around the limit region of point , system (2) undergoes a Neimark-Sacker bifurcation at the point . If , an attracting invariant closed curve bifurcates from the point for .*

*4. Numerical Simulation*

*4. Numerical Simulation*

*In this section, we will give the numerical simulations for model (2), which not only illustrate our theoretical analysis but also present the complex, rich, and surprising dynamical behavior. The bifurcation parameters are considered for the following two cases:(1)varying in the range and fixing , , , , , , and ;(2)varying in the range and fixing , , , , , , and .*

*Case 1. *For Case 1 , , , , , , and with fixed point at the initial value point . By a simple calculation, we get , , , which established Theorem 5. From bifurcation diagram (see Figure 1(c)) and the largest Lyapunov exponent of model (2) (see Figure 1(a)), the fixed point of system (2) is stable when and unstable with and there exists flip bifurcation when . Largest Lyapunov exponents (Figure 1(a)) correspond to Figure 1(c).