DDNS Discrete Dynamics in Nature and Society 1607-887X 1026-0226 Hindawi Publishing Corporation 627060 10.1155/2014/627060 627060 Research Article Pinning Synchronization of One-Sided Lipschitz Complex Networks http://orcid.org/0000-0002-0430-4381 Liu Fang 1 http://orcid.org/0000-0001-9635-1355 Song Qiang 1 http://orcid.org/0000-0003-3133-7119 Cao Jinde 2, 3 http://orcid.org/0000-0003-4423-6034 Lu Jianquan 2 Zuo Zhiqiang 1 School of Information Engineering Huanghuai University Henan 463000 China huanghuai.edu.cn 2 Department of Mathematics Southeast University Nanjing 210096 China seu.edu.cn 3 Department of Mathematics Faculty of Science King Abdulaziz University Jeddah 21589 Saudi Arabia kau.edu.sa 2014 1342014 2014 23 01 2014 19 03 2014 13 4 2014 2014 Copyright © 2014 Fang Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper studies the pinning synchronization in complex networks with node dynamics satisfying the one-sided Lipschitz condition which is less conservative than the well-known Lipschitz condition. Based on M-matrix theory and Lyapunov functional method, some simple pinning conditions are derived for one-sided Lipschitz complex networks with full-state and partial-state coupling, respectively. A selective pinning scheme is further provided to address the selection of pinned nodes and the design of pinning feedback gains for one-sided Lipschitz complex networks with general topologies. Numerical results are given to illustrate the effectiveness of the theoretical analysis.

1. Introduction

A complex network is composed of a set of coupled dynamical systems where each system updates its state based on the local neighboring information such that the whole network might exhibit collective behaviors under some conditions. Nowadays, many natural and man-made systems in our daily life, such as biological networks, generic networks, electrical power grids, and the World Wide Web, can be described by complex networks.

The synchronization problem for complex networks has attracted much attention from various disciplines over the past decade [1, 2]. By proposing the master stability function, Pecora and Carroll  developed a unified approach to study the local synchronization in a network of dynamical systems. Wu  showed that the synchronization in a complex network with a sufficiently large coupling strength can always be reached if the interaction digraph of the network contains a directed spanning tree. Lu and Chen  presented a systematic framework to investigate the synchronization in linearly coupled complex networks with general topologies. Gómez-Gardeñes et al.  deeply studied the effects of coupling strength and network topology on the synchronizability of complex networks.

The synchronization of complex networks is usually achieved by using full-state variables of network nodes which may not always be available in practice. To resolve this difficulty, some control algorithms have been proposed to synchronize complex networks by utilizing the observed states of network nodes. Based on the state observer approach, Jiang et al.  formulated a complex network model and then derived some criteria to discuss the local synchronization in the network. Wu and Jiao  extended the work of Jiang et al.  to address the synchronization in complex networks with asymmetric coupling.

If a complex network cannot achieve synchronization by itself, one can design some appropriate controllers to force the network to synchronize onto a homogenous trajectory. However, for a network consisting of a large number of nodes, the control cost will be very high when the control actions are applied to all network nodes. Fortunately, one can adopt the pinning control strategy  to achieve synchronization in complex networks by placing control injections onto a subset of network nodes. Current studies have shown that an undirected complex network can be synchronized to some desired state by specifically or randomly pinning some nodes . However, due to asymmetric coupling, the pinning control problem of a directed network is much more difficult than that of an undirected network. In recent years, some progress has been made in the pinning synchronization of directed networks. Chen et al.  and Lu et al.  have shown that a directed network can be pinned to the isolated node by using a minimum number of controllers, even a single controller. More recently, some researchers  proposed some techniques based on M-matrix theory  to study the pinning control problem for networked systems.

It is well-known that each node in complex networks is usually described by a nonlinear dynamical system. In most existing results on pinning control of complex networks, the node dynamics is assumed to satisfy the QUAD condition [12, 13, 15, 18]. Some authors considered the pinning control of complex networks with node dynamics satisfying Lipschitz condition [11, 14] and sector-restriction condition [22, 23]. In this paper, we further consider the pinning synchronization of complex networks composed of a set of nonlinear dynamical systems satisfying the one-sided Lipschitz condition  which is less conservative than the classical Lipschitz condition.

The main contribution of this paper is threefold. First, the paper studies the pinning control problem for both full-state and partial-state coupled complex networks with one-sided Lipschitz-type node dynamics. To the best of our knowledge, the synchronization in one-sided Lipschitz complex networks has not been addressed up to date. Therefore, the pinning control results for one-sided Lipschitz complex networks in this paper fill in this gap in time. Second, by using the properties of M-matrices, some simple pinning conditions in terms of low-dimensional linear matrix inequalities (LMIs) are established for both full-state and partial-state coupled complex networks. With the derived stability criteria, the pinning control problem of a large-scale network can be reduced to the test of a linear matrix inequality whose dimension is the same as that of a single network node. Third, we discuss the selection of pinned nodes and the design of pinning feedback gains for one-sided Lipschitz complex networks with both directed and undirected topologies based on M-matrix and algebraic graph theories.

The rest of this paper is organized as follows. In Section 2, some mathematical preliminaries are provided. Section 3 formulates the pinning control problem for complex networks with one-sided Lipschitz-type node dynamics. Sections 4 and 5 present some easily verified conditions for reaching pinning control in complex networks with full-state and partial-state coupling, respectively. Section 6 proposes a selective pinning scheme to satisfy the derived pinning conditions. In Section 7, numerical results are given to validate the theoretical analysis. Finally, some concluding remarks are stated in Section 8.

2. Preliminaries

This section provides some mathematical preliminaries to derive the main results of this paper.

2.1. Notations

The notations in this paper are quite standard. Let and represent the real number and complex number sets, respectively. Denote Re ( z ) as the real part of a complex number z . Let I N be the N -dimensional identity matrix and let 1 n n be the vector of all ones. For a square matrix A n × n , let A s = ( A + A T ) / 2 and 𝒥 ( A ) = min 1 i n Re ( λ i ( A ) ) with λ i ( A ) being the i th eigenvalue. For a symmetric real matrix X , let λ min ( X ) and λ max ( X ) be its minimum and maximum eigenvalues, respectively, and write X > 0 ( X 0 ) if X is positive definite (positive semidefinite). For a vector x n , let x be its Euclidean norm. Given two matrices A and B , denote A B as their Kronecker product .

2.2. One-Sided Lipschitz Dynamical System

A nonlinear dynamical system which does not satisfy Lipschitz condition may satisfy the so-called one-sided Lipschitz condition [24, 25].

Definition 1 (see [<xref ref-type="bibr" rid="B24">24</xref>]).

A nonlinear function f ( t , x ) is said to satisfy one-sided Lipschitz condition if there exist a positive definite matrix P and a constant γ such that P ( f ( t , x ) - f ( t , z ) ) , x - z γ x - z 2 holds for any x , z n , where · , · denotes the inner product of two vectors.

Remark 2.

The one-sided Lipschitz condition can be rewritten as ( x - z ) T P ( f ( t , x ) - f ( t , z ) ) γ x - z 2 . The constant γ is called the one-sided Lipschitz constant. It is worth mentioning that the one-sided Lipschitz constant can be any real number even a nonpositive number .

Remark 3.

It is well-known that a nonlinear function f ( t , x ) is said to satisfy Lipschitz condition if there exists a positive number θ > 0 such that f ( t , x ) - f ( t , z ) θ x - z holds for all x , z n . Therefore, one can easily obtain the Lipschitz constant θ by some simple calculations. However, for a one-sided Lipschitz function, it is quite challenging to find the one-sided Lipschitz constant γ because one has to design an appropriate structure for the matrix P in advance. In Section 7, some detailed examples are given to discuss the determination of the one-sided Lipschitz constant.

Remark 4.

From Definition 1 and Remark 3, one can see that a one-sided Lipschitz function automatically satisfies Lipschitz condition, but the converse may be not true.

2.3. M-Matrix Theory

Some properties of M-matrices are important to study the pinning control of complex networks.

Lemma 5 (see [<xref ref-type="bibr" rid="B21">21</xref>]).

For a nonsingular matrix A = ( a i j ) N × N N × N with a i j 0 ( i j ) , the following statements are equivalent:

A    is an M-matrix;

all entries of A - 1 are nonnegative;

all eigenvalues of A have positive real parts; that is, Re ( λ i ( A ) ) > 0 for all i = 1 , , N ;

there exists a positive definite diagonal matrix Ξ = diag ( ξ 1 , , ξ N ) > 0 such that Ξ A + A T Ξ > 0 .

3. Model Description and Problem Formulation

Consider a complex network composed of N identical coupled nodes, in which each node is a nonlinear dynamical system including linear and nonlinear terms: (1) x ˙ i ( t ) = A x i ( t ) + f ( t , x i ( t ) ) - σ j = 1 N l i j Γ x j ( t ) + u i ( t ) , i = 1 , , N , where x i = ( x i 1 , , x i n ) T is the state variable of the i th node, A n × n , f ( t , x i ( t ) ) n is a vector-valued function, σ > 0 is the coupling strength, l i j is the ( i , j ) th entry of Laplacian matrix L N × N defined as l i j < 0 if there is a directed link from node j to node i ( i j ) and l i j = 0 , otherwise, l i i = - k = 1 , k i N l i k , Γ n × n is the inner coupling matrix, and u i is the control input for node i .

The leader node (or isolated node) for complex network (1) is given by (2) x ˙ 0 ( t ) = A x 0 ( t ) + f ( t , x 0 ( t ) ) , where x 0 = ( x 01 , , x 0 n ) T .

For complex network (1), let 𝒱 = { 1 , , N } denote the set of all network nodes, where each node can only access the information of its neighbors. For the isolated node (2) which is not affected by any other node, we call its neighbors pinned nodes and let 𝒱 pin = { i 1 , , i l } 𝒱 denote the set of pinned nodes, where 1 l < N . Consider the following pinning control algorithm for complex network (1): (3) x ˙ i ( t ) = A x i ( t ) + f ( t , x i ( t ) ) - σ j = 1 N l i j Γ x j ( t ) - σ d i Γ ( x i ( t ) - x 0 ( t ) ) , i = 1 , , N , where the pinning feedback gains are defined as follows: (4) d i > 0 if i 𝒱 pin , d i = 0 if i 𝒱 𝒱 pin .

The node dynamics of complex network (1) is assumed to satisfy the one-sided Lipschitz condition.

Assumption 6.

Suppose that a positive definite matrix P and a constant γ can be found such that the nonlinear function f in complex network (1) satisfies one-sided Lipschitz condition; that is, (5) ( x - z ) T P ( f ( t , x ) - f ( t , z ) ) γ x - z 2 , x , z n .

4. Pinning Criteria for Network with Full-State Coupling

In this section, we consider the synchronization in complex network (3) with full-state coupling and derive some simple conditions for reaching pinning synchronization in the network.

Letting Γ = I n in complex network (3), we have the following pinning-controlled network: (6) x ˙ i ( t ) = A x i ( t ) + f ( t , x i ( t ) ) - σ j = 1 N l i j x j ( t ) - σ d i ( x i ( t ) - x 0 ( t ) ) , i = 1 , , N .

Remark 7.

We see that complex network (6) is full-state coupled with all state variables of network nodes being used to achieve pinning synchronization.

Let e i ( t ) = x i ( t ) - x 0 ( t ) , i = 1 , , N . From (2) and (6), we obtain the following error system: (7) e ˙ i ( t ) = A e i + f ( t , x i ) - f ( t , x 0 ) - σ j = 1 N l i j e j - σ d i e i , i = 1 , , N .

Denote x ( t ) = ( x 1 T ( t ) , , x N T ( t ) ) T , e ( t ) = ( e 1 T ( t ) , , e N T ( t ) ) T , and F ( t , x ) = ( f T ( t , x 1 ) , , f T ( t , x N ) ) T . Rewrite (7) in the matrix form as (8) e ˙ ( t ) = ( I N A ) e + F ( t , x ) - 1 N f ( t , x 0 ) - ( σ ( L + D ) I n ) e , where D = diag ( d 1 , , d N ) .

Based on M-matrix and algebraic graph theories, Song et al. [18, 19] showed that the quantity 𝒥 ( L + D ) = min 1 i N Re ( λ i ( L + D ) ) plays a fundamental role in analyzing the pinning control of networked systems. Now, we develop some pinning criteria for complex network (6).

Theorem 8.

Suppose that there exists a positive definite matrix P such that Assumption 6 holds and the following condition is satisfied: (9) P A + A T P - 2 α P + 2 γ I n 0 , where γ is the one-sided Lipschitz constant and α is a positive constant subject to (10) α ( 0 , σ 𝒥 ( L + D ) ) , where 𝒥 ( L + D ) = min 1 i N Re ( λ i ( L + D ) ) > 0 .

Then, the pinning-controlled network (6) globally asymptotically synchronizes to leader node (2); that is, x i ( t ) x 0 ( t ) , i = 1 , , N , as t , for any initial condition.

Proof.

Let λ i be the i th eigenvalue of L + D . It is easy to verify that σ λ i - α is the i th eigenvalue of σ ( L + D ) - α I N . It follows from condition (10) that Re ( σ λ i - α ) > 0 holds for all i = 1 , , N . Then, by the definition of L and Lemma 5, we know that σ ( L + D ) - α I N is an M-matrix and there exists a positive definite diagonal matrix Ξ = diag ( ξ 1 , , ξ N ) > 0 such that (11) [ Ξ ( σ ( L + D ) - α I N ) ] s > 0 .

Construct the following Lyapunov function candidate: (12) V ( t ) = 1 2 e T ( t ) ( Ξ P ) e ( t ) , where P satisfies LMI condition (9).

By Assumption 6 and the properties of Kronecker product , the time derivative of V ( t ) along the trajectory of error system (8) yields (13) V ˙ ( t ) = i = 1 N ξ i e i T P A e i + i = 1 N ξ i e i T P ( f ( t , x i ) - f ( t , x 0 ) ) - e T [ ( σ Ξ ( L + D ) ) P ] e 1 2 i = 1 N ξ i e i T ( P A + A T P ) e i + i = 1 N ξ i γ e i T e i - e T [ ( σ Ξ ( L + D ) ) P ] e = 1 2 i = 1 N ξ i e i T ( P A + A T P - 2 α P + 2 γ I n ) e i - e T ( [ Ξ ( σ ( L + D ) - α I N ) ] s P ) e .

It follows from (9), (11), and (13) that V ˙ ( t ) < 0 holds for all e 0 N n . Hence, the error system (8) is globally asymptotically stable at the origin. Then, complex network (6) globally asymptotically synchronizes to leader node (2).

Remark 9.

From condition (10), we see that the quantity 𝒥 ( L + D ) = min 1 i N Re ( λ i ( L + D ) ) is important to study the pinning synchronization in network (6). To ensure 𝒥 ( L + D ) > 0 , Song et al. [18, 19] have shown that the leader node should have a directed path to every other network node. Section 6 will provide more details to discuss the selection of pinned nodes and the design of pinning feedback gains.

Remark 10.

Note that the dimension of LMI (9) is the same as that of a single network node. By choosing the constant α according to condition (10), one can solve LMI (9) to find feasible solution of P by using some toolboxes such as YALMIP .

In Theorem 8, the interaction digraph 𝒢 of complex network is assumed to be general. When 𝒢 is undirected, it is easy to obtain L = L T and one has the following result from Theorem 8.

Corollary 11.

Assume that L is symmetric. There exists a positive definite matrix P such that Assumption 6 holds and (14) P A + A T P - 2 σ λ 1 P + 2 γ I n 0 , where λ 1 = λ min ( L + D ) > 0 . Then, complex network (6) globally asymptotically synchronizes to leader node (2).

Proof.

Obviously, L + D is symmetric with all eigenvalues being real numbers when L = L T . By Lemma 5, we know that matrix L + D is positive definite if λ min ( L + D ) > 0 .

Take the Lyapunov function candidate as follows: (15) V ( t ) = 1 2 e T ( t ) ( I N P ) e ( t ) , where P satisfies LMI condition (14).

In view of the proof of Theorem 8, the time derivative of V ( t ) along the trajectory of error system (8) satisfies (16) V ˙ ( t ) 1 2 i = 1 N e i T ( P A + A T P + 2 γ I n ) e i - e T [ ( σ ( L + D ) ) P ] e 1 2 i = 1 N e i T ( P A + A T P - 2 σ λ 1 P + 2 γ I n ) e i , where the second inequality is obtained by using L + D λ 1 I N .

From (14) and (16), we can show that the pinning-controlled complex network (6) globally asymptotically synchronizes to leader node (2) by using Lyapunov stability theory.

Remark 12.

When L = L T , the topology of network (3) is undirected. Then, the pinned nodes can be randomly or specifically chosen, and one can first pin the most highly connected nodes to achieve larger λ min ( L + D ) [9, 10].

5. Pinning Criteria for Network with Partial-State Coupling

In the previous section, complex network (6) is a full-state coupled network. However, in many practical cases, the full states of network nodes may not always be available. By utilizing the output states of network nodes, this section derives some pinning criteria for partial-state coupled complex networks.

For complex network (3), assume that a matrix C m × n ( m < n ) can be freely chosen such that the matrix pair ( A , C ) is detectable. For complex network (3), let Γ = P - 1 C T C , where P n × n is a positive definite matrix to be designed. Consider the following pinning-controlled network: (17) x ˙ i ( t ) = A x i ( t ) + f ( t , x i ) - σ j = 1 N l i j P - 1 C T C x j ( t ) - σ d i P - 1 C T C ( x i ( t ) - x 0 ( t ) ) , i = 1 , , N .

Remark 13.

Let y i ( t ) = C x i ( t ) , i = 0 , , N , be the output states of network nodes and let F = P - 1 C T n × m be the output-feedback gain matrix such that A - F C is Hurwitz. Then, complex network (17) can be rewritten as (18) x ˙ i ( t ) = A x i ( t ) + f ( t , x i ) - σ j = 1 N l i j F y j ( t ) - σ d i F ( y i ( t ) - y 0 ( t ) ) , i = 1 , , N .

Note that the output states are actually adopted to reach pinning synchronization in complex network (17), which means that network (17) is partial-state coupled rather than full-state coupled.

Remark 14.

If the matrix pair ( A , C ) is detectable, one can always find a matrix F n × m to ensure that A - F C is Hurwitz. Moreover, there always exist a positive definite matrix P and a scalar τ > 0 such that P A + A T P - τ C T C < 0 .

Theorem 15.

Suppose that there exists a positive definite matrix P such that Assumption 6 and the following condition hold: (19) P A + A T P - 2 β C T C + 2 γ I n < 0 , where β is a positive constant satisfying (20) 0 < β < σ 𝒥 ( L + D ) .

Then, the pinning-controlled network (17) globally asymptotically synchronizes to leader node (2).

Proof.

Following the similar line in the proof of Theorem 8, from condition (20), we can show that σ ( L + D ) - β I N is an M-matrix and there exists a positive definite diagonal matrix Ξ = diag ( ξ 1 , , ξ N ) > 0 such that (21) [ Ξ ( σ ( L + D ) - β I N ) ] s > 0 .

From (2) and (17), we obtain the following error system: (22) e ˙ ( t ) = ( I N A ) e + F ( t , x ) - 1 N f ( t , x 0 ) - ( σ ( L + D ) ( P - 1 C T C ) ) e .

Consider the following Lyapunov function candidate: (23) V ( t ) = 1 2 e T ( t ) ( Ξ P ) e ( t ) , where P satisfies Assumption 6 and condition (19).

Calculate the time derivative of V ( t ) along the trajectory of the error system (22) as follows: (24) V ˙ ( t ) = i = 1 N ξ i e i T P A e i + i = 1 N ξ i e i T P ( f ( t , x i ) - f ( t , x 0 ) ) - e T [ ( σ Ξ ( L + D ) ) ( C T C ) ] e 1 2 i = 1 N ξ i e i T ( P A + A T P ) e i + i = 1 N ξ i γ e i T e i - e T [ ( σ Ξ ( L + D ) ) ( C T C ) ] e = 1 2 i = 1 N ξ i e i T ( P A + A T P - 2 β C T C + 2 γ I n ) e i - e T ( [ Ξ ( σ ( L + D ) - β I N ) ] s ( C T C ) ) e .

Obviously, C T C is positive semidefinite; that is, C T C 0 . Then, it follows from (19), (21), and (24) that V ˙ < 0 holds for all e 0 N n . Hence, error system (22) is globally asymptotically stable at the origin, which indicates that network (17) globally asymptotically synchronizes to leader node (2).

Remark 16.

Note that the dimension of LMI (19) is the same as that of a single network node.

Remark 17.

In Remark 14, we have pointed out that one can always find a positive definite matrix P n × n and a scalar τ > 0 such that P A + A T P - τ C T C < 0 if the matrix pair ( A , C ) is detectable. To utilize this important result related to system detectability, we intentionally introduce the parameter β in Theorem 15 as a transitional variable to derive condition (19). Treating β as a scalar matrix, we can solve LMIs (19) and (20) to obtain appropriate parameters β and P .

Corollary 18.

Suppose that Assumption 6 holds and L is symmetric. There exists a positive definite matrix P such that (25) P A + A T P - 2 σ λ 1 C T C + 2 γ I n < 0 , where λ 1 = λ min ( L + D ) > 0 . Then, the complex network (17) globally asymptotically synchronizes to leader node (2).

Proof.

Take the following Lyapunov function candidate as (26) V ( t ) = 1 2 e T ( t ) ( I N P ) e ( t ) , where P satisfies LMI condition (25).

Then, the time derivative of V ( t ) along the trajectory of the error system (22) is obtained as follows: (27) V ˙ ( t ) 1 2 i = 1 N e i T ( P A + A T P ) e i + i = 1 N γ e i T e i - e T [ ( σ ( L + D ) ) ( C T C ) ] e 1 2 i = 1 N e i T ( P A + A T P - 2 σ λ 1 C T C + 2 γ I n ) e i .

By (25) and (27), we can show that complex network (17) globally asymptotically synchronizes to leader node (2).

6. Selective Pinning Scheme

In this section, we discuss the selection of pinned nodes and the design of pinning feedback gains. A selective pinning scheme for complex networks with general topologies is proposed to satisfy the pinning conditions in Theorems 8 and 15.

For complex network (1), let 𝒢 denote the interaction graph of the network and suppose that p is the multiplicity of the zero eigenvalue of the Laplacian matrix L . Song et al. [18, 19] showed that the digraph of the network can be partitioned into p components with each component having a directed tree. To ensure that 𝒥 ( L + D ) > 0 holds, at least one root node in each component should be pinned such that the leader node (2) has a directed path to every other network node. Moreover, Song et al.  have shown that 𝒥 ( L + D ) monotonically increases with respect to the number of pinned nodes or pinning feedback gains. However, 𝒥 ( L + D ) 𝒥 ( L ( 𝒱 𝒱 pin ) ) always holds even if the pinning feedback gains are sufficiently large. Therefore, it is better to take relatively lower pinning feedback gains. For directed complex networks, Song and Cao  pointed out that the nodes whose out-degrees are bigger than their in-degrees should be chosen as pinned candidates.

Now, we present the following selective pinning scheme for complex network (1) such that the conditions in Theorem 8 (or Theorem 15) can be satisfied for reaching pinning synchronization.

Partition the digraph 𝒢 of complex network (1) into p components where each component contains a directed tree. Initialize 𝒱 pin with p root nodes of these components.

Rearrange the remaining N - p nodes in descending order according to the differences of their out-degrees and in-degrees and choose appropriate pinning feedback gains.

Solve LMI (9) (or LMIs (19) and (20)). If no feasible solution is found, continually add more network nodes to 𝒱 pin or increase pinning feedback gains until the conditions in Theorem 8 (or Theorem 15) are satisfied.

7. Numerical Examples

In this section, some numerical results are provided to illustrate the effectiveness of our theoretical analysis. A digraph with ten nodes is shown in Figure 1, in which the element of the Laplacian matrix L is taken as l i j = - 1 if there is a directed link from node j to node i .

A digraph with ten nodes.

Note that the digraph in Figure 1 has a directed tree with node 1 being the unique root. According to the selective pinning scheme in Section 6, node 1 should be pinned. Let σ = 10 be the coupling strength and let d = 5 be the pinning feedback gain. If node 1 is pinned, we have (28) 𝒥 ( L + D ) = 0.2451 , σ 𝒥 ( L + D ) = 10 × 0.2451 = 2.451 .

Consider a complex network with one-sided Lipschitz-type node dynamics in the form of (1) where (29) A = ( 0 1 0 0 - 1 - 1.5 1 0 0 0 0 1 2 0 - 2 - 3 ) , and the nonlinear function is given by (30) f ( t , x i ) = ( 0 , 8 sin ( t ) , 0 , ( 1 + exp ( - x i 4 ) ) 1 / 3 ) T . It is not difficult to check that f ( t , x i ) is not a Lipschitz-type function. Now, we use some techniques in  to show that f ( t , x i ) satisfies one-sided Lipschitz condition. For any x i = ( x i 1 , x i 2 , x i 3 , x i 4 ) T , s = ( s 1 , s 2 , s 3 , s 4 ) T , and P = ( p 11 p 12 p 13 0 p 12 p 22 p 23 0 p 13 p 23 p 33 0 0 0 0 p 44 ) > 0 , applying the well-known mean-value theorem yields  (31) P ( f ( t , x i ) - f ( t , s ) ) , x i - s = - p 44 3 exp ( - ς ) ( 1 + exp ( - ς ) ) - 2 / 3 | x i 4 - s 4 | 2 0 , where ς ( min { x i 4 , s 4 } , max { x i 4 , s 4 } ) . Hence, f is a one-sided Lipschitz function with γ = 0 .

Case 1 (full-state coupling).

Taking α = 0.75 , we solve LMI (9) by using YALMIP toolbox and obtain (32) P = ( 1.0332 0.1460 - 0.1177 0 0.1460 0.4814 0.0942 0 - 0.1177 0.0942 1.0094 0 0 0 0 0.2483 ) .

Applying pinning control action to node 1 with algorithm (6), we depict the time variations of network states in Figure 2. Note that the network is successfully pinned to a homogenous state.

State evolutions of complex network with full-state coupling.

Case 2 (partial-state coupling).

Choosing C = ( 1 0 0 0 ) , it is easy to verify that the matrix pair ( A , C ) is detectable. The following parameters are obtained to satisfy conditions (19) and (20) of Theorem 15:

(33) β = 0.8669 , P = ( 0.9320 0.0464 - 0.4827 0 0.0464 0.5944 - 0.1730 0 - 0.4827 - 0.1730 0.8855 0 0 0 0 0.2972 ) T .

Using algorithm (17), we apply pinning control to the network by utilizing the output states of the network. The time variations of network states are shown in Figure 3, from which it is easy to see that the network achieves pinning synchronization.

State evolutions of complex network with partial-state coupling.

8. Conclusions

In this paper, the pinning synchronization problem for one-sided Lipschitz complex networks has been investigated by using M-matrix and algebraic graph theories. Some simple pinning criteria in terms of low-dimensional linear matrix inequalities have been established for full-state and partial-state coupled complex networks, respectively. In particular, the output states of network nodes are utilized to implement the distributed pinning control algorithm. A selective pinning scheme has been proposed to satisfy the derived pinning conditions for one-sided Lipschitz complex networks with general topologies. In the near future, it would be of interest to study the pinning control problem for one-sided Lipschitz-type complex networks with dynamically switching topologies and time delays.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was jointly supported by the National Science Foundation of China under Grants 61273218, 61304172, 61272530, and 61175119 and the Natural Science Foundation of Henan Province of China under Grants 122102210027, 122300410220 and 12B480005.

Wu C. W. Synchronization in Complex Networks of Nonlinear Dynamical Systems 2007 Singapore World Scientific Arenas A. Díaz-Guilera A. Kurths J. Moreno Y. Zhou C. Synchronization in complex networks Physics Reports 2008 469 3 93 153 10.1016/j.physrep.2008.09.002 MR2477097 Pecora L. M. Carroll T. L. Master stability functions for synchronized coupled systems Physical Review Letters 1998 80 10 2109 2112 2-s2.0-4243241129 Wu C. W. Synchronization in networks of nonlinear dynamical systems coupled via a directed graph Nonlinearity 2005 18 3 1057 1064 10.1088/0951-7715/18/3/007 MR2134084 ZBL1089.37024 Lu W. Chen T. New approach to synchronization analysis of linearly coupled ordinary differential systems Physica D. Nonlinear Phenomena 2006 213 2 214 230 10.1016/j.physd.2005.11.009 MR2201200 ZBL1105.34031 Gómez-Gardeñes J. Moreno Y. Arenas A. Synchronizability determined by coupling strengths and topology on complex networks Physical Review E. Statistical, Nonlinear, and Soft Matter Physics 2007 75 6 11 066106 10.1103/PhysRevE.75.066106 MR2362021 Jiang G.-P. Tang W. K.-S. Chen G. A state-observer-based approach for synchronization in complex dynamical networks IEEE Transactions on Circuits and Systems. I. Regular Papers 2006 53 12 2739 2745 10.1109/TCSI.2006.883876 MR2370459 Wu J. Jiao L. Observer-based synchronization in complex dynamical networks with nonsymmetric coupling Physica A. Statistical Mechanics and Its Applications 2007 386 469 480 Wang X. F. Chen G. Pinning control of scale-free dynamical networks Physica A. Statistical Mechanics and Its Applications 2002 310 3-4 521 531 10.1016/S0378-4371(02)00772-0 MR1946327 ZBL0995.90008 Li X. Wang X. F. Chen G. Pinning a complex dynamical network to its equilibrium IEEE Transactions on Circuits and Systems. I. Regular Papers 2004 51 10 2074 2087 10.1109/TCSI.2004.835655 MR2096915 Porfiri M. di Bernardo M. Criteria for global pinning-controllability of complex networks Automatica 2008 44 12 3100 3106 10.1016/j.automatica.2008.05.006 MR2531411 ZBL1153.93329 Yu W. Chen G. J. On pinning synchronization of complex dynamical networks Automatica 2009 45 2 429 435 10.1016/j.automatica.2008.07.016 MR2527339 ZBL1158.93308 Chen T. Liu X. Lu W. Pinning complex networks by a single controller IEEE Transactions on Circuits and Systems. I. Regular Papers 2007 54 6 1317 1326 10.1109/TCSI.2007.895383 MR2370589 Lu J. Ho D. W. C. Wang Z. Pinning stabilization of linearly coupled stochastic neural networks via minimum number of controllers IEEE Transactions on Neural Networks 2009 20 10 1617 1629 2-s2.0-70350347096 10.1109/TNN.2009.2027810 Song Q. Cao J. On pinning synchronization of directed and undirected complex dynamical networks IEEE Transactions on Circuits and Systems. I. Regular Papers 2010 57 3 672 680 10.1109/TCSI.2009.2024971 MR2729966 Su H. S. Wang X. F. Pinning Control of Complex Networked Systems: Synchronization, Consensus and Flocking of Networked Systems via Pinning 2013 Springer, Shanghai Jiaotong University Press 10.1007/978-3-642-34578-4 MR3014429 Su H. S. Rong Z. Chen M. Z. Q. Wang X. F. Chen G. Wang H. Decentralized adaptive pinning control for cluster synchronization of complex dynamical networks IEEE Transactions on Cybernetics 2013 43 1 394 399 Song Q. Liu F. Cao J. Yu W. Pinning-controllability analysis of complex networks: an M-matrix approach IEEE Transactions on Circuits and Systems. I. Regular Papers 2012 59 11 2692 2701 10.1109/TCSI.2012.2190573 MR2996246 Song Q. Liu F. Cao J. Yu W. M-matrix strategies for pinning-controlled leader-following consensus in multiagent systems with nonlinear dynamics IEEE Transactions on Cybernetics 2013 43 6 1688 1697 Wen G. Duan Z. Chen G. Yu W. Consensus tracking of multi-agent systems with Lipschitz-type node dynamics and switching topologies IEEE Transactions on Circuits and Systems. I: Regular Papers 2014 61 2 499 511 Horn R. A. Johnson C. R. Topics in Matrix Analysis 1991 Cambridge, UK Cambridge University Press 10.1017/CBO9780511840371 MR1091716 Li Z. Duan Z. Chen G. Global synchronised regions of linearly coupled Lur'e systems International Journal of Control 2011 84 2 216 227 10.1080/00207179.2010.546882 MR2783574 ZBL1222.93085 Song Q. Liu F. Cao J. Lu J. Some simple criteria for pinning a Lur'e network with directed topology IET Control Theory and Applications 2014 8 2 131 138 Hu G.-D. Observers for one-sided Lipschitz non-linear systems IMA Journal of Mathematical Control and Information 2006 23 4 395 401 10.1093/imamci/dni068 MR2271168 Zhao Y. Tao J. Shi N.-Z. A note on observer design for one-sided Lipschitz nonlinear systems Systems & Control Letters 2010 59 1 66 71 10.1016/j.sysconle.2009.11.009 MR2650441 ZBL1186.93017 Zhang W. Su H.-S. Liang Y. Han Z.-Z. Non-linear observer design for one-sided Lipschitz systems: an linear matrix inequality approach IET Control Theory & Applications 2012 6 9 1297 1303 10.1049/iet-cta.2011.0386 MR2986750 Löfberg J. YALMIP: a toolbox for modeling and optimization in MATLAB Proceedings of the IEEE International Symposium on Computer Aided Control System Design September 2004 Taipei, Taiwan 284 289 2-s2.0-20344396845 Boyd S. El Ghaoui L. Feron E. Balakrishnan V. Linear Matrix Inequalities in System and Control Theory 1994 15 Philadelphia, Pa, USA Society for Industrial and Applied Mathematics (SIAM) SIAM Studies in Applied Mathematics 10.1137/1.9781611970777 MR1284712