Abstract

Based on the stability and orthogonal polynomial approximation theory, the ordinary, dislocated, enhancing, and random feedback control methods are used to suppress the Neimark-Sacker bifurcation to fixed point in this paper. It is shown that the convergence rate of enhancing feedback control and random feedback control can be faster than those of dislocated and ordinary feedback control. The random feedback control method, which does not require any adjustable control parameters of the model, just only slightly changes the random intensity. Finally, numerical simulations are presented to verify the effectiveness of the proposed controllers.

1. Introduction

The studies of biological models gradually become one of hot spots in nonlinear dynamics. The biological models have great research background and actual significance; therefore, a growing number of researchers have shown great interests in the research of biological models. In many biological models and practical problems, bifurcation and chaos are undesirable behaviors. Thus, we need to control them. In 1976, a population model, is given by Ecologist May for the first time. A one-dimensional deterministic delayed population model, is investigated by Sun et al. [1].

Recently, the Hopf bifurcation has been given much attention, and those works about bifurcation mainly include the validated existence of bifurcation and its control [25]. The aim of bifurcation control is to design a controller to modify the bifurcation properties of a given nonlinear system and then achieve the other desirable dynamical behaviors. OGY feedback control method is studied by Ott et al. [6]. Chen et al. have investigated the feedback control in continuous-time systems [7, 8]. The control of Hopf bifurcation in time-delayed neural network system is investigated by Zhou et al. [9]. Bifurcation analysis and tracking control of an epidemic model with nonlinear incidence rate are investigated by Yi et al. [10]. Wen and Xu studied feedback control of Hopf–Hopf interaction bifurcation with development of torus solutions in high-dimensional maps [11]. Feedback control of bifurcation and chaos in dynamical systems is investigated by Abed and Wang [12]. The Hopf bifurcation control via dynamic state-feedback control is studied by Nguyen and Hong in [13]. Amplitude control of limit cycle from Hopf bifurcation is studied in [14, 15]. Hopf bifurcation control of the system based on washout filter controller is investigated by Wu and Sun [16]. Liu and Xiao have studied complex dynamic behaviors of a discrete-time predator-prey system [17].

However, owing to the uncertain factors of external environment, manufacture, material, and installation, some parameters in practical model are not constant and will be characterized as bound random parameters [18]. The stochastic system can accurately represent the original system better. Therefore the study of stochastic system is more meaningful than deterministic systems. The Hopf bifurcation control is investigated in stochastic system with random parameter [1820]. It is of interest to examine the stochastic method in biological system and explore its implications.

The rest of this letter is organized as follows. In Section 2, the conditions for the emergence of Neimark-Sacker bifurcation are reviewed. In Section 3, the ordinary, dislocated, enhancing, and random feedback controls for controlling Neimark-Sacker bifurcation are proposed. And numerical simulations are presented to verify the effectiveness of the proposed bifurcation control methods. Finally, conclusions are given in Section 4.

2. Neimark-Sacker Bifurcation

Let us consider the logistic population model [1] for a single species:where stands for the population size at time and is the growth rate. In the real environment, the population size is determined not only by the current population size but also by its size in the past. So, we consider where stands for the population size at time and is the growth rate. If we introduce in model (4), a two-dimensional discrete-time dynamical model [2] can be rewritten as

By a simple computation with mathematical software, it is straightforward to obtain the following proposition.

Proposition 1. (a) For all parameter values, model (5) has one fixed point, .
(b) If , then model (5) has, additionally, a nontrivial positive fixed point, , where .
The Jacobian matrix of model (5) evaluated at the fixed point is given byand the characteristic equation of Jacobian matrix of model (5) can be written aswhere , .

Next, according to the point of view of biology, we study the stability of the nonzero fixed points. Note that the local stability of a fixed point is determined by the modules of eigenvalues of the characteristic equation at the fixed point. From the mathematical software and Lemma 2.2 [17], the following proposition shows the local stability of the fixed point .

Proposition 2. (a) is a sink if .
(b) is a source if .
(c) is not hyperbolic if .
When the term (c) of Proposition 2 holds, we can obtain that the eigenvalues of the matrix at the fixed point are a pair of conjugate complex numbers, the modules of which are one. The condition in term (c) of Proposition 2 can be written as the setthe fixed point can undergo Neimark-Sacker bifurcation when parameters vary in the small neighborhood of . By a simple computation, all eigenvalues of (7) arewhen , we can obtain eigenvaluesObviously, the transversality condition, the nondegeneracy condition, and the additional nondegeneracy condition of Neimark-Sacker bifurcation hold (see [3]). Thus, the nontrivial fixed point loses stability in the small neighborhood of . The bifurcation diagram and phase portrait for model (5) are depicted in Figure 1.

3. Neimark-Sacker Bifurcation Control

Without control, model (5) undergoes Neimark-Sacker bifurcations at the point corresponding to the value of bifurcation parameter as , as shown in Figure 1. To control the Neimark-Sacker bifurcation to the fixed point, the ordinary, dislocated, enhancing, and random feedback control methods are introduced as shown in the following.

3.1. Ordinary Feedback Control

For the ordinary feedback control, the system variable is often multiplied by a coefficient as the feedback gain, and the feedback gain is added to the right-hand side of the corresponding equation. Let the feedback control input , and the controlled model is given bywhere is the feedback coefficient.

Theorem 3. The necessary and sufficient condition for the controlled delayed species model (11) to be asymptotically stable at fixed point is , , and .

Proof . The Jacobi matrix of model (11) isand the characteristic equation of Jacobi matrix isaccording to Lemma 2.2 [17], the eigenvalues lie inside unit circle if and only if , , and . By a simple computation, we can obtain the following conditions: (a) , (b) , and (c) . Thus, when (a), (b), and (c) hold, the controlled delayed species model (11) will gradually converge to the fixed point.
Numerical simulations are used to investigate the controlled delayed species model (11). From Theorem 3 we conclude that our model (11) will gradually converge to the point for , when . The feedback coefficient is given by . The initial values in model (11) are taken as . The behaviors of the states of the controlled delayed species model (11) with time are displayed in Figure 2, respectively.

3.2. Dislocated Feedback Control

For the dislocated feedback control, a system variable multiplied by a coefficient is added to the right-hand side of another equation. Then, this method is called dislocated feedback control. Let feedback control input , and the controlled model is given bywhere is the feedback coefficient.

Theorem 4. The necessary and sufficient condition for the controlled delayed species model (14) to be asymptotically stable at fixed point is , , and .

Proof . The Jacobi matrix of model (14) isand the characteristic equation of Jacobi matrix isaccording to Lemma 2.2 [17], the eigenvalues lie inside unit circle if and only if , , and . By a simple computation, we can obtain the following conditions: (a) ; (b) ; (c) . Thus, when (a), (b), and (c) hold, the controlled delayed species model (14) will gradually converge to the fixed point.
Numerical simulations are used to investigate the controlled delayed species model (14). From Theorem 3, we conclude that our model (14) will gradually converge to the point for , when . The feedback coefficient is given by . The initial values in model (14) are taken as . The behaviors of the states of the controlled delayed species model (14) with time are displayed in Figure 3, respectively.

3.3. Enhancing Feedback Control

For the enhancing feedback control, it is difficult for a complex system to be controlled by only one feedback variable, and in such cases the feedback gain is always very large. So we consider using multiple variables multiplied by a coefficient as the feedback gain. This method is called enhancing feedback control. Let feedback control inputs , , and the controlled model is given bywhere is the feedback coefficient.

Theorem 5. The necessary and sufficient condition for the controlled delayed species model (17) to be asymptotically stable at equilibrium point is , , and .

Proof. The Jacobi matrix of system (17) isand the characteristic equation of Jacobi matrix isaccording to Lemma 2.2 [17], the eigenvalues lie inside unit circle if and only if , , and . By a simple computation, we can obtain the following conditions: (a) ; (b) ; (c) . Thus, when (a), (b), and (c) hold, the controlled delayed species model (17) will gradually converge to the fixed point.
Numerical simulations are used to investigate the controlled delayed species model (17). From Theorem 3 we conclude that our model (17) will gradually converge to the point for , when . The feedback coefficient is given by . The initial values in model (17) are taken as . The behaviors of the states of the controlled delayed species model (17) with time are displayed in Figure 4, respectively.

3.4. Random Feedback Control

To achieve the control objectives, we need to adjust the control gains for the feedback control methods above. Thus, we use the random feedback control method to control them. If a system variable multiplied by a random coefficient is added to the right-hand side of equation, then this method is called random feedback control. Let , ; after applying the coordinate transformation, the fixed point is converted to the origin ; then we have the following system:Let feedback control input ; is a random parameter; . Taking this controller into the right side of the second equation in (20), the controlled delayed species model can be written aswhere is the statistic parameter of , and are the input direction of state variable and random variable in the controller, is random intensity, and is the random variable defined on nonnegative set integer with the probability density function . According to the orthogonal polynomial approximation [1820] of discrete random function in the Hilbert space, the response of (21) can be expressed by the following series: where is the th Charlier orthogonal polynomial and is the largest order of the polynomials we have taken. Substituting (22) into (21), we haveUsing the orthogonal polynomial approximation theory, the nonlinear term in the right side of (23) can be further reduced into a linear combination of related single polynomials. It is written asWe obtainedThe recurrence formula for Charlier polynomial isIn order to facilitate the numerical analysis of this paper, we take , . Based on the orthogonality of orthogonal polynomials, we can finally obtain the equivalent deterministic equation:And the ensemble mean response of model (21) isHere we will first discuss the influence of random feedback controller on the Neimark-Sacker bifurcation. Next we will analyze the Neimark-Sacker bifurcation control with random feedback controller based on deterministic controlled model (27). We suppose that the random intensity of the random controller is from 0 to 0.1. In this section, let , , , and , respectively. Owing to the small quantities of intensity , the same initial conditions for model (5) and model (27) are given; namely, , and , . The time history diagrams of the ensemble mean response of controlled delayed species model are shown in Figure 5. By using the same strategy we can discuss the Neimark-Sacker bifurcation control for , and , cases (not reported here).

For the random feedback control, which does not require any adjustable feedback control parameters, we need to slightly change the random intensity of random feedback controller as the random intensity is very small. By means of numerical simulations, we find that the random feedback method to control the Neimark-Sacker bifurcation is available. Next we discuss the influence of the initial values for random feedback control. Taking , as random initial values covering , according to numerical simulations, we can find that random feedback control has robustness for the random initial values. The time history diagrams of the ensemble mean response of controlled delayed species model (21) are shown in Figure 6.

4. Conclusions

In this paper, the ordinary, dislocated, enhancing, and random feedback control strategies are studied for controlling the Neimark-Sacker bifurcation in the delayed species model. It is found that the convergence rate of variables of enhancing feedback control and random feedback control can be faster than those of dislocated feedback control and ordinary feedback control. From Figure 5 we can find that the effect of the random feedback controller on controlling the Neimark-Sacker bifurcation is available. From Figure 6 we can find that the random feedback controller has robustness for random initial values. Furthermore, numerical simulations are presented to verify the effectiveness of the proposed controllers. Those methods proposed in this paper can be extended to consider other dynamical behaviors such as the double Neimark-Sacker bifurcation and chaos.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the Guizhou Science and Technology Cooperation Plan (Qian Ke He LH zi 7032); National Natural Science Foundation of China (Grants nos. 71461027 and 71471158); Science and Technology Talent Training Object of Guizhou Province Outstanding Youth (Qian ke he ren zi 06); Guizhou Province Natural Science Foundation in China (Qian Jiao He KY 295); 2013, 2014, and 2015 Zunyi 15851 Talents Elite Project Funding; Zunyi Innovative Talent Team (Zunyi KH (2015)38); Project of Teaching Quality and Teaching Reform of Higher Education in Guizhou Province (Qian Jiao gaofa 446,337); Zunyi Innovative Talent Team (Zunyi KH (2015)38); Project of Teaching Quality and Teaching Reform of Higher Education in Guizhou Province (Qian Jiao gaofa 446,337); Guizhou Science and Technology Cooperation Plan (Qian Ke He LH zi 7028); Guizhou Science and Technology Cooperation Plan (Qian Ke He LH zi 7047); and Doctoral Scientific Research Foundation of Zunyi Normal College (BS09).