Next Article in Journal
Fractional Order Mathematical Model of Serial Killing with Different Choices of Control Strategy
Next Article in Special Issue
Fixed-Time Fractional-Order Global Sliding Mode Control for Nonholonomic Mobile Robot Systems under External Disturbances
Previous Article in Journal
A Chebyshev Collocation Approach to Solve Fractional Fisher–Kolmogorov–Petrovskii–Piskunov Equation with Nonlocal Condition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Analysis and Bifurcation Study on Fractional-Order Tri-Neuron Neural Networks Incorporating Delays

1
School of Mathematics and Statistics, Henan University of Science and Technology, Luoyang 471023, China
2
Guizhou Key Laboratory of Economics System Simulation, Guizhou University of Finance and Economics, Guiyang 550004, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(3), 161; https://doi.org/10.3390/fractalfract6030161
Submission received: 7 February 2022 / Revised: 23 February 2022 / Accepted: 8 March 2022 / Published: 15 March 2022

Abstract

:
In this manuscript, we principally probe into a class of fractional-order tri-neuron neural networks incorporating delays. Making use of fixed point theorem, we prove the existence and uniqueness of solution to the fractional-order tri-neuron neural networks incorporating delays. By virtue of a suitable function, we prove the uniformly boundedness of the solution to the fractional-order tri-neuron neural networks incorporating delays. With the aid of the stability theory and bifurcation knowledge of fractional-order differential equation, a new delay-independent condition to guarantee the stability and creation of Hopf bifurcation of the fractional-order tri-neuron neural networks incorporating delays is established. Taking advantage of the mixed controller that contains state feedback and parameter perturbation, the stability region and the time of onset of Hopf bifurcation of the fractional-order trineuron neural networks incorporating delays are successfully controlled. Software simulation plots are displayed to illustrate the established key results. The obtained conclusions in this article have important theoretical significance in designing and controlling neural networks.

1. Introduction

After the classical work of Hopfield [1] on neural networks, many scholars pay much attention to the study of the dynamics of different types of neural networks since neural networks have displayed extraordinary application value in various fields such as associate memory, automatic control, image processing, biological engineering, pattern recognition and so on [2,3,4,5]. Generally speaking, delay often occurs in the signal transmission among different neurons due to the finite propagation velocity in neural networks [6]. A natural problem arises: What is the impact of time delay on the dynamical behavior of neural networks? During the past decades, in order to reveal the influence of time delay on various dynamics of neural networks, a great deal of researchers from mathematics and engineering have made great efforts to explore the dynamics of delayed neural networks and lots of valuable fruits have been achieved. For instance, Wang et al. [7] explored the fixed-time synchronization problem of delayed complex-valued BAM neural networks by using pinning control and adaptive pinning control. Gan et al. [8] discussed the anti-synchronization of a class of BAM neural networks involving the Markov scheduling protocol. Syed Ali et al. [9] set up a novel sufficient condition to ensure the global Mittag–Leffler stability of impulsive fractional-order delayed complex-valued BAM neural networks. Chen et al. [10] studied the finite-time stabilization of fractional-order fuzzy quaternion-valued BAM neural networks by applying direct quaternion approach. Making use of basic theories on fractional calculus, inequality skills of fuzzy logic and reduction to absurdity, Chen et al. [10] established the sufficient criteria that guarantee the finite-time stabilization of the fractional-order fuzzy quaternion-valued BAM neural networks. Cong et al. [11] dealt with the global exponential stability of discrete-time delayed BAM neural networks. By virtue of the mathematical induction method, a new exponential stability criterion for discrete-time delayed BAM neural networks is built. For more specific publications on this aspect, one can see [12,13,14,15,16].
Here we notice that most of the literature (e.g., [1,2,3,4,5,6,7,8,9,11,12,13,14,15,16]) only involves integer-order delayed neural networks, not considering the fractional-order cases. Since the introduction of fractional-order derivative by Leibniz and L’Hosptical [17] three hundred years ago, fractional calculus gradually arouse much attention from many scholars. The fractional derivative has more greater advantages than integer-order ones. It can adequately display the global correlation and better description of long-term development process and memory [18]. Recent research shows that fractional calculus has been widely used in many areas such as electrical engineering, wave theory in physics, electronic information, biological technique, control engineering, neural systems and so on [19,20,21]. Thus, the investigation of the dynamics of fractional-order dynamical systems has attracted much interest from numerous scholars and fruitful results have been achieved. For instance, Zhang et al. [22] explored the multistability of fractional-order competitive neural networks with delays. Li et al. [23] made a detailed analysis on the complete and finite-time synchronization of fractional-order fuzzy neural networks by virtue of nonlinear feedback control. Padmaja and Balasubramaniam [24] set up a new delay and order-dependent passivity criteria for impulsive fractional-order neural networks involving proportional delays. Ke [25] proded into the Mittag–Leffler stability and asymptotic ω -periodicity of fractional-order delayed inertial neural networks. In details, one can see [26,27,28,29,30,31,32,33,34,35,36,37].
In fractional-order differential equations, time delay is a vital factor that affects the dynamical behavior of systems. In particular, delay-induced Hopf bifurcation is an important topic in fractional-order differential systems. In order to design, optimize and control neural networks to serve humanity, we need to explore the impact of time delay on the dynamics of delayed fractional-order neural networks, especially, we need to investigate the impact of time delay on Hopf bifurcation phenomenon of delayed fractional-order neural networks. By adjusting the value of time delay, we can enlarge or narrow the stability region of neural network systems, then postpone or advance the time of onset of Hopf bifurcation of neural networks. However, many works on delay-induced Hopf bifurcation mainly focus on the integer-order delayed neural networks and there are few works on the fractional-order delayed neural networks. Recently, there has been some literature on delay-induced Hopf bifurcation of fractional-order delayed dynamical systems. For example, Xu et al. [38] established a sufficient condition to ensure the stability and the onset of Hopf bifurcation of fractional-order six-neuron BAM neural networks with multi-delays. Huang et al. [39] explored the Hopf bifurcation of fractional-order quaternion-valued neural networks. Xiao et al. [40] investigated the fractional-order PD control of Hopf bifurcations in fractional-order small-world networks with delays. For more related studies, see Refs. [41,42,43,44,45,46,47]. Although some works on Hopf bifurcation of fractional-order delayed neural networks have been carried out, many challenging problems are expected to be solved. Stimulated by the idea, in this present work, we are to deal with the following three aspects: (i) Explore the existence and uniqueness, the boundedness of solution of the involved neural networks with single delay; (ii) Set up a sufficient criterion ensuring the stability and and the creation of Hopf bifurcation of the involved neural networks with single delay; (iii) Make use of a suitable mixed controller that contains state feedback and parameter perturbation to control the time of onset of Hopf bifurcation for the involved neural networks with single delay.
In 2006, Yan [41] investigated the following tri-neuron network system:
w ˙ 1 ( t ) = ν w 1 ( t ) + h ( w 1 ( t ) ) + h 12 ( w 2 ( t ρ ) ) + h 13 ( w 3 ( t ρ ) ) , w ˙ 2 ( t ) = ν w 2 ( t ) + h ( w 2 ( t ) ) + h 21 ( w 1 ( t ρ ) ) + h 23 ( w 3 ( t ρ ) ) , w ˙ 3 ( t ) = ν w 3 ( t ) + h ( w 3 ( t ) ) + h 31 ( w 1 ( t ρ ) ) + h 32 ( w 2 ( t ρ ) ) ,
where w 1 ( t ) , w 2 ( t ) , w 3 ( t ) stand for the voltage on the input of the first neuron at time t ( t > 0 ) , the voltage on the input of the second neuron at time t ( t > 0 ) , the voltage on the input of the third neuron at time t ( t > 0 ) , respectively; ν > 0 stands for the rate with which the three neurons will reset its potential to the resetting state in isolation when disconnected from the network and external inputs. h , h 12 , h 13 , h 21 , h 23 , h 31 , h 32 are the sigmoidal activation functions of the different neurons. ρ > 0 represents the transmission delay of the signal along the axon of the neuron. In details, one can see [41,42,43]. By analyzing the characteristic equation of system (1) and selecting the time delay as bifurcation parameter, Yan [41] established the delay-independent stability and Hopf bifurcation condition ensuring the stability and the creation of Hopf bifurcation of system (1). In addition, making use of the normal form and center manifold theorem, Yan [41] explored the Hopf bifurcation peculiarities.
On the basis of the analysis above, we will deal with the delay-induced Hopf bifurcation of fractional order neural networks. Inspired by the above viewpoint, we modify model (1) as the following fractional-order form
d r w 1 ( t ) d t r = ν w 1 ( t ) + h ( w 1 ( t ) ) + h 12 ( w 2 ( t ρ ) ) + h 13 ( w 3 ( t ρ ) ) , d r w 2 ( t ) d t r = ν w 2 ( t ) + h ( w 2 ( t ) ) + h 21 ( w 1 ( t ρ ) ) + h 23 ( w 3 ( t ρ ) ) , d r w 3 ( t ) d t r = ν w 3 ( t ) + h ( w 3 ( t ) ) + h 31 ( w 1 ( t ρ ) ) + h 32 ( w 2 ( t ρ ) ) ,
where r ( 0 , 1 ] is a constant. We give the initial value of system (2) as follows:
w 1 ( ϕ ) = w 1 ϕ , ϕ [ ρ , t 0 ] , w 2 ( ϕ ) = w 2 ϕ , ϕ [ ρ , t 0 ] , w 3 ( ϕ ) = w 1 ϕ , ϕ [ ρ , t 0 ] ,
where t 0 > 0 is a constant. In order to achieve our goal in this paper, we prepare the following necessary assumptions:
  • ( Q 1 ) h , h 12 , h 13 , h 21 , h 23 , h 31 , h 32 C 1 , h ( 0 ) = h 12 ( 0 ) = h 13 ( 0 ) = h 21 ( 0 ) = h 23 ( 0 ) = h 31 ( 0 ) = h 32 ( 0 ) = 0 .
  • ( Q 2 ) There exist positive constants H , H 12 , H 13 , H 21 , H 23 , H 31 , H 32 such that
    | h ( t 1 ) h ( t 2 ) | H | t 1 t 2 | , | h 12 ( t 1 ) h 12 ( t 2 ) | H 12 | t 1 t 2 | , | h 13 ( t 1 ) h 13 ( t 2 ) | H 13 | t 1 t 2 | , | h 21 ( t 1 ) h 21 ( t 2 ) | H 21 | t 1 t 2 | , | h 23 ( t 1 ) h 23 ( t 2 ) | H 23 | t 1 t 2 | , | h 31 ( t 1 ) h 31 ( t 2 ) | H 31 | t 1 t 2 | , | h 32 ( t 1 ) h 32 ( t 2 ) | H 32 | t 1 t 2 | .
  • ( Q 3 ) There exist positive constants H , H 12 , H 13 , H 21 , H 23 , H 31 , H 32 such that
    | h ( t ) | H , | h 12 ( t ) | H 12 , | h 13 ( t ) | H 13 , | h 21 ( t ) | H 21 , | h 23 ( t ) | H 23 , | h 31 ( t ) | H 31 , | h 32 ( t ) | H 32 ,
    for all t R .
We arrange the basic structure of the article as follows. Section 2 lists some related basic theories on fractional-order dynamical system. Section 3 gives the proof existence and uniqueness of solution of system (2). Section 4 shows the proof of boundedness of solution of system (2). Section 5 establishes the delay-independent sufficient condition to ensure the stability and the creation of Hopf bifurcation of system (2). Section 6 seeks a suitable mixed controller to control the Hopf bifurcation of system (2). Section 7 carries out software simulations to support the derived key conclusions. Section 8 completes this work.

2. Prerequisite Theory

In this segment, we list some related theory on fractional-order systems. Let R + = { x | x 0 , x R } .
Definition 1
([44]). Define the Caputo fractional order derivative as follows:
D r κ ( ϱ ) = 1 Γ ( l r ) ϱ 0 ϱ κ ( l ) ( u ) ( ϱ u ) r l + 1 d u ,
where κ ( ϱ ) ( [ ϱ 0 , ) , R ) , Γ ( ϱ ) = 0 ϱ s 1 e ϱ d ϱ ,   ϱ ϱ 0 and l Z + , r ( r 1 , r ) .
The Laplace transform of Caputo-type fractional order derivative is given by
L { D r g ( t ) ; s } = s r G ( s ) j = 0 n 1 s r l 1 g ( j ) ( 0 ) , r ( n 1 , n ) , n Z + ,
where G ( s ) = L { g ( t ) } . If g ( j ) ( 0 ) = 0 , j = 1 , 2 , , n , then L { D r g ( t ) ; s } = s r G ( s ) .
Definition 2
([45]). ( w 1 , w 2 , w 3 ) is an equilibrium point of model (2) provided that
ν w 1 + h ( w 1 ) + h 12 ( w 2 ) + h 13 ( w 3 ) = 0 , ν w 2 + h ( w 2 ) + h 21 ( w 1 ) + h 23 ( w 3 ) = 0 , ν w 3 + h ( w 3 ) + h 31 ( w 1 ) + h 32 ( w 2 ) = 0 .
Lemma 1
([46]). Give the fractional-order system as follows:
d r u ( t ) d t r = g ( t , u ( t ) ) , u ( t 0 ) = u t 0 , t 0 > 0 ,
where t 0 R , r ( 0 , 1 ] , g : [ t 0 , ) × φ R + n , φ R + n , then system (5) has a unique solution defined on [ t 0 , ) if g ( t , u ) satisfies the local Lipschitz condition with respect to u.
Lemma 2
([47]). Assume that u ( t ) C [ t 0 , ) and satisfies
D r u ( t ) ζ 1 u ( t ) + ζ 2 , u ( t 0 ) = u t 0 ,
where 0 < r < 1 , ζ 1 , ζ 2 R , ζ 1 0 , t 0 0 . Then,
u ( t ) u ( t 0 ) ζ 2 ζ 1 E r [ ζ 1 ( t t 0 ) r ] + ζ 2 ζ 1 .
Lemma 3
([48]). Consider the following fractional-order model:
d r u ( t ) d t r = g ( t , u ( t ) ) , u ( 0 ) = u 0 ,
where r ( 0 , 1 ] , g ( t , u ( t ) ) : R + × R m R m . Denote u the equilibrium point of system (8). Then u is locally asymptotically stable provided that every eigenvalue ϖ of g ( t , u ) u | u = u satisfies | a r g ( ϖ ) | > r π 2 .
Lemma 4
([49]). Give the following fractional-order system as follows:
d r 1 P 1 ( t ) d t r 1 = τ 11 P 1 ( t ρ 11 ) + τ 12 P 2 ( t ρ 12 ) + + τ 1 m P m ( t ρ 1 m ) , d r 2 P 2 ( t ) d t r 2 = τ 21 P 1 ( t ρ 21 ) + τ 22 P 2 ( t ρ 22 ) + + τ 2 m P m ( t ρ 2 m ) , d r m P m ( t ) d t r m = τ m 1 P 1 ( t ρ m 1 ) + τ n 2 P 2 ( t ρ m 2 ) + + τ m m P m ( t ρ m m ) ,
where r i ( 0 , 1 ) ( i = 1 , 2 , , m ) . Set
Δ ( σ ) = σ r 1 τ 11 e σ ρ 11 τ 12 e σ ρ 12 τ 1 m e σ ρ 1 m τ 21 e σ ρ 12 σ r 2 τ 22 e σ ρ 22 τ 2 m e σ ρ 2 m τ m 1 e σ ρ m 1 τ m 2 e σ ρ m 2 σ r m τ m m e σ ρ m m . .
Then the equilibrium point of system (9) is asymptotically stable in Lyapunov sense if every root of det ( Δ ( σ ) ) = 0 has negative real parts.

3. Existence and Uniqueness

Theorem 1.
Denote Ψ = { ( w 1 , w 2 , w 3 ) R 3 : max { | w 1 | , | w 2 | , | w 3 | } < W } , where W > 0 stands for a constant. Then for arbitrary W t 0 = ( w 1 t 0 , w 2 t 0 , w 3 t 0 ) Ψ and for all t t 0 , then system (2) with the initial value W t 0 has a unique solution W ( t ) Ψ .
Proof. 
Let W = ( w 1 , w 2 , w 3 ) and W ˜ = ( w ˜ 1 , w ˜ 2 , w ˜ 3 ) and Π ( W ) = ( Π 1 ( W ) , Π 2 ( W ) , Π 3 ( W ) ) , where
Π 1 ( V ) = ν w 1 ( t ) + h ( w 1 ( t ) ) + h 12 ( w 2 ( t ρ ) ) + h 13 ( w 3 ( t ρ ) ) , Π 2 ( V ) = ν w 2 ( t ) + h ( w 2 ( t ) ) + h 21 ( w 1 ( t ρ ) ) + h 23 ( w 3 ( t ρ ) ) , Π 3 ( V ) = ν w 3 ( t ) + h ( w 3 ( t ) ) + h 31 ( w 1 ( t ρ ) ) + h 32 ( w 2 ( t ρ ) ) .
W , W ˜ Ψ , by virtue of ( Q 2 ) and ( Q 3 ) , we have
| | Π ( W ) Π ( W ˜ ) | | = j = 1 3 | Π j ( W ) Π j ( W ˜ ) | = | [ ν w 1 ( t ) + h ( w 1 ( t ) ) + h 12 ( w 2 ( t ρ ) ) + h 13 ( w 3 ( t ρ ) ) ] [ ν w ˜ 1 ( t ) + h ( w ˜ 1 ( t ) ) + h 12 ( w ˜ 2 ( t ρ ) ) + h 13 ( w ˜ 3 ( t ρ ) ) ] | + | [ ν w 2 ( t ) + h ( w 2 ( t ) ) + h 21 ( w 1 ( t ρ ) ) + h 23 ( w 3 ( t ρ ) ) ] [ ν w ˜ 2 ( t ) + h ( w ˜ 2 ( t ) ) + h 21 ( w ˜ 1 ( t ρ ) ) + h 23 ( w ˜ 3 ( t ρ ) ) ] | + | [ ν w 3 ( t ) + h ( w 3 ( t ) ) + h 31 ( w 1 ( t ρ ) ) + h 32 ( w 2 ( t ρ ) ) ] [ ν w ˜ 3 ( t ) + h ( w ˜ 3 ( t ) ) + h 31 ( w ˜ 1 ( t ρ ) ) + h 32 ( w ˜ 2 ( t ρ ) ) ] | ν | w 1 ( t ) w ˜ 1 ( t ) | + H | w 1 ( t ) w ˜ 1 ( t ) | + H 12 | w 2 ( t ρ ) w ˜ 2 ( t ρ ) | + H 13 | w 3 ( t ρ ) w ˜ 3 ( t ρ ) | + ν | w 2 ( t ) w ˜ 2 ( t ) | + H | w 2 ( t ) w ˜ 2 ( t ) | + H 21 | w 1 ( t ρ ) w ˜ 1 ( t ρ ) | + H 23 | w 3 ( t ρ ) w ˜ 3 ( t ρ ) | + ν | w 3 ( t ) w ˜ 3 ( t ) | + H | w 3 ( t ) w ˜ 3 ( t ) | + H 31 | w 1 ( t ρ ) w ˜ 1 ( t ρ ) | + H 32 | w 2 ( t ρ ) w ˜ 2 ( t ρ ) | A 1 | w 1 ( t ) w ˜ 1 ( t ) | + A 2 | w 2 ( t ) w ˜ 2 ( t ) | + A 3 | w 3 ( t ) w ˜ 3 ( t ) | ,
where
A 1 = ν + H + H 21 + H 31 , A 2 = ν + H + H 12 + H 32 , A 3 = ν + H + H 13 + H 23 .
By (12), we obtain
| | Π ( W ) Π ( W ˜ ) | | A | | W W ˜ | ,
where
A = max { A 1 , A 2 , A 3 } .
By virtue of Lemma 1, we easily know that Theorem 1 is true. The proof completes. □
Remark 1.
According to Lemma 1, we know that Π ( W ) satisfies the local Lipschitz condition with respect to W . By virtue of fixed point theorem, we can conclude that Theorem 1 is correct.

4. Boundedness

In this segment, we will prove the boundedness of the solution to system (2). Denote Ψ + = { ( w 1 , w 2 , w 3 ) Ψ : w 1 , w 2 , w 3 R + } .
Theorem 2.
All solutions of system (2) that start with Ψ + are uniformly bounded.
Proof. 
Define
V ( t ) = w 1 ( t ) + w 2 ( t ) + w 3 ( t ) .
It follows from (2) that:
d r V ( t ) d t r = d r w 1 ( t ) d t r + d r w 2 ( t ) d t r + d r w 3 ( t ) d t r = ν w 1 ( t ) + h ( w 1 ( t ) ) + h 12 ( w 2 ( t ρ ) ) + h 13 ( w 3 ( t ρ ) ) ν w 2 ( t ) + h ( w 2 ( t ) ) + h 21 ( w 1 ( t ρ ) ) + h 23 ( w 3 ( t ρ ) ) ν w 3 ( t ) + h ( w 3 ( t ) ) + h 31 ( w 1 ( t ρ ) ) + h 32 ( w 2 ( t ρ ) ) ν w 1 ( t ) + H + H 12 + H 13 ν w 2 ( t ) + H + H 21 + H 23 ν w 3 ( t ) + H + H 31 + H 32 ν V ( t ) + H 0 ,
where
H 0 = 3 H + H 12 + H 13 + H 21 + H 23 + H 31 + H 32 .
According to Lemma 2, we get
V ( t ) V t 0 H 0 ν E r [ ν ( t t 0 ) r ] + H 0 ν H 0 ν , t .
Therefore, all solutions to system (2) which begin with Ψ + are uniformly bounded. The proof finishes. □

5. Bifurcation Study

In this segment, we will analyze the stability and the creation of Hopf bifurcation of system (2). By ( Q 1 ) , we can easily know that system (2) has a unique zero equilibrium point W 0 ( 0 , 0 , 0 ) . The linear system of model (2) at W 0 ( 0 , 0 , 0 ) is:
d r w 1 ( t ) d t r = α w 1 ( t ) + a 12 w 2 ( t ρ ) + a 13 w 3 ( t ρ ) , d r w 2 ( t ) d t r = α w 2 ( t ) + a 21 w 1 ( t ρ ) + a 23 w 3 ( t ρ ) , d r w 3 ( t ) d t r = α w 3 ( t ) + a 31 w 1 ( t ρ ) + a 32 w 2 ( t ρ ) ,
where α = ν h ( 0 ) , a 12 = h 12 ( 0 ) , a 13 = h 13 ( 0 ) , a 21 = h 21 ( 0 ) , a 23 = h 23 ( 0 ) , a 31 = h 31 ( 0 ) , a 32 = h 32 ( 0 ) . The corresponding characteristic equation for (20) is given by:
det s r + α a 12 e s ρ a 13 e s ρ a 21 e s ρ s r + α a 23 e s ρ a 31 e s ρ a 32 e s ρ s r + α = 0 .
It follows from (21) that:
C 1 ( s ) + C 2 ( s ) e 2 s ρ + C 3 ( s ) e 3 s ρ = 0 ,
where
C 1 ( s ) = s 3 r + ϱ 1 s 2 r + ϱ 2 s r + ϱ 3 , C 2 ( s ) = ϱ 4 s r + ϱ 5 , C 3 ( s ) = ϱ 6 ,
where
ϱ 1 = 3 α , ϱ 2 = 3 α 2 , ϱ 3 = α 3 , ϱ 4 = ( a 12 a 21 + a 13 a 31 + a 23 a 32 ) , ϱ 5 = α ( a 12 a 21 + a 13 a 31 + a 23 a 32 ) , ϱ 6 = ( a 12 a 23 a 31 + a 13 a 21 a 32 ) .
We rewrite (22) as:
C 1 ( s ) e s ρ + C 2 ( s ) e s ρ + ϱ 6 e 2 s ρ = 0 .
Let s = i θ = θ cos π 2 + i sin π 2 be the root of Equation (25) and let C l R ( s ) and C l I ( s ) denote the real part and the imaginary part of C l ( s ) ( j = 1 , 2 ) . Then one has
C l R ( s ) = θ 3 r cos 3 r π 2 + ϱ 1 θ 2 r cos r π + ϱ 2 θ r cos r π 2 + ϱ 3 , C l I ( s ) = θ 3 r sin 3 r π 2 + ϱ 1 θ 2 r sin r π + ϱ 2 θ r sin r π 2 , C 2 R ( s ) = ϱ 4 θ r cos r π 2 + ϱ 5 , C 2 I ( s ) = ϱ 4 θ r sin r π 2 .
Let
d 1 = cos 3 r π 2 , d 2 = ϱ 1 cos r π , d 3 = ϱ 2 cos r π 2 , d 4 = ϱ 3 , d 5 = sin 3 r π 2 , d 6 = ϱ 1 sin r π , d 7 = ϱ 2 sin r π 2 , d 8 = ϱ 4 cos r π 2 , d 9 = ϱ 5 , d 10 = ϱ 4 sin r π 2 ,
then (26) becomes
C l R ( s ) = d 1 θ 3 r + d 2 θ 2 r + d 3 θ r + d 4 , C l I ( s ) = d 5 θ 3 r + d 6 θ 2 r + d 7 θ r , C 2 R ( s ) = d 8 θ r + d 9 , C 2 I ( s ) = d 10 θ r .
By virtue of (25) and (28), we have:
[ ( d 1 θ 3 r + d 2 θ 2 r + d 3 θ r + d 4 ) + i ( d 5 θ 3 r + d 6 θ 2 r + d 7 θ r ) ] ( cos θ ρ + i sin θ ρ ) + [ ( d 8 θ r + d 9 ) + i d 10 θ r ] ( cos θ ρ i sin θ ρ ) + ϱ 6 ( cos 2 θ ρ i sin 2 θ ρ ) = 0 .
It follows from (29) that:
D 1 cos θ ρ D 2 sin θ ρ = ϱ 6 cos 2 θ ρ , D 3 cos θ ρ + D 4 cos θ ρ = ϱ 6 sin 2 θ ρ ,
where
D 1 = d 1 θ 3 r + d 2 θ 2 r + ( d 3 + d 8 ) θ r + d 4 + d 9 , D 2 = d 5 θ 3 r + d 6 θ 2 r + ( d 7 + d 10 ) θ r , D 3 = d 5 θ 3 r + d 6 θ 2 r + ( d 7 d 10 ) θ r , D 4 = d 1 θ 3 r + d 2 θ 2 r + ( d 3 d 8 ) θ r + d 4 d 9 .
In view of (30), we get:
[ D 1 cos θ ρ D 2 sin θ ρ ] 2 + [ D 3 cos θ ρ + D 4 cos θ ρ ] 2 = ϱ 6 2 ,
which leads to:
( D 1 2 + D 3 2 ) cos 2 θ ρ + ( D 2 2 + D 4 2 ) sin 2 θ ρ + 2 ( D 3 D 4 D 1 D 2 ) cos θ ρ sin θ ρ = ϱ 6 2 .
By (33), we have:
( D 1 2 + D 3 2 ) cos 2 θ ρ + ( D 2 2 + D 4 2 ) sin 2 θ ρ ϱ 6 2 = 2 ( D 3 D 4 D 1 D 2 ) cos θ ρ sin θ ρ ,
which leads to
[ ( D 1 2 + D 3 2 ) cos 2 θ ρ + ( D 2 2 + D 4 2 ) sin 2 θ ρ ϱ 6 2 ] 2 = 4 ( D 3 D 4 D 1 D 2 ) 2 cos 2 θ ρ sin 2 θ ρ .
It follows from (35) that:
[ ( D 1 2 + D 3 2 ) cos 2 θ ρ + ( D 2 2 + D 4 2 ) sin 2 θ ρ ϱ 6 2 ] 2 = 4 ( D 3 D 4 D 1 D 2 ) 2 cos 2 θ ρ sin 2 θ ρ .
Namely,
[ ( D 1 2 + D 3 2 ) cos 2 θ ρ + ( D 2 2 + D 4 2 ) sin 2 θ ρ ϱ 6 2 ] 2 = 4 ( D 3 D 4 D 1 D 2 ) 2 cos 2 θ ρ sin 2 θ ρ .
Then
B 1 cos 4 θ ρ + B 2 cos 2 θ ρ + B 3 = 0 ,
where
B 1 = ( D 1 2 + D 3 2 D 2 2 D 4 2 ) 2 + 4 ( D 3 D 4 D 1 D 2 ) 2 , B 2 = 2 ( D 1 2 + D 3 2 ) cos 2 θ ρ + ( D 2 2 + D 4 2 ) ( D 2 2 + D 4 2 ϱ 6 2 ) , B 3 = ( D 2 2 + D 4 2 ϱ 6 2 ) 2 .
Assume that:
( Q 4 ) B 2 < 0 , B 2 2 4 B 1 B 3
holds, then by (38), we get
cos 2 θ ρ = B 2 + B 2 2 4 B 1 B 3 2 B 1 ,
or
cos 2 θ ρ = B 2 B 2 2 4 B 1 B 3 2 B 1 .
According to (40) and (41), we get:
cos θ ρ = ± B 2 + B 2 2 4 B 1 B 3 2 B 1 1 2 ,
or
cos θ ρ = ± B 2 B 2 2 4 B 1 B 3 2 B 1 1 2 .
Suppose that (38) has twelve real roots θ j ( j = 1 , 2 , , 6 ) . By (42) and (43), we have:
ρ j 1 l = 1 ρ j arccos B 2 + B 2 2 4 B 1 B 3 2 B 1 1 2 + 2 l π ,
ρ j 2 l = 1 ρ j arccos B 2 + B 2 2 4 B 1 B 3 2 B 1 1 2 + 2 l π ,
ρ j 3 l = 1 ρ j arccos B 2 B 2 2 4 B 1 B 3 2 B 1 1 2 + 2 l π ,
ρ j 4 l = 1 ρ j arccos B 2 B 2 2 4 B 1 B 3 2 B 1 1 2 + 2 l π ,
where l = 0 , 1 , 2 , , j = 1 , 2 , , 12 . Let
ρ 0 = min k = 1 , 2 , 3 , 4 { ρ j k 0 } , θ 0 = θ | ρ = ρ 0 .
In the sequel, the hypothesis is given:
  • ( Q 5 ) S 11 S 21 + S 12 S 22 > 0 , where
    S 11 = 3 r θ 0 3 r 1 cos ( 3 r 1 ) π 2 + 2 r ϱ 1 θ 0 2 r 1 cos ( 2 r 1 ) π 2 + r ϱ 2 θ 0 r 1 cos ( r 1 ) π 2 + r ϱ 4 θ 0 r 1 cos ( r 1 ) π 2 cos 2 θ 0 ρ 0 + r ϱ 4 θ 0 r 1 sin ( r 1 ) π 2 sin 2 θ 0 ρ 0 , S 12 = 3 r θ 0 3 r 1 sin ( 3 r 1 ) π 2 + 2 r ϱ 1 θ 0 2 r 1 sin ( 2 r 1 ) π 2 + r ϱ 2 θ 0 r 1 sin ( r 1 ) π 2 r ϱ 4 θ 0 r 1 cos ( r 1 ) π 2 sin 2 θ 0 ρ 0 + r ϱ 4 θ 0 r 1 sin ( r 1 ) π 2 cos 2 θ 0 ρ 0 , S 21 = ρ 0 θ 0 r + 1 sin r π 2 cos 2 θ 0 ρ 0 + θ 0 ϱ 4 θ 0 r cos r π 2 + ϱ 5 sin 2 θ 0 ρ 0 + 3 ρ 0 θ 0 sin 3 θ 0 ρ 0 , S 22 = ρ 0 θ 0 r + 1 sin r π 2 sin 2 θ 0 ρ 0 + θ 0 ϱ 4 θ 0 r cos r π 2 + ϱ 5 cos 2 θ 0 ρ 0 + 3 ρ 0 θ 0 cos 3 θ 0 ρ 0 .
Lemma 5.
Let s ( ρ ) = ι 1 ( ρ ) + i ι 2 ( ρ ) be the root of Equation (22) at ρ = ρ 0 such that ι 1 ( ρ 0 ) = 0 , ι 2 ( ρ 0 ) = ϖ 0 , then R e d s d ρ | ρ = ρ 0 , θ = θ 0 > 0 .
Proof. 
It follows from Equation (22) that:
3 r s 3 r 1 + 2 r ϱ 1 s 2 r 1 + r ϱ 2 s r 1 d s d ρ + r ϱ 4 s r 1 e 2 s ρ d s d ρ 2 e 2 s ρ d s d ρ ρ + s ϱ 4 s r + ϱ 5 3 ϱ 6 e 3 s ρ d s d ρ ρ + s = 0 .
Then
d s d ρ 1 = S 1 ( s ) S 2 ( s ) ρ s ,
where
S 1 ( s ) = 3 r s 3 r 1 + 2 r ϱ 1 s 2 r 1 + r ϱ 2 s r 1 + r ϱ 4 s r 1 e 2 s ρ , S 2 ( s ) = 2 s e 2 s ρ ( ϱ 4 s r + ϱ 5 ) + 3 s ϱ 6 e 3 s ρ .
Then
Re d s d ρ | ρ = ρ 0 , θ = θ 0 = Re S 1 ( s ) S 2 ( s ) | ρ = ρ 0 , θ = θ 0 = S 11 S 21 + S 12 S 22 S 21 2 + S 22 2 .
By ( Q 5 ) , one gets:
Re d s d ρ 1 | ρ = ρ 0 , θ = θ 0 > 0 ,
which completes the proof. □
Now we give the following hypothesis:
  • ( Q 6 ) The following inequalities hold:
    G 1 = ϱ 1 > 0 , G 2 = det ϱ 1 1 ϱ 3 + ϱ 5 + ϱ 6 ϱ 2 + ϱ 4 > 0 , G 3 = ( ϱ 3 + ϱ 5 + ϱ 6 ) G 2 > 0 .
Lemma 6.
If ( Q 6 ) is satisfied, then the zero equilibrium point of system (2) with ρ = 0 is locally asymptotically stable.
Proof. 
When ρ = 0 , then Equation (22) takes the form:
λ 3 + ϱ 1 λ 2 + ( ϱ 2 + ϱ 4 ) λ + ϱ 3 + ϱ 5 + ϱ 6 = 0 .
In view of ( Q 6 ) , one knows that all roots λ j of (54) satisfy | arg ( λ j ) | > r π 2 ( j = 1 , 2 , 3 ) . Thus Lemma 6 is true. The proof ends. □
Based on the analysis above, we can easily derive the following conclusion.
Theorem 3.
Assume that ( Q 1 ) ( Q 6 ) are fulfilled, then the zero equilibrium point W 0 ( 0 , 0 , 0 ) of system (2) is locally asymptotically stable if 0 ρ < ρ 0 and a Hopf bifurcation arises near W 0 ( 0 , 0 , 0 ) if ρ = ρ 0 .

6. Bifurcation Control

In this segment, we will explore the control issue of the Hopf bifurcation of system (2) via mixed controller including state feedback and parameter perturbation. Following the idea of Yuan et al. [50], we establish the following controlled system:
d r w 1 ( t ) d t r = ε 1 [ ν w 1 ( t ) + h ( w 1 ( t ) ) + h 12 ( w 2 ( t ρ ) ) + h 13 ( w 3 ( t ρ ) ) ] + ε 2 w 1 ( t ) , d r w 2 ( t ) d t r = ε 1 [ ν w 2 ( t ) + h ( w 2 ( t ) ) + h 21 ( w 1 ( t ρ ) ) + h 23 ( w 3 ( t ρ ) ) ] + ε 2 w 2 ( t ) , d r w 3 ( t ) d t r = ε 1 [ ν w 3 ( t ) + h ( w 3 ( t ) ) + h 31 ( w 1 ( t ρ ) ) + h 32 ( w 2 ( t ρ ) ) ] + ε 2 w 3 ( t ) ,
where ε 1 and ε 2 are feedback gain parameters and real constants. It is easy to see that system (57) has a unique equilibrium point W 0 ( 0 , 0 , 0 ) . The linear system of model (57) at W 0 ( 0 , 0 , 0 ) is:
d r w 1 ( t ) d t r = β w 1 ( t ) + b 12 w 2 ( t ρ ) + b 13 w 3 ( t ρ ) , d r w 2 ( t ) d t r = β w 2 ( t ) + b 21 w 1 ( t ρ ) + b 23 w 3 ( t ρ ) , d r w 3 ( t ) d t r = β w 3 ( t ) + b 31 w 1 ( t ρ ) + b 32 w 2 ( t ρ ) ,
where β = ε 1 ( ν h ( 0 ) ) , b 12 = ε 1 h 12 ( 0 ) , b 13 = ε 1 h 13 ( 0 ) , b 21 = ε 1 h 21 ( 0 ) , b 23 = ε 1 h 23 ( 0 ) , b 31 = ε 1 h 31 ( 0 ) , b 32 = ε 1 h 32 ( 0 ) . The corresponding characteristic equation for (58) is:
det s r + β b 12 e s ρ b 13 e s ρ b 21 e s ρ s r + β b 23 e s ρ b 31 e s ρ b 32 e s ρ s r + β = 0 .
It follows from (58) that
H 1 ( s ) + H 2 ( s ) e 2 s ρ + H 3 ( s ) e 3 s ρ = 0 ,
where
H 1 ( s ) = s 3 r + δ 1 s 2 r + δ 2 s r + δ 3 , H 2 ( s ) = δ 4 s r + δ 5 , H 3 ( s ) = δ 6 ,
where
δ 1 = 3 β , δ 2 = 3 β 2 , δ 3 = β 3 , δ 4 = ( b 12 b 21 + b 13 b 31 + b 23 b 32 ) , δ 5 = α ( b 12 b 21 + b 13 b 31 + b 23 b 32 ) , δ 6 = ( b 12 b 23 b 31 + b 13 b 21 b 32 ) .
We rewrite (60) as:
H 1 ( s ) e s ρ + H 2 ( s ) e s ρ + δ 6 e 2 s ρ = 0 .
Let s = i ϑ = ϑ cos π 2 + i sin π 2 be the root of Equation (63) and let H l R ( s ) and H l I ( s ) denote the real part and the imaginary part of H l ( s ) ( l = 1 , 2 ) . Then one has
H l R ( s ) = ϑ 3 r cos 3 r π 2 + δ 1 ϑ 2 r cos r π + δ 2 ϑ r cos r π 2 + δ 3 , H l I ( s ) = ϑ 3 r sin 3 r π 2 + δ 1 ϑ 2 r sin r π + δ 2 ϑ r sin r π 2 , H 2 R ( s ) = δ 4 ϑ r cos r π 2 + δ 5 , H 2 I ( s ) = δ 4 ϑ r sin r π 2 .
Let
e 1 = cos 3 r π 2 , e 2 = δ 1 cos r π , e 3 = δ 2 cos r π 2 , e 4 = δ 3 , e 5 = sin 3 r π 2 , e 6 = δ 1 sin r π , e 7 = δ 2 sin r π 2 , e 8 = δ 4 cos r π 2 , e 9 = δ 5 , e 10 = δ 4 sin r π 2 ,
then (64) becomes
H l R ( s ) = e 1 ϑ 3 r + e 2 ϑ 2 r + e 3 ϑ r + e 4 , H l I ( s ) = e 5 ϑ 3 r + e 6 ϑ 2 r + e 7 ϑ r , H 2 R ( s ) = e 8 ϑ r + e 9 , H 2 I ( s ) = e 10 ϑ r .
By virtue of (63) and (66), we have:
[ ( e 1 ϑ 3 r + e 2 ϑ 2 r + e 3 ϑ r + e 4 ) + i ( e 5 ϑ 3 r + e 6 ϑ 2 r + e 7 ϑ r ) ] ( cos ϑ ρ + i sin ϑ ρ ) + [ ( e 8 ϑ r + e 9 ) + i e 10 ϑ r ] ( cos ϑ ρ i sin ϑ ρ ) + δ 6 ( cos 2 ϑ ρ i sin 2 ϑ ρ ) = 0 .
It follows from (67) that
I 1 cos ϑ ρ I 2 sin ϑ ρ = δ 6 cos 2 ϑ ρ , I 3 cos ϑ ρ + I 4 cos ϑ ρ = δ 6 sin 2 ϑ ρ ,
where
I 1 = e 1 ϑ 3 r + e 2 ϑ 2 r + ( e 3 + e 8 ) ϑ r + e 4 + e 9 , I 2 = e 5 ϑ 3 r + e 6 ϑ 2 r + ( e 7 + e 10 ) ϑ r , I 3 = e 5 ϑ 3 r + e 6 ϑ 2 r + ( e 7 e 10 ) ϑ r , I 4 = e 1 ϑ 3 r + e 2 ϑ 2 r + ( e 3 e 8 ) ϑ r + e 4 e 9 .
In view of (69), we get:
[ I 1 cos ϑ ρ I 2 sin ϑ ρ ] 2 + [ I 3 cos ϑ ρ + I 4 cos ϑ ρ ] 2 = δ 6 2 ,
which leads to:
( I 1 2 + I 3 2 ) cos 2 ϑ ρ + ( I 2 2 + I 4 2 ) sin 2 ϑ ρ + 2 ( I 3 I 4 I 1 I 2 ) cos ϑ ρ sin ϑ ρ = δ 6 2 .
By (71), we have:
( I 1 2 + I 3 2 ) cos 2 ϑ ρ + ( I 2 2 + I 4 2 ) sin 2 ϑ ρ δ 6 2 = 2 ( I 3 I 4 I 1 I 2 ) cos ϑ ρ sin ϑ ρ ,
which leads to:
[ ( I 1 2 + I 3 2 ) cos 2 ϑ ρ + ( I 2 2 + I 4 2 ) sin 2 ϑ ρ δ 6 2 ] 2 = 4 ( I 3 I 4 I 1 I 2 ) 2 cos 2 ϑ ρ sin 2 ϑ ρ .
It follows from (73) that:
[ ( I 1 2 + I 3 2 ) cos 2 ϑ ρ + ( I 2 2 + I 4 2 ) ( 1 cos 2 ϑ ρ ) δ 6 2 ] 2 = 4 ( I 3 I 4 I 1 I 2 ) 2 cos 2 ϑ ρ ( 1 cos 2 ϑ ρ ) .
Namely,
[ ( I 1 2 + I 3 2 I 2 2 I 4 2 ) cos 2 ϑ ρ + I 2 2 + I 4 2 δ 6 2 ] 2 = 4 ( I 3 I 4 I 1 I 2 ) 2 cos 2 ϑ ρ ( 1 cos 2 ϑ ρ ) .
Then,
J 1 cos 4 ϑ ρ + J 2 cos 2 ϑ ρ + J 3 = 0 ,
where
J 1 = ( I 1 2 + I 3 2 I 2 2 I 4 2 ) 2 + 4 ( I 3 I 4 I 1 I 2 ) 2 , J 2 = 2 ( I 1 2 + I 3 2 ) cos 2 ϑ ρ + ( I 2 2 + I 4 2 ) ( I 2 2 + I 4 2 δ 6 2 ) , J 3 = ( I 2 2 + I 4 2 δ 6 2 ) 2 .
Assume that
( Q 7 ) J 2 < 0 , J 2 2 4 J 1 J 3
holds, then by (76), we get:
cos 2 ϑ ρ = J 2 + J 2 2 4 J 1 J 3 2 J 1 ,
or
cos 2 ϑ ρ = J 2 J 2 2 4 J 1 J 3 2 J 1 .
According to (78) and (79), we can get:
cos ϑ ρ = ± J 2 + J 2 2 4 J 1 J 3 2 J 1 1 2 ,
or
cos ϑ ρ = ± J 2 J 2 2 4 J 1 J 3 2 J 1 1 2 .
Suppose that (76) has twelve real roots ϑ j ( j = 1 , 2 , , 6 ) . By (80) and (81), we have:
ρ i 1 l = 1 ρ i arccos J 2 + J 2 2 4 J 1 J 3 2 J 1 1 2 + 2 l π ,
ρ i 2 l = 1 ρ i arccos J 2 + J 2 2 4 J 1 J 3 2 J 1 1 2 + 2 l π ,
ρ i 3 l = 1 ρ i arccos J 2 J 2 2 4 J 1 J 3 2 J 1 1 2 + 2 l π ,
ρ i 4 l = 1 ρ i arccos J 2 J 2 2 4 J 1 J 3 2 J 1 1 2 + 2 l π ,
where l = 0 , 1 , 2 , , i = 1 , 2 , , 6 . Let
ρ = min l = 1 , 2 , 3 , 4 { ρ i l 0 } , ϑ 0 = ϑ | ρ = ρ .
In the sequel, the hypothesis is given:
  • ( Q 8 ) N 11 N 21 + N 12 N 22 > 0 , where
    N 11 = 3 r ϑ 0 3 r 1 cos ( 3 r 1 ) π 2 + 2 r δ 1 ϑ 0 2 r 1 cos ( 2 r 1 ) π 2 + r δ 2 θ 0 r 1 cos ( r 1 ) π 2 + r δ 4 ϑ 0 r 1 cos ( r 1 ) π 2 cos 2 ϑ 0 ρ + r δ 4 ϑ 0 r 1 sin ( r 1 ) π 2 sin 2 ϑ 0 ρ , N 12 = 3 r ϑ 0 3 r 1 sin ( 3 r 1 ) π 2 + 2 r δ 1 ϑ 0 2 r 1 sin ( 2 r 1 ) π 2 + r δ 2 ϑ 0 r 1 sin ( r 1 ) π 2 r δ 4 ϑ 0 r 1 cos ( r 1 ) π 2 sin 2 ϑ 0 ρ + r δ 4 ϑ 0 r 1 sin ( r 1 ) π 2 cos 2 ϑ 0 ρ , N 21 = ρ ϑ 0 r + 1 sin r π 2 cos 2 ϑ 0 ρ + ϑ 0 δ 4 ϑ 0 r cos r π 2 + δ 5 sin 2 δ 0 ρ + 3 ρ ϑ 0 sin 3 ϑ 0 ρ , N 22 = ρ ϑ 0 r + 1 sin r π 2 sin 2 ϑ 0 ρ + ϑ 0 δ 4 ϑ 0 r cos r π 2 + δ 5 cos 2 ϑ 0 ρ + 3 ρ ϑ 0 cos 3 ϑ 0 ρ .
Lemma 7.
Let s ( ρ ) = φ 1 ( ρ ) + i φ 2 ( ρ ) be the root of Equation (60) at ρ = ρ such that φ 1 ( ρ ) = 0 , φ 2 ( ρ ) = ϑ 0 , then R e d s d ρ | ρ = ρ , ϑ = ϑ 0 > 0 .
Proof. 
It follows from Equation (60) that
3 r s 3 r 1 + 2 r δ 1 s 2 r 1 + r δ 2 s r 1 d s d ρ + r δ 4 s r 1 e 2 s ρ d s d ρ 2 e 2 s ρ d s d ρ ρ + s δ 4 s r + δ 5 3 δ 6 e 3 s ρ d s d ρ ρ + s = 0 .
Then
d s d ρ 1 = N 1 ( s ) N 2 ( s ) ρ s ,
where
N 1 ( s ) = 3 r s 3 r 1 + 2 r δ 1 s 2 r 1 + r δ 2 s r 1 + r δ 4 s r 1 e 2 s ρ , N 2 ( s ) = 2 s e 2 s ρ ( δ 4 s r + δ 5 ) + 3 s δ 6 e 3 s ρ .
Then,
Re d s d ρ | ρ = ρ , ϑ = ϑ 0 = Re N 1 ( s ) N 2 ( s ) | ρ = ρ , ϑ = ϑ 0 = N 11 N 21 + N 12 N 22 N 21 2 + N 22 2 .
By ( Q 8 ) , one gets
Re d s d ρ 1 | ρ = ρ , ϑ = ϑ 0 > 0 ,
which completes the proof. □
Now we give the following hypothesis:
  • ( Q 9 ) The following inequalities hold:
    F 1 = δ 1 > 0 , F 2 = det δ 1 1 δ 3 + δ 5 + δ 6 δ 2 + δ 4 > 0 , F 3 = ( δ 3 + δ 5 + δ 6 ) F 2 > 0 .
Lemma 8.
If ( Q 9 ) is satisfied, then the zero equilibrium point of system (57) with ρ = 0 is locally asymptotically stable.
Proof. 
When ρ = 0 , then Equation (60) takes the form:
λ 3 + δ 1 λ 2 + ( δ 2 + δ 4 ) λ + δ 3 + δ 5 + δ 6 = 0 .
In view of ( Q 9 ) , one knows that all roots λ j of (94) satisfy | arg ( λ j ) | > r π 2 ( j = 1 , 2 , 3 ) . Thus Lemma 8 is true. The proof ends. □
Based on the discussion above, we can easily derive the following conclusion.
Theorem 4.
Assume that ( Q 1 ) , ( Q 2 ) , ( Q 3 ) , ( Q 7 ) , ( Q 8 ) , ( Q 9 ) are fulfilled, then the zero equilibrium point W 0 ( 0 , 0 , 0 ) of system (57) is locally asymptotically stable if 0 ρ < ρ and a Hopf bifurcation arises near W 0 ( 0 , 0 , 0 ) if ρ = ρ .
Remark 2.
In 2006, Yan [41] investigated the Hopf bifurcation of integer-order delayed tri-neuron neural networks. In this present research, on the basis of Yan [41], we build a class of new fractional-order delayed tri-neuron neural networks. The existence and uniqueness, boundedness of the solution, stability and the onset of Hopf bifurcation, Hopf bifurcation control of fractional-order delayed tri-neuron neural networks is discussed in detail. The research approach of Yan [41] can be be applied to the fractional-order delayed tri-neuron neural networks to establish the key results of this article. Based on this standpoint, we think that our work replenishes the work of Yan [41].

7. Numerical Simulations

Example 1.
Give the fractional order system as follows:
d 0.95 w 1 ( t ) d t r = 2 w 1 ( t ) + tanh ( w 1 ( t ) ) + tanh ( w 2 ( t ρ ) ) + tanh ( w 3 ( t ρ ) ) , d 0.95 w 2 ( t ) d t r = 2 w 2 ( t ) + tanh ( w 2 ( t ) ) + tanh ( w 1 ( t ρ ) ) + tanh ( w 3 ( t ρ ) ) , d 0.95 w 3 ( t ) d t r = 2 w 3 ( t ) + tanh ( w 3 ( t ) ) + tanh ( w 1 ( t ρ ) ) + tanh ( w 2 ( t ρ ) ) .
Obviously, system (95) owns a unique equilibrium point W 0 ( 0 , 0 , 0 ) . By virtue of Matlab software, we derive θ 0 = 0.9223 and ρ 0 0.6 . Furthermore, we can verify that all the hypotheses ( Q 1 ) ( Q 6 ) of Theorem 3 are fulfilled. So the zero equilibrium point W 0 ( 0 , 0 , 0 ) of system (95) is locally asymptotically stable if ρ [ 0 , 0.6 ) and system (95) will generates Hopf bifurcation phenomenon near W 0 ( 0 , 0 , 0 ) when ρ crosses ρ 0 0.6 . Now we select ρ = 0.5 is less than ρ 0 0.6 , then we get the simulation results that are displayed in Figure 1, which shows that the three variables w 1 ( t ) , w 2 ( t ) , w 3 ( t ) gradually tend to zero when time t increases. Then we select ρ = 0.7 is bigger than ρ 0 0.6 , then we get the simulation results that are displayed in Figure 2, which shows that the three variables w 1 ( t ) , w 2 ( t ) , w 3 ( t ) keep a periodic oscillation (namely, Hopf bifurcation) near W 0 ( 0 , 0 , 0 ) when time t increases. The bifurcation figures are given in Figure 3, Figure 4 and Figure 5. Figure 3 shows the relation between ρ and w 1 , Figure 4 shows the relation between ρ and w 2 and Figure 5 shows the relation between ρ and w 3 . From Figure 3, Figure 4 and Figure 5, we can see that the bifurcation value ρ 0 0.6 . In addition, the relation of time delay ρ 0 and the amplitude θ 0 is presented in Table 1.
Example 2.
Give the fractional order system as follows:
d 0.95 w 1 ( t ) d t r = 0.6 [ 2 w 1 ( t ) + tanh ( w 1 ( t ) ) + h 12 ( w 2 ( t ρ ) ) + tanh ( w 3 ( t ρ ) ) ] + 0.4 w 1 ( t ) , d 0.95 w 2 ( t ) d t r = 0.6 [ 2 w 2 ( t ) + tanh ( w 2 ( t ) ) + h 21 ( w 1 ( t ρ ) ) + tanh ( w 3 ( t ρ ) ) ] + 0.4 w 2 ( t ) , d 0.95 w 3 ( t ) d t r = 0.6 [ 2 w 3 ( t ) + tanh ( w 3 ( t ) ) + h 31 ( w 1 ( t ρ ) ) + tanh ( w 2 ( t ρ ) ) ] + 0.4 w 3 ( t ) .
Obviously, system (96) owns a unique equilibrium point W 0 ( 0 , 0 , 0 ) . By virtue of Matlab software, we derive ϑ 0 = 1.9070 and ρ 0.3 . Furthermore, we can verify that all the hypotheses ( Q 1 ) , ( Q 2 ) , ( Q 3 ) , ( Q 7 ) , ( Q 8 ) , ( Q 9 ) of Theorem 4 are fulfilled. So the zero equilibrium point W 0 ( 0 , 0 , 0 ) of system (96) is locally asymptotically stable if ρ [ 0 , 0.3 ) and system (96) will generate the Hopf bifurcation phenomenon near W 0 ( 0 , 0 , 0 ) when ρ crosses ρ 0.3 . Now, we select ρ = 0.24 as less than ρ 0.6 , then we get the simulation results that are displayed in Figure 6, which show that the three variables w 1 ( t ) , w 2 ( t ) , w 3 ( t ) gradually tend to zero when time t increases. Then we select ρ = 0.35 as bigger than ρ 0.3 , then we get the simulation results that are displayed in Figure 7, which show that the three variables w 1 ( t ) , w 2 ( t ) , w 3 ( t ) keep a periodic oscillation (namely, Hopf bifurcation) near W 0 ( 0 , 0 , 0 ) when time t increases. The bifurcation figures are given in Figure 8, Figure 9 and Figure 10. Figure 8 shows the relation between ρ and w 1 , Figure 9 shows the relation between ρ and w 2 and Figure 10 shows the relation between ρ and w 3 . From Figure 8, Figure 9 and Figure 10, we can see that the bifurcation value ρ 0.3 . In addition, the relation of time delay ρ 0 and the amplitude ρ are presented in Table 2.
Remark 3.
In system (95), the bifurcation value ρ 0 0.6 , in system (96) (namely, the controlled system of system (95)), the bifurcation value ρ 0 0.3 . we can clearly know that the stability region of system (95) is narrowed and the time of generation of Hopf bifurcation in advance.

8. Conclusions

The investigation of neural networks has become a very important topic in mathematical and engineering areas. Revealing the intrinsic dynamic characteristics of neural networks has aroused the interest of various scholars. In this article, based on the earlier studies, we built a new fractional-order tri-neuron neural network incorporating delays. The existence and uniqueness, boundedness of the solution of the established fractional-order tri-neuron neural networks are investigated. The delay-independent condition to ensure the stability and the onset of the Hopf bifurcation of the involved fractional-order tri-neuron neural networks is derived. The bifurcation value is determined. By virtue of a mixed controller, which includes state feedback and parameter perturbation, the stability region and the time of onset of the Hopf bifurcation of the considered fractional-order tri-neuron neural networks are adjusted. The research shows that time delay is a crucial parameter which affects the stability region and the time of onset of the Hopf bifurcation of the studied fractional-order tri-neuron neural networks and its controlled neural networks. The derived results have some theoretical value in designing and controlling neural networks. Meanwhile, the exploration idea of this article can also be utilized to study a lot of fractional-order delayed neural networks and some related delayed dynamical systems. In addition, in many neural networks, there exist leakage delays. In this paper, we do not involve the Hopf bifurcation control issue of neural networks with leakage delays. We will focus on this topic in near future.

Author Contributions

Methodology, P.L. and C.X.; Software, C.X.; Formal analysis, P.L. and C.X.; Data curation, Y.S.; writing—original draft preparation; P.L. and C.X.; Writing—review and editing, P.L., J.Y. and C.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (No. 62062018) and and Project of High-level Innovative Talents of Guizhou Province ([2016]5651).

Data Availability Statement

No data were used to support this study.

Acknowledgments

This work is supported by National Natural Science Foundation of China (No. 62062018) and and Project of High-level Innovative Talents of Guizhou Province ([2016]5651). The authors would like to thank the referees and the editor for helpful suggestions incorporated into this paper.

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Hopfield, J. Neurons with graded response have collective computational properties like those of two-state neurons. Proc. Natl. Acad. Sci. USA 1984, 81, 3088–3092. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Xu, C.J.; Liao, M.X.; Li, P.L.; Liu, Z.X.; Yuan, S. New results on pseudo almost periodic solutions of quaternion-valued fuzzy cellular neural networks with delays. Fuzzy Set. Syst. 2021, 411, 25–47. [Google Scholar] [CrossRef]
  3. Thoiyab, N.M.; Muruganantham, P.; Zhu, Q.X.; Gunasekaran, N. Novel results on global stability analysis for multiple time-delayed BAM neural networks under parameter uncertainties. Chaos Solitons Fractals 2021, 152, 111441. [Google Scholar] [CrossRef]
  4. Cao, Y.; Ramajayam, S.; Sriraman, R.; Samidurai, R. Leakage delay on stabilization of finite-time complex-valued BAM neural network: Decomposition approach. Neurocomputing 2021, 463, 505–513. [Google Scholar] [CrossRef]
  5. Xu, C.J.; Tang, X.H.; Liao, M.X. Frequency domain analysis for bifurcation in a simplified tri-neuron BAM network model with two delays. Neural Netw. 2010, 23, 872–880. [Google Scholar] [CrossRef]
  6. Wei, J.J.; Velarde, M.G. Bifurcation analysis and existence of periodic solutions in a simple neural network with delays. Chaos 2004, 14, 940–953. [Google Scholar] [CrossRef]
  7. Wang, S.Z.; Zhang, Z.Y.; Lin, C.; Chen, J. Fixed-time synchronization for complex-valued BAM neural networks with time-varying delays via pinning control and adaptive pinning control. Chaos Solitons Fractals 2021, 153, 111583. [Google Scholar] [CrossRef]
  8. Gan, Y.T.; Liu, C.; Peng, H.; Liu, F.; Rao, H.X. Anti-synchronization for periodic BAM neural networks with Markov scheduling protocol. Neurocomputing 2020, 417, 585–592. [Google Scholar] [CrossRef]
  9. Syed Ali, M.; Narayanan, G.; Shekher, V.; Alsaedi, A.; Ahmad, B. Global Mittag–Leffler stability analysis of impulsive fractional-order complex-valued BAM neural networks with time varying delays. Commun. Nonlinear Sci. Numer. Simul. 2020, 83, 105088. [Google Scholar] [CrossRef]
  10. Chen, S.L.; Li, H.L.; Kao, Y.G.; Zhang, L.; Hu, C. Finite-time stabilization of fractional-order fuzzy quaternion-valued BAM neural networks via direct quaternion approach. J. Franklin Inst. 2021, 358, 7650–7673. [Google Scholar] [CrossRef]
  11. Cong, E.Y.; Han, X.; Zhang, X. Global exponential stability analysis of discrete-time BAM neural networks with delays: A mathematical induction approach. Neurocomputing 2020, 379, 227–235. [Google Scholar] [CrossRef]
  12. Duan, L.Y.; Li, J.M. Fixed-time synchronization of fuzzy neutral-type BAM memristive inertial neural networks with proportional delays. Inf. Sci. 2021, 576, 522–541. [Google Scholar] [CrossRef]
  13. Yan, M.Y.; Jian, J.G.; Zheng, S. Passivity analysis for uncertain BAM inertial neural networks with time-varying delays. Neurocomputing 2021, 435, 114–125. [Google Scholar]
  14. Ayachi, M. Measure-pseudo almost periodic dynamical behaviors for BAM neural networks with D operator and hybrid time-varying delays. Neurocomputing 2021. [Google Scholar] [CrossRef]
  15. Kumar, R.; Das, S. Exponential stability of inertial BAM neural network with time-varying impulses and mixed time-varying delays via matrix measure approach. Commun. Nonlinear Sci. Numer. Simul. 2020, 81, 105016. [Google Scholar] [CrossRef]
  16. Syed Ali, M.; Yogambigai, J.; Saravanan, S.; Elakkia, S. Stochastic stability of neutral-type Markovian-jumping BAM neural networks with time varying delays. J. Comput. Appl. Math. 2019, 349, 142–156. [Google Scholar]
  17. Bruce, W.J. Colloquium: Fractional calculus view of complexity: A tutorial. Rev. Modern Phys. 2014, 86, 1169–1186. [Google Scholar]
  18. Shi, J.P.; He, K.; Fang, H. Chaos, Hopf bifurcation and control of a fractional-order delay financial system. Math. Comput. Simul. 2022, 194, 348–364. [Google Scholar] [CrossRef]
  19. Zhou, B.; Song, Q.K. Stability and Hopf bifurcation analysis of a tri-neuron BAM neural network with distributed delay. Neurocomputing 2012, 82, 69–83. [Google Scholar] [CrossRef]
  20. Pratap, A.; Raja, R.; Cao, J.; Rihan, F.A.; Seadawy, A.R. Quasi-pinning synchronization and stabilization of fractional order BAM neural networks with delays and discontinuous neuron activations. Chaos Solitons Fractals 2020, 131, 109491. [Google Scholar] [CrossRef]
  21. Petráš, I. Control of fractional order Chua’s system. J. Elect. Eng. 2002, 53, 219–222. [Google Scholar]
  22. Zhang, F.H.; Huang, T.W.; Wu, Q.J.; Zeng, Z.G. Multistability of delayed fractional-order competitive neural networks. Neural Netw. 2021, 140, 325–335. [Google Scholar] [CrossRef] [PubMed]
  23. Li, H.L.; Hu, C.; Zhang, L.; Jiang, H.J.; Cao, J.D. Complete and finite-time synchronization of fractional-order fuzzy neural networks via nonlinear feedback control. Fuzzy Set. Syst. 2021. [Google Scholar] [CrossRef]
  24. Padmaja, N.; Balasubramaniam, P. New delay and order-dependent passivity criteria for impulsive fractional-order neural networks with switching parameters and proportional delays. Neurocomputing 2021, 454, 113–123. [Google Scholar] [CrossRef]
  25. Ke, L. Mittag–Leffler stability and asymptotic ω-periodicity of fractional-order inertial neural networks with time-delays. Neurocomputing 2021, 465, 53–62. [Google Scholar] [CrossRef]
  26. Du, F.F.; Lu, J.G. Finite-time stability of fractional-order fuzzy cellular neural networks with time delays. Fuzzy Set. Syst. 2021. [Google Scholar] [CrossRef]
  27. Xu, C.J.; Zhang, W.; Liu, Z.X.; Yao, L.Y. Delay-induced periodic oscillation for fractional-order neural networks with mixed delays. Neurocomputing 2021. [Google Scholar] [CrossRef]
  28. Udhayakumar, K.; Rihan, F.A.; Rakkiyappan, R.; Cao, J.D. Fractional-order discontinuous systems with indefinite LKFs: An application to fractional-order neural networks with time delays. Neural Netw. 2022, 145, 319–330. [Google Scholar] [CrossRef]
  29. Wang, X.H.; Wang, Z.; Xia, J.W. Stability and bifurcation control of a delayed fractional-order eco-epidemiological model with incommensurate orders. J. Franklin Inst. 2019, 356, 8278–8295. [Google Scholar] [CrossRef]
  30. Li, C.P.; Dao, X.H.; Guo, P. Fractional derivatives in complex planes. Nonlinear Anal. TMA 2009, 71, 1857–1869. [Google Scholar] [CrossRef]
  31. Guariglia, E. Riemann zeta fractional derivative-functional equation and link with primes. Adv. Diff. Equat. 2019, 261, 15. [Google Scholar] [CrossRef]
  32. Ortigueira, M.D.; Coito, F. From differences to derivatives. Frac. Calc. Appl. Anal. 2004, 7, C459–C471. [Google Scholar]
  33. Guariglia, E.; Silvestrov, S. A functional equation for the Riemann zeta fractional derivative. AIP Conf. Proceed. 2017, 1798, 020063. [Google Scholar]
  34. Wu, G.C.; Lee, E.W.M. Fractional variational iteration method and its application. Phys. Lett. A 2010, 374, 2506–2509. [Google Scholar] [CrossRef]
  35. Xu, C.J.; Zhang, W.; Aouiti, C.; Liu, Z.X.; Liao, M.X.; Li, P.L. Further investigation on bifurcation and their control of fractional-order BAM neural networks involving four neurons and multiple delays. Math. Meth. Appl. Sci. 2021. [Google Scholar] [CrossRef]
  36. Xu, C.J.; Liu, Z.X.; Aouiti, C.; Li, P.L.; Yao, L.Y.; Yan, J.L. New exploration on bifurcation for fractional-order quaternion-valued neural networks involving leakage delays. Cogn. Neurodyn. 2022. [Google Scholar] [CrossRef]
  37. Xu, C.J.; Liao, M.X.; Li, P.L.; Yuan, S. Impact of leakage delay on bifurcation in fractional-order complex-valued neural networks. Chaos Solitons Fractals 2021, 142, 110535. [Google Scholar] [CrossRef]
  38. Xu, C.J.; Liu, Z.X.; Yao, L.Y.; Aouiti, C. Further exploration on bifurcation of fractional-order six-neuron bi-directional associative memory neural networks with multi-delays. Appl. Math. Comput. 2021, 410, 126458. [Google Scholar] [CrossRef]
  39. Huang, C.D.; Nie, X.B.; Zhao, X.; Song, Q.K.; Cao, J.D. Novel bifurcation results for a delayed fractional-order quaternion-valued neural network. Neural Netw. 2019, 117, 67–93. [Google Scholar] [CrossRef]
  40. Xiao, M.; Zheng, W.X.; Lin, J.X.; Jiang, G.P.; Zhao, L.D.; Cao, J.D. Fractional-order PD control at Hopf bifurcations in delayed fractional-order small-world networks. J. Franklin Inst. 2017, 354, 7643–7667. [Google Scholar] [CrossRef]
  41. Yan, Y.P. Hopf bifurcation and stability for a delayed tri-neuron neural networks. J. Comput. Appl. Math. 2006, 196, 579–595. [Google Scholar] [CrossRef] [Green Version]
  42. Xu, C.J.; Zhang, W.; Liu, Z.X.; Li, P.L.; Yao, L.Y. Bifurcation study for fractional-order three-layer neural networks involving four time delays. Cogn. Comput. 2022. [Google Scholar] [CrossRef]
  43. Song, Y.L.; Han, M.A.; Wei, J.J. Stability and Hopf bifurcation on a simplified BAM neural network with delays. Phys. D 2005, 200, 185–204. [Google Scholar] [CrossRef]
  44. Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  45. Bandyopadhyay, B.; Kamal, S. Stabliization and Control of Fractional Order Systems: A Sliding Mode Approach; Springer: Heidelberg, Germany, 2015; Volume 317. [Google Scholar]
  46. Li, Y.; Chen, Y.Q.; Podlubny, I. Stability of fractional-order nonlinear dynamic systems: Lyapunov direct method and generalized Mittag–Leffler stability. Comput. Math. Appl. 2009, 59, 1810–1821. [Google Scholar] [CrossRef] [Green Version]
  47. Li, H.L.; Zhang, L.; Hu, C.; Jiang, Y.L.; Teng, Z.D. Dynamical analysis of a fractional-order predator-prey model incorporating a prey refuge. J. Appl. Math. Comput. 2017, 54, 435–449. [Google Scholar] [CrossRef]
  48. Matignon, D. Stability results for fractional differential equations with applications to control processing. Comput. Eng. Syst. Appl. 1996, 2, 963–968. [Google Scholar]
  49. Deng, W.H.; Li, C.P.; Lü, J.H. Stability analysis of linear fractional differential system with multiple time delays. Nonlinear Dyn. 2007, 48, 409–416. [Google Scholar] [CrossRef]
  50. Yuan, J.; Zhao, L.Z.; Huang, C.D.; Xiao, M. A novel hybrid control technique for bifurcation in an exponential RED algorithm. Int. J. Circ. Theory Appl. 2020, 48, 1476–1492. [Google Scholar] [CrossRef]
Figure 1. Computer simulation figures of system (95) with ρ = 0.5 < ρ 0 0.6 . The three variables w 1 ( t ) 0 , w 2 ( t ) 0 , w 3 ( t ) 0 when time t .
Figure 1. Computer simulation figures of system (95) with ρ = 0.5 < ρ 0 0.6 . The three variables w 1 ( t ) 0 , w 2 ( t ) 0 , w 3 ( t ) 0 when time t .
Fractalfract 06 00161 g001aFractalfract 06 00161 g001b
Figure 2. Computer simulation figures of system (95) with ρ = 0.7 > ρ 0 0.6 . The Hopf bifurcation arises in the vicinity of W 0 ( 0 , 0 , 0 ) when time t .
Figure 2. Computer simulation figures of system (95) with ρ = 0.7 > ρ 0 0.6 . The Hopf bifurcation arises in the vicinity of W 0 ( 0 , 0 , 0 ) when time t .
Fractalfract 06 00161 g002aFractalfract 06 00161 g002b
Figure 3. Bifurcation figure of system (95): ρ versus w 1 . The bifurcation value ρ 0 0.6 .
Figure 3. Bifurcation figure of system (95): ρ versus w 1 . The bifurcation value ρ 0 0.6 .
Fractalfract 06 00161 g003
Figure 4. Bifurcation figure of system (95): ρ versus w 2 . The bifurcation value ρ 0 0.6 .
Figure 4. Bifurcation figure of system (95): ρ versus w 2 . The bifurcation value ρ 0 0.6 .
Fractalfract 06 00161 g004
Figure 5. Bifurcation figure of system (95): ρ versus w 3 . The bifurcation value ρ 0 0.6 .
Figure 5. Bifurcation figure of system (95): ρ versus w 3 . The bifurcation value ρ 0 0.6 .
Fractalfract 06 00161 g005
Figure 6. Computer simulation figures of system (96) with ρ = 0.24 < ρ 0.3 . The three variables w 1 ( t ) 0 , w 2 ( t ) 0 , w 3 ( t ) 0 when time t .
Figure 6. Computer simulation figures of system (96) with ρ = 0.24 < ρ 0.3 . The three variables w 1 ( t ) 0 , w 2 ( t ) 0 , w 3 ( t ) 0 when time t .
Fractalfract 06 00161 g006aFractalfract 06 00161 g006b
Figure 7. Computer simulation figures of system (96) with ρ = 0.35 > ρ 0.3 . The Hopf bifurcation arises in the vicinity of W 0 ( 0 , 0 , 0 ) when time t .
Figure 7. Computer simulation figures of system (96) with ρ = 0.35 > ρ 0.3 . The Hopf bifurcation arises in the vicinity of W 0 ( 0 , 0 , 0 ) when time t .
Fractalfract 06 00161 g007aFractalfract 06 00161 g007b
Figure 8. Bifurcation figure of system (96): ρ versus w 1 . The bifurcation value ρ 0.3 .
Figure 8. Bifurcation figure of system (96): ρ versus w 1 . The bifurcation value ρ 0.3 .
Fractalfract 06 00161 g008
Figure 9. Bifurcation figure of system (96): ρ versus w 2 . The bifurcation value ρ 0.3 .
Figure 9. Bifurcation figure of system (96): ρ versus w 2 . The bifurcation value ρ 0.3 .
Fractalfract 06 00161 g009
Figure 10. Bifurcation figure of system (96): ρ versus w 3 . The bifurcation value ρ 0.3 .
Figure 10. Bifurcation figure of system (96): ρ versus w 3 . The bifurcation value ρ 0.3 .
Fractalfract 06 00161 g010
Table 1. The relation of θ 0 and ρ 0 of system (95).
Table 1. The relation of θ 0 and ρ 0 of system (95).
θ 0 ρ 0
3.6893 0.15
2.3058 0.24
1.4563 0.38
1.2577 0.44
0.9380 0.59
0.9223 0.60
0.7095 0.78
0.6435 0.86
0.5705 0.97
Table 2. The relation of ϑ 0 and ρ of system (96).
Table 2. The relation of ϑ 0 and ρ of system (96).
ϑ 0 ρ
4.7675 0.12
2.4874 0.23
1.9070 0.30
1.3954 0.41
1.0216 0.56
0.9227 0.62
0.7837 0.73
0.6977 0.82
0.6086 0.94
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, P.; Yan, J.; Xu, C.; Shang, Y. Dynamic Analysis and Bifurcation Study on Fractional-Order Tri-Neuron Neural Networks Incorporating Delays. Fractal Fract. 2022, 6, 161. https://doi.org/10.3390/fractalfract6030161

AMA Style

Li P, Yan J, Xu C, Shang Y. Dynamic Analysis and Bifurcation Study on Fractional-Order Tri-Neuron Neural Networks Incorporating Delays. Fractal and Fractional. 2022; 6(3):161. https://doi.org/10.3390/fractalfract6030161

Chicago/Turabian Style

Li, Peiluan, Jinling Yan, Changjin Xu, and Youlin Shang. 2022. "Dynamic Analysis and Bifurcation Study on Fractional-Order Tri-Neuron Neural Networks Incorporating Delays" Fractal and Fractional 6, no. 3: 161. https://doi.org/10.3390/fractalfract6030161

APA Style

Li, P., Yan, J., Xu, C., & Shang, Y. (2022). Dynamic Analysis and Bifurcation Study on Fractional-Order Tri-Neuron Neural Networks Incorporating Delays. Fractal and Fractional, 6(3), 161. https://doi.org/10.3390/fractalfract6030161

Article Metrics

Back to TopTop