ABSTRACT

A three-dimensional system is presented with unknown parameters that employs two nonlinearities terms. The basic characteristics of the system are studied. The stability is measured by Characteristic equation roots, Routh stability criteria, Hurwitz stability criteria and Lapiynov function, all show that the system unstable. Then, Chaoticity is measured by maximum Lapiynov exponent of ( $L_{\max }=2.509426$ ) and “Kaplan-Yorke” dimension ( $D_L=2.22349544$ ). The system is controlled effectively and synchronized by designed adaptive controllers. Furthermore, the theoretical and graphic results of the system before and after control are compared.

Subject Areas:

Dynamical System

Keywords:

Stabilization, Dissipative System, Adaptive Control, Lapiynov Exponent

1. Introduction

Researches in recent years on chaotic phenomena have increased a lot, because of the increasing frontiers of applications of chaos in engineering and non-engineering systems. “Chaos is a phenomenon which results from the exhibits sensitivity to perturbation in the structural parameters and initial conditions of some classes of dynamical systems” [1]. “Chaotic signals have a random-like nature and broadband spectrum and are non-periodic” [2]. “For a system to be chaotic, the following conditions must be satisfied. Firstly, it must be sensitive to perturbation in its initial conditions which should lead to unpredictability of its future trajectories, secondly, it is not topologically transitive and thirdly, the chaotic orbits are dense in the phase space” [3]. “Among some evolved chaotic attractors in the literature are the Chen’s” [4], “3-D, 4-wing attractor” [5], “Sundarapandian-Pehlivan” [6], “Morphous one parameter attractor” [7], “Rabinovich system” [8]. “When chaotic attractors possess one positive Lapiynov exponent, then the system is chaotic. However, the system which has two or more positive Lapiynov exponents it is a highly chaotic system and becomes hypersensitive to small perturbations in its system dynamics” [9] [10]. “Chaos Control subject has received widespread attention of research because controllability and synchronizability of chaotic attractors are index of utility in different designs such as in secure communications and robotics” [11] [12] [13].

“In the context of stability and stabilization, the principle of Lapiynov stability continued to enjoy large applications; it can effectively stabilize the dissipative systems” [14] [15].

This paper is organized as: Section 2, present a description of 3-d system. Section 3, basic analysis such as stability, dissipativity, Lapiynov dimension ‘‘Kaplan-Yοrke dimension’’. Section 4, we designed adaptive control law of the chaotic system. Section 5, a comparison of the analysis results before and after control. Section 6, we derive results for the adaptive synchronization of identical highly chaotic system. Finally Section 7, summarization of the main results.

2. System Description

A three-dimensional dynamical system [16] consist of three ordinary differential equations with state variables $x_i$ , (i = 1, 2, 3) and four unknown parameters (ρ, α, δ and φ), employs six terms include two quadratic cross-product nonlinear terms. Given by:

$\begin{array}{l}{\stackrel{\dot{}}{x}}_1=\rho \left(x_2-x_1\right)\\ {\stackrel{\dot{}}{x}}_2=a{x}_1-\delta x_1{x}_3\\ {\stackrel{\dot{}}{x}}_3=\phi x_1{x}_2-x_3\end{array}$ (1)

The parameters values are taken as

$\rho =10,\delta =40,a=296.5,\phi =10$ (2)

3. System Analysis

In this section essential, the system (1) is invested and has the following characteristics.

3.1. Equilibrium Points

The first step to analyze a system is to find its equilibrium points, so we need to solve the nonlinear equations as follows

$-10x_1+10x_2=0$

$296.5x_1-40x_1{x}_3=0$

$10x_1{x}_2-x_3=0$

We get the following equilibrium points:

$E_0=\left(0,0,0\right)$ , $E_1=\left(\frac{\sqrt{\frac{593}{2}}}{20},\frac{\sqrt{\frac{593}{2}}}{20},\frac{593}{80}\right)$ , $E_2=\left(-\frac{\sqrt{\frac{593}{2}}}{20},-\frac{\sqrt{\frac{593}{2}}}{20},\frac{593}{80}\right)$

3.2. Stability Analysis

3.2.1. Characteristic Equation Roots

“A necessary and sufficient condition for the system to be stable is that the real parts of the characteristic equation have negative real parts”.

When the parameters values are taken as in (2), the Jacobian matrix of system (1) at $E_0=\left(0,0,0\right)$ is:

$J=\left[\begin{array}{ccc}-10& 10& 0\\ 296.5& 0& 0\\ 0& 0& -1\end{array}\right]$

$\det \left(J-\lambda I\right)=0$

$\Rightarrow {\lambda }^3+11{\lambda }^2-2955\lambda -2965=0$ (3)

By using Horner’s Ruffini method [17] we get:

${\lambda }_1=-1$ , ${\lambda }_2=49.68089$ , ${\lambda }_3=-59.68089$

Similarly, we find Jacobian matrix at $E_1$ and $E_2$ , then we obtain the eigenvalues, as shown in Table 1.

So the system (1) is unstable.

3.2.2. Routh Stability Criteria

“The system is considered stable by the Routh stability states (all poles in OLHP (Open Loop Half plane)) if and only if all elements of the first column in the Routh array are positive. In addition, number of poles not in the OLHP is equal to the number of sign changes in the first column” [18].

$a_0=-2965$

$a_1=-2955$

$a_2=11$

$a_3=1$

$b_1=a_1-\frac{a_3{a}_0}{a_2}=-2685.4545$

Since, there are two negative elements in the first column. Therefore, the system (1) is unstable.

3.2.3. Hurwitz Stability Criteria

“This criterion is applied using determinants formed from coefficients of the characteristic equation. If the small minors of the square matrix J of the system (1) are all positive then the system (1) is stable, otherwise it’s unstable” [18].

If n = 3 (n denote the degree of the square matrix)

From Equation (3) we find:

${\Delta }_1=a_2=11>0$

${\Delta }_2=\left|\begin{array}{cc}{a}_2& a_0\\ a_3& a_1\end{array}\right|=a_2{a}_1-a_3{a}_0=-29.540<0$

${\Delta }_3=\left|\begin{array}{ccc}{a}_2& a_0& 0\\ a_3& a_1& 0\\ 0& a_2& a_0\end{array}\right|=a_2{a}_1{a}_0-a_0^2{a}_3=87586100>0$

Since the values of one minors is less than zero, so the system (1) is unstable.

3.2.4. Lapiynov Function

We can use quadratic function for system (1).

We assume that

$V\left(x_1,x_2,x_3\right)=\frac{1}{2}\left(x_1^2+x_2^2+x_3^2\right)$

$\stackrel{\dot{}}{V}\left(x_1,x_2,x_3\right)=x_1{\stackrel{\dot{}}{x}}_1+x_2{\stackrel{\dot{}}{x}}_2+x_3{\stackrel{\dot{}}{x}}_3$ (4)

If $\stackrel{\dot{}}{V}\left(x_1,x_2,x_3\right)<0$ then the system is stable.

By substituting (1) in Equation (4) we get:

$\stackrel{\dot{}}{V}\left(x_1,x_2,x_3\right)=-10x_1^2+306.5x_1{x}_2-30x_1{x}_2{x}_3-x_3^2$

Since $\stackrel{\dot{}}{V}\left(x_1,x_2,x_3\right)>0$ therefore the system (1) is unstable.

3.3. Dissipativity

Let $f_1=\frac{\text{d}{x}_1}{\text{d}t}$ , $f_2=\frac{\text{d}{x}_2}{\text{d}t}$ and $f_3=\frac{\text{d}{x}_3}{\text{d}t}$ in the system (1).

Then we get for the vector field

${\left({\stackrel{\dot{}}{x}}_1,{\stackrel{\dot{}}{x}}_2,{\stackrel{\dot{}}{x}}_3\right)}^{\text{T}}={\left(f_1,f_2,f_3\right)}^{\text{T}}$

thus the divergence of the vector field V on $R^3$ yields to:

$\nabla \cdot {\left({\stackrel{\dot{}}{x}}_1,{\stackrel{\dot{}}{x}}_2,{\stackrel{\dot{}}{x}}_3\right)}^{\text{T}}=\frac{\partial f_1}{\partial x_1}+\frac{\partial f_2}{\partial x_2}+\frac{\partial f_3}{\partial x_3}=-\left(\rho +1\right)=f$

Note that $f=-\left(\rho +1\right)=-11$ , so the systеm (1) is dissipative for all positive values of $\rho$ , and an exponential rate is:

$\frac{\text{d}V}{\text{d}t}=fV\Rightarrow V\left(t\right)=V_0{\text{e}}^{ft}=V_0{\text{e}}^{-11t}$

From above equation, the volume element $V_0$ is contracted by the flow into a volume element $V_0{\text{e}}^{-11t}$ at the time t.

3.4. Numerical and Graphical Results

For the numerical solution, we use Runge-Kutta method of order 5^th to solve system (1). With initial states ${x|}_{x_1\left(0\right),x_2\left(0\right),x_3\left(0\right)}=\left[-2,7,12\right]$

3.4.1. Wave Form of the System (1)

The wave-form $x_1\left(t\right),x_2\left(t\right),x_3\left(t\right)$ for the system (1) is characteristic with non-periodic shape, shown in Figures 1(a)-(c) which is one of the basic characteristic behaviors of chaotic dynamical system.

3.4.2. Phase Portrait of the System (1)

Figures 2(a)-(d) and Figures 3(a)-(c) are shows chaotic attractor for system (1) in $\left(x_1,x_2,x_3\right)$ , $\left(x_1,x_3,x_2\right)$ , $\left(x_2,x_1,x_3\right)$ , $\left(x_3,x_1,x_2\right)$ space, and 2-D attractor of system (1) in $\left(x_1,x_3\right)$ , $\left(x_1,x_2\right)$ , $\left(x_2,x_3\right)$ plane.

The orbit is dense in each graph which means the system exhibit two-scroll hyper chaotic attractor.

3.5. Lapiynov Exponent and Lapiynov Dimension

“As a rule the Lapiynov exponents refer to the average exponential rates of divergence or convergence of nearby trajectories in the phase space. The system is chaotic if there is at least one Lapiynov exponent greater than zero”. The values of lapiynov exponents are: ( $L_1=2.509426$ , $L_2=0.132019$ and $L_3=-11.818787$ ). Therefore, the Lapiynov dimension ‘‘Kaplan-Yοrke dimension’’ is:

$D_L=2+\frac{L_1+L_2}{\left|L_3\right|}=2.22349544$

So the system (1) is Highly Chaotic System, as shown in Figure 4.

4. Adaptive Control Strategy

4.1. Theoretical Results

To stabilize highly chaotic system (1), an adaptive control law is designed with unknown parameter α.

As follows:

$\begin{array}{l}{\stackrel{\dot{}}{x}}_1=10\left(x_2-x_1\right)+{\alpha }_1\\ {\stackrel{\dot{}}{x}}_2=a{x}_1-40x_1{x}_3+{\alpha }_2\\ {\stackrel{\dot{}}{x}}_3=10x_1{x}_2-x_3+{\alpha }_3\end{array}$ (5)

when ${\alpha }_1,{\alpha }_2,{\alpha }_3$ are the feedback controllers.

The adaptive control functions are:

$\begin{array}{l}{\alpha }_1=-10\left(x_2-x_1\right)-{\mu }_1{x}_1\\ {\alpha }_2=-\hat{a}{x}_1+40x_1{x}_3-{\mu }_2{x}_2\\ {\alpha }_3=-10x_1{x}_2+x_3-{\mu }_3{x}_3\end{array}$ (6)

where the constants ${\mu }_i$ , (i = 1, 2,3) are positive , $\hat{a}$ is the parameter estimate of α.

Substituting (6) into (5), we get

$\begin{array}{l}{\stackrel{\dot{}}{x}}_1=-{\mu }_1{x}_1\\ {\stackrel{\dot{}}{x}}_2=\left(a-\hat{a}\right)x_1-{\mu }_2{x}_2\\ {\stackrel{\dot{}}{x}}_3=-{\mu }_3{x}_3\end{array}$ (7)

Let the parameter estimation error

$e_a=a-\hat{a}$ (8)

Using (8), the dynamics (7) can be written compactly as

$\begin{array}{l}{\stackrel{\dot{}}{x}}_1=-{\mu }_1{x}_1\\ {\stackrel{\dot{}}{x}}_2=e_a{x}_1-{\mu }_2{x}_2\\ {\stackrel{\dot{}}{x}}_3=-{\mu }_3{x}_3\end{array}$ (9)

The Lapiynov approach is used for derivation of update law for adjusting the parameter estimate $\hat{a}$ .

Consider the lapiynov function

$V\left(x_1,x_2,x_3\right)=\frac{1}{2}\left(x_1^2+x_2^2+x_3^2+e_a^2\right)$ (10)

Notice V is positive-definite on $R^4$ .

Also

${\stackrel{\dot{}}{e}}_a=-\stackrel{\dot{}}{\hat{a}}$ (11)

Differentiating V with substituting (9) and (11), we get:

$\stackrel{\dot{}}{V}=-{\mu }_1{x}_1^2-{\mu }_2{x}_2^2-{\mu }_3{x}_3^2+e_a\left[x_1{x}_2-\stackrel{\dot{}}{\hat{a}}\right]$ (12)

In Equation (12), we update estimated parameter by:

$\stackrel{\dot{}}{\hat{a}}=x_1{x}_2+{\mu }_4{e}_a$ (13)

where the constant ${\mu }_4$ is a positive.

Now, we substitute (13) into (12), we obtain

$\stackrel{\dot{}}{V}=-{\mu }_1{x}_1^2-{\mu }_2{x}_2^2-{\mu }_3{x}_3^2-{\mu }_4{e}_a^2$ (14)

Notice $\stackrel{\dot{}}{V}$ is negative definite on $R^4$ .

Thus, by lapiynov stability, Routh-array criteria, Eigenvalues and Hurwitz stability criteria we get the below result.

Proposition 1. The chaotic system (5) with unknown parameter is stabilized for every initial value by adaptive control (6), where the estimated parameter is obtained by (13) and ${\mu }_1,{\mu }_2,{\mu }_3,{\mu }_4$ are greater than zero.

4.2. Numerical Results

To simulate the controlled highly chaotic system (7) we take the initial values ${x|}_{x_1\left(0\right),x_2\left(0\right),x_3\left(0\right)}=\left[12,9,17\right]$ and $\left[{\mu }_1,{\mu }_2,{\mu }_3\right]=\left[30,50,30\right]$ .

Figure 5 shows Controlled trajectories of system (1).

5. System Comparison Tables before & after Control

See Tables 1-6.

6. Adaptive Synchronization Technique

6.1. Theoretical Results

We apply adaptive synchronization technique of highly chaotic system with unknown parameter α.

The drive system is

(15)

where, (i = 1, 2, 3) are state variables.

As the response system, the controlled highly chaotic dynamics given by

(16)

where are nonlinear controllers to be designed, and the state variables are, (i = 1, 2, 3).

The synchronization error is defined by

(17)

then the error dynamics is obtained as

(18)

The adaptive control functіοns define as

(19)

where the constants greater than zero, and is the estimated value of the parameter α.

Substitute (19) into (18), to obtain the error dynamics as

(20)

Now, the parameter estimation error is

(21)

By substituting (21) into (20), the error dynamics simplifies to

(22)

From Lapiynov approach we derive the updated law to adjust the estimation of the parameter.

The quadratic lapiynov function is

(23)

which be a positive definite on.

Note that

(24)

Differentiating V and substituting (22) & (24) in it, we get:

(25)

update the estimated parameter in Equation (25) by the following

(26)

where the constant is greater than zero.

From (25) and (26), we obtain:

(27)

We note that (27) is negative definite on.

Hence, by Lapiynov stability [14], “it is immediate that the parameter error and synchronization error decay exponentially to zero with time for all initial values”.

Thus, we proved the results below.

Proposition 2. The drive and response identical chaotic systems (15) and (16) with unknown parameter α are synchronized for all initial values by adaptive control law (19), where the estimated parameter given by (26) and are constants greater than zero.

6.2. Numerical Results

To get the results numerically, we used the 4th-order Runge-Kutta method to solve systems (15) & (16), and solve system (18) with adaptive control law (19).

We take and as initial states of the drive system (15) and the response system (16) respectively. Also take α = 296.5 and for i = 1, 2, 3, 4.

Figure 6 shows adaptive synchronization of the highly chaotic system.

Figure 7 shows the convergent for system (18) with controller (19).

7. Conclusion

A three-dimensional dynamical system is dealt in this paper, it has quadratic cross-product nonlinear terms. The basic characteristics are analyzed by equilibrium points, stability analysis (such as characteristic equation roots, Routh criterion, Hurwitz criterion and Lapiynov function) all methods of stability shows that the system is unstable. Then, dissipativity analysis indicates system (1) is dissipative for positive values of the parameter. Lapiynov exponent, lapiynov dimension and wave-form analysis present the hyper chaos behavior when the

parameters taken as, and the maximum values of lyapenov exponents are:, and, lapiynov dimension “Kaplan-Yοrkе dimension” of the system is, which means that the system is highly chaotic system. Moreover, to stabilize the highly chaotic system, we produced an adaptive control strategy. Finally, we proposed adaptive synchronization scheme for identical highly chaotic system with upadate law for the estimiation of system parameter. Synchronization schemes are established by Lapiynov stability. Furthermore, we compared theoretical and graphical results of the system before and after control.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Aziz, M.M. and Merie, D.M. (2020) Stability and Adaptive Control with Sychronization of 3-D Dynamical System. Open Access Library Journal, 7: e6075. https://doi.org/10.4236/oalib.1106075

References