Synchronization of complex networks of identical and nonidentical chaotic systems via model-matching control

Abstract

In this work, a synchronization scheme for networks of complex systems is presented. The proposed synchronization scheme uses a control law obtained with some definitions from graph theory and solving the Model-Matching Problem for complex networks. In particular, Rössler, Chen, Lorenz and Lü chaotic systems are used as complex chaotic systems into complex networks. Particular cases with regular and irregular networks of six identical chaotic systems are implemented, with some well-known topologies as star and ring small-world, and tree topologies. Highlighting, the obtained control law is applied to synchronize an irregular network of six different chaotic systems in a tree topology. The usefulness and advantages of the proposed synchronization scheme are highlighted performing numerical simulations of the chaotic complex networks.

About the Authors:

D. López-Mancilla

Roles Conceptualization, Investigation, Supervision, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: Centro Universitario de los Lagos, Universidad de Guadalajara (UDG), Jalisco, México

ORCID logo http://orcid.org/0000-0002-0971-3176

G. López-Cahuich

Roles Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft

Affiliation: Centro Universitario de los Lagos, Universidad de Guadalajara (UDG), Jalisco, México

C. Posadas-Castillo

Roles Investigation, Methodology, Supervision, Validation, Writing – original draft

Affiliation: Universidad Autónoma de Nuevo León (UANL), Nuevo León, México

C. E. Castañeda

Roles Conceptualization, Formal analysis, Investigation

Affiliation: Centro Universitario de los Lagos, Universidad de Guadalajara (UDG), Jalisco, México

ORCID logo http://orcid.org/0000-0002-0781-0490

J. H. García-López

Roles Formal analysis, Investigation, Supervision, Validation, Visualization

Affiliation: Centro Universitario de los Lagos, Universidad de Guadalajara (UDG), Jalisco, México

J. L. Vázquez-Gutiérrez

Roles Formal analysis, Software

Affiliation: Centro Universitario de los Lagos, Universidad de Guadalajara (UDG), Jalisco, México

E. Tlelo-Cuautle

Roles Conceptualization, Resources, Supervision

Affiliation: Instituto Nacional de Astrofísica Óptica y Electrónica (INAOE), Puebla, México

ORCID logo http://orcid.org/0000-0001-7187-4686

Introduction

The field of synchronization of networks of complex systems has received a lot of attention in the last three decades, due to the potential applications in engineering, the proliferation of computer networks, communications networks as the internet, wireless communications as cellular telephony and many others [1, 2]. Different methods and topologies have been proposed for network synchronization. For example, the authors in [3–6] introduced schemes for the synchronization in small-world dynamical networks. Serrano and co-workers made a contribution on chaotic synchronization in star coupled neural networks [7]. Model-Matching Control (MMC) has been used as a solution for some open problems. It has been recently used for output synchronization of heterogeneous and nonlinear multi-agent systems [8, 9]. In other work, the authors consider the model matching problem for a class of multiple-input systems whose inputs involve different delays [10]. In chaos synchronization area, MMC has been used in the last years, but only to synchronize pairs of identical and nonidentical chaotic oscillators [11–13]. Chaotic synchronization of regular and irregular complex networks with fractional order oscillators has been presented in [14]. Nevertheless, most of the authors show results in synchronization for only structurally identical complex systems, whereas for synchronization of nonidentical complex systems, the results are only a few. Synchronization, output-synchronization, partial-synchronization or phase-synchronization of nonidentical complex systems would have implication to understand different areas. For example in neural networks, the interaction between pre-synaptic and post-synaptic neurons, where the pre-synaptic neurons output activates the nucleus cells in post-synaptic neurons to achieve some kind of synchronization between two different neurons [15]. In communications systems, sometimes the network topology has different complex systems and it is required to synchronize at least the nodes outputs. This is necessary to encryption/decryption and transmission of information [16].

The aim of this work is also devoted to demonstrate the effectiveness of the MMC for synchronizing networks of complex/chaotic systems in continuous-time. This objective is achieved by using the model-matching approach from nonlinear control theory [17, 18] and extending the results in [11, 12]. The main contribution is the generalization of a MMC for network synchronization and highlight that the proposed scheme has the following advantages:

* It is a systematic procedure to be used as guideline or a kind of recipe that anyone can reproduce.

* It can be useful to synchronize networks of identical and nonidentical complex/chaotic systems.

* It can be used in unidirectional and bidirectional coupling between master and slave systems.

* It can be used to synchronizing chaotic and hyperchaotic systems [19].

* Any node can be chosen as a master node.

* There may be more than one master node.

The proposed MMC is tested to synchronize networks of different chaotic systems, like: Rössler, Chen, Lorenz and Lü chaotic oscillators. Two classic regular network topologies are taken as cases of study, namely a star topology with an isolated node, and a ring topology. In the same way, two irregular topologies are considered, such as:: a small-world topology and a tree topology. All topologies use six chaotic systems as nodes. The election of using six nodes can be extended to larger networks. Nevertheless, it is well known that if the number of the nodes is big enough, then complete synchronization could be guaranteed only for some topologies [20]. Readers interested in pattern formation in complex networks and some new results see [21–24].

This work is organized as follows: Section 2 states the problem formulation for synchronizing a network of chaotic systems. The model-matching problem from nonlinear control theory applied to network synchronization is presented in Section 3. In Section 4, this approach is applied to synchronize networks of identical and nonidentical chaotic systems based on Rössler, Lorenz, Lü and Chen systems, with different topologies. Finally, Section 5 summarizes the concluding remarks.

Problem statement

Consider a network of complex dynamical systems described by state equations of the form:(1)where the states the inputs and the outputs , being f (xi) and g (xi) smooth and analytical functions and with i = 1, 2, …, N, where N is the number of nodes in the network. In addition, consider another nonlinear system described by:(2)where the state the input and the output , being fM (xM) and gM (xM) also smooth and analytical functions. Assume that are equilibrium points of systems described in (1), i.e., . Similarly, is an equilibrium point of system (2). Assume that the dynamical systems of (1) and (2) under certain conditions have chaotic behavior. Then, the chaotic system (1) synchronizes with the chaotic system (2), if:(3)no matter which initial conditions xi (0) and xM (0) have, and for each node i = 1, 2, …, N and suitable input signals ui (t) and uM (t).

Note that, it is mainly considering output synchronization problem between a network of chaotic systems (1) and (2). Moreover, no matter if the chaotic systems (1) and (2) are identical or nonidentical. In the next section, it will be described how to satisfy the output synchronization condition (3) from the perspective of the model-matching problem and graph theory.

Model-matching for networks

Let us consider the dynamical systems (1) like plants Pi, and (2) like a model M, respectively. The objective is to generalize a feedback control law u (t) for several plants Pi which, irrespectively of the initial states of Pi and M, makes the outputs yi (t) converge asymptotically to the output yM (t) produced by M under an arbitrary input uM (t). This problem is the called asymptotic model-matching problem from nonlinear control theory, and it was solved for synchronizing only a pair of chaotic systems in [11, 12]. Previously, different approaches to solve the model-matching problem have been proposed in the literature, see e.g. [17, 18]. In this work, the solution proposed in [11, 12] is extended to synchronize several chaotic plants with one chaotic model. In the design of the MMC we considered the possibility of using more than one chaotic model into a network of complex chaotic systems. Then, for network synchronization purpose, N auxiliary systems are defined by the following equation:(4)with state , inputs u (t) and uM (t), andNote that the output yEi = hE (xEi) of the auxiliary system (4) is the difference between the output of the plant Pi [yi = h (xi) of Eq (1)] and the output of system M [yM = hM (xM) of Eq (2)]. The control objective of the model-matching problem is contained in the following definition:

Definition 1 (Model-matching problem): Given the plants Pi and the model M around their respective equilibrium points and and points . The model-matching problem consists of finding feedback control laws for the auxiliary system Ei (4) such that, the output yEi (t) → 0 as t → ∞.

The equilibrium points in Definition 1, are stated in the following Definition:

Definition 2 (Relative degree adapted from [18]): The single-input single-output nonlinear system (1), is said to have relative degree ri at points if:

1. for all xi in a neighborhood of and for all k < ri − 1,

2. .

In Definition 2, and . A similar definition can be given for the relative degree of model (2), rM near . It is important to mention that the model matching problem is locally solvable if, and only if [18]:(5)

Now, let the auxiliary system Ei Eq (4) be in a different coordinate frame. From definition of relative degrees ri and rM; and are sets of independent functions from Pi and M, and can be chosen as new coordinates and with q = 1, …, n, around and , respectively. Consider now the auxiliary system Ei and the new coordinates [18]:where ζ (xEi) = (ζ1,i (xEi), …, ζn,i (xEi))⊤, and q = 1, …, n.

Thus, the closed-loop auxiliary system Ei, using the following feedback control law(6)takes the form:(7)

Two subsystems can be identified in the closed-loop system (7), namely:

1. The subsystem described by:which represents the dynamics of M, and

2. The subsystem described by:withwhich represents the dynamics of yEi (t). The matrix can be a single constant matrix choosing i = 1 or it can define different matrix, specially when they are used for synchronization of non-identical chaotic systems into a complex network.

Model M is stable by assumption, since it is proposed, and the control law vi (t) can be selected such that eigenvalues of matrix A* have a negative real part. Then the closed-loop system will be exponentially stable, and output synchronization condition (3) holds.

Complex networks and model-matching

From a mathematical point of view, a complex network is defined as an interconnected set of a nodes (two or more) and can be represented by a graph, where two nodes (vertices or points) joined by a connection (edges or lines) are called adjacent nodes or neighbors. Some of the most notable features for complex systems are:

* They consist of many interacting parts (nodes).

* Each part has its own internal structure and is responsible for a specific task.

Topology or coupling mesh is the layout or how the nodes of a network are connected, while the configuration is the type of connection that determines the flow of information between nodes. Now, consider the traditional control law for network synchronization [3, 14]where N denotes the size of the network or the number of nodes, c > 0 represents the coupling strength, and is a constant matrix linking the state variables. The matrix is the coupling matrix. If there is a connection between nodes i and j then, the element ai,j = 1; otherwise ai,j = 0, i ≠ j. For i = j the diagonal elements of A are defined as:

Then, given the previous result and the above definitions, it is possible to formulate a general control law for network synchronization using Model-Matching control, and choosing ri = rM = r, as:

This is:andwhere xi and xj represents the states of the chaotic systems nodes, plants and models, respectively. Notice that with this control law, any node can be chosen as the model system M.

Network synchronization through model-matching control for identical and nonidentidal systems

In this section, we use the previous material in order to illustrate how network synchronization of N chaotic systems can be achieved. Two cases of study are considered, using identical and nonidentical chaotic systems like nodes in regular and irregular topologies. The number of chaotic nodes are limited to N = 6 nodes. Fig 1 shows the block diagram of model-matching control and Fig 2 presents the auxiliary system block diagram. The model matching control works as a similar way like previous works from one of the authors, but adapted to network synchronization [11, 12]. Thus, for identical chaotic systems, complete synchronization is guaranteed for a relative degree r = n. For nonidentical chaotic systems only output synchronization is guaranteed.

[Figure omitted. See PDF.]

Fig 1. Block diagram of model-matching control.

https://doi.org/10.1371/journal.pone.0216349.g001

[Figure omitted. See PDF.]

Fig 2. Auxiliary system block diagram.

https://doi.org/10.1371/journal.pone.0216349.g002

Table 1 presents the equations of the chaotic systems and the parameter values for Rössler, Chen, Lorenz and Lü chaotic systems, used to illustrate synchronization using Model-Matching Control. Although the proposed approach can be applied to any chaotic system that holds (5) and for all plant Pi with a strong relative degree. In Table 1, Rössler systems have complete strong relative degree r = 3, the rest of the systems have a strong relative degree r = 2. For synchronization of identical systems with relative degree r = 2, stability demonstration is not trivial. But, the first two equations of the auxiliary system become linear and asymptotically stable, and the last one presents a zero dynamic isolated from the output, depending only from the third state, obtaining the form , for some positive constant β. Then, the solution ζ3 = exp (−βt) is exponentially stable.

[Figure omitted. See PDF.]

Table 1. Chaotic systems and parameter values.

https://doi.org/10.1371/journal.pone.0216349.t001

Regular networks

Some classic regular networks are used to illustrate the synchronization of networks of chaotic systems. In particular, Star topology and Ring topology are presented.

Star topology with an isolated node.

Consider a network in a Star topology with six Chen hyperchaotic systems. Fig 3 shows the Star topology for six hyperchaotic systems and the associated coupling matrix. The temporal signals of six Chen hyperchaotic systems synchronizing are illustrated in Fig 4, whereas in Fig 5 the synchronization graphics for the outputs of six Chen hyperchaotic systems are presented. See [19] for discrete-time context.

[Figure omitted. See PDF.]

Fig 3. a) Star topology, b) coupling matrix, c) random initial conditions.

https://doi.org/10.1371/journal.pone.0216349.g003

[Figure omitted. See PDF.]

Fig 4.

a) Temporal signals for the synchronization of six Chen chaotic systems in a star topology with an isolated node. b) Zoom in window showing the convergence in time for all xi1 states.

https://doi.org/10.1371/journal.pone.0216349.g004

[Figure omitted. See PDF.]

Fig 5.

a) Synchronization graphics for six Chen chaotic systems in a star topology with an isolated node. b) Chen chaotic attractor.

https://doi.org/10.1371/journal.pone.0216349.g005

Ring topology.

The second case is a Ring topology with six Lorenz systems shown in Fig 6 and its associated coupling matrix. Fig 7 illustrates the temporal signals on the synchronization of six Lorenz systems in a Ring topology and Fig 8 shows the synchronization graphics for the outputs of six Lorenz systems.

[Figure omitted. See PDF.]

Fig 6. a) Ring topology b) coupling matrix, c) random initial conditions.

https://doi.org/10.1371/journal.pone.0216349.g006

[Figure omitted. See PDF.]

Fig 7.

a) Temporal signal on the synchronization for six Lorenz chaotic systems in a ring topology. b) Zoom in window showing the convergence in time for all xi1 states.

https://doi.org/10.1371/journal.pone.0216349.g007

[Figure omitted. See PDF.]

Fig 8.

a) Synchronization graphics for six Lorenz chaotic systems in a ring topology with an isolated node. b) Lorenz chaotic attractor.

https://doi.org/10.1371/journal.pone.0216349.g008

Irregular networks

One of the most interesting cases is making irregular networks because generate some behaviour that is not expected. It could be more complex if the elements of the irregular network are structurally different. This case is presented within the irregular Small-World topology and Tree topology. The first case is realized with identical chaotic systems and the second one with different chaotic systems.

Small-world topology.

Small-world topology is a very used topology for the effectiveness to achieve synchronization. Six Rössler chaotic systems are synchronized when they are connected as shown in Fig 9 with the corresponding coupling matrix. Fig 10 shows the temporal output signals and the synchronization in a small time. Fig 11 presents the synchronization graphics for six Rössler chaotic systems in a Small-world topology.

[Figure omitted. See PDF.]

Fig 9. a) Small-World topology b) coupling matrix, c) random initial conditions.

https://doi.org/10.1371/journal.pone.0216349.g009

[Figure omitted. See PDF.]

Fig 10.

a) Synchronization graphics for six Rössler chaotic systems in a Small-World topology. b) Zoom in window showing the convergence in time for all xi2 states.

https://doi.org/10.1371/journal.pone.0216349.g010

[Figure omitted. See PDF.]

Fig 11.

a) Synchronization graphics for six Rössler chaotic systems in a Small-World topology with an isolated node. b) Rössler chaotic attractor.

https://doi.org/10.1371/journal.pone.0216349.g011

Tree topology.

In this case, an irregular Tree topology is realized using an array presented in Fig 12 with its associated coupling matrix. Notice that it uses three different chaotic systems, two Lorenz systems, two Lü systems and two Chen systems. The main master is labeled as Lorenz 1 at the node 1. The node 2 contains the Lorenz 2 that is a slave system driven by Lorenz 1 and is another Master system. For non-identical chaotic systems, output synchronization if guaranteed in [11]. Then, the expected result is output synchronization. Fig 13 shows the temporal series for output synchronization for six nodes using three different chaotic systems in a Tree topology with two master systems. The output synchronization result is illustrated again in Fig 14 with the synchronization graphics. The complete and more interesting result is shown in Fig 15. All chaotic systems show output synchronization, whereas the other states show phase synchronization with respect to the main master node 1. Nevertheless, the states between both Chen systems and both Lü systems show complete synchronization in other new state in phase with the main master.

[Figure omitted. See PDF.]

Fig 12. a) Tree topology with different chaotic systems, b) coupling matrix, c) random initial conditions.

https://doi.org/10.1371/journal.pone.0216349.g012

[Figure omitted. See PDF.]

Fig 13.

a) Output Synchronization for six nodes using three different chaotic systems in a Tree topology with two master systems. b) Zoom in window showing the convergence in time for all xi1 states.

https://doi.org/10.1371/journal.pone.0216349.g013

[Figure omitted. See PDF.]

Fig 14. Output Synchronization graphics for different chaotic systems in a Tree topology.

https://doi.org/10.1371/journal.pone.0216349.g014

[Figure omitted. See PDF.]

Fig 15. Output Synchronization, phase synchronization and complete synchronization for different chaotic systems in a Tree topology.

https://doi.org/10.1371/journal.pone.0216349.g015

Discussion

The MMC for network synchronization presented in this work is competitive with the most popular methodology used in [1–7, 14]. In that case, an error feedback between the signals to be synchronized is used. Synchronization time depends on the coupling gains k. On the other hand, in this work, synchronization time depends on the poles location. An adequate election on the poles location can achieve robust synchronization and, at the same time, an unforced control law. Nevertheless, some limitations were observed in the process: If Initial conditions are much bigger respect to the maxima amplitude of the temporal chaotic signal, the poles of the linearized auxiliary system need to be more negative from the complex axis of complex plane to achieve synchronization. This makes the MMC be increased, generating a forced control law. Other limitation is considered in Eq (5), where model matching problem is locally solvable if, and only if ri ≤ rM. Nevertheless, it is possible to obtain output synchronization for systems with different order, if the systems hold condition (5). The main advantages for network synchronization using MMC are the following: It is a systematic procedure, it is robust, can be applied for nonidentical chaotic systems, and the initial conditions rank is bigger respect to other methodologies. They could be considered as disadvantages: Mathematical analysis for computing the control law is complex. When synchronization is more robust, a more forced control law is obtained. This is not desirable for physical implementation purpose, but it could be solved using hybrid systems (computer + physical implementation) [25]. This could be truth because the MMC does not increase the computation complexity, i.e., the computer time running the program depends only on the computer hardware and software.

Conclusions

In this work, a synchronization scheme for complex networks of identical and nonidentical chaotic systems was presented. In particular, model-matching problem from nonlinear control theory was used. The results show complete synchronization for networks of identical systems and output synchronization for an irregular network of a class of non-identical systems. Although for non-identical systems, complete synchronization was achieved between Chen’s and Lü’s systems. The advantages over other cited approaches to synchronize networks of nonidentical chaotic systems are the following: this approach is a systematic procedure, it could use unidirectionally and bidirecctionally coupled systems, gains for the controller could be small, and synchronization network is obtained after a short transient behavior. Moreover, the proposed control law allows to chose any chaotic system as a model and more than one model. Numerical simulations were done. In future works, network layers and some applications for neural network communication systems will be studied. At the same time, some studies on the rank of the initial conditions and poles location could be analyzed.

Acknowledgments

The authors acknowledge the financial support from the University of Guadalajara under the project RC/075/2018, Agreement RG/006/2018, UDG, México.

Citation: López-Mancilla D, López-Cahuich G, Posadas-Castillo C, Castañeda CE, García-López JH, Vázquez-Gutiérrez JL, et al. (2019) Synchronization of complex networks of identical and nonidentical chaotic systems via model-matching control. PLoS ONE 14(5): e0216349. https://doi.org/10.1371/journal.pone.0216349

References

1. Estrada E.’Introduction to complex networks: structure and dynamics’, in Banasiak J., Mokhtar-Kharroubi M. (Ed.) “Evolutionary equations with applications in natural sciences,” (Springer, 2015, 1st edn.), pp. 93–131.

2. Arenas A., Díaz-Guilera A., Kurths J., Moreno Y., Zhou C., “Synchronization in complex networks,” Physics Reports, 2008, vol. 469, no 3, pp. 93–153.

3. Soriano-Sánchez A.G., Posadas-Castillo C., Platas-Garza M.A. “Synchronization of Generalized Chua’s Chaotic Oscillators in Small-world Topology,” J. Engineering Sci. Technology Review, 2015, 8, (2), pp. 185–191

4. Soriano-Sánchez A.G., Posadas-Castillo C. “Smart pattern to generate small–world networks,” Chaos, Solitons & Fractals, 2018, 8, (114), pp. 415–422

5. Wang X.F. “Complex networks “Topology, dynamic and synchronization,” International Journal of Bifurcation and Chaos, 2002, vol. 12, no. 5, pp. 885–916.

6. Wang X.F., Chen G. “Synchronization in small-world dynamical networks,” International Journal of Bifurcation and Chaos, 2002, vol. 12, no. 1, pp. 187–192.

7. Serrano-Guerrero H., Cruz-Hernández C., López- Gutiérrez R.M., Posadas-Castillo C., Inzunza-González E. “Chaotic synchronization in star coupled networks of three-dimensional cellular neural networks and its application in communications,” International Journal of Nonlinear Sciences and Numerical Simulation, 2010, vol. 11, no. 8, pp. 571–580.

8. Zhu L., Chen Z., Middleton R. “Robust asymptotic model matching and its application to output synchronization of heterogeneous multi-agent systems,” Syst. Control Lett, 2016, 94, pp. 37–45

9. Seyboth, G.S., Allgower, F. “Synchronized model matching: a novel approach to cooperative control of nonlinear multi-agent systems,” Proc. 19th IFAC World Congress, Cape Town, South Africa, Aug 2014, pp. 1985–1990.

10. Su H., Jia Y. “Model matching control of multiple-input-delay systems,” IET Control Theory Appl., 2010, 4, (12), pp. 2719–2732

11. López-Mancilla D., Cruz-Hernández C. “Output Synchronization of Chaotic Systems: Model-Matching Approach with Application to Secure Communication,” Nonlinear Dyn. Syst. Theory, 2005, 5, (2), pp. 141–156

12. López-Mancilla D., Cruz-Hernández C. “Output synchronization of chaotic systems under nonvanishing perturbations,” Chaos, Solitons & Fractals, 2008, 37, pp. 1172–1186

13. Pisarchik A.N., Jaimes Reategui R., García-López J. H. “Synchronization of multistable systems,” International Journal of Bifurcation and Chaos, 2008, vol. 18, No. 6, pp. 1801–1819

14. Angulo-Guzman S., Posadas-Castillo C., Platas-Garza M.A., Díaz-Romero D.A., López-Mancilla D. “Chaotic Synchronization of Regular and Irregular Complex Networks with Fractional Order Oscillators,” Int. J. Control, Autom. and Syst., 2016, 14, (4), pp. 1114–1123.

15. Pisarchik AN, Sevilla-Escoboza R, Jaimes-Reátegui R, Huerta-Cuellar G, García-López JH, Kazantsev VB. Experimental Implementation of a Biometric Laser Synaptic Sensor, Sensors. 2013; 13: 17322–31. pmid:24351638

16. Cruz-Hernández C, Martynyuk AA. Advances in chaotic dynamics and applications. Series on stability, oscillations and optimization of systems, vol. 4. Cambridge: Cambridge Scientific Publishers Ltd; 2009.

17. Di Benedetto M.D., Grizzle J.W. “Asymptotic model matching for nonlinear systems,” IEEE Trans. Automatic Control, 1994, 39, (8), pp. 1539–1549.

18. Isidori A. “Nonlinear Control Systems,” (Springer, 1995)

19. Aguilar–Bustos A.Y., Cruz–Hernández C., López–Gutiérrez R.M. Posadas–Castillo C. “Synchronization of Different Hyperchaotic Maps for Encryption” Nonlinear Dynamics and Systems Theory, 8(3), 2008, pp. 221–236.

20. Posadas-Castillo, C. “Sincronización de sistemas complejos (Synchronization of Complex Networks),” PhD Thesis, UABC, México, 2008, pp. 19–40.

21. Qian Yu, Liu Fei, Yang Keli, Zhang Ge, Yao Chenggui, Ma Jun “Spatiotemporal dynamics in excitable homogeneous random networks composed of periodically self-sustained oscillation” Scientific Reports, 7, 11885, 2017, pp. 1–13.

22. Zheng Zhi-gang, and Qian Yu, “Dominant phase-advanced driving analysis of self-sustained oscillations in biological networks” Chin. Phys. B Vol. 27, No. 1 (2018) 018901.

23. Wanga Xin, Liu Xinzhi, She Kun, Zhong Shouming, “Pinning impulsive synchronization of complex dynamical networks with various time-varying delay sizes”, Nonlinear Analysis: Hybrid Systems, 26, 2017, pp. 307–318.

24. Yang Huilan, Wang Xin, Zhong Shouming, Shu Lan, “Synchronization of nonlinear complex dynamical systems via delayed impulsive distributed control”, Applied Mathematics and Computation, 320, 2018, pp. 75–85.

25. Chiu R., López-Mancilla D., Castañeda C.E., Orozco-López O, Villafaña-Rauda E., Sevilla-Escoboza R. “Design and implementation of a jerk circuit using a hybrid analog–digital system”, Chaos, Solitons and Fractals 119, 2019, pp. 255–262.