A Research on Active Control to Synchronize a New 3d Chaotic System

This paper presents the robust synchronization problem of a 3D chaotic system by using the active control technique. Based on the Gershgorin theorem and Routh-Hurwitz criterion, sufficient algebraic conditions are derived to design a linear controller gain matrix. The conditions are then applied for the robust stability of the synchronization error dynamics in the presence of an unknown bounded smooth external disturbance. The proposed active control strategy with a suitable computation of the linear controller gain matrix is simple in design and establishes fast convergence rates of the synchronization error signals. Numerical simulation results further verified the analytical results.


Introduction
For the last two decades, chaos synchronization has been widely explored in various research fields [1][2][3].The ever present demand for secure communication is one of the foremost form of motivation for studying chaos synchronization [4][5][6].Various linear and nonlinear control techniques have been developed to carry out synchronization behavior.Some of these include, adaptive control [7], linear error state feedback control [8], backstepping method [9], active control [10], sliding mode control [11], and nonlinear control [12] are worth citing here, among others.
Recently, active control strategies have received considerable interest among many researchers due to the simple implementation for chaotic synchronization.Using the active control algorithm, Agiza and Yaseen [13] presented the synchronization of two identical Rossler and two identical Chen chaotic systems.Lei et al. [14] reported the synchronization problem of two identical nonlinear gyros using the active control based on the Routh-Hurwitz criterion and Lyapunov stability theory.Ucar et al. [15] utilized the active control technique for chaos synchronization in RCL-shunted Josephson function.Njah and Vincent [16] investigated the synchronization and anti-synchronization of two identical extended BVP chaotic oscillators by using the generalized active control method.Based on the Routh-Hurwitz criterion and Lyapunov stability theory and using the active control strategy, Shahzad and Ahmad [17] showed an experimental study of synchronization and anti-synchronization for the spin orbit problem of Enceladus.Recently, Ahmad et al. [10], addressed the synchronization problem of two identical and two non-identical chaotic systems using the active control technique based on the Lyapunov stability theory.A modification has been made based on the Lyapunov stability theory [18] for the globally, asymptotically-stable synchronization.
However, the finding of these studies [10,[13][14][15][16][17] guarantee the asymptotic stability of the resulting synchronization error dynamics by allocating negative eigenvalues of the coefficient matrix in the closed-loop system that simply obeys the Routh-Hurwitz criterion.Convergence rates of the synchronization error signals depend on the numeric values of the controller gain coefficients.In real-life applications, selecting high controller gain coefficients may lead to automatic signal saturations and the two coupled chaotic/hyperchaotic systems may lose synchronization stability completely.In fact, the above notable results [10,[13][14][15][16][17] affect each other mutually and need a systematic approach to compute suitable linear controller gains.In addition, the chaotic systems are considered free of the unknown external disturbance.
In this article, based on the Gershgorin theorem [19] and Routh-Hurwitz criterion [20] and using the active control strategy, a generalized approach is proposed to compute a suitable linear controller gain matrix that establishes globally asymptotical synchronization under the effect of unknown external disturbance.
There are three main aims for this paper.The first aim is to study the synchronization problem of a newly reported 3D chaotic system [21] in the presence of an unknown external disturbance.The second aim is to discuss the active control method for the synchronization of two identical and two non-identical 3D chaotic systems.The third aim is to derive analytically the sufficient algebraic conditions to compute a suitable linear control gain matrix that guarantees globally asymptotical stability of the closed-loop system.Two illustrative examples are given to verify the robustness and the performance of the proposed approach; complete synchronization between two identical chaotic systems [21], and the generalized synchronization between the chaotic system [21] and the Liu-Chen chaotic system [22].
The rest of the paper is organized as follows: In Section 2, the problem statement is given and derived the generic criterion to construct a proper linear control gain matrix using the active control technique.In Section 3, description of a new 3D chaotic system [21] is briefly described and synchronization behavior of two identical chaotic systems [21] are given, followed by the analysis of synchronization between the new chaotic and the Liu-Chen chaotic systems in Section 4. The paper concludes in Section 5.

Problem Statement and a Theory for the Proposed Active Control Strategy
Let us consider a drive-response chaotic system synchronization scheme is described as follows:
where X ptq " rx 1 ptq , x 2 ptq , . . ., x n ptqs T P R n and Y ptq " ry 1 ptq , y 2 ptq , . . ., y n ptqs T P R n are the state vectors, F pX ptqq and G pY ptqq are the nonlinear bounded smooth functions, and A 1 , A 2 P R nˆn are the matrices of the chaotic drive and response systems parameters alternatively, ∆ ptq and Ω ptq are the vectors of the unknown bounded smooth external disturbances present in the drive and response systems, respectively and δ ptq P R n is the control input to be designed later.

Assumption 1.
It is assumed that the unknown external disturbances are norm-bounded [2].
where a n i , b n i and Φ i are positive constants.

Definition 1.
Let us define the synchronization error system as follows: where || ¨|| represents the Euclidean norm in R n .
From the drive-response system arrangement Equation (1), it follows that: . e ptq " Y ptq A 2 ´X ptq A 1 `G pY ptqq ´F pX ptqq `Ω ptq ´∆ ptq `δ ptq ñ .e ptq " A 3 e ptq `H pX ptq , Y ptq , e ptqq `Ω ptq ´∆ ptq `δ ptq ñ .e ptq " A 3 e ptq `H pX ptq , Y ptq , e ptqq `|Ω ptq ´∆ ptq| `δ ptq (4) where H pX ptq , Y ptq , e ptqq " G pY ptqq ´F pX ptqq ``A 2 ´A2 ˘Y ptq ´`A 1 ´A1 ˘X ptq is the function of residual terms and A 2 ´A1 " A 3 is the (n ˆn) coefficient matrix of the error system Equation (4) [23]. If The chaotic synchronization problem can be considered as a stabilization of the synchronization error at the origin by a suitable controller δ ptq P R n .This means that the controller directs the synchronization error trajectories to the origin for all t ě t 0 ě 0, where t 0 is the time of control activation.

Active Controller Design
Theorem 1.Let us define the following control function: where v ptq " ´Ke ptq " ´K pY ptq ´X ptqq " ´Ke ptq is the sub-controller function with K P R nˆn is a linear controller gain matrix that regulates the strength of the feedback controller into the response system.
Proof of Theorem 1.It is assumed that the parameters of the drive and response systems are available and measurable.Using systems of Equations ( 4) and ( 6) that yields: . e ptq " A 3 e ptq `v ptq ñ where, A 3 ´K " B P R nˆn .At this stage, the problem is reduced to show that the matrix B P R nˆn of the closed-loop Equation ( 7) is Hurwitz.
Lemma 1.The necessary and sufficient condition for the matrix B P R nˆn of the closed-loop system Equation (7) being a Hurwitz matrix, if all eigenvalues of the matrix B P R nˆn have negative real parts, then, the zero solution of the closed-loop system Equation ( 7) is globally asymptotically stable.
The proof for Lemma 1 can be found in [20].Now, we state the main theorem.
Theorem 2. The coefficient matrix B P R nˆn of the closed-loop system Equation ( 7) is Hurwitz if the linear controller gain matrix K P R nˆn is design such that coefficient matrix B P R nˆn in Equation (7) satisfies the following two conditions: . b ii ă 0, for alli " 1, 2, . . ., n Then, the closed-loop system Equation ( 7) is globally asymptotically stable.
Proof of Theorem 2. Let us assumed that B P R nˆn be a real symmetric matrix.Further, it is assumed that the coefficients of the feedback controller gain matrix K P R nˆn are designed such that the system matrix B P R nˆn satisfies the two conditions of the Theorem 2.Then, the matrix B P R nˆn is a strictly diagonally dominant (SDD) non-singular matrix.By the Gershgorin theorem [19], if a matrix is SDD and all its diagonal elements are negative, then, every eigenvalue of the system matrix B P R nˆn is also negative.Thus, by the Routh-Hurwitz criterion [20], the matrix B P R nˆn is Hurwitz.Hence, by the linear control theory [20], the closed-loop system Equation ( 7) is globally asymptotically stable.This implies that the two coupled chaotic systems Equation (1) are globally asymptotical synchronized.

Description of the New 3D Chaotic System
Recently, a new 3D chaotic system is reported [21], which is contructed by replacing a constant parameter e with a switching function in Qi four-wing chaotic attractor and generated an eight-wing chaotic attractor.The system of differential equations describing the chaotic system [21] are given as follows: .
/ - (9) where rx ptq , y ptq , z ptqs T P R 3 are the state vectors and p, q, r and s are all real positive system parameters and f ptq " ξsgnpsinωtq `λ, is a parameter function with ω as the switching frequency, ξ, and λ are the constant parameters which takes values within a certain range.It has been investigated [21] that when f ptq " 1 , originally the system [21] shows a four-wing chaotic attractor and when ξ " 9 and λ = 10, the parameter function f ptq " ξsgnpsinωtq `λ with ω " 0.04π, then, the new system displays an eight-wing chaotic attractor (mother butterfly) and the new four-wing attractor (Baby butterfly) is very close to the origin [21].The proposed eight-wing chaotic attractor has a complex dynamics and topologically different structure than the original Qi four-wing chaotic attractor as shown in Figure 1.Due to the complex properties of the system Equation ( 9), it is now significant to synchronize the chaotic system [21] for further research purposes on theoretical ground in order to implement it for hybrid image encryption, secure communications, and chaotic masking, etc.

Problem Statement
In this sub-section of the paper, we briefly describe the identical synchronization scheme for the two coupled chaotic systems [21].The drive chaotic system Equation ( 9) with three state variables denoted by the subscript 1 drives the response chaotic system having identical equations denoted by the subscript 2. Although, the initial conditions of the two systems are different.The two coupled chaotic systems [21] are described in a drive-response system synchronization scheme as follows: Drive system : Response system : , , and , , T T x y z x y z R  are the state vectors and   , , , , p q r s f t are the positive parameters of the drive and response systems Equations (10) alternatively, as the control input injected to the response system, and are the vectors of unknown bounded smooth external disturbances present in the drive and response systems respectively.
Definition 2. The error dynamical system for the synchronization scheme Equation ( 10) is described as follows: where, , , e t x x e t y y e t z z       .
Theorem 3. Let us define the following active control functions: Note.For our own simplicity, we will use rx, y, zs for rx ptq , y ptq , z ptqs alternatively.

Problem Statement
In this sub-section of the paper, we briefly describe the identical synchronization scheme for the two coupled chaotic systems [21].The drive chaotic system Equation ( 9) with three state variables denoted by the subscript 1 drives the response chaotic system having identical equations denoted by the subscript 2. Although, the initial conditions of the two systems are different.The two coupled chaotic systems [21] are described in a drive-response system synchronization scheme as follows: Drive system : x 2 " p py 2 ´x2 q `f ptq y 2 z 2 `Ω1 ptq `δ1 ptq .y 2 " qx 2 `ry 2 ´x2 z 2 `Ω2 ptq `δ2 ptq .z 2 " x 2 y 2 ´sz 2 `Ω2 ptq `δ3 ptq (10) where rx 1 , y 1 , z 1 s T and rx 2 , y 2 , z 2 s T P R 3 are the state vectors and p, q, r, s, f ptq are the positive parameters of the drive and response systems Equations (10) alternatively, δ ptq " rδ 1 ptq , δ 2 ptq , δ 3 ptqs T P R 3 as the control input injected to the response system, and ∆ i ptq and Ω i ptq P R 3 are the vectors of unknown bounded smooth external disturbances present in the drive and response systems respectively.Definition 2. The error dynamical system for the synchronization scheme Equation ( 10) is described as follows:  (11) where, e 1 ptq " x 2 ´x1 , e 2 ptq " y 2 ´y1 , e 3 ptq " z 2 ´z1 .Theorem 3. Let us define the following active control functions: Then, the two coupled identical chaotic systems Equation (10) are globally asymptotically synchronized.
Proof of Theorem 3. It is assumed that all the variables and parameters of the two chaotic systems Equation (10) are available and measureable.Using systems of Equations ( 11) and (12), that yields: e ptq " Be ptq where The closed-loop system Equation ( 15) to be controlled is a linear system with control input v 1 ptq , v 2 ptq and v 3 ptq as a function of e 1 ptq , e 2 ptq and e 3 ptq alternatively, where the constants k ij , i, j " 1, 2, 3 are the linear controller gains.As long as these controller gains stabilize the closed-loop system Equation (15), the e 1 ptq , e 2 ptq and e 3 ptq converge to zero as time goes to infinity [16], then, the zero solution of the closed-loop system Equation ( 15) is obvious.This implies that the two coupled identical chaotic systems Equation (10) are globally asymptotically synchronized.
At this stage, the problem is reduced to show that the matrix B P R 3ˆ3 of the closed-loop system Equation ( 15) is Hurwitz.
Corollary 2. If the coefficient matrix B P R 3ˆ3 in Equation ( 17) of the closed-loop system Equation ( 15) is Hurwitz, then the zero solution of the error system Equation ( 15) is globally asymptotically stable.Theorem 4. The matrix B P R 3ˆ3 is Hurwitz, if the linear controller gain matrix K P R 3ˆ3 is constructed such that the coefficient matrix B P R 3ˆ3 in Equation ( 17) satisfies the following two conditions: Proof of Theorem 4. The controller Equation ( 12) is associated to the feedback controller gain matrix K P R 3ˆ3 so that the design of K P R 3ˆ3 is to make all eigenvalues of the coefficient matrix B P R nˆn of the closed-loop system Equation (15) are in the open left-half of the complex plane.For the particular choice of the gain coefficients, which satisfies the conditions Equation ( 18) are selected as follows: With this specific choice of the controller gains and considering: a " 0.4, b " 5, c " 12, p " 14, q " ´1, r " 16 and s " 43, the closed-loop system Equation ( 15) yields: Thus, by the Routh-Hurwitz criterion, the matrix B P R nˆn is Hurwitz.This proves that the closed-loop system Equation ( 15) is globally asymptotically stable.Consequently Theorem 4 and Corollary 2 are obvious.Hence Theorem 3 is proved.

Numerical Simulation and Discussion
Numerical simulation results are furnished to validate the advantages and effectiveness of the proposed approach by using Mathematica 10.The parameters for the chaotic system [21] are set as p " 14, q " ´1, r " 16, s " 43, f ptq " 7 with initial values of the state variables of the drive and response systems being chosen as rx 1 p0q , y 1 p0q , z 1 p0qs T " r10, ´10, ´15s and rx 2 p0q , y 2 p0q , z 2 p0qs T " r´25, ´50, ´40s respectively.In the simulation, the following external disturbances are applied to the drive and response systems respectively.Subsequently, Φ 1 " 0.4, Φ 2 " 0.5 and Φ 3 " 0.6, that satisfies the assumption 1.
For the identical chaotic systems [21], the time series of the state vectors of the synchronized and unsynchronized trajectories are depicted in Figures 2-4.The state trajectories of the drive system converged to the state trajectories of the response system under the control action Equation (12).
The time series of the error system Equation ( 11) for identical synchronization are shown in Figure 5.As expected, one can notice the smoothness of the synchronized error signals while converging to the zero state quickly, when the controllers are activated at t ă 2s.This demonstrates the robustness and performance of the control action Equation ( 12) for the complete synchronization under the effect of unknown external disturbance.

Problem Formulation
To study the generalized synchronization problem for the chaotic system [21], it is assumed that the chaotic system [21] drives the Liu-Chen chaotic system [22].Thus, the drive-response system synchronization scheme is described as follows: Drive system : Response system : x p y x f t y z t y qx ry x z t z x y sz t x ax y z t t y cy x z t t z x y bz t t where , , and , , are the corresponding state vectors, , , and pqr s and , and ab c are the parameters of the drive and response systems respectively,     3 and are the vectors of unknown bounded smooth external disturbances present in the drive and response systems Equation ( 22) respectively, and         as the active controller that is to be designed yet.
Definition 3. The error dynamical system for the synchronization scheme Equation ( 22) is be described as follows: Theorem 5. Let us define the following active control functions:

Problem Formulation
To study the generalized synchronization problem for the chaotic system [21], it is assumed that the chaotic system [21] drives the Liu-Chen chaotic system [22].Thus, the drive-response system synchronization scheme is described as follows: Drive system : x 2 " ax 2 ´y2 z 2 `Ω1 ptq `δ1 ptq .y 2 " ´cy 2 `x2 z 2 `Ω2 ptq `δ2 ptq .z 2 " x 2 y 2 ´bz 2 `Ω3 ptq `δ3 ptq (22) where rx 1 , y 1 , z 1 s T P R 3 and rx 2 , y 2 , z 2 s T P R 3 are the corresponding state vectors, p, q, r and s and a, b and c are the parameters of the drive and response systems respectively, ∆ i ptq and Ω i ptq P R 3 are the vectors of unknown bounded smooth external disturbances present in the drive and response systems Equation ( 22) respectively, and δ ptq " rδ 1 ptq , δ 2 ptq , δ 3 ptqs T P R 3 as the active controller that is to be designed yet.Definition 3. The error dynamical system for the synchronization scheme Equation ( 22) is be described as follows:
Proof of Theorem 5.It is assumed that all the variables and parameters of the two systems Equation ( 22) are available and measureable.Using systems of Equations ( 23) and (24), that yields: Proof of Theorem 6.For the particular choice of the linear controller gains which satisfies the conditions Equation (28) are selected as follows: With this specific selection of the linear controller gains and considering: a " 0.4, b " 5, c " 12, p " 14, q " ´1, r " 16 and s " 43, the closed-loop system Equation ( 26 This completes the proof.
For the non-identical chaotic systems [21,22], the time series of the state vectors of the synchronized and unsynchronized trajectories are illustrated in Figures 6-8.These figures demonstrated that the state trajectories of the response system converged to that state trajectories of the drive system under the synthesized control action Equation (24), while the uncontrolled state trajectories of the response system are completely different from the drive system state trajectories.Time series of the error system Equation ( 23) to the zero state are depicted in Figure 9.One can notice the smoothness of the synchronized error trajectories while converging to the zero state with fast converging rates, which illustrates the performance of the control action Equation (24) for non-identical synchronization.This feature shows that the investigated controllers are robust to the accidental parameters mismatch and unknown external disturbances in the transmitter and receiver.

 
e t e t e t This completes the proof.

Simulation Studies
The parameters for the chaotic system [21] and Liu-Chen system [22] are set as with initial values of the drive [21] and response [22] systems being chosen as alternatively.In the simulation, the following bounded smooth disturbances are applied to the drive and response systems, respectively: Accordingly, 1 2 3 0.4, 0.6 and 0.6 For the non-identical chaotic systems [21,22], the time series of the state vectors of the synchronized and unsynchronized trajectories are illustrated in Figures 6-8.These figures demonstrated that the state trajectories of the response system converged to that state trajectories of the drive system under the synthesized control action Equation (24), while the uncontrolled state trajectories of the response system are completely different from the drive system state trajectories.Time series of the error system Equation ( 23) to the zero state are depicted in Figure 9.One can notice the smoothness of the synchronized error trajectories while converging to the zero state with fast converging rates, which illustrates the performance of the control action Equation (24) for non-identical synchronization.This feature shows that the investigated controllers are robust to the accidental parameters mismatch and unknown external disturbances in the transmitter and receiver.

Conclusions
The robust synchronization problem of two identical and two non-identical chaotic systems has been investigated.Based on the Gershgorin theorem and Routh-Hurwitz criterion and using the active control strategy, sufficient conditions are derived to construct a suitable linear gain matrix that established the globally asymptotical synchronization under the effect of an unknown external

Conclusions
The robust synchronization problem of two identical and two non-identical chaotic systems has been investigated.Based on the Gershgorin theorem and Routh-Hurwitz criterion and using the active control strategy, sufficient conditions are derived to construct a suitable linear gain matrix that established the globally asymptotical synchronization under the effect of an unknown external

Conclusions
The robust synchronization problem of two identical and two non-identical chaotic systems has been investigated.Based on the Gershgorin theorem and Routh-Hurwitz criterion and using the active control strategy, sufficient conditions are derived to construct a suitable linear gain matrix that established the globally asymptotical synchronization under the effect of an unknown external

Conclusions
The robust synchronization problem of two identical and two non-identical chaotic systems has been investigated.Based on the Gershgorin theorem and Routh-Hurwitz criterion and using the active control strategy, sufficient conditions are derived to construct a suitable linear gain matrix that established the globally asymptotical synchronization under the effect of an unknown external disturbance.Numerical simulation results show that the proposed approach works very well and can be applied to a class of chaotic as well as hyperchaotic systems successfully.
In comparison with some published related works, there are three main advantages of the proposed approach which are summarized as follows.
(a) The synchronization speed is fast as well as the amplitude of the oscillations is smaller.(b) In the proposed active synchronization control approach, the eigenvalues of the coefficient matrix of the closed-loop system can be adjusted to have a desirable synchronization time.(c) Most of the chaotic systems in real practical applications have different structures, thus, we believe that the proposed active control approach will be a helpful tool in synchronizing a class of chaotic/hyperchaotic systems.

Figure 1 .
Figure 1.(a) 2D phase portrait of the new 3D chaotic system in XY-space; (b) XZ-space; and (c) YZ-space.Note.For our own simplicity, we will use   , , x y z for

Figure 1 .
Figure 1.(a) 2D phase portrait of the new 3D chaotic system in XY-space; (b) XZ-space; and (c) YZ-space.Note.For our own simplicity, we will use rx, y, zs for rx ptq , y ptq , z ptqs alternatively.

Figure 2 .
Figure 2. Time series of the state trajectories 1

Figure 3 .
Figure 3.Time series of the state trajectories 1

Figure 4 .
Figure 4. Time series of the state trajectories 12 and zz .

Figure 2 . 8 Figure 2 .
Figure 2. Time series of the state trajectories x 1 and x 2 .

Figure 3 .
Figure 3.Time series of the state trajectories 1

Figure 4 .
Figure 4. Time series of the state trajectories 12 and zz .

Figure 3 . 8 Figure 2 .
Figure 3.Time series of the state trajectories y 1 and y 2 .

Figure 3 .
Figure 3.Time series of the state trajectories 1

Figure 4 .
Figure 4. Time series of the state trajectories 12 and zz .

Figure 4 .
Figure 4. Time series of the state trajectories z 1 and z 2 .

Figure 5 .
Figure 5.Time series of the controlled synchronization error system Equation (11).

Figure 5 .
Figure 5.Time series of the controlled synchronization error system Equation (11).

Figure 6 .
Figure 6.Time series of the state trajectories 12 and xx .

Figure 6 .
Figure 6.Time series of the state trajectories x 1 and x 2 .

Figure 7 .
Figure 7. Time series of the state trajectories 12 and yy .

Figure 8 .
Figure 8.Time series of the state trajectories 12 and zz .

Figure 9 .
Figure 9.Time series of the controlled synchronization error system Equation (26).

Figure 7 . 12 Figure 7 .
Figure 7. Time series of the state trajectories y 1 and y 2 .

Figure 8 .
Figure 8.Time series of the state trajectories 12 and zz .

Figure 9 .
Figure 9.Time series of the controlled synchronization error system Equation (26).

Figure 8 . 12 Figure 7 .
Figure 8.Time series of the state trajectories z 1 and z 2 .

Figure 8 .
Figure 8.Time series of the state trajectories 12 and zz .

Figure 9 .
Figure 9.Time series of the controlled synchronization error system Equation (26).

Figure 9 .
Figure 9.Time series of the controlled synchronization error system Equation (26).
then, X ptq and Y ptq are the states of two identical chaotic systems and if, A 1 ‰ A 2 , then, X ptq and Y ptq are the states of two different chaotic systems.