Abstract and Applied Analysis

Volume 2014 (2014), Article ID 213694, 17 pages

http://dx.doi.org/10.1155/2014/213694

Research Article

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

Cheng-Hsiung Yang and Cheng-Lin Wu

Graduate Institute of Automation and Control, National Taiwan University of Science and Technology, No. 43, Section 4, Keelung Road, Taipei 106, Taiwan

Received 2 April 2014; Accepted 25 May 2014; Published 7 August 2014

Academic Editor: Xiao He

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.

Abstract.

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 [7–12] with uncertain chaotic parameters strategy by GYC partial region stability theory is proposed [13–16]. 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.

Consider

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.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 1: Phase portraits and Poincaré maps for (a) chaotic, (b) period-1, (c) period-2, (d) period-4, (e) period-3, (f) period-6, and (g) period-12, respectively, of the four-dimensional Lorenz-Stenflo system.

(a)

(b)

Figure 2: Three dimensions phase portrait of the four-dimensional Lorenz-Stenflo system and its projection.

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.

Table 1: Eigenvalues and stability of equilibrium points.

Equilibrium point
Eigenvalues
Stable/Unstable

= 3.7, = −1, = −0.7, = −3.7
Unstable

= 2.152 − 6.27i, = −12.5777, = 2.152 + 6.27i, = −0.826
Unstable

= 2.1489 + 6.50198i, = −12.5416, = 2.1489 − 6.50198i, = −0.8561
Unstable

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.

(a)

(b)

Figure 3: Lyapunov exponents and bifurcation diagrams of the four-dimensional Lorenz-Stenflo system with parameters , , and .

(a)

(b)

Figure 4: Lyapunov exponents and bifurcation diagrams of the four-dimensional Lorenz-Stenflo system with parameters , , and .

(a)

(b)

Figure 5: Lyapunov exponents and bifurcation diagrams of the four-dimensional Lorenz-Stenflo system with parameters , , and .

(a)

(b)

Figure 6: Lyapunov exponents and bifurcation diagrams of the four-dimensional Lorenz-Stenflo system with parameters , , and .

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.

(a)

(b)

(c)

(d)

Figure 7: Lyapunov dimension diagrams of the four-dimensional Lorenz-Stenflo system.

2.6. Electronic Circuit Implementation

The four-dimensional Lorenz-Stenflo system can be modeled by an electronic circuit [18–21]. 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

Figure 8: Chaotic circuit of the four-dimensional Lorenz-Stenflo system.

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.

Table 2: List of apparatus and components in the circuit experiment.

Apparatus/Component
Model
Quantity

Multiplier
AD633
2

OP-Amplifier
LF412
4

Capacitor
2000 pF~0.47 uF
4

Resister
1 k~1.2 M
19

270 kΩ
3

38.461 kΩ
1

667 kΩ
1

1429 kΩ
1

100 kΩ
8

1 MΩ
4

Variable resistor
100 k
2

DC power supply
GPS-3303
1

(a)

(b)

(c)

(d)

(e)

(f)

Figure 9: Oscilloscope graphics for voltage signals outputs of the four-dimensional Lorenz-Stenflo system: (a) , (b) , (c) , (d) , (e) , and (f) .

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.

Figure 10: Bifurcation diagram with phase portraits of the four-dimensional Lorenz-Stenflo system.

Figure 11: Experimental setup for realization of the four-dimensional Lorenz-Stenflo system.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 12: Realization and simulation circuit compare results for voltage signals outputs of the four-dimensional Lorenz-Stenflo system: (a) , (b) , (c) , (d) , (e) , and (f) .

Figure 13: Realization circuit of the four-dimensional Lorenz-Stenflo system.

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.

Consider

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.

Figure 14: Time histories of the master system and slave system.

Figure 15: Time histories of error states.

Figure 16: Time histories of estimated parameters.

(a)

(b)

(c)

(d)

Figure 17: Time histories of parameters differences.

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.

Acknowledgment

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

References

Z. Li, Y. Soh, and C. Wen, Switched and Impulsive Systems, Analysis, Design, and Applications, vol. 313, Springer, Berlin, Germany, 2005.

P. Galajda and D. Kocur, “Chua's circuit in spread spectrum communication systems,” Radio Engineering, vol. 11, no. 2, pp. 6–10, 2002.

L. Kocarev, Z. Galias, and S. Lian, Eds., Intelligent Computing Based on Chaos, vol. 184 of Studies in Computational Intelligence, Springer, Berlin, Germany, 2009.

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.

S. H. Strogatz, Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering, Perseus Books, Cambridge, Mass, USA, 1994.

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.

E. N. Lorenz, “Deterministic nonperiodic flow,” Journal of the Atmospheric Sciences, vol. 20, pp. 130–141, 1963.

L. Stenflo, “Generalized Lorenz equations for acoustic-gravity waves in the atmosphere,” Physica Scripta, vol. 53, no. 1, pp. 83–84, 1996.

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.

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.

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.

W. Xiang, “Equilibrium points and bifurcation control for Lorenz-Stenflo system,” ICIC Express Letters, vol. 3, no. 1, pp. 61–66, 2009.

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.

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.

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.

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.

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.

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.

Y. Liu, “Circuit implementation and finite-time synchronization of the 4D Rabinovich hyperchaotic system,” Nonlinear Dynamics, vol. 67, no. 1, pp. 89–96, 2012.

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.

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.