Abstract and Applied Analysis

Abstract and Applied Analysis / 2014 / Article
Special Issue

Nonlinear Analysis of Dynamical Complex Networks 2014

View this Special Issue

Research Article | Open Access

Volume 2014 |Article ID 213694 | 17 pages | https://doi.org/10.1155/2014/213694

Nonlinear Dynamic Analysis and Synchronization of Four-Dimensional Lorenz-Stenflo System and Its Circuit Experimental Implementation

Academic Editor: Xiao He
Received02 Apr 2014
Accepted25 May 2014
Published07 Aug 2014


Recently many chaotic systems’ circuits are designed to generate phenomenon of chaos signals. The ability to synchronize chaotic circuits opens a great number of ways to use them in application signals masking. In this paper, first a new nonlinear chaotic dynamical system had be design, analyze and build circuit. Second, using GYC, partial region stability theory is applied to adaptive control for two identical chaotic systems with uncertain parameters. The results of numerical simulation are performed to verify examples of the proposed nonlinear controllers.

1. Introduction

Chaotic systems provide a rich mechanism for signal design and generation, with potential applications to communications and signal processing. Because chaotic signals are typically broadband, noise-like, and difficult to predict, they can be used in various contexts for masking information bearing waveforms. They can also be used as modulating waveforms in spread spectrum systems.

This can be useful in many practical circumstances like securing communication channels [1], masking signals [2], spreading data sequence, or for generating random signals [3]. Of course, chaos is an unwanted phenomenon in many situations. The essential problem which is faced here is how to distinguish such motion from long transient behavior. A spectrum called Lyapunov exponents (LE) is often used as a quantifier [4] of chaotic motion or regular motion. These real numbers measure the average ratio of the exponential separation of the two neighborhood trajectories. For chaos it is necessary to have at least one positive of Lyapunov exponents; the second one represents a direction of the flow and must converge to zero. The last one must be negative with the largest absolute value, since the flow is dissipative. Circuits of nonlinear dynamic system provide an excellent tool for the study of chaotic behavior. Some of these circuits treat time as a discrete variable, employing sample-and-hold subcircuits and analog multipliers to model iterated maps such as the logistic map [5, 6].

In this paper, an adaptive generalized synchronization of the four-dimensional Lorenz-Stenflo system [712] with uncertain chaotic parameters strategy by GYC partial region stability theory is proposed [1316]. By using the GYC partial region stability theory the Lyapunov function is a simple linear homogeneous function of states and the structure of controllers are simpler and have less simulation error because they are in lower order than that of traditional controllers design.

This paper is organized as follows. In Section 2, the four-dimensional Lorenz-Stenflo system analyzes characteristics of nonlinear dynamic behavior, construction realization circuit of the four-dimensional Lorenz-Stenflo system. In Section 3, adaptive generalized chaos synchronization with uncertain parameters strategy by GCY partial region stability theory is proposed. In Section 4, numerical results of adaptive generalized synchronization of the four-dimensional Lorenz-Stenflo system with uncertain chaotic parameters using GYC partial region stability. Finally, some concluding remarks are given in Section 5.

2. Chaotic Dynamic Analysis of a Four-Dimensional Lorenz-Stenflo System

2.1. Phase Portraits and Poincaré Map

The Lorenz-Stenflo system is described as [7, 8] where is Prandtl number, is rotation number, is Rayleigh number, and is geometric parameter.

The of Lorenz-Stenflo system from first item changes to third item, new four-dimensional Lorenz-Stenflo system.


The Lorenz-Stenflo system has many dynamical behaviors that are one chaotic motion and six different periodic motions as shown in Figure 1. To have a better understanding of the four-dimensional Lorenz-Stenflo system, we use an innovative technique to present the three-dimensional phase portraits shown in Figure 2. These figures plot three dimensions phase portraits and their projection simultaneously. Thus, we can comprehend the four-dimensional Lorenz-Stenflo system much easier.

2.2. Equilibrium Analysis

The equilibrium points of the system are shown as follows: The parameters of four-dimensional Lorenz-Stenflo system are , , , and , so the equilibrium points become

For the equilibrium point , the Jacobian matrix of the system is The eigenvalues of the Jacobian matrix below are obvious:

Then we linearize the system about the nonzero equilibrium points and obtaine the Jacobian matrix: The eigenvalues are

By letting , the eigenvalues of the matrix at each equilibrium point can be obtained as shown in Table 1. The eigenvalues of each equilibrium point are the same. And all the equilibrium points are unstable, since at least one eigenvalue has positive real part for each equilibrium point.

Equilibrium pointEigenvaluesStable/Unstable

= 3.7, = −1, = −0.7, = −3.7Unstable
= 2.152 − 6.27i, = −12.5777, = 2.152 + 6.27i, = −0.826Unstable
= 2.1489 + 6.50198i, = −12.5416, = 2.1489 − 6.50198i, = −0.8561Unstable

2.3. Divergence Analysis

The divergence represents the volume density of the outward flux of a vector field. The divergence of a system implies whether this system dissipative or not. For the four-dimensional Lorenz-Stenflo system, its divergence is where .

The parameters are taken as , , , and ; hence, the divergence of the system becomes So, (2) is dissipative, with an exponential contraction rate:

So, a volume element is contracted by the flow into a volume element in time. All orbits of (2) are ultimately confined to a set having zero volume and asymptotic motion settles onto an attractor. The four-dimensional Lorenz-Stenflo system has bund system.

2.4. Bifurcation Diagram

For numerical analysis, all of the initial conditions are , , , and , in this and next sections. The bifurcation diagrams of the four-dimensional Lorenz-Stenflo system are plotted in Figure 3(a) to Figure 6(a). The parameters , , , and in (2) will be varied, respectively, for plotting the bifurcation diagrams. From the bifurcation diagrams, we can learn more about the behaviors of the four-dimensional Lorenz-Stenflo system.

2.5. Lyapunov Exponent and Lyapunov Dimension

The Lyapunov exponent may be used to measure the sensitive dependence on initial conditions. Once the Lyapunov exponent is greater than zero at specific value of parameters, it is said to be chaos, otherwise, periodic solution. The Lyapunov exponents diagrams of the four-dimensional Lorenz-Stenflo system are shown in Figure 3(b) to Figure 6(b). In order to verify the chaotic behavior of the system, Lyapunov exponents diagrams and bifurcation diagrams can be compared together.

The Lyapunov dimension has been proposed by Frederickson et al. [17] to measure the complexity of the attractor. The relation between Lyapunov dimension and Lyapunov exponent is where is the largest Lyapunov exponent and is defined by the condition as The Lyapunov dimension diagrams are plotted in Figure 7.

2.6. Electronic Circuit Implementation

The four-dimensional Lorenz-Stenflo system can be modeled by an electronic circuit [1821]. The main electronic components in the circuit are integrated amplifiers and inverting amplifiers. In this circuit simulation, we use electronic circuit simulation software named Multisim. All the parameters and initial conditions are the same as numerical simulation. But here we modify the value of capacitors to decrease the output voltage within a range from −12.0 to 12.0, and all the state values reduced tenfold in circuit simulation. The chaotic circuit configuration diagram is shown in Figure 8, and the governing integral equation of the circuit can be written as

Components of the chaotic system circuit are chosen to be , , , , , , and . The voltage output signals of the chaotic circuit simulated in Multisim software are shown in Figure 9. In addition, the real circuit of the four-dimensional Lorenz-Stenflo system was constructed on a breadboard shown in Figure 13. The apparatus and the components used in the circuit are listed in Table 2. The experimental setup is shown in Figure 11 and the experimental results are shown in Figure 12.


Capacitor2000 pF~0.47 uF4
Resister1 k~1.2 M19
270 kΩ3
38.461 kΩ1
667 kΩ1
1429 kΩ1
100 kΩ8
1 MΩ4
Variable resistor100 k2
DC power supplyGPS-33031

In order to have a better concept about chaotic motions and periodic motions, we present a four-dimensional diagram, which is a combination of the phase portraits of the circuits simulation and the bifurcation diagram of variable parameter a, plotted in Figure 10.

3. Adaptive Generalized Synchronization of the Four-Dimensional Lorenz-Stenflo System with Uncertain Chaotic Parameters via GYC Partial Region Stability Theory Strategy

Consider the master system and the slave system where and denote the master state vector and slave state vector, respectively, is nonlinear vector functions, is uncertain chaotic parameter in , is estimates of uncertain chaotic parameter in , and is a control input vector.

Our goal is to design a controller so that the state vector of the slave system (16) asymptotically approaches the state vector of the master system (15) plus a given vector function which is either a regular or a chaotic function of time, and finally the synchronization will be accomplished in the sense that the limit of the error vector approaches zero: where . The are positive constants which make the error dynamics always happens in the first quadrant.

From (17) we have Introducing (15) and (16) in (18), we get

A Lyapunov function is chosen as a positive definite function in the first quadrant.

Consider where . The are positive constants which make the update laws always happen in the first quadrant.

Its derivative along the solution of (20) is where and are chosen so that , and are negative constants, and is a negative definite function of and . When the generalized synchronization is obtained.

By using the GYC partial region stability theory, the Lyapunov function is easier to find, since the terms of first degree can be used to construct the definite Lyapunov function and the controller can be designed in lower order.

4. Numerical Results of Adaptive Generalized Synchronization of the Four-Dimensional Lorenz-Stenflo System with Uncertain Chaotic Parameters via GYC Partial Region Stability Theory

The master system is where , and are uncertain chaotic parameters. Equation (23) is called four-dimensional Lorenz-Stenflo system with uncertain chaotic parameters.

Consider where , and are arbitrary positive constant and .

The chaotic signals system is where , , , and .

The slave system is where , and are estimates of uncertain chaotic parameters, respectively.

The nonlinear controllers, , and , are added to each equation in (26), respectively.


The given functional system for generalized synchronization is also the four-dimensional Lorenz-Stenflo system but with different initial conditions.

Consider The initial value of the states of the master system, of the slave system, of the chaotic signals system, and of the given functional system are taken as 0.2, 0.3, 0.7, and 0.5, respectively.

The generalized synchronization error functions are The addition of constant makes the error states always greater than zero. As Figure 15 shows, we can confirm that the error states are always greater than zero. From (29), we get the following error dynamics: where ,  , and . The addition of constant makes the update laws always greater than zero. As Figure 17 shows, we can confirm that the update laws are always greater than zero.

Choose a Lyapunov function in the form of a positive definite function: and its time derivative is

We choose the update laws for those uncertain parameters as

The initial values of estimates for uncertain parameters are . Through (32) and (33), the appropriate controllers can be designed as

Substituting (33) and (34) into (32), we obtain which is a negative-definite function of and. The numerical simulation results are shown in Figures 14, 15, 16, and 17.

5. Conclusions

In summary, the paper has studied the nonlinear dynamical behaviors of the new chaotic attractor, including some basic dynamical properties, bifurcations, Lyapunov exponents, Lyapunov dimension, periodic motion, and routes to chaos motion. Also, the four-dimensional Lorenz-Stenflo system based realization circuit has been constructed to implement various wing butterfly attractors. Furthermore, by electronic circuit implementation of the proposed four-dimensional Lorenz-Stenflo system, it is shown that the chaotic attractors do physically exist. A new strategy, the GYC partial region stability theory, is proposed to achieve adaptive chaos synchronization with uncertain chaotic parameters.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


This research was supported by the National Science Council, Republic of China, under Grant no. NSC 102-2221-E-011-034.


  1. Z. Li, Y. Soh, and C. Wen, Switched and Impulsive Systems, Analysis, Design, and Applications, vol. 313, Springer, Berlin, Germany, 2005. View at: MathSciNet
  2. P. Galajda and D. Kocur, “Chua's circuit in spread spectrum communication systems,” Radio Engineering, vol. 11, no. 2, pp. 6–10, 2002. View at: Google Scholar
  3. L. Kocarev, Z. Galias, and S. Lian, Eds., Intelligent Computing Based on Chaos, vol. 184 of Studies in Computational Intelligence, Springer, Berlin, Germany, 2009. View at: Publisher Site | MathSciNet
  4. K. Grygiel and P. Szlachetka, “Lyapunov exponents analysis of autonomous and nonautonomous sets of ordinary differential equations,” Acta Physica Polonica B, vol. 26, no. 8, pp. 1321–1331, 1995. View at: Google Scholar | Zentralblatt MATH | MathSciNet
  5. S. H. Strogatz, Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering, Perseus Books, Cambridge, Mass, USA, 1994.
  6. K. Kiers, D. Schmidt, and J. C. Sprott, “Precision measurements of a simple chaotic circuit,” American Journal of Physics, vol. 72, no. 4, pp. 503–509, 2004. View at: Publisher Site | Google Scholar
  7. E. N. Lorenz, “Deterministic nonperiodic flow,” Journal of the Atmospheric Sciences, vol. 20, pp. 130–141, 1963. View at: Google Scholar
  8. L. Stenflo, “Generalized Lorenz equations for acoustic-gravity waves in the atmosphere,” Physica Scripta, vol. 53, no. 1, pp. 83–84, 1996. View at: Google Scholar
  9. J. C. Xavier and P. C. Rech, “Regular and chaotic dynamics of the Lorenz-Stenflo system,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 20, no. 1, pp. 145–152, 2010. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  10. P. Wang, D. Li, and Q. Hu, “Bounds of the hyper-chaotic Lorenz-Stenflo system,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 9, pp. 2514–2520, 2010. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  11. H. Li, X. Liao, and X. Lei, “Two fuzzy control schemes for Lorenz-Stenflo chaotic system,” Journal of Vibration and Control, vol. 18, no. 11, pp. 1675–1682, 2012. View at: Publisher Site | Google Scholar | MathSciNet
  12. W. Xiang, “Equilibrium points and bifurcation control for Lorenz-Stenflo system,” ICIC Express Letters, vol. 3, no. 1, pp. 61–66, 2009. View at: Google Scholar
  13. C.-H. Yang, “Chaos hybrid generalized synchronization of liu-chen system by GYC partial region stability theory,” Journal of Computational and Theoretical Nanoscience, vol. 10, no. 4, pp. 825–831, 2013. View at: Publisher Site | Google Scholar
  14. Z.-M. Ge and C.-H. Yang, “The generalized synchronization of a Quantum-CNN chaotic oscillator with different order systems,” Chaos, Solitons and Fractals, vol. 35, no. 5, pp. 980–990, 2008. View at: Publisher Site | Google Scholar
  15. Z.-M. Ge, C.-W. Yao, and H.-K. Chen, “Stability on partial region in dynamics,” Journal of the Chinese Society of Mechanical Engineers, vol. 15, no. 2, pp. 140–151, 1994. View at: Google Scholar
  16. Z.-M. Ge, J.-K. Yu, and H.-K. Chen, “Three asymptotical stability theorems on partial region with applications,” Japanse Journal of Applied Physics, vol. 37, article 2762, 1998. View at: Publisher Site | Google Scholar
  17. P. Frederickson, J.-L. Kaplan, E.-D. Yorke, and J.-A. Yorke, “The Liapunov dimension of strange attractors,” Journal of Differential Equations, vol. 49, no. 2, pp. 185–207, 1983. View at: Publisher Site | Google Scholar | MathSciNet
  18. X. Wang and G. A. Chen, “A chaotic system with only one stable equilibrium,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 3, pp. 1264–1272, 2012. View at: Publisher Site | Google Scholar | MathSciNet
  19. Y. Liu, “Circuit implementation and finite-time synchronization of the 4D Rabinovich hyperchaotic system,” Nonlinear Dynamics, vol. 67, no. 1, pp. 89–96, 2012. View at: Publisher Site | Google Scholar
  20. A. G. Radwan, A. M. Soliman, and A. L. El-Sedeek, “MOS realization of the double-scroll-like chaotic equation’’,” IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications, vol. 50, no. 2, pp. 285–288, 2003. View at: Publisher Site | Google Scholar
  21. A. Abooee, H. A. Yaghini-Bonabi, and M. R. Jahed-Motlagh, “Analysis and circuitry realization of a novel three-dimensional chaotic system,” Communications in Nonlinear Science and Numerical Simulation, vol. 18, no. 5, pp. 1235–1245, 2013. View at: Publisher Site | Google Scholar | MathSciNet

Copyright © 2014 Cheng-Hsiung Yang and Cheng-Lin Wu. 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.

899 Views | 639 Downloads | 6 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19.