Neimark–Sacker Bifurcation and Controlling Chaos in a Three-Species Food Chain Model through the OGY Method

+e dynamics behavior of a discrete-time three-species food chain model is investigated. By using bifurcation theory, it is shown that the equilibrium point of the system loses its stability, and the system undergoes Neimark–Sacker bifurcation, which leads to chaos as the parameter changes.+e chaotic motion is controlled on the stable periodic period-1 orbit using the implementation of the hybrid control strategy.+e factor affecting the control time of chaos is also studied. Numerical simulations are consistent with the theoretical analysis. +e results of this research prove that the chaos control method can be extended to the higher-dimensional biological model and can be realized.


Introduction
Predation and prey behaviors are common phenomena in the ecosystem. Since the last century when Volterra and Lotka constructed the predator-prey model, the predatorprey model has been concerned by many scholars. Most scholars are deepening the original model to be more practical. ey consider many factors such as time delay, functional response, diffusion, and so on. ey also discuss the dynamic properties of the system, such as stability, persistence, and the existence of positive periodic solution, and get a lot of good results. Lotka-Volterra model describes interactions between two species in an ecosystem, a predator and a prey. e model was developed independently by Lotka [1] and Voltera [2]. Ali Khan et al. [3] propose a new harvesting strategy for controlling chaotic population in a three-species Hastings and Powell food chain model: where F i (U) � (A i U/(B i + U)), i � 1, 2. Steady state, limit cycle, and period-2 and period-4 behavior of the model are obtained in [3]. is threshold harvesting strategy will be very useful for species conservation and fishery management.
Panday et al. [4] propose a three-species food chain model, where the growth rate of middle predator is reduced due to the cost of fear of top predator, and the growth rate of prey is suppressed due to the cost of fear of middle predator: e global stability analysis of the equilibrium points and the chaotic behavior of the system (2) are investigated using Poincaré section and maximum Lyapunov exponent.
e Allee effect is an important phenomenon in population biology characterized by positive density dependence, that is, a positive correlation between population density and individual fitness. Parshad et al. [5] consider a ratio-dependent spatially explicit three-species food chain model, where the top predator is subjected to a strong Allee effect. ey show the existence of a global attractor for the following system: 3 .
Most of the research is based on continuous model. In reality, for some species with small number, short life cycle, and no cross between mother and son, it is more realistic to establish discrete population model to describe their growth process. Compared with the continuous models, the discrete-time food chain models can better analyze the dynamic behavior of the system. e bifurcation and other dynamic behaviors of the discrete system will be enriched with the change of parameters. e discrete model can make better use of chaos control strategies and computer simulation.
x n k − x n x n + w n − βx n w n , y n+1 � y n + h βθ e x n w n − cy n − θy n z n + w n , z n+1 � y n + h x n z n − y n z n − dz n , w n+1 � w n + h y n z n − (d + α)w n .
In system (4), it is shown that the system shows Neimark-Sacker bifurcation at a positive fixed point. e system shows chaotic dynamics as bifurcation parameter changes. e chaotic motion of the system is controlled through implementation of hybrid control method.
Ott et al. [13] proposed a new chaos control strategy. e method is not to use the existing dynamic control strategies or destroy the conditions for the occurrence of chaotic motion. e method is further improved by Romeiras et al. [14]. In 1992, the French physicist Pyragas proposed a timedelayed feedback control method for chaos [15]. In the following ten years, the research on chaos control and synchronization has been booming, which has rapidly become an important hotspot in the field of chaos research. For example, linear state feedback control [16], sliding mode control [17], and adaptive Lyapunov control [18] are widely used. Recently, the chaos control in biological systems has been studied. Some authors have researched that these chaos control methods can be applied to population dynamics and play nontrivial evolutionary roles in [7][8][9]19]. Other methods will change the fixed point of the control or change the original properties of the original system. OGY chaos control is a small real-time disturbance to the parameters to achieve the purpose of control and maintain the original characteristics of the original system. Feng [20,21] studied the dynamics and chaos control of two-dimensional Hassell recruitment population model and Ricker-type recruitment population model using OGY method. e research on chaos and chaos control has been greatly developed and widely used. A lot of valuable research results have been obtained in [22][23][24][25][26].
Using the Euler method, the continuous-time system is discretized. Moreover, the Poincaré map was used to discretize a continuous system, to study the Neimark-Sacker bifurcation, and also to control the system by the OGY approach in [27,28]. Recently, the design of an explicit analytical expression, the Poincaré map, and its use for the control of the chaotic system were achieved in [29][30][31].
In this paper, the chaotic dynamics of a three-species food chain model is investigated and the chaos is controlled. e dynamic behavior of this model is discussed by using the stability theory of nonlinear difference equation and bifurcation theory in [32]. F: x n+1 � f x n , y n , z n � ax n 1 − x n − bx n y n , y n+1 � g x n , y n , z n � cx n y n − by n z n , z n+1 � h x n , y n , z n � ry n z n .
is food chain model describes the insects group of three fully different insects. ese insects are a population of prey x, which is predated by a first predator with population y. A third species given by a top predator z that predates on the first predator y. Each of the species with nonoverlapping generations affects each other's population dynamics. Parameters a, b, c, d, and r are all positive. Parameter a is a reproduction rate of x; b is the predation rate of predator y to prey x and top predator z to y; c is the growth rate of 2 Discrete Dynamics in Nature and Society predators y; and r is the growth rate of predators z due to the consumption of species y. Elsadany [32] has studied the stability or bifurcation of the equilibrium, but has not studied the condition of Neimark-Sacker bifurcation or the control of chaos. In the paper, the conditions of Neimark-Sacker bifurcation are given and proved strictly. By studying chaos control of the model, we can better understand the dynamic behavior of the biological system. We can adjust the value of parameters through chaos control to make the ecological model in a state of equilibrium. Numerical simulation verifies the correctness of our theoretical derivation.
For the study of Neimark-Sacker bifurcation, most of the previous studies have used the central manifold theorem to obtain the bifurcation parameter values after tedious derivation. In the paper, a method is proposed to judge the Neimark-Sacker bifurcation ( eorem 1), which makes the operation and derivation simple and clear, and gives the bifurcation parameter values. e discrete-time three-species food chain model produces a chaotic attractor when a � 4.26, b �3.7, c �3, d � 3.5, and r � 3.8. e unstable point is controlled on the period-1 orbit. It shows the control times are different as the adjustment values to be selected are different. In biology or ecology, the complex chaotic behavior of this model shows the relationship of the different species, including the number population, reproduction rate, and survival rate, whether they can survive in a balanced state, or makes the population develop in disorder or chaos.

The Equilibrium Point and Local Stability Analysis of the Food Chain Model
In order to find out the equilibria of this system and study its dynamic properties. Equation (5) is written as (5) has four equilibria: is the axial fixed point in the absence of midlevel species y � 0 and top-level Assume E i (x * , y * , z * ), (i � 1, 2, 3, 4) is positive and denote equilibrium of equation (5). By the theorem in [20], the equilibria are classified according to the eigenvalues of linearized matrix of equation (5). To carry out linear stability analysis, the Taylor series expansion of equation (5) may be written as where Consider the matrix J � P 11 P 12 P 13 where Discrete Dynamics in Nature and Society e characteristic equation is det P 11 − λ P 12 P 13 which may be rewritten in the form where c 1 � − P 11 + P 22 + P 33 , c 2 � P 11 P 22 − P 12 P 21 + P 11 P 33 + P 22 P 33 − P 23 P 32 , c 3 � − P 11 P 22 P 33 + P 11 P 23 P 32 + P 12 P 21 P 33 .
(i) e eigenvalues of matrix (9) are λ 11 � a, λ 12 � λ 13 � 0 at E 1 , with eigenvectors, respectively, given by d 1 � (1,0,0), d 2 � (0,1,0), and d 3 � (0,0,3). When a < 1, the equilibrium point E 1 is a stable point; a > 1, the equilibrium E 1 is an unstable point; otherwise, when a � 1, E 1 is called nonhyperbolic point. (ii) e eigenvalues of matrix (9) 33 are real, and they are all inside the unit circle when e interior equilibrium will be locally asymptotically stable if the coefficients of the characteristic equation (12) Meanwhile, the equilibrium E 4 is a stable point, if and only if the coefficients c 1 , c 2 , and c 3 of equation (12) are satisfied:

Neimark-Sacker Bifurcation Analysis
Neimark-Sacker bifurcation is a critical value of a parameter where the system stability changes and periodic solution arises. In the next theorem, the existence of Neimark-Sacker bifurcation is given, and we take the reproduction rate parameter a of the prey x as a bifurcation parameter. In the second part, we have analyzed the local stability near the positive equilibrium in detail. Here, without losing generality, we give the Neimark-Sacker bifurcation conclusion around the equilibrium E 3 and prove it. Theorem 1. When the reproduction rate parameter a of the prey x crosses a critical value a 0 , then the system shows Neimark-Sacker bifurcation around the positive equilibrium E 3 if the following conditions hold: Proof. e position equilibrium E 3 is locally asymptotically stable and the system (6) will lose its stability when some parameter values are changed. Hence, we take the reproduction rate parameter a of the prey x as a bifurcation parameter. If there exists a critical value, a critical value a 0 such that When a � a 0 , the characteristic equation (14) can be of the form e above equation has three roots − c 1 (a 0 ), i ����� � c 2 (a 0 ) and − i ����� � c 2 (a 0 ). To show that Neimark-Sacker bifurcation occurs at a 0 , we need to satisfy the condition dRe(λ(a)) da For all a 0 , the roots are in general of the form Now, we will verify the condition Substituting λ k (a) � μ(a) ± i](a) into equation (18) Notice that μ(a 0 ) � 0, ](a 0 ) � ����� � λ 2 (a 0 ); then, we have Solving μ ′ (a 0 ) from equation (22), we get If (c 1 (a 0 )c 2 (a 0 )) ′ ≠ c 3 ′ (a 0 ), and c 3 (a 0 ) � − λ 1 (a 0 ), thus the conditions hold and Neimark-Sacker bifurcation occurs at a � a 0 . From numerical simulation, we observed that the Neimark-Sacker bifurcation occurs with respect to the reproduction rate parameter a at the critical value a 0 � 3.01 around the equilibrium (0.5062, 0.2632, 0.1743).

Numerical Simulations
e dynamic behavior of this system is studied by numerical analysis. e bifurcation diagrams of equation (5) are shown in Figures 1(a)-1(c) for each species x, y, and z with b � 3.7, c � 3, d � 3.5, and r � 3.8, as the parameter a changes from 2 to 4.3. We now discuss the salient features of the bifurcation diagram. Because of the similarity of the bifurcation diagrams, only Figure 1(b) is shown magnified in Figure 2 to analyze the dynamic properties.
When the reproduction rate a changes between 2 and 4.3, equation (5) generates complicated features. e equilibrium point E 3 is a positive fixed point that is asymptotically stable for 2.94 < a < 2.99. e phase portraits of various a corresponding to Figure 2(a) are plotted in Figures 3(a) and 3(b). When the parameter a changes from 2.98 to 3.01, the equilibrium point E 3 becomes unstable, and at a � 3.01, a Neimark-Sacker bifurcation appears. e size and smoothness of a closed curve around the equilibrium point change continuously with the change of parameters for 3.01 < a < 3.6, and finally the closed curve converges to the equilibrium point E 3 again at a � 3.6. e phase portraits of various a corresponding to Figure 2(a) are plotted in Figures 3(c)-3(e).
When the parameter a passes through the range (3.6, 3.87), the system has an unstable period-1 orbit, and a Neimark-Sacker bifurcation occurs again at a � 3.87. ere   Discrete Dynamics in Nature and Society  Discrete Dynamics in Nature and Society is also a closed curve whose shape and size are constantly changing for 3.87 < a < 3.904. At a � 3.904, this closed curve breaks, and the system shows period-17 orbits. With the parameter a changing from 3.904 to 3.998, the system shows a cascade of period-doubling bifurcations leading to a chaotic attractor finally. e phase portraits of various a corresponding to Figure 2

The Controlling Chaos of the Three-Species Food Chain Model
As shown in Figure 4, for b � 3.7, c � 3, d � 3.5, and r � 3.8, when a � 4.26, the dynamic behavior of the system is chaotic. At the same time, there is a chaotic attractor, which is the closure of the unstable manifolds of the saddle points. And there is an infinite number of unstable periodic orbits in the chaotic attractor. ere is an unstable period-1 orbit embedding in the chaotic attractor. e chaotic motion of the three-species food chain model is controlled on the stable periodic period-1 orbit using the method in [18]. Control parameter a is perturbed slightly with time. When the mapping point wanders to the neighborhood of the periodic-1 orbit, the control parameter is perturbed.
Assume the system is written in the following form: F is sufficiently smooth, and a is an externally adjustable real parameter. Let the control parameter a be a variable parameter near the rated value a � 4.26, at b � 3.7, c � 3, d � 3.5, and r � 3.8. By equation (6), the equilibrium E 3 is (x * , y * , z * ) � (0.5062, 0.2632, 0.1743). Now, the aim is to change the parameter such that the chaotic attractor involves almost all of the initial conditions, so that the dynamic behavior of the system converges to the desired periodic orbit. Due to the ergodicity of the chaos dynasmics, when the state trajectory enter the vicinity of the unstable periodic orbit to be stabilized, a feedback control law is applied to control the trajectory to move to the desired unstable periodic orbit.
Let W * (a) denotes the unstable point E 3 ; by first-order Taylor expansion, equation (5) can be described as where A is the derivative matrix of F(w, a) to the variable w, and w is the variable set of (x, y, z), A � D w F(w, a), and B is the derivative matrix of F(w, a) to the variable a, B � D a F(w, a), Put the equilibrium point e time-dependent control parameter a is in the form of a linear function with respect to the variable: Substituting equation (32) into equation (28), If the modulus of the matrix A − BK T eigenvalues is less than 1, then the equilibrium point W * (a) is stable. According to Ogata [34], we can find the matrix C 3×3 is controllable matrix: e solution of pole placement is found out through the matrix K T , which can stabilize the chaotic motion to a stable periodic orbit, where K T � α 3 − a 3 α 2 − a 2 α 1 − a 1 T − 1 and T � CQ, Q is a matrix of order 3.
|K T (W i − W * (a))| < δ can be obtained by |a − a| < δ and equation (35). When W i goes into this region, whose width is 2δ/|K T |, the parameter is slightly perturbed; otherwise, the parameter is not perturbed. e controlling equation for the parameter a is where u is a piecewise function and δ is a small positive number whose value affects the control time.
When δ � 0.063, the chaotic motion can be controlled on the period-1 orbit at n � 990 at α 1 � − 0.073, α 2 � − 0.073, α 3 � 0 (as shown in Figures 6(a)-6(c)). According to the numerical simulation, when different values of δare taken, there is a great difference in the time of chaos control. Discrete Dynamics in Nature and Society

Conclusion
It is shown that the equilibrium point of the system loses its stability and Neimark-Sacker bifurcation appears leading to chaos for some parameter values. e chaotic food chain is controlled on the stable periodic period-1 to obtain steadystate orbit using implementation of hybrid control strategy, which is based on feedback control methodology and parameter perturbation. Numerical simulations are presented to illustrate our results with the theoretical analysis and show the effect of the control method. e results of this research prove that the chaos control method can be extended to the higher-dimensional biological model and can be realized. In the study of bifurcation, it is easier to obtain bifurcation parameters than the central manifold theorem. Compared with other methods, chaos control will save more control time and be easier to realize chaos control.
is paper only studies the dynamic behavior of the system in mathematics. In ecology and biology, how to control the parameters a, b, c, d, and r properly to achieve an ecological balance, so as not to make a species extinct or to cause damage to the ecological environment will be a very meaningful and challenging topic, which will be paid more attention in the future research. When the parameter is in a certain range, there should be subcritical or supercritical bifuraction. Whether there will be flip bifurcation is also a question worthy of discussion. I will discuss this problem in the future research work.

Data Availability
No data, models, or codes were generated or used during the study.

Conflicts of Interest
e authors declare that they have no conflicts of interest regarding the publication of this paper.