Next Article in Journal
A Novel Image Encryption Algorithm Based on DNA Encoding and Spatiotemporal Chaos
Previous Article in Journal
Modified Legendre Wavelets Technique for Fractional Oscillation Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exponential Synchronization of Two Complex Dynamical Networks of Random Disturbance with Both Mixed Coupled and Time-Varying Delay by Pinning Control

School of Computer and Software Engineering, Shenzhen Polytechnic, Shenzhen 518055, China
Entropy 2015, 17(10), 6937-6953; https://doi.org/10.3390/e17106937
Submission received: 26 June 2015 / Revised: 19 August 2015 / Accepted: 25 August 2015 / Published: 16 October 2015

Abstract

:
In this paper, the exponentially synchronization in the mean square is investigated for two different stochastic complex networks with hybrid coupling and time-varying delay via pinning control. By utilizing the Lyapunov stability theory, stochastic analysis theory, as well as matrix analysis, the sufficient conditions are derived to guarantee the exponential synchronization for any initial values through a feedback scheme. The numerical simulation is provided to show the effectiveness of the theoretical results.

1. Introduction

Nowadays, the development of society is becoming more and more complicated and networked; thus, complex networks have emerged and aroused much attention [1]. A complex network is a structure that is made up of a large set of nodes (also called vertices, oscillators, agents) that are inter-connected to varying extents by a set of links (also called edges). Complex networks, indeed, are so ubiquitously found in nature and in the modern world that it is absolutely essential for us to have a thorough understanding of their dynamical behavior (for example, coupled biological systems, such as neural networks, and socially-interacting animal species [2]), and complex network synchronization holds particular promise for applications to many fields (for example, population dynamics [2,3], power systems [4,5] and automatic control [6]).
Synchronization is a process in which two or more dynamical systems seek to adjust a certain prescribed property of their motion to a common behavior in the limit as time tends to infinity [7]. It is a phenomenon that has been widely investigated since it was discovered by Pecora and Carroll in 1990 [8]. Many synchronization patterns have been explored (for example, complete synchronization [9,10,11], cluster synchronization [12,13], phase synchronization [14], partial synchronization [15], projective synchronization [16,17]), and synchronization can be achieved by the use of pinning control [18,19], adaptive control [20,21,22], intermittent control [23], impulsive control [24,25], fuzzy control [26], hybrid control [27] or active control [28].
In the study of synchronization, inner synchronization is a synchronous phenomenon within a network and denotes the behavior of all of the nodes among one complex network becoming synchronous. Many different cases for inner synchronization have been investigated so far [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28] On the contrary, outer synchronization [29,30] is a synchronous phenomenon between networks, that is to say it refers to the synchronization occurring between two or more coupled complex networks regardless of the occurrence of inner synchronization, and many phenomena can be explained by outer synchronization in nature, such as an infectious disease spreading between different communities and different species’ development in balance [30]. Li et al. first proposed the concept of outer synchronization and derived analytically a criterion for the synchronization of two networks that have the same topological connectivity in 2007 [30]. Additionally, outer synchronization between two coupled discrete-time networks that have the same connection topologies was derived analytically [31]. On the other hand, the problem of generalized outer synchronization between two completely different complex dynamical networks was investigated in [32]. The outer synchronization between two delay-coupled complex dynamical networks with nonidentical topological structures and a noise perturbation was regarded in [33]. However, there still exists the situation that the complex networks cannot achieve synchronization.
Stochastic perturbations and time delays are important considerations when simulating realistic complex networks, because signals traveling along real physical systems are usually randomly perturbed by the environmental elements [21]. In fact, signals transmitted between nodes of complex networks are unavoidably subject to stochastic perturbations from the environment, which may cause information contained in these signals to be lost [34]; the distance and finite speed of information transmission between pairs of nodes may cause time delays. Some results have recently appeared on the synchronization of complex networks with coupling delays [35,36,37] and stochastic perturbations [38,39,40,41].
Because of the complexity of complex networks, pinning control is a technique that applies controllers to only a small fraction of the nodes in a network, and the technique is important because it greatly reduces the number of controlled nodes for real-world complex networks (which, in most cases, is huge). In fact, pinning control can be so effective for some networks that a single pinning controller is required for synchronization; namely, for complex networks that have either a symmetric or an asymmetric coupling matrix [18]. Other pinning schemes, on the other hand, are capable of globally- and exponentially-stabilizing a network into a homogeneous state by using an optimal combination of the number of pinned nodes and the feedback control gain [19].
Motivated by the above discussions, in this paper, we study the mean square exponential synchronization of two nonidentical stochastic time-varying delayed complex networks by using the pinning control method. Some sufficient conditions are derived for mean square exponential synchronization of these networks by applying the Lyapunov–Krasovskii functional method. Our results will be applicable to many synchronization problems in different fields of science and technology.
The paper is organized as follows. In Section 2, a general model of two different stochastic complex networks with both time-varying delay dynamical nodes and time-varying delay coupling and some preliminaries are given. In Section 3, some exponential synchronization criteria for such complex dynamical networks are established. In Section 4, a numerical example for verifying the effectiveness of the theoretical results is provided. We conclude the paper in Section 5.

2. Preliminaries

2.1. Notations

Throughout this paper, R n shall denote the n-dimensional Euclidean space and R n × n the set of all n × n real matrices. The superscript T shall denote the transpose of a matrix or a vector; Tr ( · ) denotes the trace of the corresponding matrix, A S = ( A + A T ) / 2 and 1 n = ( 1 , 1 , , 1 ) T R n ; and I n denotes the n-dimensional identity matrix. For square matrices M, the notation M > 0 (respectively, < 0 ) shall mean that M is a positive-definite (respectively, negative-definite) matrix, and λ max ( A ) and λ min ( A ) shall denote the greatest and least eigenvalues of a symmetric matrix, respectively. Let ( Ω , F , { F t } t 0 , P ) be a complete probability space with a filtration { F t } t 0 that is right continuous with F 0 containing all of the P -null sets. C ( [ - τ , 0 ] ; R n ) shall denote the family of continuous functions ϕ from [ - τ , 0 ] to R n with the uniform norm | | ϕ | | 2 = sup - τ s 0 ϕ ( s ) T ϕ ( s ) and C F 0 2 ( [ - τ , 0 ] ; R n ) the family of all F 0 measurable, C ( [ - τ , 0 ] ; R n ) -valued stochastic variables ξ = { ξ ( θ ) : - τ θ 0 } , such that - τ 0 E | ξ ( s ) | 2 d s , where E stands for the correspondent expectation operator with respect to the given probability measure P .

2.2. Network Model

Consider a complex network consisting of N identical linearly- and diffusively-coupled nodes with non-delayed and time-varying-delayed linear coupling. Additionally, every node in the network is an n-dimensional dynamical unit. Then, the network model of the drive system can be denoted as:
d y j ( t ) = f ( t , y j ( t ) , y j ( t - τ ( t ) ) ) + k = 1 N a j s d Σ y s ( t ) + k = 1 N b j s d Σ y s ( t - τ c ( t ) ) d t + σ ( t , y j ( t ) , y j ( t - τ ( t ) ) , y j ( t - τ c ( t ) ) ) d w ( t ) , j = 1 , 2 , , N ,
where y j ( t ) = ( y j 1 ( t ) , y j 2 ( t ) , , y j n ( t ) ) T R n is the state vector of the i-th node of the network, f ( t , y j ( t ) , y j ( t - τ ( t ) ) ) = [ f 1 ( t , y j ( t ) , y j ( t - τ ( t ) ) ) , f 2 ( t , y j ( t ) , y j ( t - τ ( t ) ) ) , , f n ( t , y j ( t ) , y j ( t - τ ( t ) ) ) ] T is a continuous vector-valued function and Σ = diag ( ϱ 1 , ϱ 2 , , ϱ n ) is an inner coupling of the networks that satisfies ϱ k > 0 , k = 1 , 2 , , n . Here, A d = [ a j s d ] R N × N and B d = [ b j s d ] R N × N are the outer coupling matrices of the network at time t and t - τ c ( t ) , respectively, such that a j s d 0 for i k , a j j d = - k = i , k i M a j s d , b j s d 0 for i k and b j j d = - k = i , k i M b j s d . τ ( t ) is the inner time-varying delay satisfying τ τ ( t ) 0 , and τ c ( t ) is the coupling time-varying delay satisfying τ c τ c ( t ) 0 . Finally, σ ( t , y i ( t ) , y j ( t - τ ( t ) ) , y j ( t - τ c ( t ) ) ) R n × n , and w ( t ) = ( w 1 ( t ) , w 2 ( t ) , , w n ( t ) ) T R n is a bounded vector-form Wiener process, satisfying:
E w k ( t ) = 0 , E w k 2 ( t ) = 1 , E w k ( t ) w k ( s ) = 0 ( s t ) .
Compared to the drive system mentioned, the response complex network is denoted as consisting of M identical linearly- and diffusively-coupled nodes with non-delayed and time-varying-delayed linear coupling. Additionally, every node in the network is an n-dimensional dynamical unit:
d x i ( t ) = f ( t , x i ( t ) , x i ( t - τ ( t ) ) ) + s = 1 M a i k r Σ x s ( t ) + s = 1 M b i k r Σ x s ( t - τ c ( t ) ) d t + σ ( t , x i ( t ) , x i ( t - τ ( t ) ) , x i ( t - τ c ( t ) ) ) d w ( t ) , i = 1 , 2 , , M ,
where x i ( t ) = ( x i 1 ( t ) , x i 2 ( t ) , , x i n ( t ) ) T R n is the state vector of the j-th node of the response network and f ( t , x i ( t ) , x i ( t - τ ( t ) ) ) = [ f 1 ( t , x i ( t ) , x i ( t - τ ( t ) ) ) , f 2 ( t , x i ( t ) , x i ( t - τ ( t ) ) ) , , f n ( t , x i ( t ) , x i ( t - τ ( t ) ) ) ] T is a continuous vector-valued function. Here, A r = [ a i k r ] R M × M and B r = [ b i k r ] R M × M are the outer coupling matrices of the network at time t and t - τ c ( t ) , respectively, such that a i k r 0 for i k , a i i r = - k = i , k i M a i k r , b i k r 0 for i k and b i i r = - k = i , k i M b i k r . In this paper, A r and A d are assumed to be irreducible in the sense that there are no isolated nodes.
The initial conditions associated with Equation (2) are:
y j ( s ) = ξ j d ( s ) , x i ( s ) = ξ i r ( s ) , - τ ˇ s 0 , j = 1 , 2 , , N , i = 1 , 2 , , M ,
where τ ˇ = max { τ ( t ) , τ c ( t ) } , ξ j d , ξ i r C F 0 b ( [ - τ ˇ , 0 ] , R n ) with the norm | | ξ | | 2 = sup - τ ˇ s 0 ξ ( s ) T ξ ( s ) .

2.3. Some Definitions, Lemmas and Assumptions

In the following, we present some definitions and lemmas that will be required throughout this paper.
Definition 1. 
The drive network (2) and the response network (1) are said to be exponentially synchronized in the mean square if the trivial solutions of Systems (2) and (1) are such that:
i = 1 M j = 1 N E | | x i ( t , t 0 , ξ i r ) - y j ( t , t 0 , ξ j d ) | | 2 K e - κ t ,
for some K > 0 and some κ > 0 for any initial data ξ i r , ξ j d C F 0 b ( [ - τ ^ , 0 ] ; R n ) .
Definition 2. 
[21] A continuous function f ( t , x , y ) : [ 0 , + ] × R n × R n R n is said to belong to the function class QUAD, denoted by f QUAD ( P , Δ , η , ζ ) for some given matrix Σ = diag { ϱ 1 , ϱ 2 , , ϱ n } , if there exists a positive definite diagonal matrix P = diag { p 1 , p 2 , , p n } , a diagonal matrix Δ = diag { δ 1 , δ 2 , , δ n } and constants η > 0 , ζ > 0 , such that f ( · ) satisfies the condition:
( x - y ) T P ( f ( t , x , z ) - f ( t , y , w ) ) - Δ Σ ( x - y ) - η ( x - y ) T ( x - y ) + ζ ( z - w ) T ( z - w )
for all x , y , z , w R n .
Almost all of the well-known chaotic systems belong to the function class QUAD, such as the Lorenz system, the Rössler system, the Chen system, the delayed Chua’s circuit, the delayed Hopfield neural network and the logistic delayed differential system [21,23].
Lemma 1. 
Consider an n-dimensional stochastic differential equation:
d x ( t ) = f ( t , x ( t ) , x ( t - τ ) ) d t + σ ( t , x ( t ) , x ( t - τ ) ) ) d ω ( t )
Let C 2 , 1 ( R + × R n ; R + ) denote the family of all nonnegative functions V ( t , x ) on R + × R n , which are twice continuously differentiable in x and once differentiable in t. If V C 2 , 1 ( R + × R n ; R + ) , define an operator L V form R + × R n to R by:
L V ( t , x ) = V t ( t , x ) + V x ( t , x ) f ( t , x , y ) + 1 2 Tr [ σ ( t , x , y ) T V x x σ ( t , x , y ) ] ,
where V t ( t , x ) = V ( t , x ) / t , V x ( t , x ) = ( V ( t , x ) / x 1 , , V ( t , x ) / x n ) , V x x ( t , x ) = ( 2 V ( t , x ) x i x j ) n × n . If V C 2 , 1 ( R + × R n ; R + ) , then for any > t > t 0 0 ,
E V ( t , x ( t ) ) = E V ( t 0 , x ( t 0 ) ) + E t 0 t L V ( s , x ( s ) ) d s
as long as the expectations of the integrals exist [42,43].
The following assumptions will be used throughout this paper for establishing the synchronization conditions.
H1 
τ ( t ) and τ c ( t ) are bounded and continuously differentiable functions, such that 0 < τ ( t ) τ , τ ˙ ( t ) < τ ¯ , 0 < τ c ( t ) τ c and τ ˙ c ( t ) < τ ¯ c . Let τ ¯ ˇ = max { τ ¯ , τ ¯ c } .
H2 
Let σ i j ( t , e i j ( t ) , e i j ( t - τ ( t ) ) , e i j ( t - τ c ( t ) ) ) = σ ( t , x i ( t ) , x i ( t - τ ( t ) ) , x i ( t - τ c ( t ) ) ) - σ ( t , y j ( t ) , y j ( t - τ ( t ) ) , y j ( t - τ c ( t ) ) ) . Then, there exist positive definite constant matrices Υ i j 1 , Υ i j 2 and Υ i j 3 for i = 1 , 2 , , M and j = 1 , 2 , , N , such that:
Tr σ i j ( t , e i j ( t ) , e i j ( t - τ ( t ) ) , e i j ( t - τ c ( t ) ) ) T P σ i j ( t , e i j ( t ) , e i j ( t - τ ( t ) ) , e i j ( t - τ c ( t ) ) ) e i j ( t ) T Υ i j 1 e i j ( t ) + e i j ( t - τ ( t ) ) T Υ i j 2 e i j ( t - τ ( t ) ) + e i j ( t - τ c ( t ) ) T Υ i j 3 e i j ( t - τ c ( t ) ) ,
where e i j ( t ) = x i ( t ) - y j ( t ) .
Remark 1. 
The assumptions are common in the literature on stochastic complex networks with time-varying delay, such as [21,23,27,38,39,40]. Assumption H1 restricts the time delay, and Assumption H2 restricts the noise strength. We shall simply write:
Υ 1 = diag Υ 11 1 , Υ 21 1 , , Υ M 1 1 , Υ 12 1 , , Υ M 2 1 , , Υ 1 N 1 , , Υ M N 1 , Υ 2 = diag Υ 11 2 , Υ 21 2 , , Υ M 1 2 , Υ 12 2 , , Υ M 2 2 , , Υ 1 N 2 , , Υ M N 2 , Υ 3 = diag Υ 11 3 , Υ 21 3 , , Υ M 1 3 , Υ 12 3 , , Υ M 2 3 , , Υ 1 N 3 , , Υ M N 3 .

3. Main Result

In the paper, the synchronization between the drive network (1) and the response network (2) will be investigated. In order to reach the synchronization, some extra controllers are required to add on the partial nodes of the response system. Without loss of generality, let the first l nodes in the response system be controlled. Additionally, the response network with a control scheme is given by:
d x i ( t ) = f ( t , x i ( t ) , x i ( t - τ ( t ) ) ) + k = 1 M a i k r Σ x k ( t ) + k = 1 M b i k r Σ x k ( t - τ c ( t ) ) + u i ( t ) d t + σ i ( t , x i ( t ) , x i ( t - τ ( t ) ) , x i ( t - τ c ( t ) ) ) d w ( t ) , i = 1 , 2 , , M ,
where u i ( t ) ( i = 1 , 2 , , M ) are the linear state feedback controllers that are defined by:
u i ( t ) = - ε i Σ ( x i ( t ) - y 1 ( t ) ) , i = 1 , 2 , , l , 0 , i = l + 1 , l + 2 , , M ,
where ε i > 0 ( i = 1 , 2 , , l ) are the control gains, denoted by Ξ = diag { ε 1 , ε 2 , , ε l , 0 , , 0 } R M × M . Define e i j ( t ) = x i ( t ) - y j ( t ) ( i = 1 , 2 , , M , j = 1 , 2 , , N ) as the synchronization error. Then, according to the controller (7), the error system is:
d e i j ( t ) = { f ( t , x i ( t ) , x i ( t - τ ( t ) ) ) - f ( t , y j ( t ) , y j ( t - τ ( t ) ) ) + k = 1 M a i k r Σ x k ( t ) - s = 1 N a j s d Σ y s ( t ) + k = 1 M b i k r Σ x k ( t - τ c ( t ) ) - s = 1 N b j s d y s ( t - τ c ( t ) ) ) + u i ( t ) } d t + σ i ( t , x i ( t ) , x i ( t - τ ( t ) ) , x i ( t - τ c ( t ) ) ) - σ j ( t , y j ( t ) , y j ( t - τ ( t ) ) , y j ( t - τ c ( t ) ) ) d w ( t ) , i = 1 , , M , j = 1 , , N .
Remark 2. 
The pinning controllers (7) are added to a part of the node in response networks by the information of the first node in drive networks, which is easy to achieve. The drive networks (1) and response networks (2) reach outer synchronization in the mean square by pinning controllers (7) as long as error System (8) is globally and exponentially stable in the mean square.
We present the outer synchronization criteria for two stochastic complex networks with time-varying delay.
Theorem 1. 
Let Assumptions H1 and H2 be true, and let f QUAD ( P , Δ , η , ζ ) . If there exist α, β, τ ( t ) and τ c ( t ) , such that:
0 τ ¯ ˇ 1 - b + c a , and a > 0
where:
a = λ min { ( η - α - β ) I M N n - I M I N ( Δ Σ ) - ( A d I M + I N A r - Ξ I N 1 ) S ( P Σ ) - 1 2 Υ 1 } , b = λ max ζ I M N n + 1 2 Υ 2 , c = λ max { ( 1 β ( B d ) T ( B d ) I M + 1 α I N ( B r ) T ( B r ) ) ( P Σ ) T ( P Σ ) + 1 2 Υ 3 } ,
then the solutions e i j ( t ) ( i = 1 , 2 , , M , j = 1 , 2 , , N ) of System (8) are globally and exponentially stable in the mean square.
Proof. 
Define the Lyapunov–Krasovskii function:
V ( t , e ( t ) ) = 1 2 i = 1 M j = 1 N e i j ( t ) T P e i j ( t ) .
By Lemma 1, we have:
L V ( t , e ( t ) , r ) = i = 1 M j = 1 N e i j ( t ) T P { f ( t , x i ( t ) , x i ( t - τ ( t ) ) ) - f ( t , y j ( t ) , y j ( t - τ ( t ) ) ) + k = 1 M a i k r Σ x k ( t ) - s = 1 N a j s d Σ y s ( t ) + k = 1 M b i k r Σ x k ( t - τ c ( t ) ) - s = 1 N b j s d y s ( t - τ c ( t ) ) ) + u i ( t ) } + 1 2 i = 1 M j = 1 N Tr { σ i j ( t , e i j ( t ) , e i j ( t - τ ( t ) ) , e i j ( t - τ c ( t ) ) ) T × P σ i j ( t , e i j ( t ) , e i j ( t - τ ( t ) ) , e i j ( t - τ c ( t ) ) ) }
For simplicity, denote e · j ( t ) = [ e 1 j T ( t ) , e 2 j T ( t ) , , e M j T ( t ) ] T R n M ; e i · ( t ) = [ e i 1 T ( t ) , e i 2 T ( t ) , , e i N T ( t ) ] T R n N ; e ( t ) = [ e · 1 T ( t ) , e · 2 T ( t ) , , e · N T ( t ) ] T R M N n .
Because f QUAD ( P , Δ , η , ζ ) , we get:
i = 1 M j = 1 N e i j ( t ) T P f ( t , x i ( t ) , x i ( t - τ ( t ) ) ) - f ( t , y j ( t ) , y j ( t - τ ( t ) ) ) - η e T ( t ) e ( t ) + ζ e T ( t - τ ( t ) ) e ( t - τ ( t ) ) + e T ( t ) I M I N ( P Δ Σ ) e ( t )
From the definitions of A r ( k = i M a i k r = 0 ), the following equation holds:
i = 1 M j = 1 N e i j T ( t ) P k = 1 M a i k r Σ x k ( t ) = i = 1 M j = 1 N e i j P k = 1 M a i k r Σ ( x k ( t ) - y j ( t ) ) = j = 1 N e · j T ( t ) A r ( P Σ ) e · j ( t ) = e T ( t ) I N A r ( P Σ ) e ( t )
Similarly,
i = 1 M j = 1 N e i j T ( t ) P k = 1 M b i k r Σ x k ( t - τ c ( t ) ) = e T ( t ) I N B r ( P Σ ) e ( t - τ c ( t ) ) , i = 1 M j = 1 N e i j T ( t ) P s = 1 N a j k d Σ x s ( t ) = e T ( t ) A d I M ( P Σ ) e ( t ) , i = 1 M j = 1 N e i j T ( t ) P s = 1 N b j k d Σ x s ( t - τ c ( t ) ) = e T ( t ) B d I M ( P Σ ) e ( t - τ c ( t ) ) .
By using x T y c x T x + 1 c y T y for any x , y R N and c > 0 , we can obtain:
e T ( t ) I N B r ( P Σ ) e ( t - τ c ( t ) ) α e T ( t ) e ( t ) + 1 α e ( t - τ c ( t ) ) I N ( B r ) T B r ( P Σ ) T ( P Σ ) e ( t - τ c ( t ) )
and:
e T ( t ) B d I M ( P Σ ) e ( t - τ c ( t ) ) β e T ( t ) e ( t ) + 1 β e T ( t - τ c ( t ) ) ( ( B d ) T B d I M ( P Σ ) T ( P Σ ) ) e ( t - τ c ( t ) )
where α > 0 and β > 0 .
Denote I N 1 = diag { 1 , 0 , , 0 } R N × N ; we have:
i = 1 l j = 1 N e i j T ( t ) P u i ( t ) = - i = 1 l j = 1 N e i j T ( t ) P ε i Σ e i 1 ( t ) = - e T ( t ) Ξ I N 1 ( P Σ ) e ( t ) .
Using Assumption H2, we get:
i = 1 M j = 1 N Tr { σ i j ( t , e i j ( t ) , e i j ( t - τ ( t ) ) , e i j ( t - τ c ( t ) ) ) T × P σ i j ( t , e i j ( t ) , e i j ( t - τ ( t ) ) , e i j ( t - τ c ( t ) ) ) } i = 1 M j = 1 N { e i j ( t ) T Υ i j 1 e i j ( t ) + e i j ( t - τ ( t ) ) T Υ i j 2 e i j ( t - τ ( t ) ) + e i j ( t - τ c ( t ) ) T Υ i j 3 e i j ( t - τ c ( t ) ) } = e ( t ) T Υ 1 e ( t ) + e ( t - τ ( t ) ) T Υ 2 e ( t - τ ( t ) ) + e ( t - τ c ( t ) ) T Υ 3 e ( t - τ c ( t ) ) .
Substituting Inequalities (12)– (18) into Equality (11), we have:
L V ( t , e ( t ) ) = e T ( t ) { ( - η + α + β ) I M N n + I M I N ( Δ Σ ) + ( A d I M + I N A r - Ξ I N 1 ) ( P Σ ) + 1 2 Υ 1 } e ( t ) + e T ( t - τ ( t ) ) ζ I M N n + 1 2 Υ 2 e ( t - τ ( t ) ) + e T ( t - τ c ( t ) ) { ( 1 β ( B d ) T ( B d ) I M + 1 α I N ( B r ) T ( B r ) )
( P Σ ) T ( P Σ ) + 1 2 Υ 3 } e ( t - τ c ( t ) )
By Equation (10), we get:
L V ( t ) - a V ( t ) + b V ( t - τ ( t ) ) + c V ( t - τ c ( t ) ) .
Define:
W ( t ) = e γ t V ( t )
where constant γ > 0 will be decided later, and use Equation (21) to compute the operator:
L W ( t ) = e γ t [ γ V ( t ) + L V ( t ) ] e γ t [ γ V ( t ) - a V ( t ) + b V ( t - τ ( t ) ) + c V ( t - τ c ( t ) ) ] ,
which, after applying the generalized Itô’s formula, gives:
e γ t E V ( t ) = e γ t 0 E V ( t 0 ) + E t 0 t L W ( s ) d s
for any t > t 0 = 0 . Hence, we have:
e γ t E V ( t ) sup - τ ˇ s 0 E V ( s ) + E 0 t e γ s { γ V ( s ) - a V ( s ) + b V ( s - τ ( s ) ) + c V ( s - τ c ( s ) ) } d s sup - τ ˇ s 0 E V ( s ) + ( γ - a ) 0 t e γ s E V ( s ) d s + b e γ τ 0 t e γ ( s - τ ( s ) ) E V ( s - τ ( s ) ) d s + c e γ τ c 0 t e γ ( s - τ c ( s ) ) E V ( s - τ c ( s ) ) d s
which, by using the change of variables s - τ ( s ) = u , gives:
0 t e γ ( s - τ ( s ) ) E V ( s - τ ( s ) ) d s = - τ ( 0 ) t - τ ( t ) e γ u E V ( u ) d u 1 - τ ˙ ( t ) τ 1 - τ ¯ sup - τ ˇ s 0 E V ( s ) + 1 1 - τ ¯ 0 t e γ u E V ( u ) d u
and a further change of variables s - τ c ( s ) = u gives:
0 t e s - τ c ( s ) E V ( s - τ c ( s ) ) d s = - τ c ( 0 ) t - τ c ( t ) e γ u E V ( u ) d u 1 - τ ˙ c ( t ) τ c 1 - τ ¯ c sup - τ ˇ s 0 E V ( s ) + 1 1 - τ ¯ c 0 t e γ u E V ( u ) d u .
Substituting Inequalities (23) and (24) into Equality (22), we have:
e γ t E V ( t ) ( 1 + b τ 1 - τ ¯ e γ τ + c τ c 1 - τ ¯ c e γ τ c ) sup - τ ˇ s 0 E V ( s ) + ( γ - a + b 1 - τ ¯ e γ τ + c 1 - τ ¯ c e γ τ c ) 0 t e γ s E V ( s ) d s
By Condition (9), this equation:
γ - a + b 1 - τ ¯ e γ τ + c 1 - τ ¯ c e γ τ c = 0 ,
has at least one positive root. Let γ be the smallest positive root of (25). Therefore, we obtain:
E V ( t ) ( 1 + b τ 1 - τ ¯ e γ τ + c τ c 1 - τ ¯ c e γ τ c ) sup - τ ˇ s 0 E V ( s ) e - γ t ,
so that:
E | | e ( t ) | | 2 2 p ^ ( 1 + b τ 1 - τ ¯ e γ τ + c τ c 1 - τ ¯ c e γ τ c ) sup - τ ˇ s 0 E V ( s ) e - γ t ,
where p ^ = min { p 1 , p 2 , , p n } . The proof is hence complete. ☐
Remark 3. 
Inequality (9) relies on the coupled matrices, the dynamics system, the noise strength and the time delays. The linear matrix inequality method is inadequate for assessing the stability of the synchronous solution.
When the time-varying delays are constant (i.e., τ ( t ) = τ , τ c ( t ) = τ c ), we obtain the following corollary:
Corollary 1. 
Let Assumptions H1 and H2 be true, and let f QUAD ( P , Δ , η , ζ ) . If there exist τ ( t ) and τ c ( t ) , such that:
b + c - a < 0 ,
where:
a = λ min ( η - α - β ) I M N n - I M I N ( Δ Σ ) - A d I M + I N A r - Ξ I N 1 S ( P Σ ) - 1 2 Υ 1 , b = λ max ζ I M N n + 1 2 Υ 2 , c = λ max ( 1 β ( B d ) T ( B d ) I M + 1 α I N ( B r ) T ( B r ) ) ( P Σ ) T ( P Σ ) + 1 2 Υ 3 ,
then the solutions e i j ( t ) ( i = 1 , 2 , , M , j = 1 , 2 , , N ) of System (8) are globally and exponentially stable in the mean square.
When σ i ( · ) = 0 , we can get the following corollary:
Corollary 2. 
Let Assumptions H1 and H2 be true, and let f QUAD ( P , Δ , η , ζ ) . If there exist τ ( t ) and τ c ( t ) , such that:
0 τ ¯ ˇ 1 - ζ + c a , and a > 0
where:
a = λ min η I M N n - I M I N ( Δ Σ ) - A d I M + I N A r - Ξ I N 1 S ( P Σ ) , c = λ max ( 1 β ( B d ) T ( B d ) I M + 1 α I N ( B r ) T ( B r ) ) ( P Σ ) T ( P Σ ) ,
then the solutions e i j ( t ) ( i = 1 , 2 , , M , j = 1 , 2 , , N ) of System Equation (8) are globally and exponentially stable in the mean square.
When A d = A r = A , B d = B r = B and σ i ( · ) = 0 , we can get the following corollary:
Corollary 3. 
Let Assumptions H1 and H2 be true, and let f QUAD ( P , Δ , η , ζ ) . If there exist τ ( t ) and τ c ( t ) , such that:
0 τ ¯ ˇ 1 - ζ + c a , and a > 0
where:
a = λ min η I M N n - I M I N ( Δ Σ ) - A I M + I N A - Ξ I N 1 S ( P Σ ) , c = λ max ( 1 β B T B I M + 1 α I N B T B ) ( P Σ ) T ( P Σ ) ,
then the solutions e i j ( t ) ( i = 1 , 2 , , M , j = 1 , 2 , , N ) of System (8) are globally and exponentially stable in the mean square.

4. Numerical Simulation

In this section, we present some numerical simulation results that validate the theorem of the previous section.
Consider the 2D chaotic delayed neural network:
s ˙ ( t ) = f ( t , s ( t ) , s ( t - τ ( t ) ) ) = - C s ( t ) + D tanh ( s ( t ) ) + E tanh ( s ( t - τ ( t ) ) ) ,
where τ ( t ) = 1 , C = 1 0 0 1 , D = 2 - 0 . 1 - 5 4 . 5 and E = - 1 . 5 - 0 . 1 - 0 . 2 - 4 .
Taking P = diag { 1 , 1 } and Δ = diag { 5 , 11.5 } , we have η = 0.15 and ζ = 3.25 , so that Condition (3) is satisfied [21]. Thus:
d x i ( t ) = f ( t , x i ( t ) , x i ( t - τ ( t ) ) ) + j = 1 7 a i k r Σ x k ( t ) + k = 1 7 b i k r Σ x k ( t - τ c ( t ) ) - ε i Σ ( x i ( t ) - y 1 ( t ) ) d t + σ i ( t , x i ( t ) , x i ( t - τ ( t ) ) , x i ( t - τ c ( t ) ) ) d w ( t ) , i = 1 , 2 , , 7 ,
and:
d y i ( t ) = f ( t , y j ( t ) , y j ( t - τ ( t ) ) ) + s = 1 5 a j s r Σ y s ( t ) + s = 1 5 b j s r Σ y s ( t - τ c ( t ) ) d t + σ j ( t , y j ( t ) , y j ( t - τ ( t ) ) , y j ( t - τ c ( t ) ) ) d w ( t ) , j = 1 , 2 , , 5 ,
where:
A r = 203 . 5999 - 4 1 1 1 1 1 - 4 1 1 1 1 1 - 3 1 0 1 1 1 - 4 1 1 1 0 1 - 3 , B r = - 4 1 1 1 1 1 - 3 1 0 1 1 1 - 3 1 0 1 0 1 - 2 0 1 1 0 0 - 2 , A d = 212 . 8970 - 6 1 1 1 1 1 1 1 - 3 1 1 0 0 0 1 1 - 5 1 1 1 0 1 1 1 - 6 1 1 1 1 0 1 1 - 4 0 1 1 0 1 1 0 - 3 0 1 0 0 1 1 0 - 3 , B d = - 3 1 1 0 1 0 0 1 - 5 1 1 1 1 0 1 1 - 4 1 0 0 1 0 1 1 - 2 0 0 0 1 1 0 0 - 4 1 1 0 1 0 0 1 - 2 0 0 0 1 0 1 0 - 2 ,
Ξ = 718.4148 diag { 1 , 1 , 0 , 0 , 0 } , τ c ( t ) = 0.01 e t 1 + e t , σ i ( t , x i ( t ) , x i ( t - τ ( t ) ) , x i ( t - τ c ( t ) ) ) = 0.01 diag { x i 1 ( t ) + x i 1 ( t - τ ( t ) ) + x i 1 ( t - τ c ( t ) ) , x i 2 ( t ) + x i 2 ( t - τ ( t ) ) + x i 2 ( t - τ c ( t ) ) } and σ i ( t , y i ( t ) , y i ( t - τ ( t ) ) , y i ( t - τ c ( t ) ) ) = 0.01 diag { y i 1 ( t ) + y i 1 ( t - τ ( t ) ) + y i 1 ( t - τ c ( t ) ) , y i 2 ( t ) + y i 2 ( t - τ ( t ) ) + y i 2 ( t - τ c ( t ) ) } .
Computations then yield τ = 1 , τ ¯ = 0 , τ c = 0.01 , τ ¯ c = 0.0025 , Υ i j = 0.01 I 2 for i = 1 , 2 , , M , and j = 1 , 2 . Then, the solutions of Equalities (10) are: a = 5.40 , b = 3.26 , c = 2.0334 ; after solving Equation (25), we obtain γ = 0.016 , so that the conditions of Theorem 1 are satisfied.
The initial conditions for this simulation are x i j ( t 0 ) = - ( 2 i + k ) cos ( t 0 ) , i = 1 , 2 , , 7 , y j s ( t 0 ) = - ( 2 ( i + 7 ) + s ) cos ( t 0 ) , j = 1 , 2 , , 5 , k = 1 , 2 and for all t 0 [ - 1 , 0 ] . The numerical approximation to the solution of stochastic differential Equation (27) and Equation (28) are found by the Euler–Maruyama method [44] with Δ t = 0.001 , and the trajectories of the pinning control gains are shown in Figure 1, Figure 2 and Figure 3. Figure 4 shows the time evolution of the synchronization errors without pinning control. One can see that the pinning control directly affects the synchronization of two complex network, which cannot be reached without pinning control.
Figure 1. The trajectories of the state variables of x i and y j ( i = 1 , 2 , , 7 and j = 1 , , 5 ) in Systems (27) and (28) under pinning control. The blue lines show the evolution of the driver systems, and the black lines represent the trajectory of the driver systems in the figure, as they do for the following figures.
Figure 1. The trajectories of the state variables of x i and y j ( i = 1 , 2 , , 7 and j = 1 , , 5 ) in Systems (27) and (28) under pinning control. The blue lines show the evolution of the driver systems, and the black lines represent the trajectory of the driver systems in the figure, as they do for the following figures.
Entropy 17 06937 g001
Figure 2. The time evolution of x i 2 and y j 2 ( i = 1 , 2 , , 7 and j = 1 , , 5 ) in Systems (27) and (28) under pinning control.
Figure 2. The time evolution of x i 2 and y j 2 ( i = 1 , 2 , , 7 and j = 1 , , 5 ) in Systems (27) and (28) under pinning control.
Entropy 17 06937 g002
Figure 3. The time evolution of x i 2 and y j 2 ( i = 1 , 2 , , 7 and j = 1 , , 5 ) in Systems (27) and (28) under pinning control.
Figure 3. The time evolution of x i 2 and y j 2 ( i = 1 , 2 , , 7 and j = 1 , , 5 ) in Systems (27) and (28) under pinning control.
Entropy 17 06937 g003
Figure 4. The time evolution of x i and y j ( i = 1 , 2 , , 5 ) in Systems (27) and (28) without pinning control.
Figure 4. The time evolution of x i and y j ( i = 1 , 2 , , 5 ) in Systems (27) and (28) without pinning control.
Entropy 17 06937 g004

5. Conclusions

In this paper, we investigated the synchronization of two different stochastic complex networks with hybrid coupling and time-varying delay. Specifically, we achieved global exponential synchronization in the mean square by applying a pinning control scheme to a small fraction of the nodes and derived sufficient conditions for the global exponential stability of synchronization in the mean square. In addition, we considered some numerical examples that have verified the effectiveness of the proposed techniques. In the future, we will consider the analysis of pinning-controlled networks with switching topologies and drive-response delays.

Acknowledgments

This work was supported by the Shenzhen Basic Research Project (JCYJ20130331152625792). The author is very grateful to the reviewers and the editor for their valuable comments and suggestions to improve the presentation of the paper.

Conflicts of Interest

The author declares that there is no conflict of interest regarding the publication of this paper.

References

  1. Girvan, M.; Newman, M.E. Community structure in social and biological networks. Proc. Natl. Acad. Sci. USA 2002, 99, 7821–7826. [Google Scholar] [CrossRef] [PubMed]
  2. Arenas, A.; Diaz-Guilera, A.; Kurths, J.; Moreno, Y.; Zhou, C. Synchronization in complex networks. Phys. Rep. 2008, 469, 93–153. [Google Scholar] [CrossRef]
  3. Dye, C.; Williams, B.G. The population dynamics and control of tuberculosis. Science 2010, 328, 856–861. [Google Scholar] [CrossRef] [PubMed]
  4. Dorfler, F.; Bullo, F. Synchronization and transient stability in power networks and nonuniform Kuramoto oscillators. SIAM J. Control Optim. 2012, 50, 1616–1642. [Google Scholar] [CrossRef]
  5. Dörfler, F.; Chertkov, M.; Bullo, F. Synchronization in complex oscillator networks and smart grids. Proc. Natl. Acad. Sci. USA 2013, 110, 2005–2010. [Google Scholar] [CrossRef] [PubMed]
  6. Zhou, J.; Zhang, H.; Xiang, L.; Wu, Q. Synchronization of coupled harmonic oscillators with local instantaneous interaction. Automatica 2012, 48, 1715–1721. [Google Scholar] [CrossRef]
  7. Boccaletti, S.; Kurths, J.; Osipov, G.; Valladares, D.L.; Zhou, C.S. The synchronization of chaotic systems. Phys. Rep. 2002, 366, 1–101. [Google Scholar] [CrossRef]
  8. Pecora, L.M.; Carroll, T.L. Synchronization in chaotic systems. Phys. Rev. Lett. 1990, 64, 821–824. [Google Scholar] [CrossRef] [PubMed]
  9. Tönjes, R.; Blasius, B. Perturbation analysis of complete synchronization in networks of phase oscillators. Phys. Rev. E 2009, 80. [Google Scholar] [CrossRef]
  10. Lu, J.; Ho, D.W.C. Globally exponential synchronization and synchronizability for general dynamical networks. IEEE Trans. Syst. Man Cybern. 2010, 40, 350–361. [Google Scholar]
  11. Yang, X.; Cao, J.; Lu, J. Synchronization of coupled neural networks with random coupling strengths and mixed probabilistic time-varying delays. Int. J. Robust Nonlinear Control 2013, 23, 2060–2081. [Google Scholar] [CrossRef]
  12. Belykh, V.N.; Belykh, I.V.; Mosekilde, E. Cluster synchronization modes in an ensemble of coupled chaotic oscillators. Phys. Rev. E 2001, 63. [Google Scholar] [CrossRef]
  13. Wang, J.; Feng, J.; Xu, C.; Zhao, Y. Cluster synchronization of nonlinearly-coupled complex networks with nonidentical nodes and asymmetrical coupling matrix. Nonlinear Dyn. 2012, 67, 1635–1646. [Google Scholar] [CrossRef]
  14. Rosenblum, M.G.; Pikovsky, A.S.; Kurths, J. Phase Synchronization of Chaotic Oscillators. Phys. Rev. Lett. 1996, 76, 1804–1807. [Google Scholar] [CrossRef] [PubMed]
  15. Van Vreeswijk, C. Partial synchronization in populations of pulse-coupled oscillators. Phys. Rev. 1996, E54, 5522–5537. [Google Scholar] [CrossRef]
  16. Zhou, P.; Bai, R.; Zheng, J. Projective synchronization for a class of fractional-order chaotic systems with fractional-order in the (1, 2) Interval. Entropy 2015, 17, 1123–1134. [Google Scholar] [CrossRef]
  17. Zhou, X.; Xiong, L.; Cai, X. Adaptive switched generalized function projective synchronization between two hyperchaotic systems with unknown parameters. Entropy 2014, 16, 377–388. [Google Scholar] [CrossRef]
  18. Chen, T.; Liu, X.; Lu, W. Pinning complex networks by a single controller. IEEE Trans. Circuits Syst. I Regul. Pap. 2007, 54, 1317–1326. [Google Scholar] [CrossRef]
  19. Zhao, J.; Lu, J.; Wu, X. Pinning control of general complex dynamical networks with optimization. Sci. China Inf. Sci. 2010, 53, 813–822. [Google Scholar] [CrossRef]
  20. Lu, J.; Cao, J.; Ho, D.W.C. Adaptive stabilization and synchronization for chaotic Lur’e systems with time-varying delay. IEEE Trans. Circuits Syst. I Regul. Pap. 2008, 55, 1347–1356. [Google Scholar]
  21. Yang, X.; Cao, J. Adaptive pinning synchronization of complex networks with stochastic perturbations. Discrete Dyn. Nat. Soc. 2010, 2010. [Google Scholar] [CrossRef]
  22. Ma, W.; Li, C.; Wu, Y. Adaptive synchronization of fractional neural networks with unknown parameters and time delays. Entropy 2014, 16, 6286–6299. [Google Scholar] [CrossRef]
  23. Wang, J.; Feng, J.; Xu, C.; Zhao, Y. Exponential synchronization of stochastic perturbed complex networks with time-varying delays via periodically intermittent pinning. Commun. Nonlin. Sci. Numer. Simul. 2013, 18, 3146–3157. [Google Scholar] [CrossRef]
  24. Haeri, M.; Dehghani, M. Impulsive synchronization of Chen’s hyperchaotic system. Phys. Lett. A 2006, 356, 226–230. [Google Scholar] [CrossRef]
  25. Dai, A.; Zhou, W.; Feng, J.; Fang, J.; Xu, S. Exponential synchronization of the coupling delayed switching complex dynamical networks via impulsive control. Adv. Differ. Equ. 2013, 2013. [Google Scholar] [CrossRef]
  26. Ahmadlou, M.; Adeli, H. Fuzzy synchronization likelihood with application to attention-deficit/hyperactivity disorder. Clin. EEG Neurosci. 2011, 42, 6–13. [Google Scholar] [PubMed]
  27. Yang, X.; Cao, J.; Lu, J. Stochastic synchronization of complex networks with nonidentical nodes via hybrid adaptive and impulsive control. IEEE Trans. Circuits Syst. I Regul. Pap. 2012, 59, 371–384. [Google Scholar] [CrossRef]
  28. Razminia, A.; Majd, V.J.; Baleanu, D. Chaotic incommensurate fractional order Rössler system: Active control and synchronization. Adv. Differ. Equ. 2011, 2011. [Google Scholar] [CrossRef]
  29. Wu, Y.; Liu, L. Exponential outer synchronization between two uncertain time-varying complex networks with nonlinear coupling. Entropy 2015, 17, 3097–3109. [Google Scholar] [CrossRef]
  30. Li, C.; Sun, W.; Kurths, J. Synchronization between two coupled complex networks. Phys. Rev. E 2007, 76. [Google Scholar] [CrossRef]
  31. Li, C.; Xu, C.; Sun, W.; Xu, J.; Kurths, J. Outer synchronization of coupled discrete-time networks. Chaos Interdiscip. J. Nonlin. Sci. 2009, 19. [Google Scholar] [CrossRef] [PubMed]
  32. Wu, X.; Zheng, W.X.; Zhou, J. Generalized outer synchronization between complex dynamical networks. Chaos Interdiscip. J. Nonlin. Sci. 2009, 19. [Google Scholar] [CrossRef] [PubMed]
  33. Wang, G.; Cao, J.; Lu, J. Outer synchronization between two nonidentical networks with circumstance noise. Physica A 2010, 389, 1480–1488. [Google Scholar] [CrossRef]
  34. Cao, J.; Wang, Z.; Sun, Y. Synchronization in an array of linearly stochastically coupled networks with time delays. Physica A 2007, 385, 718–728. [Google Scholar] [CrossRef]
  35. Li, C.; Chen, G. Synchronization in general complex dynamical networks with coupling delays. Physica A 2004, 343, 263–278. [Google Scholar] [CrossRef]
  36. Lu, W.; Chen, T.; Chen, G. Synchronization analysis of linearly coupled systems described by differential equations with a coupling delay. Physica D 2006, 221, 118–134. [Google Scholar] [CrossRef]
  37. He, W.; Cao, J. Exponential synchronization of hybrid coupled networks with delayed coupling. IEEE Trans. Neural Netw. 2010, 21, 571–583. [Google Scholar] [PubMed]
  38. Cao, J.; Wang, Z.; Sun, Y. Synchronization in an array of linearly stochastically coupled networks with time delays. Physica A 2007, 385, 718–728. [Google Scholar] [CrossRef]
  39. Ding, X.; Gao, Y.; Zhou, W.; Tong, D.; Su, H. Adaptive almost surely asymptotically synchronization for stochastic delayed neural networks with Markovian switching. Adv. Differ. Equ. 2013, 2013, 1–12. [Google Scholar] [CrossRef]
  40. Wang, J.; Xu, C.; Feng, J.; Kwong, M.K.; Austin, F. Mean-square exponential synchronization of markovian switching stochastic complex networks with time-varying delays by pinning control. Abstr. Appl. Anal. 2012, 2012. [Google Scholar] [CrossRef]
  41. Lu, J.; Ho, D.W.C.; Wu, L. Exponential stabilization of switched stochastic dynamical networks. Nonlinearity 2009, 22, 889–911. [Google Scholar] [CrossRef]
  42. Mao, X. Stability of stochastic differential equations with Markovian switching. Stoch. Process Their Appl. 1999, 79, 45–67. [Google Scholar] [CrossRef]
  43. Yuan, C.; Mao, X. Robust stability and controllability of stochastic differential delay equations with Markovian switching. Automatica 2004, 40, 343–354. [Google Scholar] [CrossRef]
  44. Higham, D.J. An algorithmic introduction to numerical simulation of stochastic differential equations. SIAM Rev. 2001, 43, 525–546. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Wu, X. Exponential Synchronization of Two Complex Dynamical Networks of Random Disturbance with Both Mixed Coupled and Time-Varying Delay by Pinning Control. Entropy 2015, 17, 6937-6953. https://doi.org/10.3390/e17106937

AMA Style

Wu X. Exponential Synchronization of Two Complex Dynamical Networks of Random Disturbance with Both Mixed Coupled and Time-Varying Delay by Pinning Control. Entropy. 2015; 17(10):6937-6953. https://doi.org/10.3390/e17106937

Chicago/Turabian Style

Wu, Xuefei. 2015. "Exponential Synchronization of Two Complex Dynamical Networks of Random Disturbance with Both Mixed Coupled and Time-Varying Delay by Pinning Control" Entropy 17, no. 10: 6937-6953. https://doi.org/10.3390/e17106937

APA Style

Wu, X. (2015). Exponential Synchronization of Two Complex Dynamical Networks of Random Disturbance with Both Mixed Coupled and Time-Varying Delay by Pinning Control. Entropy, 17(10), 6937-6953. https://doi.org/10.3390/e17106937

Article Metrics

Back to TopTop