Computational and Mathematical Methods in Medicine

Volume 2012 (2012), Article ID 230980, 11 pages

doi:10.1155/2012/230980

Research Article

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

Muhammad Rehan¹ and Keum-Shik Hong²

¹Department of Electrical Engineering, Pakistan Institute of Engineering and Applied Sciences (PIEAS), P.O. Box 45650, Islamabad, Pakistan

²Department of Cogno-Mechatronics Engineering and School of Mechanical Engineering, Pusan National University, 30 Jangjeon-Dong, Geumjeong-Gu, Busan 609-735, Republic of Korea

Received 11 August 2012; Revised 3 October 2012; Accepted 3 October 2012

Academic Editor: Ricardo Femat

Copyright © 2012 Muhammad Rehan and Keum-Shik Hong. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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 , where denotes the Euclidian norm of . A positive definite symmetric matrix is denoted as . For with the th diagonal entry and , represents a diagonal matrix.

2. Model Description

Consider two identical uncertain coupled chaotic FHN neurons with time-delays under EES:

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 indicates the time-delay due to separation between neurons. The master and the slave neurons are under external electrical stimulation with currents and , respectively. Here, and 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

along with the parametric uncertainties

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

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 and ; (b) phase portrait of and ; (c) synchronization error versus ; (d) synchronization error 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

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 and 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 makes the controller input , 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

are satisfied, where

Proof. Using the control law (5) in (4) under Assumption 2, the overall closed-loop system becomes

Defining the synchronization errors as

the error dynamics are

The error dynamics can be written as

where , , and matrices and are given by (7). Now consider the Lyapunov functional

such that , . The derivative of (12) along (11) is

For stability, , that is, , 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

where . Moreover, the Lipschitz nonlinearities also satisfy

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 , where and are the minimum and maximum limits on (or ), respectively.

From Definition 8, it can be seen that

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 . The nonlinear control law (5) synchronizes the coupled neurons asymptotically with proper selection of parameters , , and , if the following matrix inequalities are verified:

where

and matrices and are given by (7).

Proof. Incorporating the control law (5) into (4), the overall closed-loop system becomes

Using the same procedure as discussed in the previous section, the following error dynamics model is obtained:

where

On the basis of Assumption 9 and inequality (16), we have

Combining the local and global Lipschitz constraints of (22), we have

As we know from the global Lipschitz condition,

Inequalities (23)–(25) imply

and (26) implies

Constructing the LK functional

with and .

Taking the derivative of (29) along (20)

and, using inequalities (27)-(28), we obtain

This further implies that

where

For stability, , and hence , 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 , , , and , 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

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 . Suppose the optimization problem

such that

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

such that

where

and , as already defined in (29). It has already been shown, in Theorem 10, that leads to matrix inequalities in (36), ensuring the robust asymptotic synchronization of neurons. Hence as , which completes the proof of statement (i) in Theorem 12. Now, integrating from to , and multiplying by , we obtain

From (40), we have , because , , and (given that 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

which further can be written as

where

Using the Schur complement, inequality (44) can be written as , 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 as . Note that is dependent on initial condition . Therefore, the effects of the uncertain nonlinear terms contained in and the initial-condition-dependent term 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 ) 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, 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 , suppose the optimization problem

such that

Then, the control law (5), with proper selection of parameters and along with , 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 and , respectively.

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 and have no memory feedback, because . By solving Theorem 12, the gain for is obtained. Local synchronization of the FHN models is considered for the region with . Phase portraits and synchronization error plots for 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 . (a) Phase portrait of and ; (b) phase portrait of and ; (c) synchronization error versus ; (d) synchronization error versus .

The controller , though robust, requires two control inputs in addition to measurement of the recovery potentials. We then check the performance of controller for the same region and . The gain 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 is not better than that of , though we use a higher value of gain for . This is due to the lack of a second control input in the single-input control case. Nonetheless, the results obtained using , 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 . (a) Phase portrait of and ; (b) phase portrait of and ; (c) synchronization error versus ; (d) synchronization error 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).

References

A. D. Garbo, M. Barbi, and S. Chillemi, “The synchronization properties of a network of inhibitory interneurons depend on the biophysical model,” BioSystems, vol. 88, no. 3, pp. 216–227, 2007.

M. Aqil, K.-S. Hong, M.-Y. Jeong, and S. S. Ge, “Detection of event-related hemodynamic response to neuroactivation by dynamic modeling of brain activity,” Neuroimage, vol. 63, no. 1, pp. 553–568, 2012.

A. Alvarellos-González, A. Pazos, and A. B. Porto-Pazos, “Computational models of neuron-astrocyte interactions lead to improved efficacy in the performance of neural networks,” Computational and Mathematical Methods in Medicine, vol. 2012, Article ID 476324, 10 pages, 2012.

X.-S. Hu, K.-S. Hong, and S. S. Ge, “fNIRS-based online deception decoding,” Journal of Neural Engineering, vol. 9, no. 2, pp. 1–8, 2012.

E. Kuriscak, P. Marsalek, J. Stroffek, and Z. Wünsch, “The effect of neural noise on spike time precision in a detailed CA3 neuron model,” Computational and Mathematical Methods in Medicine, vol. 2012, Article ID 595398, 16 pages, 2012.

L. H. Nguyen and K.-S. Hong, “Adaptive synchronization of two coupled chaotic Hindmarsh—Rose neurons by controlling the membrane potential of a slave neuron,” Applied Mathematical Modelling. In press

R. Femat and G. Solís-Perales, “On the chaos synchronization phenomena,” Physics Letters A, vol. 262, no. 1, pp. 50–60, 1999.

O. Cornejo-Pérez and R. Femat, “Unidirectional synchronization of Hodgkin-Huxley neurons,” Chaos, Solitons and Fractals, vol. 25, no. 1, pp. 43–53, 2005.

P. Limousin and I. M. Torres, “Deep brain stimulation for Parkinson's disease,” Neurotherapeutics, vol. 5, no. 2, pp. 309–319, 2008.

B. Jobst, “Brain stimulation for surgical epilepsy,” Epilepsy Research, vol. 89, no. 1, pp. 154–161, 2010.

J. L. Ostrem and P. A. Starr, “Treatment of dystonia with deep brain stimulation,” Neurotherapeutics, vol. 5, no. 2, pp. 320–330, 2008.

J. Wang, T. Zhang, and B. Deng, “Synchronization of FitzHugh-Nagumo neurons in external electrical stimulation via nonlinear control,” Chaos, Solitons and Fractals, vol. 31, no. 1, pp. 30–38, 2007.

T. Zhang, J. Wang, X. Fei, and B. Deng, “Synchronization of coupled FitzHugh-Nagumo systems via MIMO feedback linearization control,” Chaos, Solitons and Fractals, vol. 33, no. 1, pp. 194–202, 2007.

M. Aqil, K.-S. Hong, and M. Y. Jeong, “Synchronization of coupled chaotic FitzHugh-Nagumo systems,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 4, pp. 1615–1627, 2012.

C. J. Thompson, D. C. Bardos, Y. S. Yang, and K. H. Joyner, “Nonlinear cable models for cells exposed to electric fields I. General theory and space-clamped solutions,” Chaos, solitons and fractals, vol. 10, no. 11, pp. 1825–1842, 1999.

L. H. Nguyen and K.-S. Hong, “Synchronization of coupled chaotic FitzHugh-Nagumo neurons via Lyapunov functions,” Mathematics and Computers in Simulation, vol. 82, no. 4, pp. 590–603, 2011.

J. Wang, B. Deng, and X. Fei, “Synchronization of chaotic neurons coupled with gap junction with time delays in external electrical stimulation,” Chaos, Solitons and Fractals, vol. 35, no. 3, pp. 512–518, 2008.

N. Burić, K. Todorović, and N. Vasović, “Dynamics of noisy FitzHugh-Nagumo neurons with delayed coupling,” Chaos, Solitons and Fractals, vol. 40, no. 5, pp. 2405–2413, 2009.

B. Wang, “Pullback attractors for the non-autonomous FitzHugh-Nagumo system on unbounded domains,” Nonlinear Analysis, Theory, Methods and Applications, vol. 70, no. 11, pp. 3799–3815, 2009.

D. Q. Wei, X. S. Luo, B. Zhang, and Y. H. Qin, “Controlling chaos in space-clamped FitzHugh-Nagumo neuron by adaptive passive method,” Nonlinear Analysis, vol. 11, no. 3, pp. 1752–1759, 2010.

J. Wang, Z. Zhang, and H. Li, “Synchronization of FitzHugh-Nagumo systems in EES via H∞ variable universe adaptive fuzzy control,” Chaos, Solitons and Fractals, vol. 36, no. 5, pp. 1332–1339, 2008.

Y. Q. Che, J. Wang, S. S. Zhou, and B. Deng, “Robust synchronization control of coupled chaotic neurons under external electrical stimulation,” Chaos, Solitons and Fractals, vol. 40, no. 3, pp. 1333–1342, 2009.

M. Rehan and K.-S. Hong, “LMI-based robust adaptive synchronization of FitzHugh-Nagumo neurons with unknown parameters under uncertain external electrical stimulation,” Physics Letters A, vol. 375, no. 15, pp. 1666–1670, 2011.

M. Rehan, K.-S. Hong, and M. Aqil, “Synchronization of multiple chaotic FitzHugh-Nagumo neurons with gap junctions under external electrical stimulation,” Neurocomputing, vol. 74, no. 17, pp. 3296–3304, 2011.

R. Femat, L. Kocarev, L. Van Gerven, and M. E. Monsivais-Pérez, “Towards generalized synchronization of strictly different chaotic systems,” Physics Letters A, vol. 342, no. 3, pp. 247–255, 2005.

J. G. B. Ramírez, K. P. C. Galarza, and R. Femat, “Generalized synchronization of strictly different systems: partial-state synchrony,” Chaos, Solitons & Fractals, vol. 45, no. 3, pp. 193–204, 2012.

C. Wang and S. S. Ge, “Synchronization of two uncertain chaotic systems via adaptive backstepping,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 11, no. 6, pp. 1743–1751, 2001.

C. Wang and S. S. Ge, “Adaptive synchronization of uncertain chaotic systems via backstepping design,” Chaos, solitons and fractals, vol. 12, no. 7, pp. 1199–1206, 2001.

F. Blanchini, “The gain scheduling and the robust state feedback stabilization problems,” IEEE Transactions on Automatic Control, vol. 45, no. 11, pp. 2061–2070, 2000.

A. Iftar and U. Ozguner, “Local LQG/LTR controller design for decentralized systems,” IEEE Transactions on Automatic Control, vol. 32, no. 10, pp. 926–930, 1987.

S. Skogestad and I. Postlethwaite, Multivariable Feedback Control Analysis and Design, John Wiley & Sons, England, UK, 2nd edition, 2005.

G. Lu and D. W. C. Ho, “Full-order and reduced-order observers for Lipschitz descriptor systems: the unified LMI approach,” IEEE Transactions on Circuits and Systems II, vol. 53, no. 7, pp. 563–567, 2006.

E. Ruiz-Velázquez, R. Femat, and D. U. Campos-Delgado, “Blood glucose control for type I diabetes mellitus: a robust tracking H∞ problem,” Control Engineering Practice, vol. 12, no. 9, pp. 1179–1195, 2004.

C. P. Flores-Gutiérrez, R. Femat, and E. Ruiz-Velázquez, “On hypoglycemic levels induced by H∞ control on type i diabetes mellitus,” Applied Mathematics and Computation, vol. 218, no. 2, pp. 376–385, 2011.

M. Rehan, A. Ahmed, and N. Iqbal, “Static and low order anti-windup synthesis for cascade control systems with actuator saturation: an application to temperature-based process control,” ISA Transactions, vol. 49, no. 3, pp. 293–301, 2010.