Research Article | Open Access

# Bifurcation Analysis and Chaos Control in a Discrete Epidemic System

**Academic Editor:**Garyfalos Papashinopoulos

#### 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).**Let**If 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)*(iii)

*and ;*(ii.2)*and ;**is nonhyperbolic if satisfying(iii.1)*(iv)

*and ;*(iii.2)*and ;**is a saddle, except for that values of parameters which (i) to (iii). If (iii.1) holds, then and .*

*Let*

*If (iii.2) holds, then and are a pair of complex conjugate numbers and and .*

*Let*

#### 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

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

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

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

**(a)**

**(b)**

**(c)**

**(d)**

The phase portraits which are associated with Figures 1(a) and 1(c) are displayed in Figure 3. There are the orbits of periods 2, 4, and 8 with . Moreover, when , stable period windows and chaotic sets alternate in appearing. For example, there emerge the stable window of period-9, period-7, period-3, period-5, period-6, period-1 orbits. The phase portraits which are associated with Figures 1(b) and 1(d) are displayed in Figure 3. With increasing, chaos and period-doubling bifurcation take place by turn. Finally, the system enters into chaos.

*Case 2. *For Case 2 , , , , , , and . The initial value is , after a simple calculation for the fixed point , we get , , , , , , , (see Figure 2(b)); it established Theorem 6. From Figures 2(a) and 2(c), we observe that when , the fixed point of map (2) is stable, it loses its stability at , and there appears an invariant circle when exceeds 2.96. Figures 2(b) and 2(d) are local amplifications of Figures 2(a) and 2(c) respectively for .

**(a)**

**(b)**

**(c)**

**(d)**

The phase portraits which are associated with Figures 1(a) and 1(c) are displayed in Figure 4. Lyapunov exponents are smaller than 0 when ; that is to say, the nonchaotic region is smaller than the chaotic region . For , some Lyapunov exponents are larger than 0 and some Lyapunov exponents are smaller than 0, so there exists stable fixed point or stable period windows in the chaotic region. As we all know, Lyapunov exponents are bigger than 0 which implies the existence of chaos. With increases, circle disappears and period-53, period-34, period-19, period-23, period-19, period-42, period-15, and period-11 orbits appear, attracting chaotic sets (see Figure 4). From Figure 4, we can see that the attractive cycle is smooth in the vicinity of the bifurcation point. However, when decreases, the cycle inflates gradually and loses its smoothness, and the system enters chaos. This shows that not only the stability changes of system but also the limit cycle occur when system crosses the bifurcation point.

#### 5. Chaos Control

In this section, a feedback control law as the control force is added to (2): the following controlled form of model (2) iswhere and are the feedback gain and is the positive point of model (2). The Jacobian matrix of controlled system (46) is given bywhere are given in (17). The characteristic equation of (48) isLet and be the two roots of (49); thenThe lines of marginal stability are determined by solving the equation of and .

These conditions guarantee that the eigenvalues and have modulus less than 1. Assume that , and from (51) we have . Assume that , and from (50) and (51) we have Assume that , and from (50) and (51) we have The stable eigenvalues lie within a triangular region by lines , and (see Figure 5(a)).

**(a)**

**(b)**

**(c)**

We have performed some numerical simulations to see how the state feedback method controls the unstable fixed point. Parameter values are fixed as , ,