A Comparative Study of Synchronization Methods of Rucklidge Chaotic Systems with Design of Active Control and Backstepping Methods ()

Absana Tarammim^{*}, Musammet Tahmina Akter^{}

Department of Mathematics, Chittagong University of Engineering & Technology (CUET), Chattogram, Bangladesh.

**DOI: **10.4236/ijmnta.2022.112003
PDF
HTML XML
115
Downloads
830
Views
Citations

Department of Mathematics, Chittagong University of Engineering & Technology (CUET), Chattogram, Bangladesh.

The performance of two widely used chaos synchronization approaches, active control and backstepping control, is investigated in this study. These two methods are projected to synchronize two chaotic systems (Master/Drive of Rucklidge Systems) that are identical but have different initial conditions. The paper’s significant feature is that based on error dynamics, controllers are designed using the appropriate variable and the time synchronization between master Rucklidge and drive Rucklidge systems using both methods. The control function of the active control method is designed on the proper selection of matrices. The chaotic behavior is controlled using a recursive backstepping design based on the Lyapunov stability theory with a validated Lyapunov function. The effectiveness of the controller in eradicating the chaotic behavior from the state trajectories is also revealed using numerical simulations with Matlab. The backstepping method is superior to the active control method for synchronization of the measured pair of systems, as it takes less time to synchronize while exhausting the first one than the second one with great performance, according to numerical simulation and graphical outcomes.

Share and Cite:

Tarammim, A. and Akter, M. (2022) A Comparative Study of Synchronization Methods of Rucklidge Chaotic Systems with Design of Active Control and Backstepping Methods. *International Journal of Modern Nonlinear Theory and Application*, **11**, 31-51. doi: 10.4236/ijmnta.2022.112003.

1. Introduction

Natural sciences are founded on the basis of chaos [1], so it is applied to various disciplines in mathematics, computer science, microbiology, meteorology, biology, engineering, geology, finance, economics, algorithmic trading, politics, population dynamics, psychology, philosophy, and robotics [2] [3] [4]. A variety of engineering and natural sciences problems are modeled using chaotic dynamical systems. Dynamic systems designated by non-linear differential equations can be sensitive to initial conditions [5]. This phenomenon is recognized as deterministic chaos to mean that, even if the system’s mathematical description is deterministic, its activities demonstrate to be unpredictable. System instabilities and dynamic characteristics are generous in practice. Also, chaos is unpredictable and may lead to tremblings and exhaustion failures in mechanical systems; yet, chaos suppression is usually advantageous. Many researchers are interested in investigating mathematical models and the possibilities of chaos and its control mechanisms after discovering chaotic dynamics in the deterministic non-linear system [6] [7] [8]. The dynamic analysis of biologic and technological models has emerged as a significant research topic [9] [10] [11] [12] [13]. Many features of chaotic systems are studied, such as chaos control, chaos stability, amplitude death, chaos synchronization, pattern formation, etc. Because of its various applications in physics [14], control theory, biological networks, secure communication [15], artificial neural networks, chemical reactors [16], etc., chaos synchronization is a vital aspect in non-linear dynamical science.

Pecora and Carrol demonstrated that two chaotic systems could synchronize in 1990s [17] [18]. The concept of synchronizing two similar nonlinear chaotic systems that begin with different initial conditions and a first system (Master system) can follow the trajectories of a second one (Drive system) when a suitable control law is applied, and it appears that two chaotic systems cannot synchronize with one another. However, if the two systems share information on a regular basis, they will be able to synchronize. A wide variety of methods are projected to achieve chaos synchronization such as the linear state error feedback method [19], time-delay feedback method [20], active control approach [21] [22], impulsive method [23], backstepping approach [24] [25], and some other controlling methods are important in recent times [26] [27] [28] [29]. The above methods are applied to many practical systems such as Van der Pol Duffing oscillators [21], the Rikitake two-disc dynamo—a geophysical system [19], Chua’s circuits [30], nonlinear Bloch equations modeling nuclear magnetic resonance [20], electric circuits modeling “jerk” equation [23], complex dynamos [31], nonlinear equations of acoustic gravity waves [32] and some other chaotic system are important in recent times [24] [29] [32] [33]. Backstepping design and active control, in particular, are regarded as two powerful strategies for controlling and synchronizing chaos.

Bai and Lonngren [34] projected the uses of chaos synchronization in the active control method. In this method, the synchronization speed is fast, and the amplitude of the oscillations is less. Because of the ease of implementation of chaotic synchronization, active control schemes have piqued the interest of researchers. Lately, we have used these methods to achieve synchronization and analyzed the effect of eigenvalues [35]. These eigenvalues of the chaotic system’s coefficient matrix can be changed to obtain the required synchronization time. Using these studies, we can easily synchronize two identical chaos systems (*i.e.*, two of the same parameters) with active control. The active control process is a powerful algorithm for synchronizing two chaotic systems, both are identical or non-identical [36]. Most practical systems have non-identical components, as is widely known, for which the active control method is the more efficient procedure. If the chaotic system’s non-linearity is known, then linear active control approaches can be designed to achieve global chaos control and synchronization based on the chaotic system’s provided conditions. This approach has been used on a variety of real-world systems [34]; in the chaotic synchronization, the receptor tries to identify the chaotic signal sent from the emitter that means if two chaotic signals are asymptotically identical when the time goes to infinity so that two chaotic signals will be synchronized. There are no derivatives in the controller, and the Lyapunov exponents are not required for their implementation, and these properties provide active control methods an edge over other conventional control systems [27].

The backstepping control algorithm is a type of non-linear controller design. It is very existent for handling mismatched perturbation [31]. It operates on multiple chaotic systems regardless of whether they understand external stimulation or not, it requires just one controller to recognize synchronization between chaotic systems, and the controller has no derivatives [22]. Furthermore, the Backstepping design is a type of synthetic procedure to the controller that recursively connects the choice of a Lyapunov function [37]. The backstepping method has the advantage of being able to avoid cancellations of beneficial non-linearities acting in the system. So the backstepping method acts on stabilization and tracking other than the linearization method. It is reported [38] that, backstepping technique is employed for controlling, tracking, synchronizing many chaotic systems [39] and it can be guaranteed global stability and transient performance non-linear systems. The initiation of a non-linear function passed through by course error signal under the backstepping design framework is projected to simplify the non-linear controller design process, decrease the number of undetermined parameters, increase the strength, and reduce the energy consumption of the course-keeping controller and improved summarizing design method, while the design process of the controller is simplified to only one step. The Lyapunov function is to choose for stabilizing the system of different time steps with characterizations of the current approach and the control function is designed at the final step. Finally, this paper intends, schemes, and analyses a non-linear controller that comprehends the synchronized error fast and oscillation free convergence to zero.

To control chaos in the Rucklidge system, the Active control and Backstepping approaches are applied in this work. To provide global synchronization between two identical chaotic systems, we project active controllers and a recursive backstepping control and compare simulation results of the two strategies. In Section 1, we give the theoretical foundation of our paper. The remainder of the paper will be written as follows: Dynamical Analysis of the Rucklidge System with two Parameters, Section 2; Synchronization via Active Control is presented in Section 3, with numerical results, whereas Backstepping control design for chaos synchronization is presented in Section 4, with numerical results and a comparative research description in Section 5. Finally, Section 6 brings the paper to a close.

2. Dynamical Analysis of the Rucklidge Chaotic System

In this work, the Rucklidge chaotic system (1) is described by [40]:

$\begin{array}{l}\stackrel{\dot{}}{x}=-bx+ay-yz\\ \stackrel{\dot{}}{y}=x\\ \stackrel{\dot{}}{z}=-z+{y}^{2}\end{array}$ (1)

where $x\mathrm{,}y\mathrm{,}z$ and $a\mathrm{,}b$ are the state variables and positive constant parameters, respectively. We show in this paper that the system (1) is chaotic when the parameters are $a=6.7,b=2$. We use the chaotic system’s (8) initial values for numerical simulations as: ${x}_{1}\left(0\right)=0.001$, ${y}_{1}\left(0\right)=0.001$, ${z}_{1}\left(0\right)=0.001$.

And the chaotic system (9) as: ${x}_{2}\left(0\right)=5$, ${y}_{2}\left(0\right)=8$, ${z}_{2}\left(0\right)=4$. Figures 1-3 demonstrate the Rucklidge System’s (1) 2-D projections on $\left(x,y\right)$, $\left(y,z\right)$ and $\left(z,x\right)$ space projections, respectively.

Figure 1. On the *X*-*Y* plane, a 2-D phase portrait of the Rucklidge chaotic system.

Figure 2. On the *Y*-*Z* plane, a 2-D phase portrait of the Rucklidge chaotic system.

The time response of before synchronization states for the Master system $\left({x}_{1}\mathrm{,}{y}_{1}\mathrm{,}{z}_{1}\right)$ and the Drive system $\left({x}_{2}\mathrm{,}{y}_{2}\mathrm{,}{z}_{2}\right)$ is shown in Figure 4.

The error dynamics in the uncontrolled state are shown in Figure 5, whereas the error dynamics in the controlled state are shown in Figure 6 and Figure 12.

Figure 3. On the *X*-*Z* plane, a 2-D phase portrait of the Rucklidge chaotic system.

Figure 4. Time response of the states $x\mathrm{,}y\mathrm{,}z$ before synchronization. (a) ${x}_{1}\mathrm{,}{x}_{2}$ ; (b) ${y}_{1}\mathrm{,}{y}_{2}$ ; (c) ${z}_{1}\mathrm{,}{z}_{2}$.

Figure 5. Time response of the error states before synchronization.

Figure 6. Time response of the synchronized states of error $\left({e}_{x}\mathrm{,}{e}_{y}\mathrm{,}{e}_{z}\right)$.

3. Synchronization via Active Control

When Master and Drive systems are in the same parameters, achieving synchronization using active control is appropriate and effective. Assume the existence of a Master system as

$\stackrel{\dot{}}{x}=Mx+g\left(x\right)$ (2)

where
$x={\left({x}_{1},{y}_{1},\cdots ,{z}_{1}\right)}^{\text{T}}\in {R}^{n}$,
$M\in {R}^{n}\times {R}^{n}$, *M* is a constant system matrix and
$g\left(x\right)$ is a sequence function that is nonlinear.

The Drive system is

$\stackrel{\dot{}}{y}=My+g\left(y\right)+u\left(t\right)$ (3)

where $\stackrel{\dot{}}{y}={\left({x}_{2},{y}_{2},\cdots ,{z}_{2}\right)}^{\text{T}}\in {R}^{n}$, $u\left(t\right)=\left({u}_{i}\left(t\right)\right)\in {R}^{n}$, where $i=\left(1,2,\cdots ,n\right)$.

3.1. Definition

If an appropriate controller $u\left(t\right)$ exists that satisfies $\forall x\mathrm{,}y\mathrm{,}e\in {R}^{n}$, ${\mathrm{lim}}_{t\to 0}\Vert y-x\Vert ={\mathrm{lim}}_{t\to 0}\Vert e\Vert =0$ the Master and Drive systems are then designed to be synchronized.

As a result, the error is defined as $e=y-x$.

Then the error dynamics is

$\stackrel{\dot{}}{e}=\stackrel{\dot{}}{y}-\stackrel{\dot{}}{x}=Me+G\left(x,y\right)+u\left(t\right)$ (4)

where
$G\left(x,y\right)=g\left(y\right)-g\left(x\right)$. Without *e* of system (4), controller
$u\left(t\right)$ may eliminate non-linear section. That is

$u\left(t\right)=v\left(t\right)-G\left(x,y\right)$ (5)

where $v\left(t\right)=Ke$ represents a linear section with error variables. From Equation (5) and Equation (4) we get

$\stackrel{\dot{}}{e}=Me+v\left(t\right)$ (6)

Equation (6) becomes when $v\left(t\right)$ is a linear section with error variables and $v\left(t\right)=ke$ is a constant matrix.

$\stackrel{\dot{}}{e}=\left(M+K\right)e$ (7)

3.2. Proposition [7]

Satisfing the requirements needed for diagonal matrix $\left(M+K\right)$ is ${\lambda}_{i}\le 0$, where ${\lambda}_{i}$ is the eigenvalue of matrix $\left(M+K\right)$, state vectors of system (7) is asymptotically converge to zero, as a result Master system (2) and Drive system (3) asymptotically synchronize. In this context, the Master and Drive systems are defined as follows:

The Master Rucklidge system is defined by

$\begin{array}{l}{\stackrel{\dot{}}{x}}_{1}=-b{x}_{1}+a{y}_{1}-{y}_{1}{z}_{1}\\ {\stackrel{\dot{}}{y}}_{1}={x}_{1}\\ {\stackrel{\dot{}}{z}}_{1}=-{z}_{1}+{y}_{1}^{2}\end{array}$ (8)

where *a* and *b* are the system parameters.

The Drive Rucklidge system is defined by

$\begin{array}{l}{\stackrel{\dot{}}{x}}_{2}=-b{x}_{2}+a{y}_{2}-{y}_{2}{z}_{2}+{u}_{1}\\ {\stackrel{\dot{}}{y}}_{2}={x}_{2}+{u}_{2}\\ {\stackrel{\dot{}}{z}}_{2}=-{z}_{2}+{y}_{2}^{2}+{u}_{3}\end{array}$ (9)

${u}_{i}=\left(1,2,3\right)$ are active control functions for the Master system. In Equation (9), we presented three control functions: ${u}_{1}\left(t\right)$, ${u}_{2}\left(t\right)$, and ${u}_{3}\left(t\right)$ which are must be regulated. We subtract Equation (8) from Equation (9) to estimate the control functions. Consider the state errors between the controlling system (9) and the master system (8) that needs to be controlled by using

${e}_{x}={x}_{2}-{x}_{1};{e}_{y}={y}_{2}-{y}_{1};{e}_{z}={z}_{2}-{z}_{1};$ (10)

Applying the active control design methods, we subtract Equation (8) from Equation (9) and use the definitions in Equation (10) to derive the error dynamics equation:

$\begin{array}{l}{\stackrel{\dot{}}{e}}_{x}=-b{e}_{x}+a{e}_{y}-{y}_{2}{z}_{2}+{y}_{1}{z}_{1}+{u}_{1}\\ {\stackrel{\dot{}}{e}}_{y}={e}_{x}+{u}_{2}\\ {\stackrel{\dot{}}{e}}_{z}=-{e}_{z}+{e}_{y}\left({y}_{2}+{y}_{1}\right)+{u}_{3}\end{array}$ (11)

tracking the re-describing control functions as:

$\begin{array}{l}{v}_{1}={u}_{1}-{y}_{2}{z}_{2}+{y}_{1}{z}_{1}\\ {v}_{2}={u}_{2}\\ {v}_{3}={u}_{3}\end{array}$ (12)

the error dynamics Equation (11) takes on a new state

$\begin{array}{l}{\stackrel{\dot{}}{e}}_{x}=-b{e}_{x}+a{e}_{y}+{v}_{1}\\ {\stackrel{\dot{}}{e}}_{y}={e}_{x}+{v}_{2}\\ {\stackrel{\dot{}}{e}}_{z}=-{e}_{z}+{e}_{y}\left({y}_{2}+{y}_{1}\right)+{v}_{3}\end{array}$ (13)

The controllable state of error system (11) is a linear system with control effort
$\left({v}_{1}\left(t\right)\mathrm{,}{v}_{2}\left(t\right)\mathrm{,}{v}_{3}\left(t\right)\right)$ as a function of error states
$\left({e}_{x}\mathrm{,}{e}_{y}\mathrm{,}{e}_{z}\right)$. The error states
$\left({e}_{x}\mathrm{,}{e}_{y}\mathrm{,}{e}_{z}\right)$ converge to zero as time
$t\to \infty $ passes if these stabilize the system. This means that active control is used to synchronize the Master and Drive systems. For the control
$\left({v}_{1}\left(t\right)\mathrm{,}{v}_{2}\left(t\right)\mathrm{,}{v}_{3}\left(t\right)\right)$, there are various possible sets. We created a constant matrix *A* that will govern the error dynamics (11) in accordance with the active control approach, such that

${\left[{v}_{1},{v}_{2},{v}_{3}\right]}^{\text{T}}=A{\left[{e}_{x},{e}_{y},{e}_{z}\right]}^{\text{T}}$ (14)

where *A* is a constant matrix. To stabilize the state of the error system, the components of the matrix *A* must be chosen in such a way that the feedback system has all eigenvalues with negative real portions. In the tracking form, choose the matrix *A*:

$A=\left(\begin{array}{ccc}-b& -a& 0\\ -1& -1& 0\\ 0& -\left({y}_{2}+{y}_{1}\right)& 0\end{array}\right)$

Bring Equation (14) into Equation (13), we may find

$\stackrel{\dot{}}{e}=Ke$

$\left(\begin{array}{c}{\stackrel{\dot{}}{e}}_{x}\\ {\stackrel{\dot{}}{e}}_{y}\\ {\stackrel{\dot{}}{e}}_{z}\end{array}\right)=\left(\begin{array}{ccc}-b& a& 0\\ 1& 0& 0\\ 0& {y}_{2}+{y}_{1}& -1\end{array}\right)\left(\begin{array}{c}{e}_{1}\\ {e}_{2}\\ {e}_{3}\end{array}\right)+A\left(\begin{array}{c}{e}_{1}\\ {e}_{2}\\ {e}_{3}\end{array}\right)$

$K=\left(\begin{array}{ccc}-2b& 0& 0\\ 0& -1& 0\\ 0& 0& -1\end{array}\right)$

In this specific selection, the eigenvalues of the system (13) are
$\left(-2b\mathrm{,}-\mathrm{1,}-1\right)$. Due to the linearity system’s stability theory, this option will result in the convergence of error states
$\left({e}_{x}\mathrm{,}{e}_{y}\mathrm{,}{e}_{z}\right)$ to zero as time *t* approaches infinity, and therefore the synchronization of the two systems is achieved under the control system (15)

$\begin{array}{l}{u}_{1}=-b{e}_{x}-a{e}_{y}+{y}_{2}{z}_{2}-{y}_{1}{z}_{1}\\ {u}_{2}=-{e}_{x}-{e}_{y}\\ {u}_{3}=-\left({y}_{2}+{y}_{1}\right){e}_{y}\end{array}$ (15)

3.3. Simulation and Results

To get numerical results, we apply MatLab Simulink and the Runge-Kutta algorithm of 4th order with a time grid of 0.05. We continue with the initial conditions of the Master system Equation (8) as follows:

${x}_{1}\left(0\right)=0.001,{y}_{1}\left(0\right)=0.001,{z}_{1}\left(0\right)=0.001$ (16)

and the Drive system’s initial conditions (9) as follows:

${x}_{2}\left(0\right)=5,{y}_{2}\left(0\right)=8,{z}_{2}\left(0\right)=4$ (17)

The time response of the error states obtained from numerical simulations of synchronization systems (11) is shown in Figure 6 under the controller (15) stimulated at time *t* = 100 seconds, thereby guaranteeing the synchronization of system (8) and (9). It is clear that after control signals are activated, the error vectors converge to zero quickly. We calculate the synchronization measure for verifying the synchronization act from Figure 7, the average error on the system state variables are given by

$e=\sqrt{{e}_{x}^{2}+{e}_{y}^{2}+{e}_{z}^{2}}$

and Figures 8-10 demonstrate the response time of the synchronization of the driven system (9) to the driver (8), control signals are activated at the time *t* = 0, and Figure 11 show the response time of action control (15) for achieve chaos synchronization between the two chaotic systems.

Figure 7. The response time of the synchronized states of average error.

Figure 8. The response time of the synchronized states $\left({x}_{1}\mathrm{,}{x}_{2}\right)$.

Figure 9. Time response of the synchronized states $\left({y}_{1}\mathrm{,}{y}_{2}\right)$.

Figure 10. Time response of the synchronized states $\left({z}_{1}\mathrm{,}{z}_{2}\right)$.

The purpose of this section Figures 8-10 is to exhaust active control procedure to synchronize two chaotic chaotic systems (8) and (9) by defining a controller that the Drive system capability to trail the Master System and the states of two chaotic systems (8) and (9) show similar activities for all future states [21] [41].

4. Backstepping Control Design for Chaos Synchronization

This section focuses on Backstepping Control systems because this class offers various examples of chaotic circuits and systems [42] we explore the backstepping control design for a general system of the form:

$\stackrel{\dot{}}{x}=f\left(x\right)+g\left(x\right){\mu}_{1}$ (18)

be a responsive affine nonlinear system, with
${\mu}_{1}\in R$ as the control input,
$x\in {R}^{n}$ as the state, and *f* and *g* as nonlinear functions, with
$f\left(0\right)=0$. The control system 18 can use the usual backstepping design. To stabilize systems that are in strict feedback form, we can use recursive application of backstepping control design as follows.

Figure 11. Time response of the synchronized states ${u}_{i}$. (a) ${u}_{1}$ ; (b) ${u}_{2}$ ; (c) ${u}_{3}$.

$\begin{array}{l}\stackrel{\dot{}}{x}=f\left(x\right)+g\left(x\right){\mu}_{1}\\ {\stackrel{\dot{}}{\mu}}_{1}={f}_{1}\left(x,{\mu}_{1}\right)+{g}_{1}\left(x,{\mu}_{1}\right){\mu}_{2}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\vdots \\ {\stackrel{\dot{}}{\mu}}_{k-1}={f}_{k-1}\left(x,{\mu}_{1},\cdots ,{\mu}_{k-1}\right)+{g}_{k-1}\left(x,{\mu}_{1},\cdots ,{\mu}_{k-1}\right){\mu}_{k}\\ {\stackrel{\dot{}}{\mu}}_{k}={f}_{k}\left(x,{\mu}_{1},\cdots ,{\mu}_{k}\right)+{g}_{k}\left(x,{\mu}_{1},\cdots ,{\mu}_{k}\right)u\end{array}$ (19)

where
${\mu}_{1}\mathrm{,}{\mu}_{2}\mathrm{,}\cdots \mathrm{,}{\mu}_{k}$ and *u* are scalars [43].

Backstepping design is a Lyapunov-based control approach that is well-organized. Backstepping control design of systems (18) is a recursive strategy that ensures the system’s global asymptotic stability. The fundamental concept is to extend Lyapunov’s method by breaking down the full system design model (19) into a series of design challenges for lower-order systems. Depending on the state variables, control parameters, and stabilizing functions, incorporate new variables into transformation processes. The *i*th subsystem may be stabilized with respect to a specific Lyapunov function
${V}_{i},i=1,2,\cdots ,n$ by utilizing the backstepping control design at the *i*th step. By treating the variable
${\mu}_{1}$ as a virtual control input, the approach is used to stabilize the first equation. Similarly, the second equation is stabilized by continuing to use the variable
${\mu}_{2}$ as the virtual control, when
${\mu}_{1}$ is designed and so on. So the design of the final control input *u*, which generally depends on *x* and
${\mu}_{1}\mathrm{,}{\mu}_{2}\mathrm{,}\cdots \mathrm{,}{\mu}_{k}$, is steadily achieved in *n* steps [43] [44]. So it is proved that backstepping design is suitable for controlling chaos, stabilization, and tracking problems.

Consider the following Master system of Rucklidge systems:

$\begin{array}{l}{\stackrel{\dot{}}{x}}_{1}=-b{x}_{1}+a{y}_{1}-{y}_{1}{z}_{1}\\ {\stackrel{\dot{}}{y}}_{1}={x}_{1}\\ {\stackrel{\dot{}}{z}}_{1}=-{z}_{1}+{y}_{1}^{2}\end{array}$ (20)

in relation to the Drive system as

$\begin{array}{l}{\stackrel{\dot{}}{x}}_{2}=-b{x}_{2}+a{y}_{2}-{y}_{2}{z}_{2}+u\\ {\stackrel{\dot{}}{y}}_{2}={x}_{2}\\ {\stackrel{\dot{}}{z}}_{2}=-{z}_{2}+{y}_{2}^{2}\end{array}$ (21)

where *u* is a to-be-defined control function. We required only one controller in backstepping method. By using the error states definition (10) and subtracting Equation (20) from (21), we obtain

$\begin{array}{l}{\stackrel{\dot{}}{e}}_{x}=-b{e}_{x}+a{e}_{y}-{z}_{2}{e}_{y}-{y}_{1}{e}_{z}+u\\ {\stackrel{\dot{}}{e}}_{y}={e}_{x}\\ {\stackrel{\dot{}}{e}}_{z}=-{e}_{z}+{e}_{y}\left({y}_{2}+{y}_{1}\right)\end{array}$ (22)

to re-arrange error system (22) into system (23) according to sequence of full system (19)

$\begin{array}{l}{\stackrel{\dot{}}{e}}_{z}=-{e}_{z}+{e}_{y}\left({y}_{2}+{y}_{1}\right)\\ {\stackrel{\dot{}}{e}}_{y}={e}_{x}\\ {\stackrel{\dot{}}{e}}_{x}=-b{e}_{x}+a{e}_{y}-{z}_{2}{e}_{y}-{y}_{1}{e}_{z}+u\end{array}$ (23)

The previous system (23) is an equilibrium point
$\left(\mathrm{0,0,0}\right)$ in the absence of all control *u*. If *u* is chosen, then the equilibrium point residues are unchanged, the challenge of synchronization between the Master and Drive systems can be transformed into a problem of asymptotical system stabilization (23). The goal of this project is to find a control law *u* that will stabilize the system’s error variables (23) at the origin. Since the system comprises three states or three non-linear differential equations, a recursive strategy will have three steps. We split down the system (23) into three subsystems, each with a single input and output. We initiate the formation with the first subsystem of the first non-linear differential equation, and one continues until the last subsystem. In the procedure of formation a change of coordinates
${w}_{i}=\left({e}_{x},{e}_{y},{e}_{z}\right)$ is done.

4.1. First Step

To begin, we use
${e}_{y}$ as a virtual controller to stabilize the first equation in (23), the object of the control is to drive
$\left({e}_{x},{e}_{y},{e}_{z}\right)=\left(0,0,0\right)$ *i.e..*
$\left({x}_{1}\mathrm{;}{y}_{1}\mathrm{;}{z}_{1}\right)=\left({x}_{2}\mathrm{;}{y}_{2}\mathrm{;}{z}_{2}\right)$, and define the first new virtual variable of the Backstepping design as
${w}_{1}={e}_{z}$ and Considering another new virtual variable as
${w}_{2}={e}_{y}-{\alpha}_{1}$. This last term will not be used in the first step, but it is needs to
${w}_{2}$ to join the first subsystem on
${w}_{1}$ to the next subsystem on
${w}_{2}$ which will be considered in the second step,
${\alpha}_{1}$ is the stabilization function. To find the stabilization function
${\alpha}_{1}$, consider the Lyapunov function [37], given by:

${V}_{1}\left({w}_{1}\right)=\frac{1}{2}{\left({w}_{1}\right)}^{2}$ (24)

substituting the derivative of ${w}_{1}$ for the time derivative of Equation (24) yields:

${\stackrel{\dot{}}{V}}_{1}=-{\left({w}_{1}\right)}^{2}+{w}_{1}{w}_{2}\left({y}_{2}+{y}_{1}\right)+{w}_{1}\left({y}_{2}+{y}_{1}\right){\alpha}_{1}$ (25)

If the estimated stabilisation function ${\alpha}_{1}$ is chosen as ${\alpha}_{1}=0$, as a result, ${\stackrel{\dot{}}{V}}_{1}$ is indefinitely negative. The second term in Equation (25) will be eradicated in the next step.

4.2. Second Step

We stabilize the second equation in (23) for the stability of the second subsystem and its dynamics are calculated by putting the second equation of the system (23) into the second virtual variable with derivatives:

${\stackrel{\dot{}}{w}}_{2}={\stackrel{\dot{}}{e}}_{y}-{\stackrel{\dot{}}{\alpha}}_{1}$

The stabilization function ${\alpha}_{2}$ for the second subsystem is preferred such as the new virtual state variable ${w}_{3}={e}_{x}-{\alpha}_{2}$. Let Lyapunov function for the second subsystem as:

${V}_{2}\left({w}_{1},{w}_{2}\right)={V}_{1}\left({w}_{1}\right)+\frac{1}{2}{\left({w}_{2}\right)}^{2}$ (26)

Substituting the derivative of ${w}_{2}$ in time derivative of ${\stackrel{\dot{}}{V}}_{2}$ then we get:

${\stackrel{\dot{}}{V}}_{2}=-{w}_{1}^{2}-{w}_{2}^{2}+{w}_{2}{w}_{3}+{w}_{2}\left({w}_{2}+{w}_{1}\left({y}_{2}+{y}_{1}\right)+{\alpha}_{2}\right)$ (27)

The stabilization function ${\alpha}_{2}$ is chosen as:

${\alpha}_{2}=-{w}_{2}-{w}_{1}\left({y}_{2}+{y}_{1}\right)$ (28)

And by substituting ${\alpha}_{2}\mathrm{,}{\stackrel{\dot{}}{V}}_{2}$ converted:

${\stackrel{\dot{}}{V}}_{2}=-{w}_{1}^{2}-{w}_{2}^{2}+{w}_{2}{w}_{3}$ (29)

then ${\stackrel{\dot{}}{V}}_{2}$ is negative definite, the third term in Equation (29) will be eradicated in the next step.

4.3. Third Step

To get the third virtual variable state, the time derivative of both the error on the third coordinate state and the stabilization function ${\alpha}_{2}$. The equation of the third subsystem become:

${\stackrel{\dot{}}{w}}_{3}={\stackrel{\dot{}}{e}}_{x}-{\stackrel{\dot{}}{\alpha}}_{2}$

${\stackrel{\dot{}}{w}}_{3}=-b{e}_{x}+a{e}_{y}-{z}_{2}{e}_{y}-{y}_{1}{e}_{z}+u-{\stackrel{\dot{}}{\alpha}}_{2}$ (30)

We choose the Lyapunov function for the third subsystem is as:

${V}_{3}\left({w}_{2},{w}_{3}\right)={V}_{2}\left({w}_{2}\right)+\frac{1}{2}{\left({w}_{3}\right)}^{2}$ (31)

It’s time derivative gives:

$\begin{array}{l}{\stackrel{\dot{}}{V}}_{3}=-{w}_{1}^{2}-{w}_{2}^{2}+{w}_{2}{w}_{3}+{w}_{3}{\stackrel{\dot{}}{w}}_{3}\\ {\stackrel{\dot{}}{V}}_{3}=-{w}_{1}^{2}-{w}_{2}^{2}-b{w}_{3}^{2}+{w}_{3}\left({w}_{2}-b{\alpha}_{2}+\left(a-{z}_{2}\right){e}_{y}-{y}_{1}{e}_{z}-{\stackrel{\dot{}}{\alpha}}_{2}+u\right)\end{array}$ (32)

If the control law *u* is preferred as follows

$u=-{w}_{2}+b{\alpha}_{2}-\left(a-{z}_{2}\right){e}_{y}+{y}_{1}{e}_{z}+{\stackrel{\dot{}}{\alpha}}_{2}$ (33)

Then,

${\stackrel{\dot{}}{V}}_{3}=-{w}_{1}^{2}-{w}_{2}^{2}-b{w}_{3}^{2}$ (34)

is negative definite. The final system of the time derivative of the error model in $\left({w}_{1}\mathrm{,}{w}_{2}\mathrm{,}{w}_{3}\right)$ coordinates in (35), the error dynamics $\left({e}_{x}\mathrm{,}{e}_{y}\mathrm{,}{e}_{z}\right)$ will converge to zero as $t\to \infty $, whereas the equilibrium $\left(\mathrm{0,0,0}\right)$ residues are asymptotically stable, according to the LaSalle-Yoshizawa theorem [43]. As a result, the master-drive system has been synchronized. We now study the complete space of $\left({\stackrel{\dot{}}{w}}_{1}\mathrm{,}{\stackrel{\dot{}}{w}}_{2}\mathrm{,}{\stackrel{\dot{}}{w}}_{3}\right)$ :

$\begin{array}{l}{\stackrel{\dot{}}{w}}_{1}=-{w}_{1}+{w}_{2}\left({y}_{2}+{y}_{1}\right)\\ {\stackrel{\dot{}}{w}}_{2}={e}_{x}\\ {\stackrel{\dot{}}{w}}_{3}=-b{e}_{x}+a{e}_{y}-{z}_{2}{e}_{y}-{y}_{1}{e}_{z}+u-{\stackrel{\dot{}}{\alpha}}_{2}\end{array}$ (35)

The equilibrium
$\left(\mathrm{0,0,0}\right)$ of system (35) is universally asymptotically stable, as illustrated in Figure 13, according to LaSalle-Yoshizawa theorem [26] [27]. That is, the control law *u* has no effect on the system’s equilibrium of (22), *i.e.*
$\left(\mathrm{0,0,0}\right)$ is still the equilibrium. As a result, the new chaotic system (35) is stabilized at the origin under the direction of the controller (33), and the two Rucklidge chaotic systems are synchronized.

4.4. Simulation and Results

We solved system (22) with the controllers defined in (33) using the fourth-order Runge-Kutta algorithm with initial conditions
${x}_{1}\left(0\right)=0.001$,
${y}_{1}\left(0\right)=0.001$,
${z}_{1}\left(0\right)=0.001$ and
${x}_{2}\left(0\right)=5$,
${y}_{2}\left(0\right)=8$,
${z}_{2}\left(0\right)=4$ with a time step of 0.01, and fixing the parameter values of *a* and *b* as shown in Figures 1-3. To demonstrate the efficiency of the control law, we show numerical findings. The time response of the error state
$\left({e}_{x}\mathrm{,}{e}_{y}\mathrm{,}{e}_{z}\right)$ is shown in Figure 12, which was

Figure 12. Tracking error ${e}_{x}\mathrm{,}{e}_{y}\mathrm{,}{e}_{z}$.

obtained from numerical simulations of the synchronization systems (20) and (21) under the controller (33).

From Figure 13 to examine the equilibrium $\left(\mathrm{0,0,0}\right)$ of system (35) and the average error propagation on the system state variables is used to calculate the synchronization quantity this are shown in Figure 14 given by

$e=\sqrt{{e}_{x}^{2}+{e}_{y}^{2}+{e}_{z}^{2}}$

After an initial transcience of around $t=6\text{\hspace{0.17em}}\text{s}$, the Master-Drive system is globally synchronized. The exponential convergence of the synchronization quality described by error propagation on average error states also confirms this in Figure 17(b).

Figure 15 shows the time response of the control law.

The persistence of this part Figure 16 is to synchronize two chaotic systems (20) and (21) shattering Backstepping control process by defining a controller, and for all future states, the effectiveness of the Drive system to track the Master System and the states of two chaotic systems (20) and (21) indicate similar characteristics.

5. Active Control and Backstepping Approaches Are Compared

We obtained a relationship between synchronization time to investigate and compare the synchronization performance of the two strategies. Figure 17 illustrates the synchronization error (*e*) for the two techniques when controls are activated at *t* = 0. So, for backstepping control, synchronization was achieved at
$t=6\text{\hspace{0.17em}}\text{s}$, and for active control, synchronization was achieved at (
$t=22\text{\hspace{0.17em}}\text{s}$ ), with a time delay of 16s. It is demonstrated that the error signals converge to the origin extremely smoothly with a low decay rate and sufficient synchronization speed, indicating that the researched controllers.

The active control and backstepping control techniques are used to create a controller for a three-dimensional autonomous chaotic system, with backstepping control reducing the number of controllers from three to one, and active control works according to dimension. So Backstepping control performs significantly, reducing the controller complexity and cost than the active control method.

Figure 13. Final system of the error model after synchronization.

Figure 14. Average error *e*.

Figure 15. Control effort *u*.

Figure 16. Time response of the states $x\mathrm{,}y\mathrm{,}z$ after synchronization. (a) Time response of $\left({x}_{1}\mathrm{,}{x}_{2}\right)$ ; (b) Time response of $\left({y}_{1}\mathrm{,}{y}_{2}\right)$ ; (c) Time response of $\left({z}_{1}\mathrm{,}{z}_{2}\right)$.

6. Conclusions

The performance of two control techniques for chaos synchronization, active control and recursive backstepping control, was investigated in this work. The two techniques are proven to have outstanding synchronization performance, with the active control slightly outperforming the backstepping. The performance of theoretically designed nonlinear controllers was verified by numerical simulations that confirmed the proposed controller’s effectiveness. The results are presented in graphical style, together with a time history (Figures 1-16). The summaries are as follows:

• Three controllers are found in active control design, and one controller is originated for the Backstepping design;

Figure 17. Comparison of synchronization times for (
$0\le t\le 100$ ) between (a) active control and (b) backstepping approach with controller activated. (a) Zoom of the time response of average error *e* in Active control; (b) Zoom of the time response of average error *e* in backstepping control.

• In two techniques, the error dynamics converge to zero as $t\to \infty $, hence the equilibrium point $\left(\mathrm{0,0,0}\right)$ remains asymptotically stable;

• Similar activities are shown in both ways for the states of two chaotic systems (Master and Drive);

• The results reveal that the backstepping strategy converges to zero faster than the active control technique for momentary error dynamics.

Conflicts of Interest

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

[1] |
Lorenz, E.N. (1963) Deterministic Nonperiodic Flow. Journal of Atmospheric Sciences, 20, 130-141. https://doi.org/10.1175/1520-0469(1963)020<0130:DNF>2.0.CO;2 |

[2] |
Pourmahmood, M., Khanmohammadi, S. and Alizadeh, G. (2011) Synchronization of Two Different Uncertain Chaotic Systems with Unknown Parameters Using a Robust Adaptive Sliding Mode Controller. Communications in Nonlinear Science and Numerical Simulation, 16, 2853-2868. https://doi.org/10.1016/j.cnsns.2010.09.038 |

[3] |
Kyrtsou, C. and Labys, W.C. (2006) Evidence for Chaotic Dependence between Us Inflation and Commodity Prices. Journal of Macroeconomics, 28, 256-266. https://doi.org/10.1016/j.jmacro.2005.10.019 |

[4] |
Kyrtsou, C. and Labys, W.C. (2007) Detecting Positive Feedback in Multivariate Time Series: The Case of Metal Prices and Us Inflation. Physica A: Statistical Mechanics and Its Applications, 377, 227-229. https://doi.org/10.1016/j.physa.2006.11.002 |

[5] |
Alligood, K.T., Sauer, T.D. and Yorke, J.A. (1996) Chaos. Springer, Berlin. https://doi.org/10.1007/b97589 |

[6] |
Luo, R., Su, H. and Zeng, Y. (2017) Chaos Control and Synchronization via Switched Output Control Strategy. Complexity, 2017, Article ID: 6125102. https://doi.org/10.1155/2017/6125102 |

[7] |
Lei, Y., Xu, W. and Xie, W. (2007) Synchronization of Two Chaotic Four-Dimensional Systems Using Active Control. Chaos, Solitons & Fractals, 32, 1823-1829. https://doi.org/10.1016/j.chaos.2005.12.014 |

[8] |
Matouk, A. and Elsadany, A. (2016) Dynamical Analysis, Stabilization and Discretization of a Chaotic Fractional-Order GLV Model. Nonlinear Dynamics, 85, 1597-1612. https://doi.org/10.1007/s11071-016-2781-6 |

[9] |
Zhang, H., Chen, D., Xu, B. and Wang, F. (2015) Nonlinear Modeling and Dynamic Analysis of Hydro-Turbine Governing System in the Process of Load Rejection Transient. Energy Conversion and Management, 90, 128-137. https://doi.org/10.1016/j.enconman.2014.11.020 |

[10] |
Zhang, H., Chen, D., Wu, C., Wang, X., Lee, J.-M. and Jung, K.-H. (2017) Dynamic Modeling and Dynamical Analysis of Pump-Turbines in S-Shaped Regions during Runaway Operation. Energy Conversion and Management, 138, 375-382. https://doi.org/10.1016/j.enconman.2017.01.053 |

[11] |
Zhang, H., Chen, D., Xu, B., Patelli, E. and Tolo, S. (2018) Dynamic Analysis of a Pumped-Storage Hydropower Plant with Random Power Load. Mechanical Systems and Signal Processing, 100, 524-533. https://doi.org/10.1016/j.ymssp.2017.07.052 |

[12] |
Zhang, H., Chen, D., Wu, C. and Wang, X. (2018) Dynamics Analysis of the Fast-Slow Hydro-Turbine Governing System with Different Time-Scale Coupling. Communications in Nonlinear Science and Numerical Simulation, 54, 136-147. https://doi.org/10.1016/j.cnsns.2017.05.020 |

[13] |
Zhang, H., Chen, D., Guo, P., Luo, X. and George, A. (2018) A Novel Surface-Cluster Approach towards Transient Modeling of Hydro-Turbine Governing Systems in the Start-Up Process. Energy Conversion and Management, 165, 861-868. https://doi.org/10.1016/j.enconman.2018.03.097 |

[14] |
Hamidzadeh, S. and Esmaelzadeh, R. (2014) Control and Synchronization Chaotic Satellite Using Active Control. International Journal of Computer Applications, 94, 29-33. https://doi.org/10.5120/16380-5887 |

[15] |
Shahverdiev, E. and Shore, K. (2009) Impact of Modulated Multiple Optical Feedback Time Delays on Laser Diode Chaos Synchronization. Optics Communications, 282, 3568-3572. https://doi.org/10.1016/j.optcom.2009.05.068 |

[16] |
Robert, B., Feki, M. and Iu, H.H. (2006) Control of a PWM Inverter Using Proportional plus Extended Time-Delayed Feedback. International Journal of Bifurcation and Chaos, 16, 113-128. https://doi.org/10.1142/S0218127406014629 |

[17] |
Pecora, L.M. and Carroll, T.L. (1991) Driving Systems with Chaotic Signals. Physical Review A, 44, 2374-2383. https://doi.org/10.1103/PhysRevA.44.2374 |

[18] |
Pecora, L.M. and Carroll, T.L. (1990) Synchronization in Chaotic Systems. Physical Review Letters, 64, 821-823. https://doi.org/10.1103/PhysRevLett.64.821 |

[19] |
Olusola, O., Vincent, U. and Njah, A. (2010) Synchronization, Multistability and Basin Crisis in Coupled Pendula. Journal of Sound and Vibration, 329, 443-456. https://doi.org/10.1016/j.jsv.2009.09.030 |

[20] |
Idowu, B.A., Olusola, O.I., Onma, O.S., Vaidyanathan, S., Ogabi, C.O. and Adejo, O.A. (2019) Chaotic Financial System with Uncertain Parameters—Its Control and Synchronisation. International Journal of Nonlinear Dynamics, and Control, 1, 271-286. https://doi.org/10.1504/IJNDC.2019.098682 |

[21] |
Njah, A. and Vincent, U. (2008) Chaos Synchronization between Single and Double Wells Duffing-Van der Pol Oscillators Using Active Control. Chaos, Solitons & Fractals, 37, 1356-1361. https://doi.org/10.1016/j.chaos.2006.10.038 |

[22] |
Vincent, U. (2008) Chaos Synchronization Using Active Control and Backstepping Control: A Comparative Analysis. Nonlinear Analysis: Modelling and Control, 13, 253-261. https://doi.org/10.15388/NA.2008.13.2.14583 |

[23] |
Li, C., Liao, X. and Wong, K.-W. (2005) Lag Synchronization of Hyperchaos with Application to Secure Communications. Chaos, Solitons & Fractals, 23, 183-193. https://doi.org/10.1016/j.chaos.2004.04.025 |

[24] |
Peng, C.-C., Hsue, A.W.-J. and Chen, C.-L. (2011) Variable Structure Based Robust Backstepping Controller Design for Nonlinear Systems. Nonlinear Dynamics, 63, 253-262. https://doi.org/10.1007/s11071-010-9801-8 |

[25] |
Yu, Y. and Zhang, S. (2004) Adaptive Backstepping Synchronization of Uncertain Chaotic System. Chaos, Solitons & Fractals, 21, 643-649. https://doi.org/10.1016/j.chaos.2003.12.067 |

[26] |
Yan, J.-J., Yang, Y.-S., Chiang, T.-Y. and Chen, C.-Y. (2007) Robust Synchronization of Unified Chaotic Systems via Sliding Mode Control. Chaos, Solitons & Fractals, 34, 947-954. https://doi.org/10.1016/j.chaos.2006.04.003 |

[27] |
Yao, Q. (2021) Synchronization of Second-Order Chaotic Systems with Uncertainties and Disturbances Using Fixed-Time Adaptive Sliding Mode Control. Chaos, Solitons & Fractals, 142, Article ID: 110372. https://doi.org/10.1016/j.chaos.2020.110372 |

[28] |
Zhang, H., Huang, W., Wang, Z. and Chai, T. (2006) Adaptive Synchronization between Two Different Chaotic Systems with Unknown Parameters. Physics Letters A, 350, 363-366. https://doi.org/10.1016/j.physleta.2005.10.033 |

[29] | Elabbasy, E., Agiza, H. and El-Dessoky, M. (2006) Global Chaos Synchronization for Four-Scroll Attractor by Nonlinear Control. Scientific Research and Essays, 1, 65-71. |

[30] |
Wang, X. and Wang, Y. (2011) Adaptive Control for Synchronization of a Four-Dimensional Chaotic System via a Single Variable. Nonlinear Dynamics, 65, 311-316. https://doi.org/10.1007/s11071-010-9893-1 |

[31] | Olusola, O.I., Vincent, E., Njah, A.N. and Ali, E. (2011) Control and Synchronization of Chaos in Biological Systems via Backsteping Design. International Journal of Nonlinear Science, 11, 121-128. |

[32] |
Stenflo, L. (1996) Generalized Lorenz Equations for Acoustic-Gravity Waves in the Atmosphere. Physica Scripta, 53, 83-84. https://doi.org/10.1088/0031-8949/53/1/015 |

[33] |
Mata-Machuca, J., Martinez-Guerra, R. and Aguilar-Lopez, R. (2010) An Exponential Polynomial Observer for Synchronization of Chaotic Systems. Communications in Nonlinear Science and Numerical Simulation, 15, 4114-4130. https://doi.org/10.1016/j.cnsns.2010.01.040 |

[34] |
Bai, E.-W. and Lonngren, K.E. (1997) Synchronization of Two Lorenz Systems Using Active Control. Chaos, Solitons & Fractals, 8, 51-58. https://doi.org/10.1016/S0960-0779(96)00060-4 |

[35] |
Ho, M.-C., Hung, Y.-C. and Chou, C.-H. (2002) Phase and Anti-Phase Synchronization of Two Chaotic Systems by Using Active Control. Physics Letters A, 296, 43-48. https://doi.org/10.1016/S0375-9601(02)00074-9 |

[36] |
Codreanu, S. (2003) Synchronization of Spatiotemporal Nonlinear Dynamical Systems by an Active Control. Chaos, Solitons & Fractals, 15, 507-510. https://doi.org/10.1016/S0960-0779(02)00128-5 |

[37] |
Jankovic, M., Sepulchre, R. and Kokotovic, P.V. (1996) Constructive Lyapunov Stabilization of Nonlinear Cascade Systems. IEEE Transactions on Automatic Control, 41, 1723-1735. https://doi.org/10.1109/9.545712 |

[38] |
Yassen, M. (2006) Chaos Control of Chaotic Dynamical Systems Using Backstepping Design. Chaos, Solitons & Fractals, 27, 537-548. https://doi.org/10.1016/j.chaos.2005.03.046 |

[39] |
Tan, X., Zhang, J. and Yang, Y. (2003) Synchronizing Chaotic Systems Using Backstepping Design. Chaos, Solitons & Fractals, 16, 37-45. https://doi.org/10.1016/S0960-0779(02)00153-4 |

[40] |
Tee, L.S. and Salleh, Z. (2013) Dynamical Analysis of a Modified Lorenz System. Journal of Mathematics, 2013, Article ID: 820946. https://doi.org/10.1063/1.4801119 |

[41] |
Vincent, U. (2008) Synchronization of Identical and Non-Identical 4-D Chaotic Systems Using Active Control. Chaos, Solitons & Fractals, 37, 1065-1075. https://doi.org/10.1016/j.chaos.2006.10.005 |

[42] |
Yang, T., Li, X.F. and Shao, H.H. (2001) Chaotic Synchronization Using Backstepping Method with Application to the Chua’s Circuit and Lorenz System. Proceedings of the 2001 American Control Conference, Arlington, VA, USA, 25-27 June 2001, Vol. 3, 2299-2300. https://doi.org/10.1109/ACC.2001.946094 |

[43] | Krstic, M., Kokotovic, P.V. and Kanellakopoulos, I. (1995) Nonlinear and Adaptive Control Design. John Wiley & Sons, Inc., Hoboken. |

[44] |
Khalil, H.K. (1996) Adaptive Output Feedback Control of Nonlinear Systems Represented by Input-Output Models. IEEE Transactions on Automatic Control, 41, 177-188. https://doi.org/10.1109/9.481517 |

Journals Menu

Contact us

customer@scirp.org | |

+86 18163351462(WhatsApp) | |

1655362766 | |

Paper Publishing WeChat |

Copyright © 2023 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.