SciELO - Scientific Electronic Library Online

vol.21 número3Numerical Investigation of the Droplet Behavior in Cascades Using a Finite Volume MethodPerformance Analysis of Simulated Annealing Cooling Schedules in the Context of Dense Image Matching índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados




Links relacionados

  • Não possue artigos similaresSimilares em SciELO


Computación y Sistemas

versão impressa ISSN 1405-5546

Comp. y Sist. vol.21 no.3 México Jul./Set. 2017 

Regular Articles

Trajectory Tracking of Complex Dynamical Network for Chaos Synchronization Using Recurrent Neural Network

Jose P. Perez1 

Joel Perez P.1 

Angel Flores H.1 

Martha S. Lopez de la Fuente2 

1 Universidad Autonoma de Nuevo Leon (UANL), Facultad de Ciencias Fisico Matematicas, Monterrey, Mexico

2 Universidad de Monterrey (UDEM), Monterrey, Nuevo Leon, Mexico


In this paper the problem of trajectory tracking is studied. Based on the Lyapunov theory, a control law that achieves the global asymptotic stability of the tracking error between a recurrent neural network and a complex dynamical network is obtained. To illustrate the analytic results we present a tracking simulation of a dynamical network with each node being just one Lorenz´s dynamical system and three identical Chen’s dynamical systems.

Keywords: Trajectory tracking; recurrent neural net-work; complex dynamical network; Lyapunov analysis

Corresponding author is Jose P. Perez.

1 Introduction

Since the most famous random graph model was proposed by Erdös and Rényi [1], the complex network has attracted much attention in many fields of research, such as biology, physics, computer networks, the World Wide Web (WWW) [2], and so on. Network synchronization has obvious advantages, it has great application value in practice. Therefore, Atay et al. [3] studied synchronization of complex network when delays exist among the nodes; Motter et al. [4] studied the influence of coupling strength on the synchronizing ability of a complex network; Timme et al. [5] studied the web synchronization law of pulse-coupled dynamical systems; Checco et al. [6] studied the synchronization of random web. Lü et al. [7] constructed general complex dynamical networks and studied the synchronization; Lu and Chen [8] studied synchronization analysis of linearly coupled networks of discrete time systems; Han and Lu [9] studied the changes of synchronization ability of coupled networks from ring networks to chain networks; He and Yang [10] studied adaptive synchronization in nonlinearly coupled dynamical networks.

The analysis and control of complex behavior in complex networks, which consist of dynamical nodes, has become a point of great interest in recent studies, [11, 12, 13]. The complexity in networks comes from their structure and dynamics but also from their topology, which often axoects their function. Recurrent neural networks have been widely used in the fields of optimization, pattern recognition, signal processing and control systems, among others. They have to be designed in such a way that there is one equilibrium point that is globally asymptotically stable. In biological and artificial neural networks, time delays arise in the processing of information storage and transmission. Also, it is known that these delays can create oscillatory or even unstable trajectories. Trajectory tracking is a very interesting problem in the field of theory of systems control; it allows the implementation of important tasks for automatic control such as: high speed target recognition and tracking, real-time visual inspection, and recognition of context sensitive and moving scenes, among others.

The motivation in this paper lies in the complex network synchronization and chaos control importance. Network synchronization is one of the most practical and valuable issues. A synchronization of network means the situation in which the output of all nodes in the study of the complex network is consistent with any given external input signal under a certain condition. Numerical simulations are used to verify the effectiveness of the proposed techniques. We present the results of the design of a control law that guarantees the tracking of general complex dynamical networks.

2 Mathematical Models

2.1 General Complex Dynamical Network

Consider a network consisting of N linearly and diffusively coupled nodes, with each node being an n-dimensional dynamical system, described by


where x i = (x i1 , x i2 , . . . , x in ) T ∈ ℝ n are the state vectors of node i, f i : ℝ n 7 ↦ ℝ n represents the self-dynamics of node i, constants c ij > 0 are the coupling strengths between node i and node j, with i, j = 1, 2, . . . , N. Γ = (τij) ∈ ℝ n×n is a constant internal matrix that describes the way of linking the components in each pair of connected node vectors (x j − x i ): that is to say for some pairs (i, j) with 1 ≤ i, j ≤ n and τij ≠ 0 the two coupled nodes are linked through their ith and jth sub-state variables, respectively, while the coupling matrix A = (a ij ) ∈ ℝ N×N denotes the coupling configuration of the entire network: that is to say if there is a connection between node i and node j(ij), then a ij = a ji = 1; otherwise a ij = a ji = 0.

2.2 Recurrent Neural Network

Consider a recurrent neural network in the following form:


where xni = (xni, xni2,,xnin) T ∈ ℝ n is the state vector of neural network i, Uni ∈ ℝ n is the input of neural network i, Ani = -λni I n×n , i = 1, 2, . . . , N, is the state feedback matrix, with λni being a positive constant, Wni ∈ ℝ n×n is the connection weight matrix with i = 1, 2, . . . , N, and σ(·) ∈ ℝ n is a Lipschitz sigmoid vector function [14], [15], such that σ(xni ) = 0 only at xni = 0, with Lipschitz constant Lσi , i = 1, 2, . . . , N and neuron activation functions σ i (·) = tanh(·), i = 1, 2, . . . , n.

3 Trajectory Tracking

The objetive is to develop a control law such that the ith neural network (2) tracks to the trajectory of the ith dynamical system (1). We define the tracking error as e i = x ni − x i , i = 1, 2, . . . , N whose time derivative is


Substituting (1) and (2) in (3), we obtain


Adding and substracting Wni σ(x i ), α i (t), i = 1, 2, . . . , N, to (4), where α i to be determined below, and taking into account that xni = e i + x i , i = 1, 2, . . . , N, then


In order to guarantee that the ith neural network (2) tracks the ith reference trajectory (1), the following assumption has to be satisfied:

Assumption 1. There exist functions ρ i (t) and α i (t), i = 1, 2, . . . , N, such that


Let define


Considering (6) and (7), the equation (5) is reduced to


Rewriting the summations as


also taking into account that cninj = c ij and aninj = a ij , then, using the equations above (8) becomes


It is clear that e i = 0, i = 1, 2, . . . , N is an equilibrium point of (10), when , i = 1, 2, . . . , N. In this way, the tracking problem can be restated as a global asymptotic stabilization problem for system (10).

4 Tracking Error Stabilization and Control Design

In order to establish the convergence of (10) to ei = 0, i = 1, 2, . . . , N, which ensures the desired tracking, first, we propose the following candidate Lyapunov function


The time derivative of (11), along the trajectories of (10), is


Reformulating (12), we get


Next, let consider the following inequality, proved in [16, 17]:


which holds for all matrices X, Y ∈ ℝ n×k and Λ ∈ ℝ n×n with Λ = Λ T > 0. Applying (14) with Λ = I n×n to the term e i T W in φ σ (e i , x i ), i = 1, 2, . . . , N we get


Then we have that, and taking into account that φ σ is Lypchitz:


with Lipschitz constant . Applying (16) to we obtain


By simplifying (15), we obtain


Then we have that:


We define e, i = 1, 2, . . . , N, then (19) becomes in


Now, we propose to use the following control law:


then for all e ≠ 0.

This means that the proposed control law (21) can globally and asymptotically stabilize the ith error system (10), thereby ensuring the tracking of (1) by (2). Finally, the control action driving the recurrent neural networks is given by


5 Simulations

In order to illustrate the applicability of the discussed results, we consider a dynamical network with just one Lorentz’s node and three identical Chen’s nodes. The single Lorentz´s system is described by


and the Chen’s oscillator is described by Eq. 24:


If the system parameters are selected as p 1 = 35, p 2 = 3, p 3 = 28, then the Lorentz’s system and Chen’s system are shown in Fig. 1 and Fig.2 respectively. In this set of system parameters, one unstable equilibrium point of the oscillator (24) is x = (7 : 9373; 7 : 9373; 21)T .

Fig. 1 Sub-State of Lorentz’s attractor with initial condition X 1(0) = (10, 0, 10) T  

Fig. 2 Sub-States of Chen’s attractor with initial condition X2;3;4(0) = (10, 0, 37) T  

Suppose that each pair of two connected Lorentz and Chen’s oscillators are linked together through their identical sub-state variables, i.e., Γ = diag(1, 1, 1), and the coupling strengths are c 12 = c 21 = π, c 13 = c 31 = π, c 23 = c 32 = π, c 14 = c 41 = 2π, c 24 = c 42 = 2π, c 34 = c 43 = 2π. Fig. 3 visualizes this entire dynamical network:

Fig. 3 Structure of the network with each node being a Lorentz and Chen’s system 

The neural network was selected as


Theorem 1 For the unknown nonlinear sys-tem modeled by (1), the on-line learning law tr {W T W} = −e T W σ(x) and the control law (22) ensure the tracking of to the nonlinear reference model (2).

Remark 2 From (20) we have, ∀ e ≠ 0, ∀W, and therefore V is decreasing and bounded from below by V (0). Since , then we conclude that e, WL 1; this means that the weights remain bounded.

The experiment is performed as follows. Both systems, the recurrent neural network (2) and the dynamical networks (24) and (25), evolve independently; at that time, the proposed control law (23) is incepted. Simulation results are presented in Fig. 4 - Fig. 6 for sub-sates of node 1. As can be seen, tracking is successfully achieved and error is asymptotically stable, as it is shown in Fig. 7 -Fig. 9 for sub-states of node 4.

Fig. 4 Time evolution for sub-states 1 with initial state X n1 (0) = (10, 0, 10) T  

Fig. 5 Time evolution for sub-states 1 with initial state X n1 (0) = (10, 0, 10) T  

Fig. 6 Time evolution for sub-states 2 with initial state X n1 (0) = (10, 0, 10) T  

Fig. 7 Time evolution for sub-states 4 with initial state X n4 (0) = (20, 20, 10) T  

Fig. 8 Time evolution for sub-states 4 with initial state X n4 (0) = (20, 20, 10) T  

Fig. 9 Time evolution for sub-states 4 with initial state X n4 (0) = (20, 20, 10) T  

6 Conclusions

We have presented the controller design for trajectory tracking determined by a general complex dynamical network. This framework is based on dynamic recurrent neural networks and the methodology is based on Lyapunov theory. The proposed control is applied to a dynamical network with each node being a Lorenz and Chen’s dynamical system, respectively, being able to also stabilize in asymptotic form the tracking error between two systems. The results of the simulation shows clearly the desired tracking. In future work, we will consider the stochastic case for the complex dynamical network.


The authors thank the support of CONACYT and the Dynamical Systems group of the Facultad de Ciencias Fisico-Matematicas, Universidad Autonoma de Nuevo Leon, Mexico.


1. Erdös, P., & Rényi, A. (1961). On the evolution of random graphs. Publication of the Mathematical Institute of the Hungarian Academy of Sciences, Vol. 5, pp. 17. [ Links ]

2. Newman, M. E. J., Strogatz, S. H., & Watts, D. J. (2001). Random graphs with arbitrary degree distributions and their applications. Physical ReviewE, Vol. 64, Article ID26118, 17 pp. [ Links ]

3. Atay, F. M., Jost, J., & Wende, A. (2004). Delays, connection topology, and synchronization of coupled chaotic maps. Physical Review Letters, Vol. 92, No. 14, pp. 144101-1. [ Links ]

4. Motter, A. E., Zhou, C., & Kurths, J. (2005). Network synchronization, diffusion, and the paradox of heterogeneity. Physical Review E, Vol. 71, Article ID016116, 9 pp. [ Links ]

5. Timme, M., Wolf, F., & Geisel, T. (2004). Topological speed limits to network synchronization. Physical Review Letters, Vol. 92, Article ID 074101, 4 pp. [ Links ]

6. Checco, P., Biey, M., & Kocarev, L. (2008). Synchronization in random networks with given expected degree sequences. Chaos, Solitons & Fractals, Vol. 35, No. 3, pp. 562-577. [ Links ]

7. Lü, J., Yu, X., & Chen, G. (2004). Chaos synchronization of general complex dynamical networks. Physica A, Vol. 334, No. 1-2, pp. 281-302. [ Links ]

8. Lu, W., & Chen, T. (2004.) Synchronization analysis of linearly coupled networks of discrete time systems. Physica D, Vol. 198, No. 1-2, pp. 148-168. [ Links ]

9. Han, X., & Lu, J. (2007). The changes on synchronizing ability of coupled networks from ring networks to chain networks. Science in China F, Vol. 50, No. 4, pp. 615-624. [ Links ]

10. He, G., & Yang, J. (2008). Adaptive synchronization in nonlinearly coupled dynamical networks. Chaos, Solitons and Fractals, Vol. 38, No. 5, pp. 1254-1259. [ Links ]

11. Wang, X (2002). Complex networks: Topology, dynamics and synchronization. Int. J. Bifurcation Chaos, Vol. 12, No. 5, pp. 885-916. [ Links ]

12. Wu, C (2002). Synchronization in Coupled Chaotic Circuits and Systems. World Scientific, Singapore. [ Links ]

13. Newman, M (2003). The structure and function of complex networks. SIAM Review, Vol. 45, pp. 167- 256. [ Links ]

14. Khalil, H (1996). Nonlinear System Analysis. 2nd. Ed. Prentice Hall, Upper Saddle River, NJ, USA. [ Links ]

15. Krstic, M., & Deng, H. (1998). Stabilization of Nonlinear Uncertain Systems. Springer Verlag, USA. [ Links ]

16. Sanchez, E., & Perez, P. (1999). Input-to-state stability analysis for dynamic neural networks. IEEE Trans. on Circuits Syst. I, Vol. 46, pp. 1395-1398. [ Links ]

17. Sanchez, E., Perez, J. P., & Perez, J. (2006). Trajectory Tracking for Delayed Recurrent Neural Networks. American Control Conference, pp. 14-16 [ Links ]

18. Li, X., Wang, X., & Chen, G. (2004). Pinning a complex dynamical network to its equilibrium. IEEE Transactions on Circuits and Systems I: Regular Papers, Vol. 51, No. 10, pp. 2074-2087. [ Links ]

Received: January 13, 2015; Accepted: April 17, 2017

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License