Linear matrix inequality approach for synchronization control of fuzzy cellular neural networks with mixed time delays

Chin. Phys. B Vol. 21, No. 4 (212 4842 Linear matrix inequality approach for synchronization control of fuzzy cellular neural networks with mixed time delays P. Balasubramaniam a, M. Kalpana a, and R. Rakkiyappan b a Department of Mathematics, Gandhigram Rural Institute Deemed University, Gandhigram-624 32, Tamilnadu, India b Department of Mathematics, Bharathiar University, Coimbatore-641 46, Tamilnadu, India (Received 26 August 211; revised manuscript received 3 August 211 Fuzzy cellular neural networks (FCNNs are special kinds of cellular neural networks (CNNs. Each cell in an FCNN contains fuzzy operating abilities. The entire network is governed by cellular computing laws. The design of FCNNs is based on fuzzy local rules. In this paper, a linear matrix inequality (LMI approach for synchronization control of FCNNs with mixed delays is investigated. Mixed delays include discrete time-varying delays and unbounded distributed delays. A dynamic control scheme is proposed to achieve the synchronization between a drive network and a response network. By constructing the Lyapunov Krasovskii functional which contains a triple-integral term and the free-weighting matrices method an improved delay-dependent stability criterion is derived in terms of LMIs. The controller can be easily obtained by solving the derived LMIs. A numerical example and its simulations are presented to illustrate the effectiveness of the proposed method. Keywords: asymptotic stability, chaos, fuzzy cellular neural networks, linear matrix inequalities, synchronization PACS: 84.35.i, 5.45.Gg, 5.45.Xt, 7.5.Mh DOI: 1.188/1674-156/21/4/4842 1. Introduction Cellular neural networks (CNNs are locally connected nonlinear networks. They originally stemmed from cellular automata and artificial neural networks. Local connectedness is the most significant property of CNN. The local connectedness restricts the ability of a CNN to solve many global problems that cannot be decomposed into local components. However, the local property has its advantages, such as easy implementation using VLSI technology and efficiency for solving local problems, as proposed by Chua and Yang in Refs. 1] and 2]. Fuzzy set theory provides mathematical support to the capture of uncertainties associated with human cognitive processes. A fuzzy cellular neural network (FCNN is a generalized case of the CNN structure proposed by Yang in Refs. 3] and 4]. In the past decade, synchronization of chaotic systems has been extensively studied due to its potential applications in many different areas including secure communication, chemical and biological systems, information science, optics, and so on. Carroll and Pecora 5] have introduced the drive-response concept, and used the output of the drive system to control the response system so that state synchronization was achieved. The study of the dynamical properties of FC- NNs has mainly concentrated on stability analysis. 6,7] Time delays in the FCNNs make the dynamic behaviors become more complicated, and may destabilize the stable equilibria and admit periodic oscillation, bifurcation, and chaos. Therefore, considerable attention has been paid to the study of delay systems in control theory. 8] In fact, as special complex networks, delayed neural networks (DNNs have been found to exhibit complicated dynamics and even chaotic behavior if the parameters and time delays are appropriately chosen. 9 18] However, CNNs and FCNNs have also been found to exhibit complicated dynamics and even chaotic behavior if the parameters and time delays are appropriately chosen. 19 24] For instance, Xia et al. 22] investigated the synchronization schemes for coupled identical Yang Yang type FCNNs. To the best of the authors knowledge, there are Project supported by No. DST/INSPIRE Fellowship/21/293]/dt. 18/3/211. Corresponding author. E-mail: balugru@gmail.com 212 Chinese Physical Society and IOP Publishing Ltd http://iopscience.iop.org/cpb http://cpb.iphy.ac.cn 4842-1

Chin. Phys. B Vol. 21, No. 4 (212 4842 no results on the linear matrix inequality (LMI approach for synchronization control of FCNNs with discrete and unbounded distributed delays. Based on the Lyapunov Krasovskii functional, which contains a triple-integral term and free-weighting matrices method, an improved delay-dependent stability criterion is derived in terms of linear matrix inequalities (LMIs. For the synchronization conditions given in Refs. 21] and 24], the time delays are constant delays or time-varying delays that are differentiable such that their derivatives are not greater than one or finite. A new result is presented by relaxing these restrictive conditions on time-varying delays, which means that our presented results have a wider adaptive range. Finally, a numerical example and numerical simulations are given to show the effectiveness of our result. This paper is organized as follows. In Section 2, we introduce some model formulation and preliminaries. In Section 3, the main results for asymptotical synchronization of chaotic FCNNs are presented. In Section 4, a numerical example is given in order to demonstrate the effectiveness of our results. Finally, in Section 5, some conclusions are summarized. Notations R n denotes the n-dimensional Euclidean space. For any matrix A = a ij ] n n, let A T and A 1 denote the transpose and the inverse of A, respectively. A = a ij ] n n. Let A > (A < denote the positive-definite (negative-definite symmetric matrix, respectively. I denotes the identity matrix of appropriate dimension and Λ = {1, 2,..., n}. Symbol denotes the symmetric terms in a symmetric matrix. 2. Model formulation and preliminaries Most of the synchronization methods belong to the master-slave (drive-response type, which means that the two systems are coupled so that the behavior of the second one is influenced by the behavior of the first one, but the behavior of the first one is independent of the second one. The first system will be called the master system or drive system, and the second system will be the slave system or response system. The main objective of this paper is to design a controller to let the slave system synchronize with the master system. Now let us consider the following chaotic FCNN with discrete and unbounded distributed delays ẋ i (t = d i x i (t a ij f j (x j (t b ij f j (x j (t τ(t I i k j (t sf j (x j (sds k j (t sf j (x j (sds, i Λ, x i (s = ϕ i (s, s (, ], (1 where ϕ i ( C((, ], R; and are the elements of fuzzy feedback MIN template, fuzzy feedback MAX template, respectively; a ij and b ij are the elements of feedback template; and denote the fuzzy AND and fuzzy OR operations, respectively; x i and I i denote the state and external input of the i-th neuron, respectively; d i is a diagonal matrix, d i represents the rates with which the i-th neuron will reset their potential to the resting state in isolation when disconnected from the networks and external inputs; f j represents the neuron activation function; k j ( is the delay kernel function. In this section, the following assumptions are needed: A1 The neuron activation functions f j ( are continuously bounded and satisfy l j f j(u f j (v u v l j, for any u, v R, u v, j Λ, where l j and l j are some real constants and they may be positive, zero or negative. A2 The transmission delay τ(t is a time varying delay, and it satisfies τ(t τ, where τ is a positive constant. A3 The delay kernels k j are some real-valued non-negative continuous functions defined in, satisfying k j : = ˆk j : = for some η >. k j (sds > and k j (s e ηs ds <, j Λ, 4842-2

Chin. Phys. B Vol. 21, No. 4 (212 4842 Moreover, considering that a chaotic FCNN system depends extremely on initial values, the initial condition of the response system is defined to be different from that of the drive system. Therefore, a response system with the same form and parameters of drive system (1 is introduced as follows: ẏ i (t = d i y i (t a ij f j (y j (t b ij f j (y j (t τ(t I i k j (t sf j (y j (sds (2 k j (t sf j (y j (sds u i (t, i Λ, y i (s = φ i (s, s (, ], where φ i ( C((, ], R; d i, a ij, b ij,,, I i, and f j are the same as system (1 and u i (t is the appropriate control input that will be designed in order to obtain a certain control objective. Let e(t = ( e 1 (t, e 2 (t,..., e n (t : = x(t y(t be the error state. Then, the error dynamical system between system (1 and system (2 is given by ė i (t = d i e i (t a ij g j (e j (t b ij g j (e j (t τ(t e i (s = ϕ i (s φ i (s, s (, ], k j (t sf j (x j (sds k j (t sf j (x j (sds where g j (e j ( = f j (x j ( f j (y j (, j Λ. The appropriate control input is k j (t sf j (y j (sds k j (t sf j (y j (sds u i (t, i Λ, u i (t : = h ij e i (t, i, j Λ. (4 (3 Substituting Eq. (4 into system (3, it is easy to have ė i (t = ( d i h ij e i (t a ij g j (e j (t b ij g j (e j (t τ(t e i (s = ϕ i (s φ i (s, s (, ]. k j (t sf j (x j (sds k j (t sf j (x j (sds k j (t sf j (y j (sds k j (t sf j (y j (sds, i Λ, (5 Remark 1 Note that the assumptions A1 and A2 on activation function and time-varying delay in this paper are weaker than those generally used in the literature, 21,24] namely, the boundedness of the activation function f j. Further, the differentiability of the time varying delay τ(t is not required in this paper. Now, we state the following a few lemmas which will be used in the following discussion. matrix inequality Q(x S(x >, S T (x R(x where Q(x = Q T (x, R(x = R T (x is equivalent to R(x > and Q(x S(xR 1 (xs T (x >. Lemma 2 For any x, y R n, ϵ >, and positivedefinite matrix Q R n n, the following matrix inequality holds: Lemma 1 (Schur complement 25] The linear 2x T y ϵx T Qx ϵ 1 y T Q 1 y. 4842-3

Chin. Phys. B Vol. 21, No. 4 (212 4842 Lemma 3 26] Let z, z be two states of system (1, then we have f j (z f j (z f j (z f j (z, f j (z f j (z f j (z f j (z. Lemma 4 27] For any x R n, any constant matrix A = a ij ] n n with a ij, the following matrix inequality holds: x T A T Ax nx T A T s A s x, where ( n A s = diag a i1, a i2,..., a in. Lemma 5 28] Given any real matrix M = M T > of appropriate dimension, a scalar η >, and a vector function ω( : a, b] R n, such that the integrations concerned are well defined, then T b ] b ω(sds] M ω(sds a (b a b 3. Main results a a ω T (smω(sds. Theorem 1 Assume that assumptions A1 A3 hold. The error dynamical system (3 is globally asymptotically stable, if there exist n n positive diagonal matrices P, R, U 1, and U 2, some n n positive-definite symmetric matrices Q, W, M 1, and M 2, ( two scalars µ 1 >, µ 2 >, and a 2n 2n matrix T11 T 12 T 22 > such that the following LMI has feasible solution: Ω i,j Γ1 T Γ2 T Ω = µ 1 n 1 I <, (6 µ 2 n 1 I where i, j = 1, 2,..., 1 with Ω 1,1 = 2P D Q T 1 τ M 1 2M 2 U 1 Σ 1, Ω 1,2 = 1 τ M 1, Ω 1,3 = T T 12, Ω 1,4 = D T P T Q T, Ω 1,5 =, Ω 1,6 = 2 τ M 2, Ω 1,7 = 2 τ M 2, Ω 1,8 = P A U 1 Σ 2, Ω 1,9 = P B, Ω 1,1 =, Ω 2,2 = Q 1 τ M 1, Ω 2,3 =, Ω 2,4 =, Ω 2,5 =, Ω 2,6 =, Ω 2,7 =, Ω 2,8 =, Ω 2,9 =, Ω 2,1 =, Ω 3,3 = τt 11 2T T 12 U 2 Σ 1, Ω 3,4 =, Ω 3,5 =, Ω 3,6 =, Ω 3,7 =, Ω 3,8 =, Ω 3,9 = U 2 Σ 2, Ω 3,1 =, Ω 4,4 = 2P W τm 1 τ 2 2 M 2 τt 22, Ω 4,5 =, Ω 4,6 =, Ω 4,7 =, Ω 4,8 = P A, Ω 4,9 = P B, Ω 4,1 =, Ω 5,5 = W, Ω 5,6 =, Ω 5,7 =, Ω 5,8 =, Ω 5,9 =, Ω 5,1 =, Ω 6,6 = 2 τ 2 M 2, Ω 6,7 = 2 τ 2 M 2, Ω 6,8 =, Ω 6,9 =, Ω 6,1 =, Ω 7,7 = 2 τ 2 M 2, Ω 7,8 =, Ω 7,9 =, Ω 7,1 =, Ω 8,8 = RK U 1, Ω 8,9 =, Ω 8,1 =, Ω 9,9 = U 2, Ω 9,1 =, Ω 1,1 = µ 1 I µ 2 I R, K = diag { k1ˆk1,..., k } nˆkn, { n } α s = diag α i1, α i2,..., α in, { n } β s = diag β i1, β i2,..., β in, S = α s β s, Γ T 1 = Γ T 2 = (W S T ] T, (W S T ] T. Moreover, the estimation gain H = P 1 Q. Proof Consider the following Lyapunov Krasovskii functional 7 V (t = V i (t, (7 where V 1 (t = e T (tp e(t = V 2 (t = V 3 (t = p i e 2 i (t, e T (sqe(sds r j kj ė T (sw ė(sds, k j (θ gj 2 (e j (sdsdθ, t θ 4842-4

Chin. Phys. B Vol. 21, No. 4 (212 4842 V 4 (t = V 5 (t = tθ θ ė T (sm 1 ė(sdsdθ, tλ ė T (sm 2 ė(sdsdλdθ, θ e(θ τ(θ V 6 (t = θ(θ ė(s e(θ τ(θ dsdθ, ė(s V 7 (t = tθ T ė T (st 22 ė(sdsdθ. T 11 T 12 T 22 From Lemma 3, we obtain = k j (t sf j (x j ds k j (t sf j (y j ds ( k j (t s f j (x j f j (y j ds k j (t sg j (e j (sds. By calculating the time derivation of V i (t along the trajectory of system (5, we obtain V 1 (t = 2 { p i e i (t ( d i h ij e i (t 2e T (tp a ij g j (e j (t k j (t sf j (x j ds k j (t sf j (x j ds b ij g j (e j (t τ(t ] ( D He(t Ag(e(t Bg(e(t τ(t k j (t sf j (y j ds } k j (t sf j (y j ds 2e T (tp ( α β 2e T (tp De(t 2e T (tqe(t 2e T (tp Ag(e(t 2e T (tp Bg(e(t τ(t µ 1 1 net (tp ( α s β s ( α s β s T P e(t ( T ( µ 1 k(t sg(e(sds k(t sg(e(sds ] 2ė T (tp ė(t ė(t 2e T (tp De(t 2e T (tqe(t 2e T (tp Ag(e(t 2e T (tp Bg(e(t τ(t µ 1 1 net (tp ( α s β s ( α s β s T P e(t ( T ( µ 1 k(t sg(e(sds k(t sg(e(sds 2ė T (tp ė(t 2ė T (tp De(t 2ė T (tqe(t 2ė T (tp Ag(e(t 2ė T (tp Bg(e(t τ(t k(t sg(e(sds µ 1 2 nėt (tp ( α s β s ( α s β s T P ė(t ( T ( µ 2 k(t sg(e(sds k(t sg(e(sds, (8 V 2 (t = e T (tqe(t e T (t τqe(t τ ė T (tw ė(t ė T (t τw ė(t τ, (9 V 3 (t = r j kj k j (θgj 2 (e j (tdθ r j kj k j (θgj 2 (e j (t θdθ g T (e(trkg(e(t g T (e(trkg(e(t r j k j (θdθ ( r j k j (θg 2 j (e j (t θdθ 2 k j (θg j (e j (t θdθ 4842-5

( = g T (e(trkg(e(t Chin. Phys. B Vol. 21, No. 4 (212 4842 V 4 (t = τė T (tm 1 ė(t ė T (sm 1 ė(sds V 5 (t = τ 2 V 6 (t = TR ( k(t sg(e(sds k(t sg(e(sds, (1 τė T (tm 1 ė(t 1 τ et (tm 1 e(t 2 τ et (tm 1 e(t τ 1 τ et (t τm 1 e(t τ, (11 2 ėt (tm 2 ė(t τ 2 2 ėt (tm 2 ė(t 2 τ 2 ( tθ = τ 2 2 ėt (tm 2 ė(t 2 τ 2 ( τe(t = τ 2 2 ėt (tm 2 ė(t 2 τ 2 ( τe(t ( τe(t (t e(sds ė T (sm 2 ė(sdsdθ tθ (t (t TM2 ( ė(sdsdθ TM2 ( e(sds τe(t e(sds e(sds = τ 2 2 ėt (tm 2 ė(t 2e T (tm 2 e(t 2 τ et (tm 2 2 (t τ et (tm 2 e(sds 2 τ 2 τ 2 (t (t (t e(sds 2 τ 2 (t e T (sdsm 2 e(t τ(t ė(s T (t (t (t tθ ė(sdsdθ e(sds e(sds TM2 e(sds e T (sdsm 2 e(t 2 τ 2 (t e T (sdsm 2 e(sds 2 τ (t T 11 T 12 T 22 e(sds 2 (t τ 2 e(t τ(t ds ė(s (t (t e T (sdsm 2 e T (sdsm 2 e(t (t e T (sdsm 2 e(sds, (12 = τ(te T (t τ(tt 11 e(t τ(t 2e T (tt T 12e(t τ(t 2e T (t τ(tt T 12e(t τ(t (t ė T (st 22 ė(sds ] e T (t τ(t τt 11 2T12 T e(t τ(t 2e T (tt12e(t T τ(t V 7 (t = τė T (tt 22 ė(t ė T (t θt 22 ė(t θdθ = τė T (tt 22 ė(t ė T (st 22 ė(sds, (13 ė T (st 22 ė(sds. (14 In addition, for any n n diagonal matrices U 1 >, U 2 >, the following inequality holds by the methods proposed in Ref. 29]: e(t g(e(t T U 1Σ 1 U 1 Σ 2 U 1 e(t e(t τ(t g(e(t g(e(t τ(t T U 2Σ 1 U 2 Σ 2 U 2 e(t τ(t. (15 g(e(t τ(t Hence, from expressions (7 (15 we have where V (t ξ T (t Ω i,j Γ T 1 µ 1 1 nγ 1 Γ T 2 µ 1 2 nγ 2] ξ(t, = ξ T (tωξ(t, (16 4842-6

ξ(t = Chin. Phys. B Vol. 21, No. 4 (212 4842 e T (t, e T (t τ, e T (t τ(t, ė T (t, ė T (t τ, T g T (e(t, g T (e(t τ(t, k(t sg (e(sds] T, Ω = Ω i,j Γ T 1 µ 1 1 nγ 1 Γ T 2 µ 1 2 nγ 2. (t e T (sds, (t e T (sds, By expression (6, it yields V (t ξ T (t Ω ξ(t, t >, where Ω = Ω >, which guarantees globally asymptotically stable of the error dynamical system (3 by Lyapunov stability theory. As a result, the chaotic FCNN with discrete and distributed delays (2 is globally synchronized with the master FCNN (1. Remark 2 Yu and Jiang 23] have investigated the global exponential synchronization of fuzzy cellular neural networks with delays and reaction diffusion terms. Feng et al. 24] have investigated the global exponential synchronization of delayed fuzzy cellular neural networks with impulsive effects. The authors in the previous literature have tried to solve the synchronization of FCNNs with constant delays. So far no results have been reported on the synchronization of FCNNs with unbounded distributed delays using the LMI approach. It is challenging and complicated to synchronize the FCNNs model with unbounded distributed delays. By constructing a suitable Lyapunov krasovskii functional and simple control, we derive a simple and efficient criterion in terms of LMIs for synchronization. The results obtained in this paper are simple and easy to apply in practical applications. Let us consider the master FCNN described by the following state equation: ẋ i (t = d i x i (t a ij f j (x j (t b ij f j (x j (t τ(t I i f j (x j (t τ(t f j (x j (t τ(t, i Λ, x i (s = ϕ i (s, s (, ]. The corresponding response system of system (17 is given by ẏ i (t = d i y i (t a ij f j (y j (t b ij f j (y j (t τ(t I i f j (y j (t τ(t f j (y j (t τ(t u i (t, i Λ, y i (s = φ i (s, s (, ], (17 (18 where ϕ i (, φ i (, d i, a ij, b ij,,, I i, f j, and u i (t described in systems (17 and (18 are the same as those in systems (1 and (2. Let e(t = ( e 1 (t, e 2 (t,..., e n (t : = x(t y(t be the error state. Then, the error dynamical system between systems (17 and (18 is given by ė i (t = d i e i (t a ij g j (e j (t b ij g j (e j (t τ(t f j (x j (t τ(t f j (x j (t τ(t e i (s = ϕ i (s φ i (s, s (, ], f j (y j (t τ(t f j (y j (t τ(t u i (t, i Λ, (19 4842-7

Chin. Phys. B Vol. 21, No. 4 (212 4842 where g j (e j ( = f j (x j ( f j (y j (, j Λ. Theorem 2 Assume that assumptions A1 and A2 hold. The error dynamical system (19 is globally asymptotically stable, if there exist n n positive diagonal matrices P, U 1, U 2, and some n n positive-definite symmetric matrices Q, W, M 1, M 2, two scalars µ 1 >, µ 2 >, and a 2n 2n matrix ( T11 T 12 T 22 > such that the following LMI has feasible solution: Ω i,j Γ1 T Γ2 T Ω = µ 1 n 1 I <, (2 µ 2 n 1 I where i, j = 1, 2,..., 9 with Ω 1,1 = 2P D Q T 1 τ M 1 2M 2 U 1 Σ 1, Ω 1,2 = 1 τ M 1, Ω 1,3 = T T 12, Ω 1,4 = D T P T Q T, Ω 1,5 =, Ω 1,6 = 2 τ M 2, Ω 1,7 = 2 τ M 2, Ω 1,8 = P A U 1 Σ 2, Ω 1,9 = P B, Ω 1,1 =, Ω 2,2 = Q 1 τ M 1, Ω 2,3 =, Ω 2,4 =, Ω 2,5 =, Ω 2,6 =, Ω 2,7 =, Ω 2,8 =, Ω 2,9 =, Ω 2,1 =, Ω 3,3 = τt 11 2T T 12 U 2 Σ 1, Ω 3,4 =, Ω 3,5 =, Ω 3,6 =, Ω 3,7 =, Ω 3,8 =, Ω 3,9 = U 2 Σ 2, Ω 3,1 =, Ω 4,4 = 2P W τm 1 τ 2 2 M 2 τt 22, Ω 4,5 =, Ω 4,6 =, Ω 4,7 =, Ω 4,8 = P A, Ω 4,9 = P B, Ω 4,1 =, Ω 5,5 = W, Ω 5,6 =, Ω 5,7 =, Ω 5,8 =, Ω 5,9 =, Ω 5,1 =, Ω 6,6 = 2 τ 2 M 2, Ω 6,7 = 2 τ 2 M 2, Ω 6,8 =, Ω 6,9 =, Ω 6,1 =, Ω 7,7 = 2 τ 2 M 2, Ω 7,8 =, Ω 7,9 =, Ω 7,1 =, Ω 8,8 = U 1, Ω 8,9 =, Ω 8,1 =, Ω 9,9 = µ 1 I µ 2 I U 2, Ω 9,1 =, { n } α s = diag α i1, α i2,..., α in, { n } β s = diag β i1, β i2,..., β in, ] T S = α s β s, Γ1 T = (W S T, ] T Γ2 T = (W S T. Moreover, the estimation gain H = P 1 Q. Proof Consider the following Lyapunov Krasovskii functional where V (t = V 1 (t = e T (tp e(t = V 2 (t = V 3 (t = V 4 (t = 6 V i (t, (21 p i e 2 i (t, e T (sqe(sds tθ θ ė T (sm 1 ė(sdsdθ, tλ ė T (sw ė(sds, ė T (sm 2 ė(sdsdλdθ, θ e(θ τ(θ V 5 (t = θ(θ ė(s e(θ τ(θ dsdθ, ė(s V 6 (t = tθ T ė T (st 22 ė(sdsdθ. T 11 T 12 T 22 The proof of this Theorem is immediately follows from Theorem 1. 4. Numerical example Example 1 Consider the chaotic FCNN model ẋ i (t = d i x i (t a ij f j (x j (t b ij f j (x j (t τ(t I i k j (t sf j (x j (sds k j (t sf j (x j (sds, i Λ, x i (s = ϕ i (s, s (, ] (22 4842-8

with parameters defined as Chin. Phys. B Vol. 21, No. 4 (212 4842 ϕ(s = (1,.5, 1 T, s (, ], f j (x j = 1 2( xj 1 x j 1, j = 1, 2, 3, k(s =.17 e s, I i =, i = 1, 2, 3, τ(t =.4 sin(t, 1.25 3.21 3.2 4.3 7.5 3 1 A = 3.2 1.1 4.4 B = 3 1.2 5 D = 1 3.2 4.4 1 3.2 4.5 2.3 1 1/32 1/32 1/32 α = β = 1/32 1/32 1/32. 1/32 1/32 1/32 The corresponding response system is designed as ẏ i (t = d i y i (t a ij f j (y j (t b ij f j (y j (t τ(t I i k j (t sf j (y j (sds u i (t, i Λ, y i (s = φ i (s, s (, ], k j (t sf j (y j (sds (23 where u is given by Eq. (4 and the initial condition is φ(s = ( 4.5, 3, 4.8 T, s (, ]. By using the Matlab LMI toolbox to solve the LMI (6 in Theorem 1, it can be found that the LMI is feasible and.27 P =.39.28 1.533 R = 1.5291 1.5443.27..9 Q =..33.9.9.9.36.439.492.2845 W = 1 3.492.5644.2946.2845.2946.781.97.12.7 M 1 =.12.13.73.7.73.165.9241.53.184 M 2 =.53 1.913.215.184.215.8319 3.786.35.1892 T 11 =.35 4.73.2625.1892.2625 3.3565.3225.7.189 T 12 =.21.4482.261.139.43.327.327.6.79 T 22 =.6.493.88.79.88.365.6533 U 1 =.71.5919.284 U 2 =.444.255 µ 1 =.5531, µ 2 =.5576. Consequently, the controller gain matrix H is designed as follows: 4842-9

Chin. Phys. B Vol. 21, No. 4 (212 4842.9949.17.3447 H = P 1 Q =.76.8578.2225. (24.3422.314 1.2978 By Theorem 1, models (22 and (23 are asymptotically synchronized. The simulation results are illustrated in Figs. 1(a 1(d in which the controller designed in Eq. (24 is applied. 2 (a 4 1-1 -2 1 2 3 4 5 2-2 (b -4 1 2 3 4 5 1 5-5 (c -1 1 2 3 4 5 2 1-1 (d -2 1 2 3 4 5 Fig. 1. State trajectories and error trajectories of drive system (22 and response system (23 with control input (24: (a the state x 1 and its estimation, (b the state x 2 and its estimation, (c the state x 3 and its estimation, (d the error states. 5 1-1 -5 (a -2 1 2 3 4 5 (b -1 1 2 3 4 5 15 1 (c 3 2 (d 5 1-5 -1-1 1 2 3 4 5-2 1 2 3 4 5 Fig. 2. State trajectories and error trajectories of drive system (22 and response system (23 without control input: (a the state x 1 and its estimation, (b the state x 2 and its estimation, (c the state x 3 and its estimation, (d the error states. 4842-1

Remark 3 The simulation results can be described as follows. Figures 2(a 2(d provide the state trajectories and the error trajectories between the drive system (22 and the response system (23 without control input. One may observe that the drive system (22 and the response system (23 without control input cannot be synchronized. Figures 3(a and 3(b describe the chaotic behavior in phase space of the drive system (22 and the response system (23 without control input, respectively. Figures 1(a 1(d show that the state trajectories and error trajectories of the drive system (22 and the response system (23 with control input (24 to be asymptotically synchronized. In the simulations, we choose the time step size h =.1 and time segment T = 5. In addition, we should point out that for simplicity of our computer simulations, the delay kernel k is used as k(s =.17 e s for s, 4] and k(s = for s > 4. x3 t y3 t 1-1 4 2-2 4 (a 2 x 2 t (b 2 Chin. Phys. B Vol. 21, No. 4 (212 4842-2 5 1-4 x 1 t 15-4-2-15 -1-5 -2 y 2 t y 1 t Fig. 3. (a The chaotic behavior of drive system (22 in phase space with the initial condition ϕ(s = (1,.5, 1 T, s (, ]; (b The chaotic behavior of response system (23 in phase space without control input with the initial condition φ(s = ( 4.5, 3, 4.8 T, s (, ]. 5. Conclusion In this paper, the synchronization for FCNNs with mixed delays was considered. Mixed delays include discrete time-varying delays and unbounded distributed delays. A dynamic control scheme was proposed to achieve the synchronization between a drive network and a response network. By constructing the Lyapunov Krasovskii functional which contains a triple-integral term and the free-weighting matrices method, we derived a simple and efficient criterion in terms of LMIs for synchronization. The controller can be easily obtained by solving the derived LMIs. The assumptions A1 and A2 on activation function and time-varying delay in this paper are weaker than those generally used in the previous literature, 21,24] namely, the boundedness of the activation function f j. Further the differentiability of the time varying delay τ(t is not required in this paper. A numerical example and its simulations are presented to illustrate the effectiveness of the proposed method. References 1] Chua L O and Yang L 1988 IEEE Trans. Circuits Syst. 35 1257 2] Chua L O and Yang L 1988 IEEE Trans. Circuits Syst. 35 1273 3] Yang T, Yang L B, Wu C W and Chua L O 1996 Proceedings of the IEEE International Workshop on Cellular Neural Networks and Applications (Singapore: World Scientific p. 181 4] Yang T, Yang L B, Wu C W and Chua L O 1996 Proceedings of the IEEE International Workshop on Cellular Neural Networks and Applications (Singapore: World Scientific p. 225 5] Carroll T L and Pecora L M 1991 IEEE Trans. Circuits Syst. I 38 453 6] Sun J, Liu G P, Chen J and Rees D 29 Phys. Lett. A 373 342 7] Balasubramaniam P, Kalpana M and Rakkiyappan R Circuits Systems Signal Process. 3 1595 8] Blythe S, Mao X and Liao X 21 J. Franklin Inst. 338 481 9] Gao X, Zhong S and Gao F 29 Nonlinear Anal. 71 23 1] Liu Y, Wang Z and Liu X 28 Int. J. Comput. Math. 85 1299 11] Zhang Y, Xu S and Chu Y 211 Int. J. Comput. Math. 88 249 12] Zhang C, He Y and Wu M 21 Neurocomputing 74 265 13] Yu W, Cao J and Lu W 21 Neurocomputing 73 858 14] Li T, Fei S, Zhu Q and Cong S 28 Neurocomputing 71 35 15] Li T, Song A, Fei S and Guo Y 29 Nonlinear Anal. 71 2372 16] Li X and Bohner M 21 Math. Comput. Modelling 52 643 17] Huang H and Feng G 29 Neural Networks 22 869 18] Li X and Fu X 211 Commun. Nonlinear Sci. Numer. Simul. 1 885 19] Park Ju H 29 Chaos, Solitons and Fractals 42 1299 2] Wang K, Teng Z and Jiang H 21 Math. Comput. Modelling 52 12 21] Tang Y and Fang J 29 Neurocomputing 72 3253 22] Xia Y, Yang Z and Han M 29 Commun. Nonlinear Sci. Numer. Simul. 14 3645 23] Yu F and Jiang H 211 Neurocomputing 74 59 24] Feng X, Zhang F and Wang W 211 Chaos, Solitons and Fractals 44 9 25] Boyd S, Ghaoui L E, Feron E and Balakrishnan V 1994 Linear Matrix Inequalities in System and Control Theory (Philadelphia, Pennsylvania: SIAM 26] Yang T and Yang L B 1996 IEEE Trans. Circuits Syst. I 43 88 27] Liu Z, Zhang H and Wang Z 29 Neurocomputing 72 156 28] Gu K 2 Proceedings of the 39th IEEE Conference on Decision and Control December 2 Sydney, Australia p. 285 29] Li T, Fei S and Zhu Q 29 Nonlinear Anal. Real World Appl. 1 1229 4842-11