Next Article in Journal
Efficient Jacobian Computations for Complex ECT/EIT Imaging
Previous Article in Journal
Efficient Maintenance of Minimum Spanning Trees in Dynamic Weighted Undirected Graphs
Previous Article in Special Issue
Optimal Control Strategy for SLBRS with Two Control Inputs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Clustering Component Synchronization of Nonlinearly Coupled Complex Networks via Pinning Control

1
College of Electrical and Information Engineering, Lanzhou University of Technology, Lanzhou 730050, China
2
Department of Applied Mathematics, Lanzhou University of Technology, Lanzhou 730050, China
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(7), 1022; https://doi.org/10.3390/math12071022
Submission received: 17 January 2024 / Revised: 14 March 2024 / Accepted: 27 March 2024 / Published: 28 March 2024
(This article belongs to the Special Issue Control Problem of Nonlinear Systems with Applications, 2nd Edition)

Abstract

:
In this paper, the problem of clustering component synchronization of nonlinearly coupled complex networks with nonidentical nodes and asymmetric couplings is investigated. A pinning control strategy is designed to achieve the clustering component synchronization with respect to the specified components. Based on matrix analysis and stability theory, clustering component synchronization criteria are established. Two numerical simulations are also provided to show the effectiveness of the theoretical results.

1. Introduction

Synchronization, a ubiquitous and fascinating collective behavior in complex networks, has been extensively studied through the past decades [1,2,3,4,5,6,7,8]. In reality, the complex network may be split into a few subnetworks called clusters due to some partition laws. When the dynamical nodes within a cluster synchronize with each other but desynchronize among different clusters, cluster synchronization appears. Due to its potential applications in neural networks, secure communication, etc., cluster synchronization has been intensively investigated and has yielded fruitful research results [9,10,11,12,13,14,15,16,17,18,19,20,21,22]. For example, in [11], cluster synchronization on multiple nonlinearly coupled dynamical subnetworks of complex networks with nonidentical nodes and stochastic perturbations was considered by Zhou et al., and the criteria that can ensure the network synchronizes to the reference state exponentially were established. In [12], Liu et al., discussed the asymptotic and finite-time cluster synchronization for coupled fractional-order neural networks with time delay, derived sufficient conditions for both synchronizations, respectively, and estimated the upper bound of the settling time in the case of finite-time cluster synchronization. In 2020, Della Rossa et al. [13] discussed the problem of symmetries and cluster synchronization in arbitrary multilayer networks, defined and computed the symmetries of the networks, and investigated the cluster synchronization that emerged analytically and experimentally. In 2023, Lu et al. [17] studied the cluster synchronization problem of a directed reaction-diffusion complex network with stochastic noise and Markovian switching, proposed switched constant gain, centralized adaptive and decentralized adaptive pinning control strategy, respectively, to realize the cluster synchronization in mean square of the network.
The pinning control strategy, which drives a complex network to the selected synchronization state only by controlling a small portion of the nodes in the network, has been widely adopted to synchronize the network since it was first introduced by Wang and Chen [23] in 2002. Much work has been dedicated to the investigation of the synchronization of complex networks via pinning control [24,25,26,27,28,29,30,31,32,33,34]. For instance, in [24], by adding linear negative feedback controllers to the last node of each cluster, Wu et al., successfully drove a general linearly coupled network to a specified cluster synchronization pattern for any initial value and delivered the pinning cluster synchronization criteria. Feng et al. [25] extended the results of [24] to nonlinearly coupled networks with symmetric coupling and tamed the network to prescribed cluster synchronization patterns by pinning the first l u nodes of each cluster. In [26], the problem of adaptive cluster synchronization for directed networks with nonidentical nodes was addressed by Wang et al., by controlling the root nodes of all the clusters.
The synchronization mentioned above refers to the convergence of all components of the node’s state variables in the complex network. However, in the synchronization study of complex networks, especially when the complex network is composed of high-dimensional chaotic systems, the requirement of convergence of all state components is sometimes difficult or even impossible due to the sensitivity of dynamic behavior to initial conditions. Moreover, in some instances, one may only focus on the convergence of some components of the node’s state variables rather than all components. These raise the necessity of studying partial component synchronization (consensus). Ma and his team worked a lot on this issue (one can see [35,36,37,38,39] and references therein). For example, Li et al., addressed the partial component synchronization problem on connected complex networks and derived sufficient conditions ascertaining the synchronization in [36]. Moreover, in [37], the authors studied the clustering component synchronization problem of unconnected complex networks. With the help of pinning control, the first m components of the node’s state variables in each cluster of the network were synchronized to the presetting synchronous trajectories.
It is worth noting that the clustering component synchronization studied in [37] is on a linearly coupled network, with the assumption that the coupling matrix is symmetric. In reality, however, most networks are directed and the state variables cannot be observed directly. Therefore, the study of clustering component synchronization for nonlinearly coupled complex networks with asymmetric couplings is necessary in both theoretical research and applications. Enlightened by the literature mentioned above, in this paper, we investigate the clustering component synchronization problem for nonlinearly coupled complex networks via pinning control. The novelty of this paper is that the synchronization studied herein is the convergence on any m specified components of the node’s state variables of the network rather than all components. For the purpose that the synchronization behavior of the specified components can be investigated, we move these components to the front by rearranging the components of error variables through a suitable transformation. Hence, the investigation of the clustering component synchronization with respect to these specified components transforms into the study of the stability of partial variables. Then, the stability theory of partial variables is applied, and sufficient conditions ensuring the realization of the clustering component synchronization are established. The main contributions of this paper can be summarized as follows:
(1)
The inner coupling of the complex networks studied in this paper is nonlinear, which is more realistic since the observed data are usually a nonlinear function of the state variable rather than itself.
(2)
The coupling matrix of the complex network investigated in this paper is asymmetric, which means that the complex network is directed. Moreover, the intrinsic dynamics of the nodes are uniform within a cluster but not the same in different ones. Therefore, the complex network under investigation in this paper is more in line with reality.
(3)
The synchronization investigated in this paper is clustering component synchronization, which is concerned with the asymptotic convergence of specified components of the node’s state variables in each cluster of the network rather than all components. When the specified components are all state components, clustering component synchronization becomes cluster synchronization in general. Hence, compared with traditional cluster synchronization investigated in [24,26], clustering component synchronization has more potential applications.
The rest of this paper is arranged as follows. In Section 2, some notations, definitions and lemmas that will be used later are presented. In Section 3, the criteria for clustering component synchronization are obtained. Two numerical simulations are provided to verify the effectiveness of the theoretical results in Section 4, and the paper is concluded in Section 5.

2. Preliminaries

In this section, we provide some notations, definitions and lemmas that will be used later. Throughout this paper, R + , R n and R n × n denote the sets of nonnegative real numbers, n-dimensional real column vectors with Euclidean norm · and n × n real matrices, respectively. I denotes the identity matrix with appropriate dimension, diag ( d 1 , d 2 , , d n ) denotes the diagonal matrix with diagonal entries d 1 to d n . The superscript “T” stands for the transpose of a matrix and the notation ⊗ represents the Kronecker product. A s means the symmetric part of a square matrix A which is defined as A s = 1 2 ( A + A T ) . For symmetric matrix B, B < 0 means that B is negative definite and λ max ( B ) means the maximum eigenvalue of B.
What follows are some basic definitions and lemmas.
Definition 1 
([9]). An irreducible matrix A = ( a i j ) R q × q is said to belong to class A , denoted by A A , if the elements besides the diagonal entries are nonnegative and the diagonal elements a i i = j = 1 , j i q a i j , i = 1 , 2 , , q .
Definition 2 
([9]). For an L × L matrix A, which can be written as
A = A 11 A 12 A 1 k A 21 A 22 A 2 k A k 1 A k 2 A k k ,
where A u v R ( q u q u 1 ) × ( q v q v 1 ) , u , v = 1 , 2 , , k , here, q 0 = 0 , q k = L , is said to belong to class M , denoted by A M , if each row-sum of A u v is zero and A u u A , u , v = 1 , 2 , , k .
Definition 3 
([1]). A nonlinear function h : R R is said to belong to the acceptable nonlinear coupling function class, denoted by h N C F ( β , δ ) , if there are nonnegative scalars β and δ, such that h ( υ ) β υ satisfies the following Lipschitz condition:
| h ( υ 1 ) h ( υ 2 ) β ( υ 1 υ 2 ) | δ | υ 1 υ 2 | ,
where υ 1 , υ 2 R .
Definition 4 
([37,40]). A function ϕ : R + R + is said to belong to the K class function, denoted by ϕ K , if ϕ is continuous, strictly monotone increasing and ϕ ( 0 ) = 0 .
Consider the s-dimensional non-autonomous system
x ˙ = f ( t , x ) , t R + ,
where x = ( y T , z T ) T , here y = ( x 1 , x 2 , , x l ) T and z = ( x l + 1 , x l + 2 , , x s ) T , f C [ R + × R s , R s ] and f ( t , 0 ) 0 for t R + .
Assume that the existence and uniqueness of solutions to (1) subject to x ( t 0 ) = x 0 as well as their dependence on initial values are guaranteed.
Lemma 1 
([37,40]). Let ϕ , ψ , α K . If there is a Lyapunov function V : R + × R s R + with V ( t , 0 ) = 0 for t R + , such that
ϕ ( y ) V ( t , x ) ψ ( y ) for ( t , x ) ( R + , R n ) ,
and its derivative along the trajectories of (1) meets
d V d t ( 1 ) α ( y ( t ) ) , t R + ,
then the trivial solution of the system (1) is asymptotically stable with respect to the variable y.
Lemma 2. 
Let B = ( b i j ) R L × L and C = ( c i j ) R n × n . Then, for any permutation ( p 1 , p 2 , , p n ) of ( 1 , 2 , , n ) , there exist orthogonal matrices P R L n × L n and Q R n × n such that the equality
P ( B C ) P T = ( Q C Q T ) B
holds.
Proof. 
For any permutation ( p 1 , p 2 , , p n ) of ( 1 , 2 , , n ) , we define
P = ( θ p 1 ) T , , ( θ ( L 1 ) n + p 1 ) T , ( θ p 2 ) T , , ( θ ( L 1 ) n + p 2 ) T , , ( θ p n ) T , , ( θ ( L 1 ) n + p n ) T T
and
Q = ( ϑ p 1 ) T , ( ϑ p 2 ) T , , ( ϑ p n ) T T ,
where θ l ( l = 1 , 2 , , L n ) is an L n -dimensional row vector whose lth element is 1 and all the other elements are 0, while ϑ s ( s = 1 , 2 , , n ) is an n-dimensional row vector whose sth element is 1 and all the other elements are 0. Then, it is obvious that P and Q are orthogonal matrices, and (4) is true after direct calculation. □
Lemma 3 
([41]). For any ξ , η R L and positive definite matrix K R L × L ,
ξ T η 1 2 ξ T K ξ + 1 2 η T K 1 η .
Lemma 4 
([42]). For symmetric matrix W = ( w i j ) R L × L , if each row-sum of W is zero, then for any two vectors ξ = ( ξ 1 , ξ 2 , , ξ L ) T and η = ( η 1 , η 2 , , η L ) T ,
ξ T W η = i = 1 L j = 1 L ξ i w i j η j = j > i w i j ( ξ j ξ i ) ( η j η i ) .
Lemma 5 
([43]). If A , B R L × L are symmetric, ξ R L is a nonzero vector, a , b , c R and a , b > 0 , then
(1)
c A is symmetric;
(2)
A ± B is symmetric;
(3)
ξ T A ξ ξ T ξ λ max ( A ) ;
(4)
λ max ( a A + b B ) a λ max ( A ) + b λ max ( B ) .

3. Main Results

Consider a nonlinearly coupled complex network composed of L nodes with index set L = { 1 , 2 , , L } . Suppose these nodes are split into k ( 1 k < L ) nonempty clusters. Without loss of generality, let { C 1 , C 2 , , C k } be the partition, where C l = { q l 1 + 1 , q l 1 + 2 , , q l } , here q 0 = 0 , q k = L , q l 1 < q l , l = 1 , 2 , , k . Then the network can be described as
x ˙ i ( t ) = f i ¯ ( x i ( t ) ) + j = 1 L c i j a i j h ( x j ( t ) ) , t R + , i L ,
where x i ( t ) = x 1 i ( t ) , x 2 i ( t ) , , x n i ( t ) T R n is the state variable of the node i at time t; i ¯ denotes the index of the cluster that the node i belongs to; f i ¯ : R n R n is a continuous function which describes the intrinsic dynamic of the nodes in the i ¯ th cluster and f i ¯ f j ¯ for i ¯ j ¯ ; h : R n R n is the nonlinear coupling function, which is defined by h ( ζ ) = h 1 ( ζ 1 ) , h 2 ( ζ 2 ) , , h n ( ζ n ) T for ζ = ( ζ 1 , ζ 2 , , ζ n ) T R n and h satisfies h ( 0 ) = 0 . A = ( a i j ) R L × L denotes the topological structure of the complex network, which is defined as follows: if there is a connection from node j to node i( i j ), then a i j 0 ; otherwise, a i j = 0 and a i i = j = 1 , j i L a i j , i L . Here, A can be asymmetric or reducible. c i j > 0 represents the coupling strength between node i and node j, which is defined as follows: if i ¯ = j ¯ = l , then c i j = c l > 0 ; if i ¯ j ¯ , then c i j = 1 , i , j L . Denote G : = ( c i j a i j ) R L × L . According to the partition, G can be written as:
G = G 11 G 12 G 1 k G 21 G 22 G 2 k G k 1 G k 2 G k k = c 1 A 11 A 12 A 1 k A 21 c 2 A 22 A 2 k A k 1 A k 2 c k A k k ,
where G u v , A u v R ( q u q u 1 ) × ( q v q v 1 ) , u , v = 1 , 2 , , k . In this paper, we always assume that A M .
Suppose π 1 ( t ) , π 2 ( t ) , , π k ( t ) are the target trajectories of the k clusters, satisfying π ˙ l ( t ) = f l π l ( t ) , and π l ( t ) π j ( t ) for l j , t R + , l , j = 1 , 2 , , k . Here, we assume that the existence and uniqueness of solutions to π ˙ l ( t ) = f l π l ( t ) subject to π l ( t 0 ) = π l 0 and their dependence on initial values are guaranteed. Select any m ( 1 m n ) components of the node’s state variables as the components that are required to be synchronized. Denote these m components as p 1 , p 2 , , p m , and the remaining components as p m + 1 , p m + 2 , , p n . The objective of this paper is to design an appropriate pinning control strategy such that the network (5) can achieve clustering component synchronization with respect to the specified components. To this end, we add controllers on the first r u ( 1 r u q u q u 1 ) nodes of the uth cluster, u = 1 , 2 , , k .
In what follows, we consider the pinning controlled network
x ˙ i ( t ) = f i ¯ x i ( t ) + j = 1 L c i j a i j h x j ( t ) + u i ( t ) , t R + , i J , x ˙ i ( t ) = f i ¯ x i ( t ) + j = 1 L c i j a i j h x j ( t ) , t R + , i J ,
where J = { 1 , , r 1 , q 1 + 1 , , q 1 + r 2 , , q k 1 + 1 , , q k 1 + r k } , u i ( t ) is the controller defined by
u i ( t ) = c i ¯ d i ¯ h x i ( t ) h ( π i ¯ t ) , t R + , i J ,
with d i ¯ > 0 . For convenience, if i J , then let d i ¯ = 0 . Hence, (6) can be rewritten as
x ˙ i ( t ) = f i ¯ x i ( t ) + j = 1 L c i j a i j h x j ( t ) + u i ( t ) , t R + , i L .
Define the error variable e i ( t ) = x i ( t ) π i ¯ ( t ) , t R + , i L . Define
e ^ p l ( t ) = e p l 1 ( t ) , e p l 2 ( t ) , , e p l L ( t ) T , t R + , l = 1 , 2 , , n .
Then, we have the following definition of clustering component synchronization with respect to the specified components p 1 , p 2 , , p m :
Definition 5. 
If lim t + l = 1 m e ^ p l ( t ) = 0 , then the pinning controlled network (8) is said to achieve clustering component synchronization with respect to the specified components p 1 , p 2 , , p m .
In order to derive the synchronization criteria for the pinning controlled network (8), we make the following assumptions:
(A1) There exist constants ϖ u > 0 such that for any ξ , η R n , the inequalities
( ξ η ) T Λ f u ( ξ ) f u ( η ) ϖ u ( ξ η ) T Λ ( ξ η ) , u = 1 , 2 , , k
hold, where Λ = diag ( λ 1 , λ 2 , , λ n ) , here
λ j = 1 , j { p 1 , p 2 , , p m } , 0 , o t h e r w i s e ;
(A2) For j { p 1 , p 2 , , p m } , there exist constants β j δ j 0 such that h j N C F ( β j , δ j ) .
Now, we give a general sufficient condition for the pinning controlled network (8) to achieve clustering component synchronization with respect to the specified components, which implies the relationship between the synchronization of the specified components and the topological structure of the network, strength of couplings, self-dynamics of the isolated node, inner couplings and the control.
Theorem 1. 
Suppose that (A1) and (A2) hold. If the following conditions are satisfied:
( ϖ u + ( 1 1 q u q u 1 ) β p l ( k 1 ) + ( 1 1 q u q u 1 ) k δ p l 2 ) I + 1 2 β p l v = 1 , v u k G u v G u v T + 1 2 v = 1 k G u v G u v T + β p l ( G u u s D u ) + δ p l D u < 0 , u = 1 , 2 , , k , l = 1 , 2 , , m ,
where D u = diag ( c u d u , , c u d u r u , 0 , , 0 q u q u 1 r u ) , then the pinning controlled network (8) achieves clustering component synchronization with respect to the specified components p 1 , p 2 , , p m .
Proof. 
In view of A M , it follows
j = 1 L c i j a i j h x j ( t ) = u = 1 k j C u c i j a i j h x j ( t ) h π j ¯ ( t ) + u = 1 k j C u c i j a i j h π j ¯ ( t ) = j = 1 L c i j a i j h x j ( t ) h π j ¯ ( t ) .
Therefore,
e ˙ i ( t ) = f i ¯ x i ( t ) f i ¯ π i ¯ ( t ) + j = 1 L c i j a i j h x j ( t ) h π j ¯ ( t ) c i ¯ d i ¯ h x i ( t ) h ( π i ¯ t ) , t R + , i L .
Denote E ( t ) = e 1 ( t ) T , e 2 ( t ) T , , e L ( t ) T T , D = diag ( D 1 , D 2 , , D k ) , where D u = diag ( c u d u , , c u d u r u , 0 , , 0 q u q u 1 r u ) , u = 1 , 2 , , k . Then, (10) can be rewritten as the following compact form
E ˙ ( t ) = F x 1 ( t ) , x 2 ( t ) , , x L ( t ) F π 1 ¯ ( t ) , π 2 ¯ ( t ) , , π L ¯ ( t ) + ( G D ) I H x 1 ( t ) , x 2 ( t ) , , x L ( t ) H π 1 ¯ ( t ) , π 2 ¯ ( t ) , , π L ¯ ( t ) , t R + ,
where F and H are defined by
F ( ζ 1 , ζ 2 , , ζ L ) = f 1 1 ¯ ( ζ 1 ) , f 2 1 ¯ ( ζ 1 ) , , f n 1 ¯ ( ζ 1 ) , f 1 2 ¯ ( ζ 2 ) , f 2 2 ¯ ( ζ 2 ) , , f n 2 ¯ ( ζ 2 ) , , f 1 L ¯ ( ζ L ) , f 2 L ¯ ( ζ L ) , , f n L ¯ ( ζ L ) T
and
H ( ζ 1 , ζ 2 , , ζ L ) = h 1 ( ζ 1 1 ) , h 2 ( ζ 2 1 ) , , h n ( ζ n 1 ) , h 1 ( ζ 1 2 ) , h 2 ( ζ 2 2 ) , , h n ( ζ n 2 ) , , h 1 ( ζ 1 L ) , h 2 ( ζ 2 L ) , , h n ( ζ n L ) T
for ζ i = ( ζ 1 i , ζ 2 i , , ζ n i ) T , i L , respectively.
Define E ^ ( t ) = e ^ p 1 ( t ) T , e ^ p 2 ( t ) T , , e ^ p n ( t ) T T , t R + . It is obvious that E ^ ( t ) = P E ( t ) , where P is defined as in Lemma 2, and
E ^ ˙ ( t ) = F ^ x 1 ( t ) , x 2 ( t ) , , x L ( t ) F ^ π 1 ¯ ( t ) , π 2 ¯ ( t ) , , π L ¯ ( t ) + I ( G D ) H ^ x 1 ( t ) , x 2 ( t ) , , x L ( t ) H ^ π 1 ¯ ( t ) , π 2 ¯ ( t ) , , π L ¯ ( t ) , t R + ,
where F ^ and H ^ are defined by
F ^ ( ζ 1 , ζ 2 , , ζ L ) = f ^ p 1 ( ζ 1 , ζ 2 , , ζ L ) T , , f ^ p n ( ζ 1 , ζ 2 , , ζ L ) T T = f p 1 1 ¯ ( ζ 1 ) , f p 1 2 ¯ ( ζ 2 ) , , f p 1 L ¯ ( ζ L ) , , f p n 1 ¯ ( ζ 1 ) , f p n 2 ¯ ( ζ 2 ) , , f p n L ¯ ( ζ L ) T
and
H ^ ( ζ 1 , ζ 2 , , ζ L ) = h ^ p 1 ( ζ 1 , ζ 2 , , ζ L ) T , , h ^ p n ( ζ 1 , ζ 2 , , ζ L ) T T = h p 1 ( ζ p 1 1 ) , h p 1 ( ζ p 1 2 ) , , h p 1 ( ζ p 1 L ) , , h p n ( ζ p n 1 ) , h p n ( ζ p n 2 ) , , h p n ( ζ p n L ) T
for ζ i = ( ζ 1 i , ζ 2 i , , ζ n i ) T , i L , respectively.
Now, we construct the following Lyapunov function
V ( t , x ) = 1 2 x T ( Γ I ) x for ( t , x ) R + × R L n ,
where x = ( y T , z T ) T , here y = ( x 1 , x 2 , , x L m ) T and z = ( x L m + 1 , x L m + 2 , , x L n ) T ; Γ = diag ( 1 , , 1 m , 0 , , 0 n m ) .
Obviously, V ( t , x ) = 1 2 x T ( Γ I ) x = 1 2 i = 1 L m ( x i ) 2 = 1 2 y T y = 1 2 y 2 . Thus, V : R + × R L n R + and V ( t , 0 ) = 0 for t R + . Then, for chosen functions ϕ ( ν ) = 1 4 ν 2 and ψ ( ν ) = ν 2 , ν R + , V ( t , x ) meets (2) of Lemma 1.
Differentiating V ( t , x ) along the trajectories of (11), we obtain
d V d t ( 11 ) = E ^ ( t ) T ( Γ I ) E ^ ˙ ( t ) = E ^ ( t ) T ( Γ I ) ( F ^ x 1 ( t ) , x 2 ( t ) , , x L ( t ) F ^ π 1 ¯ ( t ) , π 2 ¯ ( t ) , , π L ¯ ( t ) + I ( G D ) H ^ x 1 ( t ) , x 2 ( t ) , , x L ( t ) H ^ π 1 ¯ ( t ) , π 2 ¯ ( t ) , , π L ¯ ( t ) ) = V 1 ( t ) + V 2 ( t ) + V 3 ( t ) , t R + ,
where
V 1 ( t ) = E ^ ( t ) T ( Γ I ) F ^ x 1 ( t ) , x 2 ( t ) , , x L ( t ) F ^ π 1 ¯ ( t ) , π 2 ¯ ( t ) , , π L ¯ ( t ) , t R + ,
V 2 ( t ) = E ^ ( t ) T ( Γ G ) H ^ x 1 ( t ) , x 2 ( t ) , , x L ( t ) H ^ π 1 ¯ ( t ) , π 2 ¯ ( t ) , , π L ¯ ( t ) , t R +
and
V 3 ( t ) = E ^ ( t ) T ( Γ D ) H ^ x 1 ( t ) , x 2 ( t ) , , x L ( t ) H ^ π 1 ¯ ( t ) , π 2 ¯ ( t ) , , π L ¯ ( t ) , t R + .
For the sake of convenience, denote e ^ p l u ( t ) = e p l q u 1 + 1 ( t ) , e p l q u 1 + 2 ( t ) , , e p l q u ( t ) T , t R + , u = 1 , 2 , , k , l = 1 , 2 , , n . Then, by (A1), we have
V 1 ( t ) = E ^ ( t ) T ( Γ I ) F ^ x 1 ( t ) , x 2 ( t ) , , x L ( t ) F ^ π 1 ¯ ( t ) , π 2 ¯ ( t ) , , π L ¯ ( t ) = l = 1 m e ^ p l ( t ) T f ^ p l x 1 ( t ) , x 2 ( t ) , , x L ( t ) f ^ p l π 1 ¯ ( t ) , π 2 ¯ ( t ) , , π L ¯ ( t ) = l = 1 m u = 1 k i C u e p l i ( t ) f p l i ¯ x i ( t ) f p l i ¯ π i ¯ ( t ) = u = 1 k i C u j = 1 n e j i ( t ) λ j f j i ¯ x i ( t ) f j i ¯ π i ¯ ( t ) = u = 1 k i C u e i ( t ) T Λ f i ¯ x i ( t ) f i ¯ π i ¯ ( t ) u = 1 k i C u ϖ i ¯ e i ( t ) T Λ e i ( t ) = l = 1 m u = 1 k ϖ u e ^ p l u ( t ) T e ^ p l u ( t ) , t R + .
Let W u = ( w i j ) = I 1 q u q u 1 1 · 1 T , where 1 = ( 1 , 1 , , 1 q u q u 1 ) T , u = 1 , 2 , , k . It is easy to verify that W is symmetric and that each row-sum of W is zero. Moreover, for any matrix G u v , it follows G u v W v = G u v , u , v = 1 , 2 , , k . Denote h ^ p l u ( ζ 1 , ζ 2 , , ζ q u q u 1 ) = h p l ( ζ p l 1 ) , h p l ( ζ p l 2 ) , , h p l ( ζ p l q u q u 1 ) T for ζ i = ( ζ 1 i , ζ 2 i , , ζ n i ) T , i = 1 , 2 , , q u q u 1 , u = 1 , 2 , , k , l = 1 , 2 , , n .
In consideration of Lemmas 3, 4 and (A2), we can obtain
V 2 ( t ) = E ^ ( t ) T ( Γ G ) H ^ x 1 ( t ) , x 2 ( t ) , , x L ( t ) H ^ π 1 ¯ ( t ) , π 2 ¯ ( t ) , , π L ¯ ( t ) = l = 1 m e ^ p l ( t ) T G h ^ p l x 1 ( t ) , x 2 ( t ) , , x L ( t ) h ^ p l π 1 ¯ ( t ) , π 2 ¯ ( t ) , , π L ¯ ( t ) = l = 1 m u = 1 k β p l e ^ p l u ( t ) T G u u e ^ p l u ( t ) + l = 1 m u = 1 k v = 1 , v u k β p l e ^ p l u ( t ) T G u v W v e ^ p l v ( t ) + l = 1 m u = 1 k v = 1 k e ^ p l u ( t ) T G u v W v ( h ^ p l v x q v 1 + 1 ( t ) , x q v 1 + 2 ( t ) , , x q v ( t ) h ^ p l v π v ( t ) , π v ( t ) , , π v ( t ) β p l e ^ p l v ( t ) ) l = 1 m u = 1 k β p l e ^ p l u ( t ) T G u u s e ^ p l u ( t ) + l = 1 m u = 1 k v = 1 , v u k ( 1 2 β p l e ^ p l u ( t ) T G u v G u v T e ^ p l u ( t ) + 1 2 β p l e ^ p l v ( t ) T W v T W v e ^ p l v ( t ) ) + l = 1 m u = 1 k v = 1 k ( 1 2 e ^ p l u ( t ) T G u v G u v T e ^ p l u ( t ) + 1 2 ( h ^ p l v x q v 1 + 1 ( t ) , x q v 1 + 2 ( t ) , , x q v ( t ) h ^ p l v π v ( t ) , π v ( t ) , , π v ( t ) β p l e ^ p l v ( t ) ) T W v T W v ( h ^ p l v x q v 1 + 1 ( t ) , x q v 1 + 2 ( t ) , , x q v ( t ) h ^ p l v π v ( t ) , π v ( t ) , , π v ( t ) β p l e ^ p l v ( t ) ) ) = l = 1 m u = 1 k e ^ p l u ( t ) T β p l G u u s + 1 2 β p l v = 1 , v u k G u v G u v T + 1 2 v = 1 k G u v G u v T e ^ p l u ( t ) + 1 2 l = 1 m u = 1 k v = 1 , v u k i , j C v , j > i β p l 1 q v q v 1 e p l j ( t ) e p l i ( t ) 2 + 1 2 l = 1 m u = 1 k v = 1 k i , j C v , j > i 1 q v q v 1 ( h p l x p l j ( t ) h p l π p l j ¯ ( t ) β p l e p l j ( t ) h p l x p l i ( t ) h p l π p l i ¯ ( t ) β p l e p l i ( t ) ) 2 l = 1 m u = 1 k e ^ p l u ( t ) T β p l G u u s + 1 2 β p l v = 1 , v u k G u v G u v T + 1 2 v = 1 k G u v G u v T e ^ p l u ( t ) + l = 1 m u = 1 k v = 1 , v u k i C v β p l ( 1 1 q v q v 1 ) e p l i ( t ) 2 + l = 1 m u = 1 k v = 1 k i C v 1 q v q v 1 δ p l 2 e p l j ( t ) 2 + e p l i ( t ) 2 = l = 1 m u = 1 k e ^ p l u ( t ) T ( β p l G u u s + 1 2 β p l v = 1 , v u k G u v G u v T + 1 2 v = 1 k G u v G u v T + ( 1 1 q u q u 1 ) β p l ( k 1 ) I + ( 1 1 q u q u 1 ) k δ p l 2 I ) e ^ p l u ( t ) , t R + .
Similarly, with (A2), we have
V 3 ( t ) = E ^ ( t ) T ( Γ D ) H ^ x 1 ( t ) , x 2 ( t ) , , x L ( t ) H ^ π 1 ¯ ( t ) , π 2 ¯ ( t ) , , π L ¯ ( t ) = l = 1 m e ^ p l ( t ) T D h ^ p l x 1 ( t ) , x 2 ( t ) , , x L ( t ) h ^ p l π 1 ¯ ( t ) , π 2 ¯ ( t ) , , π L ¯ ( t ) = l = 1 m u = 1 k β p l e ^ p l u ( t ) T D u e ^ p l u ( t ) l = 1 m u = 1 k i C u e p l i ( t ) c i ¯ d i ¯ h p l x p l i ( t ) h p l π p l i ¯ ( t ) β p l e ^ p l u ( t ) l = 1 m u = 1 k β p l e ^ p l u ( t ) T D u e ^ p l u ( t ) + l = 1 m u = 1 k i C u c i ¯ d i ¯ δ p l e p l i ( t ) 2 = l = 1 m u = 1 k ( β p l δ p l ) e ^ p l u ( t ) T D u e ^ p l u ( t ) , t R + .
Substituting inequalities (13), (14), (15) into (12), we obtain
d V d t ( 11 ) l = 1 m u = 1 k e ^ p l u ( t ) T ( ϖ u + ( 1 1 q u q u 1 ) β p l ( k 1 ) + ( 1 1 q u q u 1 ) k δ p l 2 I + 1 2 β p l v = 1 , v u k G u v G u v T + 1 2 v = 1 k G u v G u v T + β p l ( G u u s D u ) + δ p l D u ) e ^ p l u ( t ) , t R + ,
which together with (9) implies that
d V d t ( 11 ) l = 1 m u = 1 k λ max ( ϖ u + ( 1 1 q u q u 1 ) β p l ( k 1 ) + ( 1 1 q u q u 1 ) k δ p l 2 I + 1 2 β p l v = 1 , v u k G u v G u v T + 1 2 v = 1 k G u v G u v T + β p l ( G u u s D u ) + δ p l D u ) e ^ p l u ( t ) T e ^ p l u ( t ) h l = 1 m u = 1 k e ^ p l u ( t ) T e ^ p l u ( t ) , t R + ,
where h = max 1 l m , 1 u k λ max ϖ u + ( 1 1 q u q u 1 ) β p l ( k 1 ) + ( 1 1 q u q u 1 ) k δ p l 2 I + 1 2 β p l v = 1 , v u k G u v G u v T + 1 2 v = 1 k G u v G u v T + β p l ( G u u s D u ) + δ p l D u .
Choose α ( ν ) = h ν 2 , ν R + . Obviously, α K and (3) of Lemma 1 is satisfied. Then, from Lemma 1, the trivial solution of the network (11) is asymptotically stable with respect to the variable y. Therefore, lim t + l = 1 m e ^ p l ( t ) = 0 . Thus, we have the conclusion that the pinning controlled network (8) achieves clustering component synchronization with respect to the specified components p 1 , p 2 , , p m . □
Although Theorem 1 provides a very general sufficient condition, it is not convenient to use. To make Theorem 1 more applicable, we give the following corollary.
Corollary 1. 
Suppose that (A1) and (A2) hold. If
ϖ u + ( 1 1 q u q u 1 ) β p l ( k 1 ) + ( 1 1 q u q u 1 ) k δ p l 2 + λ max ( 1 2 β p l v = 1 , v u k G u v G u v T + 1 2 v = 1 k G u v G u v T + β p l ( G u u s D u ) + δ p l D u ) < 0 , u = 1 , 2 , , k , l = 1 , 2 , , m ,
then the pinning controlled network (8) achieves clustering component synchronization with respect to the specified components p 1 , p 2 , , p m .
Remark 1. 
If the specified components are the first m components of the node’s state variables in each cluster of the complex network, then the clustering component synchronization we discussed turns into clustering component synchronization with respect to the first m components in [37]. In fact, for a linearly coupled complex network with identical nodes and symmetric couplings, we can easily obtain the result of [37] from Corollary 1, which means that the clustering component synchronization we discuss is a generalization of that in [37].
Remark 2. 
If the specified components are all state components of the node, then the clustering component synchronization with respect to the specified components becomes cluster synchronization in general.
If the specified components are all state components, then (A1) and (A2) turn into (A1′) There exist constants ϖ u > 0 such that for any ξ , η R n , the inequalities
( ξ η ) T f u ( ξ ) f u ( η ) ϖ u ( ξ η ) T ( ξ η ) , u = 1 , 2 , , k
hold.
(A2′) There exist constants β j δ j 0 such that h j N C F ( β j , δ j ) , j = 1 , 2 , , n .
Therefore, we have the following corollary about cluster synchronization from Theorem 1.
Corollary 2. 
Suppose that (A1′), (A2′) hold, and the following conditions are satisfied:
( ϖ u + ( 1 1 q u q u 1 ) β j ( k 1 ) + ( 1 1 q u q u 1 ) k δ j 2 ) I + 1 2 β j v = 1 , v u k G u v G u v T + 1 2 v = 1 k G u v G u v T + β j ( G u u s D u ) + δ j D u < 0 , u = 1 , 2 , , k , j = 1 , 2 , , n .
Then, the pinning controlled network (8) achieves cluster synchronization.

4. Numerical Simulations

In this section, two numerical simulations are provided to show the effectiveness of the theoretical results.
Example 1. 
Consider a nonlinearly coupled complex network composed of nine nodes with index set L = { 1 , 2 , , 9 } . Suppose these nodes are split into two clusters: C 1 = { 1 , 2 , , 6 } and C 2 = { 7 , 8 , 9 } . Therefore, the network can be described as
x ˙ i ( t ) = f i ¯ x i ( t ) + j = 1 9 c i j a i j h x j ( t ) , t R + , i L ,
where x i ( t ) = x 1 i ( t ) , x 2 i ( t ) , x 3 i ( t ) T , i L . f i ¯ and h are defined by
f 1 ( ζ ) = ( 0.1 ζ 1 2 ζ 1 3 tanh ζ 2 ζ 3 , 1.2 ζ 1 + ζ 2 , 0.15 ζ 1 0.1 ζ 3 + 0.01 ) T ,
f 2 ( ζ ) = ( 0.5 ζ 1 2 ζ 1 3 sin ζ 2 ζ 3 , 2.3 ζ 2 + 1 , 0.11 ζ 1 0.1 ζ 3 + 0.02 ) T
and
h ( ζ ) = ( 10 ζ 1 + sin ζ 1 , ζ 2 + tanh ζ 2 , 10 ζ 3 + sin ζ 3 , ) T
for ζ = ( ζ 1 , ζ 2 , ζ 3 ) T , respectively. c 1 = 3 , c 2 = 2.5 and the coupling configuration matrix
A = A 11 A 12 A 21 A 22 = 1.9 0 0.8 0 1.1 0 0 0 0 1.3 1.3 0 0 0 0 0 0 0 0 1.8 1.8 0 0 0 0 0 0 0 0 0.8 1.4 0 0.6 0 0 0 0.8 0 0 0.5 1.3 0 0 0 0 1.3 0 0 0 0 1.3 0 0 0 0 0 0 1.4 0 1.4 1 0 1 0 0 0 0 0 0 0.9 2 1.1 0 0 0 0 0 0 0.8 1.2 2 .
Let π 1 ( t ) and π 2 ( t ) be the target trajectories of two clusters with initial values π 1 ( 0 ) = ( 1 , 1 , 1 ) T and π 2 ( 0 ) = ( 1 , 2 , 3 ) T , respectively. Now, we investigate the pinning controlled network
x ˙ i ( t ) = f i ¯ x i ( t ) + j = 1 9 c i j a i j h x j ( t ) + u i ( t ) , t R + , i L ,
where u i ( t ) is defined by (7) with D = diag ( 30 , 30 , 0 , 0 , 0 , 0 , 20 , 0 , 0 ) .
Let p 1 = 1 and p 2 = 3 . In what follows, we will verify that the pinning controlled network (17) can achieve clustering component synchronization with respect to the specified components p 1 and p 2 .
Obviously, A M . It is easy to verify that (A1) holds with ϖ 1 = 3.9 and ϖ 2 = 4.3 , (A2) holds with β p l = 10 and δ p l = 1 , l = 1 , 2 . Moreover, ϖ u + ( 1 1 q u q u 1 ) β p l ( k 1 ) + ( 1 1 q u q u 1 ) k δ p l 2 + λ max 1 2 β p l v = 1 , v u k G u v G u v T + 1 2 v = 1 k G u v G u v T + β p l ( G u u s D u ) + δ p l D u is equal to 0.2283 and 0.1739 for u = 1 and u = 2 , respectively, l = 1 , 2 . All the requirements of Corollary 1 are met. Therefore, it follows Corollary 1 that the network (17) achieves clustering component synchronization with respect to the specified components p 1 = 1 and p 2 = 3 . Figure 1, Figure 2 and Figure 3 show the time history of the components of the error variables corresponding to the pinning controlled network (17). It can be seen that the first and third components of the error variables tend to 0 when t + , respectively, while the second component of the error variables does not. All these indicate that the network (17) achieves clustering component synchronization with respect to the specified components p 1 and p 2 .
Example 2. 
Consider a nonlinearly coupled complex network composed of 10 nodes with index set L = { 1 , 2 , , 10 } . Suppose these nodes are split into two clusters: C 1 = { 1 , 2 , , 5 } and C 2 = { 6 , 7 , , 10 } , and the node dynamics of the two clusters, respectively, are the 3-D neural network and the Chua’s circuit, which are described by
ζ ˙ = f 1 ( ζ ) = I ζ + B g ( ζ ) , ζ C 1
and
ζ ˙ = f 2 ( ζ ) = T ζ + w ( ζ ) , ζ C 2 .
where ζ = ( ζ 1 , ζ 2 , ζ 3 ) T , g ( ζ ) = 1 2 ( | ζ 1 + 1 | | ζ 1 1 | ) , 1 2 ( | ζ 2 + 1 | | ζ 2 1 | ) , 1 2 ( | ζ 3 + 1 | | ζ 3 1 | ) T , w ( ζ ) = 27 14 ( | ζ 1 + 1 | | ζ 1 1 | ) , 0 , 0 T , and matrices B and T are, respectively, defined by
B = 1.25 3.2 3.2 3.2 1.1 4.4 3.2 4.4 1 , T = 18 7 9 0 1 1 1 0 100 7 0 .
Choose c 1 = 2 , c 2 = 2.2 , the coupling configuration matrix
A = A 11 A 12 A 21 A 22 = 3.7 1.3 0 0.8 1.6 0.3 0.4 0 0.7 0 0 1.2 0.5 0 0.7 0 0 0 0 0 1.5 0 1.5 0 0 0 0 0 0 0 1.6 0 0 1.6 0 0 0 0 0 0 0 2.1 0 0 2.1 0 0 0 0 0 0 0 0 0 0 2.9 0.7 0.8 0 1.4 0 0 0 0 0 0 0.5 0 0.5 0 0 0 0 0 0 1.2 0 1.2 0 0 0 0 0 0 0 0 0.8 1 2.9 1.1 0 0 0 0 0 1.3 0 0 1.2 2.5 ,
and the nonlinear coupling function h ( ζ ) = ( 9 ζ 1 + tanh ζ 1 , 9 ζ 2 + tanh ζ 2 , 9 ζ 3 + tanh ζ 3 , ) T for ζ = ( ζ 1 , ζ 2 , ζ 3 ) T .
Let π 1 ( t ) and π 2 ( t ) be the target trajectories with initial values π 1 ( 0 ) = ( 0.1 , 0.2 , 0.3 ) T and π 2 ( 0 ) = ( 1 , 0.1 , 1 ) T , respectively. Now, we investigate the pinning controlled network
x ˙ i ( t ) = f i ¯ x i ( t ) + j = 1 10 c i j a i j h x j ( t ) + u i ( t ) , t R + , i L ,
where u i ( t ) is defined by (7) with D = diag ( 12 , 12 , 12 , 0 , 0 , 22 , 22 , 22 , 0 , 0 ) .
Let p 1 = 1 , p 2 = 2 and p 3 = 3 . In what follows, we will verify that the pinning controlled network (18) can achieve cluster synchronization.
It is easy to verify that A M , (A1) holds with ϖ 1 = 5.6 and ϖ 2 = 10 , and (A2) holds with β p l = 9 and δ p l = 1 , l = 1 , 2 , 3 . Moreover, ϖ u + ( 1 1 q u q u 1 ) β j ( k 1 ) + ( 1 1 q u q u 1 ) k δ j 2 I + 1 2 β j v = 1 , v u k G u v G u v T + 1 2 v = 1 k G u v G u v T + β j ( G u u s D u ) + δ j D u < 0 for u = 1 , 2 and l = 1 , 2 , 3 . All the conditions of Corollary 2 are fulfilled. Therefore, it follows Corollary 2 that the network (18) achieves cluster synchronization. Figure 4, Figure 5 and Figure 6 show the time history of the components of the error variables corresponding to the pinning controlled network (18). It can be seen that all the components of the error variables tend to 0 when t + . These indicate that the network (18) achieves cluster synchronization.

5. Conclusions

In this paper, the problem of clustering component synchronization of nonlinearly coupled complex networks with nonidentical nodes and asymmetric couplings is investigated. With the use of matrix analysis and stability theory, some criteria for clustering component synchronization are established. Two numerical simulation examples are also provided to show the effectiveness of the theoretical results obtained.

Author Contributions

Conceptualization, J.-P.S. and J.L.; methodology, J.L.; software, J.L.; validation, J.L.; formal analysis, J.-P.S. and J.L.; investigation, J.L.; data curation, J.L.; writing—original draft preparation, J.L.; writing—review and editing, J.-P.S. and J.L.; visualization, J.L.; supervision, J.-P.S.; project administration, J.-P.S.; funding acquisition, J.-P.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China under Grant 12361039.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Liu, X.; Chen, T. Synchronization analysis for nonlinearly-coupled complex networks with an asymmetrical coupling matrix. Phys. A Stat. Mech. Appl. 2008, 387, 4429–4439. [Google Scholar] [CrossRef]
  2. Zhang, L.; Lei, Y.; Wang, Y.; Chen, H. Generalized outer synchronization between non-dissipatively coupled complex networks with different-dimensional nodes. Appl. Math. Model. 2018, 55, 248–261. [Google Scholar] [CrossRef]
  3. Li, Y.; Zhang, J.; Lu, J.; Lou, J. Finite-time synchronization of complex networks with partial communication channels failure. Inf. Sci. 2023, 634, 539–549. [Google Scholar] [CrossRef]
  4. Abdurahman, A.; Abudusaimaiti, M.; Jiang, H. Fixed/predefined-time lag synchronization of complex-valued BAM neural networks with stochastic perturbations. Appl. Math. Comput. 2023, 444, 127811. [Google Scholar] [CrossRef]
  5. Syed Ali, M.; Stamov, G.; Stamova, I.; Ibrahim, T.F.; Dawood, A.A.; Osman Birkea, F.M. Global asymptotic stability and synchronization of fractional-order reaction–diffusion fuzzy BAM neural networks with distributed delays via hybrid feedback controllers. Mathematics 2023, 11, 4248. [Google Scholar] [CrossRef]
  6. Liang, Y.; Deng, Y.; Zhang, C. Outer synchronization of two muti-layer dynamical complex networks with intermittent pinning control. Neural Netw. 2023, 11, 3543. [Google Scholar] [CrossRef]
  7. Kumar, A.; Das, S.; Singha, S.; Rajeev. Quasi-projective synchronization of inertial complex-valued recurrent neural networks with mixed time-varying delay and mismatched parameters. Chaos Solitons Fractals 2023, 166, 112948. [Google Scholar] [CrossRef]
  8. Wei, C.; Wang, X.; Hui, M.; Zeng, Z. Quasi-synchronization of fractional multiweighted coupled neural networks via aperiodic intermittent control. IEEE Trans. Cybern. 2024, 54, 1671–1684. [Google Scholar] [CrossRef]
  9. Liu, X.; Chen, T. Cluster synchronization in directed networks via intermittent pinning control. IEEE Trans. Neural Netw. 2011, 22, 1009–1020. [Google Scholar] [CrossRef]
  10. Li, K.; Zhou, J.; Yu, W.; Small, M.; Fu, X. Adaptive cluster synchronization in networks with time-varying and distributed coupling delays. Appl. Math. Model. 2014, 38, 1300–1314. [Google Scholar] [CrossRef]
  11. Zhou, L.; Wang, C.; Du, S.; Zhou, L. Cluster synchronization on multiple nonlinearly coupled dynamical subnetworks of complex networks with nonidentical nodes. IEEE Trans. Neural Netw. Learn. Syst. 2017, 28, 570–583. [Google Scholar] [CrossRef]
  12. Liu, P.; Zeng, Z.; Wang, J. Asymptotic and finite-time cluster synchronization of coupled fractional-order neural networks with time delay. IEEE Trans. Neural Netw. Learn. Syst. 2020, 31, 4956–4967. [Google Scholar] [CrossRef]
  13. Della Rossa, F.; Pecora, L.; Blaha, K.; Shirin, A.; Klickstein, I.; Sorrentino, F. Symmetries and cluster synchronization in multilayer networks. Nat. Commun. 2020, 11, 3179. [Google Scholar] [CrossRef]
  14. Wang, Y.; Wu, Z. Cluster synchronization in variable-order fractional community network via intermittent control. Mathematics 2021, 9, 2596. [Google Scholar] [CrossRef]
  15. Wang, X.; Zhai, S.; Luo, G.; Huang, T. Cluster synchronization in a network of nonlinear systems with directed topology and competitive relationships. Appl. Math. Comput. 2022, 421, 126931. [Google Scholar] [CrossRef]
  16. Tong, L.; Liang, J.; Liu, Y. Generalized cluster synchronization of Boolean control networks with delays in both the states and the inputs. J. Frankl. Inst. 2022, 359, 206–223. [Google Scholar] [CrossRef]
  17. Lu, B.; Jiang, H.; Hu, C.; Abdurahman, A.; Liu, M. Adaptive pinning cluster synchronization of a stochastic reaction-diffusion complex network. Neural Netw. 2023, 166, 524–540. [Google Scholar] [CrossRef]
  18. Zhai, S.; Wang, X.; Zheng, Y. Cluster synchronization of a nonlinear network with fixed and switching topologies. IEEE Syst. J. 2023, 17, 3752–3761. [Google Scholar] [CrossRef]
  19. Liu, L.; Zhou, W.; Huang, C. Finite/Prescribed-time cluster synchronization of complex dynamical networks with multiproportional delays and asynchronous switching. IEEE Trans. Syst. Man Cybern. Syst. 2023, 53, 3683–3694. [Google Scholar] [CrossRef]
  20. Zhai, S.; Zheng, W.X. Stability conditions for cluster synchronization in directed networks of diffusively coupled nonlinear systems. IEEE Trans. Circuits Syst. I Regul. Pap. 2023, 70, 413–423. [Google Scholar] [CrossRef]
  21. Lin, L.; Cao, J.; Liu, X.; Lu, G.; Abdel-Aty, M. Cluster synchronization of finite-field networks. IEEE Trans. Cybern. 2024, 54, 635–640. [Google Scholar] [CrossRef]
  22. Hui, M.; Liu, X.; Zhu, S.; Cao, J. Event-triggered impulsive cluster synchronization of coupled reaction–diffusion neural networks and its application to image encryption. Neural Netw. 2024, 170, 46–54. [Google Scholar] [CrossRef]
  23. Wang, X.F.; Chen, G. Pinning control of scale-free dynamical networks. Phys. A Stat. Mech. Appl. 2002, 310, 521–531. [Google Scholar] [CrossRef]
  24. Wu, W.; Zhou, W.; Chen, T. Cluster synchronization of linearly coupled complex networks under pinning control. IEEE Trans. Circuits Syst. I Regul. Pap. 2009, 56, 829–839. [Google Scholar] [CrossRef]
  25. Feng, J.; Wang, J.; Xu, C.; Austin, F. Cluster synchronization of nonlinearly coupled complex networks via pinning control. Discret. Dyn. Nat. Soc. 2011, 2011, 262349. [Google Scholar] [CrossRef]
  26. Wang, Y.; Ma, Z.; Cao, J.; Alsaedi, A.; Alsaadi, F.E. Adaptive cluster synchronization in directed networks with nonidentical nonlinear dynamics. Complexity 2016, 21, 380–387. [Google Scholar] [CrossRef]
  27. Li, L.; Ho, D.W.C.; Cao, J.; Lu, J. Pinning cluster synchronization in an array of coupled neural networks under event-based mechanism. Neural Netw. 2016, 76, 1–12. [Google Scholar] [CrossRef]
  28. Zhang, X.; Zhou, W.; Karimi, H.R.; Sun, Y. Finite- and fixed-time cluster synchronization of nonlinearly coupled delayed neural networks via pinning control. IEEE Trans. Neural Netw. Learn. Syst. 2021, 32, 5222–5231. [Google Scholar] [CrossRef]
  29. Yang, Q.; Wu, H.; Cao, J. Pinning exponential cluster synchronization for fractional-order complex dynamical networks with switching topology and mode-dependent impulses. Neurocomputing 2021, 428, 182–194. [Google Scholar] [CrossRef]
  30. Liu, L.; Ding, X.; Zhou, W. Prescribed-time cluster synchronization of uncertain complex dynamical networks with switching via pinning control. Neurocomputing 2021, 419, 136–147. [Google Scholar] [CrossRef]
  31. Wang, J.; Liu, X. Cluster synchronization for multi-weighted and directed complex networks via pinning control. IEEE Trans. Circuits Syst. II Express Briefs 2022, 69, 1347–1351. [Google Scholar] [CrossRef]
  32. Fan, H.; Shi, K.; Zhao, Y. Pinning impulsive cluster synchronization of uncertain complex dynamical networks with multiple time-varying delays and impulsive effects. Phys. A Stat. Mech. Appl. 2022, 587, 126534. [Google Scholar] [CrossRef]
  33. Wei, J.; Zhang, C.; Guo, Y.; Wang, F. Cluster synchronization of stochastic two-layer delayed neural networks via pinning impulsive control. J. Appl. Math. Comput. 2024, 2024, 1–18. [Google Scholar] [CrossRef]
  34. Chen, Z.; Sun, Y.; Li, J.; Liu, M. Distributed fault-tolerant pinning control for cluster synchronization of multiple unmanned surface vehicles. Ocean Eng. 2024, 296, 116882. [Google Scholar] [CrossRef]
  35. Wu, B.; Ma, Z.; Wang, Y. Partial component consensus of leader-following multi-agent systems. Acta Phys. Sin. 2017, 66, 060201. [Google Scholar] [CrossRef]
  36. Li, F.; Ma, Z.; Duan, Q. Partial component synchronization on chaotic networks. Phys. A Stat. Mech. Appl. 2019, 515, 707–714. [Google Scholar] [CrossRef]
  37. Li, F.; Ma, Z.; Duan, Q. Clustering component synchronization in a class of unconnected networks via pinning control. Phys. A Stat. Mech. Appl. 2019, 525, 394–401. [Google Scholar] [CrossRef]
  38. Hu, W.; Zhang, W.; Ma, Z.; Li, K. Partial component consensus analysis of second-order and third-order nonlinear multi-agent systems. Phys. A Stat. Mech. Appl. 2022, 593, 126857. [Google Scholar] [CrossRef]
  39. Wang, Y.; Song, H.; Chen, G.; Ma, Z.; Cao, J. p components of cluster-lag consensus for second-order multiagent systems with adaptive controller on cooperative-competitive networks. IEEE Trans. Cybern. 2023, 53, 2852–2863. [Google Scholar] [CrossRef]
  40. Liao, X. Mathematical Theory of Stability and Its Application; Central China Normal University Press: Wuhan, China, 2001. [Google Scholar]
  41. Wu, J.; Jiao, L. Synchronization in complex delayed dynamical networks with nonsymmetric coupling. Phys. A Stat. Mech. Appl. 2007, 386, 513–530. [Google Scholar] [CrossRef]
  42. Wu, C.W.; Chua, L.O. Synchronization in an array of linearly coupled dynamical systems. IEEE Trans. Circuits Syst. I Fundam. Theory Appl. 1995, 42, 430–447. [Google Scholar] [CrossRef]
  43. Lütkepohl, H. Handbook of Matrices; Wiley: New York, NY, USA, 1996. [Google Scholar]
Figure 1. The time history of e 1 i ( t ) , i = 1 , 2 , , 9 .
Figure 1. The time history of e 1 i ( t ) , i = 1 , 2 , , 9 .
Mathematics 12 01022 g001
Figure 2. The time history of e 2 i ( t ) , i = 1 , 2 , , 9 .
Figure 2. The time history of e 2 i ( t ) , i = 1 , 2 , , 9 .
Mathematics 12 01022 g002
Figure 3. The time history of e 3 i ( t ) , i = 1 , 2 , , 9 .
Figure 3. The time history of e 3 i ( t ) , i = 1 , 2 , , 9 .
Mathematics 12 01022 g003
Figure 4. The time history of e 1 i ( t ) , i = 1 , 2 , , 10 .
Figure 4. The time history of e 1 i ( t ) , i = 1 , 2 , , 10 .
Mathematics 12 01022 g004
Figure 5. The time history of e 2 i ( t ) , i = 1 , 2 , , 10 .
Figure 5. The time history of e 2 i ( t ) , i = 1 , 2 , , 10 .
Mathematics 12 01022 g005
Figure 6. The time history of e 3 i ( t ) , i = 1 , 2 , , 10 .
Figure 6. The time history of e 3 i ( t ) , i = 1 , 2 , , 10 .
Mathematics 12 01022 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, J.; Sun, J.-P. Clustering Component Synchronization of Nonlinearly Coupled Complex Networks via Pinning Control. Mathematics 2024, 12, 1022. https://doi.org/10.3390/math12071022

AMA Style

Liu J, Sun J-P. Clustering Component Synchronization of Nonlinearly Coupled Complex Networks via Pinning Control. Mathematics. 2024; 12(7):1022. https://doi.org/10.3390/math12071022

Chicago/Turabian Style

Liu, Jie, and Jian-Ping Sun. 2024. "Clustering Component Synchronization of Nonlinearly Coupled Complex Networks via Pinning Control" Mathematics 12, no. 7: 1022. https://doi.org/10.3390/math12071022

APA Style

Liu, J., & Sun, J. -P. (2024). Clustering Component Synchronization of Nonlinearly Coupled Complex Networks via Pinning Control. Mathematics, 12(7), 1022. https://doi.org/10.3390/math12071022

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop