Adaptive Pinning Synchronization of Complex Networks with Stochastic Perturbations

The adaptive pinning synchronization is investigated for complex networks with nondelayed and delayed couplings and vector-form stochastic perturbations. Two kinds of adaptive pinning controllers are designed. Based on an Lyapunov-Krasovskii functional and the stochastic stability analysis theory, several sufficient conditions are developed to guarantee the synchronization of the proposed complex networks even if partial states of the nodes are coupled. Furthermore, three examples with their numerical simulations are employed to show the effectiveness of the theoretical results.


Introduction
Recently, synchronization of all dynamical nodes in a network is one of the hot topics in the investigation of complex networks.It is well known that there are many useful synchronization phenomena in real life, such as the synchronous transfer of digital or analog signals in communication networks.Adaptive feedback control has witnessed its effectiveness in synchronizing a complex network 1-4 .By using the adaptive feedback control scheme, Chen and Zhou 1 studied synchronization of complex nondelayed networks, Cao et al. 2 investigated the complete synchronization in an array of linearly stochastically coupled identical networks with delays.In 3 , Zhou et al. considered complex dynamical networks with uncertain couplings.In 4 , Lu et al. studied the synchronization in arrays of delay-coupled neural networks.However, it is assumed that all the nodes need to be controlled in 1-4 .As we know, the real-world complex networks normally have a large number of nodes; it is usually impractical and impossible to control a complex network by adding the controllers to all nodes.To overcome this difficulty, pinning control, in which controllers are only applied to a small fraction of nodes, has been introduced in recent years 5-9 .By using adaptive pinning control method, Zhou et al. 10 studied local and global synchronizations of complex networks without delays; authors of 11, 12 considered the global synchronizations of the complex networks with nondelayed and delayed couplings.Note that, in most of existing results of complex networks' synchronization, all of the states are coupled for connected nodes.However, adaptive pinning synchronization results in which only partial states of the nodes are coupled are few.Hence, in this paper, we consider two different adaptive pinning controllers, which synchronize complex networks with partial or complete couplings of the nodes' states.
In the process of studying synchronization of complex networks, two main factors should be considered: time delays and stochastic perturbations.Time delays commonly exist in the real world and even vary according to time.In the subsystems, time delay can give rise to chaos, such as delayed neural network and delayed Chua's circuit system in Section 4 of this paper.Moreover, time-delayed couplings between subsystems cannot be ignored, such as in long-distance communication and traffic congestions, and so forth.In 8 , authors investigated synchronization of complex networks with delayed subsystems, while in 4, 12, 13 , authors investigated synchronization of complex networks with coupling delays.On the other hand, in real world, due to random uncertainties, such as stochastic forces on the physical systems and noisy measurements caused by environmental uncertainties, a stochastic behavior should be produced instead of a deterministic one.In fact, signals transmitted between subsystems of complex networks are unavoidably subject to stochastic perturbations from environment, which may cause information contained in these signals to be lost 2 .Therefore, transmitted signals may not be fully detected and received by other subsystems.This can have a great influence on the behavior of complex networks.There are some works in the field of synchronization of complex networks 2, 14-16 .Noise perturbations in 2, 14, 15 are all one-dimensional, which means that the signal transmitted by subsystems is influenced by the same noise.In 16 , Yang and Cao considered stochastic synchronization of coupled neural networks with intermittent control, in which noise perturbations are vector forms.Vector-form perturbation means that different subsystem is influenced by different noise, and hence is practical in the real world.
Based on the above analysis, in this paper, we study the synchronization of complex networks with time-varying delays and vector-form stochastic elements.Two different adaptive pinning controllers according to the different properties of inner couplings are considered.To obtain our main results, we first formulate a new complex network with nondelayed and delayed couplings and vector-form Wiener processes.Then, by virtue of an Lyapunov-Krasovskii functional and the stochastic stability analysis theory, we develop several theoretical results guaranteeing the synchronization of the new complex networks.Our adaptive pinning controllers are simple.The coupling matrices can be symmetric or asymmetric.Numerical simulations testify the effectiveness of our theoretical results.

Notations
In the sequel, if not explicitly stated, matrices are assumed to have compatible dimensions.I N denotes the N × N identity matrix.The Euclidean norm in R n is denoted as • ; accordingly, for vector x ∈ R n , x x T x, where T denotes transposition.G g ij N×N denotes a matrix of N-dimension, G λ max G T G , where λ max A respectively, λ min A means the largest respectively, smallest eigenvalue of matrix A. G s 1/2 G G T .G > 0 or G < 0 denotes that the matrix G is symmetric and positive or negative definite matrix, respectively.G l denotes the matrix of the first l row-column pairs of G. G N−l denotes the minor matrix of matrix G by removing all the first l row-column pairs of G. E • is the mathematical expectation.
The rest of this paper is organized as follows.In Section 2, new model of delayed complex networks with nondelayed and delayed couplings and vector-form stochastic Weiner processes is presented.Some necessary assumptions, definitions, and lemmas are also given in this section.Our main results and their rigorous proofs are described in Section 3. In Section 4, two directed networks and one Barabási-Albert BA network 17 with their numerical simulations are employed to show the effectiveness of our results.Finally, in Section 5, conclusions are given.

Preliminaries
Consider complex a network consisting of N identical nodes with nondelayed and delayed linear couplings and vector-form stochastic perturbations, which is described as ≥ 0 for i / j, g ii − N j 1,j / i g ij , and u ij ≥ 0 for i / j, u ii − N j 1,j / i u ij , i 1, 2, . . ., N. Matrix h i x 1 t , x 2 t , . . ., x N t ∈ R n×n describes the unknown coupling of the complex networks satisfying h i x t , x t , . . ., x t 0 zero matrix of n dimension .ω i t ω i1 t , ω i2 t , . . ., ω in t T is a bounded vector-form Weiner process.In this paper, we always assume that ω i t and ω j t are independent processes of one another for i / j, and matrix G g ij N×N is irreducible in the sense that there is no isolated node.
The initial condition of system 2.1 is given in the following form: where < ∞, and R −τ, 0 , R n denotes the family of all continuous R n -valued functions ϕ i s on −τ, 0 with the norm ϕ i sup −τ≤s≤0 ϕ T i s ϕ i s .We always assume that system 2.1 has a unique solution with respect to initial condition.
Our objective of synchronization is to control network 2.1 to the trajectory s t ∈ R n of the uncoupled system: ds t f t, s t , s t − τ t dt, 2.3 where s t s 1 t , s 2 t , . . ., s n t T can be any desired state: equilibrium point, a nontrivial periodic orbit, or even a chaotic orbit.To achieve this goal, some adaptive pinning controllers will be added to part of its nodes.Without loss of generality, rearrange the order of the nodes in the network, and let the first l nodes be controlled.Thus, the pinning controlled network can be described by where R i , i 1, 2, . . ., l, are control inputs.
For convenience of writing, in the sequel, we denote h i t, x 1 t , x 2 t , . . ., x N t as h i t, x t .Accordingly, denote h i t, s t , s t , . . ., s t as h i t, s t .Let e i t x i t − s t , e t e T 1 t , e T 2 t , . . ., e T N t T , i 1, 2, . . ., N. Subtracting 2.3 from 2.4 yields the following error dynamical system:

where F t, e i t , e i t − τ t f t, x i t , x i t − τ t − f t, s t , s t − τ t , and h i t, e t h i t, x t − h i t, s t .
Then the objective here is to find some appropriate adaptive pinning controllers R i , i 1, 2, . . ., l, such that the trivial solution of error system 2.5 is globally asymptotically stable in the mean square, that is, 0, i 1, 2, . . ., N.
Remark 2.2.The function class QUAD1 P, Δ, η, θ includes almost all the well-known chaotic systems with delays or without delays such as Lorenz system, R össler system, Chen system, delayed Chua's circuit as well as logistic delay differential system, delayed Hopfield neural networks and delayed CNNs, and so on.In fact, all those systems mentioned above satisfy the following condition 8, 18 : there exist positive k ij > 0 such that for for any x t x 1 t , x 2 t , . . ., x n t T ∈ R n , y t y 1 t , y 2 t , . . ., y n t T ∈ R n .
From condition 2.8 , we get Discrete Dynamics in Nature and Society where

Main Results
In this section, two different adaptive pinning feedback controllers corresponding to two kinds of properties of the inner coupling Γ are designed, and several sufficient conditions for synchronization of the complex networks 2.1 are derived.
For the the inner coupling matrix Γ, there exist the following two cases: a γ j ≥ 0, j 1, 2, . . ., n, and there exists at least one j 0 ∈ {j 1, 2, . . ., n} such that Obviously, case b is the special case of case a .Case b means that all the states of connected nodes are coupled, while in case a this is not necessary.However, most of existing papers only consider case b ; for example, see references 7, 12, 13, 19 .In this paper, we will consider both of the two cases.We first study the general case a .

Theorem 3.1. Suppose that the assumptions H
then the trivial solution of 2.5 is globally asymptotically stable in the mean square with the adaptive pinning controllers R i −cε i e i t , i 1, 2, . . ., l, εi v i e T i t Pe i t , i 1, 2, . . ., l.

3.1
Proof.According to H 1 and H 2 , the origin is an equilibrium point of the error system 2.5 .We define the Lyapunov-Krasovskii functional as where In view of Lemma 2.5, Π < 0 is equivalent to where G s * is matrix with appropriate dimension.Since 2cγ G s N−l W N−l < 0, obviously, there exist positive constants k 1 , k 2 , . . ., k l such that Again, from Lemma 2.5 we obtain Υ < 0. Hence, 3.8 holds, which means that Π < 0. Taking the mathematical expectation of both sides of 3.6 , we have

3.11
where λ min −Π is positive.In view of the LaSalle invariance principle of stochastic differential equation, which was developed in 21 , we have lim t → ∞ E N i 1 e i t 2 0, which in turn illustrates that lim t → ∞ E e i t 2  0, and at the same time, lim t → ∞ ε i t ε i constants .This completes the proof.
Next, we consider the case b .Obviously, the adaptive pinning controllers 3.1 can also synchronize 2.1 for this case.Here we use another kind of controllers, which is related to the inner matrix Γ, to synchronize the complex networks 2.1 .
Differentiating both sides of 3.13 along the solution of 2.
In view of Lemma 2.5, Π j < 0 is equivalent to where ΞG s * is the matrix with appropriate dimension.Since 2cp j γ j ΞG s

3.18
Again, from Lemma 2.5 we obtain Υ j < 0. Hence, 3.16 holds, which implies Π j < 0. The rest steps are the same as those in the proof of Theorem 3.1.This completes the proof.
Remark 3.3.If there is no noise perturbation in 2.1 , then Φ 0 in Theorem 3.1 and Φ 0 in Theorem 3.2, respectively.Hence, if 2.1 is synchronized with the first l nodes to be controlled by adding controllers 3.1 or 3.12 , then 2.1 with no noise perturbation can sure be synchronized by using adapting pinning controllers 3.1 or 3.12 with no more than l nodes to be controlled.Remark 3.4.Theorems 3.1 and 3.2 are applicable to directed networks and undirected networks.Although they do not point out which nodes should be controlled first; however, according to 22, 23 , for the graph corresponding to G, nodes whose out-degrees are larger than their in-degrees 24 should be controlled first for directed networks, while the controlled nodes can be randomly selected and it is better to successively pin the most-highly connected nodes for undirected networks.

Numerical Examples
In this section, we provide three examples to illustrate the effectiveness of the results obtained above.
Example 4.1.Consider the following chaotic delayed neural networks: where f s tanh s , τ 1, In the case that the initial condition is chosen as s 1 t 0.4, s 2 t 0.6, for all t ∈ −1, 0 , the chaotic attractor can be seen in Figure 1.
In order to verify our new results, consider the following coupled networks: where Γ and D are 2-dimensional identity matrices, c 4, d 2, and G and U are asymmetric and zero-row sum matrices as the following: Through simple computation, we get that the left eigenvector of G corresponding to eigenvalue 0 is 0.5494, 0.8242, 0.1374 T and ΞU 1.3901.Let h i x According to Theorem 3.2, the complex networks 4.3 can be controlled to the state of 4.1 under the adaptive pinning controllers 3.12 with the first two nodes to be controlled.
The initial conditions of the numerical simulations are as follows:  where c 1.7, d 0.
The initial conditions of the numerical simulations are as follows: 4 shows the time evolutions of synchronization errors with adaptive pinning control.The trajectories of the adaptive pinning control gains are shown in Figure 5. where s t s 1 t , s 2 t , s 3 t T , g and δ 10, ξ 19.53, ω 0.1636, a −1.4325, b −0.8731, ν 0.5, 0.2, τ 0.02.In the case that the initial condition is chosen as s t t 2; s 2 t 0.2; s 3 t 0.3, for all t ∈ −0.02, 0 , the chaotic attractor can be seen in Figure 6.Now we construct a complex network, which obeys the scale-free distribution of the Barabási-Albert model 17 .The parameters in the process of constructing are the following: initial graph is complete with m 0 5 nodes, m 5 edges are added to the network when a new node is introduced, and the final number of nodes is N 100.See Figure 7 for the BA scare-free network.
Consider the following complex networks:  where G U g ij 100×100 is the graph laplacian of the BA scale-free network, h i x 1 t , x 2 t , . . ., x 100 t diag x i,1 t − x i 1,1 t , x i,2 t − x i 1,2 t , x i,3 t − x i 1,3 t with x 101 t x 1 t , i 1, 2, . . ., 100.Without loss of generality, we take the first 10 nodes to be controlled.The initial conditions of the numerical simulations are as follows: x i s −8 i, −5 i, −10 i T , s ∈ −0.02, 0 , 1 ≤ i ≤ 100, ε j 0 1, v j 0.04.For the synchronization errors e ij t x ij t − s j t , 1 ≤ i ≤ 100, j 1, 2, 3 and trajectories of the control gains; see Figures 8 and 9.
Remark 4.4.In Example 4.2, only the first state of all nodes is coupled; however, the complex network 4.6 is synchronized by adaptive pinning of the first two nodes, which verifies the effectiveness of our theoretical results.In 7,12,13,19 , however, all states of the nodes should be coupled.Therefore, our results extend some of existing results.On the other hand, Examples 4.1 and 4.2 are directed networks, while Example 4.3 is undirected network with 100 nodes.Examples with numerical simulations show that our theoretical results are applicable to directed and undirected networks, even networks of large size.

Conclusions
In this paper, we have studied the adaptive pinning synchronization of the complex networks with delays and vector-form stochastic perturbations, in which dynamical behaviors are more realistic and complicated.Since the new condition QUAD1, our proof is very simple.Via two simple adaptive pinning feedback control schemes, several sufficient conditions guaranteeing the synchronization of the proposed complex networks are obtained.Our results are valid even if only partial states of nodes are coupled.Numerical simulations verified the effectiveness of our results.Our results improve and extend some existing results.Models and results in this paper provide possible new applications for network designers.

Figure 2 :
Figure 2: Time evolutions of synchronization errors with the first two nodes being controlled.

Figure 7 :
Figure 7: BA scale-free network initial graph is complete with m 0 5 nodes; m 5 edges are added in the network when a new node is introduced and the final number of nodes is N 100.

5 cFigure 8 :
Figure 8: Synchronization errors e i1 t a , e i2 t b , and e i3 t c 1 ≤ i ≤ 100 of 4.10 with the first 10 nodes being controlled.
and 2ab ≤ a 2 b 2 is used in the above inequality.Hence, condition 2.7 is satisfied.Note that function class QUAD1 P, Δ, η, θ is more general than the usual function class QUAD P, Δ, η see 13, 19 for chaotic systems and includes it as a special class when θ 0.Before giving our main results, we present the following lemmas, which are needed in the next section.
are positive constants, and Q and k i , i 1, 2, . .., l, are to be determined.Differentiating V 1 t along the solution of 2.5 and using the properties of vector-form Weiner process 16 , we get 1 t , e 2 t , . . ., e N t T , Q is a positive definite matrix, k i , i 1, 2, . . ., l, T h i t, e t dt N i 1 e T i t P F t, e i t , e i t − τ t − Δe i t T t Qζ t − 1 − τ t ζ T t − τ t Qζ t − τ t ≤ ζ T t Qζ t − 1 − h ζ T t − τ t Qζ t − τ t dt.where K diag k 1 , k 2 , . . ., k l , 0, . . ., 0 , T t ζ T t , ζ T t − τ t T , and 1 t , x 2 t , x 3 t 0.2 diag x i,1 t − x i 1,1 t , x i,2 t − x i 1,2 t , i 1, 2, 3, where .Hence, the assumption conditions H 1 and H 2 are satisfied.Selecting the first two nodes to be controlled, we have 2cp 1 γ 1 ΞG s 04. Figure2shows the time evolutions of synchronization errors e ij tx ij t − s j t , i 1, 2, 3, j 1, 2, with the first two nodes being controlled, which verify Theorem 3.2 perfectly.The trajectories of the gains are shown in Figure3.