Outer Synchronization between Fractional-Order Complex Networks: A Non-Fragile Observer-based Control Scheme

Meichun Zhao <sup>1</sup> and Junwei Wang <sup>2</sup>, \*


*Received: 25 February 2013; in revised form: 4 April 2013 / Accepted: 8 April 2013 / Published: 15 April 2013*

Abstract: This paper addresses the global outer synchronization problem between two fractional-order complex networks coupled in a drive-response configuration. In particular, for a given fractional-order complex network composed of Lur'e systems, an observer-type response network with non-fragile output feedback controllers is constructed. Both additive and multiplicative uncertainties that perturb the control gain matrices are considered. Then, using the stability theory of fractional-order systems and eigenvalue distribution of the Kronecker sum of matrices, we establish some sufficient conditions for global outer synchronization. Interestingly, the developed results are cast in the format of linear matrix inequalities (LMIs), which can be efficiently solved via the MATLAB LMI Control Toolbox. Finally, numerical simulations on fractional-order networks with nearest-neighbor and small-world topologies are given to support the theoretical analysis.

Keywords: outer synchronization; fractional-order derivative; complex network; observer control; LMI

#### 1. Introduction

Most real systems in nature, society and engineering can be properly described by models of complex networks of interacting dynamical units with diverse topologies [1]. In complex networks, various collective behaviors can emerge through the interaction of the units, among which synchronization represents one of the most interesting ones. Since the first observation of synchronization phenomenon of two pendulum clocks by Huygens in 1665, this phenomenon has been discovered in many biological and physical systems, such as pacemaker cells in the heart and nervous systems, synchronously flashing fireflies, networks of neurons in the circadian pacemaker [2,3]. Another topic that is closely related to the synchronization of complex networks is the consensus of multi-agent systems, which means that a team of agents reaches an agreement on certain quantities of interest through local communication. The study of consensus problem not only helps us understand natural phenomena (e.g., schooling of fish, flocking of birds and swarming of bees), but also has a variety of engineering applications (e.g., cooperative control of unmanned aerial vehicles, rendezvous of mobile robots and communication among sensor networks) [4]. Recently, the relation between synchronization of complex networks and consensus of multi-agent systems has been discussed [5,6]. The past few years have witnessed dramatic advances concerned with synchronization of complex networks and consensus of multi-agent systems (for more details, see [7–17] and references therein).

Notwithstanding the vast technical literature on synchronization, the great majority of research efforts has focused on complex networks of coupled integer-order systems, whose dynamics are described by integer-order differential equations. However, it has been recognized that many physical systems are more suitable to be modeled by fractional-order differential equations (*i.e*., differential equations involving fractional-order derivatives) rather than the classic integer-order ones [18,19]. Fractional-order derivatives provide an excellent instrument for the description of memory and hereditary properties of various materials and processes. Moreover, they include traditional integer-order derivatives as a special case. In addition, many natural collective behaviors can be explained by the complex networks with fractional-order dynamics: for example, the synchronized motion of agents in fractional circumstances, such as macromolecule fluids and porous media [20], CaoY2010, ShenJ2012. Therefore, it is meaningful to study synchronization problem in complex networks of coupled fractional-order systems. For convenience, we here call networks composed of integer-order and fractional-order systems "integer-order complex networks (ICNs)" and "fractional-order complex networks (FCNs)", respectively.

In recent years, synchronization in FCNs has begun to receive research attention within the scientific community [23–28], with the first systematic studies on synchronization of FCNs emerging in [29]. A common underlying assumption in the above mentioned literature is that they describe the synchronization behavior inside a single network, which has been termed "inner synchronization". For this case, it has been shown that stability of the synchronized state depends on the details of the underlying network topology. The route to inner synchronization differs from the synchronization taking place between two coupled networks. In the latter case, also known as "outer synchronization", the corresponding nodes of two coupled networks will achieve synchronization [30]. In the real world, there are a great many examples about relationships between different networks, such as the original spreading of infectious diseases between two communities, the balance of beneficial bacteria and pathogenic bacteria in our digestive systems, predator-prey interactions in ecological systems. This shows the great importance and challenge to study the dynamics between two coupled networks.

Despite the fact that some advances have been made for outer synchronization [31–38], these literatures almost exclusively report outer synchronization between coupled ICNs. For outer synchronization between two FCNs, the related research work is just beginning and the only two attempts is the work in [39,40]. In [39], the authors treated outer synchronization problem between two different bidirectionally coupled FCNs. However, the stability condition in that work depends on the eigenvalues of a large system matrix. From the viewpoint of computational complexity, the eigenvalue computation of such large size matrix is difficult and even prohibitive. Another work in [40] considered the robust outer synchronization in a setup consisting of two FCNs coupled unidirectionally by way of an open-plus-closed-loop scheme. Although their approach avoids the need to compute eigenvalues of a large system matrix, the outer synchronization behavior is achieved in the "*local*" sense, not in the "*global*" sense. When initial conditions of two FCNs lie far away from each other, the two FCNs in [40] would fail to achieve outer synchronization. Therefore, whether outer synchronization behavior between two FCNs can be achieved globally still remains an open problem, which motivates the research of this work. Another motivation comes from concerns on controller gain variations. In practical applications, uncertainties or inaccuracies do occur in the controller implementation stage due to finite word length, round-off errors in numerical computations and finite resolution measuring instruments [41]. Consequently, even though the designed controllers are robust with respect to system uncertainties, they may be very fragile to their own uncertainties. This brings the so-called non-fragile controller design problem, *i.e*., how to design a controller such that the controller is insensitive (or resilient) to some uncertainties in its coefficients. To the best of authors' knowledge, the non-fragile control problem for global outer synchronization in coupled FCNs has never been reported in the literature. The above situation is exactly what concerns and interests us.

The aim of this paper is to discuss the global outer synchronization problem between coupled FCNs by designing appropriate non-fragile controllers. We will first construct two drive-response coupled FCNs from the viewpoint of observer theory. Then, using the stability theory of fractional-order systems and the characteristics of the eigenvalue distribution of Kronecker sum of two matrices, we present a basic theorem for outer synchronization. Based on this basic theorem, two sufficient conditions for outer synchronization in the LMI format are derived for the additive and multiplicative controller gain perturbations, respectively. Finally, the effectiveness and feasibility of the designed control strategy for outer synchronization are demonstrated by numerical simulations on chaotic drive-response FCNs with nearest-neighbor or small-world topologies.

The organization of this work is as follows. In Section 2, the preliminaries and problem statement are introduced. The outer synchronization analysis for two designed coupled FCNs is studied and synchronization criteria are proposed in Section 3. In Section 4, corresponding numerical simulations

*Notations.* In the sequel, if not explicitly stated, matrices are assumed to have compatible dimensions for algebraic operations. R<sup>n</sup> denotes the n-dimensional Euclidean space, R<sup>m</sup>×<sup>n</sup> is the set of all <sup>m</sup> <sup>×</sup> <sup>n</sup> real matrices; *<sup>I</sup>*<sup>n</sup> represents the <sup>n</sup> <sup>×</sup> <sup>n</sup> identity matrix, *<sup>I</sup>* means an identity matrix of appropriate order; the shorthand diag{· · · } represents a block diagonal matrix; and the superscript "T" stands for matrix transposition. <sup>λ</sup>min(*A*) and <sup>λ</sup>max(*A*) denote the smallest and the largest eigenvalues of a real symmetric matrix *A*, respectively. The operator Sym(*A*) denotes *A*+*A*<sup>T</sup> . The notation *<sup>P</sup>* <sup>&</sup>gt; 0 (<sup>≥</sup> <sup>0</sup>, < <sup>0</sup>, <sup>≤</sup> 0) means that <sup>P</sup> is positive (semi-positive, negative, semi-negative) definite, and the symbol denotes the elements below the main diagonal of a symmetric block matrix. For *<sup>A</sup>* <sup>∈</sup> <sup>R</sup><sup>m</sup>×<sup>n</sup> and *<sup>B</sup>* <sup>∈</sup> <sup>R</sup><sup>p</sup>×<sup>q</sup>, *<sup>A</sup>*⊗*<sup>B</sup>* <sup>∈</sup> <sup>R</sup>mn×pq denotes the Kronecker product of the two matrices.

#### 2. Preliminaries and Network Model

are presented. Finally, conclusions are drawn in Section 5.

#### *2.1. Basic Concepts and Lemmas*

In this section, we first recall some definitions related to fractional-order derivatives that will be used throughout this paper. Note that there are different definitions for fractional derivatives, such as Grünwald–Letnikov, Riemann–Liouville, and Caputo definitions [18]. Here we make use of the following Caputo fractional derivative, because its Laplace transform allows utilization of initial conditions of classical integer-order derivative that has known physical interpretations.

*Definition 2.1.* Let <sup>m</sup> <sup>−</sup> <sup>1</sup> < α <sup>≤</sup> <sup>m</sup>, with <sup>m</sup> <sup>∈</sup> <sup>N</sup>, the Caputo fractional derivative is defined as

$$D\_t^\alpha h(t) = J\_t^{m-\alpha} D\_t^m h(t) = \begin{cases} \frac{1}{\Gamma(m-\alpha)} \int\_0^t (t-\tau)^{m-\alpha-1} h^{(m)}(\tau) d\tau, \ \alpha \neq m\\ h^{(m)}(t), \ \alpha = m \end{cases}$$

where h(m) (t) = <sup>d</sup>mh(t) dt<sup>m</sup> is the ordinary derivative of integer order m.

To prove the main results in the next section, we need the following lemmas. Lemma 2.2 ([42]). For the linear fractional-order system

$$D\_t^\alpha x(t) = A x(t), \text{with } x(0) = x\_0 \tag{1}$$

where <sup>0</sup> <α< <sup>1</sup>, *<sup>x</sup>* <sup>∈</sup> <sup>R</sup><sup>n</sup> and *<sup>A</sup>* <sup>∈</sup> <sup>R</sup><sup>n</sup>×<sup>n</sup>. System (1) is asymptotically stable if and only if

$$\left| \arg(\text{spec}(\mathbf{A})) \right| > \frac{\alpha \pi}{2} \tag{2}$$

where spec(*A*) is the spectrum of system matrix *A*.

Lemma 2.3 ([43]). Let *<sup>A</sup>* <sup>∈</sup> <sup>R</sup><sup>n</sup>×<sup>n</sup> have eigenvalues <sup>λ</sup><sup>i</sup> (1 <sup>≤</sup> <sup>i</sup> <sup>≤</sup> <sup>n</sup>), and let *<sup>B</sup>* <sup>∈</sup> <sup>R</sup><sup>m</sup>×<sup>m</sup> have eigenvalues <sup>μ</sup><sup>i</sup> (1 <sup>≤</sup> <sup>i</sup> <sup>≤</sup> <sup>m</sup>). Then the Kronecker sum *<sup>A</sup>* <sup>⊕</sup> *<sup>B</sup>* - (*I*<sup>m</sup> <sup>⊗</sup> *<sup>A</sup>*)+(*<sup>B</sup>* <sup>⊗</sup> *<sup>I</sup>*n) <sup>∈</sup> <sup>R</sup>mn×mn has mn eigenvalues

$$\left(\lambda\_1 + \mu\_1, \dots, \lambda\_1 + \mu\_m, \lambda\_2 + \mu\_1, \dots, \lambda\_2 + \mu\_m, \dots, \lambda\_n + \mu\_1, \dots, \lambda\_n + \mu\_m\right)$$

Lemma 2.4 ([44]). Let *<sup>A</sup>* <sup>∈</sup> <sup>R</sup><sup>n</sup>×<sup>n</sup> and <sup>0</sup> <α< <sup>1</sup>. The fractional-order system <sup>D</sup><sup>α</sup> <sup>∗</sup> *<sup>x</sup>*(t) = *Ax*(t) is asymptotically stable if and only if there exist two real symmetric positive definite matrices *<sup>P</sup>*k<sup>1</sup> <sup>∈</sup> <sup>R</sup><sup>n</sup>×<sup>n</sup> (<sup>k</sup> = 1, 2), and two skew-symmetric matrices *<sup>P</sup>*k<sup>2</sup> <sup>∈</sup> <sup>R</sup><sup>n</sup>×<sup>n</sup> (<sup>k</sup> = 1, 2), such that

$$\sum\_{i=1}^{2} \sum\_{j=1}^{2} \text{Sym}\{\Theta\_{ij} \otimes (\mathbf{AP}\_{ij})\} < 0 \tag{3}$$

$$
\begin{bmatrix}
\mathbf{P}\_{11} & \mathbf{P}\_{12} \\
\end{bmatrix} > 0,
\begin{bmatrix}
\mathbf{P}\_{21} & \mathbf{P}\_{22} \\
\end{bmatrix} > 0\tag{4}
$$

where

$$\begin{aligned} \Theta\_{11} &= \begin{bmatrix} \sin \theta & -\cos \theta \\ \cos \theta & \sin \theta \end{bmatrix}, \Theta\_{12} = \begin{bmatrix} \cos \theta & \sin \theta \\ -\sin \theta & \cos \theta \end{bmatrix} \\ \Theta\_{21} &= \begin{bmatrix} \sin \theta & \cos \theta \\ -\cos \theta & \sin \theta \end{bmatrix}, \Theta\_{22} = \begin{bmatrix} -\cos \theta & \sin \theta \\ -\sin \theta & -\cos \theta \end{bmatrix}, \theta = \frac{q\pi}{2} \end{aligned}$$

Lemma 2.5 ([45]). For matrices *X* and *Y* with appropriate dimensions, the following inequality holds for ε > 0

$$\mathbf{X}^T \mathbf{Y} + \mathbf{Y}^T \mathbf{X} < \varepsilon \mathbf{X}^T \mathbf{X} + \frac{1}{\varepsilon} \mathbf{Y}^T \mathbf{Y} \tag{5}$$

Lemma 2.6 ([46]). For a given symmetric matrix *S* = *S*<sup>T</sup> , the following assertions are equivalent


#### *2.2. Network Model*

In this paper, we will focus on the outer synchronization problem of two coupled fractional-order complex networks with the drive-response (or master-slave, unidirectional) coupling structure, in which the drive network does not receive any information from the response network.

The drive network with each node being a n-dimensional fractional-order differential system in the form of Lur'e type is described as

$$\begin{aligned} D\_t^\alpha x\_i &= \mathbf{A}x\_i + \mathbf{f}(Hx\_i, t) + \sum\_{j=1}^N c\_{ij} x\_j\\ y\_i &= \mathbf{H}x\_i, \ i = 1, 2, \cdots, N\end{aligned} \tag{6}$$

where <sup>0</sup> <α< <sup>1</sup> is the fractional commensurate order, *<sup>x</sup>*<sup>i</sup> <sup>∈</sup> <sup>R</sup><sup>n</sup> and <sup>y</sup><sup>i</sup> <sup>∈</sup> <sup>R</sup> denote the state and scalar output vectors of the <sup>i</sup>-th node, respectively. The nonlinear vector-valued function *<sup>f</sup>*(·, ·) : <sup>R</sup> <sup>×</sup> <sup>R</sup> −→ <sup>R</sup><sup>n</sup> is continuously differentiable. The constant matrix *<sup>A</sup>* <sup>∈</sup> <sup>R</sup><sup>n</sup>×<sup>n</sup> combing with *<sup>f</sup>*(·, ·) describes the dynamics of individual nodes. *<sup>H</sup>* <sup>∈</sup> <sup>R</sup><sup>1</sup>×<sup>n</sup> is the observer matrix. The matrix *<sup>C</sup>* <sup>=</sup> (cij ) <sup>∈</sup> <sup>R</sup><sup>N</sup>×<sup>N</sup> is the outer coupling matrix of the drive network, which denotes the network topology, and is defined as follows: if there is a connection between node i and node j (i = j), then cij = <sup>c</sup>ji <sup>&</sup>gt; <sup>0</sup>; otherwise, <sup>c</sup>ij <sup>=</sup> <sup>c</sup>ji = 0 (<sup>j</sup> <sup>=</sup> <sup>i</sup>), and the diagonal elements of matrix *<sup>C</sup>* are defined by <sup>c</sup>ii <sup>=</sup> <sup>−</sup><sup>N</sup> j=1 j=i <sup>c</sup>ij , i = 1, <sup>2</sup>, ··· , N. It is noted that the configuration matrix *<sup>C</sup>* is symmetric with non-positive real eigenvalues and not necessarily irreducible.

Based on the design idea of observer in the control theory, the response network is described as a nonlinear observer

$$\begin{aligned} D\_t^\alpha \hat{x}\_i &= \mathbf{A} \hat{x}\_i + \mathbf{f}(y\_i, t) + \sum\_{j=1}^N c\_{ij} \hat{x}\_j + \mathbf{U}\_i(y\_i, \hat{y}\_i) \\ \hat{y}\_i &= \mathbf{H} \hat{x}\_i, \ i = 1, 2, \cdots, N \end{aligned} \tag{7}$$

where *<sup>x</sup>***ˆ**<sup>i</sup> <sup>∈</sup> <sup>R</sup><sup>n</sup> and <sup>y</sup>ˆ<sup>i</sup> <sup>∈</sup> <sup>R</sup> denote the state and scalar output vectors of node <sup>i</sup> of the response network, respectively. *<sup>A</sup>*, *<sup>f</sup>*(·, ·), *<sup>H</sup>* and *<sup>C</sup>* are the same as in system Equation (6). *<sup>U</sup>*i(·, ·) : <sup>R</sup> <sup>×</sup> <sup>R</sup> −→ <sup>R</sup><sup>n</sup> is the control input to be designed.

Considering the control gain perturbations, the actual implemented control input is assumed to be

$$U\_i(y\_i, \hat{y}\_i) = (\mathbf{L} + \Delta \mathbf{L}(t))(y\_i - \hat{y}\_i) \tag{8}$$

where *<sup>L</sup>* <sup>∈</sup> <sup>R</sup><sup>n</sup> is the nominal control gain matrix and the term <sup>Δ</sup>*L*(t) represents the control gain variations. In this paper, the following two classes of the control gain variations are considered: Type 1: Δ*L*(t) is with the norm-bounded additive form:

$$
\Delta L(t) = \Delta\_1(t) = \mathbf{M}\_1 \mathbf{F}\_1(t) \mathbf{N}\_1 \tag{9}
$$

Type 2: Δ*L*(t) is with the norm-bounded multiplicative form:

$$
\Delta L(t) = \Delta\_2(t) = \mathbf{M}\_2 \mathbf{F}\_2(t) \mathbf{N}\_2 \mathbf{L} \tag{10}
$$

where *M*1, *N*1, *M*<sup>2</sup> and *N*<sup>2</sup> are known matrices with appropriate dimensions, *F*1(t) and *F*2(t) are unknown time-varying matrices satisfying the relation

$$\mathbf{F}\_1(t)\mathbf{F}\_1^T(t) \le \mathbf{I}, \; \mathbf{F}\_2(t)\mathbf{F}\_2^T(t) \le \mathbf{I} \tag{11}$$

Let us define the errors of outer synchronization *<sup>e</sup>*<sup>i</sup> <sup>=</sup> *<sup>x</sup>*<sup>i</sup> <sup>−</sup> *<sup>x</sup>*ˆi, then the following error dynamics of the outer synchronization can be obtained from Equations (6) and (7):

$$D\_t^\alpha \mathbf{e}\_i = (\mathbf{A} - (\mathbf{L} + \Delta \mathbf{L})\mathbf{H})\mathbf{e}\_i + \sum\_{j=1}^N c\_{ij}\mathbf{e}\_j, \ i = 1, 2, \cdots, N \tag{12}$$

In this paper, we aim at establishing computable criteria in the LMI format to find the control gain *L* such that the coupled fractional-order complex networks Equations (6) and (7) achieve global outer synchronization in the following sense

$$\lim\_{t \to \infty} \|\mathbf{e}\_i(t)\| = \lim\_{t \to \infty} \|\mathbf{x}\_i(t) - \hat{\mathbf{x}}\_i(t)\| = 0, i = 1, 2, \cdots, N \tag{13}$$

for any initial conditions *<sup>x</sup>*i(0) and *<sup>x</sup>*ˆi(0), where · refers to the Euclidean norm.

Using the Kronecker product, the error system Equation (12) can be compactly rewritten equivalently as

$$\mathbf{d}\,D\_t^\alpha \mathbf{e} = \left(\mathbf{I}\_N \otimes \mathbf{A}\_L + \mathbf{C} \otimes \mathbf{I}\_n\right) \mathbf{e}(t) \tag{14}$$

where *e* = \* *e*T <sup>1</sup> , *e*<sup>T</sup> <sup>2</sup> , ··· , *<sup>e</sup>*<sup>T</sup> N +T <sup>∈</sup> *<sup>R</sup>*nN and *<sup>A</sup><sup>L</sup>* <sup>=</sup> *<sup>A</sup>* <sup>−</sup> (*<sup>L</sup>* + Δ*L*)*H*. It is implied from the above representation that the global outer synchronization problem between the FCNs Equations (6) and (7) is converted into an equivalent global asymptotical stability problem of the linear error system Equation (14).

#### 3. Global Outer Synchronization Analysis

From Lemma 2.2, the error system Equation (14) is asymptotically stable if the spectrum of system matrix *<sup>I</sup>*<sup>N</sup> <sup>⊗</sup> *<sup>A</sup><sup>L</sup>* <sup>+</sup> *<sup>C</sup>* <sup>⊗</sup> *<sup>I</sup>*<sup>n</sup> satisfy the inequality (2). It can been seen that the system matrix *<sup>I</sup>*<sup>N</sup> <sup>⊗</sup> *<sup>A</sup><sup>L</sup>* <sup>+</sup> *<sup>C</sup>* <sup>⊗</sup> *<sup>I</sup>*<sup>n</sup> is the Kronecker sum of *<sup>A</sup><sup>L</sup>* and *<sup>C</sup>*, then its eigenvalues can be expressed by the sum of eigenvalues of *A<sup>L</sup>* and eigenvalues of *C* using Lemma 2.3. Hence, we obtain the following basic theorem.

Theorem 3.1. The fractional-order complex networks Equations (6) and (7) with the designed non-fragile controllers (8) will achieve outer synchronization behavior globally, if the following condition is satisfied

$$\left| \arg(\text{spec}(\mathbf{A}\_L)) \right| > \frac{\alpha \pi}{2} \tag{15}$$

where spec(*AL*) is the spectrum of system matrix *<sup>A</sup><sup>L</sup>* <sup>=</sup> *<sup>A</sup>* <sup>−</sup> (*<sup>L</sup>* + Δ*L*)*H*.

*Remark 3.2.* The importance of this theorem lies in the fact that it converts the outer synchronization problem between coupled FCNs Equations (6) and (7) into the eigenvalue distribution of the uncertain matrix *A<sup>L</sup>* with the same dimension as a single node, thereby significantly reducing the computational complexity. In addition, it should be noted that in previous work [40] on outer synchronization of FCNs, the linear error systems are often derived through a suitable linearization of the system's nonlinear functions. This approach often implies "*local*" stability of the outer synchronization manifold. However, in many applications, global stability of the outer synchronization manifold is very desirable but difficult to achieve. To circumvent this difficulty, here an observer-based design procedure is performed instead of linearizing approximation of the nonlinear function *<sup>f</sup>*(·, ·) in Equations (6) and (7). Thus the above proposed condition theoretically guarantees the "*global*" outer synchronization between FCNs Equations (6) and (7).

*Remark 3.3.* For a given pair (*H*, *A*), whether the observer-type response network Equation (7) synchronizes the drive network Equation (6) globally depends on whether the spectrum of system matrix *A<sup>L</sup>* satisfies the inequality (15). In case when the control gain *L* is deterministic without uncertainty, *i.e*., <sup>Δ</sup>*<sup>L</sup>* = 0, the control gain *<sup>L</sup>* <sup>∈</sup> <sup>R</sup><sup>n</sup> may be chosen such that the inequality (15) is satisfied. As is known from the control theory [47], if the pair (*H*, *A*) is observable, *i.e*., if rank[*H*<sup>T</sup> , *<sup>A</sup>*T*H*<sup>T</sup> , ··· ,(*A*<sup>T</sup> )<sup>n</sup>−<sup>1</sup>*H*<sup>T</sup> ] = <sup>n</sup>, then there exists *<sup>L</sup>* providing the matrix *<sup>A</sup><sup>L</sup>* <sup>=</sup> *<sup>A</sup>* <sup>−</sup> *LH* with any given eigenvalues. Particularly, all eigenvalues of *A<sup>L</sup>* can be designed to locate the region defined by inequality (15). For other general pair (*H*, *A*), the control gain *L* is usually determined through a considerable amount of trial and error. However, with control gain uncertainties (9) or (10), the condition (15) is not easy to be checked because there is infinite number of eigenvalues [48]. To effectively avoid this difficulty, in the following we will develop an LMI-based design method.

Now we are in a position to present an LMI-type solvability condition for the outer synchronization problem of coupled FCNs (6) and (7) with the control gain perturbations (9) or (10). Theorem 3.4. Consider the fractional-order complex networks (6) and (7) with the control gain perturbation Δ1(t) in Equation (9). Then the outer synchronization between networks (6) and (7) will be achieved globally, if there exist a constant ε > 0, a symmetric positive definite matrix *P* > 0 and a matrix *Q*, satisfying the following LMI:

$$
\begin{bmatrix}
\Psi & I\_2 \otimes (N\_1HP)^T \\
\star & -\frac{1}{2}\varepsilon I
\end{bmatrix} < 0 \tag{16}
$$

where

$$\Psi = \sum\_{i=1}^{2} \text{Sym}\left\{\Theta\_i \otimes (\mathbf{AP} - \mathbf{Q})\right\} + 2\varepsilon \left(\mathbf{I}\_2 \otimes \mathbf{M}\_1 \mathbf{M}\_1^T\right)^2$$

with

$$\begin{aligned} \Theta\_1 = \begin{bmatrix} \sin \theta & -\cos \theta \\ \cos \theta & \sin \theta \end{bmatrix}, \Theta\_2 = \begin{bmatrix} \sin \theta & \cos \theta \\ -\cos \theta & \sin \theta \end{bmatrix}, \theta = \frac{\alpha \pi}{2} \end{aligned}$$

Moreover, the non-fragile control gain *L* is given by

$$L = QP^{-1}H^T \left(HH^T\right)^{-1} \tag{17}$$

*Proof.* Setting *P*<sup>11</sup> = *P*<sup>21</sup> = *P* , *P*<sup>12</sup> = *P*<sup>22</sup> = **0** in Lemma 2.4, we have that if there exists real symmetric positive definite matrix *P* such that

$$\sum\_{i=1}^{2} \text{Sym}\{\Theta\_i \otimes (A\_L P)\} < 0 \tag{18}$$

then <sup>|</sup>arg(spec(*AL*))<sup>|</sup> <sup>&</sup>gt; απ <sup>2</sup> , where spec(*AL*) is the spectrum of system matrix *<sup>A</sup>L*.

With *<sup>A</sup><sup>L</sup>* <sup>=</sup> *<sup>A</sup>* <sup>−</sup> (*<sup>L</sup>* + Δ*L*)*H*, the left hand side of Equation (18) can be rewritten as

$$\begin{aligned} \sum\_{i=1}^{2} \text{Sym} \{ \Theta\_i \otimes (\mathbf{A}\_L \mathbf{P}) \} &= \sum\_{i=1}^{2} \text{Sym} \{ \Theta\_i \otimes (\mathbf{AP} - LHP - \Delta LHP) \} \\ &= \sum\_{i=1}^{2} \text{Sym} \{ \Theta\_i \otimes (\mathbf{AP} - LHP) \} + \sum\_{i=1}^{2} \text{Sym} \{ \Theta\_i \otimes (-M\_1 \mathbf{F}\_1(t) N\_1 HP) \} \end{aligned} (19)$$

From Equation (11), one has

$$\begin{aligned} (I\_2 \otimes F\_1)(I\_2 \otimes F\_1)^T &= (I\_2 \otimes F\_1)(I\_2 \otimes F\_1^T) \\ &= I\_2 \otimes (F\_1 F\_1^T) \le I \end{aligned} \tag{20}$$

By Equation (20) and Θi<sup>1</sup>Θ<sup>T</sup> <sup>i</sup><sup>1</sup> = *I*2, it directly follows from Lemma 2.5 that for any real scalar ε > 0

$$\begin{aligned} &\sum\_{i=1}^{2} \text{Sym}\{\Theta\_{i}\otimes(-M\_{1}\mathbf{F}\_{1}(t)N\_{1}HP)\} \\ &= \sum\_{i=1}^{2} \text{Sym}\{- (\Theta\_{i}\otimes M\_{1})(I\_{2}\otimes \mathbf{F}\_{1})(I\_{2}\otimes (N\_{1}HP))\} \\ &\leq \sum\_{i=1}^{2} \left\{ \varepsilon (\Theta\_{i}\otimes M\_{1})(I\_{2}\otimes \mathbf{F}\_{1})(I\_{2}\otimes \mathbf{F}\_{1})^{T}(\Theta\_{i}\otimes M\_{1})^{T} + \frac{1}{\varepsilon}(I\_{2}\otimes N\_{1}HP)^{T}(I\_{2}\otimes N\_{1}HP) \right\} \\ &\leq 2\varepsilon \left(I\_{2}\otimes M\_{1}M\_{1}^{T}\right) + \frac{2}{\varepsilon}(I\_{2}\otimes N\_{1}HP)^{T}(I\_{2}\otimes N\_{1}HP) \end{aligned} \tag{21}$$

Substituting Equation (21) into Equation (19), one has

$$\begin{aligned} \sum\_{i=1}^{2} \text{Sym}\left\{\Theta\_{i} \otimes (\mathbf{A}\_{L}\mathbf{P})\right\} &= \sum\_{i=1}^{2} \text{Sym}\left\{\Theta\_{i} \otimes (\mathbf{AP} - \mathbf{L}HP - \Delta LHP)\right\} \\ &\leq \sum\_{i=1}^{2} \text{Sym}\left\{\Theta\_{i} \otimes (\mathbf{AP} - \mathbf{L}HP)\right\} \\ &+ 2\varepsilon \left(\mathbf{I}\_{2} \otimes \mathbf{M}\_{1}\mathbf{M}\_{1}^{T}\right) + \frac{2}{\varepsilon} (\mathbf{I}\_{2} \otimes \mathbf{N}\_{1}\mathbf{HP})^{T} (\mathbf{I}\_{2} \otimes \mathbf{N}\_{1}\mathbf{HP}) \end{aligned} (22)$$

Let

$$Q = LHP \tag{23}$$

and following from Equation (22), the inequality Equation (18) holds if

$$\sum\_{i=1}^{2} \text{Sym}\left\{\Theta\_i \otimes (AP - \mathbf{Q})\right\} + 2\varepsilon \left(I\_2 \otimes M\_1 M\_1^T\right) + \frac{2}{\varepsilon} (I\_2 \otimes N\_1 HP)^T (I\_2 \otimes N\_1 HP) < 0 \quad (24)$$

Using Lemma 2.6, it is easily seen that Equation (24) is in turn equivalent to the linear matrix inequality Equation (16), which is the condition stated in the theorem. Therefore, <sup>|</sup>arg(spec(*AL*))<sup>|</sup> <sup>&</sup>gt; απ <sup>2</sup> , which implies that outer synchronization between the fractional-order networks (6) and (7) will occur globally by using Theorem 3.1. This completes the proof.

For the FCNs (6) and (7) with the control gain perturbation defined in Equation (10), we have the following results.

Theorem 3.5. Consider the fractional-order complex networks (6) and (7) with the control gain perturbation Δ2(t) in Equation (10). Then the outer synchronization between networks (6) and (7) will be achieved globally, if there exist a constant ε > 0, a symmetric positive definite matrix *P* > 0 and a matrix *Q*, satisfying the following LMI:

$$
\begin{bmatrix}
\Psi & I\_2 \otimes (\mathbf{N}\_2 \mathbf{Q})^T \\
\star & -\frac{1}{2} \varepsilon I
\end{bmatrix} < 0 \tag{25}
$$

where

$$\Psi = \sum\_{i=1}^{2} \text{Sym}\left\{\Theta\_i \otimes (AP - Q)\right\} + 2\varepsilon \left(I\_2 \otimes M\_2 M\_2^T\right)^2$$

with

$$\begin{aligned} \Theta\_1 = \begin{bmatrix} \sin \theta & -\cos \theta \\ \cos \theta & \sin \theta \end{bmatrix}, \Theta\_2 = \begin{bmatrix} \sin \theta & \cos \theta \\ -\cos \theta & \sin \theta \end{bmatrix}, \theta = \frac{\alpha \pi}{2} \end{aligned}$$

Moreover, the non-fragile control gain *L* is given by

$$L = QP^{-1}H^T \left(HH^T\right)^{-1} \tag{26}$$

*Proof.* It is similar to that of Theorem 3.4, the details are then omitted.

A particular case of Theorems 3.4 and 3.5 is Δ*L* = 0. In such a case, we can still provide a sufficient condition for global outer synchronization.

Corollary 3.6. Consider the fractional-order complex networks (6) and (7) without the control gain perturbation (*i.e*., Δ*L* = 0). Then the outer synchronization between networks (6) and (7) will be achieved globally, if there exist a symmetric positive definite matrix *P* > 0 and a matrix *Q*, satisfying the following LMI:

$$
\Theta\_1 \otimes (\mathbf{AP} - \mathbf{Q}) + \Theta\_1^T \otimes (\mathbf{AP} - \mathbf{Q})^T + \Theta\_2 \otimes (\mathbf{AP} - \mathbf{Q}) + \Theta\_2^T \otimes (\mathbf{AP} - \mathbf{Q})^T < 0 \tag{27}
$$

where

$$\begin{aligned} \Theta\_1 = \begin{bmatrix} \sin \theta & -\cos \theta \\ \cos \theta & \sin \theta \end{bmatrix}, \Theta\_2 = \begin{bmatrix} \sin \theta & \cos \theta \\ -\cos \theta & \sin \theta \end{bmatrix}, \theta = \frac{\alpha \pi}{2} \end{aligned}$$

Moreover, the non-fragile control gain *L* is given by

$$L = QP^{-1}H^T \left(HH^T\right)^{-1} \tag{28}$$

*Proof.* It follows directly from Theorems 3.4 and 3.5, and the details are therefore omitted here. *Remark 3.7.* The above theorems and corollary present sufficient conditions for the solvability of non-fragile outer synchronization problem for coupled fractional-order complex networks, which are related to the solutions to LMIs. In this case, these LMIs can be solved efficiently by resorting to some standard numerical algorithms [46].

#### 4. Numerical Simulations

In this section, numerical examples are given to verify the effectiveness of the above design scheme. For the control gain perturbations in additive form, the coupled fractional-order jerk model based on nearest-neighbor topology is first utilized to demonstrate the main results of Theorem 3.4. Then, a small-world complex network consisted of fractional-order Duffing oscillators is introduced to illustrate the correctness of Theorem 3.5 for the control gain perturbations in multiplicative form.

#### *4.1. Outer Synchronization between Two FCNs with Nearest-Neighbor Network Topology*

This subsection considers two coupled FCNs with N = 10 nodes each and nearest-neighbor network topology. The dynamics of the nodes in two networks is determined by the following fractional-order jerk model [49]:

$$\begin{cases} D\_t^\alpha x\_1 = x\_2 \\ D\_t^\alpha x\_2 = x\_3 \\ D\_t^\alpha x\_3 = -x\_2 - \rho x\_3 + \varphi(x\_1) \end{cases} \tag{29}$$

with nonlinear characteristic

$$
\varphi(x\_1) = -1.2x\_1 + 2\text{sgn}(x\_1).
$$

and the measured output

$$y(t) = x\_1(t)$$

where x1, x<sup>2</sup> and x<sup>3</sup> are, respectively, the position, velocity, and acceleration of the object, ρ > 0 is the control parameter. This model in its integer-order version (*i.e*., α = 1) is used to determine the time derivative of acceleration of an object and is known to give chaos for ρ = 0.6. For system (29), we show that the chaotic behavior is preserved in the fractional-order case, as shown in Figure 1 for α = 0.95. In Lur'e form, the fractional-order jerk model (29) can be represented with

$$\mathbf{A} = \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 0 & -1 & -\rho \end{pmatrix}, \mathbf{f}(y) = \begin{pmatrix} 0 \\ 0 \\ \varphi(y) \end{pmatrix}, \mathbf{H} = \{1, 0, 0\} \tag{30}$$

Figure 1. Chaotic behavior of the fractional-order jerk model (29). The fractional orders are: (a) α = 1 and (b) α = 0.95.

The additive control gain perturbations Δ*L*(t) can be described by (9) with

$$\mathbf{M}\_1 = \begin{pmatrix} 0.05 & -0.13 & 0.1 \\ -0.24 & 0.54 & 0.2 \\ -0.15 & -0.2 & 0.16 \end{pmatrix}, \mathbf{F}\_1(t) = \begin{pmatrix} \cos(3t) & 0 & 0 \\ 0 & \cos(0.1t) & 0 \\ 0 & 0 & \sin(2t) \end{pmatrix}, \mathbf{N}\_1 = \begin{pmatrix} 0.1 \\ 0.2 \\ 0.15 \end{pmatrix} \tag{31}$$

For the nearest-neighbor coupling structure, the coupling matrix is given by

$$C = \begin{pmatrix} -2k & k & 0 & \cdots & 0 & k \\ k & -2k & k & 0 & \cdots & 0 \\ & & \ddots & \ddots & \ddots & \\ 0 & \cdots & 0 & k & -2k & k \\ k & 0 & \cdots & 0 & k & -2k \end{pmatrix} \tag{32}$$

where k > 0 denotes the coupling strength of the whole network.

Since the above Δ*L*(t) is in the form of Type 1, Theorem 3.4 is used to design a non-fragile observer-based control (8). Using the MATLAB LMI Control Toolbox, we find that the LMI (16) in Theorem 3.4 is feasible. A feasible solution is presented as follows:

$$\varepsilon = 85.0407, \; \mathbf{P} = \begin{pmatrix} 41.2019 & -0.0000 & 0.0000 \\ -0.0000 & 53.1504 & -0.0000 \\ 0.0000 & -0.0000 & 53.1504 \end{pmatrix}, \; \mathbf{Q} = 10^3 \times \begin{pmatrix} 0.0146 & 2.9380 & -0.0124 \\ -2.8902 & 0.0299 & -0.4075 \\ 0.0154 & 0.4040 & -0.0148 \end{pmatrix}$$

Figure 2. Synchronization errors between the FCNs (6) and (7), where each node is a chaotic fractional-order jerk model (29). (a) The time evolutions of ei<sup>1</sup>(t) = xi<sup>1</sup>(t) − xˆi<sup>1</sup>(t); (b) the time evolutions of ei<sup>2</sup>(t) = xi<sup>2</sup>(t) − xˆi<sup>2</sup>(t); (c) the time evolutions of ei<sup>3</sup>(t) = xi<sup>3</sup>(t) − xˆi<sup>3</sup>(t) (i = 1, 2, ··· , 10).

Therefore, based on the Equation (17) in Theorem 3.4, the nominal control gain is given by:

$$L = QP^{-1}H^T \left(HH^T\right)^{-1} = \begin{pmatrix} 0.3539\\ -70.1471\\ 0.3734 \end{pmatrix} \tag{33}$$

With the aforementioned control gain matrix and k = 1, the simulation results for synchronization errors eij (i = 1, 2, ··· , 10, j = 1, 2, 3) of networks (6) and (7) are given in Figure 2, where the initial conditions *x*i(0) and *x*ˆi(0) are randomly chosen. As seen in Figure 2, the trajectories of the synchronization errors approach zero, which imply outer synchronization between complex networks (6) and (7) with fractional-order jerk models as nodes' dynamics.

#### *4.2. Outer Synchronization between Two FCNs with Small-World Network Topology*

In this simulation, two small-world FCNs of N = 100 Duffing oscillators are constructed. A single fractional-order Duffing oscillator [50] is described by:

$$\begin{cases} D\_t^\alpha x\_1 = x\_2 \\ D\_t^\alpha x\_2 = -p\_1 x\_2 - p\_2 x\_1 - p\_3 x\_1^3 + q \cos(\omega t) \end{cases} \tag{34}$$

with the measured output

$$y(t) = H\mathbf{z}(t) = (1, 0)(x\_1(t), x\_2(t))^T \tag{35}$$

where p1, p2, p<sup>3</sup> and q are system parameters. The system can be represented in Lur'e form with

$$\mathbf{A} = \begin{pmatrix} 0 & 1 \\ -p\_2 & -p\_1 \end{pmatrix}, \mathbf{f}(y, t) = \begin{pmatrix} 0 \\ -p\_3 y^3 + q \cos(\omega t) \end{pmatrix} \tag{36}$$

For the parameters p<sup>1</sup> = 1/25, p<sup>2</sup> = −1/5, p<sup>3</sup> = 8/15, q = 2/5 and ω = 0.2, the fractional-order Duffing oscillator (34) exhibits chaotic behavior for α = 0.98 (see Figure 3).

The additive controller uncertainties in Equation (10) are considered through

$$\mathbf{M}\_2 = \begin{pmatrix} -0.5 & -0.24 \\ 0.1 & -0.35 \end{pmatrix}, \ F\_2(t) = \begin{pmatrix} \sin(0.1t) & 0 \\ 0 & \cos(5t) \end{pmatrix}, \ \mathbf{N}\_2 = \begin{pmatrix} 0.2 & 0.16 \\ -0.05 & 0.32 \end{pmatrix} \tag{37}$$

Since the control gain perturbations Δ*L*(t) is in the form of Type 2, Theorem 3.5 is used to design the non-fragile observer-based control (8). Using the MATLAB LMI Control Toolbox, a feasible solution to the LMI (25) in Theorem 3.5 is given by

$$\varepsilon = 71.6499, \; \mathbf{P} = \begin{pmatrix} 40.9209 & -1.7321 \\ -1.7321 & 31.6650 \end{pmatrix}, \; \mathbf{Q} = 10^3 \times \begin{pmatrix} 0.0212 & -7.5497 \\ 7.5921 & 0.0151 \end{pmatrix}.$$

Therefore, based on the Equation (26) in Theorem 3.5, the nominal control gain is given by:

$$L = QP^{-1}H^T \left(HH^T\right)^{-1} = \begin{pmatrix} -9.5952\\ 185.9816 \end{pmatrix} \tag{38}$$

Figure 3. Chaotic behavior of the fractional-order Duffing oscillator (34). The fractional orders are: (a) α = 1 and (b) α = 0.98.

The zero-row-sum coupling matrix *C* is generated from the known Watts–Strogatz small-world network model [51] with N = 100, m = 3 and p = 0.1 . According to Theorem 3.5, the outer synchronization between two fractional-order complex networks of coupled Duffing oscillators with small-world topology will be achieved globally. Figure 4 shows the changes in synchronization errors eij (i = 1, 2, ··· , 100, j = 1, 2), respectively. From these simulation results, it can be seen the designed drive-response networks achieve outer synchronization globally and the effectiveness of the theoretical analysis is demonstrated.

Figure 4. Synchronization errors between the FCNs (6) and (7), where each node is a chaotic fractional-order Duffing oscillator (34). (a) The time evolutions of ei<sup>1</sup>(t) = xi<sup>1</sup>(t) − xˆi<sup>1</sup>(t); (b) the time evolutions of ei<sup>2</sup>(t) = xi<sup>2</sup>(t) − xˆi<sup>2</sup>(t) (i = 1, 2, ··· , 100).

#### 5. Conclusions

In this paper, we have proposed a novel observer-based control scheme for outer synchronization between two complex networks with fractional-order derivatives. The designed controllers have the following two features: (i) they use only scalar output signals to couple two FCNs in a drive-response manner; (ii) they are non-fragile for both additive and multiplicative control gain perturbations. Therefore, this is more practical and economical for real network applications, such as communication networks. Taking advantage of the eigenvalue distribution of Kronecker sum of two matrices, we presented a basic theorem for outer synchronization of coupled FCNs. Then, two sufficient conditions in the form of LMI for outer synchronization of FCNs are further provided. Compared with previous results, the proposed conditions can not only ensure outer synchronization to be achieved in the "global" sense but also facilitate it with the help of MATLAB LMI Control Toolbox.

#### Acknowledgments

The work is supported by the National Natural Science Foundation of China (Grant Nos. 61104138, 11271139), Guangdong Natural Science Foundation (Grant No. S2011040001704), the Foundation for Distinguished Young Talents in Higher Education of Guangdong, China (Grant No. LYM10074), and the Business Intelligence Key Team of Guangdong University of Foreign Studies (TD1202).

#### References


Reprinted from *Entropy*. Cite as: Chen, L.; Qu, J.; Chai, Y.; Wu, R.; Qi, G. Synchronization of a Class of Fractional-Order Chaotic Neural Networks. *Entropy* 2013, *15*, 3355–3366.

### *Article*
