Next Article in Journal
Stability Properties of Multi-Term Fractional-Differential Equations
Previous Article in Journal
A Miniaturized Tri-Wideband Sierpinski Hexagonal-Shaped Fractal Antenna for Wireless Communication Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Hopf–Hopf Interactions Induced by Multiple Delays for Inertial Hopfield Neural Models

College of Mathematics and Information Science, Henan University of Economics and Law, Zhengzhou 450046, China
Fractal Fract. 2023, 7(2), 116; https://doi.org/10.3390/fractalfract7020116
Submission received: 29 November 2022 / Revised: 18 January 2023 / Accepted: 19 January 2023 / Published: 26 January 2023
(This article belongs to the Section Complexity)

Abstract

:
The investigation of dynamic behaviors of inertial neural networks depicted by second-order delayed differential equations has received considerable attention. Substantial research has been performed on the transformed first-order differential equations using traditional variable substitution. However, there are few studies on bifurcation dynamics using direct analysis. In this paper, a multi-delay Hopfield neural system with inertial couplings is considered. The perturbation scheme and non-reduced order technique are firstly combined into studying multi-delay induced Hopf–Hopf singularity. This combination avoids tedious computation and overcomes the disadvantages of the traditional variable-substitution reduced-order method. In the neighbor of Hopf–Hopf interaction points, interesting dynamics are found on the plane of self-connected delay and coupled delay. Multiple delays can induce the switching of stable periodic oscillation and periodic coexistence. The explicit expressions of periodic solutions are obtained. The validity of theoretical results is shown through consistency with numerical simulations.

1. Introduction

It is well-known that time delays are unavoidable due to the finite transmission speed of signals, system reaction, and processing time between synapses. Therefore, it is meaningful and interesting to focus on how time delays influence the dynamic behavior of natural systems. In the past few decades, numerous excellent results on dynamic analysis have been obtained for a natural system of first-order delayed differential equations such as stability [1,2,3,4], synchronization [5,6,7,8,9], bifurcation [10,11,12,13,14], and chaotic activity [15]. Artificial neural networks are very complex nonlinear dynamical systems. Hopf bifurcation is studied in fractional-order delayed neural models [16]. Recently, inertial neural models depicted using second-order delayed differential equations have received considerable attention because of their theoretical and practical significance [17,18,19]. The authors considered complex codimension-two bifurcations and coexistence dynamics for delayed neural networks with inertia terms [20,21,22,23,24,25,26]. The criteria for global synchronization and dissipativity are established on inertial neural networks with discrete-time or time-varying delays [27,28,29,30].
It is worth mentioning that some weaknesses need to be addressed in the existing research results. Most discussions of single delay-induced bifurcation dynamics employ central manifold analysis [31]. Despite being multiple-delay systems, authors usually convert them into single-delay systems under the limitation of delays for analysis convenience, which cannot well reveal the real and complex dynamic characteristics of a nonlinear neural system. To the best of our knowledge, for a system of second-order differential equations, studies are also always made on the rewritten first-order equations by introducing proper variable transformations. Some authors [32] pointed out that this substitution doubles the dimension of the addressed systems and greatly increases the complexity of theoretical analysis and the difficulty of studied results. Therefore, there is an urgent need for a new technique to avoid and overcome those questions arising from the traditional variable-substitution reduced-order method. To be more concise, a natural idea is to perform a dynamic analysis directly on the second-order differential systems. Very recently, authors began to develop non-reduced order strategies to reveal the global exponential stability or dissipativity, or state estimation for inertial network models with different kinds of delays [33,34,35,36,37]. Limit cycles have been investigated for second-order differential models with a single delay [38]. However, there are few studies on bifurcation dynamics based on the non-reduced method at present.
Motivated by the above discussion, and due to the high precision of artificial neural networks, we incorporate inertial couplings into a multi-delay neural system in this paper. We attempt to study the joint influences of multiple time delays on Hopf–Hopf interactions with the aid of the non-reduced order technique. The main works of this paper are as follows:
(1)
Diverse time delays are two major parameters of our concern. Compared with the previous related works, investigating the joint influences of diverse delays on inertial neural systems is more realistic and meaningful.
(2)
The perturbation scheme and non-reduced order technique are first combined into the study of Hopf–Hopf interactions. In contrast with the traditional reduced-order method, it is simple and valid with less computation.
(3)
The search for analytical bifurcating solutions is converted to the problem of solving four algebraic equations.
(4)
A developed model with inertial couplings and multiple delays is investigated and theoretical results demonstrate its validity.
This paper is arranged as follows. In Section 2, a perturbation scheme involving the non-reduced order technique is demonstrated to investigate Hopf–Hopf singularity induced by multiple delays. In Section 3, the Hopf–Hopf bifurcation of a multi-delay Hopfield inertial neural network is discussed using the method proposed in Section 2 and the numerical results are displayed to verify the accuracy of the theoretical analysis. Finally, a brief conclusion is given in Section 4.

2. Methodology Formulation

In this section, a perturbation scheme with the aid of the non-reduced order technique is presented to investigate Hopf–Hopf singularities near the bifurcation point when considering two delays as control parameters for a general second-order nonlinear differential system.

2.1. Hopf–Hopf Bifurcation Point

Consider a general delayed system of n -dimension second-order nonlinear equations with two discrete delays as follows
Z ¨ t = L Z ˙ t + E 0 Z t + E 1 Z t τ 1 + E 2 Z t τ 2 + ε F Z t , Z t τ 1 , Z t τ 2 ,
where Z t R n , Z ¨ t = d 2 Z t / d t 2 , Z ˙ t = d Z t / d t , L , E 0 , E 1 , and E 2 are real square matrices; ε is a nonlinear coupling strength; τ 1 and τ 2 are discrete delays; and a nonlinear smooth function F satisfies F 0 , 0 , , 0 = 0 .
The equilibrium point of Equation (1) is the zero solution 0 , 0 , , 0 . Linearization for Equation (1) near 0 , 0 , , 0 produces the following linear system
Z ¨ t = L Z ˙ t + E 0 Z t + E 1 Z t τ 1 + E 2 Z t τ 2 .
The corresponding characteristic equation is
det λ 2 I λ L E 0 e λ τ 1 E 1 e λ τ 2 E 2 = 0 .
Our main goals concern the joint effects of two diverse delays on Hopf–Hopf interactions in Equation (1). Therefore, to obtain our main results, the following assumptions first need to be made.
(a1)
Diverse delays τ 1 and τ 2 are chosen as two control parameters.
(a2)
Increasing two control delays, Equation (1) undergoes weak resonant Hopf–Hopf bifurcations at the critical values τ 1 = τ 1 c and τ 2 = τ 2 c . That is, the other roots of Equation (2) have negative real parts except for two distinct pairs of roots, ± i k 1 ω and ± i k 2 ω   ω > 0 for τ 1 = τ 1 c and τ 2 = τ 2 c where k 1 and k 2 are positive real constants.
The critical values τ 1 c and τ 2 c of Hopf–Hopf interactions can be obtained by studying the distribution of the roots of the associated characteristic Equation (2).

2.2. Bifurcation Sets and Periodic Solutions

In this section, the non-reduced order technique combined with a perturbation scheme is constructed to derive the explicit expressions of stable and unstable periodic solutions and sets of dynamical classifications in the neighborhood of the critical point τ 1 c , τ 2 c .
For control parameters τ 1 and τ 2 , there always exists a small perturbation such as
τ 1 = τ 1 c + δ 1 ε , τ 2 = τ 2 c + δ 2 ε .
Accordingly, Equation (1) is equivalent to being transformed as
Z ¨ = L Z ˙ + E 0 Z + E 1 Z τ 1 c + E 2 Z τ 2 c + F ^ ,
where
Z = Z t , Z τ 1 c = Z t τ 1 c , Z τ 2 c = Z t τ 2 c ,
F ^ = E 1 Z τ 1 Z τ 1 c + E 2 Z τ 2 Z τ 2 c + ε F Z , Z τ 1 , Z τ 2 , ε .
For ε = 0 in Equation (3), periodic solutions of the linearized equations are expressed by
Z 0 t = m = 1 2 a m 1 a m 2   a m n cos k m ω t + b m 1 b m 2   b m n sin k m ω t ,
where an integer n is the dimension of Equation (1).
For ε = 0 , taking Equation (4) to (3), we can solve the values of a m j and b m j m = 1 , 2 ;   j = 1 , 2 , , n in Equation (4) by satisfying the following equations
M m b m 1 b m 2   b m n = N m a m 1 a m 2   a m n , M m a m 1 a m 2   a m n = N m b m 1 b m 2   b m n ,
where
M m = k m ω L E 1 sin k m ω τ 1 c E 2 sin k m ω τ 2 c , N m = k m 2 ω 2 I E 0 E 1 cos k m ω τ 1 c E 2 cos k m ω τ 2 c .
For a small value ε   > 0 , the bifurcation solutions of the addressed Equation (3) can be seen as a small perturbation of periodic solutions in Equation (4), which are expressed in a polar coordinate by the following
Z t = m = 1 2 r m ε cos k m ω + σ m ε t + θ m z 2 t z n t ,
where r m 0 = r m , and θ m are determined by initial values.
The general expressions of z 2 t ,   , and z n t in the above, Equation (7), are too complicated and omitted here. For a specific differential system, the expressions of solutions are easy to simplify. Furthermore, it is not difficult to derive the values of the coefficients r m ε and perturbations σ m ε   m = 1 , 2 based on the adjoint operator. That is, firstly, according to the definition of the adjoint operator, one can obtain the adjoint system (8) in Lemma 1 of the original linear system (see Appendix A). Secondly, based on Lemma 1, r m ε and σ m ε   m = 1 , 2 can be solved from Theorem 1.
Lemma 1.
If Y t is the periodic solution of the following linear system
Y ¨ = L T Y ˙ + E 0 T Y + E 1 T Y t + τ 1 c + E 2 T Y t + τ 2 c ,
then the analytical expression Y t is given by
Y t = y 1 t y 2 t y n t = m = 1 2 p m 1 p m 2   p m n cos k i ω t + q m 1 q m 2   q m n sin k i ω t ,
with Y t = Y t + 2 π ω , where the unknown coefficients p m j and q m j m = 1 , 2 ;   j = 1 , 2 , , n are computed by the following equations
M m T p m 1 p m 2   p m n = N m T q m 1 q m 2   q m n , M m T q m 1 q m 2   q m n = N m T p m 1 p m 2   p m n ,
and M m and N m are expressed by Equation (6).
Theorem 1.
Based on Lemma 1, the values of r m ε and σ m ε m = 1 , 2 in perturbation solutions in Equation (7) are derived by
Y T 0 Z ˙ 2 π / ω Z ˙ 0 Y ˙ T 0 Z 2 π / ω Z 0 Y T 0 L Z 2 π / ω Z 0 τ 1 c 0 Y T t + τ 1 c E 1 Z t Z t + 2 π / ω d t τ 2 c 0 Y T t + τ 2 c E 2 Z t Z t + 2 π / ω d t 0 2 π / ω Y T t F ^ d t = 0 .
Proof.
Multiplying both sides of Equation (3) by Y T t and integrating on t , we have
0 2 π / ω Y T t Z ¨ d t = 0 2 π / ω Y T t L Z ˙ + E 0 Z + E 1 Z τ 1 c + E 2 Z τ 2 c + F ^ d t .
Due to Y t = Y t + 2 π / ω , we have the following equations through combination with the partial integral method
0 2 π / ω Y T t Z ¨ t d t = Y T 0 Z ˙ 2 π / ω Z ˙ 0 Y ˙ T 0 Z 2 π / ω Z 0 0 2 π / ω Y ˙ T t L Z d t + 0 2 π / ω Y T t E 0 Z d t + 0 2 π / ω Y T t + τ 1 c E 1 Z d t + 0 2 π / ω Y T t + τ 2 c E 2 Z d t ,
0 2 π / ω Y T t L Z ˙ d t = Y T 0 L Z 2 π / ω Z 0 0 2 π / ω Y ˙ T t L Z d t ,
0 2 π / ω Y T t E 1 Z τ 1 c d t = 0 2 π / ω Y T t + τ 1 c E 1 Z d t + τ 1 c 0 Y T t + τ 1 c E 1 Z Z t + 2 π / ω d t ,
0 2 π / ω Y T t E 2 Z τ 2 c d t = 0 2 π / ω Y T t + τ 2 c E 2 Z d t + τ 2 c 0 Y T t + τ 2 c E 2 Z Z t + 2 π / ω d t .
Equation (11) holds from Equations (13)–(16). □
Remark 1.
Equation (8) is the adjoint system of linearized system for the addressed Equation (1), which is obtained in Appendix A.
Remark 2.
The perturbation scheme and non-reduced order technique are first combined to discuss bifurcation dynamics. The incorporation of the non-reduced order technique can be considered as a complement of bifurcation investigation on the second-order differential model with the traditional reduced-order method [14,20].
Remark 3.
In contrast to the perturbation scheme with the traditional reduced order technique in [14], the main difference is the increase in two items in the key Equation (11) in Theorem 1, i.e., Y T 0 Z ˙ 2 π / ω Z ˙ 0 and Y ˙ T 0 Z 2 π / ω Z 0 . However, the two terms are very simple to calculate because the expressions of Z , Z ˙ and Y , Y ˙ are obtained easily by the perturbation solution (7) and periodic solution (9), respectively. The non-reduced order technique introduced in this paper avoids and overcomes the above-mentioned problem of twice raising the dimension arising from the traditional reduced order technique. Furthermore, it has a very clear procedure. Therefore, it can be easier to program the computation of the periodic solution, using a direct analysis, than the reduced order technique in [14].
Remark 4.
Equation (11) can yield four algebraic equations for r 1 ε , r 2 ε , σ 1 ε , and σ 2 ε . When bifurcation parameters are in the neighborhood of Hopf–Hopf bifurcation points (i.e., δ 1 ε and δ 2 ε are very small), it is not difficult to analytically derive periodic solutions in Equation (7) from them.

3. Multiple-Delay Neural Networks with Inertial Couplings

In this section, we employ the methodology in Section 2 to investigate multiple delay-induced Hopf–Hopf bifurcations for a pair of neurons’ multi-delay systems with inertial couplings.

3.1. Existence of Hopf–Hopf Bifurcation Point

Due to the high-precision description of artificial neural networks, we incorporate inertial couplings into the dynamical model [26]. The developed model with inertial couplings and multiple delays can be taken into account
M 1 x ¨ t = m 1 x ˙ t K 1 x t τ 1 + c 1 tanh y t τ 2 , M 2 y ¨ t = m 2 y ˙ t K 2 y t τ 1 + c 2 tanh x t τ 2 ,
where M i and m i   i = 1 , 2 are positive real constants; K i i = 1 , 2 denotes self-connection weights from one neuron to itself; c i i = 1 , 2 are the cross-interaction weights through two neurons; and discrete delays τ 1 and τ 2 are non-negative real constants due to some physical meanings.
The characteristic equation of the linearization system for (17) around the trivial equilibrium point is determined by
λ 4 + λ 3 m 1 M 1 + m 2 M 2 + λ 2 K 1 M 1 e λ τ 1 + K 2 M 2 e λ τ 2 + m 1 m 2 M 1 M 2 + λ K 2 m 1 e λ τ 1 M 1 M 2 + K 1 m 2 e λ τ 2 M 1 M 2 + K 1 K 2 e 2 λ τ 1 M 1 M 2 c 1 c 2 e 2 λ τ 2 M 1 M 2 = 0 .
For K 1 K 2 c 1 c 2 (i.e., λ 0 ), substituting λ = ± i ω ω > 0 into the above characteristic Equation (18) and separating real and imaginary parts yields the following equations
cos 2 ω τ 2 = K 1 K 2 cos 2 ω τ 1 + K 2 m 1 + K 1 m 2 ω sin ω τ 1 c 1 c 2 + ω 4 M 1 M 2 ω 2 m 1 m 2 K 2 M 1 + K 1 M 2 ω 2 cos ω τ 1 c 1 c 2 , sin 2 ω τ 2 = K 1 K 2 sin 2 ω τ 1 K 2 m 1 + K 1 m 2 ω cos ω τ 1 c 1 c 2 + ω 3 m 2 M 1 + m 1 M 2 K 2 M 1 + K 1 M 2 ω 2 sin ω τ 1 c 1 c 2 .
Eliminating τ 2 from cos 2 2 ω τ 2 + sin 2 2 ω τ 2 = 1 in (19), we obtain the following equation on ω
k 1 2 k 2 2 c 1 2 c 2 2 + ω 8 M 1 2 M 2 2 2 ω k 1 k 2 m 1 k 2 + k 1 m 2 sin ω τ 1 + ω 2 2 cos 2 ω τ 1 k 1 k 2 m 1 m 2 + k 2 m 1 + k 1 m 2 2 2 cos ω τ 1 k 1 k 2 k 2 M 1 + k 1 M 2 2 ω 3 sin ω τ 1 m 1 m 2 k 2 m 1 + k 1 m 2 + sin 2 ω τ 1 k 1 k 2 m 2 M 1 + m 1 M 2 + ω 4 2 cos ω τ 1 k 1 m 2 2 M 1 + 2 cos 2 ω τ 1 k 1 k 2 M 1 M 2 + ω 4 k 2 M 1 + k 1 M 2 2 + m 1 2 m 2 2 2 cos ω τ 1 k 2 M 2 + ω m 1 2 M 2 2 ω 5 + ω 5 2 k 1 sin ω τ 1 m 1 M 2 2 2 k 1 ω cos ω τ 1 M 1 M 2 2 + ω m 2 2 M 1 2
2 k 2 M 1 2 sin ω τ 1 m 2 + ω cos ω τ 1 M 2 ω 5 = 0 .
The transversality condition on Hopf bifurcations occurring should be satisfied. That is, Re λ τ 2 λ = ± i ω 0 , where
λ τ 2 = 2 c 1 c 2 λ e 2 τ 1 λ / G ,
and
G = 2 c 1 c 2 τ 2 e 2 τ 1 λ + K 2 m 1 + K 1 m 2 e 2 τ 2 + τ 1 λ 2 K 1 K 2 τ 1 e 2 τ 2 λ + 2 m 1 m 2 λ e 2 τ 2 + τ 1 λ + 2 K 2 m 1 + K 1 m 2 λ e 2 τ 2 + τ 1 λ K 2 m 1 + K 1 m 2 τ 1 λ e 2 τ 2 + τ 1 λ + 3 m 2 M 1 λ 2 e 2 τ 2 + τ 1 λ K 2 m 1 + K 1 m 2 τ 1 λ 2 e 2 τ 2 + τ 1 λ + 4 M 1 M 2 λ 3 e 2 τ 2 + τ 1 λ .
If Equation (18) has two pairs of purely imaginary roots λ 1 = ± i ω 1 and λ 2 = ± i ω 2 ω 1 > ω 2 > 0 , then the critical delays are expressed by
τ 2 1 , j = ϕ 1 + 2 j π 2 ω 1 , j = 0 , 1 , 2 , ,   ϕ 1 0 , 2 π
where
cos ϕ 1 = K 1 K 2 cos 2 ω 1 τ 1 + K 2 m 1 + K 1 m 2 ω 1 sin ω 1 τ 1 c 1 c 2 +   ω 1 4 M 1 M 2 K 2 M 1 + K 1 M 2 ω 1 2 cos ω 1 τ 1 ω 1 2 m 1 m 2 c 1 c 2 , sin ϕ 1 = K 1 K 2 sin 2 ω 1 τ 1 K 2 m 1 + K 1 m 2 ω 1 cos ω 1 τ 1 c 1 c 2 + ω 1 3 m 2 M 1 + m 1 M 2 K 2 M 1 + K 1 M 2 ω 1 2 sin ω 1 τ 1 c 1 c 2 .
and
τ 2 2 , j = ϕ 2 + 2 j π 2 ω 2 , j = 0 , 1 , 2 , ,   ϕ 2 0 , 2 π
where
cos ϕ 1 = K 1 K 2 cos 2 ω 2 τ 1 + K 2 m 1 + K 1 m 2 ω 2 sin ω 2 τ 1 c 1 c 2 +   ω 2 4 M 1 M 2 K 2 M 1 + K 1 M 2 ω 2 2 cos ω 2 τ 1 ω 2 2 m 1 m 2 c 1 c 2 , sin ϕ 1 = K 1 K 2 sin 2 ω 2 τ 1 K 2 m 1 + K 1 m 2 ω 2 cos ω 2 τ 1 c 1 c 2 + ω 2 3 m 2 M 1 + m 1 M 2 K 2 M 1 + K 1 M 2 ω 2 2 sin ω 2 τ 1 c 1 c 2 .
It is well-known that double-Hopf bifurcation points are the intersection points of two Hopf bifurcation curves. Due to the existence of multiple delays in the neural system (17), double-Hopf bifurcation points cannot be derived analytically. However, for given parameter values, bifurcation points can be easily obtained.
For example, some values of parameters are chosen as follows
M 1 = M 2 = 1.0 ,   K 1 = K 2 = 1.0 ,   m 1 = 1.0 ,   m 2 = 1.6 ,   c 1 = 0.8 ,   c 2 = 0.2 .
Then, Equation (18) has only two pairs of simple purely imaginary roots ± i 8 ω and ± i 11 ω , with ω = 0.080575 in red and all the other eigenvalues with negative real parts in green, in Figure 1. Therefore, we can obtain that Hopf–Hopf bifurcation occurs at τ 1 = τ 1 c = 1.1383 and τ 2 = τ 2 c = 4.2847 .

3.2. Bifurcation Sets and Numerical Simulations

To apply the mentioned non-reduced order method, let some perturbations on bifurcation parameters
τ 1 = τ 1 c + δ 1 ε 2 ,   τ 2 = τ 2 c + δ 2 ε 2 ,
and variable transformation
x t ε x t ,   y t ε y t .
System (17) becomes Equation (1)
L = m 1 0 0 m 2 , E 0 = 0 0 0 0 , E 1 = 1 0 0 1 ,
E 2 = 0 c 1 c 2 0 , F = ε 3 c 1 y 3 t τ 1 + h . o . t ε 3 c 2 x 3 t τ 2 + h . o . t .
Based on Equations (5) and (6), the perturbation solution of system (17) close to Hopf–Hopf bifurcation points can be expressed as
Z t = ε r 1 cos 8 ω + ε 2 σ 1 t + θ 1 + ε r 2 cos 11 ω + ε 2 σ 2 t + θ 2 0.368219 ε r 1 cos 8 ω + ε 2 σ 1 t + θ 1 0.180657 ε r 1 sin 8 ω + ε 2 σ 1 t + θ 1 + 0.280884 ε r 2 cos 11 ω + ε 2 σ 2 t + θ 2 0.152974 ε r 2 sin 11 ω + ε 2 σ 2 t + θ 2 .
Using Lemma 1, the adjoint periodic solutions in Equation (9) is written as
Y t = p 11 cos 8 ω t + q 11 sin 8 ω t + p 21 cos 11 ω t + q 21 sin 11 ω t p 12 cos 8 ω t + q 12 sin 8 ω t + p 22 cos 11 ω t + q 22 sin 11 ω t ,
where
p 11 = 0.547394 p 12 + 0.268365 q 12 , q 11 = 0.268365 p 12 0.547394 q 12 ,  
p 21 = 0.686302 p 22 + 0.373805 q 22 , q 21 = 0.373805 p 22 + 0.686302 q 22 .
Putting Equations (20) and (21) into Equation (11), omitting the higher order on, four algebraic equations on ε 2 r 1 , ε 2 r 2 , ε 2 σ 1 , and ε 2 σ 2 are derived as
ε 2 r 1 0.843989 r 1 2 1.59277 r 2 2 + 18.5664 δ 1 9.33851 δ 2 + 69.1896 σ 1 = 0 , ε 2 r 1 2.11533 r 1 2 + 3.99186 r 2 2 13.773 δ 1 3.72481 δ 2 + 94.2979 σ 1 = 0 , ε 2 r 2 3.60921 r 1 2 1.70279 r 2 2 32.3505 δ 1 8.43129 δ 2 + 144.122 σ 2 = 0 , ε 2 r 2 2.77835 r 1 2 + 1.31079 r 2 2 11.3497 δ 1 10.9523 δ 2 + 7.97716 σ 2 = 0 ,
where
ε 2 σ 1 = 0.122374 δ 1 ε 2 + 0.101342 δ 2 ε 2 , ε 2 σ 2 = 0.304847 δ 1 ε 2 + 0.146674 δ 2 ε 2 .
Substituting the perturbed expressions of time delays τ 1 = τ 1 c + ε 2 δ 1 , τ 2 = τ 2 c + ε 2 δ 2 in Equation (22), except for the solution 0 , 0 , we can easily obtain the other three solutions r 1 , r 2 depending on two parameters τ 1 and τ 2 , as follows
1.66036 0.65625 + 4.34064 τ 1 τ 2 , 0   for   τ 2 < 0.65625 + 4.34064 τ 1 ,
0 , 2.73183 5.32241 + 0.911631 τ 1 + τ 2   for   τ 2 > 5.32241 0.911631 τ 1 ,
2.36928 4.34369 + 0.0518228 τ 1 + τ 2 ,   2.10601 2.69693 + 1.39486 τ 1 τ 2
for   4.34369 0.0518228 τ 1 < τ 2 < 2.69693 + 1.39486 τ 1 .
In addition, four critical bifurcation lines are given as follows:
l 1 :   τ 2 = 4.34064 τ 1 0.65625 , l 2 :   τ 2 = 5.32241 0.911631 τ 1 ,
l 3 :   τ 2 = 4.34369 0.0518228 τ 1 , l 4 :   τ 2 = 1.39486 τ 1 + 2.69693 .
The bifurcation diagram on the plane τ 1 , τ 2 is given close to the Hopf–Hopf interaction points in Figure 2, where the bifurcation point (1.1383,4.2847) is indicated by the intersection black dotted point. The delayed parameter plane is divided into six regions denoted by I, II, III, IV, V, and VI, which are divided by the line l i   i = 1 , 2 , 3 , 4 in red and their dynamic classifications depend on two delays. The interesting phenomena are exhibited involving stable equilibrium, the switching of periodic motions, and the coexistence of periodic solutions.
(1)
τ 1 , τ 2 region (I), the trivial equilibrium of Equation (17) is asymptotically stable. In addition, Equation (17) undergoes a Hopf bifurcation at the line l 1 . Region (I) is a stability zone.
(2)
τ 1 , τ 2 region (II), 0 , 0 loses its stability and a stable periodic oscillation emerges.
(3)
τ 1 , τ 2 region (III), when τ 1 , τ 2 crosses the line l 2 , there exist two periodic solutions and the trivial equilibrium point. A periodic solution has high frequency and the trivial equilibrium point is unstable and the other periodic solution with low frequency is stable.
(4)
τ 1 , τ 2 region (IV), when τ 1 , τ 2 crosses the line l 3 , there exist three periodic solutions and the trivial equilibrium point. The new emerging periodic solution and the trivial equilibrium are unstable.
(5)
τ 1 , τ 2 region (V), when τ 1 , τ 2 crosses the line l 4 , the unstable periodic solution disappears. The periodic solution with high frequency is stable and the other periodic solution and the trivial equilibrium point are unstable.
(6)
τ 1 , τ 2 region (VI), there exists a periodic solution with high frequency and a trivial equilibrium point. The periodic solution bifurcating from the trivial equilibrium point is stable.
Next, we plot some stable-state solutions using the software Matlab. Firstly, the comparison of analytical periodic solutions from the reduced-order method [14] and the non-reduced order method from this paper is displayed in Figure 3, using time-response trajectories (left) and phase portraits (right) with 1.15 , 4.2 case   ( II ) .
Secondly, for τ 1 , τ 2 = 1.0 , 4.3 region   ( I ) , time histories and all the eigenvalues of the characteristic equation (roots with negative real parts denoted by green) are given to verify a stable trivial equilibrium point in Figure 4. Figure 5 shows the multi-stability coexistence on periodic solutions, where τ 1 , τ 2 = 1.25 , 4.31 region ( IV ) . Finally, Figure 6 exhibits the stability switching of period oscillation from low frequency to high frequency corresponding to four bifurcation sets, where τ 1 , τ 2 = 1.15 , 4.2 region   ( II ) , 1.25 , 4.25 region   ( II I ) , 1.25 , 4.35 region   ( V ) and 1.1 , 4.45 region   ( VI ) , respectively.
It is clear to see that our analytical results given by the non-reduced order technique are in excellent agreement with those displayed using the reduced-order method [14]. All the numerical results also show that the perturbation scheme in combination with the non-reduced order method is effective and easier for discussing bifurcation analytical solutions and dynamical classification in general second-order differential systems.

4. Conclusions

In this article, in order to clearly reveal the real and complex dynamic characteristics of neural systems, a multi-delay neural model with inertial couplings is considered. Inertia and multiple delays can produce richer nonlinear dynamics, including the switching of periodic oscillation and periodic coexistence. The perturbation scheme and non-reduced order technique are first combined to analyze the joint influences of diverse delays on Hopf–Hopf interactions. Results from this combination are in good agreement with those from the traditional reduced-order method. Furthermore, it has a very clear procedure and can easily be programmed to derive dynamical sets and explicit bifurcation solutions in closed form in the neighborhood of the critical point.
The inclusion of the non-reduced order technique avoids and overcomes the disadvantages of the traditional variable substitution reduced order method. These advantages have been verified in comparison with the previous related studies. The key research findings of our work have significant theoretical guiding value in dominating and optimizing networks.
To the best of our knowledge, bifurcation dynamics involving the non-reduced order technique are discussed for the first time in this paper. Further research will consider absolute synchronization based on the non-reduced order technique for a system of second-order delayed differential equations.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 11872175), the Henan Natural Science Foundation for Outstanding Youth (Grant No. 212300410021), and Young Talents Fund of HUEL.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are very grateful to the editors and anonymous reviewers for their constructive comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The adjoint operator K * of K is given by
K Z , Y = Z , K * Y
where
Z , Y = 0 b Y T Z d t ,   b is   a   real   constant .
For Equation (3), let
K Z = Z ¨ L Z ˙ E 0 Z E 1 Z τ 1 c E 2 Z τ 2 c .
Choosing b = 2 π / ω , and substituting Equation (A2) into the left of Equation (A1), one can obtain
K Z , Y = 0 2 π / ω Y T t Z ¨ L Z ˙ E 0 Z E 1 Z τ 1 c E 2 Z τ 2 c d t .
According to Y t = Y t + 2 π / ω , we have the following
0 2 π / ω Y T Z ¨ d t = 0 2 π / ω Y ¨ T Z d t ,   0 2 π / ω Y T L Z ˙ d t = 0 2 π / ω Y ˙ T L Z d t , 0 2 π / ω Y T E 1 Z τ 1 c d t = 0 2 π / ω Y T t + τ 1 c E 1 Z d t , 0 2 π / ω Y T E 2 Z τ 2 c d t = 0 2 π / ω Y T t + τ 2 c E 2 Z d t .
Substituting Equation (A4) into Equation (A3), we obtain
K Z , Y = 0 2 π / ω Y ¨ t + L T Y ˙ t E 0 T Y t E 1 T Y t + τ 1 c E 2 T Y t + τ 2 c T Z d t .
Based on Equation (A1), Equation (8) is derived.

References

  1. Wang, F.; Liu, M. Global exponential stability of high-order bidirectional associative memory (BAM) neural networks with time delays in leakage terms. Neurocomputing 2016, 177, 515–528. [Google Scholar] [CrossRef]
  2. Dong, Y.; Takeuchi, Y.; Nakaoka, S. A mathematical model of multiple delayed feedback control system of the gut microbiota-Antibiotics injection controlled by measured metagenomic data. Nonlinear Anal. Real 2018, 43, 1–17. [Google Scholar] [CrossRef]
  3. Yang, J.; Liu, C.; Mei, M. Global solutions for bistable degenerate reaction-diffusion equation with time-delay and nonlocal effect. Appl. Math. Lett. 2022, 125, 107726. [Google Scholar] [CrossRef]
  4. Mao, X.; Wang, Z. Stability, bifurcation, and synchronization of delay-coupled ring neural networks. Nonlinear Dyn. 2016, 84, 1063–1078. [Google Scholar] [CrossRef]
  5. Shen, H.; Hu, X.; Wang, J.; Cao, J.; Qian, W. Non-fragile H∞ synchronization for Markov jump singularly perturbed coupled neural networks subject to double-layer switching regulation. IEEE Trans. Neural Netw. Learn Syst. 2022, 1–11. [Google Scholar] [CrossRef] [PubMed]
  6. Cao, Y.; Zhao, L.; Wen, S.; Huang, T. Lag H∞ synchronization of coupled neural networks with multiple state couplings and multiple delayed state couplings. Neural Netw. 2022, 151, 143–155. [Google Scholar] [CrossRef] [PubMed]
  7. Ge, J.; Xu, J. Computation of synchronized periodic solution in a BAM network with two delays. IEEE Trans. Neural Netw. Learn Syst. 2010, 21, 439–450. [Google Scholar]
  8. Shi, K.; Tang, Y.; Liu, X.; Zhong, S. Non-fragile sampled-data robust synchronization of uncertain delayed chaotic Lurie systems with randomly occurring controller gain fluctuation. ISA T. 2017, 66, 185–199. [Google Scholar] [CrossRef]
  9. Shi, K.; Liu, X.; Zhu, H.; Zhong, S.; Zeng, Y.; Yin, C. Novel delay-dependent master-slave synchronization criteria of chaotic Lurie systems with time-varying delay feedback control. App. Math. Comput. 2016, 282, 137–154. [Google Scholar] [CrossRef]
  10. He, X.; Li, C.; Huang, T.; Li, C. Bogdanov-Takens Singularity in Tri-Neuron Network with Time Delay. IEEE Trans. Neural Netw. Learn Syst. 2013, 24, 1001–1007. [Google Scholar]
  11. Gupta, P.; Majee, N.; Roy, A. Stability and Hopf bifurcation analysis of delayed BAM neural network under dynamic thresholds. Nonlinear Anal. Model. 2009, 14, 435–461. [Google Scholar] [CrossRef] [Green Version]
  12. Ma, S. Hopf bifurcation of a type of neuron model with multiple time delays. Int. J. Bifurcat. Chaos 2020, 29, 1950163. [Google Scholar] [CrossRef]
  13. Pei, L.; Wang, S. Double Hopf bifurcation of differential equation with linearly state-dependent delays via MMS. Appl. Math. Comput. 2019, 341, 256–276. [Google Scholar] [CrossRef]
  14. Ge, J.; Xu, J. An analytical method for studying double Hopf bifurcations induced by two delays in nonlinear differential systems. Sci. China Technol. Sci. 2020, 63, 597–602. [Google Scholar] [CrossRef]
  15. Yu, D.; Zhou, X.; Wang, G.; Ding, Q.; Li, T.; Jia, Y. Effects of chaotic activity and time delay on signal transmission in FitzHugh-Nagumo neuronal system. Cogn. Neurodyn. 2022, 16, 887–897. [Google Scholar] [CrossRef]
  16. Zhao, L.; Huang, C.; Cao, J. Effects of double delays on bifurcation for a fractional-order neural network. Cogn. Neurodyn. 2022, 16, 1189–1201. [Google Scholar] [CrossRef]
  17. He, X.; Li, C.; Huang, T.; Yu, J. Bifurcation behaviors of an Euler discretized inertial delayed neuron model. Sci. China Technol. Sci. 2016, 59, 418–427. [Google Scholar] [CrossRef]
  18. He, X.; Huang, T.; Yu, J.; Li, C.; Li, C. An inertial projection neural network for solving variational inequalities. IEEE Trans. Cybern. 2017, 47, 809–814. [Google Scholar] [CrossRef]
  19. Ge, J.; Xu, J. Stability and Hopf bifurcation on four-neuron neural networks with inertia and multiple delays. Neurocomputing 2018, 287, 34–44. [Google Scholar] [CrossRef]
  20. Song, Z.; Wang, C.; Zhen, B. Codimension-two bifurcation and multistability coexistence in an inertial two-neuron system with multiple delays. Nonlinear Dyn. 2016, 85, 2099–2113. [Google Scholar] [CrossRef]
  21. Song, Z.; Xu, J. Self-/mutual-symmetric rhythms and their coexistence in a delayed half-center oscillator of the CPG neural system. Nonlinear Dyn. 2022, 108, 2595–2609. [Google Scholar] [CrossRef]
  22. Ge, J.; Xu, J. Stability switches and fold-Hopf bifurcations in an inertial four-neuron network model with coupling delay. Neurocomputing 2013, 110, 70–79. [Google Scholar] [CrossRef]
  23. Song, Z.; Zhen, B.; Hu, D. Multiple bifurcations and coexistence in an inertial two-neuron system with multiple time delays. Cogn. Neurodyn. 2020, 14, 359–374. [Google Scholar] [CrossRef]
  24. Ge, J.; Xu, J. Weak resonant double Hopf bifurcations in an inertial four-neuron model with time delay. Int. J. Neural Syst. 2012, 22, 63–75. [Google Scholar] [CrossRef]
  25. Yao, S.; Ding, L.; Song, Z.; Xu, J. Two bifurcation routes to multiple chaotic coexistence in an inertial two-neural system with time delay. Nonlinear Dyn. 2019, 95, 1549–1563. [Google Scholar] [CrossRef]
  26. Ge, J. Multi-delay-induced bifurcation singularity in two-neuron neural models with multiple time delays. Nonlinear Dyn. 2022, 108, 4357–4371. [Google Scholar] [CrossRef]
  27. Lakshmanan, S.; Prakash, M.; Lim, C.; Rakkiyappan, R.; Balasubramaniam, P.; Nahavandi, S. Synchronization of an inertial neural network with time-varying delays and its application of secure communication. IEEE Trans. Neural Netw. Learn Syst. 2018, 29, 195–207. [Google Scholar] [CrossRef] [PubMed]
  28. Guo, Z.; Gong, S.; Yang, S.; Huang, T. Global exponential synchronization of multiple coupled inertial memristive neural networks with time-varying delay via nonlinear coupling. Neural Netw. 2018, 108, 260–271. [Google Scholar] [CrossRef] [PubMed]
  29. Chen, X.; Lin, D.; Lan, W. Global dissipativity of delayed discrete-time inertial neural networks. Neurocomputing 2020, 390, 131–138. [Google Scholar] [CrossRef]
  30. Cui, Q.; Li, L.; Cao, J. Stability of inertial delayed neural networks with stochastic delayed impulses via matrix measure method. Neurocomputing 2022, 471, 70–78. [Google Scholar] [CrossRef]
  31. Hassard, B.D.; Kazarinoff, N.D.; Wan, Y. Theory and Applications of Hopf bifurcation; Cambridge University Press: Cambridge, UK, 1981. [Google Scholar]
  32. Li, X.; Li, X.; Hu, C. Some new results on stability and synchronization for delayed inertial neural networks based on non-reduced order method. Neural Netw. 2017, 96, 91–100. [Google Scholar] [CrossRef]
  33. Huang, C.; Liu, B. New studies on dynamic analysis of inertial neural networks involving non-reduced order method. Neurocomputing 2019, 325, 283–287. [Google Scholar] [CrossRef]
  34. Wu, K.; Jian, J. Non-reduced order strategies for global dissipativity of memristive neural-type inertial neural networks with mixed time-varying delays. Neurocomputing 2021, 436, 174–183. [Google Scholar] [CrossRef]
  35. Chen, S.; Jiang, H.; Lu, B.; Yu, Z.; Li, L. Pinning bipartite synchronization for inertial coupled delayed neural networks with signed digraph via non-reduced order method. Neural Netw. 2020, 129, 392–402. [Google Scholar] [CrossRef] [PubMed]
  36. Wang, J.; Wang, X.; Wang, Y.; Zhang, X. Non-reduced order method to global h-stability criteria for proportional delay high-order inertial neural networks. Appl. Math. Comput. 2021, 407, 126308. [Google Scholar] [CrossRef]
  37. Tu, Z.; Dai, N.; Wang, L.; Yang, X.; Wu, Y.; Li, N.; Cao, J. H state estimation of quaternion-valued inertial neural networks: Non-reduced order method. Cogn. Neurodyn. 2022. [Google Scholar] [CrossRef]
  38. Liu, P.; Kalmár-Nagy, T. High-dimensional harmonic balance analysis for second-order delay-differential equations. J. Vib. Control 2010, 16, 1189–1208. [Google Scholar] [CrossRef]
Figure 1. The eigenvalues of characteristic equation for networks (17).
Figure 1. The eigenvalues of characteristic equation for networks (17).
Fractalfract 07 00116 g001
Figure 2. Bifurcation sets on τ 1 versus τ 2 and dynamics near Hopf–Hopf interaction point of inertial models (17). The black intersection dotted point is the bifurcation point (1.1383,4.2847) in the portrait. The delayed plane is divided into six regions by the red line l i   i = 1 , 2 , 3 , 4 . (I) A stable equilibrium; (II) and (III) periodic motion with low frequency; (IV) coexistence of periodic oscillations; (V) and (VI) periodic motion with high frequency.
Figure 2. Bifurcation sets on τ 1 versus τ 2 and dynamics near Hopf–Hopf interaction point of inertial models (17). The black intersection dotted point is the bifurcation point (1.1383,4.2847) in the portrait. The delayed plane is divided into six regions by the red line l i   i = 1 , 2 , 3 , 4 . (I) A stable equilibrium; (II) and (III) periodic motion with low frequency; (IV) coexistence of periodic oscillations; (V) and (VI) periodic motion with high frequency.
Fractalfract 07 00116 g002
Figure 3. Comparison of the analytical periodic solutions derived from non-reduced order via reduced order technique.
Figure 3. Comparison of the analytical periodic solutions derived from non-reduced order via reduced order technique.
Fractalfract 07 00116 g003
Figure 4. Time−response trajectories of X and Y versus t and the eigenvalues of characteristic equation for networks (17) with 1.0 , 4.3 case   ( I ) corresponding to the stable zero solution.
Figure 4. Time−response trajectories of X and Y versus t and the eigenvalues of characteristic equation for networks (17) with 1.0 , 4.3 case   ( I ) corresponding to the stable zero solution.
Fractalfract 07 00116 g004
Figure 5. Multi−stability periodic coexistence of network (17) with different initial conditions exhibited by time-response trajectories of X and Y versus t and three phase portraits for 1.25 , 4.31 case   ( IV ) .
Figure 5. Multi−stability periodic coexistence of network (17) with different initial conditions exhibited by time-response trajectories of X and Y versus t and three phase portraits for 1.25 , 4.31 case   ( IV ) .
Fractalfract 07 00116 g005
Figure 6. Time−response trajectories of X versus t for networks (17) corresponding to four bifurcation sets of Figure 1 denoted by II, III, V and VI where time delays are chosen as 1.15 , 4.2 case   ( II ) , 1.25 , 4.25 case   ( II I ) , 1.25 , 4.35 case   ( V ) , and 1.1 , 4.45 case   ( VI ) .
Figure 6. Time−response trajectories of X versus t for networks (17) corresponding to four bifurcation sets of Figure 1 denoted by II, III, V and VI where time delays are chosen as 1.15 , 4.2 case   ( II ) , 1.25 , 4.25 case   ( II I ) , 1.25 , 4.35 case   ( V ) , and 1.1 , 4.45 case   ( VI ) .
Fractalfract 07 00116 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ge, J. Analysis of Hopf–Hopf Interactions Induced by Multiple Delays for Inertial Hopfield Neural Models. Fractal Fract. 2023, 7, 116. https://doi.org/10.3390/fractalfract7020116

AMA Style

Ge J. Analysis of Hopf–Hopf Interactions Induced by Multiple Delays for Inertial Hopfield Neural Models. Fractal and Fractional. 2023; 7(2):116. https://doi.org/10.3390/fractalfract7020116

Chicago/Turabian Style

Ge, Juhong. 2023. "Analysis of Hopf–Hopf Interactions Induced by Multiple Delays for Inertial Hopfield Neural Models" Fractal and Fractional 7, no. 2: 116. https://doi.org/10.3390/fractalfract7020116

APA Style

Ge, J. (2023). Analysis of Hopf–Hopf Interactions Induced by Multiple Delays for Inertial Hopfield Neural Models. Fractal and Fractional, 7(2), 116. https://doi.org/10.3390/fractalfract7020116

Article Metrics

Back to TopTop