# Controlling Neimark-Sacker Bifurcation in Delayed Species Model Using Feedback Controller

Authors: Jie Ran; Yanmin Liu; Jun He; Xiang Li
Journal: Advances in Mathematical Physics
 (2016)

---

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


---

## Body

## 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,(1)xn+1=rxn1-xn,r∈0,4,x∈0,1,is given by Ecologist May for the first time. A one-dimensional deterministic delayed population model, (2)xn+1=rxn1-xn-1,r∈0,2.28,x∈0,1,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 [2–5]. 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 [18–20]. 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 Section2, 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:(3)xn+1=rxn1-xn,where xn stands for the population size at time n and r 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 (4)xn+1=rxn1-xn-1,where xn-1 stands for the population size at time n-1 and r is the growth rate. If we introduce yn=xn-1 in model (4), a two-dimensional discrete-time dynamical model [2] can be rewritten as(5)xn+1=rxn1-yn,yn+1=xn.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, O0,0.
(b) Ifr>1, then model (5) has, additionally, a nontrivial positive fixed point, O∗x∗,y∗, where x∗=y∗=1-1/r.
The Jacobian matrix of model (5) evaluated at the fixed point x,y is given by(6)J=r-ry-rxx0,and the characteristic equation of Jacobian matrix of model (5) can be written as(7)fλ=λ2+px,yλ+qx,y=0,where px,y=ry-r, qx,y=-rx.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 O∗x∗,y∗.Proposition 2.
(a)O∗x∗,y∗ is a sink if 1<r<2.
(b)O∗x∗,y∗ is a source if r>2.
(c)O∗x∗,y∗ is not hyperbolic if r=2.
When the term (c) of Proposition2 holds, we can obtain that the eigenvalues of the matrix J at the fixed point O∗x∗,y∗ 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 set(8)Hc=r:r=2;the fixed point O∗x∗,y∗ can undergo Neimark-Sacker bifurcation when parameters vary in the small neighborhood of Hc. By a simple computation, all eigenvalues of (7) are(9)λ1,2=12±54-r;when r=rc=2, we can obtain eigenvalues(10)λ1,2=12±32i.
Obviously, 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 Hc. The bifurcation diagram and phase portrait for model (5) are depicted in Figure 1.Figure 1
(a) Bifurcation diagram of model (5) with bifurcation parameter r covering [1.8,2.25] and (b) phase portrait corresponding to bifurcation parameter r=2.01. The initial values are x0=0.2,y0=0.2.
(a) (b)


## 3. Neimark-Sacker Bifurcation Control

Without control, model (5) undergoes Neimark-Sacker bifurcations at the point 0.5025,0.5025 corresponding to the value of bifurcation parameter as r=2.01, 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 inputu=k1yn-y∗, and the controlled model is given by(11)xn+1=rxn1-yn,yn+1=xn+k1yn-y∗,where k1 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 r-1>0, 2k1+r+1>0, and k1+r-2<0.Proof .
The Jacobi matrix of model (11) is(12)J=1-r+11k1,and the characteristic equation of Jacobi matrix J is(13)fλ=λ2-k1+1λ+k1+r-1=0;according to Lemma 2.2 [17], the eigenvalues lie inside unit circle if and only if f1>0, f-1>0, and qx∗,y∗<1. By a simple computation, we can obtain the following conditions: (a) r-1>0, (b) 2k1+r+1>0, and (c) k1+r-2<0. 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 0.5025,0.5025 for k1∈-1.505,-0.01, when r=2.01. The feedback coefficient is given by k1=-0.05. The initial values in model (11) are taken as x0=0.2,y0=0.2. The behaviors of the states xn,yn of the controlled delayed species model (11) with time n are displayed in Figure 2, respectively.Figure 2
Time history diagram of the controlled delayed species model (11) for k1=-0.05 and r=2.01; the initial values are x0=0.2,y0=0.2.
(a) (b)


### 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 inputu=k2xn-x∗, and the controlled model is given by(14)xn+1=rxn1-yn,yn+1=xn+k2xn-x∗,where k2 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 r+rk2-k2-1>0, rk2+r-k2+1>0, and rk2+r-k2-2<0.Proof .
The Jacobi matrix of model (14) is(15)J=1-r+11+k20,and the characteristic equation of Jacobi matrix J is(16)fλ=λ2-λ+rk2+r-k2-1=0;according to Lemma 2.2 [17], the eigenvalues lie inside unit circle if and only if f1>0, f-1>0, and qx∗,y∗<1. By a simple computation, we can obtain the following conditions: (a) r+rk2-k2-1>0; (b) rk2+r-k2+1>0; (c) rk2+r-k2-2<0. 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 0.5025,0.5025 for k2∈-1,-0.00990099, when r=2.01. The feedback coefficient is given by k2=-0.03. The initial values in model (14) are taken as x0=0.2,y0=0.2. The behaviors of the states xn,yn of the controlled delayed species model (14) with time n are displayed in Figure 3, respectively.Figure 3
Time history diagram of the controlled delayed species model (14) for k2=-0.03 and r=2.01. The initial values are x0=0.2,y0=0.2.
(a) (b)


### 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 inputsu1=k3xn-x∗, u2=k3yn-y∗, and the controlled model is given by(17)xn+1=rxn1-yn+k3xn-x∗,yn+1=xn+k3yn-y∗,where k3 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 -1-k3+k32+r+(1/r)k3-k32>0, 1+3k3+k32+r+(1/r)k3+k32>0, and k3+k32+r-(1/r)k32-2<0.Proof.
The Jacobi matrix of system (17) is(18)J=1+k31-1r-r+11k3,and the characteristic equation of Jacobi matrix J is(19)fλ=λ2+-1-2k3+k3rλ+k3+k32+r-k32r-1=0;according to Lemma 2.2 [17], the eigenvalues lie inside unit circle if and only if f1>0, f-1>0, and qx∗,y∗<1. By a simple computation, we can obtain the following conditions: (a) -1-k3+k32+r+(1/r)k3-k32>0; (b) 1+3k3+k32+r+(1/r)k3+k32>0; (c) k3+k32+r-(1/r)k32-2<0. 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 0.5025,0.5025 for k3∈-1.980048,-0.010050, when r=2.01. The feedback coefficient is given by k3=-0.41. The initial values in model (17) are taken as x0=0.2,y0=0.2. The behaviors of the states xn,yn of the controlled delayed species model (17) with time n are displayed in Figure 4, respectively.Figure 4
Time history diagram of the controlled delayed species model (17) for k3=-0.41 and r=2.01. The initial values are x0=0.2,y0=0.2.
(a) (b)


### 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. Letun=xn-x∗, vn=yn-y∗; after applying the coordinate transformation, the fixed point is converted to the origin O0,0; then we have the following system:(20)un+1=un-runvn-r-1un+vn,vn+1=un.Let feedback control input u=kun+vn; k is a random parameter; k=αk4+βδξ. Taking this controller into the right side of the second equation in (20), the controlled delayed species model can be written as(21)un+1=un-runvn-r-1un+vn,vn+1=un+αk4+βδξun+αk4+βδξvn,where k4 is the statistic parameter of k, α 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 pξ. According to the orthogonal polynomial approximation [18–20] of discrete random function in the Hilbert space, the response of (21) can be expressed by the following series: (22)un,ξ=∑i=0MuinCiξ,vn,ξ=∑i=0MvinCiξ,uin=∑k=0Npξun,ξCiξ,vin=∑k=0Npξvn,ξCiξ,where Cik is the ith Charlier orthogonal polynomial and M is the largest order of the polynomials we have taken. Substituting (22) into (21), we have(23)∑i=0Muin+1Ciξ=∑i=0MuinCiξ-r∑i=0MuinCiξ∑i=0MvinCiξ-r-1∑i=0MuinCiξ+∑i=0MvinCiξ,∑i=0Mvin+1Ciξ=∑i=0MuinCiξ+αk4+βδξ∑i=0MuinCiξ+αk4+βδξ∑i=0MvinCiξ.Using 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 as(24)∑i=0MuinCik∑i=0MvinCik=∑i=02MXinCik.We obtained(25)∑i=0Muin+1Ciξ=∑i=0MuinCiξ-r∑i=02MXinCiξ-r-1∑i=0MuinCiξ+∑i=0MvinCiξ,∑i=0Mvin+1Ciξ=∑i=0MuinCiξ+αk4+βδξ∑i=0MuinCiξ+αk4+βδξ∑i=0MvinCiξ.The recurrence formula for Charlier polynomial is(26)ξCiξ=Ci+1ξ+i+λCiξ+iλCi-1ξ.In order to facilitate the numerical analysis of this paper, we take M=1, λ=1. Based on the orthogonality of orthogonal polynomials, we can finally obtain the equivalent deterministic equation:(27)u0n+1=u0n-ru0nv0n+u1nv1n-r-1u0n+v0n,v0n+1=u0n+αk4u0n+αk4v0n+βδu0n+βδv0n,u1n+1=u1n-ru0nv1n+u1nv0n+u1nv1n-r-1u1n+v1n,v1n+1=u1n+αk4v1n+αk4u1n+βδ2u1n+u0n+βδ2v1n+v0n.And the ensemble mean response of model (21) is(28)Eun,k=∑i=01uinEPik=u0n,Evn,k=∑i=01vinEPik=v0n.Here 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 α=-1, β=1, k4=0.1, and δ=0.001, respectively. Owing to the small quantities of intensity δ, the same initial conditions for model (5) and model (27) are given; namely, u0=u0=0.2, v0=v0=0.2 and U0=0.2,0T, V0=0.2,0T. 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 α=1, β=±1 and α=-1, β=-1 cases (not reported here).Figure 5
Time history diagrams of the controlled delayed species model (21) for δ=0.001, k4=0.1. The initial values are [u(0)=0.2,v(0)=0.2].
(a) (b)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. Takingu0, v0 as random initial values covering 0,0.5, 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.Figure 6
Time history diagrams of the controlled delayed species model (21) for δ=0.001, k4=0.1. The initial values u(0) and v(0) take random values covering [0,0.5].
(a) (b)


## 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 inputu=k1yn-y∗, and the controlled model is given by(11)xn+1=rxn1-yn,yn+1=xn+k1yn-y∗,where k1 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 r-1>0, 2k1+r+1>0, and k1+r-2<0.Proof .
The Jacobi matrix of model (11) is(12)J=1-r+11k1,and the characteristic equation of Jacobi matrix J is(13)fλ=λ2-k1+1λ+k1+r-1=0;according to Lemma 2.2 [17], the eigenvalues lie inside unit circle if and only if f1>0, f-1>0, and qx∗,y∗<1. By a simple computation, we can obtain the following conditions: (a) r-1>0, (b) 2k1+r+1>0, and (c) k1+r-2<0. 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 0.5025,0.5025 for k1∈-1.505,-0.01, when r=2.01. The feedback coefficient is given by k1=-0.05. The initial values in model (11) are taken as x0=0.2,y0=0.2. The behaviors of the states xn,yn of the controlled delayed species model (11) with time n are displayed in Figure 2, respectively.Figure 2
Time history diagram of the controlled delayed species model (11) for k1=-0.05 and r=2.01; the initial values are x0=0.2,y0=0.2.
(a) (b)


## 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 inputu=k2xn-x∗, and the controlled model is given by(14)xn+1=rxn1-yn,yn+1=xn+k2xn-x∗,where k2 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 r+rk2-k2-1>0, rk2+r-k2+1>0, and rk2+r-k2-2<0.Proof .
The Jacobi matrix of model (14) is(15)J=1-r+11+k20,and the characteristic equation of Jacobi matrix J is(16)fλ=λ2-λ+rk2+r-k2-1=0;according to Lemma 2.2 [17], the eigenvalues lie inside unit circle if and only if f1>0, f-1>0, and qx∗,y∗<1. By a simple computation, we can obtain the following conditions: (a) r+rk2-k2-1>0; (b) rk2+r-k2+1>0; (c) rk2+r-k2-2<0. 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 0.5025,0.5025 for k2∈-1,-0.00990099, when r=2.01. The feedback coefficient is given by k2=-0.03. The initial values in model (14) are taken as x0=0.2,y0=0.2. The behaviors of the states xn,yn of the controlled delayed species model (14) with time n are displayed in Figure 3, respectively.Figure 3
Time history diagram of the controlled delayed species model (14) for k2=-0.03 and r=2.01. The initial values are x0=0.2,y0=0.2.
(a) (b)


## 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 inputsu1=k3xn-x∗, u2=k3yn-y∗, and the controlled model is given by(17)xn+1=rxn1-yn+k3xn-x∗,yn+1=xn+k3yn-y∗,where k3 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 -1-k3+k32+r+(1/r)k3-k32>0, 1+3k3+k32+r+(1/r)k3+k32>0, and k3+k32+r-(1/r)k32-2<0.Proof.
The Jacobi matrix of system (17) is(18)J=1+k31-1r-r+11k3,and the characteristic equation of Jacobi matrix J is(19)fλ=λ2+-1-2k3+k3rλ+k3+k32+r-k32r-1=0;according to Lemma 2.2 [17], the eigenvalues lie inside unit circle if and only if f1>0, f-1>0, and qx∗,y∗<1. By a simple computation, we can obtain the following conditions: (a) -1-k3+k32+r+(1/r)k3-k32>0; (b) 1+3k3+k32+r+(1/r)k3+k32>0; (c) k3+k32+r-(1/r)k32-2<0. 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 0.5025,0.5025 for k3∈-1.980048,-0.010050, when r=2.01. The feedback coefficient is given by k3=-0.41. The initial values in model (17) are taken as x0=0.2,y0=0.2. The behaviors of the states xn,yn of the controlled delayed species model (17) with time n are displayed in Figure 4, respectively.Figure 4
Time history diagram of the controlled delayed species model (17) for k3=-0.41 and r=2.01. The initial values are x0=0.2,y0=0.2.
(a) (b)


## 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. Letun=xn-x∗, vn=yn-y∗; after applying the coordinate transformation, the fixed point is converted to the origin O0,0; then we have the following system:(20)un+1=un-runvn-r-1un+vn,vn+1=un.Let feedback control input u=kun+vn; k is a random parameter; k=αk4+βδξ. Taking this controller into the right side of the second equation in (20), the controlled delayed species model can be written as(21)un+1=un-runvn-r-1un+vn,vn+1=un+αk4+βδξun+αk4+βδξvn,where k4 is the statistic parameter of k, α 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 pξ. According to the orthogonal polynomial approximation [18–20] of discrete random function in the Hilbert space, the response of (21) can be expressed by the following series: (22)un,ξ=∑i=0MuinCiξ,vn,ξ=∑i=0MvinCiξ,uin=∑k=0Npξun,ξCiξ,vin=∑k=0Npξvn,ξCiξ,where Cik is the ith Charlier orthogonal polynomial and M is the largest order of the polynomials we have taken. Substituting (22) into (21), we have(23)∑i=0Muin+1Ciξ=∑i=0MuinCiξ-r∑i=0MuinCiξ∑i=0MvinCiξ-r-1∑i=0MuinCiξ+∑i=0MvinCiξ,∑i=0Mvin+1Ciξ=∑i=0MuinCiξ+αk4+βδξ∑i=0MuinCiξ+αk4+βδξ∑i=0MvinCiξ.Using 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 as(24)∑i=0MuinCik∑i=0MvinCik=∑i=02MXinCik.We obtained(25)∑i=0Muin+1Ciξ=∑i=0MuinCiξ-r∑i=02MXinCiξ-r-1∑i=0MuinCiξ+∑i=0MvinCiξ,∑i=0Mvin+1Ciξ=∑i=0MuinCiξ+αk4+βδξ∑i=0MuinCiξ+αk4+βδξ∑i=0MvinCiξ.The recurrence formula for Charlier polynomial is(26)ξCiξ=Ci+1ξ+i+λCiξ+iλCi-1ξ.In order to facilitate the numerical analysis of this paper, we take M=1, λ=1. Based on the orthogonality of orthogonal polynomials, we can finally obtain the equivalent deterministic equation:(27)u0n+1=u0n-ru0nv0n+u1nv1n-r-1u0n+v0n,v0n+1=u0n+αk4u0n+αk4v0n+βδu0n+βδv0n,u1n+1=u1n-ru0nv1n+u1nv0n+u1nv1n-r-1u1n+v1n,v1n+1=u1n+αk4v1n+αk4u1n+βδ2u1n+u0n+βδ2v1n+v0n.And the ensemble mean response of model (21) is(28)Eun,k=∑i=01uinEPik=u0n,Evn,k=∑i=01vinEPik=v0n.Here 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 α=-1, β=1, k4=0.1, and δ=0.001, respectively. Owing to the small quantities of intensity δ, the same initial conditions for model (5) and model (27) are given; namely, u0=u0=0.2, v0=v0=0.2 and U0=0.2,0T, V0=0.2,0T. 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 α=1, β=±1 and α=-1, β=-1 cases (not reported here).Figure 5
Time history diagrams of the controlled delayed species model (21) for δ=0.001, k4=0.1. The initial values are [u(0)=0.2,v(0)=0.2].
(a) (b)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. Takingu0, v0 as random initial values covering 0,0.5, 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.Figure 6
Time history diagrams of the controlled delayed species model (21) for δ=0.001, k4=0.1. The initial values u(0) and v(0) take random values covering [0,0.5].
(a) (b)


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



---

*Source: 2028037-2016-09-04.xml*