# Chaotic Dynamics and Chaos Control of Hassell-Type Recruitment Population Model

Authors: Guo Feng
Journal: Discrete Dynamics in Nature and Society
 (2020)

---

## Abstract

For certain parameters, the mapping of a Hassell-type recruitment population model has a chaotic attractor. The control parameter is disturbed slightly with time by the improvement OGY method. When the mapping point wanders to the neighborhood of the periodic point, the control parameter is perturbed. The chaotic motion is controlled on the stable periodic period-1 point and period-2 orbits, and the influence of different control parameter ranges on the control average time is analyzed. When the selected regulator poles are different, the number of iterations used to control chaotic motion on a stable periodic orbit is different. Numerical simulations are presented to illustrate our results with the theoretical analysis and show the effect of the control method.


---

## Body

## 1. Introduction

Leslie models [1] with nonlinear fertility and mortality can have complicated dynamical behavior. Ugarcovici and Weiss [2] studied the Ricker model. This model is described by the two-dimensional mapping Ra,b:ℝ+2⟶ℝ+2:(1)Ra,bx,y=ax+γaye−λx+y,bx,where x and y stand for the density of the first age group and the second age group. a and γa are the group’s initial fertility rates (a, γ > 0), b is the survival rate from the first age group to the second one, and λ is the decay index, λ > 0. In equation (1), the fertility rate monotonically decreases as a function of the total population size, and the fertility decay is exponential. The other model is Hassell model. It is described by the two-dimensional mapping Ha,b:ℝ+2⟶ℝ+2:(2)Ha,bx,y=ax+γay1+x+y−β,bx,where a and γa are the group’s initial fertility rates, b is the survival rate from the first age group to the second one, and β is the decay index, β > 1. In equation (2), the fertility rate monotonically decreases as a function of the total population size, and the fertility decay is polynomial.For some parameter values, these models admit an ergodic attractor which supports a unique physical probability measure. This physical measure satisfies in the strongest possible sense the population biologist’s requirement for ergodicity in their population models. The Lebesgue measure basin is the whole planeℝ+2; Hassell [3], Wikan and Mjølhus [4], and Ugarcovici and Weiss [2] showed that Ricker mapping and Hassell mapping produce Hénon-like chaotic attractors.Here, we present a procedure to contend chaotic dynamics in a Hassell-type recruitment population model. Asa is between 0 and 60, b = 0.7, γ=0.8,andβ=22, Hassell mapping shows rich chaotic dynamic behavior. In biology or ecology, the complex chaotic behavior of this mapping shows the relationship between the number, birth rate, and survival rate in a population, whether it survives in a balanced state or makes the population develop in disorder or chaos. This research can provide theoretical basis and help for the research in biology or ecology. For example, it is applied in marine fishing or in the reproduction and population growth of a certain species in nature.The study of chaos control began in the late 1980s. The early research idea is to use the existing dynamic control strategies and destroy the conditions for the occurrence of chaotic motion. Ott et al. [5] proposed a method to control chaos. The basic idea of controlling chaos can be understood by considering the following one-dimensional logistic map, one of the best studied chaotic systems:(3)xn+1=fxn,r=rxn1−xn,where x is restricted to the unit interval [0 1] and r is a control parameter. The logistic map in the neighborhood of a periodic orbit can be approximated by a linear equation expanded around the periodic orbit. Denote the target period-m orbit to be controlled as x(i), i = 1, 2, …, m, where x(i + 1) = f(x(i)) and x(m + 1) = x(1). Assume that, at time n, the trajectory falls into the neighborhood of component i of the period-m orbit. The linearized dynamics in the neighborhood of component (i+1) is(4)xn+1−xi+1=∂f∂xxn−xi+∂f∂rΔrn=r01−2xixn−xi+xi1−xiΔrn,where the partial derivatives are evaluated at x = x (i) and r=r0. To make xn+1 stay in the neighborhood of x (i + 1), let xn+1−xi+1=0, which gives(5)Δrn=r02xi−1xn−xixi1−xi.Equation (5) holds only when the trajectory xn enters a small neighborhood of the period-m orbit, i.e., when xn−xi≪1, and hence, the required parameter perturbation Δrn is small. Let the length of a small interval defining the neighborhood around each component of the period-m orbit be 2ε. In general, the required maximum parameter perturbation δ is proportional to ε. Since ε can be chosen to be arbitrarily small, δ can also be made arbitrarily small. The average transient time before a trajectory enters the neighborhood of the target periodic orbit depends on ε (or δ). When the trajectory is outside the neighborhood of the target periodic orbit, any parameter perturbation is not applied, so the system evolves at its nominal parameter value r0. Hence, we set Δrn=0 when Δrn>δ. The parameter perturbation Δrn depends on xn and is time-dependent.Pyragas [6], Shinbrot et al. [7], Pyragas [8], and Kocarev and Parlitz [9] have put forward different improvement measures according to various situations and further developed the OGY method, which has laid a good foundation for the application of chaos. Flipe et al. [10] adopt the pole assignment technique in system control to further improve the OGY method. The chaos of high periodic states and high-dimensional dynamic systems in chaotic attractors is controlled by the improved OGY method.By the OGY method, Guo et al. [11, 12] studied the chaos control of two-dimensional Lauwerier mapping(6)xn+1=bxn1−2yn+yn,yn+1=ayn1−ynand a piecewise linear Lozi mapping(7)xn+1=1−axn+yn,yn+1=bxn.The chaotic motions of equations (6) and (7) are both controlled on periodic-1 and periodic-2 orbits. Guo et al. [13] studied the chaos control of two-degree-of-freedom collision vibration system by the OGY method. A four-dimensional discontinuous system is studied. The chaotic motion is controlled on periodic-1 and periodic-2 orbits.The Hassell mapping produces a chaotic attractor asa = 31, b = 0.7, γ=0.8,andβ=22. The unstable period-1 point and unstable period-2 points are controlled on the periodic orbits.


## 2. The Hassell Recruitment Population Model

The Hassell recruitment population model has two generations [3]. It is described by the two-dimensional mapping Ha,b:ℝ+2⟶ℝ+2:(8)Ha,bx,y:fx,y=ax+γay1+x+y−β,gx,y=bx,where x and y stand for the density of the first age group and the second one. a and γa are the fertility rates of the group at time zero, b is the survival rate from the first age group to the second one, and β is the decay index, β > 1. When β = 1, the case corresponds to the model of Beverton and Holt [14]. The attractors are fixed points.For the parametersγ=0.8andβ=22 and fixed survival rate b = 0.7, when the initial fertility rate a changes between 0 and 35, equation (8) generates complicated features, shown in Figure1(a)–1(e). For 0 < a < 0.455, (0, 0) is a global attractor. There exists a positive fixed point that is asymptotically stable for 0.455 < a < 7. For a between 7 and 10, unstable coexistence between the first age group and the second one can be observed, with chaotic bands having narrow and wide periodic windows. Subsequently, a cascade of period-halving is shows in Figure 1(b). For a between 10 and 15, the system goes through a quasi-periodic region (Hopf bifurcation occurs at a = 12.5, including frequency lockings which appear as a collapse of the invariant circle to a periodic orbit), beyond which equation (8) coexists in a stable manner for a large range of the parameter a (as shown in Figure 1(c)). The system enters quasi-periodicity again when a = 15. For a between 15 and 17, it coexists in an unstable manner again, with chaotic bands having narrow and wide periodic windows. For a between 17 and 25, Hopf bifurcation occurs again (as shown in Figure 1(d)). For a between 25 and 29, the bifurcation diagram is a periodic window with a periodic attractor, but that it includes ambiguous dynamic behavior. Several attractors coexist in this region (as shown in Figure 1(e)). For a greater than 31, a chaotic attractor appears. A Hénon chaotic attractor [15] coexists with an attracting invariant curve (as shown in Figure 2). For a between 0 and 60, the global bifurcation is shown in Figure 1(a).Figure 1
Bifurcation diagrams of the Hassell mapping (γ=0.8,β=22, b = 0.7, and 0 < a < 35).
(a)
(b)
(c)
(d)
(e)Figure 2
The chaotic attractor of the Hassell mapping with unstable periodic orbits of period-1(.) and period-2(+).
(a)
(b)The largest Lyapunov exponent measures the average exponential rate of divergence or convergence of nearby orbits in the phase space. Whenever two nearby orbits move too far apart, one of the orbits has to be moved back to the vicinity of the other along the line of separation. For a chaotic attractor, the largest Lyapunov exponent must be positive. If the largest Lyapunov exponent is negative, this implies a stable state or a periodic attractor. The largest Lyapunov exponents of the Hassell mapping have been calculated and plotted in Figure3. When 0 < a < 15, the largest Lyapunov exponent is negative. When 15 < a < 30, the largest Lyapunov exponent can be negative or positive, and the system has Hopf bifurcation and enters quasi-periodicity and at last generates the chaotic state. When a > 30, the largest Lyapunov exponent is positive. The dynamics behavior is chaotic. And the chaotic attractor is found as LE max = 2.1392e0.3.Figure 3
Lyapunov exponent of Hassell mapping (γ=0.8,β=22, b = 0.7, and 0 < a < 40).
(a)
(b)
(c)


## 3. Pole Assignment Technique for Controlling Chaos

The system is first written in the form of a discrete time [10]:(9)Zi+1=FZi,a,Zi∈R2,a∈R.F is sufficiently smooth and a is an externally adjustable real parameter. That is, it requires a−a¯<δ at some time, and a¯ is a rated value. It is assumed that there is a chaotic attractor of equation (9) for a=a¯. Now, the aim is to change the parameters 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 in the attractor. By the OGY method, due to the ergodicity of the chaos dynamics, when the state trajectory enters 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.IfZ∗a¯ is the unstable fixed point, by first-order Taylor expansion, equation (9) can be written as(10)Zi+1−Z∗a¯=AZi−Z∗a¯+Ba−a¯.Find out the values of the matricesA and B at Z=Z∗a¯ and a=a¯, where A is the partial derivative matrix of F (z, a) to z, A=DzFz,a, and B is the partial derivative matrix of F (z, a) to a, B=DaFz,a. The time-dependent control parameter a is in the form of a linear function with respect to the variable(11)a−a¯=−KTZi−Z∗a¯.Replace equations (11) into (10):(12)Zi+1−Z∗a¯=A−BKTZi−Z∗a¯.So, as long as the matrixA−BKT is asymptotically stable, that is, if the modulus of its eigenvalues is less than 1, the fixed point Z∗a¯ is stable. The following key questions is how to determine the matrix KT, which can stabilize the chaotic motion at a stable periodic point. The pole assignment is solved according to [16]. The matrix Cn×n is a controllable matrix, and the rank is n:(13)C=BABA2B⋅⋅⋅An−1B.The solution of pole assignment is given byKT=αn−an,…,α1−a1T−1, where T=CW, W is a matrix of order n:(14)W=an−1an−2…a11an−2an−3…10……………a11…0010…00,where aii=1,…,n are the coefficients of the characteristic polynomial of the matrix A, that is,(15)detλI−A=λn+a1λn−1+⋅⋅⋅+an,and α1,α2,…,αn are the coefficients of the characteristic polynomial detA−BKT, that is,(16)detsI−A−BKT=sn+α1sn−1+⋅⋅⋅+αn.After working outKT,KTZi−Z∗a¯<δ is obtained by a−a¯<δ and equation (12).


## 4. The Chaos Control of Hassell Mapping

As shown in Figures2(a) and 2(b), for γ=0.8,β=22,andb=0.7, when a = 31, the dynamic has period-1 point and period-2 points. The dynamic behavior of the mapping is chaotic at a¯=31. 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. There are an unstable fixed point and two unstable period-2 points embedding in the chaotic attractor.


### 4.1. Control of Period-1 Point of Hassell Mapping

Let the control parametera be a variable near the rated value a¯=31 at γ=0.8,β=22,andb=0.7, By fx,y=x and gx,y=y, the fixed point is x∗,y∗=0.46,0.11. Bring the fixed point into the matrixes A and B:(17)A=a1+1x−βx+y−βγy1+x+y−β−1aγ+γx−βx+γy−γyβ1+x+y−β−1b0x∗,y∗=−0.010−0.0111.50,B=x+γy1+x+y−β0x∗,y∗=1.440.A controllable matrix is a matrix with rank 2 as(18)C=BAB=1.44−0.14402.16.The solution of the pole assignment problem is given byKT=α2−a2,α1−a1T−1, where T=CW, W=a1110, where aii=1,2 are the coefficients of the characteristic polynomial of the matrix A.(19)detλI−A=detλ+0.0100.011−1.5λ=λ2+0.010λ+0.017.So,a1=0.010anda2=0.017.(20)T=CW=1.44−0.14402.160.010110=−0.131.442.160,T−1=00.46300.69440.0418.The characteristic roots ofA which can also be obtained at the fixed point x∗,y∗=0.46,0.11 are λs=−0.005−0.13i,λu=−0.005+0.13i. α1,α2 are the coefficients of the characteristic polynomial of the matrix A−BKT. The assumed characteristic roots μ1,μ2 are called the adjustment values, that is,(21)detsI−A−BKT==s2−μ1+μ2s+μ1μ2.The relationships between the roots and the coefficients are obtained as(22)α1=−μ1+μ2,α2=μ1μ2.Whenμ1=1, α1=−1−α2; when μ1=−1, α1=1+α2; and when μ1μ2=1, α2=1 according to the range of the values α1 and α2 determined by equation (22), as in the shaded area in Figure 4.Figure 4
Choice of regulator pole area (triangle regionMNQ).From Ogata [16] and equation (9), the matrix KT is not unique. As long as the matrix KT is obtained with the values α1 and α2 in the triangle region MNQ as Figure 4, it can make the matrix A−BKT be asymptotically stable, that is, the modulus of its eigenvalues is less than 1, so we can take μ1=0,μ2=λs. As α1,α2=−λs,0, we can obtain(23)KT=0−a2,−λs−a1T−1=−0.017,−0.12.WhenKT is found, KTZi−Z∗a¯<δ is obtained by a−a¯<δ and equation (12). There is a region whose width is 2δ/KT, and when Zi is in this region, the parameter is controlled; otherwise, the parameter is not controlled. The control rate is given by the following formula:(24)a−a¯=−KTZi−Z∗a¯×uδ−KTZi−Z∗a¯,where u is a step function:(25)uα=0,α<0,1,α>0.As shown in Figure5, we get the different values in the triangle region of Figure 4. When we choose α1=1.75andα2=0.74, the chaotic motion can be controlled on the period-1 orbit at n = 5760 (as shown in Figures 5(a) and 5(b)). When we choose α1=−0.5andα2=0.6, the chaotic motion can be controlled on the period-1 orbit at n = 3300 (as shown in Figures 5(c) and 5(d)). When the values of α1 and α2 are different from the previous values, the number of equation (8) iterations is different to control chaos.Figure 5
Control of period-1 of Hassell mapping.
(a)
(b)
(c)
(d)


### 4.2. Control of Period-2 of Hassell Mapping

By iterating equation (8), we can obtain(26)HHa,bx,y:Fx,y,a,Gx,y,a,where(27)Fx,y,a=a2x+γy1+x+y−β+γay1+ax+γay1+x+y−β+y−β,Gx,y,a=b2x.Whena¯=31, the points x1,y1 and x2,y2 of period two are (0.15, 0.26) and (0.34, 0.056), which are obtained as follows: Fx,y,a=x,Gx,y,a=y.According to(28)A1=∂Fx,y,a∂x∂Fx,y,a∂y∂Gx,y,a∂x∂Gx,y,a∂yx1,y1,A2=∂Fx,y,a∂x∂Fx,y,a∂y∂Gx,y,a∂x∂Gx,y,a∂yx2,y2,B1=∂Fx,y,a∂a∂Gx,y,a∂ax1,y1,B2=∂Fx,y,a∂a∂Gx,y,a∂ax2,y2,where(29)∂Fx,y,a∂x=1+ax+γay1+x+y−β+y−βa21+x+γy+a2x+γy1+x+y−β+γaya+ax+γay1+x+y−β,(30)∂Fx,y,a∂y=1+ax+γay1+x+y−β+y−βa21+x+y−β1+x+γy+γaa2x+γy1+x+y−β+γay1+x+y−βa+ax+γay+1,(31)∂Gx,y,a∂x=b2,∂Gx,y,a∂y=0,(32)∂Fx,y,a∂a=1+ax+γay1+x+y−β+y−β2ax+γy1+x+y−β+γya2x+γy1+x+y−β+γay1+x+y−βx+γy,(33)∂Gx,y,a∂a=0.Bringx1,y1=0.15,0.26 and x2,y2=0.34,0.056 into equation (29)–(33):(34)A1=−1.54−0.980−3.235,A2=0.080.280−3.28,B1=00.34,B2=00.44.The controllable matrices are obtained as follows:(35)C1=B1A1B1=0−0.330.34−1.10,C2=B2A2B2=00.120.44−1.44,and aii=1,2 are the coefficients of the characteristic polynomial of the matrix A1, a1=4.77,a2=4.97. At the same time, the characteristic roots are λ1s=−0.54,λ1u=−3.23.The solution of the pole assignment problem is given by(36)K1T=α2−a2,α1−a1T1−1,where α1andα2 are the coefficients of the characteristic polynomial of the matrix A1−B1K1T. According to the discussion of the values α1,α2, we take α1,α2=−λ1s,0, where T1=C1W1, so(37)K1T=0−a2,−λ1s−a1T1−1=−4.97,−4.23,aii=3,4 are the coefficients of the characteristic polynomial of the matrix A2−B2K2T, a3=3.2,a4=−0.26, and the characteristic roots are λ2s=0.08,λ2u=−3.28. We take α3,α4=−λ2s,0, and we obtain K2T=−0.46,4.39. The control rate can be given by the following equation:(38)a−a¯=−KiTZn−Zia¯×uδ−KiTZn−Zia¯,i=1,2.As shown in Figure6, we get the different values in the triangle region MNQ of Figure 4. When we choose α1=−0.37,α2=0.3,α3=−0.25,andα4=−0.54, the chaotic motion can be controlled on the period-1 orbit at n = 6554 (as shown in Figures 6(a) and 6(b)). When we choose α1=0.73,α4=0.24,α3=0.34,andα4=0.44, the chaotic motion can be controlled on the period-1 orbit at n = 2700 (as shown in Figures 6(c) and 6(d)).Figure 6
Control of period-2 of Hassell mapping.
(a)
(b)
(c)
(d)


## 4.1. Control of Period-1 Point of Hassell Mapping

Let the control parametera be a variable near the rated value a¯=31 at γ=0.8,β=22,andb=0.7, By fx,y=x and gx,y=y, the fixed point is x∗,y∗=0.46,0.11. Bring the fixed point into the matrixes A and B:(17)A=a1+1x−βx+y−βγy1+x+y−β−1aγ+γx−βx+γy−γyβ1+x+y−β−1b0x∗,y∗=−0.010−0.0111.50,B=x+γy1+x+y−β0x∗,y∗=1.440.A controllable matrix is a matrix with rank 2 as(18)C=BAB=1.44−0.14402.16.The solution of the pole assignment problem is given byKT=α2−a2,α1−a1T−1, where T=CW, W=a1110, where aii=1,2 are the coefficients of the characteristic polynomial of the matrix A.(19)detλI−A=detλ+0.0100.011−1.5λ=λ2+0.010λ+0.017.So,a1=0.010anda2=0.017.(20)T=CW=1.44−0.14402.160.010110=−0.131.442.160,T−1=00.46300.69440.0418.The characteristic roots ofA which can also be obtained at the fixed point x∗,y∗=0.46,0.11 are λs=−0.005−0.13i,λu=−0.005+0.13i. α1,α2 are the coefficients of the characteristic polynomial of the matrix A−BKT. The assumed characteristic roots μ1,μ2 are called the adjustment values, that is,(21)detsI−A−BKT==s2−μ1+μ2s+μ1μ2.The relationships between the roots and the coefficients are obtained as(22)α1=−μ1+μ2,α2=μ1μ2.Whenμ1=1, α1=−1−α2; when μ1=−1, α1=1+α2; and when μ1μ2=1, α2=1 according to the range of the values α1 and α2 determined by equation (22), as in the shaded area in Figure 4.Figure 4
Choice of regulator pole area (triangle regionMNQ).From Ogata [16] and equation (9), the matrix KT is not unique. As long as the matrix KT is obtained with the values α1 and α2 in the triangle region MNQ as Figure 4, it can make the matrix A−BKT be asymptotically stable, that is, the modulus of its eigenvalues is less than 1, so we can take μ1=0,μ2=λs. As α1,α2=−λs,0, we can obtain(23)KT=0−a2,−λs−a1T−1=−0.017,−0.12.WhenKT is found, KTZi−Z∗a¯<δ is obtained by a−a¯<δ and equation (12). There is a region whose width is 2δ/KT, and when Zi is in this region, the parameter is controlled; otherwise, the parameter is not controlled. The control rate is given by the following formula:(24)a−a¯=−KTZi−Z∗a¯×uδ−KTZi−Z∗a¯,where u is a step function:(25)uα=0,α<0,1,α>0.As shown in Figure5, we get the different values in the triangle region of Figure 4. When we choose α1=1.75andα2=0.74, the chaotic motion can be controlled on the period-1 orbit at n = 5760 (as shown in Figures 5(a) and 5(b)). When we choose α1=−0.5andα2=0.6, the chaotic motion can be controlled on the period-1 orbit at n = 3300 (as shown in Figures 5(c) and 5(d)). When the values of α1 and α2 are different from the previous values, the number of equation (8) iterations is different to control chaos.Figure 5
Control of period-1 of Hassell mapping.
(a)
(b)
(c)
(d)


## 4.2. Control of Period-2 of Hassell Mapping

By iterating equation (8), we can obtain(26)HHa,bx,y:Fx,y,a,Gx,y,a,where(27)Fx,y,a=a2x+γy1+x+y−β+γay1+ax+γay1+x+y−β+y−β,Gx,y,a=b2x.Whena¯=31, the points x1,y1 and x2,y2 of period two are (0.15, 0.26) and (0.34, 0.056), which are obtained as follows: Fx,y,a=x,Gx,y,a=y.According to(28)A1=∂Fx,y,a∂x∂Fx,y,a∂y∂Gx,y,a∂x∂Gx,y,a∂yx1,y1,A2=∂Fx,y,a∂x∂Fx,y,a∂y∂Gx,y,a∂x∂Gx,y,a∂yx2,y2,B1=∂Fx,y,a∂a∂Gx,y,a∂ax1,y1,B2=∂Fx,y,a∂a∂Gx,y,a∂ax2,y2,where(29)∂Fx,y,a∂x=1+ax+γay1+x+y−β+y−βa21+x+γy+a2x+γy1+x+y−β+γaya+ax+γay1+x+y−β,(30)∂Fx,y,a∂y=1+ax+γay1+x+y−β+y−βa21+x+y−β1+x+γy+γaa2x+γy1+x+y−β+γay1+x+y−βa+ax+γay+1,(31)∂Gx,y,a∂x=b2,∂Gx,y,a∂y=0,(32)∂Fx,y,a∂a=1+ax+γay1+x+y−β+y−β2ax+γy1+x+y−β+γya2x+γy1+x+y−β+γay1+x+y−βx+γy,(33)∂Gx,y,a∂a=0.Bringx1,y1=0.15,0.26 and x2,y2=0.34,0.056 into equation (29)–(33):(34)A1=−1.54−0.980−3.235,A2=0.080.280−3.28,B1=00.34,B2=00.44.The controllable matrices are obtained as follows:(35)C1=B1A1B1=0−0.330.34−1.10,C2=B2A2B2=00.120.44−1.44,and aii=1,2 are the coefficients of the characteristic polynomial of the matrix A1, a1=4.77,a2=4.97. At the same time, the characteristic roots are λ1s=−0.54,λ1u=−3.23.The solution of the pole assignment problem is given by(36)K1T=α2−a2,α1−a1T1−1,where α1andα2 are the coefficients of the characteristic polynomial of the matrix A1−B1K1T. According to the discussion of the values α1,α2, we take α1,α2=−λ1s,0, where T1=C1W1, so(37)K1T=0−a2,−λ1s−a1T1−1=−4.97,−4.23,aii=3,4 are the coefficients of the characteristic polynomial of the matrix A2−B2K2T, a3=3.2,a4=−0.26, and the characteristic roots are λ2s=0.08,λ2u=−3.28. We take α3,α4=−λ2s,0, and we obtain K2T=−0.46,4.39. The control rate can be given by the following equation:(38)a−a¯=−KiTZn−Zia¯×uδ−KiTZn−Zia¯,i=1,2.As shown in Figure6, we get the different values in the triangle region MNQ of Figure 4. When we choose α1=−0.37,α2=0.3,α3=−0.25,andα4=−0.54, the chaotic motion can be controlled on the period-1 orbit at n = 6554 (as shown in Figures 6(a) and 6(b)). When we choose α1=0.73,α4=0.24,α3=0.34,andα4=0.44, the chaotic motion can be controlled on the period-1 orbit at n = 2700 (as shown in Figures 6(c) and 6(d)).Figure 6
Control of period-2 of Hassell mapping.
(a)
(b)
(c)
(d)


## 5. Conclusion

The dynamic complexities of the Hassell recruitment population model have been analyzed. These results demonstrate that the reproductive rate and the resulting population growth can also act as a strong destabilizing factor, leading to rich dynamic behaviors. By improving the method of the OGY, the chaotic control of the Hassell mapping is controlled, the perturbation quantity of the control parameter is selected by the pole configuration method in the linear control theory, and unstable period-1 and unstable period-2 are controlled to be stable period orbits. At the same time, the adjustment values to be selected are different, and the number of the mapping iterations is different to control chaos. By numerical simulation, the effectiveness of the method is demonstrated. This research can provide theoretical basis and help for the research in biology or ecology. According to the research results, the birth ratea, survival rate b, or decay coefficient β can be properly adjusted to achieve an ecological balance.



---

*Source: 8148634-2020-04-09.xml*