Small amplitude periodic solution of Hopf Bifurcation Theorem for fractional differential equations of balance point in group competitive martial arts

Shuang Chen; Yuanjin Ren

Détails du journal et du numéro

Applied Mathematics and Nonlinear Sciences

Détails du magazine

Format: Magazine
eISSN: 2444-8656
Première parution: 01 Jan 2016
Périodicité: 2 fois par an
Langues: Anglais

Differential equation modelling was earlier used to discover better and understand various biological phenomena and social problems. We hope to understand the stability of the system and the Hopf bifurcation based on the characteristic roots of the linear system. Because group competitive sports require participants to have certain competitive skills, those who do not have sports skills but want to develop into activities must receive training and specific training. Therefore, based on the research background, the article proposes a time-lag group competitive martial arts activity model with a time lag effect. Through delay differential equation theory and Hopf bifurcation theory, the stability of the equilibrium point and the existence of periodic solutions generated by the Hopf bifurcation caused by the ‘instability’ of the equilibrium point are discussed. Finally, the theoretical results are simulated and verified with the help of MATLAB software.

Keywords

group competitive martial arts activities
balance point
stability
time lag
Hopf bifurcation
numerical simulation

MSC 2010

34A08

Article

Introduction

Differential equation modelling can be traced back to the early population model of Malthus and the predator–prey model of Lotka and Volterra. These models were once used to discover and understand various biological phenomena and social problems [1]. However, an overly simple model is difficult to accurately reflect the observed complex dynamic behaviours, such as periodic solutions. To this end, the model needs to be continuously improved. One method increases the equation's dimensionality continuously, but the cost is that it is difficult to estimate the parameters multiplied by actual data. Another way is to consider the time lag effect. The time lag effect is common in real problems. Time lag can correspond to the incubation period, delivery delay and response delay of the disease.

Moreover, simple time-delay differential systems often contain a wealth of complex dynamic behaviours. Some scholars have studied the gene regulation model of the time lag effect. Some scholars have conducted Hopf bifurcation research on Volterra predator–predator system with time delay [2]. Until the 1960s, the research on delayed differential systems mainly focused on stability, boundedness, asymptotics and equilibrium, periodic solutions, and oscillations of almost periodic solutions. Compared with ordinary differential systems, there is relatively little research on branch theory.

Hall was the first to study the local bifurcation of time-delay differential systems. He studied the existence of the central flow pattern of time-delay differential systems and the Hopf bifurcation theorem. However, Hall's theory is difficult to apply to practical problems. For time-delay differential systems with finite time delays, we hope to understand the stability of the system and the Hopf bifurcation based on the characteristic roots of the linear system. Since the characteristic equation of a linear system is a function of time delay, the characteristic root is also a function of time delay [3].

Moreover, the stability of the singularity will change as the time delay changes so that Hopf branches will occur near some critical values. Part of the time-delay differential equations will alternately appear from stable to unstable and then to stable as the singularity of the time-delay changes. This phenomenon is the so-called stable switching phenomenon.

Some scholars have proposed an ordinary differential equation model for group competitive sports activities [4]. Competitive activities here mean that participants must have certain skills rather than professional competitive sports activities. The article divides the human population into three categories. He set the total number of three types of people to 1, then the ordinary differential equation model of group competitive sports activities is as follows: ${\begin{array}{l} \frac{d p_{1} (t)}{d t} = (r_{1} - b) p_{1} (t) - α p_{1} (t) p_{2} (t) + a \\ \frac{d p_{2} (t)}{d t} = (r_{2} - b - β) p_{2} (t) + α p_{1} (t) p_{2} (t) \end{array}$ \left\{ {\matrix{ {{{d{p_1}(t)} \over {dt}} = ({r_1} - b){p_1}(t) - \alpha {p_1}(t){p_2}(t) + a} \hfill \cr {{{d{p_2}(t)} \over {dt}} = ({r_2} - b - \beta ){p_2}(t) + \alpha {p_1}(t){p_2}(t)} \hfill \cr } } \right. where α represents the conversion rate of category I individuals into category II individuals under the influence of category II individuals. β represents the conversion rate of the II individual into the III individual. It is assumed here that the III individual has no influence. a, b stands for birth rate and death rate, respectively. r_i(i = 1, 2), respectively, represents the migration rate of the I, II individual.

Because group competitive sports require participants to have certain competitive skills, those who do not have sports skills but want to develop into activities must receive training and specific training. It must take a certain amount of time [5]. The above model does not consider the time lag effect. This paper discusses the stability of the equilibrium point and the existence of periodic solutions generated by the Hopf bifurcation caused by the ‘instability’ of the equilibrium point through the theory of delay differential equations and the Hopf bifurcation theory. Finally, the theoretical results are simulated and verified with the help of MATLAB software.

Improved model

Considering that it takes a certain amount of time for an individual who participates in sports skills to transform into an individual who can participate in sports competitions, it is set as in Ref. [6]. Then the improved delay differential equation model can be expressed as: (1) ${\begin{array}{l} \frac{d x}{d t} = (r_{1} - b) x (t) - α x (t) y (t) + a \\ \frac{d y}{d t} = (r_{2} - b - β) y (t) + α x (t - τ) y (t - τ) \end{array}$ \left\{ {\matrix{ {{{dx} \over {dt}} = ({r_1} - b)x(t) - \alpha x(t)y(t) + a} \hfill \cr {{{dy} \over {dt}} = ({r_2} - b - \beta )y(t) + \alpha x(t - \tau )y(t - \tau )} \hfill \cr } } \right. where x, y, z represents the total number of individuals of class I, class II, and class III at time t; α represents the conversion rate of class I individuals into class II individuals under the influence of class II individuals; and β represents the conversion rate of the II individual into the III type individual. It is assumed here that the category III individuals have no influence. a, b represents the birth rate and death rate, respectively. r_i(i = 1, 2, 3) respectively represents the migration rate of the I, II, III individual [7]. Since the sum of the total number of people in the three groups remains unchanged, only the above-mentioned two-dimensional system composed of I, II type individuals must be considered.

The stability of the equilibrium point and the existence of the Hopf branch

First calculate the balance point of the system. We make the right end of the system (1) zero, that is, $\frac{d x}{d t} = 0$ {{dx} \over {dt}} = 0 , $\frac{d y}{d t} = 0$ {{dy} \over {dt}} = 0 . Then the balance points are $E_{0} (x_{0}, 0) = (\frac{a}{b - r_{1}}, 0)$ {E_0}({x_0},0) = \left( {{a \over {b - {r_1}}},0} \right) and E_*(x^*, y^*), where $x^{*} = \frac{b - r_{2} + β}{a}$ {x^*} = {{b - {r_2} + \beta } \over a} , $y^{*} = \frac{α a + (r_{1} - b) (b - r_{2} + β)}{a (b - r_{2} + β)}$ {y^*} = {{\alpha a + ({r_1} - b)(b - {r_2} + \beta )} \over {a(b - {r_2} + \beta )}} .

3.1

Asymptotic stability of boundary equilibrium point E₀

Boundary equilibrium point $E_{0} (x_{0}, 0) = (\frac{a}{b - r_{1}}, 0)$ {E_0}({x_0},0) = \left( {{a \over {b - {r_1}}},0} \right) , because the number of people in each type of group in competitive activities needs to be a positive value [8]. Assuming that system (1) satisfies H₁: b − r₁ > 0, the boundary equilibrium point is a non-negative equilibrium point: (2) ${\begin{array}{l} \frac{d x}{d t} = (r_{1} - b) x (t) - \frac{α a}{b - r_{1}} y (t) \\ \frac{d y}{d t} = (r_{2} - b - β) y (t) + \frac{α a}{b - r_{1}} y (t - τ) \end{array}$ \left\{ {\matrix{ {{{dx} \over {dt}} = \left( {{r_1} - b} \right)x(t) - {{\alpha a} \over {b - {r_1}}}y(t)} \hfill \cr {{{dy} \over {dt}} = \left( {{r_2} - b - \beta } \right)y(t) + {{\alpha a} \over {b - {r_1}}}y(t - \tau )} \hfill \cr } } \right. For the convenience of calculation, we make a series of variable substitutions. Assuming m₁ = r₁ − b, m₂ = r₂ − b − β, $m_{3} = \frac{α a}{b - r_{1}}$ {m_3} = {{\alpha a} \over {b - {r_1}}} , then formula (2) can be expressed as: (3) ${\begin{array}{l} \frac{d x}{d t} = m_{1} x (t) - m_{3} y (t) \\ \frac{d y}{d t} = m_{2} y (t) + m_{3} y (t - τ) \end{array}$ \left\{ {\matrix{ {{{dx} \over {dt}} = {m_1}x(t) - {m_3}y(t)} \hfill \cr {{{dy} \over {dt}} = {m_2}y(t) + {m_3}y(t - \tau )} \hfill \cr } } \right. The characteristic equation corresponding to system (3) at equilibrium point $E_{0} (x_{0}, 0) = (\frac{a}{b - r_{1}}, 0)$ {E_0}({x_0},0) = \left( {{a \over {b - {r_1}}},0} \right) is: (4) $(λ - m_{1}) (λ - m_{2} - m_{3} e^{- λ τ}) = 0$ \left( {\lambda - {m_1}} \right)\left( {\lambda - {m_2} - {m_3}{e^{ - \lambda \tau }}} \right) = 0 According to the hypothesis H₁, λ = m₁ = r₁ − b is always the negative root of Eq. (4). So let us study the second factor λ − m₂ − m₃e^−λτ of the characteristic Eq. (4).

Hypothesis (5) $G (λ) = λ - m_{2} - m_{3} e^{- λ τ}$ G(\lambda ) = \lambda - {m_2} - {m_3}{e^{ - \lambda \tau }} If m₃ > −m₂ is true, you can get: $G (0) = - m_{2} - m_{3} < 0, lim_{λ \to + \infty} G (λ) = + \infty$ G(0) = - {m_2} - {m_3} < 0,\;\mathop {\lim }\limits_{\lambda \to + \infty } G(\lambda ) = + \infty . Therefore G(λ) = 0 has at least one positive real root. When $m_{2} + m_{3} = r_{2} - b - β + \frac{α a}{b - r_{1}} > 0$ {m_2} + {m_3} = {r_2} - b - \beta + {{\alpha a} \over {b - {r_1}}} > 0 , the boundary equilibrium point E₀(x₀, 0) is unstable. If m₃ < −m₂ is true, when τ = 0, we know that the system is asymptotically stable at the boundary equilibrium point. Consider the case of τ > 0 below.

Lemma 1

When τ > 0, if m₃ < −m₂ holds, the characteristic Eq. (4) has no pure imaginary roots.

Proof

Suppose λ = iω(ω > 0) is a pure imaginary root of the characteristic Eq. (4). We substitute λ = iω into Eq. (4) and separate the real and imaginary parts from obtaining the following equation: (6) ${\begin{array}{l} m_{2} = - m_{3} cos ω τ \\ ω = - m_{3} sin ω τ \end{array}$ \left\{ {\matrix{ {{m_2} = - {m_3}\cos \omega \tau } \hfill \cr {\omega = - {m_3}\sin \omega \tau }\hfill \cr } } \right. Add the squares to get $ω^{2} = m_{3}^{2} - m_{2}^{2} < 0$ {\omega ^2} = m_3^2 - m_2^2 < 0 This contradicts the assumption. There is no pure imaginary root for any τ ≥ 0 Eq. (4). In other words, the zero points of the characteristic equation will not appear on the imaginary axis or cross the imaginary axis [9]. Based on the above analysis, when the system satisfies the assumption H₁, the conclusions that can be obtained for any τ ≥ 0 are as follows: a) If m₃ > −m₂, E₀(x₀, 0) is unstable; b) If m₃ < −m₂, E₀(x₀, 0) is locally asymptotically stable.

3.2

Asymptotic stability of the positive equilibrium point E_*(x^*, y^*)

Lemma 2

Assuming that H₂: b − r₂ + β > 0, αa + (r₁ − b)(b − r₂ + β) > 0 holds, then system (1) has the only positive equilibrium point E_*(x^*, y^*): System (1) The linearised system at the positive balance point E_* (x^*, y^*) is: x′(t) = Ax(t) + Bx(t − τ). Among them: $A = (\begin{matrix} r_{1} - b - a y^{*} & - x y^{*} \\ 0 & r_{2} - b - β \end{matrix})$ A = \left( {\matrix{ {{r_1} - b - a{y^*}} & { - x{y^*}} \cr 0 & {{r_2} - b - \beta } \cr } } \right) , $B = (\begin{matrix} 0 & 0 \\ a y^{*} & a x^{*} \end{matrix})$ B = \left( {\matrix{ 0 & 0 \cr {a{y^*}} & {a{x^*}} \cr } } \right) .

According to the above formula, the linearisation system at the positive equilibrium point E_* (x^*, y^*) can be obtained as: (7) ${\begin{array}{l} \frac{d x}{d t} = (r_{1} - b - a y^{*}) x (t) - a x^{*} y (t) \\ \frac{d y}{d t} = a y^{*} x (t - τ) + (r_{2} - b - β) \\ y (t) + a x^{*} y (t - τ) \end{array}$ \left\{ {\matrix{ {{{dx} \over {dt}} = ({r_1} - b - a{y^*})x(t) - a{x^*}y(t)} \hfill \cr {{{dy} \over {dt}} = a{y^*}x(t - \tau ) + ({r_2} - b - \beta )}\hfill \cr {y(t) + a{x^*}y(t - \tau )} \hfill \cr } } \right. The characteristic equations corresponding to the system (7) can be sorted out: (8) $λ^{2} + (p_{0} λ + p_{1}) e^{- λ τ} + q_{0} λ + q_{1} = 0$ {\lambda ^2} + \left( {{p_0}\lambda + {p_1}} \right){e^{ - \lambda \tau }} + {q_0}\lambda + {q_1} = 0 Among them: $\begin{matrix} p_{0} = - a x^{*}, q_{0} = (b - r_{1} + a y^{*}) - (r_{2} - b - β), \\ p_{1} = a x^{*} (r_{1} - b), q_{1} = (r_{1} - b - a y^{*}) (r_{2} - b - β) \end{matrix}$ \matrix{ {{p_0} = - a{x^*},\quad {q_0} = \left( {b - {r_1} + a{y^*}} \right) - \left( {{r_2} - b - \beta } \right),} \cr {{p_1} = a{x^*}\left( {{r_1} - b} \right),\quad {q_1} = \left( {{r_1} - b - a{y^*}} \right)\left( {{r_2} - b - \beta } \right)} \cr } .

Lemma 3

Assume that condition H₂ holds and there is H₃: q₁ − p₁ = αa − (r₁ − b)(b − r₂ + β) < 0. Then there is τ₀ when τ increases from zero.

Proof

First, prove that the characteristic equation has a pair of purely imaginary roots ±iω₀, ω₀ > 0 at the positive equilibrium point E_*.

When τ = 0 is the formula (8), we can get: (9) $λ^{2} + (p_{0} + q_{0}) λ + (p_{1} + q_{1}) = 0$ {\lambda ^2} + \left( {{p_0} + {q_0}} \right)\lambda + \left( {{p_1} + {q_1}} \right) = 0 For the above condition without time delay, it is known that when τ = 0 is the positive equilibrium point E_* is asymptotically stable [10]. The case of τ > 0 is discussed below. When τ > 0, if λ = iω(ω > 0) is the pure imaginary root of the characteristic Eq. (8). We substitute it and separate the real part and the imaginary part to get: (10) ${\begin{array}{l} ω^{2} - q_{1} = p_{0} ω sin ω τ + p_{1} cos ω τ \\ q_{0} ω = p_{1} sin ω τ - p_{0} ω cos ω τ \end{array}$ \left\{ {\matrix{ {{\omega ^2} - {q_1} = {p_0}\omega \sin \omega \tau + {p_1}\cos \omega \tau } \hfill \cr {{q_0}\omega = {p_1}\sin \omega \tau - {p_0}\omega \cos \omega \tau } \hfill \cr } } \right. Eliminate the trigonometric function from Eq. (10) to get: (11) $ω^{4} + (q_{0}^{2} - p_{0}^{2} - 2 q_{1}) ω^{2} + (q_{1}^{2} - p_{1}^{2}) = 0$ {\omega ^4} + \left( {q_0^2 - p_0^2 - 2{q_1}} \right){\omega ^2} + \left( {q_1^2 - p_1^2} \right) = 0 In: (12) $q_{0}^{2} - p_{0}^{2} - 2 q_{1} = {(r_{1} - b - a y^{*})}^{2} > 0$ q_0^2 - p_0^2 - 2{q_1} = {\left( {{r_1} - b - a{y^*}} \right)^2} > 0 If another H₃: q₁ − p₁ = αa − (r₁ − b)(b − r₂ + β) < 0 is established. At the same time, it is known from (12) that Eq. (11) has a unique positive real root ω₀, then Eq. (8) has a pair of purely imaginary roots ±iω₀ that are conjugate to each other when τ = τ_j, j = 0, 1, 2,··· is.

According to Eq. (10), we can get: (13) $ω_{0} = \sqrt{\frac{1}{2} [(p_{0}^{2} - q_{0}^{2} + 2 q_{1}) + \sqrt{{(q_{0}^{2} - p_{0}^{2} - 2 q_{1})}^{2} - 4 (q_{1}^{2} - p_{1}^{2})}]}$ {\omega _0} = \sqrt {{1 \over 2}\left[ {(p_0^2 - q_0^2 + 2{q_1}) + \sqrt {{{(q_0^2 - p_0^2 - 2{q_1})}^2} - 4(q_1^2 - p_1^2)} } \right]} Combining (10)–(12), define (14) $τ_{j} = \frac{1}{ω_{0}} (arccos \frac{φ}{φ} + 2 j π), j = 0, 1, 2, \dots$ {\tau _j} = {1 \over {{\omega _0}}}\left( {\arccos {\varphi \over \varphi } + 2j\pi } \right),\quad j = 0,1,2, \cdots where $φ = (p_{1} - p_{0} q_{0}) ω_{0}^{2} - p_{1} q_{1}$ \varphi = \left( {{p_1} - {p_0}{q_0}} \right)\omega _0^2 - {p_1}{q_1} , $\tilde{ω} = p_{0}^{2} ω_{0}^{2} + p_{1}^{2}$ \widetilde \omega = p_0^2\omega _0^2 + p_1^2 . Then (τ, λ) = (τ_j, iω₀) satisfies Eq. (8), with the definition (15) $λ (τ) = a (τ) + i ω (τ)$ \lambda (\tau ) = a(\tau ) + i\omega (\tau ) It is only necessary to prove that the transversal condition $(\frac{d R e λ (τ)}{d (τ)} |_{τ = τ_{j}}) > 0$ \left( {{{dRe\lambda (\tau )} \over {d(\tau )}}\left| {_{\tau = {\tau _j}}} \right.} \right) > 0 holds. When τ > τ_j, there is at least one characteristic root whose real part is >0 so that the condition of Hopf bifurcation is satisfied and a periodic solution is generated near E_*. Substituting λ = λ (τ) into Eq. (8) and deriving the derivation of both sides concerning τ at the same time: (16) ${(\frac{d λ}{d τ})}^{- 1} = \frac{2 λ + p_{0} e^{- λ τ} + q_{0}}{p_{0} λ^{2} e^{- λ τ} + p_{1} λ e^{- λ τ}} - \frac{τ}{λ}$ {\left( {{{d\lambda } \over {d\tau }}} \right)^{ - 1}} = {{2\lambda + {p_0}{e^{ - \lambda \tau }} + {q_0}} \over {{p_0}{\lambda ^2}{e^{ - \lambda \tau }} + {p_1}\lambda {e^{ - \lambda \tau }}}} - {\tau \over \lambda } Hence ${(\frac{dRe λ (τ)}{d τ} |_{τ = τ_{j}})}^{- 1} = Re {(\frac{2 λ + p_{0} e^{- λ τ} + q_{0}}{p_{0} λ^{2} e^{- λ τ} + p_{1} e^{- λ τ}})}_{τ = τ_{j}} = \frac{M X + N Y}{X^{2} + Y^{2}}$ {\left( {{{dRe\lambda (\tau )} \over {d\tau }}\left| {_{\tau = {\tau _j}}} \right.} \right)^{ - 1}} = \mathop {{\rm{Re}}}{\left( {{{2\lambda + {p_0}{e^{ - \lambda \tau }} + {q_0}} \over {{p_0}{\lambda ^2}{e^{ - \lambda \tau }} + {p_1}{e^{ - \lambda \tau }}}}} \right)_{\tau = {\tau _j}}} = {{MX + NY} \over {{X^2} + {Y^2}}} . Its M = p₀ cos ωτ, N = 2ω − p₀ sin ωτ, X = −p₀ cos ωτ + p₁ω sin ωτ, Y = p₀ω² sin ωτ + p₁ω cos ωτ. According to formula (13), the following formula can be deduced to hold $\frac{M X + N Y}{X^{2} + Y^{2}} = \frac{(2 ω^{2} - 2 q_{1} - p_{0}^{2})}{(p_{0}^{2} ω^{2} + p_{1}^{2})} > 0$ {{MX + NY} \over {{X^2} + {Y^2}}} = {{(2{\omega ^2} - 2{q_1} - p_0^2)} \over {(p_0^2{\omega ^2} + p_1^2)}} > 0 . This proves that the transversal condition holds. Therefore the Hopf branch occurs at ω = ω₀, τ = τ_j.

3.3

Hopf branch of the system

The stability of the time-delay system at τ = τ_j changes. Therefore, according to the Hopf bifurcation theory, the time-lag group competitive sports activity model (1) has a Hopf bifurcation at its positive balance point E_*(x^*, y^*) [11].

Numerical simulation

We select a set of parameters as follows: $a = 0.012, a = 0.38, r_{1} - b = 0.21, r_{2} - b - β = - 0.126$ {a = 0.012,\quad a = 0.38,\quad {r_1} - b = 0.21,\quad {r_2} - b - \beta = - 0.126} According to the previous analysis and formula, it is calculated that $\begin{matrix} p_{0} = - 0.1260, q_{0} = 0.1622, p_{1} = 0.0265, q_{1} = 0.0046 \\ E_{*} (x^{*}, y^{*}) = (0.3316, 0.6479), ω^{2} = 0.1594; ω = 0.0254; τ_{0} = 1.5222 \end{matrix}$ \matrix{ {{p_0} = - 0.1260,\quad {q_0} = 0.1622,\quad {p_1} = 0.0265,\quad {q_1} = 0.0046} \cr {{E_*}({x^*},{y^*}) = (0.3316,0.6479),\quad {\omega ^2} = 0.1594;\quad \omega = 0.0254;\quad {\tau _0} = 1.5222} \cr } Drawing through MATLAB software, System (1) obtains the numerical simulation result at time lag τ = 1.3 > τ₀ and obtains the waveform diagram and the plane phase diagram of the numerical solution of the system (Figures 1 and 2). The numerical simulation results at time lag τ = 1.56 > τ₀ are the same as above (see Figures 3 and 4).

MATLAB numerical simulation shows that the stability of the positive equilibrium has changed with the change of time delay [12].

Waveform graph of the percentage of the I and II population in the total population versus time

Plane phase diagrams of I and II groups of people

Waveform diagram of the percentages of the I and II groups in the total number of people over time

Planar phase diagram of the number of people in category I and category II

Conclusion

There are certain critical points in the process of the time delay increasing from 0. At these critical points, the stability of the system equilibrium point always changes from stable to unstable; from zero to increase every time the critical value is passed from the positive equilibrium point. Expenditure periodic solution. Second, a numerical simulation of the time-lag group competitive sports activity model was carried out by MATLAB software. The theoretical results of the above model were verified by selecting a set of parameters.

Figures et tableaux

Références

[1] Mboupda Pone, J. R., Kingni, S. T., Kol, G. R., & Pham, V. T. Hopf bifurcation, antimonotonicity and amplitude controls in the chaotic Toda jerk oscillator: analysis, circuit realization and combination synchronization in its fractional-order form. Automatika., 2019; 60(2): 149–161. Mboupda PoneJ. R. KingniS. T. KolG. R. PhamV. T. Hopf bifurcation, antimonotonicity and amplitude controls in the chaotic Toda jerk oscillator: analysis, circuit realization and combination synchronization in its fractional-order form Automatika 2019 60 2 149 161 10.1080/00051144.2019.1600109 Search in Google Scholar

[2] El-Saka, H. A. A., Lee, S., & Jang, B. Dynamic analysis of fractional-order predator–prey biological economic system with Holling type II functional response. Nonlinear Dynamics., 2019; 96(1): 407–416. El-SakaH. A. A. LeeS. JangB. Dynamic analysis of fractional-order predator–prey biological economic system with Holling type II functional response Nonlinear Dynamics 2019 96 1 407 416 10.1007/s11071-019-04796-y Search in Google Scholar

[3] Xu, C., & Aouiti, C. Comparative analysis on Hopf bifurcation of integer-order and fractional-order two-neuron neural networks with delay. International Journal of Circuit Theory and Applications., 2020; 48(9): 1459–1475. XuC. AouitiC. Comparative analysis on Hopf bifurcation of integer-order and fractional-order two-neuron neural networks with delay International Journal of Circuit Theory and Applications 2020 48 9 1459 1475 10.1002/cta.2847 Search in Google Scholar

[4] Shi, M., Yu, Y., & Xu, Q. Hopf bifurcation of the fractional-order Hindmarsh–Rose neuron model with time-delay. Rocky Mountain Journal of Mathematics., 2020; 50(6): 2213–2222. ShiM. YuY. XuQ. Hopf bifurcation of the fractional-order Hindmarsh–Rose neuron model with time-delay Rocky Mountain Journal of Mathematics 2020 50 6 2213 2222 10.1216/rmj.2020.50.2213 Search in Google Scholar

[5] Tang, B. Dynamics for a fractional-order predator-prey model with group defense. Scientific reports., 2020; 10(1): 1–17. TangB. Dynamics for a fractional-order predator-prey model with group defense Scientific reports 2020 10 1 1 17 10.1038/s41598-020-61468-3707823032184437 Search in Google Scholar

[6] Naim, M., Lahmidi, F., & Namir, A. Mathematical analysis of a fractional order SIS epidemic model with double diseases, Beddington-DeAngelis functional response and time delay. International Journal of Nonlinear Science., 2020;29(1): 47–59. NaimM. LahmidiF. NamirA. Mathematical analysis of a fractional order SIS epidemic model with double diseases, Beddington-DeAngelis functional response and time delay International Journal of Nonlinear Science 2020 29 1 47 59 Search in Google Scholar

[7] Wang, X., Wang, Z., & Xia, J. Stability and bifurcation control of a delayed fractional-order eco-epidemiological model with incommensurate orders. Journal of the Franklin Institute., 2019; 356(15): 8278–8295. WangX. WangZ. XiaJ. Stability and bifurcation control of a delayed fractional-order eco-epidemiological model with incommensurate orders Journal of the Franklin Institute 2019 356 15 8278 8295 10.1016/j.jfranklin.2019.07.028 Search in Google Scholar

[8] Shi, R., Ren, J., & Wang, C. Stability analysis and Hopf bifurcation of a fractional order mathematical model with time delay for nutrient-phytoplankton-zooplankton. Mathematical Biosciences and Engineering., 2020; 17(4): 3836–3868. ShiR. RenJ. WangC. Stability analysis and Hopf bifurcation of a fractional order mathematical model with time delay for nutrient-phytoplankton-zooplankton Mathematical Biosciences and Engineering 2020 17 4 3836 3868 10.3934/mbe.202021432987557 Search in Google Scholar

[9] Xu, C., Liao, M., & Li, P. Bifurcation control of a fractional-order delayed competition and cooperation model of two enterprises. Science China Technological Sciences., 2019; 62(12): 2130–2143. XuC. LiaoM. LiP. Bifurcation control of a fractional-order delayed competition and cooperation model of two enterprises Science China Technological Sciences 2019 62 12 2130 2143 10.1007/s11431-018-9376-2 Search in Google Scholar

[10] Sulaiman, T., Bulut, H. & Baskonus, H. On the exact solutions to some system of complex nonlinear models. Applied Mathematics and Nonlinear Sciences., 2021; 6(1): 29–42. SulaimanT. BulutH. BaskonusH. On the exact solutions to some system of complex nonlinear models Applied Mathematics and Nonlinear Sciences 2021 6 1 29 42 10.2478/amns.2020.2.00007 Search in Google Scholar

[11] Hu, X., Li, J. & Aram, Research on style control in planning and designing small towns. Applied Mathematics and Nonlinear Sciences., 2021; 6(1): 57–64. HuX. LiJ. Aram Research on style control in planning and designing small towns Applied Mathematics and Nonlinear Sciences 2021 6 1 57 64 10.2478/amns.2020.2.00077 Search in Google Scholar

[12] Moustafa, M., Mohd, M. H., Ismail, A. I., & Abdullah, F. A. Stage structure and refuge effects in the dynamical analysis of a fractional order Rosenzweig-MacArthur prey-predator model. Prog. Fract. Differ. Appl., 2019; 5(1): 49–64. MoustafaM. MohdM. H. IsmailA. I. AbdullahF. A. Stage structure and refuge effects in the dynamical analysis of a fractional order Rosenzweig-MacArthur prey-predator model Prog. Fract. Differ. Appl. 2019 5 1 49 64 10.18576/pfda/050106 Search in Google Scholar