Bifurcation and Chaos of a Discrete-Time Population Model

A Leslie population model for two generations is investigated by qualitative analysis and numerical simulation. For the different parameters a and b in the model, the dynamics of the system are studied, respectively. It shows many complex dynamic behavior, including several types of bifurcations leading to chaos, such as period-doubling bifurcations and Neimark–Sacker bifurcations. With the change of parameters, attractor crises and chaotic bands with periodic windows appear.-e largest Lyapunov exponents are numerically computed and can verify the rationality of the theoretical analysis.


Introduction
In 1964, Hénon took KAM theorem as the background and found that there was deterministic random behavior in twodimensional nonintegrable Hamiltonian system, that is claimed Hénon attractor. Ruelle and Takens proposed the concept of "strange attractor," which promoted the research of Smale horseshoe attractor [1]. Lorenz [2] pointed out that there must be a connection between the inexact recurrence of climate and the inability of long-term weather forecast. It is found that the chaos phenomenon is "extremely sensitive to the initial conditions." e word "chaos" is formally used by Li and Yorke [3], which is considered as the first formal expression of chaos theory. In [1], an example of horseshoe mapping is given, which opened the mathematical method of studying chaos. Lasota [4] studied the initial value problem of the first order nonlinear partial differential equation of the wormhole model. Brunovsky [5] gave the definition of chaotic mapping.
A single species continuous-time population model is studied by [6] and later studied further by Pearl and Reed [7].
May [8] studied several discrete-time models with a singlespecies displaying chaotic behavior for certain parameter values. Beddington et al. [9] studied that discrete-time host, parasitoid models show complex dynamic behavior.
Leslie [10,11] introduced an age-structured linear population model. Later, this kind of discrete-time models contain linear or nonlinearity forms, which are called Leslie models. Ugarcovici and Weiss [12] studied the dynamical behavior of Ricker model. is model is described by the two-dimensional mapping R a,b : R 2 where x and y stand for the density of the first age group and the second age group. a and ca are the group's initial fertility rates (a, c > 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. e other mold is the Hassell model. It is described by the two-dimensional mapping H a,b : R 2 where a and ca 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. Guo [13] studied the chaos control of the population model.
For some parameter values, these models admit an ergodic attractor which supports a unique physical probability measure.
is physical measure satisfies in the strongest possible sense the population biologist's requirement for ergodicity in their population models. Wikan and Mjølhus [14] and Ugarcovici and Weiss [12] showed that Ricker mapping and Hassell mapping produce Hénon-like chaotic attractors.
Here, for the different parameters a and b in the model, the dynamics of the system is studied, respectively. It shows many complex dynamic behavior, including several types of bifurcations leading to chaos, such as period-doubling bifurcations and Neimark-Sacker bifurcations. With the change of parameters, attractor crises and chaotic bands with periodic windows appear. e largest Lyapunov exponents are numerically computed and can verify the rationality of the theoretical analysis. Zhu and Zhao [15] studied the dynamic complexity of a host, parasitoid ecological model with the Hassell growth function for the host. Liu and Xiao [16] studied the complex dynamic behaviors of a discrete-time predator-prey system.
In biology or ecology, the complex chaotic behavior of this mapping shows the relationship among 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. is research can provide theoretical basis and help for the research in biology or ecology.
rough the study of these complex dynamic properties of the population model, it can guide people to make corresponding meaningful work to keep the ecological balance. For example, it is applied in the marine fishing, or in the reproduction and population growth of a certain species in nature.

The Nonlinear Hassell Population Model
Hassell model is described by the two-dimensional mapping H a,b : R 2 where where x t and y t are the density of the first age group and the second one. a and ca 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 (3), the fertility rate monotonically decreases as a function of the total population size and the fertility decay is polynomial.

Stability Analysis of Hassell
Population Model e equation (3) has one equilibrium point: E(x * , y * ), where (x * , y * ) are positive and satisfy bx * � y * . (6) e Taylor series expansion of equation (3) at the equilibrium point: E(x * , y * ) is written as where Let the matrix zF/zx zF/zy zG/zx zG/zy (x * ,y * ) be A � P 11 P 12 where e characteristic equation of equation (9) is Equation (11) can be rewritten as where 2 Discrete Dynamics in Nature and Society e roots of equation (12) are Both eigenvalues are real numbers and |λ 1,2 | < 1, if Equation (15) can be solved as e eigenvalues λ 1,2 become complex numbers and are inside the unit circle in the complex λ-plane for at is, e positive equilibrium point E(x * , y * ) is stable under conditions (17) and (18).

Bifurcation Analysis
Equation (3) cannot be solved analytically, and therefore its long-term behavior must be studied by numerical simulation. Figure 1 shows the bifurcation diagrams of equation (3) for the density of the first age group x with c � 0.8, β � 22, and b � 0.7, as the parameter a increases. Figure 1 is magnified as Figure 2.
When the initial fertility rate a changes between 0 and 60, equation (3) shows complicated features. For 0 < a < 0.75, (0, 0) is a global attractor. ere exists a positive fixed point that is asymptotically stable for 0.75 < a < 15.8. At a � 15.8, it appears a cascade of period-3 orbit. For a between 15.8 and 24, it exists a period-3 window which is embedded in the strange attractor for a long time. en, each of the period-3 orbits begins to undergo a flip bifurcation, leading to chaos. e phase portraits of various a corresponding to Figure 2(a) are plotted in Figures 3(a)-3(b). Tianyan-Li and Yorke [3] proved that if a system has period-3 point, it has all periodic points. is is also verified in the following analysis as a increases.
When a passes through the range (24.4, 27), the bifurcation diagram in Figure 2(b) shows that this window is not a periodic window with a cascade of periodic attractors, but that it includes more complex dynamic patterns. It appears that several attractors coexist in this region: period-3 and period-5 attractors. ere is a very wide chaotic band with period-6 orbit meanwhile (see Figure 3(b)).
When a � 27.2, the chaotic band suddenly disappears in a crisis and the system enters a periodic window with a cascade of period-doubling bifurcations leading to a chaotic attractor with periodic windows. It appears period-11 and starts to appear period-9 when a � 27.24. e phase portraits which are associated with Figure 2(b) are disposed in Figures 3(b)-3(c). When a passes through the range (27.24, 30), the detailed bifurcation diagram in Figure 2(c). As a > 28.26, there is a cascade of period-doubling bifurcation leading to a wide chaotic region, characterized by tangent bifurcation and attractor crises. A very narrow periodic window appears again when a � 29. When a increases beyond 28.5, a chaotic attractor abruptly appears and the periodic attractor disappears. e chaotic band has changed into four chaotic bands.
is kind of transilience of the chaotic state caused by the continuous change of parameters is called chaos crisis by Grebogi et al. [17]. Because the unstable periodic orbit meets the secondary chaotic band, the orbits in the chaotic band is filled between all levels of orbits, leading to the emergence of chaos crisis. e phase portraits of various a corresponding to Figure 2(c) are plotted in Figures 3(d)-3(f ). When a passes through the range (30.8, 31.7) in Figure 2(d), a large periodic window appears; then, chaos appears because of period-doubling bifurcation. When a � 49, it appears as period-4, and it appears as period-9 at a � 51.2, at last, equation (6) appears as a chaotic attractor. e phase portraits which are associated with Figure 2(e) are disposed in Figures 3(g) and 3(h).
From the above analysis, it can be seen that the twodimensional nonlinear mapping is a process of intermittently breaking up in accordance with periodic behavior and chaos phenomenon so that the system appears chaotic motion state, which is called Pomeau-Manneville path (through intermittence chaos) by Eckmann [18]. In the process of parameter change, when the periodic window appears, the one of the roots of equation (11)  In Figure 4(a), as b increases from 0 to 0.03, a narrow chaotic band appears. When b � 0.03, the chaotic regime suddenly disappears in a crisis and the system appears a period-2 orbit for b ∈ (0.03, 0.435), and a stable fixed point for b ∈ (0.435, 0.5526). A flip bifurcation occurs at b � 0.435. Figure 5(a) presents the phase portrait of period-2. As b ∈ (0.5526, 0.61) in Figure 4(c), equation (3) goes through a quasi-periodic region (including tangent bifurcation, narrow and wide periodic windows, and frequency lockings which appear as a collapse of the invariant circle to a periodic orbit). We can see that there is a stable fixed point for b < 0.5526, and the fixed point loses its stability as b increases. A Neimark-Sacker bifurcation occurs at b � 0.5526 and an attracting invariant cycle bifurcates from the fixed point. In this case, equation (11) has two conjugate complex roots λ and λ, satisfying |λλ| � 1.
As b is increasing, the smooth invariant loop becomes larger and no longer smooth. At last, the invariant loop breaks. As b � 0.61, there is a cascade of period-doubling bifurcation leading to a wide chaotic region, characterized by tangent bifurcation and attractor crises. Furthermore, when b ∈ (0.565, 0.57) and b ∈ (0.615, 0.619), we can observe the period-5, 8 windows within the chaotic regions, respectively. e phase portraits of various b corresponding to Figure 4(c) are plotted in Figures 5(b)-5(f ).
According to the above analysis, when b is increasing, chaos appears by the Ruelle Takens Newhouse scheme (through Neimark-Sacker bifurcation) and the Pomeau-Manneville path (through intermittence chaos). In the larger region of the parameter space, the path leading to chaos is related to the Neimark-Sacker bifurcation. In these ways, the locking phase and quasiperiodic motion can be observed.

The Largest Lyapunov Exponent
In this section, the largest Lyapunov exponent, which has proven to be the most useful dynamic diagnostic tool for chaotic systems, is considered. is quantity represents the average exponential rate of divergence or convergence of nearby orbits in phase space [19]. For a chaotic attractor, the largest Lyapunov exponent LE max must be positive. If LE max is negative, this implies a stable state or a periodic attractor.
e largest Lyapunov exponents corresponding to the cases shown in Figures 1 and 4 have been calculated and plotted in Figures 6(a) and 6(b), respectively. e existence of chaotic regions in the parametric space is clearly visible in these figures. e largest Lyapunov exponents for the strange attractors were found to be LE max � 20.02,31.81, corresponding to the parameters a and b, respectively. In Figure 6(a), we can easily see that the maximum Lyapunov exponents are negative for a ∈ (0, 25.5), that is to say, equation (3) has no chaotic region. For a ∈ (25.5, 35), some Lyapunov exponents are bigger than 0, some are smaller than 0, so there exists stable fixed point or stable periodic windows in the chaotic region. In Figure 6(b), the  We find that the values of the maximum Lyapunov exponents are consistent with the dynamic properties of the system under different parameters. windows appear. e largest Lyapunov exponents are numerically computed and can verify the rationality of the theoretical analysis. Numerical simulation results not only show the consistence with the theoretical analysis but also display the new and interesting dynamical behaviors, including different periodic orbits in chaotic regions and boundary crisis. When the parameter c changes in different value ranges, what more complex dynamic properties of the model will be a problem worthy of further study.

Conclusion
In biology or ecology, the complex chaotic behavior of this mapping shows the relationship among 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. is research can provide theoretical basis and help for the research in biology or ecology. For example, it is applied in the marine fishing, or in the reproduction and population growth of a certain species in nature.

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

Conflicts of Interest
e authors declare that they have no conflicts of interest.