Identifying the interactions in a colored dynamical network*
Wu Zhao-Yan, Gong Xiao-Li
College of Mathematics and Information Science, Jiangxi Normal University, Nanchang 330022, China

Corresponding author. E-mail: zhywu7981@gmail.com

*Project supported by the National Natural Science Foundation of China (Grant No. 61463022), the Natural Science Foundation of Jiangxi Educational Committee, China (Grant No. GJJ14273), and the Graduate Innovation Fund of Jiangxi Normal University, China (Grant No. YJS2014061).

Abstract

The interactions of a colored dynamical network play a great role in its dynamical behaviour and are denoted by outer and inner coupling matrices. In this paper, the outer and inner coupling matrices are assumed to be unknown and need to be identified. A corresponding network estimator is designed for identifying the unknown interactions by adopting proper adaptive laws. Based on the Lyapunov function method and Barbalat’s lemma, the obtained result is analytically proved. A colored network coupled with chaotic Lorenz, Chen, and Lü systems is considered as a numerical example to illustrate the effectiveness of the proposed method.

Keyword: 02.40.Pc; 64.60.aq; 05.45.Xt; colored network; interaction; identification
1. Introduction

Many large-scale physical systems consisting of interactive individuals are usually modelled by complex dynamical networks, in which the nodes denote the individuals and the edges denote the interactions among them.[110] Generally, the interactions are denoted by the outer and inner coupling matrices. The outer coupling matrix denotes the topology and the coupling strength of the network, from which one can know whether a pair of nodes affect each other. The inner coupling matrix denotes the interactive manner between a pair of nodes. Take a network coupled with a chaotic Lorenz system as an example, if the first, second, and third components of a pair of nodes affect each other correspondingly, the inner coupling matrix between them is an identity matrix. Due to the complexity of real dynamical networks, the node dynamics and the interactions may be different. Take a dynamical network coupled with chaotic Lorenz and Chen systems as an example, one pair of nodes interact only in the first component, and another pair of nodes interact only in the last component. For better describing this kind of network, namely a colored graph in mathematics, a colored network consisting of colored nodes and colored edges is proposed.[9, 10] Two nodes with different colors have different node dynamics, and two edges with different colors denote the two pairs of nodes have different interactions.

As we know, interactions of a network play a great role in its dynamical behaviour. For example, a dynamical network with different outer and inner coupling matrices can display a totally different synchronization region, [11] even discontinuous synchronization regions.[12] Moreover, the literature has studied the relationship between the network synchronizability and network structural parameters.[13, 14] However, in many real networks, the interactions are not exactly known beforehand. Therefore, how to identify the interactions in a network becomes a hot topic and has been widely studied.[1524] Yu et al. suggested a method for estimating the topology of a network for the first time.[15] Wu et al. identified the topologies of general weighted networks with time-varying coupling delays.[16] Mei et al. considered the finite-time structure identification of a dynamical network.[22] In all the above-mentioned literature, the inner coupling matrix is assumed to be known, i.e., only the outer coupling matrix is unknown and to be identified. Naturally, if the outer and inner coupling matrix are both unknown, how to identify them is a challenging and important problem. In particular, the inner matrices of a colored network may be totally different, and how to identify them deserves further studies.

Motivated by the above discussion, this paper considers the identification of interactions in a colored network with the same dimension node dynamics. By introducing proper adaptive laws, an effective network estimator is designed for identifying the unknown interactions. According to the Lyapunov function method and Barbalat’ s lemma, the proof is completed. Section 2 introduces the network model and some preliminaries. Section 3 gives the main result. Section 4 performs a numerical simulation to verify the effectiveness of the derived theoretical result. Finally, Section 5 concludes the paper.

2. Model description and preliminaries

Consider a colored network consisting of N nodes, which can be described by

where xi(t) = (xi1(t), xi2(t), … , xin(t))TRn is the state variable of node i, fi : RnRn is a nonlinear vector-valued function and denotes the local dynamics. Matrix A = (aij) ∈ RN× N is the outer coupling matrix and denotes the coupling strength and topology. If there is a connection from node j to node i (ij), then aij ≠ 0; otherwise, aij = 0. The diagonal entries of A are zero, i.e., aii = 0, i = 1, 2, … , N. Matrix is the inner coupling matrix. If the change rate of the k-th component of node i (i.e., ik(t)) is affected by the difference between the l-th components of nodes i and j (i.e., xjl(t) − xil(t)), then ; otherwise, .

In this paper, both the outer and inner coupling matrices are assumed to be unknown and need to be identified. Let and , then the problem of identifying aij and Hij becomes the problem of identifying Cij.

Define . Then, the network (1) can be written as

For identifying the matrices Cij in network (2), we design the following network estimator:

where i, j = 1, 2, … , N, k, l = 1, 2, … , n, yi(t) = (yi1(t), yi2(t), … , yin(t))TRn is the state variable of node i in the network estimator, , is the estimation of is the adaptive feedback controller applied on node i, d(t) is the adaptive feedback gain, η and δ are positive constants.

Let ei(t) = yi(t) − xi(t) be the synchronization error, we have the following error system:

Assumption 1 Suppose that there exists a positive L such that (yx)T (fi(y) − fi(x)) ≤ L(yx)T(yx), i = 1, 2, … , N, holds for any x, yRn.

Assumption 2 Suppose that are linearly independent for t > 0.

Barbalat’ s Lemma[25] If Φ : RR+ is a uniformly positive function for t ≥ 0 and if the limit of the integral exists and is finite, then limt→ ∞ Φ (t) = 0.

3. Main result

Theorem 1 Suppose that Assumptions 1 and 2 hold. Then the unknown matrices of the colored network (2) can be identified by the estimated matrices of the network estimator (3).

Proof Consider the following Lyapunov function candidate:

where d* is an arbitrary positive constant to be determined.

Then the derivative of V(t) with respect to t along the trajectories of Eq. (4) is

where , INn is Nn-dimensional identity matrix, and the matrix C is defined as

Let λ C be the largest eigenvalue of matrix (C + CT)/2, then one can choose d* = L + λ C + 1 such that

Then, one has

and supt≥ 0V(t) ≤ V(0), i.e., V(t) is bounded. That is to say, e(t) ∈ L2 and e(t) are bounded, i.e., E(t) ∈ L. From error system (4), Ė (t) exists and is bounded for t ∈ [0, + ∞ ). According to Barbalat’ s lemma, we obtain limt→ ∞ e(t) = 0, followed by limt→ ∞ ė (t) = 0. Then, as t → ∞ , converges to the following largest invariant set:

which is equivalent to

From Assumption 2, we have . That is, the coupling matrices of network (1) are identified when the synchronization between networks (2) and (3) is achieved. Thus, the proof is completed.

Remark 1 In colored network (1), the inner coupling matrices Hij need not be diagonal. If the inner coupling matrices Hij are diagonal, i.e., for kl, then the adaptive laws about in (3) can be replaced by for kl and , k = 1, 2, … , n, i, j = 1, 2, … , N.

4. Numerical simulations

In this section, a simple colored network consisting of 3 colored nodes and 3 colored edges is considered in a numerical example to verify the effectiveness of the obtained result, which is shown in Fig. 1. The first, second, and third node dynamics are chosen as the Lorenz system[26]

Chen system[27]

and Lü system[28]

Fig. 1. A simple colored network consisting of 3 colored nodes and 3 colored edges. The red, green, and blue nodes are chosen as the Lorenz, Chen, and Lü systems, respectively.

The outer and inner coupling matrices in Fig. 1 are chosen as follows:

In numerical simulations, choose η = 5, δ = 5, the initial values , and d(0) = 10. Moreover, choose the initial values of xi(t) and yi(t) randomly. From Fig. 2, one can find , and converge to 0.2, and others converge to 0 as t → ∞ , which yields a12 = 0.2, . That is, the interactions between nodes 1 and 2 are identified. Similarly, from Figs. 3 and 4, we have a13 = 0.1, a23 = 0.1, , and . That is, all the interactions in the colored network are identified.

Fig. 2. The orbits of , k, l = 1, 2, 3. Clearly, , and converge to 0.2 and others converge to 0 as t → ∞ , which yields a12 = 0.2 and . That is, the interactions between nodes 1 and 2 are identified.

Remark 2 In the above numerical example, for simplicity, a colored network consisting of 3 colored nodes and 3 colored edges is considered. From the network model (1) and the proof of Theorem 1, a colored network considered in this paper consists of N nodes. That is, the proposed method in this paper can handle the networks with more than 3 nodes.

Fig. 3. The orbits of , k, l = 1, 2, 3. Clearly, , , converge to 0.1 and others converge to 0 as t → ∞ , which yields a13 = 0.1, . That is, the interactions between nodes 1 and 3 are identified.
Fig. 4. The orbits of , k, l = 1, 2, 3. Clearly, , , and converge to 0.3 and others converge to 0 as t → ∞ , which gives a23 = 0.1, . That is, the interactions between nodes 2 and 3 are identified.

5. Conclusions and discussion

In this paper, the identification problem of a colored dynamical network is considered. Compared with the existing results about the identification problem of dynamical networks, not only the outer coupling matrix but also the inner coupling matrices are unknown and need to be identified. The inner coupling matrices can be totally different. For identifying the unknown interactions (outer and inner coupling matrices), an effective network estimator with proper adaptive laws is designed. According to the Lyapunov function method and Barbalat’ s lemma, the result is rigorously proved. Finally, a numerical example is provided to illustrate the effectiveness of the obtained result.

Though the node dynamics of the colored network considered in this paper can be different, they have the same dimension. In real networks, the node dynamics may have different dimensions. How to design an effective network estimator for identifying the interactions of a colored network with different-dimensional node dynamics is a challenging problem and deserves further studies. On the other hand, if the node dynamics of a colored network are also unknown, how to extend the topology identification methods from observed dynamical time series, such as Granger causality[20] and compressive sensing[8, 24] method, to identify the unknown interactions of a colored network will be considered in our future work.

Reference
1 Albert R and Barabási A L 2002 Rev. Modern Phys. 74 47 DOI:10.1103/RevModPhys.74.47 [Cited within:1]
2 Arenas A, Díaz-Guilera A, Kurths J, Moreno Y and Zhou C S 2008 Phys. Rep. 469 93 DOI:10.1016/j.physrep.2008.09.002 [Cited within:1]
3 Boccaletti S, Latora V, Moreno Y, Chavezf M and Hwang D U 2006 Phys. Rep. 424 175 DOI:10.1016/j.physrep.2005.10.009 [Cited within:1]
4 Li K Z, He E, Zeng Z R and Tse C K 2013 Chin. Phys. B 22 070504 DOI:10.1088/1674-1056/22/7/070504 [Cited within:1]
5 Zhang X Y, Hu X, Kurths J and Liu Z H 2013 Phys. Rev. E 88 010802 DOI:10.1103/PhysRevE.88.010802 [Cited within:1]
6 Wang S G and Yao H X 2012 Chin. Phys. B 21 050508 DOI:10.1088/1674-1056/21/5/050508 [Cited within:1]
7 Liu D F, Wu Z Y and Ye Q L 2014 Nonlinear Dyn. 75 209 DOI:10.1007/s11071-013-1059-5 [Cited within:1]
8 Wang W X, Yang R, Lai Y C, Kovanis V and Harrison M A F 2011 Europhys. Lett. 94 48006 DOI:10.1209/0295-5075/94/48006 [Cited within:1]
9 Wu Z Y, Xu X J, Chen G R and Fu X C 2012 Chaos 22 043137 DOI:10.1063/1.4769991 [Cited within:1]
10 Sun M, Li D D, Han D and Jia Q 2013 Chin. Phys. Lett. 30 090202 DOI:10.1088/0256-307X/30/9/090202 [Cited within:2]
11 Duan Z S, Chen G R and Huang L 2008 Phys. Lett. A 372 3741 DOI:10.1016/j.physleta.2008.02.056 [Cited within:1]
12 Duan Z S, Chen G R and Huang L 2009 IEEE Trans. Auto. Control 54 845 DOI:10.1109/TAC.2008.2009690 [Cited within:1]
13 Dadashi M, Barjasteh I and Jalili M 2010 Chaos 20 043119 DOI:10.1063/1.3515170 [Cited within:1]
14 Shi D H, Chen G R, Thong W W K and Yan X Y 2013 IEEE Circ. Syst. Mag. 13 66 DOI:10.1109/MCAS.2012.2237145 [Cited within:1]
15 Yu D C, Righero M and Kocarev L 2006 Phys. Rev. Lett. 97 188701 DOI:10.1103/PhysRevLett.97.188701 [Cited within:2]
16 Wu X Q 2008 Physica A 387 997 DOI:10.1016/j.physa.2007.10.030 [Cited within:1]
17 Wu Z Y and Fu X C 2009 Chin. Phys. Lett. 26 070201 DOI:10.1088/0256-307X/26/7/070201 [Cited within:1]
18 Liu H, Lu J A, J H and Hill D J 2009 Automatica 45 1799 DOI:10.1016/j.automatica.2009.03.022 [Cited within:1]
19 Xu Y H, Zhou W N, Fang J A and Sun W 2011 Commun. Nonlinear Sci. Numer. simulat. 16 3337 DOI:10.1016/j.cnsns.2010.11.003 [Cited within:1]
20 Wu X Q, Wang W H and Zheng W X 2012 Phys. Rev. E 86 046106 DOI:10.1103/PhysRevE.86.046106 [Cited within:1]
21 Wu Z Y and Fu X C 2013 IET Contr. Theory Appl. 7 1269 DOI:10.1049/iet-cta.2013.0201 [Cited within:1]
22 Mei J, Jiang M H and Wang J 2013 Commun. Nonlinear Sci. Numer. simulat. 18 999 DOI:10.1016/j.cnsns.2012.08.039 [Cited within:1]
23 Liu D F, Wu Z Y and Ye Q L 2014 Chin. Phys. B 23 040504 DOI:10.1088/1674-1056/23/4/040504 [Cited within:1]
24 Li G J, Wu X Q, Liu J, Lu J A and Guo C 2015 Chaos 25 043102 DOI:10.1063/1.4916788 [Cited within:2]
25 Khalil H K 2002 Nonlinear Systems Englewood Cliffs, NJ Prentice-Hall [Cited within:1]
26 Lorenz E N 1963 J. Atmos. Sci. 20 130 DOI:10.1175/1520-0469(1963)020<0130:DNF>2.0.CO;2 [Cited within:1]
27 Chen G R and Ueta T 1999 Int. J. Bifurc. Chaos 9 1465 DOI:10.1142/S0218127499001024 [Cited within:1]
28 J H and Chen G R 2002 Int. J. Bifurc. Chaos 12 659 DOI:10.1142/S0218127402004620 [Cited within:1]