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

and Applied Analysis 3 0 10 20 0 20 0 10 20 30 40 50


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][8][9][10][11][12] with uncertain chaotic parameters strategy by GYC partial region stability theory is proposed [13][14][15][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.

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  4 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.

Equilibrium Analysis.
The equilibrium points of the system are shown as follows: 0 (0, 0, 0, 0) ) . ( The parameters of four-dimensional Lorenz-Stenflo system are  = 3.7,  = 1.5,  = 26, and  = 0.7, so the equilibrium points become  0 (0, 0, 0, 0) For the equilibrium point  0 , 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 det |  − | = 0, 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
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  = 3.7,  = 1.5,  = 26, and  = 0.7; hence, the divergence of the system becomes div Φ = 3.7 − 1 + 0.7 − 3.7 = −9.1. ( So, ( 2) is dissipative, with an exponential contraction rate: So, a volume element Π 0 is contracted by the flow into a volume element Φ 0  −(2++1) in time.All orbits of (2) are ultimately confined to a set having zero volume and asymptotic motion settles onto an attractor.The fourdimensional Lorenz-Stenflo system has bund system.

Bifurcation Diagram.
For numerical analysis, all of the initial conditions are  1 (0) = 20,  2 (0) = 20,  3 (0) = 20, and  4 (0) = 20, 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.

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 [18][19][20][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 Components of the chaotic system circuit are chosen to be   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.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.
A Lyapunov function (, Ã) is chosen as a positive definite function in the first quadrant. Consider where () and Ȧ are chosen so that V =  +  Ã, and  are negative constants, and V is a negative definite function of  1 ,  2 , . . .,   and Ã1 , Ã2 , . . ., Ã .When lim 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.
The slave system is where Â1 (), Â2 (), Â3 (), and Â4 () are estimates of uncertain chaotic parameters, respectively.The nonlinear controllers,  1 ,  2 ,  3 , and  4 , 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 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 its time derivative is We choose the update laws for those uncertain parameters as The initial values of estimates for uncertain parameters are Â1 (0) = Â2 (0) = Â3 (0) = Â4 (0) = 0. Through (32) and (33), the appropriate controllers can be designed as

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.

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

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

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

Figure 16 :
Figure 16: Time histories of estimated parameters.

Figure 17 :
Figure 17: Time histories of parameters differences.

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