Abstract
The cluster synchronization of linearly coupled complex networks with identical and nonidentical nodes is studied. Without assuming symmetry, we proved that these linearly coupled complex networks could achieve cluster synchronization under certain pinning control schemes. Sufficient conditions guaranteeing cluster synchronization for any initial values are derived by using Lyapunov function methods. Moreover, the adaptive feedback algorithms are proposed to adjust the control strength. Several numerical examples are given to illustrate our theoretical results.
1. Introduction
Recently, an increasing interest has been devoted to the study of complex networks. Among them, synchronization is the most interesting. In fact, synchronization of complex networks has been found to be a universal phenomenon in nature and has important potential applications in real-world dynamical systems. Great interests and attentions have been received for the synchronization of complex networks in many research and application fields including secure communication, seismology, parallel image processing, chemical reaction, and others [1–7].
There are many widely studied synchronization patterns, such as complete synchronization [8], lag synchronization [9], cluster synchronization [10], phase synchronization [11], and partial synchronization [12] Among them, the studies on cluster synchronization have received more and more attentions. 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, which could describe the behaviors of the complex network in the real world. For instance, the metabolic, neural, or software networks containing some different function communities. Thus, it is a natural idea to consider the cluster synchronization of such community networks.
The complex network we considered in this paper is the linearly coupled ordinary differential equations (LCODEs). In fact, LCODEs are a large class of dynamical systems with continuous time and state, as well as discrete space, which are widely used to describe coupling oscillators. Nowadays, cluster synchronization of different kinds of LCODEs has been widely studied, and many results have already exist on the various properties of such problem. For instance, Ma et al. [13] constructed a novel coupling scheme with cooperative and competitive weight couplings that guarantees the cluster synchronization of any connected networks with identical nodes. The authors also derived a sufficient condition for the global stability of cluster synchronization. In [14], Wu et al. have discussed the problem of driving linearly coupled networks to an arbitrarily selected cluster synchronization pattern via pinning control. They introduced a single negative feedback controller for each cluster to pin the coupled system to the assigned cluster synchronization pattern for any initial values. However, in most cases, couplings between nodes are not the same even if the diffusive condition is still satisfied. Nodes usually receive instantaneous information as well as delayed information from their neighbors. Thus, it is a nature idea to study the synchronization of the networks with both delayed and nondelayed coupling. In this paper, we would investigate cluster synchronization of LCODEs with both delayed and nondelayed coupling under pinning control scheme. First, we assume all the node in LCODEs are identical. By utilizing the Lyapunov stability method, the global stability of cluster synchronization in networks is investigated, and several sufficient conditions for the global stability are given. Furthermore, we propose an adaptive feedback algorithms to adjust the control strength for LCODEs.
Since in the real world, many networks contain some different function communities 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 [15]. One method to solve such problem is to consider the cluster synchronization of community networks with nonidentical nodes.There have been already some papers focused on sufficient conditions for the global stability of cluster synchronization of some related networks. A number of sufficient conditions were similarly obtained by Lu et al. [16] for the cluster synchronization of dynamical networks with community structure and nonidentical nodes in the presence of time delays by using a certain feedback control scheme. Lu et al. [17] studied the cluster synchronization of general bi-directed networks with nonidentical clusters and derived sufficient conditions for achieving local cluster synchronization of networks. The authors also discovered a relationship between the cluster synchronizability of a network and its intra-to-intercluster link ratio with the help of numerical examples. Recently, Wang et al. [18] considered the cluster synchronization of dynamical networks with community structure and nonidentical nodes and with identical local dynamics for all individual nodes in each community by using pinning control schemes. In this paper, we also investigate cluster synchronization of LCODEs with nonidentical nodes under pinning control scheme. By utilizing the Lyapunov-Krasovskii stability method, the global stability of cluster synchronization in networks is investigated, and several sufficient conditions for the global stability are given. Compared with [16–18], the complex network model we considered in this paper is more general. And moreover, the coupling matrices with and without time delay are asymmetric.
The paper is organized as follows. In Section 2, some necessary and useful definitions and lemmas are given. In Section 3, we study the global cluster synchronization of LCODEs with identical nodes and give a sufficient condition for it. Then, the adaptive feedback algorithms on control strength are proposed to achieve cluster synchronization in the complex network. In Section 4 cluster synchronization of LCODEs with nonidentical nodes is investigated and sufficient conditions are derived to achieve cluster synchronization. And the adaptive feedback algorithms on control strength are also proposed. In Section 5, numerical simulation are presented. We conclude the paper in Section 6.
2. Preliminaries
First, we introduce the mathematical definition of cluster synchronization.
Definition 2.1 (see [14]). Let be a partition of the set into nonempty subsets, that is, and . For , let denote the subscript of the subset in which the number is, that is, . A network with identical oscillators is said to realize -cluster synchronization with the partition if, for any initial values, the state variables of the oscillators satisfy for and for .
For convenience of the statement to our main results, we now make some definitions for a class of functions and a class of matrices.
Definition 2.2 (see [19–21]). Suppose that is a class of continuous functions . Let be a positive definite diagonal matrix and be a diagonal matrix. if and only if holds for some , where and .
Definition 2.3 (see [14]). For matrix
where , . If each block is a zero-row-sum matrix, then we say that .
For an asymmetric matrix with zero-row-sums, we have the following.
Lemma 2.4 (see [19]). Let and be two symmetric matrices, and matrix has suitable dimension. Then if and only if both and .
3. Cluster Synchronization of LCODEs with Identical Nodes
The complex network we considered in this section can be described as where is the networks size, is the state vector of the th oscillator, is a continues map, and is the coupling strength. is the time delay. and are the coupling configuration matrix with zero-sum rows. It represents the topological structure of the network, in which if there is a connection from node to node () and is zero otherwise. has the same properties. Here, and need not be symmetric since asymmetric topological structures are most common in the real world.
Without loss of generality, we set the partition of nodes , , , , where , and . Let be the special solutions of the homogenous system , which satisfy for . By using similar pinning control method in [14], we let the controlled oscillator set be , which means we only put the control on the last node of each partition . We use linear negative feedback controllers and the LCODEs (3.1) become where with are the control strengths.
In this section, sufficient conditions are derived for the attainment of cluster synchronization for any initial value by control pinning, that is, by making
Theorem 3.1. Suppose that the coupling matrices and in (3.2) satisfy and . Let be a positive definite diagonal matrix and be a diagonal matrix such that . Define . For , if where satisfies Then, for any initial values, the solution , , , of the system (3.1) under the control (3.2) can achieve cluster synchronization and satisfies (3.3).
Proof. We define , where . Denote for . Since , which means that hold for all and . It is readily seen that
And similarly, . Thus, for satisfies the following differential equation:
Choose a Lyapunov function as
Note that . We differentiate (3.8) along (3.7) and have
where is defined in (3.5).
If (3.4) is satisfied, from Lemma 2.4, it can be easily seen that , and thus (3.2) could achieve cluster synchronization.
By using adaptive adjustments, we can find relatively small control strength to realize cluster synchronization. We regard the control strength of the network functions varying with time. Then, we could design the adaptive control strength. Then, we have the following result.
Theorem 3.2. Suppose that the coupling matrix and in (3.2) satisfy and . The control strength in (3.2) is defined as where is a positive constant. Let be a positive definite diagonal matrix and be a diagonal matrix such that . For , if where . is a constant and satisfies if . Then, for any initial values, the solution , , , of the system (3.1) under the control (3.2) can achieve cluster synchronization and satisfies (3.3).
Proof. Choose a Lyapunov function as Note that . We differentiate (3.12) along (3.7) and have By using similar calculations in (3.9), we have Noticing the inequalities in (3.11), we could obtain . Thus, we get , and . By Cauchy convergence principle, converges.
4. Cluster Synchronization of LCODEs with Nonidentical Nodes
The complex network considered in this section is where are the state variables of node . The complex network (4.1) has nodes and communities with . If node belongs to the th community, then we let . We denote by the set of all nodes in the th community and let , which is the subset of , be the index set of all nodes in the th community having direct connections to other communities in . And by the similar way, we can define . The function describes the local dynamics of nodes in the th community, which is differentiable and capable of performing abundant dynamical behaviors. For any pair of indices and , if , which means that node and node belong to different communities, then . The constant denotes the coupling strength. is the time delay in couplings. For and are the coupling configuration matrices with zero-sum rows, which represent the topological structure of the network. Take for an example and has the same properties. if there is a connection and is zero otherwise. Also, and need not be symmetric.
Let , , be the control inputs, then the controlled dynamical network with respect to (4.1) can be described by We define the error variables by for , where satisfies , which describes the identical local dynamics for the nodes in the th community. The nodes are said to achieve cluster synchronization if which means that the nodes within th community are fully synchronized to dynamic state , while nodes in separate communities behave independently. Define a set as the cluster synchronization manifold for network (4.2). In fact, condition (4.3) implies that the manifold is stable.
According to the above definition of error variables, we can write the corresponding error system with respect to (4.2) as
Since and are zero-row-sum matrices, we have Let stands for the feedback control strength. We design a local feedback control as where for and for , which means we put control on the nodes that have communications with nodes in other different clusters. By intuition, the terms and in (4.6) is used to synchronize all nodes in the same cluster, while the remainder terms and in the controller is to weaken the mutual influences among clusters at the intersection nodes. It is easy to verify that the manifold is an invariant manifold for the network (4.2).
Theorem 4.1. Let be a positive definite diagonal matrix and be a diagonal matrix such that for . For , if where and . is a given positive constant. Then, for any initial values, the solution , , , of the system (4.1) under the control (4.6) can achieve cluster synchronization and satisfies (4.3).
Proof. Denote , . Define a Lyapunov-Krasovskii function as
Differentiating (4.8) along the solution of error system (4.4) under the control (4.6) gives
Noticing the inequalities of (4.7), from Lemma 2.4, we obtain
Thus, we get .
5. Numerical Simulation
In this section, we give numerical simulations to verify the theorems obtained in Section 3 and Section 4. Consider the unified system, The system is chaotic for all . If , the system (5.1) is Lorenz’s attractor. The ultimate bound and positively invariant set for system (5.1) is given in [22]. For , and , where Thus, we have , and . Let and , we obtain Let , we have Set , , and . Let , and respectively in (5.1) since the nodes are nonidentical in different communities. If , we have and . If , we have and . If , we have and .
5.1. Simulation to Cluster Synchronization with Identical Nodes
In this simulation, we consider a network with 30 nodes and 3 communities. It is too high and we do not show it out. We show the topology structure in Figure 1.
(a)
(b)
(c)
(d)
Let , , be the solution of the uncoupling system with initial values , , and . Define . If , the complex network achieves cluster synchronization.
For the controlled dynamic network (4.6). Let , , , , in the Theorem 3.1. It can be easily verified that such parameters fits (4.7).
Figures 2(a), 2(b), and 2(c) give the behavior of , as well as Figure 2(d) which shows how evolve in pinning united chaotic attractor complex network with initial values chosen randomly in the interval .
(a)
(b)
(c)
(d)
For the controlled dynamic network (4.6) with adaptive control strength , let and keep other parameters the same as selected above, which makes (3.11) correct.
Figure 3 shows the behavior of and how evolve in pinning Lorenz chaotic attractor complex network with adaptive control strength and the initial values are chosen randomly in the interval .
(a)
(b)
(c)
(d)
(e)
5.2. Simulation to Cluster Synchronization with Nonidentical Nodes
In this simulation, we consider a network with 30 nodes and 3 communities. It is too high and we do not show it out. We give the topology structure Figure 4
(a)
(b)
(c)
(d)
Let , , be the solution of the uncoupling system with initial values , , and . Define . If , the complex network achieves cluster synchronization.
For the controlled dynamic network (4.6). Let , , , , and in the Theorem 4.1. It can be easily verified that such parameters fits (4.7).
Figures 5(a), 5(b), and 5(c) give the behavior of , as well as Figure 5(d) which shows how evolve in pinning united chaotic attractor complex network with initial values chosen randomly in the interval .
(a)
(b)
(c)
(d)
6. Conclusion
In the paper, we have investigated the cluster synchronization on pinning control of LCODEs with identical and nonidentical nodes. We give a sufficient condition to make the complex network achieve cluster synchronization. Moreover, adaptive feedback control techniques are used to adjust control strength. Finally, some numerical examples are given, which is essential to verify our theoretical analysis.
Acknowledgments
The authors are grateful to the editors and the reviewers for their valuable suggestions and comments. This work was supported by Guangdong Education University Industry Cooperation Projects (2009B090300355), Shenzhen Basic Research Project (JC201006010743A), and the 2011 Foundation for Distinguished Young Talents in Higher Education of Guangdong (LYM11115).