This paper investigates the cluster synchronization of impulsive complex networks with stochastic perturbation and time-varying delays. Besides, the nodes in the complex networks are nonidentical. By utilizing the Lyapunov stability theory, stochastic analysis theory, and linear matrix inequalities (LMI), sufficient conditions are derived to guarantee the cluster synchronization. The numerical simulation is provided to show the effectiveness of the theoretical results.

1. Introduction

Over the past few years, the study of complex networks has become an important issue. The complex networks are constructed with individual units called nodes, connected by links that exhibit complex topological properties, such as coupled biological and chemical system, neural networks, social interacting species, the Internet, and the World Wide Web [1]. Thereafter, the problem of control and synchronization in complex networks has attracted increasing attention for its wide applications in real world.

From the literature, there are three common phenomena in many evolving networks: delay effects, stochastic effects, and impulsive effects. First of all, considering finite switching speed of amplifiers and finite signal propagation time, time delay is ubiquitous in the implementation of electronic networks. There are two kinds of time delays in complex dynamical networks. The inner delay time causes chaos, such as delayed neural networks and delayed Chua’s circuit system and the outer coupling delay time also exists widely, such as in communication and traffic congestion. Therefore, time delays cannot be ignored in order to simulate realistic networks; see [25]. Second, the nodes in dynamical networks are often subject to various types of noise and uncertainty, which can have a great influence on the behavior of dynamical networks. For instance, signals transmitted between nodes of complex networks are unavoidably subject to stochastic perturbations from environment, which may cause information contained in these signals to be lost; see [6, 7]. Therefore, stochastic perturbations should be included in the complex network model. Third, impulsive exchange is a common phenomenon in many evolving networks. The state of nodes is often subject to instantaneous perturbations and experience abrupt change at certain instants which may be caused by switching phenomena, frequency change, or sudden noise. For example, the states of electronic networks and biological networks are often subject to instantaneous disturbances and experience abrupt changes at certain instants, which may be caused by switching phenomenon, frequency change, or other sudden noises, which could be expressed as impulsive effects. Thus, impulsive dynamical networks, which involve sudden changes at certain discrete times, are receiving more and more attention of researchers for various fields, that is, [810] and the references therein. Since time delays, stochastic perturbations, and impulsive effects can heavily affect the dynamical behaviors of the networks, it is necessary to investigate the time delay, stochastic perturbations, and impulsive effects on the synchronization of dynamical networks. Actually, impulsive stochastic perturbed complex networks have attracted the interest of many researchers for their various applications in information science, economic systems, automated control systems, and so forth.

Synchronization is one of the hot topics in the investigation of complex networks. Generally speaking, synchronization is the 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 either by virtue of coupling or by forcing [11]. Synchronization of complex network has observed wide potential applications in many different areas such as population dynamic, power system, chemical process simulation, automatic control, and many others [1218]. Synchronization patterns are classified as complete synchronization [19], lag synchronization [20], cluster synchronization [21], phase synchronization [22], and partial synchronization [23]. Among them, the cluster synchronization has received more and more attention. The cluster synchronization requires that the coupled oscillators split into subgroups called clusters, such that the oscillators synchronize with one another in the same cluster, but there is no synchronization among different clusters. The clusters can be classified by different function communities contained in the networks of real world, and the local dynamics between two function communities are different. For instances, in metabolic, neural, or software community networks, the individual nodes in each community can be viewed as the identical functional units, whereas the nodes in different communities are different since they have different functions [24]. There have been already some papers focused on sufficient conditions for the global stability of cluster synchronization of some related networks; see [25, 26] and the references therein. In [27], Wang et al. investigated the exponential synchronization of stochastic perturbed complex networks with time-varying delays via periodically intermittent pinning. Compared with the work in [27], in this paper, we would focus on the cluster synchronization of complex networks with stochastic perturbations, time-varying delays, and impulsive effects. It should also be mentioned that synchronization problems could be seen as an application of stability. There are many papers that are concerned with the stability of complex systems with time delays, such as [2831], which inspire the studies of synchronization problems.

As everyone knows, the real-world networks normally have a large number of nodes, and it is usually impractical to control a complex network by adding the controllers to all nodes. Pinning control, in which controllers are only applied to a small fraction of nodes, is an effective way to reduce the number of controlled nodes. Chen et al. [32] proved that the complex networks with symmetric or asymmetric coupling matrix achieved synchronization by a single pinning controller. In [33], Zhao et al. proposed a pinning scheme to solve the problems how to stabilize a network onto a homogeneous state globally exponentially and how to select an optimal combination between the number of pinned nodes and the feedback control gain.

Based on the above analysis, in this paper, we study the cluster synchronization of impulsive complex networks with time-varying delays and stochastic perturbations by adding feedback controllers and impulsive controllers on a fraction of selected nodes. To obtain our main results, we first formulate a new complex network nondelayed and time-varying delayed linear coupling and vector-form stochastic perturbations. By using the Lyapunov functional method, the stochastic stability analysis theory, and linear matrix inequality technique (LMI), some novel sufficient conditions are derived to guarantee cluster synchronization of the complex networks.

Remark 1. The main tools to derive the results in this paper are Lyapunov functional method, Itô formula, and linear matrix inequality technique (LMI). It is well known that Lyapunov functional method is one of the most useful tools to handle the stability problems. Since stochastic perturbations are considered in the model of this paper, Itô formula is employed to deal with the stochastic differential equations arising in the analysis of synchronization. Linear matrix inequality technique also plays an important role in the proof of our main theorem.

The paper is organized as follows. In Section 2, a general model of impulsive complex network with both time-varying delays and stochastic perturbations is given, and some preliminaries are stated. In Section 3, some cluster synchronization criteria for such complex dynamical network 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

For the facility of statements, we give some definitions of notations and lemmas, which would be used in the analysis of the next sections. Similar definitions can also be found in [27]. Suppose thatanddenote the-dimensional Euclidean space and the set of allreal matrices, respectively. The superscriptdenotes the transpose of a matrix or vector,denotes the trace of the corresponding matrix, anddenotes the identity matrix withdimensions. For square matrices, the notation() denotes thatis positive-definite (negative-definite) matrix. Suppose thatis a symmetric matrix. Denoteas the largest eigenvalue of, and denoteas the smallest eigenvalue of. Letbe a completed probability space with a filtrationwhich is right continuous andcontains all-null sets. Letbe an-dimensional Brownian motion defined on probability space. Denote bythe family of continuous functionfromtowith the uniform norm. Denote bythe family of allmeasurable,-valued stochastic variablessuch that, wherestands for the correspondent expectation operator with respect to the given probability measure.

Consider an impulsive complex network consisting ofnodes andcommunities () with nondelayed and time-varying delayed linear coupling and vector-form stochastic perturbations, which is described as whererepresents the state vector of theth node of the network;,is a continuous vector-form function. If nodebelongs to theth community, then we let. We denote bythe set of all nodes in theth community and, which is the subset of, is the index set of all nodes in theth community having direct connections to other communities in. And by the similar way, we can define. The functiondescribes the local dynamics of nodes in theth community, which is differentiable and capable of performing abundant dynamical behaviors. For any pair of indicesand, if, which means nodeand nodebelong to different communities, then.is an inner coupling of the networks, satisfying,.andare outer coupling matrices of the networks at timeand, respectively, which are both symmetric, satisfyingfor,, andfor,;is the inner time-varying delay satisfying, andis the coupling time-varying delay satisfying;,,is the noise intensity matrix, andis a bounded vector-form Weiner process, satisfying In this paper, we always assume thatis irreducible in the sense that there is no isolated node and, which means that corresponding impulsive effects are synchronizing.

The initial conditions associated with (1) are where,with the norm. We definethat satisfies which describes the identical local dynamics for the nodes in theth community. Obviously, the synchronization state in the same communityis uniform. Therefore, we assume that.

In order to achieve the cluster synchronization objective, the feedback controllers as well as impulsive controllers are added to part of its nodes. Whenfor and whenfor whererepresents control gain and by definingforotherwise, we can denote. Define the synchronization error as(), according to the controllers (5) and (6), the error system is derived as follows:

Remark 2. Since both state coupling and impulsive coupling are considered in (1), it is reasonable to expect that both the state-feedback controller (5) and the impulsive controller (6) will play important roles in the stabilizing process of the complex network (1). Thus, the controller designed in this paper consists of both the state-feedback controller (5) and the impulsive controller (6).

Definition 3. The complex network (4) is said to be exponentially synchronized if the trivial solution of system (9) satisfies whereand, for any initial data.

Definition 4 (see [6]). A continuous functionis said to be in thefunction class, denoted as, for given matrix, if there exist positive definite diagonal matrices, diagonal matricesand constants,, such thatsatisfies the following condition: for all.

We will simply introduce the following notations:

Definition 5 (see [34]). The average impulsive interval of the impulsive sequenceis less than, if there exist a positive integerand a positive number, such that wheredenotes the number of impulsive times of the impulsive sequencein the time interval.

The following assumptions will be used throughout this paper in establishing our synchronization condition.(H1).(H2)andare abounded and continuously differentiable functions, satisfying,,, and. Denote(H3)Denote,. There exist appropriate dimensional positive definite constant matrices,, andforsuch that

Lemma 6 (see [35] (Itô formula)). Consider an-dimensional stochastic differential equation Letdenote the family of all nonnegative functionson, which are twice continuously differentiable inand once differentiable in. If   , define an operatorfromtoby where/If; then for any, as long as the expectations of the integrals exist.

3. Main Results

We would give some sufficient conditions for the cluster synchronization of impulsive complex networks with time-varying delays and stochastic perturbations in this section.

Theorem 7. Suppose the assumptions (H1), (H2), and (H3) hold and. If there exist positive constantsand, such that whereis the smallest root of the equation then the solutionof system (4) could reach cluster synchronization under controller (5) when the average impulsive interval is less than, and.

Proof. Let,. Define a Lyapunov function When, according to Lemma 6, we have From condition (16), we obtain
Define whereis the largest root of (18). For, we calculate the operator Therefore, by the generalized Itô formula, we have that, for any, For, we have By changing variable, we have Substituting variablegives Taking (27) and (28) into (26) and by condition (18), we obtain When, we haveand Similarly, we can prove that, when, we have If, we get where The proof is then completed.

Remark 8. In Theorem 7, we can see that the speed of stabilizing process depends on the dynamical structure and stochastic perturbations of the complex network, as well as the impulsive effectsand the average impulsive interval. The smallerandare, the higher the speed of the stabilizing process will be.

When there is no time-varying delay coupling, the network (1) can be changed as

Corollary 9. Suppose that the assumptions (H1), (H2), and (H3) hold and. If there exist positive constantssuch that whereis the smallest root of the equation then the solutionof system (34) could reach cluster synchronization under the controller (5) when the average impulsive interval is less thanand.

When the time-varying delays are constant (i.e.,,), we get the following corollary.

Corollary 10. Suppose the assumptions (H1), (H2), and (H3) hold and. If there exist positive constantsand, such that whereis the smallest root of the equation then the solutionof system (4) could reach cluster synchronization under controller (5) when the average impulsive interval is less thanand.

4. Numerical Simulation

In this section, we give numerical simulation to verify the theorem given in the previous section.

Consider the following chaotic delayed neural networks: where,, Takeandand we have,. Hence the condition (9) can be satisfied for(see [6]).

In order to verify our results, we consider the following complex network: where In this simulation, we consider a network with the topological structure of Figure 1 for 20 nodes with 2 communities. The nodes numberedbelong to the first community, while the rest of the nodes belong to the second community. By the pinning scheme above, we pin the nodes,,,,,,,,, and.

Through computation, we get,,,, andforand. Let the control strength. By using the Matlab LMI toolbox, we can obtain the following solution for the conditions of Theorem 7:,,,,,, and.

The initial conditions of the numerical simulations are as follows:,,, for all, whereare chosen randomly in. By using the Euler-Maruyama method, the simulation is carried out with. The trajectories of the state variables ofandin system (41) under pinning control are shown in Figure 2. And Figure 3 exhibits the time evolutions of synchronization errors under pinning control.

5. Conclusion

In this paper, we investigated the cluster synchronization of impulsive complex networks with time-varying delay coupling and stochastic perturbations. Specifically, we achieved global exponential synchronization by applying pinning control scheme to a small fraction of nodes and derived sufficient conditions for the global exponential stability of synchronization. Finally, for clarity of exposition, a numerical example was considered to illustrate the theoretical analysis by using Matlab.


The authors thank the referees and the editor for their valuable comments on this paper. This work is supported by the Natural Science Foundation of China (Grant no. 61273220), Guangdong Education University Industry Cooperation Projects (Grant no. 2009B090300355), the Shenzhen Basic Research Project (JC201006010743A, JCYJ20120613105730482), and 2011 Foundation for Distinguished Young Talents in Higher Education of Guangdong (LYM11115).