Control and Synchronization of Chaotic and Hyperchaotic Lorenz Systems via Extended Backstepping Techniques

Wepropose novel controllers for stabilization and tracking of chaotic andhyperchaotic Lorenz systems using extended backstepping techniques. Based on the proposed approach, generalized weighted controllers were designed to control chaotic behaviour as well as to achieve synchronization in chaotic and hyperchaotic Lorenz systems. The effectiveness and feasibility of the proposed weighted controllers were verified numerically and showed their robustness against noise.


Introduction
Chaos theory has found application in many areas of studies; these include mathematics, physics, biology, engineering, economics, and politics [1][2][3].One of the most successful applications of chaos theory has been in ecology, where dynamical systems have been used to show how population growth under density dependence can lead to chaotic dynamics.Chaos theory is also currently being applied to medical studies of epilepsy, specifically to the prediction of seemingly random seizures by observing initial conditions [4].Furthermore, a related field of physics called quantum chaos theory investigates the relationship between chaos and quantum mechanics [5].In addition, another field called relativistic chaos [6] has emerged to describe systems that follow the laws of general relativity.
Chaotic phenomenon could be beneficial in some applications; however, it is undesirable in many engineering and other physical applications and should therefore be controlled in order to improve the system performance.Chaos control is concerned with using some designed control inputs to modify the characteristics of a parameterized nonlinear system so that the system becomes stable at a chosen position or tracks a desired trajectory [7].Several techniques have been deviced for chaos control but mostly are for development of two basic approaches: the OGY (Ott, Grebogi, and Yorke) method and Pyragas continuous control.Both methods require a previous determination of the unstable periodic orbits of the chaotic system before the controlling algorithm can be designed.Experimental control of chaos by one or both of these methods has been achieved in a variety of systems, including turbulent fluids, oscillating chemical reactions, magneto-mechanical oscillators, and cardiac tissues [8].
Since the idea of synchronization of chaotic systems was proposed by Pecora and Carroll in 1990, [9], chaos synchronization has received an increasing attention due to its theoretical challenge and its potential applications in secure communications, chemical reactions, biological systems, information science, and plasma technologies [10].Chaos synchronization involves the coupling or forcing of two systems so that both systems achieve identical dynamics asymptotically with time.A wide variety of approaches have been proposed to achieve chaos synchronization in the coupled/forced chaotic systems such as linear state error feedback method [11], time-delay feedback method [12], active control approach [13,14], impulsive method [15], adaptive control [16,17], and backstepping approach [18].The backstepping approach has been widely applied to achieve synchronization in chaotic and hyperchaotic systems.It has the ability to achieve global stability tracking transient performance for a board class of strict-feedback nonlinear system (see [7] and the references therein).Also the method requires less control effort in comparison to the differential geometric method [19].However, control function designed via this method has been shown to be difficult for practical implementation because of its complexity and high signal strength.
In order to reduce the complexity in controllers and energy consumption, the present work is focused on the design of generalized weighted controllers for stabilization and tracking of chaotic and hyperchaotic systems using extended backstepping method.Based on this approach, controller complexity is minimized and controller is singularity free from nonlinear quadratic terms [7,11].The proposed approach is suitable for practical implementation and this shall be a focus in future work.
To the best of our knowledge, this approach has not been used to control as well as synchronize chaotic and hyperchaotic Lorenz systems.Also, the approach may be suitable for practical implementation in some real systems.

Systems Description
In this work, two nonlinear systems are studied, namely, chaotic Lorenz system and hyperchaotic Lorenz system.In the first instance, we consider the chaotic Lorenz system of the form where , , and  are the constant parameters and , , and  are the state variables of the system.The dynamics of system (1) are chaotic for parameter values with the chaotic attractor shown in Figure 1.
The second system that we studied is the hyperchaotic Lorenz system describable by the following differential equations: where , , ,  1 , and  2 are the control parameters and , , , and  are the state variables of the system under investigation.

Extended Controller for Tracking Control of Chaotic System
3.1.Design of Recursive Backstepping Controller.Let us recast system (1) as follows: where   (),  = 1, 2, 3, are the control inputs to be determined such that the state variables , , and  of system (3) can take the desired values   ,   , and   , respectively.We define the error states between the state variable and the desired values as To design general control functions   (),  = 1, 2, 3, that can control (3) to track any trajectory () that is a smooth function of time, we let where   ,  = 1, 2, 3, are the arbitrary control parameter to be chosen appropriately.By substituting ( 5) into (4) and differentiating the resulting equation, we have To stabilize the error system (6), we consider a Lyapunov function of the form where   ,   , and   are positive constants coefficients and obtain its time derivative as To satisfy the condition for asymptotic stability of the error dynamics ( 6) necessary for tracking, we substitute ( 6) into ( 8) and choose   () ( = 1, 2, 3) such that V = − ∑    2  < 0,  = , , and , as follows: We observed that system (3) is effectively controlled with  1 =  3 = 1 and  2 = 0 which simplify the controllers in (9) to By introducing weights to the above controllers, we obtain the following: where   ,  = 1, 2, 3, are the weights added to the usual controllers.
The results obtained show that the state variables move chaotically with time when the controller is switched off and when the controller is activated at time  = 50 the state variables are controlled to track the desired functions () = 35 cos 0.57.The result is shown in Figure 3.With the value of   = 0.5 (Figure 4) where  = 1, 2, 3 we observed that the controllers ( 10) and ( 11) are capable of controlling the dynamics of the chaotic system (3) to track any desired smooth function, (), and to stabilize it at any position  (case () = ) and when  = 0 the system becomes stabilized at the equilibrium point (0, 0, 0).The results in Figure 4 showed that the proposed extended backstepping techniques reduced the control strength and complexity of the designed controllers by 50%, thereby producing economic controllers with low energy consumption which are suitable for practical implementation.

Extended Controller for Tracking Control of the Hyperchaotic System
where   ,  = 1, 2, 3, 4, are the control inputs to be determined later such that the state variables , , , and  of system (13) can take the desired values   ,   ,   , and   , respectively.We define the error states between the state variable and the desired values as To design general control functions,   (), that can control system (12) to track any trajectory () that is a smooth function of time, we let The parameters of the system are as shown in Figure 1.
where   ( = 1, 2, . . .6) are arbitrary control parameters to be chosen appropriately.Substitution of ( 14) into ( 13) and carrying out some algebraic manipulations yield To stabilize the error system (15), we consider a Lyapunov function of the form where  ,   ,   and   , are positive constants coefficients and obtain its time derivative as To satisfy the condition for asymptotic stability of the error dynamics ( 6) necessary for tracking, we substitute (15) into (17) and choose   () ( = 1, 2, 3, 4) such that V = − ∑    2  < 0,  = , , , and , as follows: We have observed from the results of computations that system ( 12) is effectively controlled with only  1 = 1 and so we set  2 =  3 =  4 =  5 =  6 = 0 which simplify the controllers (18) to By introducing weights to the above controllers, we obtain the following: where   ,  = 1, 2, 3, 4, are the weights added to the usual controllers.

Numerical Simulation.
Using the fourth-order Runge-Kutta algorithm with initial conditions (, , , ) = (0.0, 0.5, 0.6, 0.1), a time step of 0.001, and fixing the parameter values of , , ,  1 , and  2 as in Figure 2 to ensure hyperchaotic dynamics of the state variables, we solved system (12) with the controllers   (),  = 1, 2, 3, 4, as defined in ( 19) and (20).The results obtained show that the state variables move hyperchaotically with time when the controller is deactivated and when the controller is switched on at time  = 50 the state variables are controlled to track the desired functions () = 35 cos 0.57 as shown in Figure 5.With the choice of the weights,  1 =  2 = 0.7 and  3 =  4 = 0.8 in (20), we displayed the result of our simulation in Figure 6.We observed that the controllers (19) and ( 20) are capable of controlling the dynamics of the hyperchaotic system (12) to track any desired smooth function, (), and to stabilize it at any position  (case () = ) and when  = 0 the system becomes stabilized at the equilibrium point (0, 0, 0, 0).By comparing the results in Figure 5 (recursive backstepping controllers) with those of Figure 6 (extended backstepping method), it is obvious that when the weights were added to the controllers, control strength and complexity were reduced significantly.Thus, the proposed extended backstepping controllers can produce economic controllers with low energy consumption which may be suitable for practical implementation.

Design of the Extended Backstepping Controller for Chaotic
Lorenz System.In order to derive controllers for synchronization in chaotic Lorenz system, we rewrite (1) in the form Let system (21) be the drive system and let the following be the response system: where   ,  = 1, 2, 3, are the control functions.Subtract (21) from ( 22) and use the notations The error system is obtained as The objective is to find control functions,   (), that can stabilize the error states in (24) at origin.Firstly, we stabilize the first equation in (24) by regarding  2 as a controller, choosing a Lyapunov function  1 ( 1 ) = (1/2) 2 1 , and differentiating it with respect to time to have We estimate that the controller  2 =  1 ( 1 ); then (25) can be written as; if  1 = 0 and estimative function  1 ( 1 ) = 0, then V1 = − 2 1 is negative definite (since  > 0), and hence the  1 subsystem is stabilized.The error  2 between  2 and  1 ( 1 ) is Substitution of ė 2 and  2 from (24) and (26), respectively, into the time derivative of (26) yields We now stabilize the ( 1 ,  2 ) subsystem given by (27) as follows.We choose a Lyapunov function  2 ( 1 ,  2 ) =  1 ( 1 )+ (1/2) 2 2 and obtain its time derivative as The parameters of the system are as shown in Figure 1.

Numerical Simulation.
Using the fourth-order Runge-Kutta routine with initial conditions ( 1 ,  2 ,  3 ) = (0.8 − 0.8, −1.0), ( 1 ,  2 ,  3 ) = (1.0,−1.0, −0.1), a time step of 0.001, and fixing the parameter values as in Figure 1 to ensure chaotic dynamics of the state variables, we solved systems (21) and ( 22) with the control function as defined in (37) and (38).The result obtained shows that the error state variable moves chaotically with time when the controllers are switched off and when the controllers are switched on at  = 50 as shown in Figures 7 and 8 the error state variables converge to zero, thereby guaranteeing the synchronization of systems (21) and (22).This is confirmed by the synchronization quality  given by Remarkably, we notice that the control strength and its complexity are reduced by about 70% when the proposed extended backstepping techniques is applied (Figure 8) compared with the usual backstepping approach.Thus, the new approach proposed in this study produces economic controllers with low energy consumption which may be of vital importance for practical application.Let system (40) be the drive system and let the following be the response system:

Figure 3 :
Figure 3: Recursive backstepping tracking control of chaotic Lorenz system for () = 35 cos 0.57 when the controller is activated at  = 50.The parameters of the system are as shown in Figure 1.

Figure 4 :
Figure4: Extended recursive backstepping tracking control of chaotic Lorenz system for () = 35 cos 0.57 when the controller is activated at  = 50.The parameters of the system are as shown in Figure1with   = 0.5, where  = 1, 2, 3.

Figure 5 :
Figure 5: Recursive backstepping tracking control of hyperchaotic Lorenz system for () = 35 cos 0.57 when the controller is activated at  = 50.The parameters of the system are as shown in Figure 2.

Figure 6 :
Figure 6: Extended recursive backstepping tracking control of hyperchaotic Lorenz system for () = 35 cos 0.57 when the controller is activated at  = 50.The parameters of the system are as shown in Figure 2 with  1 =  2 = 0.7 and  3 =  4 = 0.8.

Figure 7 :
Figure 7: Error dynamics between the two chaotic Lorenz systems with the controller deactivated for 0 <  < 50 and activated for  ≥ 50.The parameters of the system are as shown in Figure 1.

Figure 8 :
Figure 8:  Error dynamics between the two chaotic Lorenz systems with the controller deactivated for 0 <  < 50 and activated for  ≥ 50.The parameters of the system are as shown in Figure1with   = 0.3, where  = 1, 2, 3.

Figure 9 :
Figure 9: Error dynamics between the two hyperchaotic Lorenz systems with the controller deactivated for 0 <  < 50 and activated for  ≥ 50.The parameters of the system are as shown in Figure 2.

Figure 10 :
Figure 10: Error dynamics between the two hyperchaotic Lorenz systems with the controller deactivated for 0 <  < 50 and activated for  ≥ 50.The parameters of the system are as shown in Figure 2 with   = 0.3, where  = 1, 2, 3, 4.