
This study focused on the dynamical behavior analysis of a discrete fractional Leslie-Gower model incorporating antipredator behavior and a Holling type Ⅱ functional response. Initially, we analyzed the existence and stability of the model's positive equilibrium points. For the interior positive equilibrium points, we investigated the parameter conditions leading to period-doubling bifurcation and Neimark-Sacker bifurcation using the center manifold theorem and bifurcation theory. To effectively control the chaos resulting from these bifurcations, we proposed two chaos control strategies. Numerical simulations were conducted to validate the theoretical results. These findings may contribute to the improved management and preservation of ecological systems.
Citation: Yao Shi, Zhenyu Wang. Bifurcation analysis and chaos control of a discrete fractional-order Leslie-Gower model with fear factor[J]. AIMS Mathematics, 2024, 9(11): 30298-30319. doi: 10.3934/math.20241462
[1] | Yi Tian . Approximate solution of initial boundary value problems for ordinary differential equations with fractal derivative. Mathematical Modelling and Control, 2022, 2(2): 75-80. doi: 10.3934/mmc.2022009 |
[2] | Abduljawad Anwar, Shayma Adil Murad . On the Ulam stability and existence of Lp-solutions for fractional differential and integro-differential equations with Caputo-Hadamard derivative. Mathematical Modelling and Control, 2024, 4(4): 439-458. doi: 10.3934/mmc.2024035 |
[3] | Ravindra Rao, Jagan Mohan Jonnalagadda . Existence of a unique solution to a fourth-order boundary value problem and elastic beam analysis. Mathematical Modelling and Control, 2024, 4(3): 297-306. doi: 10.3934/mmc.2024024 |
[4] | Zhen Xin, Yuhe Yang, Qiaoxia Li . Controllability of nonlinear ordinary differential equations with non-instantaneous impulses. Mathematical Modelling and Control, 2022, 2(1): 1-6. doi: 10.3934/mmc.2022001 |
[5] | K. Venkatachalam, M. Sathish Kumar, P. Jayakumar . Results on non local impulsive implicit Caputo-Hadamard fractional differential equations. Mathematical Modelling and Control, 2024, 4(3): 286-296. doi: 10.3934/mmc.2024023 |
[6] | Iman Malmir . Novel closed-loop controllers for fractional nonlinear quadratic systems. Mathematical Modelling and Control, 2023, 3(4): 345-354. doi: 10.3934/mmc.2023028 |
[7] | Hui Li, Nana Jin, Yu Zhang . Existence of nonoscillatory solutions for higher order nonlinear mixed neutral differential equations. Mathematical Modelling and Control, 2024, 4(4): 417-423. doi: 10.3934/mmc.2024033 |
[8] | Mrutyunjaya Sahoo, Dhabaleswar Mohapatra, S. Chakraverty . Wave solution for time fractional geophysical KdV equation in uncertain environment. Mathematical Modelling and Control, 2025, 5(1): 61-72. doi: 10.3934/mmc.2025005 |
[9] | Kexin Ouyang, Xinmin Qu, Huiqin Lu . Sign-changing and signed solutions for fractional Laplacian equations with critical or supercritical nonlinearity. Mathematical Modelling and Control, 2025, 5(1): 1-14. doi: 10.3934/mmc.2025001 |
[10] | Zhenguo Luo, Liping Luo . New criteria for oscillation of damped fractional partial differential equations. Mathematical Modelling and Control, 2022, 2(4): 219-227. doi: 10.3934/mmc.2022021 |
This study focused on the dynamical behavior analysis of a discrete fractional Leslie-Gower model incorporating antipredator behavior and a Holling type Ⅱ functional response. Initially, we analyzed the existence and stability of the model's positive equilibrium points. For the interior positive equilibrium points, we investigated the parameter conditions leading to period-doubling bifurcation and Neimark-Sacker bifurcation using the center manifold theorem and bifurcation theory. To effectively control the chaos resulting from these bifurcations, we proposed two chaos control strategies. Numerical simulations were conducted to validate the theoretical results. These findings may contribute to the improved management and preservation of ecological systems.
Nonlinear equations arise always in electroanalytical chemistry with complex and esoteric nonlinear terms[1,2], though there are some advanced analytical methods to deal with nonlinear problems, for examples, the Gamma function method[3], Fourier spectral method[4], the reproducing kernel method[5], the perturbation method[6], the homotopy perturbation method[7,8], He's frequency formulation[9,10,11] and the dimensional method[12], chemists are always eager to have a simple one step method for nonlinear equations. This paper introduces an ancient Chinese algorithm called as the Ying Buzu algorithm[13] to solve nonlinear differential equations.
We first introduce the Taylor series method[14]. Considering the nonlinear differential equation:
d2udx2+F(u)=0. | (0.1) |
The boundary conditions are
dudx(a)=α, | (0.2) |
u(b)=β. | (0.3) |
If u(a) is known, we can use an infinite Taylor series to express the exact solution[14]. We assume that
u(a)=c. | (0.4) |
From (0.1), we have
u″(a)=−F(u(a))=−F(c), |
u‴(a)=−∂F(c)∂uu′(a)=−α∂F(c)∂u. |
Other higher order derivatives can be obtained with ease, and its Taylor series solution is
u(x)=u(a)+(x−a)u′(a)+12!(x−a)2u″(a)+13!(x−a)3u‴(a)+...+1N!(x−a)Nu(N)(a), |
the constant c can be determined by the boundary condition of (0.3).
The Ying Buzu algorithm[15,16] was used to solve differential equations in 2006[13], it was further developed to He's frequency formulation for nonlinear oscillators[13,17,18,19,20,21,22,23] and Chun-Hui He's algorithm for numerical simulation[24].
As c in (0.4) is unknown, according to the Ying Buzu algorithm[13,15,16], we can assume two initial guesses:
u1(a)=c1,u2(a)=c2. | (0.5) |
where c1 and c2 are given approximate values.
Using the initial conditions given in (0.2) and (0.5), we can obtain the terminal values:
u(b,c1)=β1,u(b,c2)=β2. |
According to the Ying Buzu algorithm[6,7,8,9,10,11,12], the initial guess can be updated as
u(a)est=c3=c1(β−β2)−c2(β−β1)(β−β2)−(β−β1), |
and its terminal value can be calculated as
u(b,c3)=β3. |
For a given small threshold, ε, |β−β3|≤ε, we obtain u(a)=c3 as an approximate solution.
Here, we take Michaelis Menten dynamics as an example to solve the equation. Michaelis Menten reaction diffusion equation is considered as follows[25,26]:
d2udx2−u1+u=0. | (0.6) |
The boundary conditions of it are as follows:
dudx(0)=0,u(1)=1. | (0.7) |
We assume
u(0)=c. |
From (0.6), we have
u″(0)=c1+c, |
u‴(0)=0, | (0.8) |
u(4)=c(1+c)3. |
The 2nd order Taylor series solution is
u(x)=u(0)+u′(0)1!x+u″(0)2!x2=c+c2(1+c)x2. |
In view of the boundary condition of (0.7), we have
u(1)=c+c2(1+c)=1, | (0.9) |
solving c from (0.9) results in
c=0.7808. |
So we obtain the following approximate solution
u(x)=0.7808+0.2192x2. |
Similarly the fourth order Taylor series solution is
u(x)=c+c2!(1+c)x2+c4!(1+c)3x4. |
Incorporating the boundary condition, u(1)=1, we have
c+c2!(1+c)+c4!(1+c)3=1. | (0.10) |
We use the Ying Buzu algorithm to solve c, and write (0.10) in the form
R(c)=c+c2(1+c)+c24(1+c)3−1. |
Assume the two initial solutions are
c1=0.8,c2=0.5. |
We obtain the following residuals
R1(0.8)=0.0279,R2(0.5)=−0.3271. |
By the Ying Buzu algorithm, c can be calculated as
c=R2c1−R1c2R2−R1=0.0279×0.5+0.3271×0.80.0279+0.3271=0.7764. |
The exact solution of (0.10) is
c=0.7758. |
The 4th order Taylor series solution is
u(x)=0.7758+0.2192x2+0.0057x4. |
Figure 1 shows the Taylor series solutions, which approximately meet the requirement of the boundary condition at x=1.
Now we use the Ying Buzu algorithm by choosing two initial guesses
u1(0)=0.5,u2(0)=1, |
which lead to u1=0.6726 and u2=1.2550, respectively, see Figure 2 (a) and (b).
It is obvious that the terminal value at x=1 deviates from u(1)=1 for each guess, according to the Ying Buzu algorithm, the initial guess can be updated as
u3(0)=0.5×(1−1.2550)−1×(1−0.6726)(1−1.2550)−(1−0.6726)=0.7810. | (0.11) |
The shooting process using (0.11) results in
u3(1)=1.0058, |
which deviates the exact value of u(1)=1 with a relative error of 0.5%, see Figure 3.
We can continue the iteration process to obtain a higher accuracy by using two following two guesses u1(0)=0.5, u3(0)=0.7810:
u4(0)=0.5×(1−1.0058)−0.7810×(1−0.6726)(1−1.0058)−(1−0.6726)=0.7761. |
Using this updated initial value, the shooting process leads to the result
u(1)=1.0001, |
so the approximate u(0)=0.7761 has only a relative error of 0.01%.
The above solution process couples the numerical method, and the ancient method can also be solved independently.
We assume that solution is
u(x)=c+(1−c)x2. | (0.12) |
Equation (0.12) meets all boundary conditions.
The residual equation is
R(x)=d2udx2−u1+u. |
It is easy to find that
R(0)=2(1−c)−c1+c. |
We choose two guesses:
c1=0.5,c2=1. |
We obtain the following residuals
R1(0)=2(1−0.5)−0.51+0.5=23, |
R2(0)=2(1−1)−11+1=−12. |
The Ying Buzu algorithm leads to the updated result:
c=c2R1(0)−c1R2(0)R1(0)−R2(0)=23×1+12×0.523+12=0.7857. |
The relative error is 1.2%, and the process can continue if a higher accuracy is still needed.
The ancient Chinese algorithm provides a simple and straightforward tool to two-point boundary value problems arising in chemistry, and it can be used for fast insight into the solution property of a complex problem.
The authors declare that they have no conflicts of interest to this work.
[1] |
P. H. Leslie, Some further notes on the use of matrices in population mathematics, Biometrika, 35 (1948), 213–245. https://doi.org/10.2307/2332342 doi: 10.2307/2332342
![]() |
[2] |
P. H. Leslie, J. Gower, The properties of a stochastic model for the predator-prey type of interaction between two species, Biometrika, 47 (1960), 219–234. https://doi.org/10.2307/2333294 doi: 10.2307/2333294
![]() |
[3] |
S. Hsu, T. Huang, Global stability for a class of predator-prey systems, SIAM J. Appl. Math., 55 (1995), 763–783. https://doi.org/10.1137/S0036139993253201 doi: 10.1137/S0036139993253201
![]() |
[4] |
Y. J. Li, M. X. He, Z. Li, Dynamics of a ratio-dependent Leslie-Gower predator-prey model with Allee effect and fear effect, Math. Comput. Simulation, 201 (2022), 417–439. https://doi.org/10.1016/j.matcom.2022.05.017 doi: 10.1016/j.matcom.2022.05.017
![]() |
[5] |
M. He, Z. Li, Dynamics of a Leslie-Gower predator-pery model with square root response function and generalist predator, Appl. Math. Lett., 157 (2024), 109193. https://doi.org/10.1016/j.aml.2024.109193 doi: 10.1016/j.aml.2024.109193
![]() |
[6] |
X. Wang, L. Zanette, X. Zou, Modelling the fear effect in predator-prey interactions, J. Math. Biol., 73 (2016), 1179–1204. https://doi.org/10.1007/s00285-016-0989-1 doi: 10.1007/s00285-016-0989-1
![]() |
[7] |
S. K. Sasmal, Population dynamics with multiple Allee effects induced by fear facotrs-A mathematical study on prey-predator interactions, Appl. Math. Model., 64 (2018), 1–14. https://doi.org/10.1016/j.apm.2018.07.021 doi: 10.1016/j.apm.2018.07.021
![]() |
[8] |
H. Zhang, Y. Cai, S. Fu, W. Wang, Impact of the fear effect in a prey-predator model incorporating a prey refuge, Appl. Math. Comput., 356 (2019), 328–337. https://doi.org/10.1016/j.amc.2019.03.034 doi: 10.1016/j.amc.2019.03.034
![]() |
[9] |
Y. Xue, Impact of both-density-dependent fear effect in a Leslie-Gower predator-prey model with Beddington-DeAngelis functional response, Chaos Soliton. Fract., 185 (2024), 115055. https://doi.org/10.1016/j.chaos.2024.115055 doi: 10.1016/j.chaos.2024.115055
![]() |
[10] |
R. K. Ghaziani, J. Alidousti, A. B. Eshkaftaki, Stability and dynamics of a fractional order Leslie-Gower prey-predator model, Appl. Math. Model., 40 (2016), 2075–2086. https://doi.org/10.1016/j.apm.2015.09.014 doi: 10.1016/j.apm.2015.09.014
![]() |
[11] |
C. Maji, Impact of fear effect in a fractional-order predator-prey system incorporating constant prey refuge, Nonlinear Dyn., 107 (2022), 1329–1342. https://doi.org/10.1007/s11071-021-07031-9 doi: 10.1007/s11071-021-07031-9
![]() |
[12] |
G. R. Kumar, K. Ramesh, A. Khan, K. Lakshminarayan, T. Abdeljawad, Dynamical study of fractional order Leslie-Gower model of predator-prey with fear, Allee effect, and inter-species rivalry, Res. Control Optim., 14 (2024), 100403. https://doi.org/10.1016/j.rico.2024.100403 doi: 10.1016/j.rico.2024.100403
![]() |
[13] |
H. Zhang, A. Muhammadhaji, Dynamics of a delayed fractional-order predator-prey model with cannibalism and disease in prey, Fractal Fract., 8 (2024), 333. https://doi.org/10.3390/fractalfract8060333 doi: 10.3390/fractalfract8060333
![]() |
[14] |
M. Awadalla, J. Alahmadi, K. R. Cheneke, S. Qureshi, Fractional optimal control model and bifurcation analysis of human syncytial respiratory virus transmission dynamics, Fractal Fract., 8 (2024), 44. https://doi.org/10.3390/fractalfract8010044 doi: 10.3390/fractalfract8010044
![]() |
[15] | I. Podlubny, Fractional differential equations, London: Academic Peress, 1999. |
[16] |
Y. Shi, Y. Q. Ma, X. Ding, Dynamical behaviors in a discrete fractional-order predator-prey system, Filomat, 32 (2018), 5857–5874. https://doi.org/10.2298/FIL1817857S doi: 10.2298/FIL1817857S
![]() |
[17] |
B. Wang, X. Li, Modeling and dynamical analysis of a fractional-order predator-prey system with anti-predator behavior and a Holling type Ⅳ functional response, Fractal Fract., 7 (2023), 722. https://doi.org/10.3390/fractalfract7100722 doi: 10.3390/fractalfract7100722
![]() |
[18] |
A. Singh, V. S. Sharma, Bifurcations and chaos control in a discrete-time prey-predator model with Holling type-Ⅱ functional response and prey refuge, J. Comput. Appl. Math., 418 (2023), 114666. https://doi.org/10.1016/j.cam.2022.114666 doi: 10.1016/j.cam.2022.114666
![]() |
[19] |
M. Berkal, M. B. Almatrafi, Bifurcation and stability of two-dimensional activator-inhibitor model with fractional-order derivative, Fractal Fract., 7 (2023), 344. https://doi.org/10.3390/fractalfract7050344 doi: 10.3390/fractalfract7050344
![]() |
[20] |
R. Saadeh, A. Abbes, A. Al-Husban, A. Ouannas, G. Grassi, The fractional discrete predator-prey model: Chaos, control and synchronization, Fractal Fract., 7 (2023), 120. https://doi.org/10.3390/fractalfract7020120 doi: 10.3390/fractalfract7020120
![]() |
[21] |
Q. Din, R. A. Naseem, M. S. Shabbir, Predator-prey interaction with fear effects: Stability, bifurcation and two-parameter analysis incorporating complex and fractal behavior, Fractal Fract., 8 (2024), 221. https://doi.org/10.3390/fractalfract8040221 doi: 10.3390/fractalfract8040221
![]() |
[22] |
Q. Din, Complexity and chaos control in a discrete-time prey-predator model, Commun. Nonlinear Sci., 49 (2017), 113–134. https://doi.org/10.1016/j.cnsns.2017.01.025 doi: 10.1016/j.cnsns.2017.01.025
![]() |
[23] |
Q. Din, Bifurcation analysis and chaos control in discrete-time glycolysis models, J. Math. Chem., 56 (2018), 904–931. https://doi.org/10.1007/s10910-017-0839-4 doi: 10.1007/s10910-017-0839-4
![]() |
[24] |
Q. Din, W. Ishaque, M.A. Iqbal, U. Saeed, Modification of Nicholson-Bailey model under refuge effects with stability, bifurcation, and chaos control, J. Vib. Control, 28 (2022), 3524–3538. https://doi.org/10.1177/10775463211034021 doi: 10.1177/10775463211034021
![]() |
[25] |
W. Ishaque, Q. Din, K. Khan, R. Mabela, Dynamics of predator-prey model based on fear effect with bifurcation analysis and chaos control, Qual. Theory Dyn. Syst., 23 (2024), 26. https://doi.org/10.1007/s12346-023-00878-w doi: 10.1007/s12346-023-00878-w
![]() |
[26] | S. Wiggins, Introduction to applied nonlinear dynamical systems and chaos, 2 Eds., New York: Springer New York, 2003. https://doi.org/10.1007/b97481 |
[27] |
H. Fan, J. Tang, K. Shi, Y. Zhao, Hybrid impulsive feedback control for drive-response synchronization of fractional-order multi-link Memristive neural networks with multi-delays, Fractal Fract., 7 (2023), 495. https://doi.org/10.3390/fractalfract7070495 doi: 10.3390/fractalfract7070495
![]() |
[28] |
K. Ding, Q. Zhu, Intermittent static output feedback control for stochastic delayed-switched positive systems with only partially measurable information, IEEE Trans. Autom. Control, 68 (2023), 8150–8157. https://doi.org/10.1109/TAC.2023.3293012 doi: 10.1109/TAC.2023.3293012
![]() |
1. | Yi Tian, Approximate solution of initial boundary value problems for ordinary differential equations with fractal derivative, 2022, 2, 2767-8946, 75, 10.3934/mmc.2022009 | |
2. | Guang-Qing Feng, Jing-Yan Niu, The analysis for the dynamic pull-in of a micro-electromechanical system, 2022, 1461-3484, 146134842211455, 10.1177/14613484221145588 | |
3. | YONG-GANG KANG, SHUAI-JIA KOU, SI-REN SONG, YU-ZHEN CHANG, AN-YANG WANG, YONG-GANG CHEN, VIBRATION ANALYSIS OF HEAVY WEAPONS IN TRANSIT BY AIRCRAFT IN FRACTAL SPACE CONSIDERING LOCATION DEVIATION, 2024, 32, 0218-348X, 10.1142/S0218348X2450018X | |
4. | Jing-Yan Niu, Guang-Qing Feng, A mini-review on ancient mathematics’ modern applications with an emphasis on the old Babylonian mathematics for MEMS systems, 2024, 12, 2296-424X, 10.3389/fphy.2024.1532630 |