Analysis and control of the fractional chaotic Hopfield neural network

The fractional Hopfield neural network (HNN) model is studied here analyzing its symmetry, uniqueness of the solution, dissipativity, fixed points etc. A Lyapunov and bifurcation analysis of the system is done for specific as well as variable fractional order. Since a very long time ago, HNN has been carefully studied and applied in various fields. Because of the exceptional non-linearity of the neuron activation function, the HNN system is stoutly non-linear. Chaos control using adaptive SMC considering disturbances and uncertainties is done about randomly chosen points by designing suitable controllers. Numerical simulations performed in MATLAB verify the efficacy of the designed controllers.


Introduction
Various mathematical models have been proposed to understand the various phenomena better. Proposing biological models [1] has become most important as it helps scientists in better understanding and bringing new insight into it, such as epidemic modeling that may help in control of the epidemic. Most of the biological models [2,3] that have been proposed fall in the category of non-linear biological systems. Chaos, antimonotonicity, different types of bifurcations, multi-stability are the many properties [4] which have been visualized in many biological models. In order to analyze the proposed nonlinear dynamical models [5][6][7][8], one should use non-linear methods such as Lyapunov exponents, stagnation points analysis, basins of attraction, and bifurcation diagrams. In order to model the changing dynamics, the considered parameter values can be changed over a wide range.
Since a very long time ago, the Hopfield neural network (HNN) has been carefully studied and applied in various fields such as in image encryption, data storage, information processing, and associative memory. Because of the exceptional non-linearity of the neuron activation function, the HNN model is stoutly non-linear [9]. Like the chaotic time delay systems, Chua circuit, and coupled HR neuron circuits [10], the HNN model generates very complex behaviors like hyper-chaos, periodic chaos and quasi-period [11], which implies that the HNN model simulates the prominent chaotic behavior like that of the brain [12,13]. Therefore, one must investigate these systems as they have a theoretical signifi-cance as well as practical significance. Up to now various HNN models have been studied and proposed in the literature; such as the two dimensional neuronal model [14], three dimensional neuronal model [15] and four dimensional neuronal model [16]. For these systems a variety of studies has been done, such as multi-stability, coexisting attractors, circuit implementation, and synchronization [17][18][19][20][21][22][23].
Many chaos control methods [24][25][26][27] have been developed in the recent past to tame chaos and increase its application across various disciplines. Some of the popular methods used are active control, tracking control, sliding mode, and the parameter estimation method. Chaos synchronization is also a way to contain chaos between master and slave systems. Some popular synchronization methods use the anti-synchronization, projective synchronization, difference synchronization, and matrix synchronization [28] methods [29][30][31][32].
We study fractional [33,34] HNN model varying parameter values and fractional orders. Some basic dynamic methods [35][36][37] have been used such as Lyapunov exponents [38,39], bifurcation diagrams [40], phase portraits, and time series. We have also discussed the existence of the solution of the considered system. As the considered system did not have any stagnation point, we have used chaos control using adaptive sliding mode about two arbitrarily chosen desired points considering external disturbances and uncertainties. The disturbances have been estimated and the error converging to zero has been achieved, which have been plotted using MATLAB software.

The fractional Hopfield neural network
The fractional HNN model is [41]: For P = 0.8 and I.C. (0, 0.01, 0), system (1) shows chaotic behavior. Here we introduce the fractional version of HNN, perform its thorough dynamical analysis and chaos control. Caputo's derivative is used in the paper: The fractional HNN model is given as where Z = (Z 1 , Z 2 , Z 3 ) T ∈ R 3 are the state variables and P ∈ R is a parameter value.

Dynamics of fractional HNN model
The dynamics of the fractional HNN with tanhyperbolic terms is explored here, studying the symmetry, dissipative, uniqueness of the solution, bifurcation and Lyapunov dynamics, fixed point analysis, etc.

Symmetry, stagnation point analysis and dissipativity
System (2) in matrix form can be written as The fractional HNN model (2) does not remain invariant under the Z i → -Z i , Z j → -Z j , Z k → Z k transformation i.e. the system possesses asymmetric behavior about all the axes. However, the system is symmetric about the origin as under the transformation Z 1 → -Z 1 , Z 2 → -Z 2 , Z 3 → -Z 3 , the system remains invariant. The divergence of G is Therefore (2) is dissipative. Equating G i (Z 1 , Z 2 , Z 3 ) for i = 1, 2, 3 to 0, the system can be explored for stagnation points i.e.
For P = 0.8 we obtain no stagnation points. The absence of stagnation points hints at the complex nature of the chaotic system.

Solution of HNN model
Theorem The I. V.P. of system (2) , then G(Z) is continuous and bounded on [Z 0 -, Z 0 + ] for any > 0, therefore Lipschitz continuity over [Z 0 -, Z 0 + ] proves the existence and uniqueness of the solution.
We have We have Thus ) and W (t) ∈ R 3 . Hence system (2) possesses a unique solution.  From the formula

Lyapunov dynamics and bifurcation
where p is such that

Controlling chaos
Chaos in fractional HNN model exposed to uncertainties and disturbances is controlled using an adaptive SMC technique. Suitably designed controllers are constructed to stabilize chaos in the trajectories of the system about an arbitrarily chosen point (p 1 , p 2 , p 3 ). The fractional HNN model exposed to uncertainties and external disturbances is where H i are uncertainties and D i are disturbances, v i are controllers designed about desired point. Figures 6 and 7 give the trajectories and plots of the exposed system. Consider | H i | and D i to be bounded by positive values C i and F i withĈ i ,F i being their estimates.  The sliding surface is defined as To have (5) in sliding mode, the necessary condition is Differentiating (6): Then from (7), we have Equation (9) is stable using Matignon's theorem [42]. The designed controllers are with sign(·), the signum function. Parameter update conditions arė with c i , f i > 0 are constants.

Figure 7
Chaotic attractors of exposed system (3) Proof The proof is based on Lyapunov's direct method, defining [43] the Lyapunov function by where Differentiating (13): From (8), we havė Substituting D q e i ,Ċ i andḞ i in (15): Finally, Thus ∃ a real T ≥ 0 so that From Lyapunov stability theory s i → 0 as t → ∞. Hence the errors converge to s i = 0 implying stability about the desired point.

Conclusion
Dynamical properties of fractional order HNN model are studied in the paper by analyzing the system's symmetry, uniqueness of the solution, dissipativity and fixed points.
Lyapunov dynamics and bifurcations of the system are studied for specific order as well as variable order. The chaos in a fractional system exposed to uncertainties and disturbances is contained by designing suitable controllers based on an adaptive sliding mode control technique about two arbitrarily chosen desired points. The simulations performed in Matlab have been displayed and discussed. Synchronization of the fractional HNN model with some other system involves future scope of work in this direction. Also studying the system as regards its hidden attractors as well as its electronic circuit would be interesting.