Robust Synchronization of Delayed Chaotic FitzHugh-Nagumo Neurons under External Electrical Stimulation

Abstract

Synchronization of chaotic neurons under external electrical stimulation (EES) is studied in order to understand information processing in the brain and to improve the methodologies employed in the treatment of cognitive diseases. This paper investigates the dynamics of uncertain coupled chaotic delayed FitzHugh-Nagumo (FHN) neurons under EES for incorporated parametric variations. A global nonlinear control law for synchronization of delayed neurons with known parameters is developed. Based on local and global Lipschitz conditions, knowledge of the bounds on the neuronal states, the Lyapunov-Krasovskii functional, and the 𝐿₂ gain reduction, a less conservative local robust nonlinear control law is formulated to address the problem of robust asymptotic synchronization of delayed FHN neurons under parametric uncertainties. The proposed local control law guarantees both robust stability and robust performance and provides the 𝐿₂ bound for uncertainty rejection in the synchronization error dynamics. Separate conditions for single-input and multiple-input control schemes for synchronization of a wide class of FHN systems are provided. The results of the proposed techniques are verified through numerical simulations.

1. Introduction

The neuron is the fundamental unit of the functioning brain [1]. Its dynamical investigation, for the purpose of measuring brain activity and understanding how the neural system transmits electrochemical signals to the muscles, is one of the most significant challenges facing brain researchers [2–8]. Neural-system malfunctions can cause potentially fatal motor-function-impairment diseases such as Parkinson’s, Huntington’s, and epilepsy [9–11]. Synchronization of chaotic neurons with gap junctions under external electrical stimulation (EES), one of the most fundamental research issues, has been relentlessly studied in order to improve therapy-based treatments of neurodegenerative disorders [12, 13]. In this respect, the FitzHugh-Nagumo (FHN) model, which has been applied in other fields (e.g., chemical reaction kinetics [14]) as well, is one of the most pertinent neural models utilized in synchronization studies [15–18].

Given the FHN model’s wide applicability, many of its significant and complex dynamical aspects, including chaos, bifurcation, synchronization, control, noise effects and filtering, coupling, and medium effects, not to mention disturbance rejection, have already been reviewed extensively in the literature [12–20]. Researchers have applied nonlinear, adaptive, fuzzy, neural-network-, and observer-based, as well as robust control methodologies to the synchronization of FHN neurons under EES [12–14, 16, 21–24]. However, these traditional conservative synchronization controller design techniques [12–14, 16, 21–24], which in mathematical models ignore the time-delay arising due to the separation between coupled neurons, cannot synchronize distant FHN neurons. The dynamics of coupled chaotic delayed FHN neurons with gap junctions under EES recently have been investigated [17, 18], which can be accounted for synchronization studies.

In the present work, we examined, preparatory to a brief numerical simulation study, the dynamics of coupled delayed FHN neurons in consideration of the model parametric variations. This parametric-variation-based model, with its separate visualization of each uncertain component, offers better insight into dynamical uncertainty in actual neurons; as such, it is a superior means of control law formulation for neuronal synchronization. Control law derivation for synchronization of delayed FHN neurons separated by gap junctions under EES has remained extremely rare to this date. In the interests of filling this research gap, we propose a global synchronization control strategy for identical neurons with known parameters. To that end, we also propose a novel local robust control law that guarantees asymptotic convergence of synchronization errors to zero under time-delays and parametric uncertainties. This regional control methodology, which is based on local and global Lipschitz constraints on nonlinear and uncertain components of neuronal dynamics, knowledge of state bounds, and the Lyapunov-Krasovskii (LK) functional, is less conservative in its performance within the desired locality. Additionally, we developed conditions under which robust control law performance bounds, which distinguish our work and enable the choice of a suitable robust controller, can be determined. The proposed computationally simple control strategy, with its easy design procedure, ensures both robust stability and robust performance in neuronal synchronization. Details on the robust single-input and multiple-input control strategies are presented herein to facilitate their application to a wide class of FHN systems. Finally, the proposed control schemes are successfully validated by numerical simulations. The main contributions of this paper are as follows. (i)The dynamics of coupled delayed FHN neurons under parametric uncertainties are studied in order to provide better insight into uncertain coupled neurons.(ii)Global synchronization control of delayed FHN neurons, with guaranteed convergence of synchronization errors to zero, is achieved. (iii)A less conservative nonlinear control law for local robust synchronization of delayed FHN neurons under parametric uncertainties is developed that ensures asymptotic convergence of synchronization errors to zero.(iv)A robust performance assessment tool, in terms of the robustness bound, is provided for evaluation of the performance of a local controller.(v)Single-input and multiple-input synchronization control laws that select specific and different objective functions for robust performance in their respective case, are derived in order to broaden the scope of the proposed schemes.

This paper is organized as follows. Section 2 provides a brief overview of two identical coupled chaotic FHN neurons with gap junctions, delays, and parametric uncertainties. Section 3 proposes a nonlinear control law for global synchronization of neurons with delays. Section 4 addresses the issue of local robust synchronization of delayed uncertain FHN neurons. Section 5 provides and discusses the simulation results for the proposed schemes. Section 6 draws conclusions.

Standard notations are used in this paper. The 𝐿₂ norm of a vector 𝑧 is defined as ‖𝑧‖₂=(∫^∞₀‖𝑧‖²𝑑𝑡)^1/2, where ‖𝑧‖ denotes the Euclidian norm of 𝑧. A positive definite symmetric matrix 𝑋 is denoted as 𝑋>0. For 𝑥_𝑖 with the 𝑖th diagonal entry and 𝑖=1,2,…,𝑛, diag(𝑥₁,𝑥₂,…,𝑥_𝑛) represents a diagonal matrix.

2. Model Description

Consider two identical uncertain coupled chaotic FHN neurons with time-delays under EES: 𝑑𝑥₁𝑑𝑡=𝑥₁𝑥₁−11−(𝑟+Δ𝑟)𝑥₁−𝑦₁−𝑔+Δ𝑔𝑥₁−𝑥₂(𝑡−𝜏)+𝑠₁(𝑡),𝑑𝑦₁=()𝑥𝑑𝑡𝑏+Δ𝑏₁−(𝑣+Δ𝑣)𝑦₁,𝑑𝑥₂𝑑𝑡=𝑥₂𝑥₂−11−(𝑟+Δ𝑟)𝑥₂−𝑦₂−𝑔+Δ𝑔𝑥₂−𝑥₁(𝑡−𝜏)+𝑠₂(𝑡),𝑑𝑦₂=()𝑥𝑑𝑡𝑏+Δ𝑏₂−(𝑣+Δ𝑣)𝑦₂, (1) where 𝑥 and 𝑦 represent the states of a neuron in terms of activation potential and recovery voltage, respectively; (𝑥₁,𝑦₁) and (𝑥₂,𝑦₂) represent the states of the master and slave FHN neurons, respectively; 𝑔 describes the strength of the gap junctions between neurons; the parameter 𝜏>0 indicates the time-delay due to separation between neurons. The master and the slave neurons are under external electrical stimulation with currents 𝑠₁(𝑡)=(𝑎₁/𝜔₁)cos𝜔₁𝑡 and 𝑠₂(𝑡)=(𝑎₂/𝜔₂)cos𝜔₂𝑡, respectively. Here, 𝜔₁=2𝜋𝑓₁ and 𝜔₂=2𝜋𝑓₂ are dimensionless angular frequencies, 𝑎₁/𝜔₁ and 𝑎₁/𝜔₁ are stimulation amplitudes, and 𝑡 denotes time.

In biological models, we often know only the estimated or nominal values of parameters, not the exact or true values. Thus, in synchronizing identical neurons, the exact values of the model parameters, unlike the cases in scenarios examined in previous studies [12, 13, 21, 22], are unknown. In contrast to the literature [12–18, 21, 22], the terms Δ𝑔, Δ𝑟, Δ𝑏, and Δ𝑣, representing the parametric variations in 𝑔, 𝑟, 𝑏, and 𝑣, respectively, are added to the FHN neurons (1). Another reality of biological models is that the two neurons cannot be at all identical. Certainly, synchronization of identical coupled chaotic systems [25, 26], such as delayed FHN neurons under parametric uncertainties, has remained as a complex, challenging, and nontrivial problem. This problem can be resolved by utilization of adaptive control that is computationally complex [27, 28], and indeed especially complex when adaptation laws are required for a number of parametric variations. In order to resolve this problem, we now address the issue of robust nonlinear controller design. The model parameters selected are 𝑟=10,𝑔=0.1,𝑓₁=𝑓₂𝑎=0.138,𝑏=1,₁=𝑎₂=0.1,𝜏=20,𝑣=0.1 (2) along with the parametric uncertainties Δ𝑔=0.1,Δ𝑟=1,Δ𝑏=−0.1,Δ𝑣=−0.01. (3)

Figure 1 shows phase portraits and synchronization error plots for FHN neurons. These chaotic neurons are not synchronous, as the synchronization errors in Figures 1(a) and 1(b) do not converge to zero. We develop a control strategy utilizing two control inputs, 𝑢₁ and 𝑢₂, for synchronization of the slave neuron with the master neuron. The FHN model with control inputs, then, is given by 𝑑𝑥₁𝑑𝑡=𝑥₁𝑥₁−11−(𝑟+Δ𝑟)𝑥₁−𝑦₁−𝑔+Δ𝑔𝑥₁−𝑥₂(𝑡−𝜏)+𝑠₁(𝑡),𝑑𝑦₁=()𝑥𝑑𝑡𝑏+Δ𝑏₁−(𝑣+Δ𝑣)𝑦₁,𝑑𝑥₂𝑑𝑡=𝑥₂𝑥₂−11−(𝑟+Δ𝑟)𝑥₂−𝑦₂−𝑔+Δ𝑔𝑥₂−𝑥₁(𝑡−𝜏)+𝑠₂(𝑡)+𝑢₁,𝑑𝑦₂=()𝑥𝑑𝑡𝑏+Δ𝑏₂−(𝑣+Δ𝑣)𝑦₂+𝑢₂. (4) Now, the design of control inputs 𝑢₁ and 𝑢₂ for synchronization of the FHN neurons will be addressed.

(a)

(b)

(c)

(d)

Figure 1 

Chaotic behavior of nonsynchronous coupled delayed FHN neurons with parametric variations under EES. (a) Phase portrait of 𝑥1 and 𝑦1; (b) phase portrait of 𝑥2 and 𝑦2; (c) synchronization error 𝑒1=𝑥1−𝑥2 versus 𝑡; (d) synchronization error 𝑒2=𝑦1−𝑦2 versus 𝑡.

3. Global Nonlinear Control

The traditional control techniques have not clearly addressed synchronization of delayed FHN neurons, not even in the absence of uncertainty. Indeed, delayed neuronal models recently have been analyzed [17, 18], and their synchronization remains an open problem. We therefore develop a global nonlinear control strategy to address this issue. The proposed nonlinear control law for the synchronization of both FHN neurons is given by 𝑢₁=𝐶₀𝑥₁−𝑥₂+𝐶₁𝑥₁(𝑡−𝜏)−𝑥₂()−()𝑥𝑡−𝜏1+𝑟²₂−𝑟𝑥³₂+(1+𝑟)𝑥²₁−𝑟𝑥³₁+𝑠₁(𝑡)−𝑠₂(𝑡),𝑢₂=𝐶₂𝑦₁−𝑦₂, (5) where 𝐶₀, 𝐶₁, and 𝐶₂ are the controller gains. The term 𝐶₀(𝑥₁−𝑥₂) in the control law is required for convergence for the activation potential of the slave neuron to the corresponding state of the master neuron. The delayed term 𝐶₁(𝑥₁(𝑡−𝜏)−𝑥₂(𝑡−𝜏)) provides an extra degree of freedom for convergence of the activation potential. Identical behavior of recovery variables of the two neurons can be achieved by application of 𝐶₂(𝑦₁−𝑦₂). The nonlinear components −((1+𝑟)𝑥²₂−𝑟𝑥³₂) and ((1+𝑟)𝑥²₁−𝑟𝑥³₁) are used to cancel the effect of known nonlinear parts in the synchronization error dynamics, which is required to simplify the controller design (see, e.g, [23, 24] and references therein). To overcome the effect of difference in stimulation signals for the master and the slave FHN neurons, the term 𝑠₁(𝑡)−𝑠₂(𝑡) is used in the control law (5).

Remark 1. The proposed control law is computationally simpler than the so-called traditional synchronization techniques [12, 13, 21–23], despite the fact that it deals with FHN neuron dynamics that are more complex owing to delays and parametric uncertainties (discussed below). Moreover, to further reduce control workload and memory utilization, the parameters 𝐶₁ and 𝐶₂ can be selected as zero. These parameters are included in the control law to make it more general and, thus, applicable to other complex chaotic systems as well. Selection of 𝐶₂=0 makes the controller input 𝑢₂=0, which is desired for synchronization of FHN neurons. The techniques presented in this study, then, are useful for designing either a single-input or a multiple-input controller according to the given requirement.

To address the synchronization of neurons in the absence of uncertainties, we make the following assumption.

Assumption 2. The parametric uncertainties Δ𝑔, Δ𝑟, Δ𝑏, and Δ𝑣 are zero.

The following theorem provides the sufficient condition for global synchronization of delayed chaotic FHN neurons.

Theorem 3. Suppose that the FHN neurons (4) satisfying Assumption 2, the nonlinear control law given by (5), synchronizes the coupled neurons asymptotically with proper selection of parameters 𝐶₀, 𝐶₁, and 𝐶₂, if the matrix inequalities 𝐴𝑃>0,𝑄>0,Ω=^𝑇𝑃+𝑃𝐴+𝑄𝑃𝐴₁∗−𝑄<0 (6) are satisfied, where −𝐴=1+𝐶₀+𝑔−1𝑏−𝑣+𝐶₂,𝐴₁=−𝐶₁0.+𝑔00 (7)

Proof. Using the control law (5) in (4) under Assumption 2, the overall closed-loop system becomes 𝑑𝑥₁𝑑𝑡=𝑥₁𝑥₁−11−𝑟𝑥₁−𝑦₁−𝑔𝑥₁−𝑥₂(𝑡−𝜏)+𝑠₁(𝑡),𝑑𝑦₁𝑑𝑡=𝑏𝑥₁−𝑣𝑦₁,𝑑𝑥₂𝑑𝑡=−𝑟𝑥³₁+(1+𝑟)𝑥²₁−𝑥₂−𝑦₂−𝑔𝑥₂−𝑥₁(𝑡−𝜏)+𝑠₁(𝑡)+𝐶₀𝑥₁−𝑥₂+𝐶₁𝑥₁(𝑡−𝜏)−𝑥₂(𝑡−𝜏),𝑑𝑦₂𝑑𝑡=𝑏𝑥₂−𝑣𝑦₂+𝐶₂𝑦₁−𝑦₂. (8) Defining the synchronization errors as 𝑒₁=𝑥₁−𝑥₂,𝑒₂=𝑦₁−𝑦₂, (9) the error dynamics are 𝑑𝑒₁𝑑𝑡=−1+𝐶₀+𝑔𝑒₁−𝑒₂−𝐶₁+𝑔𝑒₁(𝑡−𝜏),𝑑𝑒₂𝑑𝑡=𝑏𝑒₁−𝑣+𝐶₂𝑒₂. (10) The error dynamics can be written as 𝑑𝑒𝑑𝑡=𝐴𝑒+𝐴₁𝑒(𝑡−𝜏), (11) where 𝑒=[𝑒₁𝑒₂]^𝑇, 𝑒(𝑡−𝜏)=[𝑒₁(𝑡−𝜏)𝑒₂(𝑡−𝜏)]^𝑇, and matrices 𝐴 and 𝐴₁ are given by (7). Now consider the Lyapunov functional 𝑉(𝑒,𝑡)=𝑒^𝑇𝑃𝑒+^𝑡_𝑡−𝜏𝑒^𝑇(𝜁)𝑄𝑒(𝜁)𝑑𝜁, (12) such that 𝑃>0, 𝑄>0. The derivative of (12) along (11) is 𝑉̇(𝑒,𝑡)=𝑒^𝑇𝐴^𝑇𝑃+𝑃𝐴+𝑄𝑒+𝑒^𝑇(𝑡−𝜏)𝐴^𝑇₁𝑃𝑒+𝑒^𝑇𝑃𝐴₁𝑒(𝑡−𝜏)−𝑒^𝑇(𝑡−𝜏)𝑄𝑒(𝑡−𝜏). (13) For stability, 𝑉̇(𝑒,𝑡)<0, that is, Ω<0, which completes the proof of Theorem 3.

Remark 4. Traditional studies on control and synchronization of FHN models for a variety of applications ignored the factor of time-delays, which can produce unrealistic outcomes. In the present study, the consideration of time-delays due to separation between FHN systems such as neurons makes the synchronization problem more realistic than the classical schemes [12, 13, 21–23].

Remark 5. The proposed nonlinear control strategy utilizing matrix inequalities, derived by employing the LK functional, guarantees global asymptotic synchronization of delayed FHN neurons. These inequalities can be resolved, for a selection of control parameters, by using linear matrix inequality (LMI)-based tools. This makes the computation of matrices 𝑃 and 𝑄 easier and hence helpful for selection of control parameters. It is also possible to incorporate additional LMI constraints for better performance, for instance to fix the upper bound on the synchronization error convergence rate.

4. Local Nonlinear Control

In the previous section, all parameters of FHN neurons were assumed to be known, with zero variations. In this section, we derive a sufficient condition for the synchronization of both FHN neurons under time-delays and parametric variations. Additionally, we provide sufficient conditions for the robust performance of single-input and multiple-input controllers. It is not always necessary to synchronize two identical oscillators globally. If the regional bounds on the states of the oscillators are known, a local controller can be a better choice. In reality, local controllers are less conservative, due to easy management of performance, robustness, and computation reduction for a given specific region (rather than emphasizing the whole space) [14, 29, 30]. Before delving into design methodology, we will review some basic definitions from the literature [31, 32].

Definition 6. The 𝐿₂ gain of a system from signal 𝑑 to 𝑒 is said to be less than a positive scalar 𝛾 if ‖𝑒‖₂<𝛾‖𝑑‖₂+𝛽, where 𝛽 is a small positive constant [31].

Definition 7. A function 𝑓(𝑥) is said to be Lipschitz if it satisfies the Lipschitz condition 𝑓(𝑥)−𝑓(𝑥)≤‖𝐿(𝑥−𝑥)‖, (14) where 𝑥,𝑥,𝐿∈𝑅. Moreover, the Lipschitz nonlinearities also satisfy 𝜕𝑓(𝑥≤‖𝐿‖,𝜕𝑥) (15) which is an inequality useful in determining 𝐿, by application of numerical algorithms.

Definition 8. A function 𝑓(𝑥) is said to be locally Lipschitz for 𝑥∈𝑅 if it satisfies the conditions (14)-(15) locally for a bounded region {𝑥,𝑥∈[𝑥_min,𝑥_max]∶𝑥_min,𝑥_max∈𝑅}, where 𝑥_min and 𝑥_max are the minimum and maximum limits on 𝑥 (or 𝑥), respectively.
From Definition 8, it can be seen that −𝑥³₁+𝑥²₁+𝑥³₂−𝑥²₂≤𝐿_𝑎𝑥₁−𝑥₂,∀𝑥₁,𝑥₂∈𝑥_min,𝑥_max, (16) where 𝐿_𝑎 is the Lipschitz constant for the local region, which can be selected by solving (15). To address the issue of the synchronization of neurons under parametric uncertainties, we make the following assumption.

Assumption 9. The parametric uncertainties are bounded by ‖Δ𝑔‖≤𝑔_𝑚, ‖Δ𝑏‖≤𝑏_𝑚, ‖Δ𝑟‖≤𝑟_𝑚, and ‖Δ𝑣‖≤𝑣_𝑚.
It is clear from (16) that the nonlinearity present in the FHN model satisfies the local Lipschitz condition. It has been reported previously [14, 29, 32] that a local controller can be designed for locally Lipschitz nonlinear systems if the bounds on the states of the system are known. Thanks to case studies, it is in fact well known that the states of a real neuron are always bounded in terms of the limits on activation potential and recovery voltage. This fact can be observed also in the simulation results shown in Figure 1 (see [12, 13, 17, 18, 21–24] as well). Therefore, by incorporating the knowledge of the minimum and maximum values of the states of neurons, the idea of bounds 𝑔_𝑚, 𝑏_𝑚, 𝑟_𝑚, and 𝑣_𝑚, and by noting the fact that the nonlinear part of the dynamics is locally Lipschitz, a regional robust controller can be designed for synchronization of FHN neurons.

Theorem 10. Suppose that the FHN neurons (4) satisfy Assumption 9 with states bounded by a region 𝑥₁,𝑥₂∈[𝑥_min,𝑥_max]. The nonlinear control law (5) synchronizes the coupled neurons asymptotically with proper selection of parameters 𝐶₀, 𝐶₁, and 𝐶₂, if the following matrix inequalities are verified: 𝑃_𝑠>0,𝑄_𝑠Γ>0,_𝑠=𝐴^𝑇𝑃+𝑃𝐴+𝑄+𝐹₁𝑃𝐴₁𝑃𝑃𝑃∗−𝑄+𝐹₂000∗∗−𝐼00∗∗∗−𝐼0∗∗∗∗−𝐼<0, (17) where 𝐹₁=𝑟²_𝑚𝐿²_𝑎+𝑔²_𝑚+𝑏²_𝑚00𝑣²_𝑚,𝐹₂=𝑔²_𝑚000, (18) and matrices 𝐴 and 𝐴₁ are given by (7).

Proof. Incorporating the control law (5) into (4), the overall closed-loop system becomes 𝑑𝑥₁𝑑𝑡=𝑥₁𝑥₁−11−(𝑟+Δ𝑟)𝑥₁−𝑦₁−𝑔+Δ𝑔𝑥₁−𝑥₂(𝑡−𝜏)+𝑠₁(𝑡),𝑑𝑦₁=()𝑥𝑑𝑡𝑏+Δ𝑏₁−(𝑣+Δ𝑣)𝑦₁,𝑑𝑥₂𝑑𝑡=−𝑟𝑥³₁+(1+𝑟)𝑥²₁−𝑥₂−𝑦₂−×𝑥𝑔+Δ𝑔₂−𝑥₁(𝑡−𝜏)+𝑠₁(𝑡)+𝐶₀𝑥₁−𝑥₂+𝐶₁𝑥₁(𝑡−𝜏)−𝑥₂(𝑡−𝜏)−Δ𝑟𝑥³₂+Δ𝑟𝑥²₂,𝑑𝑦₂=()𝑥𝑑𝑡𝑏+Δ𝑏₂−(𝑣+Δ𝑣)𝑦₂+𝐶₂𝑦₁−𝑦₂. (19) Using the same procedure as discussed in the previous section, the following error dynamics model is obtained: 𝑑𝑒𝑑𝑡=𝐴𝑒+𝐴₁𝑒(𝑡−𝜏)+Ψ+Θ+Φ, (20) where Ψ=𝜓𝑥₁,𝑦₁−𝜓𝑥₂,𝑦₂,Θ=𝜃𝑥₁−𝜃𝑥₂,Φ=𝜙𝑥₁(𝑡−𝜏)−𝜙𝑥₂(),𝜓𝑥𝑡−𝜏₁,𝑦₁=−Δ𝑟𝑥³₁−𝑥²₁−Δ𝑔𝑥₁−Δ𝑣𝑦₁,𝜓𝑥₂,𝑦₂=−Δ𝑟𝑥³₂−𝑥²₂−Δ𝑔𝑥₂−Δ𝑣𝑦₂,𝜃𝑥₁=0Δ𝑏𝑥₁,𝜃𝑥₂=0Δ𝑏𝑥₂,𝜙𝑥₁(𝑡−𝜏)=−Δ𝑔𝑥₁()0,𝜙𝑥𝑡−𝜏₂(𝑡−𝜏)=−Δ𝑔𝑥₂()0.𝑡−𝜏 (21) On the basis of Assumption 9 and inequality (16), we have Δ𝑟−𝑥³₁+𝑥²₁+𝑥³₂−𝑥²₂≤𝑟_𝑚𝐿_𝑎𝑥₁−𝑥₂,∀𝑥₁,𝑥₂∈𝑥_min,𝑥_max,Δ𝑔−𝑥₁+𝑥₂≤𝑔_𝑚𝑥₁−𝑥₂. (22) Combining the local and global Lipschitz constraints of (22), we have Δ𝑟−𝑥³₁+𝑥²₁+𝑥³₂−𝑥²₂+Δ𝑔−𝑥₁+𝑥₂≤𝑟²_𝑚𝐿²_𝑎+𝑔²_𝑚𝑥₁−𝑥₂,∀𝑥₁,𝑥₂∈𝑥_min,𝑥_max. (23) As we know from the global Lipschitz condition, Δ𝑏𝑥₁−𝑥₂≤𝑏_𝑚𝑥₁−𝑥₂,Δ𝑣−𝑦₁+𝑦₂≤𝑣_𝑚𝑦₁−𝑦₂,Δ𝑔−𝑥₁(𝑡−𝜏)+𝑥₂()≤𝑔𝑡−𝜏_𝑚𝑥₁(𝑡−𝜏)−𝑥₂(𝑡−𝜏). (24)(25)(26) Inequalities (23)–(25) imply Ψ^𝑇Ψ+Θ^𝑇Θ≤𝑒^𝑇𝐹₁𝑒, (27) and (26) implies Φ^𝑇Φ≤𝑒^𝑇(𝑡−𝜏)𝐹₂𝑒(𝑡−𝜏). (28) Constructing the LK functional 𝐸(𝑒,𝑡)=𝑒^𝑇𝑃_𝑠𝑒+^𝑡_𝑡−𝜏𝑒^𝑇(𝜁)𝑄_𝑠𝑒(𝜁)𝑑𝜁, (29) with 𝑃_𝑠>0 and 𝑄_𝑠>0.
Taking the derivative of (29) along (20) 𝐸̇(𝑒,𝑡)=𝑒^𝑇𝐴^𝑇𝑃_𝑠+𝑃_𝑠𝐴+𝑄_𝑠𝑒+𝑒^𝑇(𝑡−𝜏)×𝐴^𝑇₁𝑃_𝑠𝑒+𝑒^𝑇𝑃_𝑠𝐴₁𝑒(𝑡−𝜏)−𝑒^𝑇(𝑡−𝜏)×𝑄_𝑠𝑒(𝑡−𝜏)+Ψ^𝑇𝑃_𝑠𝑒+𝑒^𝑇𝑃_𝑠Ψ+Θ^𝑇𝑃_𝑠𝑒+𝑒^𝑇𝑃_𝑠Θ+Φ^𝑇𝑃_𝑠𝑒+𝑒^𝑇𝑃_𝑠Φ, (30) and, using inequalities (27)-(28), we obtain 𝐸̇(𝑒,𝑡)≤𝑒^𝑇𝐴^𝑇𝑃_𝑠+𝑃_𝑠𝐴+𝑄_𝑠𝑒+𝑒^𝑇(𝑡−𝜏)𝐴^𝑇₁𝑃_𝑠𝑒+𝑒^𝑇𝑃_𝑠𝐴₁𝑒(𝑡−𝜏)−𝑒^𝑇(𝑡−𝜏)𝑄_𝑠𝑒(𝑡−𝜏)+Ψ^𝑇𝑃_𝑠𝑒+𝑒^𝑇𝑃_𝑠Ψ+Θ^𝑇𝑃_𝑠𝑒+𝑒^𝑇𝑃_𝑠Θ+Φ^𝑇𝑃_𝑠𝑒+𝑒^𝑇𝑃_𝑠Φ−Ψ^𝑇Ψ−Θ^𝑇Θ−Φ^𝑇Φ+𝑒^𝑇𝐹₁𝑒+𝑒^𝑇(𝑡−𝜏)𝐹₂𝑒(𝑡−𝜏). (31) This further implies that 𝐸̇(𝑒,𝑡)≤𝜉^𝑇Γ_𝑠𝜉, (32) where 𝜉^𝑇=𝑒^𝑇𝑒^𝑇(𝑡−𝜏)Ψ^𝑇Θ^𝑇Φ^𝑇𝑇. (33) For stability, 𝐸̇(𝑒,𝑡)<0, and hence Γ_𝑠<0, which completes the proof of Theorem 10.

Remark 11. It is notable that Theorem 10 also guarantees global asymptotic synchronization in the absence of parametric uncertainties. By ignoring Φ, Ψ, and Θ in (31), corresponding to Δ𝑔=0, Δ𝑟=0, Δ𝑏=0, and Δ𝑣=0, one can obtain the matrix inequalities in (6).

Theorem 10 provides a sufficient condition for the local robust asymptotic synchronization of uncertain delayed FHN neurons, ensuring zero synchronization errors in the steady state. The other pertinent issue is robust performance in terms of the 𝐿₂ gain reduction from the uncertain nonlinearities Ψ, Θ, and Φ to the error 𝑒. By selecting a controller with a smaller size of the error 𝑒 with respect to the uncertainties Ψ, Θ, and Φ, the required robust performance can be achieved. For this purpose, again consider system (20), but in an alternative form given by 𝑑𝑒𝑑𝑡=𝐴𝑒+𝐴₁𝑒()+Ψ𝑡−𝜏𝐼𝐼𝐼𝑑,with𝑑=^𝑇Θ^𝑇Φ^𝑇𝑇, (34) where 𝐼 represents the identity matrix of appropriate dimensions. Although the asymptotic convergence of error 𝑒 to zero under parametric uncertainties can be ensured by Theorem 10, the performance of the synchronization control can be improved for robustness with the help of additional constraints addressing the minimization of the effects of uncertainties in 𝑑 at error 𝑒 (see also [33–35]). To that end, we provide a sufficient condition for robust asymptotic synchronization of FHN neurons with robustness bound 𝛾 in terms of the 𝐿₂ gain from the uncertain nonlinearities to the error.

Theorem 12. Consider FHN neurons (4) satisfying Assumption 9 with states bounded by a region 𝑥₁,𝑥₂∈[𝑥_min,𝑥_max]. Suppose the optimization problem min𝛾, (35) such that 𝑃_𝑠>0,𝑄_𝑠>0,Γ_𝑠𝑃<0,_𝑟>0,𝑄_𝑟Γ>0,𝛾>0,_𝑟=𝐴^𝑇𝑃_𝑟+𝑃_𝑟𝐴+𝑄_𝑟𝑃_𝑟𝐴₁𝑃_𝑟𝑃_𝑟𝑃_𝑟𝐼∗−𝑄_𝑟0000∗∗−𝛾𝐼000∗∗∗−𝛾𝐼00∗∗∗∗−𝛾𝐼0∗∗∗∗∗−𝛾𝐼<0. (36)(37) Then, the control law (5), with proper selection of parameters 𝐶₀, 𝐶₁, and 𝐶₂ ensures the following.(i)asymptotic synchronization of neurons with zero steady state synchronization error;(ii)the 𝐿₂ gain from the nonlinear uncertainties 𝑑 to the error 𝑒 less than 𝛾.

Proof. Consider the objective function min𝛾 (38) such that 𝐸̇()𝐽()=𝐸̇𝑒,𝑡<0,𝑒,𝑡_𝑟(𝑒,𝑡)+1/𝛾𝑒^𝑇𝑒−𝛾𝑑^𝑇𝑑<0, (39) where 𝐸_𝑟(𝑒,𝑡)=𝑒^𝑇𝑃_𝑟𝑒+^𝑡_𝑡−𝜏𝑒^𝑇(𝜁)𝑄_𝑟𝑒(𝜁)𝑑𝜁>0,with𝑃_𝑟>0,𝑄_𝑟>0, (40) and 𝐸(𝑒,𝑡)>0, as already defined in (29). It has already been shown, in Theorem 10, that 𝐸̇(𝑒,𝑡)<0 leads to matrix inequalities in (36), ensuring the robust asymptotic synchronization of neurons. Hence 𝑒(𝑡)→0 as 𝑡→∞, which completes the proof of statement (i) in Theorem 12. Now, integrating 𝐽 from 𝑡=0 to 𝑡→∞, and multiplying by 𝛾, we obtain ‖𝑒‖²₂<𝛾²‖𝑑‖²₂+𝛾𝐸_𝑟(𝑒,0)−𝐸_𝑟(𝑒,∞). (41) From (40), we have 𝐸_𝑟(𝑒,0)>0, because 𝑃_𝑟>0, 𝑄_𝑟>0, and 𝐸_𝑟(𝑒,∞)=0 (given that 𝑒(𝑡)→0 as 𝑡→∞). Accordingly, (41) shows that the 𝐿₂ gain from 𝑑 to 𝑒 is less than 𝛾. Taking the derivative of (40) along (34) and incorporating it into (39), we obtain 𝐽(𝑒,𝑡)=𝑒^𝑇𝐴^𝑇𝑃_𝑟+𝑃_𝑟𝐴+𝑄_𝑟𝑒+𝑒^𝑇(𝑡−𝜏)𝐴^𝑇₁𝑃_𝑟𝑒+𝑒^𝑇𝑃_𝑟𝐴₁𝑒(𝑡−𝜏)−𝑒^𝑇(𝑡−𝜏)𝑄_𝑟𝑒(𝑡−𝜏)+Ψ^𝑇𝑃_𝑟𝑒+𝑒^𝑇𝑃_𝑟Ψ+Θ^𝑇𝑃_𝑟𝑒+𝑒^𝑇𝑃_𝑟Θ+Φ^𝑇𝑃_𝑟𝑒,+𝑒^𝑇𝑃_𝑟Φ+1/𝛾𝑒^𝑇𝑒−𝛾Ψ^𝑇Ψ+Θ^𝑇Θ+Φ^𝑇Φ (42) which further can be written as 𝐽(𝑒,𝑡)=𝜁^𝑇Π𝜁<0, (43) where Π=𝐴^𝑇𝑃_𝑟+𝑃_𝑟𝐴+𝑄_𝑟+1/𝛾𝐼𝑃_𝑟𝐴₁𝑃_𝑟𝑃_𝑟𝑃_𝑟∗−𝑄_𝑟000∗∗−𝛾𝐼00∗∗∗−𝛾𝐼0∗∗∗∗−𝛾𝐼<0. (44) Using the Schur complement, inequality (44) can be written as Γ_𝑟<0, which completes the proof of statement (ii) in Theorem 12.

Remark 13. The sufficient condition for synchronization of FHN neurons, provided by Theorem 12, ensures ‖𝑒‖²₂<𝛾²‖𝑑‖²₂+𝛾(𝐸_𝑟(𝑒,0)) as 𝐸_𝑟(𝑒,∞)=0. Note that 𝐸_𝑟(𝑒,0) is dependent on initial condition 𝑒(0). Therefore, the effects of the uncertain nonlinear terms contained in 𝑑 and the initial-condition-dependent term 𝐸_𝑟(𝑒,0) are minimized by minimizing 𝛾. Hence, the proposed condition, by managing a single parameter 𝛾, ensures robustness against both the uncertainties and the initial condition.

Remark 14. Recently, global synchronization of FHN models with unknown parameters and in the absence of time-delays, by application of a nonlinear robust adaptive control methodology, was presented [23]. Such a control scheme is computationally complex due to the utilization of adaptation laws. In the present work, synchronization of delayed FHN models under parametric uncertainties is addressed by designing a noncomplex locally robust controller.

In neuronal synchronization, only measurements of the activation potentials and control input 𝑢₁ are available, though the robust synchronization control addressed by Theorem 10 is general for selection of a control law with a single control input 𝑢₁ (by taking 𝐶₂=0) or multiple control inputs 𝑢₁ and 𝑢₂. However, Theorem 12, providing robust synchronization along with robust performance, is better for the two-control-input case in which the second control input requires measurement of recovery potentials. For a single control input, it is better to ensure the minimization of the 𝐿₂ gain from 𝑑 to 𝑒₁ than from 𝑑 to 𝑒, for three reasons. First, the most relevant uncertain neuronal state is the activation potential due to uncertainty in 𝑟 and 𝑔. Because 𝑒₁=𝑥₁−𝑥₂, the above-noted criterion is helpful when dealing with variations in activation potentials. The second reason is that we have no control input to handle the uncertainties Δ𝑏 and Δ𝑣 (because, in the case of neurons, 𝑢₂=0 with a single-input controller). And third, we can still ensure the robustness of 𝑒₂ by minimizing the 𝐿₂ gain from 𝑑 to 𝑒₁, because the dynamical equation 𝑒₂̇=(𝑏+Δ𝑏)𝑒₁−(𝑣+Δ𝑣)𝑒₂ also contains 𝑒₁. This indicates that the robustness of 𝑒₁ will, somehow, ensure the robustness of 𝑒₂. Although obtaining matrix inequalities for this performance criterion is a straightforward extension of Theorem 12, this important issue remains unaddressed in the literature; we prefer therefore to apply the results for this case to the following Theorem.

Theorem 15. Given FHN neurons (4) satisfying Assumption 9 with states bounded by a region 𝑥₁,𝑥₂∈[𝑥_min,𝑥_max], suppose the optimization problem min𝛾, (45) such that 𝑃_𝑠>0,𝑄_𝑠>0,Γ_𝑠𝑃<0,_𝑟>0,𝑄_𝑟>0,𝛾₁𝐴>0,_𝑇𝑃_𝑟+𝑃_𝑟𝐴+𝑄_𝑟𝑃_𝑟𝐴₁𝑃_𝑟𝑃_𝑟𝑃_𝑟diag(1,0)∗−𝑄_𝑟0000∗∗−𝛾₁𝐼000∗∗∗−𝛾₁𝐼00∗∗∗∗−𝛾₁𝐼0∗∗∗∗∗−𝛾₁𝐼<0. (46) Then, the control law (5), with proper selection of parameters 𝐶₀ and 𝐶₁ along with 𝐶₂=0, ensures(i)asymptotic synchronization of neurons with zero steady state synchronization error;(ii)the 𝐿₂ gain from the nonlinear uncertainties 𝑑 to the error 𝑒₁ less than 𝛾₁.

Proof. The proof is similar to that of Theorem 12.

Remark 16. The proposed techniques provide both “robust stability” (in terms of asymptotic convergence of synchronization errors to zero) and “robust performance” (in terms of uncertainty rejection) as addressed by Theorems 10–15. In contrast to the traditional schemes, the robust performance of the proposed synchronization schemes, addressed in Theorems 12–15, is an extra feature of the proposed controller (in addition to asymptotic stabilization).

Remark 17. The present work, in contrast to traditional synchronization techniques for FHN neurons, is novel in many respects. Our techniques consider time-delay between interlinked neurons. The easy LMI-based means of control parameter selection is another key characteristic, as is the exceptional idea of a less conservative local robust nonlinear controller for local synchronization of uncertain delayed FHN neurons. Our work furthers the design of both single- and multiple-input controllers by addressing their robust performances.

Remark 18. The synchronization techniques proposed by Theorems 3–15 for FHN models can be used for a number of purposes. These techniques can be applied to estimate the control signal 𝑢₁ responsible for synchronization of two uncertain delayed neurons, which can be helpful in future for measuring brain activity and for improving stimulation-therapy-based treatments for brain disorders. The proposed schemes for synchronization of two delayed neurons can be generalized to deal with a delayed neural network. Further, such methodologies can be used for biomimetic systems in order to develop artificial neural networks, which can be useful for humanoid robotic applications. Furthermore, synchronization studies, owing to the capability of the FHN model to represent complex processes such as the reaction-diffusion system, can be applied to control (or synchronize) (the chemical kinematics of) industrial plants under time-delays.

Theorem 12 facilitates the selection of a suitable multiple-input robust controller with guaranteed robust stability, ensured by constraints (36), as well as robust performance in terms of 𝐿₂ gain 𝛾, ensured by constraints (37). Theorem 15 provides further results for the special and important case of a single-input robust controller. The traditional techniques, with the help of additional constraints, have not yet adequately addressed robust stability and robust performance simultaneously for synchronization of FHN neurons. In future, the ideas presented in this work will be extended to address the issue of robust synchronization of different uncertain delayed neurons.

5. Simulation Results

To confirm the validity of the proposed schemes, we choose multiple- and single-input controllers 𝐾_I and 𝐾_II, respectively. 𝐾_I=𝐶₀=10,𝐶₁𝐶=0,₂𝐾=10._II=𝐶₀=17,𝐶₁𝐶=0,₂=0. (47)(48) Note that these controllers are global in the absence of parametric uncertainties and local otherwise, due to the feasibility of both Theorems 3 and 10. Moreover, both controllers 𝐾_I and 𝐾_II have no memory feedback, because 𝐶₁=0. By solving Theorem 12, the 𝐿₂ gain 𝛾=0.141 for 𝐾_I is obtained. Local synchronization of the FHN models is considered for the region [𝑥_min,𝑥_max]=[−0.5,1] with 𝐿_𝑎=2. Phase portraits and synchronization error plots for 𝐾_I are shown in Figure 2. Clearly, the FHN neurons are synchronized with zero steady-state synchronization errors, demonstrating suitable robust performance under the parametric uncertainties.

(a)

(b)

(c)

(d)

Figure 2 

Synchronization of delayed coupled FHN neurons with parametric uncertainties using robust multiple-input controller 𝐾I. (a) Phase portrait of 𝑥1 and 𝑦1; (b) phase portrait of 𝑥2 and 𝑦2; (c) synchronization error 𝑒1=𝑥1−𝑥2 versus 𝑡; (d) synchronization error 𝑒2=𝑦1−𝑦2 versus 𝑡.

The controller 𝐾_I, though robust, requires two control inputs in addition to measurement of the recovery potentials. We then check the performance of controller 𝐾_II for the same region [𝑥_min,𝑥_max]=[−0.5,1] and 𝐿_𝑎=2. The 𝐿₂ gain 𝛾₁=0.508 is obtained by solving Theorem 15. Figure 3 provides phase portraits and synchronization error plots validating the results obtained by Theorem 15. The performance of controller 𝐾_II is not better than that of 𝐾_I, though we use a higher value of gain 𝐶₀ for 𝐾_II. This is due to the lack of a second control input 𝑢₂ in the single-input control case. Nonetheless, the results obtained using 𝐾_II, which offers additional computational simplicity with no requirement for multiple control inputs or measurement of recovery potentials, are reasonable.

(a)

(b)

(c)

(d)

Figure 3 

Synchronization of delayed coupled FHN neurons with parametric uncertainties using robust single-input controller 𝐾II. (a) Phase portrait of 𝑥1 and 𝑦1; (b) phase portrait of 𝑥2 and 𝑦2; (c) synchronization error 𝑒1=𝑥1−𝑥2 versus 𝑡; (d) synchronization error 𝑒2=𝑦1−𝑦2 versus 𝑡.

6. Conclusions

This paper provided a brief look at coupled delayed FHN neurons with various parametric variations under EES. A global nonlinear control law was developed for the asymptotic synchronization of delayed FHN neurons. By integrating the ideas of neuronal state bounds, local and global Lipschitz conditions for the nonlinear and uncertain components of the dynamics of delayed neurons, the LK functional, and 𝐿₂ gain reduction, a less conservative regional robust synchronization control was developed that ensures both robust stability and robust performance. Computational simplicity, a simple design procedure, guaranteed zero steady-state synchronization error, a computed robustness bound and applicability to both multiple- and single-input controllers are additional distinguishing features of the proposed schemes. Simulations of the uncertain coupled chaotic delayed FHN neuronal synchronization validated the proposed methodology.

Acknowledgment

This research was supported by the World Class University Program funded by the Ministry of Education, Science, and Technology through the National Research Foundation of Korea (Grant no. R31-20004).