Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2009 (2009), Article ID 174546, 14 pages
Research Article

Diffusive Synchronization of Hyperchaotic Lorenz Systems

Department of Electrical Engineering, School of Engineering at São Carlos, University of São Paulo, 13566-590 São Carlos, SP, Brazil

Received 6 January 2009; Accepted 2 March 2009

Academic Editor: José Roberto Castilho Piqueira

Copyright © 2009 Ruy Barboza. 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.


The synchronizing properties of two diffusively coupled hyperchaotic Lorenz 4D systems are investigated by calculating the transverse Lyapunov exponents and by observing the phase space trajectories near the synchronization hyperplane. The effect of parameter mismatch is also observed. A simple electrical circuit described by the Lorenz 4D equations is proposed. Some results from laboratory experiments with two coupled circuits are presented.

1. Introduction

Coupled oscillators are currently studied in physics, chemistry, biology, neural networks, and other fields. A large number of coupled oscillators form a complex system for which the investigation of coherence and synchronization is important and in the last decade the discovery that chaotic systems can synchronize added interest to this topic [14]. An approach for contributing to the investigation of large networks is studying the properties of a small number of coupled oscillators—this is the approach used in the present work. Chaos synchronization is also of interest in secure communication systems. For such application, hyperchaos has drawn attention for providing more complex waveforms than simply chaotic systems, thus improving the masking process. This is because hyperchaos is characterized by at least two positive Lyapunov exponents, while simple chaos shows a single one. Related to this subject is the question on how many variables are necessary to be coupled in order to obtain synchronization. Although chaotic systems can synchronize by a single variable coupling, it was for some time believed that in the case of hyperchaos the minimum number of coupling variables had to be equal to the number of positive Lyapunov exponents [5]. It was later demonstrated in [6] that it is not true, and some hyperchaotic systems can achieve synchronization by a single variable. Other systems, however, for example the Rössler equations for hyperchaos, are unable to synchronize by only one of its variables [6, 7]. Another problem refers to synchronization under parameter mismatch, since usually perfect synchronization is achieved only if the coupling systems are identical. In this work we are concerned with the question of whether a hyperchaotic Lorenz system can synchronize and with which variables. The effect of mismatch is also observed. We present some results from numerical and laboratory experiments on an eight-dimensional dynamical system obtained by diffusively coupling two hyperchaotic Lorenz systems.

2. Equations and Numerical Simulations

2.1. Diffusive Coupling

Consider the four-dimensional system , where . Unidirectional diffusive coupling involving two such systems is obtained by using the system to drive the response system . In this work we consider a diagonal matrix, where each element is the strength of the coupling related to the corresponding variable. We also consider , and each entry , , , or is either zero or equal to a positive constant . For example, the case , , , and means simultaneous and coupling. The two coupled systems form a compound eight-dimensional system. If the systems synchronize, the motion must remain on the hyperplane . For small we have , so we can write the variational equation where and . In the case of bidirectional diffusive coupling we have and , therefore (2.1) also applies, now with . We use (2.1) as the locally linear dynamical system associated to a fiducial trajectory [8] of the coupled system to calculate the transverse Lyapunov exponents [2, 3, 7, 9], so called because the perturbation is transverse to the synchronization hyperplane. The coupled system will remain stably synchronized if all transverse Lyapunov exponents (TLEs) are negative.

2.2. Complete Replacement

If the drive system 1 transmits the scalar component and the corresponding variable of the response system 2 is replaced by the transmitted one , this is called complete replacement [4]. As explained in [4], unidirectional diffusive coupling and complete replacements are related, since at very high values of the variable slaves . Therefore, in our numerical and experimental investigations, complete replacement of one or more variables corresponds to .

2.3. Coupled Hyperchaotic Lorenz Systems

The preceding coupling scheme will now be applied to systems described by the following Lorenz equations linearly extended to four dimensions: which shows hyperchaos and was theoretically analyzed in [10]. In the case of unidirectional coupling the eight-dimensional system is given by For this system the matrix is for unidirectional coupling. In the case of bidirectional coupling, the only modification is replacing by along the diagonal. The synchronizing properties of system (2.2) will now be numerically investigated by calculating the TLE of (2.3) as a function of the coupling strength . Parameter mismatch is also examined. In the following, only unidirectional coupling is considered.

2.3.1. Parameters

We first illustrate the synchronizing properties of (2.2) for the above parameter values. The values and are the classical, or most popular, used in studies of the original Lorenz 3D system. For the hyperchaotic Lorenz 4D system the extra parameter is included. In Figure 1, the two largest TLE for single-variable coupling are plotted as a function of , for . The synchronization thresholds are and for and couplings, respectively. The systems will never get synchronized if alone is the coupling variable. On the other hand, coupling provides a small window of stable synchronization.

Figure 1: Transversal Lyapunov exponents (TLEs) for

2.3.2. Parameters

The above values of the parameters are of interest in this work because small values of and are easier to realize with practical component values in the circuit model presented in Section 3. (At this point it is worth remembering the observation by Sparrow [11] that small leads to very complex behavior of the Lorenz equations.) Therefore, the system properties for such small parameter values will be examined with more detail in the following. In Figure 2, the Lyapunov spectrum for these parameter values is plotted as a function of , showing a broad range of hyperchaotic behavior. Also shown is the one-dimensional bifurcation diagram along the same range.

Figure 2: (a) Lyapunov spectrum ; (b) bifurcation diagram obtained from the crossings through the Poincaré surface , corresponding to maxima of . Parameter values:

The two largest TLEs for single-variable coupling, plotted as a function of , are shown in Figures 3(a)3(c). The synchronization thresholds are and for and couplings, respectively. For coupling, the system shows a window of stable synchronization ranging from to The system does not synchronize in the case of coupling. However, the nonsynchronizing variables and , when working together in the plus double-coupling scheme, provide stable synchronization above In the case of all-variable coupling, stable synchronization is achieved above

Figure 3: Transversal Lyapunov exponents (TLEs) for (a) and (b) One-variable coupling; (c) two- and four-variable couplings.

2.3.3. Parameter Mismatch

A qualitative method of investigating the hyperchaos synchronization phenomena is by observing the projections of the eight-dimensional attractor onto the planes , , , and . In these planes the straight lines , , , and correspond to the synchronization hyperplane. In the following we show only the plane, since the components and seem to be the most difficult to synchronize. For all plots we used coupling. Figure 4 refers to identical parameters (i.e., without mismatch), illustrating the inability of the systems to synchronize if is less than the threshold value obtained from Figure 3(a), while perfect synchronization is achieved above the threshold: the same alignment along the diagonal is observed in all the four projection planes.

Figure 4: Synchronization for identical parameters: below and above the threshold

In Figures 5 and 6, we observe the effect of mismatch on synchronization (for coupling). In these examples we applied the same mismatch to all parameters, that is, . In Figure 5, where , we see that for values just above the threshold, some good degree of synchronization is obtained for 1% mismatch; however, for 5% large deviations from the diagonal are observed. Figure 6 shows the effect of mismatch for 1% and 5% in the case of .

Figure 5: Effect of mismatch on synchronization for
Figure 6: Effect of mismatch on synchronization for

3. Experiments with a Simple Electrical Circuit

3.1. Circuit Description and Equations

The Lorenz system, being one of the most important paradigms of chaos, has inspired many attempts to make a physical system representing its equations, mainly in the form of an electrical circuit. Some authors have proposed replacing the cross-products of variables by discontinuities (switching circuits) as in [12, 13], or by continuous piecewise linear resistors [14], thus resulting in very simple and practical circuits, although not truly described by the Lorenz equations. More accurate realization, though more complex, is by the analog computer approach using smooth cross-product functions, as in [15], which employs 10 integrated circuits (2 multipliers and 8 op amps) and 23 passive components, therefore a total of 33 circuit components for the Lorenz 3D circuit. In the present work we are proposing a simpler easy-to-build circuit with smooth functions, aiming at encouraging more experimental approaches on hyperchaos investigation, even by those researchers not trained in electronics. As stated in [3], to facilitate experiments with coupled chaotic oscillators the circuit is required to exhibit chaos in a large range of parameters in order that the coupling will not destroy the attractor, and it is needed to be simple enough so that several practically identical oscillators can be easily constructed. The simplest possible smooth Lorenz 3D circuit appeared in [16], using only 2 integrated circuits (2 multipliers) and 7 passive components (a total of 9 components, thus about 70% smaller than the circuit by Cuomo et al. in [15]), as shown in Figure 7(a). The good performance of this circuit is illustrated in Figure 7(b), which shows the experimental attractor and examples of single-cusp and double-cusp Lorenz maps (displayed in-line by the circuit, via a Poincaré-section circuitry). In the present work we extended to 4D that simplest circuit by adding 2 op amps and 6 passive components, obtaining the hyperchaotic Lorenz circuit shown in Figure 8(a), redrawn in Figure 8(b) using circuit theoretic symbols. The following equations describe the circuit: where we used the multipliers transfer function , where on the input of the first multiplier (on the left side in Figure 8(a)), and on the input of the second one. In deriving (3.1) we assumed and . Now, by defining the new variables and parameters we obtain (2.2). Therefore, the proposed circuit realizes, exactly, the Lorenz hyperchaotic system given by (2.2), obviously with some usual practical restrictions imposed by parasitic effects and finite bandwidth, slew-rate, excursion range, and so forth.

Figure 7: (a) Simplest circuit for the 3D standard Lorenz equations. (b) Experimental results from the circuit of Figure 7(a): butterfly attractor on the plane ; Lorenz maps showing single and double cusps . The straight line in each picture is given by .
Figure 8: (a) Simplest circuit for the 4D hyperchaotic Lorenz system of (2.2). (b) Equivalent Lorenz hyperchaotic circuit.

3.2. Some Experimental Results

For our experiments, two circuits following the schematic diagram of Figure 8(a) were constructed, each one using two AD633 multipliers and two LM351 op amps, all powered with 15 V, and the following passive components: , , , , , , , , , , . Using (3.2), the corresponding parameters of (2.2) are , , We worked with , giving (Note: we verified that with an independent DC power source it is possible to use values up to 30 V, or , without waveform clipping.) The six projections of the experimental hyperchaotic attractor are shown in Figure 9(b); the calculated attractor is shown in Figure 9(a). For the experiments on synchronization, bidirectional diffusive coupling can be obtained simply by connecting a resistor linking the capacitor of the first circuit with the capacitor of the second circuit, since in this work we have tested only , or , coupling. For unidirectional coupling a voltage follower is added in series with the resistor , as sketched in Figure 10. The coupling strength is given by the relation , as can be easily verified by adding the term to the first of (3.1). Note that each variable of the second circuit is represented by the same symbol as the corresponding one of the first circuit, but with an uppercase prime. The four projections of the trajectories near the synchronization hyperplane are shown in Figures 11(a) and 11(b) for several values of the coupling resistor . The waveforms and are shown in Figure 12. All these results refer to unidirectional coupling through the variables and .

Figure 9: (a) Projections of a calculated hyperchaotic attractor described by (2.2) for , and From left to right, top: ; ; ; bottom: (b) Experimental hyperchaotic Lorenz attractor generated by the circuit of Figure 8, as projected on the oscilloscope screen. From left to right, top: bottom: The same parameter values as those of Figure 9(a)—see text.
Figure 10: Coupling schemes: (a) bidirectional; (b) unidirectional.
Figure 11: (a) Experimental trajectories near the synchronization hyperplane for several values of the unidirectional coupling resistor. A: 0 ; B: 100 ; C: 220 . (b) Continuation of Figure 11(a). D: 470 ; E: 1.0 ; F: 2.2 .
Figure 12: Experimental waveforms and corresponding to the first column of Figure 11.

4. Conclusion

In this work the synchronizing properties of diffusively coupled hyperchaotic Lorenz 4D systems, described by (2.2), have been studied both numerically and experimentally. The numerical investigation was realized by calculating the transverse Lyapunov exponents as a function of the coupling strength , and also by visually inspecting the phase space trajectories near the synchronization hyperplane. We concluded that using a single coupling variable, either or (but neither nor ), guarantees stable synchronization. Although and are not good choices for the single-variable scheme, double coupling with both and easily provides synchronization. We also verified that a small degree of parameter mismatch seems tolerable. For the laboratory work a very simple electrical circuit described by the Lorenz 4D system was proposed and described here for the first time. The experiments confirmed the qualitative behavior predicted by the numerical approach.


This work was supported by the Brazilian agency FAPESP. The author thanks the reviewers for the useful comments and suggestions.


  1. L. M. Pecora and T. L. Carroll, “Synchronization in chaotic systems,” Physical Review Letters, vol. 64, no. 8, pp. 821–824, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  2. L. M. Pecora and T. L. Carroll, “Driving systems with chaotic signals,” Physical Review A, vol. 44, no. 4, pp. 2374–2383, 1991. View at Publisher · View at Google Scholar
  3. J. F. Heagy, T. L. Carroll, and L. M. Pecora, “Synchronous chaos in coupled oscillator systems,” Physical Review E, vol. 50, no. 3, pp. 1874–1885, 1994. View at Publisher · View at Google Scholar
  4. L. M. Pecora, T. L. Carroll, G. A. Johnson, D. J. Mar, and J. F. Heagy, “Fundamentals of synchronization in chaotic systems, concepts, and applications,” Chaos, vol. 7, no. 4, pp. 520–543, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  5. K. Pyragas, “Predictable chaos in slightly perturbed unpredictable chaotic systems,” Physics Letters A, vol. 181, no. 3, pp. 203–210, 1993. View at Publisher · View at Google Scholar
  6. J. H. Peng, E. J. Ding, M. Ding, and W. Yang, “Synchronizing hyperchaos with a scalar transmitted signal,” Physical Review Letters, vol. 76, no. 6, pp. 904–907, 1996. View at Publisher · View at Google Scholar
  7. A. Tamasevicious and A. Cenis, “Synchronizing hyperchaos with a single variable,” Physical Review E, vol. 55, no. 1, pp. 297–299, 1997. View at Publisher · View at Google Scholar
  8. A. Wolf, J. B. Swift, H. L. Swinney, and J. A. Vastano, “Determining Lyapunov exponents from a time series,” Physica D, vol. 16, no. 3, pp. 285–317, 1985. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  9. J. N. Blakely and D. J. Gauthier, “Attractor bubbling in coupled hyperchaotic oscillators,” International Journal of Bifurcation and Chaos, vol. 10, no. 4, pp. 835–847, 2000. View at Google Scholar
  10. R. Barboza, “Dynamics of a hyperchaotic Lorenz system,” International Journal of Bifurcation and Chaos, vol. 17, no. 12, pp. 4285–4294, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  11. C. Sparrow, The Lorenz Equations: Bifurcations, Chaos, and Strange Attractors, vol. 41 of Applied Mathematical Sciences, Springer, New York, NY, USA, 1982. View at Zentralblatt MATH · View at MathSciNet
  12. S. Özoğuz, A. S. Elwakil, and M. P. Kennedy, “Experimental verification of the butterfly attractor in a modified Lorenz system,” International Journal of Bifurcation and Chaos, vol. 12, no. 7, pp. 1627–1632, 2002. View at Publisher · View at Google Scholar
  13. E. H. Baghious and P. Jarry, ““Lorenz attractor” from differential equations with piecewise-linear terms,” International Journal of Bifurcation and Chaos, vol. 3, no. 1, pp. 201–210, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  14. R. Tokunaga, T. Matsumoto, L. O. Chua, and S. Miyama, “The piecewise-linear Lorenz circuit is chaotic in the sense of Shilnikov,” IEEE Transactions on Circuits and Systems, vol. 37, no. 6, pp. 766–786, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  15. K. M. Cuomo, A. V. Oppenheim, and S. H. Strogatz, “Synchronization of Lorenz-based chaotic circuits with applications to communications,” IEEE Transactions on Circuits and Systems II, vol. 40, no. 10, pp. 626–633, 1993. View at Publisher · View at Google Scholar
  16. R. Barboza, “Experiments on Lorenz system,” in Proceedings of the International Symposium on Nonlinear Theory and Its Applications (NOLTA '04), vol. 1, pp. 529–532, Fukuoka, Japan, November 2004.