Enhancement of the Quality and Robustness in Synchronization of Nonlinear Lur ’ e Dynamical Networks

In order to improve the synchronous reliability and dependability of complex dynamical networks, methods need to be proposed to enhance the quality and robustness of the synchronization scheme. The present study focuses on the robust fault detection issue within the synchronization for a class of nonlinear dynamical networks composed by identical Lur’e systems. Sufficient conditions in terms of linear matrix inequalities LMIs are established to guarantee global robust H−/H∞ synchronization of the network. Under such a synchronization scheme, the error dynamical system is globally asymptotically stable, the effect of external disturbances is suppressed, and at the same time, the network is sensitive to possible faults based on a mixed H−/H∞ performance. The fault sensitivity H− index, moreover, can be optimized via a convex optimization algorithm. The effectiveness and applicability of the analytical results are demonstrated through a network example composed by the Chua’s circuit, and it shows that the quality and robustness of synchronization has been greatly enhanced.


Introduction
In daily life, many physical systems can be characterized by various complex network models whose nodes are the elements of the network and the edges represent the interactions among them 1 .Treated as typical versions of large-scale systems, the notion of complex dynamical networks has drawn more and more attentions in recent years 2, 3 .One of the interesting and significant phenomena in complex dynamical networks is the synchronization of all dynamical nodes, which is a kind of typical collective behaviors and basic motions in nature 4-9 .Aiming at deriving global synchronization conditions, attempts have been made to consider the synchronization for a special class of networks composed of nonlinear Lur'e systems 10-12 .The main reason is that, in various fields of theory and engineering applications, vast amounts of nonlinear systems can be represented as the Lur'e type, including the Chua's circuit 13 , the Goodwin model 14 , and the swarm model 15 .Primary methods of dealing with such problems, among others, are developed under the framework of absolute stability theory 16 .
In order to improve the synchronous reliability and dependability, methods have been proposed to enhance the quality and robustness of the synchronization scheme.Due to the instability and poor performance that caused by noise or disturbances, it is reasonable to take the noise phenomenon into account during the synchronization process of complex dynamical networks 17, 18 .On the other hand, research in fault diagnosis has been gaining increasing consideration worldwide in the past decades 19-23 .One of the key issues related to fault detection is concerned with its robustness.Large amounts of the relevant jobs have been done for the linear systems in order to examine the robust fault detection RFD problem see 22, 23 and the references therein .In a recent work, we have investigated the robust fault sensitive synchronization of nonlinear Lur'e systems coupled in a masterslave fashion 24 .Similarly, in complex dynamical networks, since it is inevitable for faults to happen within each of the single node, a fault-free synchronization process cannot always be guaranteed.Even though, there is a few work concentrating upon the RFD problem of large-scale nonlinear systems, and hardly there is any previous work that brought the notion "fault" into physical aspects such as synchronization of nonlinear dynamical networks.
Based on these considerations, this present study considers the fault detection and disturbance rejection problem within robust synchronization for a class of dynamical networks.The network model is composed by identical nodes with each node being a perturbed nonlinear Lur'e system, while at the same time, subject to possible faults.The main challenge in evaluating the synchronization scheme is to distinguish failures from other disturbances, and accordingly, the H − /H ∞ paradigm is introduced 25 .For the purpose of description, the robustness objectives during synchronization are considered in virtue of the H ∞ norm, while the fault sensitivity specifications are expressed by utilizing the formulation of H − index.In this manner, the closed-loop error system is asymptotically stable with the H ∞ -norm from the disturbance input to controlled output reduced to a prescribed level, and at the same time, with the H − performance index maximized.By transforming the synchronization problem of dynamical networks into absolute stability problem of corresponding error systems as well as applying Lur'e system method in control theory 16 , sufficient conditions to the global robust H − /H ∞ synchronization within nonlinear Lur'e networks are developed in terms of sets of linear matrix inequalities LMI 26 .Furthermore, the derived high-dimensional LMI condition is simplified into three groups of lower-dimensional LMIs, which are easier to handle.It should be pointed out that no linearization technique is involved through derivation of all the synchronization criteria.
The rest of the paper is organized as follows.Section 2 proposes the model to be examined in this study, and gives the mathematical formulations of the global robust H − /H ∞ synchronization problem to be solved.In Section 3, the global robust H ∞ synchronization scheme of the networks is firstly studied, based on which the criteria on H − /H ∞ synchronization are then proposed in virtue of the LMI technique.Moreover, performance analysis of the network is also discussed in this part.The dynamical network composed by ten identical Chua's circuits is adopted as a numerical example in Section 4, and Section 5 closes the paper.

Notations and Preliminaries
The notations used in this study are fairly standard.R n×n is the set of n × n real matrices.For a matrix A, A T denotes its transpose.He is the Hermit operator with He A A A T .If A is a real symmetric negative definite matrix, it is shown by A < 0. diag • implies a diagonal or block-diagonal matrix.A ⊗ B indicates the Kronecker product of an n × m matrix A and a p × q matrix B, that is, If not explicitly stated, matrices are assumed to have compatible dimensions, and the terms replaced by * of a matrix refer to the terms in a symmetric position that do not need to be written out.

Basic Knowledge on Lur'e Systems
The basic model of nonlinear Lur'e systems subject to input noise and possible faults considered in this paper is described by where x ∈ R n is the state vector and z ∈ R m represents the measurement output vector.d 0 ∈ R p is an unknown input vector including disturbance, uninterested fault as well as some norm-bounded unstructured model uncertainty belonging to L 2 0, ∞ , while f 0 ∈ R q denotes the process, sensor, or actuator fault vector to be detected and isolated.Depending on specific situations under consideration, f 0 and d 0 can be modeled as different types of signals. where and the nonlinearity ϕ y is said to be in the sector 0, Θ 0 if it satisfies 2.4 .

Mathematical Problems in Engineering
Definition 2.1.Nonlinear system 2.2 is said to be absolutely stable with respect to the sector 0, Θ 0 , if the equilibrium point x 0 is globally asymptotically stable for the nonlinearity ϕ y satisfying 2.4 .
In order to characterize the influence of the disturbance and fault input, several definitions are introduced.Definition 2.2.Consider the following transfer function d 0 → z of system 2.2 : Then its H ∞ norm is defined as , where σ represents the maximal singular value.Definition 2.3.For system 2.2 , the transfer function from the input f 0 to output z is given as: , where σ stands for the minimum singular value and ω denotes the frequency band 0, ω .
Remark 2.4.The H − index defined has been widely adopted to measure the sensitivity of residual to fault in the frequency domain.A system is said to possess a better level of RFD, if the H ∞ norm of its transfer function from the disturbances to the performance variable is small; meanwhile, the H − index of the transfer function from fault to the output variable is large 19 .Various kinds of H − /H ∞ performance criteria have been proposed to determine the RFD issue 21 , and the performance is mostly adopted as a trade-off between robustness and sensitivity.In this study, for the sake of simplicity, we will consider the case of maximizing the fault sensitivity K rf s − with disturbance attenuation K rd s ∞ being a prescribed constant.

Dynamical Networks Composed of Lur'e Nodes
Consider a class of complex dynamical network model with each node being a general Lur'e system 2.2 shown as follows: where x i ∈ R n and z i ∈ R m are the state and measurement output of the ith node, respectively.d 0 and f 0 are defined as in system 2.2 , which are supposed to be the same with respect to each node.The inner coupling matrix Γ τ ij n×n denotes the coupling pattern between two nodes.G g ij N×N is the outer coupling matrix, standing for the coupling configuration of the network.If there is a connection between node i and node j i / j , then g ij g ji 1; otherwise, g ij g ji 0 i / j .The row sums of G are zero, that is, N j 1, j / i g ij −g ii , i . ., N satisfy the following.
ii If G 0ij ≥ 0 for all 1 ≤ i, j ≤ N, i / j, then the real parts of all eigenvalues of G 0 are less than or equal to 0 and all possible eigenvalues with zero part are 0. In fact, 0 is its eigenvalue of multiplicity 1.
Assume that the network 2.7 has no isolate clusters; namely, the network is connected.Under this circumstance, the coupling matrix G is symmetric and irreducible; hence it satisfies all the properties given in Lemma 2.5.Besides, suppose that the coupling matrix G has q distinct different eigenvalues λ 1 , . . ., λ q ; then there exists a nonsingular matrix U with U T U I N such that U T GU Λ, where Λ is in the following form: 2.9 Here, λ 1 0 is the maximum eigenvalue of multiply 1 and λ i is the eigenvalue of multiply Definition 2.6.When d 0 f 0 0, the dynamical network 2.7 is said to achieve global asymptotical synchronization if which can be an equilibrium point, a periodic orbit, or even a nonperiodic orbit.
From the properties of the internal coupling matrix G, the following condition holds: Define error signals e i x i −x s and residual signals r i z i −z s for i 1, 2, . . ., N. By subtracting 2.12 from 2.7 , one arrives at the dynamics of synchronization residual error: where derive from 2.8 that for c T i e / 0, the nonlinear functions φ l c T l e i ; x s , l 1, 2, . . ., m, satisfy the following sector restrictions: which leads to φ l c T l e i ; x s φ l c T l e i ; x s − θ 1l c T l e i ≤ 0, i 1, 2, . . ., N, l 1, 2, . . ., m, 2.15 and thus φ Ce i ; x s also belongs to the sector 0, Θ 1 .
Remark 2.7.Based on the basic knowledge of synchronization, the residual error dynamics must be asymptotically stable in order for the whole process to work.Note that the dynamics of the residual error signal r depends not only on f 0 , d 0 , and φ y but also on the states of each isolated node x i .In consequence, this study aims at ensuring the residual error dynamical system to be sensitive to possible faults in the regard of H − index, but the error dynamics also remain robustly asymptotically stable to external disturbance in the H ∞ sense.Under such circumstances, the dynamical network composed of Lur'e nodes is said to achieve global synchronization with a guaranteed H − /H ∞ performance.
Reformulating system 2.13 in virtue of the Kronecker product 28 as . . . 2.17 with Θ I N ⊗ Θ 1 ∈ R Nm > 0 and Φ Ce; X s belonging to the sector 0, Θ .accordingly, the residual error dynamical system 2.16 can be treated as an Nn-dimensional nonlinear Lur'e system, and the H − /H ∞ synchronization of the nonlinear dynamical networks 2.7 can be transformed into the performance analysis and stabilization problem of the corresponding residual error dynamics 2.16 .
For system 2.16 , the transfer function B f D denotes the transfer function f → r.Specifically speaking, the main objective of this present study is to determine under what condition the residual error dynamics 2.16 could be asymptotically stable and, at the same time, satisfy the following conditions: where γ is a prescribed positive constant, and β is a constant to be optimized.By applying the well-known Parseval theorem to the frequency-domain expressions 2.18 , where the ratios H ∞ norm and H − index are presented in Definitions 2.2 and 2.3, respectively, we arrive at the equivalent statements as follows: Accordingly, the definition of robust H − /H ∞ synchronization is derived as follows.
Definition 2.8.The dynamical networks composed of nonlinear Lur'e nodes in 2.7 are said to achieve global robust H − /H ∞ synchronization with disturbance attenuation γ and fault sensitivity β over the frequency range 0, ω where ω could be both finite and infinite , if with zero disturbance and zero fault, the synchronization residual error signal 2.16 is asymptotically stable, while with zero initial condition and given constants γ > 0, β > 0, conditions 2.19 -2.20 hold.

Main Results
The intention of this part is to investigate the fault sensitivity as well as disturbance rejection ability of the complex dynamical network 2.7 .In order to quantify these two performance indices, one borrows the concept of H − index and H ∞ -norm defined in the previous section.

Global H ∞ Synchronization of Nonlinear Lur'e Networks
In this subsection, we first consider the case that there is no fault existed in the network by extending previous results on H ∞ synchronization between two identical Lur'e systems to that of nonlinear Lur'e dynamical networks.Accordingly, the network model is described by

3.2
The disturbance rejection problem within the synchronization of nonlinear dynamical network 3.1 is summarized in the following definition.
Definition 3.1.Given constant scalar γ > 0, the dynamical network 3.1 is said to achieve global robust H ∞ synchronization, if system 3.2 is globally asymptotically stable with zero disturbance, and meanwhile, the performance index 2.19 is satisfied with zero initial conditions.
The robust H ∞ synchronization can be determined in virtue of the following criterion.
Theorem 3.2.Suppose that γ > 0 is a prescribed constant.For a given scalar α, if there exist positive-definite matrices P P T > 0, diagonal matrices Δ 1 diag δ 1 , . . ., δ m > 0, Π 1 diag π 1 , . . ., π m > 0, and Ω 1 diag ω 1 , . . ., ω m > 0, and matrices Q 1 and Q 2 such that the following LMI It is thus expected that Theorem 3.2 will be less conservative than some existing results due to the increasing freedom of these slack variables 29 .With the derived H ∞ synchronization conditions on Lur'e networks, the fault detection issue will then be examined in the next subsection.However, if the number of nodes is large, condition 3.3 would become a high-dimensional LMI, which is rather tedious to verify.To this end, both of these criteria will be further simplified to the test of three groups of lower-dimensional LMIs.

Fault Detection within Global H ∞ Synchronization
The RFD within a synchronization configuration can be treated as a multiple objective design task; that is, the design objective is not only being as sensitive as possible to faults such that early detection of faults is possible, but on the other hand, the sensitivity of possible faults is maximized, also suppressing the effect of disturbances and modeling errors on the synchronization error and subsequently on the residual, in order to prevent the synchronization process from being destroyed.Next theorem gives an LMI formulation for global robust H − /H ∞ synchronization.
Proof.See Appendix B.
Theorem 3.5.Suppose that α, β > 0 and γ > 0 are given scalars.If there exist matrices W i > 0, V i , and diagonal matrices Δ 1 > 0, Π 1 > 0, and Ω 1 > 0 such that the following conditions for i 1, 2 and q hold: where Corollary 3.6.For a constant α, let β > 0 and γ > 0 be prescribed constant scalars.The dynamical network 2.7 is said to achieve global robust H − /H ∞ synchronization with disturbance attenuation γ and fault sensitivity β, if there exist matrices W i > 0, V i T i , and diagonal matrices Δ 1 > 0, Ω 1 >, Π 1 > 0 such that the LMI conditions 3.5 -3.6 hold for i 1, 2 and q (corresponding to the largest, second largest, and smallest eigenvalues, resp.).Remark 3.7.If the number of nodes N is large, the H − /H ∞ synchronization criterion of the dynamical network would become a group of LMIs with rather high dimensions.In order to tackle this problem, the synchronization of the nN × nN-dimensional network has been disposed in a lower n-dimensional space through verifying three groups of n-dimensional LMIs in Corollary 3.6, and the derived conditions are quite convenient to use.
As an immediate consequence, we arrive at the simplified criterion for global robust H ∞ synchronization of nonlinear dynamical network 3.1 summarized as in the following corollary.
Corollary 3.8.For a constant α, let β > 0 and γ > 0 be prescribed constant scalars.If there exist matrices W i > 0, V i T i , and diagonal matrices Δ 1 > 0, Ω 1 >, Π 1 > 0 such that the LMIs 3.5 for i 1, 2 and q are feasible, then the dynamical network 2.7 is said to achieve global robust H ∞ synchronization.

H − /H ∞ Performance Analysis
It comes from Corollary 3.6 that the H − /H ∞ synchronization within a dynamical network can be cast into that of three sets of independent systems whose dimensions are the same as that of each isolate node.Namely, if the following systems ėλi r λi He λi Df 0 , 3.7 satisfy 3.5 -3.6 for i 1, 2 and q, then the conditions given in Definition 2.8 will be guaranteed.Suppose the transfer function of system 3.7 from d 0 → r λi and f 0 → r λi for i 1, 2, . . ., N as K rd i and K rf i , respectively.Then denote where K rdλ and K rfλ are in the following form:

3.9
On the other hand, consider the following system: ėλ where e λ e T λ1 , . . ., e T λN T and r λ r T λ1 , . . ., r T λN T .It can be found that the transfer functions from d → r λ and f → r λ of system 3.10 are just those defined in 3.8 .Moreover, by carrying out unitary transformation, K rdλ is similar to K rd , and so do K rfλ and K rf .Recall the definition of the H ∞ norm and H − index previously stated in Definitions 2.2 and 2.3; then we arrive at the following relationships between the H ∞ norms of K rd and K rdi as well as the H − indexes of K rf and K rfi for i 1, 2, . . ., N:

3.11
Conditions 3.11 show that the H ∞ norm of the transfer function from d → r in 2.16 equals to the maximum of those of the N systems 3.7 , whilst the corresponding H − index is the minimum value within those of 3.7 .Accordingly, the RFD of the network 2.7 can be cast into those of 3.7 ; thus we have the following corollary.Corollary 3.9.For a given scalar γ > 0, the performance indexes of the dynamical network 2.7 satisfy K rd ∞ < γ, and K rf − > β, if max i 1,...,N K rdi ∞ < γ and min i 1,...,N K rfi − > β hold in the decoupled systems 3.7 for i 1, 2, . . ., N.
The following corollary presents a method of deriving the maximum value of fault sensitivity and, at the same time, suppresses the external disturbance to a prescribed level for the global robust H − /H ∞ synchronization of network 2.7 .
Corollary 3.10.The nonlinear dynamical networks 2.7 are said to achieve global synchronization with guaranteed H ∞ performance γ and the maximum fault detection sensitivity β 0 √ ρ, where ρ is the global minimum of the following generalized eigenvalue minimization problem with respect to matrices W i > 0, V i for i {1, 2, q} as well as diagonal matrices Δ 1 > 0, Π 1 > 0, and Ω 1 > 0: as well as the LMI condition 3.5 holds.Here, Υ 11 and Υ 13 are described in Theorem 3.5 with constant scalars α and γ > 0 prescribed.

Numerical Examples
A lower-dimensional dynamical network model is concerned in this part so as to demonstrate the applicability and effectiveness of the approaches proposed in the previous sections.Throughout our numerical simulations, each node of the network is supposed to be a concrete Chua's circuit, which is frequently observed in various fields of theory and engineering applications 30 .
In the first stage, it will be shown that how the results derived in Section 3.1 can be used to guarantee the global robust H ∞ synchronization of the dynamical network 2.7 .Let us take a group of ten dimensionless state equations of Chua's oscillators, for example, where one of the node system is shown as system S a in Figure 1, a 1, 2, . . ., 10: Here, R 0 and R are linear resistors.The voltages across the capacitors C 1 and C 2 are denoted by v a1 and v a2 , i a3 is the current through the inductances L, and i ad is an external disturbance Mathematical Problems in Engineering 13 current that system S a subjects to.The nonlinear characteristic g v 1 represents the current through the nonlinear resistor N R , which is a piecewise-linear function expressed as Suppose that each node of the dynamical network developed by 2.7 is a circuit in the form of 4.1 .The possible coupling between two arbitrary Chua's circuits, as shown in Figure 1, indicates that there is a connection from S b to S a but none from S a to S b , where the element F plays the role of unicommunication.Depending on different values of the controller gain, the resistor R 1 can be adjusted.It is straightforward to reformulate system 4.1 into the Lur'e form as ẋi Ax i Bϕ Cx i 10 j 1 with parameter matrix H 1 0 0 .Choose system parameters as R C 2 1, p 1/RC 1 5.5, q 1/L 7.3, s R 0 /L 4, M 0 −1/7, M 1 2/7.In the following, R 1 0.3 Ω is taken.The network topology is assumed as star-like with ten nodes; thus G has the eigenvalues as follows: Picking α 3, and prescribing disturbance attenuation γ 0.9, we arrive at the feasible solutions given in Appendix B by solving the LMIs 3.5 , which, according to Corollary 3.8, implies that the dynamical network composed of Chua's circuits has achieved the global robust H ∞ synchronization.
Simulation results also confirm the effectiveness of the design.Figure 2 depicts the time response of synchronization error of the nominal dynamical network without disturbance signal d t , and it shows that the synchronization error converges to zero exponentially.Herein, initial values are taken arbitrarily.To observe the H ∞ performance with disturbance attenuation, assume the unknown input noise disturbance d i to be as Accordingly, the time response of the output residual error of Lur'e dynamical network with the above disturbance signals and zero initial conditions are shown in Figure 3.
In what follows, let us consider the global robust H − /H ∞ synchronization of the dynamical network 4.4 in the presence of fault signal f.For the purpose of illustration, the process fault is supposed to be a faulty current flowing in the same direction as i m3 along with the leftmost branch of each of the circuits, which will be simulated as two different types.Accordingly, it leads to ẋi Ax i Bϕ Cx i Remaining γ 0.9 and picking the fault sensitivity β 0.6, we arrive at solution of the LMI 3.5 -3.6 with α 3 presented in Appendix C, which on its turn ensures that the network 4.4 has achieved global robust H − /H ∞ synchronization in the presence of possible faults and external disturbances.
As for the corresponding simulation results, first let the process fault be a pulse of unit amplitude occurred from 5s to 10s and is zero otherwise .The generated residual signals r i t , i 1, 2, . . ., 10 are depicted in Figure 4 a , from which one observes that the effect of the disturbance input d i t on the residual error signal r i t , i 1, 2, . . ., 10 has been greatly reduced, and the residuals have rather large amplitudes so that the synchronization process remains sensitive to fault.With the same disturbance d i t , then redo the simulation with a different: and the results are plotted in Figure 4 b .By solving the generalized eigenvalue problem corresponding to the minimization problem given in Corollary 3.6, we get estimates of the maximum values of fault sensitivity as β 1m 0.7961, β 2m • • • β 9m 0.8548, and β 10m 0.9524., which also guarantees in terms of Corollary 3.6 that Lur'e dynamical networks achieve H − /H ∞ synchronization with β < β 0m , where β 0m min i 1,...,10 {β im } 0.7961.

Conclusion and Future Work
Aiming at enhancing the reliability and robustness of synchronization, the global robust H − /H ∞ synchronization scheme has been introduced into a class of nonlinear dynamical networks in the existence of possible faults and external disturbances.The criterion on synchronization was developed in virtue of the LMI technique such that each of the node systems of the network is robustly synchronized as well as sensitive to faults according to a mixed H − /H ∞ performance.Since both of the external disturbance and system fault are, respectively, considered, such synchronization scheme proposed here may be more practical than the synchronization in the previous literature.Moreover, the fault sensitivity H − index could be optimized via a convex optimization algorithm.In order to demonstrate the effectiveness and applicability of the derived results, a low-dimensional dynamical network with each node being a Chua's circuit has been adopted as an example.
As for future work, it will be interesting to study the synchronization of complex networks with different disturbance from various sources.Also, it is possible to extend the present results to stochastic complex networks 31-34 .

A. Proof of Theorem 3.2
Proof.First, we will show the global asymptotical stability of the residual error dynamics 3.2 with d 0 no fault is taken into account here , and accordingly 3.2 is represented as ė Ae BΦ Ce; X s . A.1 Under the given conditions, the performance indexes J 1 in condition 2.19 are then proven to be satisfied.Choose a Lyapunov functional candidate in the form of where P > 0 and Π I N ⊗ Π 1 with Π 1 diag π 11 , . . ., π 1m > 0 need to be determined.By calculating the time derivative of V along with the trajectory of the residual error dynamics 3.2 , it yields where Ψ Ce ϕ T Ce 1 , . . ., ϕ T Ce N .Then consider the sector restrictions that nonlinearities Φ Ce; X s and Ψ Ce satisfy, namely, for any diagonal matrices Δ 1 diag δ 11 , . . ., δ 1m > 0 and Ω 1 diag ω 11 , . . ., ω 1m > 0: The results are obtained with the assumption that each subsystem has the same diagonal matrices Π 1 , Ω 1 , and Δ 1 which does not affect the feasibility of inequality 2.18 .For the sake of simplicity, denote Φ Φ Ce; X s and Ψ Ψ Ce in the following contexts.Moreover, it is known from A.1 that there exist free-weighting matrices Q 1 and Q 2 with appropriate dimensions such that

Mathematical Problems in Engineering
Incorporating formulations A.4 -A.5 into equality A.3 derives and it follows that Ξ < 0 is guaranteed by the upper left block of LMI 3.3 ; hence the synchronization residual error dynamics 3.2 is globally asymptotically stable.
In the following, we will show that the restriction on performance index J 1 given in 2.19 is satisfied under zero initial conditions for all nonzero d ∈ L 2 0, ∞ .In this case, the error dynamics 3.2 is expressed by ė Ae BΦ Ce; X s B d d, r He.

A.8
Based on A.6 and A.8 , it is not difficult to derive where η 1 e T Φ T ėT Ψ T d T T , and Ξ 1 is described in condition 3.3 with Ξ 1 < 0. It further implies that for any d / 0, r t T r t − γ 2 d T t d t V t < 0. Under zero initial condition, the Lyapunov function V defined in A.2 satisfies V 0 0 and V t ≥ 0 for t > 0, hence and 2.19 is satisfied, which completes the proof.

B. Proof of Theorem 3.4
Proof.On the basis of Theorem 3.2, it is known that if there exist solutions to LMI 3.3 , the network achieves global synchronization and robust to input disturbances.As for the condition of fault detection within the H ∞ synchronization, namely, the synchronization process should be sensitive to possible input faults, it then comes to an extra verification of condition 2.20 under zero initial conditions for all nonzero f ∈ L 2 0, ∞ .In this situation, the error dynamics is given by ė Ae BΦ Ce; X s B f f, r He Df.

B.1
Following the same line of the proof of J 1 < 0 in Theorem 3.2, we know that if holds, then the constraint 2.20 will be satisfied where V is defined in A.2 , and further, the inequality condition B.2 is guaranteed by where η 2 e T Φ T ėT Ψ T f T T with Ξ 2 < 0 given in 3.4 .Thus the performance index J 2 > 0 is satisfied, and the proof is completed.

C. Proof of Theorem 3.5
Proof.To facilitate the design of the coupling matrix Γ, we designate Q 1 S and Q 2 αS, respectively, where α is a constant scalar; also it can be seen from 3.3 that α S S T < 0, and thus S is nonsingular.
Recall that there exists a unitary matrix U such that U T GU Λ with Λ defined in 2.9 .Pre-and postmultiplying to both sides of LMIs 3. 1 that all the matrices appearing in this LMI are diagonal except for matrices V , W. To this end, suppose that there exist matrices V i and W i such that for i 1, 2, . . ., N, the N LMIs 3.5 hold; then there must exist diagonal matrices V diag V 1 , V 2 , . . ., V N Mathematical Problems in Engineering and W diag W 1 , W 2 , . . ., W N as solutions to condition C.1 , and accordingly 3.3 holds.In a similar pattern, the feasibility of LMI 3.6 means that condition 3.4 holds.Moreover, since the coupling matrix G has q distinct different eigenvalues as 2.9 , it is evident to find that the number of LMI groups to be examined in 3.5 -3.6 can be reduced from N to q.On the other hand, it is noted that due to the convex property of LMI 26 , each of the rest q − 3 groups of LMIs for i 3, . . ., q − 1 can be written as a linear combination of the two groups of LMIs corresponding to the second-maximum λ 2 and the minimum eigenvalue λ q .In this situation, the synchronization condition only requires the feasibility of three groups LMIs 3.5 -3.6 with i 1, 2 and q; thus it completes the proof.E.1

3 . 1 and
the corresponding error dynamics in form of Kronecker product is expressed as ė Ae BΦ Ce; X s B d d, r He.
then the conditions given in Theorem 3.4 are ensured.Proof.See Appendix C.

4 . 4 and the nonlinear function ϕ Cx 1 / 2 10 cFigure 2 :
Figure 2: The time responses of synchronization errors of the nominal dynamical network.

10 bFigure 4 :
Figure 4: The residual responses r i t with d i t 0.5 sin 2t and a f t f 1 t ; b f t f 2 t , i 1, 2, . . ., 10.
3 by U diag U T ⊗ I N , U T ⊗ I m , U T ⊗ I m , U T ⊗ I m and U T I N S U ⊗ I N , W U T ⊗ I N P U ⊗ I N , and A Λ I ⊗ A Λ ⊗ ΓH.It implies from C.
N, 2.10 where • 2 means the Euclidean norm.x s ∈ R n is a solution of an isolate node given by ẋs Ax s Bϕ y s , y s Cx s , and Ω I N ⊗ Ω 1 , then the dynamical network 3.1 achieves global robust H ∞ synchronization with disturbance attenuation γ.Theorem 3.2 has provided a sufficient condition for the global robust H ∞ synchronization of nonlinear Lur'e networks by introducing slack matrices Q 1 and Q 2 into LMI 3.3 .