Research Article Synchronization Analysis of Two Coupled Complex Networks with Time Delays

Transcription

1 Discrete Dynamics in ature and Society Volume, Article ID 9, pages doi:.55//9 Research Article Synchronization Analysis of Two Coupled Complex etworks with Time Delays Weigang Sun, Jingyuan Zhang, and Changpin Li School of Science, Hangzhou Dianzi University, Hangzhou 8, China Department of Mathematics, Shanghai University, Shanghai, China Correspondence should be addressed to Weigang Sun, Received 7 August ; Revised January ; Accepted January Academic Editor: Recai Kilic Copyright q Weigang Sun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. This paper studies the synchronized motions between two complex networks with time delays, which include individual inner synchronization in each network and outer synchronization between two networks. Based on the Lyapunov stability theory and the linear matrix equality LMI, a synchronous criterion for inner synchronization inside each network is derived. umerical examples are given which fit the theoretical analysis. In addition, the involved numerical results show that the delays between two networks have little effect on inner synchronization. It is also shown that synchronous motions within each network or between two networks are not enhanced if individual intranetwork connections are allowed.. Introduction We refer to network synchronization, that is, synchronizing all the nodes inside a network, as inner synchronization, which is widely studied in the science world based on the appearance of small-world and scale-free network models. Afterwards, the improved and expanded work in this respect that is, introducing weighted connections, time varying coupling matrices, nonlinear coupling function, time delays, and so forth can be found in the literature and many references cited therein. The main approaches on studying synchronization inside a network is decoupling network systems, and studying the lowdimensional systems through the master stability function, or using the tool in equality LMI toolbox. Recently outer synchronization was proposed, which aimed at studying synchronization between coupled networks. In, the authors theoretically and numerically demonstrated the possibility of synchronization between two networks. By the open-plusclosed-loop OPCL method, synchronization between two networks can be realized with the same topological structures. In reality, if network nodes are of similar properties

2 Discrete Dynamics in ature and Society same node dynamics, we regard it as one network, otherwise, as more networks. Several examples on two networks are presented in the literature,. Actions between two networks are colorful, such as through the nonlinear signals, bipartite connection, or by special nodes. Recently, Sorrentino and Ott studied network synchronization of groups and discussed the feasibility of inner synchronization of individual network. The problem of collective behavior in a network or between two networks is of broad interest. As a first example, in subway systems e.g., Shanghai, China, when the trains arrive at the platform, the outer and inner doors simultaneously open or close, which shows that both inner synchronization and outer synchronization happen. Another example is two unmanned vehicles are assigned to accomplish independent tasks, such as cooperative searches and attacks, which focuses on inner synchronization. If the communicated information between two networks is not amount to lead to desynchronization, we should adopt control skills 5 7 to make synchronization happen between them, which has received increasing attention; see 8. It is also found that present studies on synchronization inside a network and between two networks simultaneously are much less. Time delays commonly exist in natural systems, since the finite speed of signal transmission over a distance gives rise to finite time delay, which plays an important role in the stability of synchronization. For the two coupled networks, the signals from one network to another often create delays due to the distance. Therefore, motivated by the impact of topological structures and the delays on the dynamics of the networks, this paper mainly focuses on the effect of delays on inner synchronization inside each network and outer synchronization between two coupled networks. In Section, network models and synchronization analysis are presented, and numerical examples are shown in Section, including two examples for the same or different dimension of node dynamics. Finally, the discussions are included in the last section.. Model Presentation and Synchronization Analysis Consider the following network models: the node dynamics in network X is ẋ i f x i t, i,..., x,andẏ j g y j t, j,..., y, for the nodes in network Y, where f x : R n x R n x, g y : R n y R n y are continuously differential functions and x i y j is an n x -dimensional n y -dimensional state vector. The dynamical equations of the network systems with time delays are as follows: ( ) ẋ i t f x i t A ij Γ y j t τy, i,,...,x, j ẏ j t g ( y j t ) x B ji Γ x i t τ x, j,,..., y, i. where A is an x y dimensional coupling matrix, whose entries A ij represent the intensity of the direct interaction from i in network X to j in network Y, analogously the entries of B ji.matrixγ Γ R n x n y R n y n x is the inner-coupling matrix. τ x,τ y are the time delays between networks. The action sketch between two networks with time delays is shown in Figure.

3 Discrete Dynamics in ature and Society x τ x y x y x τ y y etwork X etwork Y Figure : A simple action sketch between two coupled networks with time delays. For the more general models, we allow connections in each network, the networked system then reads as ( ) x ẋ i t f x i t A ij Γ y j t τy C im Γ x m t, i,,..., x, j m ẏ j t g ( y j t ) x B ji Γ x i t τ x D jk Γ y k t, j,,..., y, i k. where C C im R x x, D D jk R y y are the coupling matrices, satisfying the sum of every row being zero. Γ R n x n x, Γ R n y n y are also inner-coupling matrices. Let us now consider the possibility that the individual networks achieve synchronization; that is, x t x x t x s t and y t y y t y s t. If there exist such synchronous states, satisfying A ij μ, i X, j x B ji μ, j Y,. i without loss of generality, we assume that μ μ. Thus, the synchronized state equations are ẋ s t f x s t Γ y s ( t τy ), ẏ s t g ( y s t ) Γ x s t τ x.. Linearizing the synchronous states around x s and y s,weget ( ) δẋ i t J t δx i t A ij Γ δy j t τy, i,,...,x, j x δẏ j t W t δy j t B ji Γ δx i t τ x, j,,..., y, i.5

4 Discrete Dynamics in ature and Society where J t Df x s t and W t Dg y s t are the Jacobians of f x t and g y t at x s and y s, respectively. We assume that the x y -independent solutions of.5 can be expressed in the form δx i Φ xi δx, i,,..., x, δy j Φ yj δy, j,,..., y, where {Φ xi }, {Φ yj } are the suitable time-independent scalars and δx, δy are the appropriate variables. If the dimension of the space vectors given by the values of Φ xi, Φ yj i,,..., x,j,,..., y is x y, we can see that this assumed form includes all possible linear solutions of.5. Substituting this assumption into.5, they become Φ xi δẋ t Φ xi J t δx t A ij Φ yj Γ δy ( ) t τ y, i,,...,x,.6 j x Φ yj δẏ t Φ yj W t δy t B ji Φ xi Γ δx t τ x, j,,..., y..7 i Thus, in order that.6 resp.,.7 is satisfied for all i resp., j, we require that j A ijφ yj η, where η is independent of i, andφ y j i B jiφ xi η, where η is Φ x i independent of j. SetΦ x Φ x,...,φ xx and Φ y Φ y,...,φ yy. From this, we get AΦ y η Φ x, BΦ x η Φ y ;thatis, ( A B )( Φx Φ y ) ( ) η Φ x..8 η Φ y Substituting.8 in.6 and.7,weobtain δẋ t J t δx t η Γ δy ( t τ y ), δẏ t W t δy t η Γ δx t τ x..9 One particular solution of.8 is derived when η η λ, then ( Φ ) ( x Φ ) ( ) x A Θ λ, Θ,. B Φ y Φ y where λ is the possibly complex eigenvalues of the matrix Θ. Set η λξ, η λ/ξ, Φ x Φ x, Φ y ξφ y, where ξ is a free parameter. Equation. becomes ( )( A Φ ) λξ Φ x x ( ) λ,. B ξφ y ξφ y ξ

5 Discrete Dynamics in ature and Society 5 which shows that the solution of. contains all the possible solutions of.8. We rewrite.9 as δẋ t J t δx t λγ δy ( t τ y ), δẏ t W t δy t λγ δx t τ x,. where λ λ,λ,...,λ x y. In, the authors gave the explicit analysis on the spectrum of Θ, and two forms of the constructed matrix Θ were shown. Here, we adopt one that matrix Θ has real eigenvalues and only study network.. In the sequel, we utilize the LMI to derive a synchronous theorem. Theorem.. Consider network model.. If there exist two positive matrices P, Q >, satisfying Ψ λpγ Ψ λqγ Ξ λγ T Q I <,. n x λγ T P I n y where Ψ J t T P PJ t I nx, Ψ W t T Q QW t I ny, then the network. asymptotically synchronizes to x s,y s defined by. for the fixed delays τ x,τ y >, respectively. Proof. Consider.. Choose a Lyapunov-Krasovskii functional as t t V t δx T t Pδx t δy T t Qδy t δx T α δx α dα t τ x δy T α δy α dα. t τ y. Therefore, V t δẋ T t Pδx t δx T t Pδẋ t δẏ T t Qδy t δy T t Qδẏ t δx T t δx t δx T t τ x δx t τ x δy T t δy t δy T( ) ( ) t τ y δy t τy δx T t Ψ δx t δy T( ) t τ y λγ T Pδx t δyt t Ψ δy t δx T t τ x λγ T Qδy t δy T t λqγ δx T t τ x δx T t λpγ δy T( ) t τ y δx T t τ x δx t τ x δy T( t τ y ) δy ( t τy ) T δx t δx t δy t δy t Ξ. δx t τ x δx t τ x δy t τ y δy ( ) t τ y.5

6 6 Discrete Dynamics in ature and Society According to the known condition Ξ < and Lyapunov-Krasovskii stability Theorem 6,the zero solutions of. are asymptotically stable, which shows that synchronization inside network X and network Y happens when t tends to. Remark.. Matrix Ξ relies upon t and τ x,τ y. The condition can be changed into Ξ < when t>tfor enough big T>. The linear matrix equality method is inadequate for assessing the stability of the synchronous solution when both intragroup and extragroup connections are allowed in the network... umerical Examples In this section, we will give two examples to illustrate our obtained results, which includes two cases: n x / n y and n x n y. For the case n x / n y, we only consider the inner synchronization inside each network. On the other hand, n x n y, we discuss the inner synchronization within each network and investigate the outer synchronization between two networks. Example.. Consider the coupled network dynamical equations below, which are in the form. : ( ) ẋ i t a x i t s x i t A ij y j t τy, i,...,x, j x ẏ j t y j t y j t B ji x i t τ x, i. ẏ j t by j t, j,,..., y, where s x m x / m m x x. Here, we take a.6, b 6., m 8/7, m 5/7, for numerical simulation. We introduce two quantities E x max i,...,x x i x s and E y max j,...,y y j y s to measure the extent to which synchronization is achieved. The initial values are randomly chosen in,. The intercoupling matrices are A, B.. When a.6, b 6., m 8/7, m 5/7, τ x τ y.8 anda, B; byusingthe Matlab LMI Toolbox, we solve the LMI. for P>, Q>, and obtain P.587 ; ( ) Q,

7 Discrete Dynamics in ature and Society Ex,Ey..5 t E x E y Figure : Inner synchronization evolution curves with τ x τ y Ex Ey.5 Ex Ey τ x τ y a b Figure : The plot shows E x E y at t with regard to τ x or τ y, a τ y.8; b τ x.8. which indicates that inner synchronization inside each network appears. The inner synchronization evolution curves are depicted in Figure. WeplotsE x E y with different values of τ x or τ y in Figure, which shows that the delays τ x,τ y have little influence on inner synchronization with such intercoupling matrices A, B. ote that the order of magnitude for E x E y is 8. ext, we change the connection topology of A, B and consider a special case globally connected, the elements of A ij / y, B ji / x, i,..., x,andj,..., y with x 5, y. The other parameters are the same as in Figure. The numerical results are given in Figures and 5. It is shown that the node numbers and connection topology influence the inner synchronization less.

8 8 Discrete Dynamics in ature and Society Ex,Ey t E x E y Figure : Inner synchronization evolution curves with globally connected topology, where x 5, y, and τ x τ y Ex Ey τ y Figure 5: The curves of E x E y at t versus τ y,withτ x.8, x 5, and y. Example.. In this example, we let n x n y and x y. Therefore, we discuss inner synchronization inside network X or Y and study synchronization between two networks, that is, outer synchronization proposed in, denoting the quantity E outer max i,..., x i y i, to demonstrate that outer synchronization happens. Consider a one-order dynamical system as the dynamical nodes of the complex networks which is described by ẋ i t x i t ρ ( ) x i t r x A ij y j t τy, j i,,...,,

9 Discrete Dynamics in ature and Society Ex,Ey,Eouter t E x E y E outer Figure 6: Inner and outer synchronization evolution curves with r x.5, r y., τ x., τ y.5, ρ., ρ.6, and 5. ẏ j t y j t ρ y j t r y B ji x i t τ x, i j,,...,.. The elements of A ij /, B ji /, i, j,...,,and 5, when ρ., ρ.6, τ x., and τ y.5, by solving the LMI. for P>, Q>, we obtain P.96, Q.98,.5 satisfying the conditions of a synchronous criterion, then the individual network achieves its own inner synchronized state. Figure 6 plots the inner and outer synchronization evolution curves. In addition, we plot the curves of quantities E x,e y,e outer with regard to r x for r y. in Figure 7. We find that the individual network easily achieves inner synchronization for the arbitrary values of r x, while only when r x., outer synchronization is realized. In the following, we numerically discuss the delay effect on inner and outer synchronization. Fix the values of ρ., ρ.6, r y., r x., and τ y.5, and let the value of τ x vary. We observe that the delay does not influence inner synchronization, and when τ x >.5, outer synchronization disappears. Figure 8 plots the curves of E x,e y,and E outer with regard to τ x. ow, we add the individual connections in network X and network Y, andthe following network systems are written as: ẋ i t x i t ρ ( ) x i t r x A ij y j t τy C im x m t, j m ẏ j t y j t ρ y j t r y B ji x i t τ x D jk y k t, i k i,,...,, j,,...,,.6

10 Discrete Dynamics in ature and Society Ex,Ey,Eouter r x E x E y E outer Figure 7: The curves of E x,e y,ande outer versus r x for r y.. The other parameters are τ x., τ y.5, ρ., and ρ Ex,Ey,Eouter τ x E x E y E outer Figure 8: The plot show E x,e y,ande outer versus τ x with τ y.5, r y., r x., ρ., and ρ.6. where C {C im }, D {D jk } are the Laplacian matrices: m C im for all i and k D jk for all j; obviously, the diffusive couplings m C im, k D jk are null in the synchronized manifold. The parameters are chosen as the same as those in Figure 8, and individual couplings C, D are taken as two random matrices. umerical results show that the adding individual connections have little effect on inner and outer synchronization.

11 Discrete Dynamics in ature and Society. Conclusions In conclusion, we have studied inner synchronization inside networks X and Y and outer synchronization between them with time delays. A synchronous criterion for the occurrence of inner synchronization has been derived. umerical examples show that the delays have little effect on inner synchronization, how to estimate the domain of delays τ x and τ y is an interesting work. It is noted that in Example., outer synchronization only happens in the coupling strength r x., and perhaps the node dynamics in networks X and Y plays a central role, or the bidirectional delayed coupling is weak. Because of the various kinds of actions between two networks, how to derive the criteria on inner and outer synchronization simultaneously is an technical challenging. The theoretical understanding is helpful to study the synchronization between two or more neural networks with appropriate couplings. We hope that such work will appear elsewhere. Acknowledgments This work was supported by the ational atural Science Foundation of China nos. 879, 56, and 67 and Research Foundation of Hangzhou Dianzi University no. KYF756. References D. J. Watts and S. H. Strogatz, Collective dynamics of small-world networks, ature, vol. 9, no. 668, pp., 998. A.-L. Barabási and R. Albert, Emergence of scaling in random networks, Science, vol. 86, no. 59, pp. 59 5, 999. C. Zhou, A. E. Motter, and J. Kurths, Universality in the synchronization of weighted random networks, Physical Review Letters, vol. 96, no., Article ID, 6. A. E. Motter, C. S. Zhou, and J. Kurths, Enhancing complex-network synchronization, Europhysics Letters, vol. 69, no., pp., 5. 5 J. Lü, X. Yu, and G. Chen, Chaos synchronization of general complex dynamical networks, Physica A, vol., no. -, pp. 8,. 6 Q. Wang, Z. Duan, G. Chen, and Z. Feng, Synchronization in a class of weighted complex networks with coupling delays, Physica A, vol. 87, no., pp , 8. 7 C. P. Li, W. Sun, and D. L. Xu, Synchronization of complex dynamical networks with nonlinear inner-coupling functions and time delays, Progress of Theoretical Physics, vol., no., pp , 5. 8 S. M. Cai, J. Zhou, L. Xiang, and Z. Liu, Robust impulsive synchronization of complex delayed dynamical networks, Physics Letters. Section A, vol. 7, no., pp , 8. 9 W. Kinzel, A. Englert, G. Reents, M. Zigzag, and I. Kanter, Synchronization of networks of chaotic units with time-delayed couplings, Physical Review E, vol. 79, no. 5, Article ID 567, 9. Y. Sun and J. Ruan, Synchronization in coupled time-delayed systems with parameter mismatch and noise perturbation, Chaos, vol. 9, no., Article ID, 9. C. P. Li, W. Sun, and J. Kurths, Synchronization between two coupled complex networks, Physical Review E, vol. 76, no., Article ID 6, 7. I. Grosu, R. Banerjee, P. K. Roy, and S. K. Dana, Design of coupling for synchronization of chaotic oscillators, Physical Review E, vol. 8, no., Article ID 6, 9. W. Sun, R. Wang, W. Wang, and J. Cao, Analyzing inner and outer synchronization between two coupled discrete-time networks with time delays, Cognitive eurodynamics, vol., no., pp. 5,. F. Sorrentino and E. Ott, etwork synchronization of groups, Physical Review E, vol. 76, no. 5, Article ID 56, p., 7. 5 Z. S. Duan, J. Wang, G. Chen, and L. Huang, Stability analysis and decentralized control of a class of complex dynamical networks, Automatica, vol., no., pp. 8 5, 8.

12 Discrete Dynamics in ature and Society 6 F. Li and J. Sun, Asymptotic stability of a genetic network under impulsive control, Physics Letters. Section A, vol. 7, no. -, pp. 77 8,. 7 J. Cao and X. Li, Stability in delayed Cohen-Grossberg neural networks: LMI optimization approach, Physica D, vol., no. -, pp. 5 65, 5. 8 H. Tang, L. Chen, J.-A. Lu, and C. K. Tse, Adaptive synchronization between two complex networks with nonidentical topological structures, Physica A, vol. 87, no., pp , 8. 9 M. Sun, C.-Y. Zeng, and L.-X. Tian, Generalized projective synchronization between two complex networks with time-varying coupling delay, Chinese Physics Letters, vol. 6, no., Article ID 5, 9. S. Zheng, Q. Bi, and G. Cai, Adaptive projective synchronization in complex networks with timevarying coupling delay, Physics Letters. Section A, vol. 7, no. 7, pp , 9. X. Wu, W. X. Zheng, and J. Zhou, Generalized outer synchronization between complex dynamical networks, Chaos, vol. 9, no., Article ID 9, p. 9, 9. V. Flunkert, S. Yanchuk, T. Dahms, and E. Schöll, Synchronizing distant nodes: a universal classification of networks, Physical Review Letters, vol. 5, no. 5, Article ID 5,. W.-G. Sun, C.-X. Xu, C.-P. Li, and J.-Q. Fang, Synchronization and bifurcation of general complex dynamical networks, Communications in Theoretical Physics, vol. 7, no. 6, pp. 7 75, 7. J. L. Liang, Z. Wang, Y. Liu, and X. Liu, Robust synchronization of an array of coupled stochastic discrete-time delayed neural networks, IEEE Transactions on eural etworks, vol. 9, no., pp. 9 9, 8.

13 Advances in Operations Research Advances in Decision Sciences Applied Mathematics Algebra Probability and Statistics The Scientific World Journal International Differential Equations Submit your manuscripts at International Advances in Combinatorics Mathematical Physics Complex Analysis International Mathematics and Mathematical Sciences Mathematical Problems in Engineering Mathematics Discrete Mathematics Discrete Dynamics in ature and Society Function Spaces Abstract and Applied Analysis International Stochastic Analysis Optimization