Tracking Control of Chaotic Systems via Optimized Active Disturbance Rejection Control

For tracking control of chaotic systems, we develop an active disturbance rejection (ADR) control method. Using the first state of the system as the only available state, a time-varying bandwidth extended state observer reconstructs the remaining states and the total disturbance. A time-varying bandwidth feedback controller forces all the states of the system to follow exactly the reference signal and its derivative. The parameters of the ADR controller are optimized using a genetic algorithm. As the objective function, we chose the weighted sum of the integral of the absolute error and the integral of the absolute control signal. Two chaotic systems—the Duffing system and the Genesio–Tesi system—are considered in computer simulation tests. Results of these simulations are presented to demonstrate the effectiveness of the ADRC method in controlling chaotic systems.


Introduction
Chaos is the complex, unpredictable, and irregular behavior of systems.The behavior is very sensitive to a slight change in the initial conditions.Chaotic systems are very interesting nonlinear systems and have been intensively studied in numerous fields.Indeed, chaos can be found in many engineering systems such as oscillators, chemical reactions, robotics, lasers, and secure communications.Controlling a chaotic system requires designing a controller that stabilizes the system towards equilibrium points or forces the system to follow a time-varying reference signal.Transforming the system's response into useful signals is important in practice.Towards this end, a variety of control methods have been proposed, ever since the pioneering work reported in [1].Examples include sliding mode control [2][3][4][5][6][7][8], adaptive control [9][10][11][12][13][14][15][16][17][18], and backstepping control [19,20].
The control methods cited above are examples from two categories.The first assumes that, for a system model described accurately by a mathematical model, all states are available, neglecting uncertainties and external disturbances.This leads to a simple control method.The second assumes that the mathematical description of the system is partially known.This leads to a complex control method.In most cases in practice, only one state is available; the parameters of the system are uncertain and the system is subject to external disturbances.These uncertainties and external disturbances affect the performance of the controller.Therefore, it is essential to develop techniques that provide reliable estimates of the other states and the external disturbances, from the available single output.This problem can be solved if an observer is used.The observer design is a branch of control theory, for which solutions exist for linear and nonlinear systems.An observer is designed to reconstruct the unmeasured states and the external disturbance.
The active disturbance rejection (ADR) control method [21][22][23] is a nonlinear control method that is efficient and easy to implement.In addition, it is an effective method in dealing with uncertain nonlinear systems subject to external disturbances.In the (ADR) control method the unknown dynamics of the system and the external disturbances are considered the "total disturbance".The extended state observer (ESO) reconstructs the missing states and the total disturbance.These outputs of the ESO are used by a feedback controller.The procedure keeps the controller relatively simple, as the plant is reduced to a chain of integrators.

Mathematical Problems in Engineering
The objective of this paper is as follows: (i) To control chaotic systems with uncertain parameters, in the presence of external disturbances using ADRC method.
(ii) To optimize the ADRC controller parameters so that the system performs well by rejecting the total disturbances acting on the system and minimizing the tracking error to improve the overall system performance.
(iii) To suppress the peaking phenomena and improve the controller performance by using the time-varying observer and controller bandwidths.
As an objective function, a weighted sum of the integral of the absolute error and the integral of the absolute control signal is used, so that good tracking and concomitantly a smaller control signal can be produced.Two well-known examples are considered: the Duffing system (DS) as a typical secondorder chaotic system and the Genesio-Tesi system (GTS) as a typical third-order chaotic system.The organization of this article is as follows.Section 2 states the problem at hand in mathematical terms.Section 3 describes the ADRC method.In Section 4, simulation experiments are performed to show the effectiveness of the proposed scheme in controlling the Duffing and Genesio-Tesi chaotic systems.Section 5 concludes the paper.

Problem Formulation
In this study, the chaotic systems to be controlled are a class of n-th order single-input and single-output (SISO) continuous nonlinear systems, described by the dynamic equation: where  = [, ẋ , . . ., The aim of a control problem is to design a controller () to force the system output () to follow the given bounded reference signal () and its derivatives, that is, to minimize the tracking error defined as () = () − (), where ‖ ⋅ ‖ is the Euclidian norm (or 2-norm) of a vector.

ADR Controller
In practice, only partial information of the mathematical model is known for most control systems; that is, the system has uncertainties and external disturbance.The main part of ADR controller is the ESO [18].The main idea of the ESO is to use an additional state  +1 = (, ) + Δ(, ) + () so that system (1) can be written in state equation form as Here  +1 = (, ) is the total disturbance or the extended state of the system.Note that the extended state can be considered as total unknown factor including linear/nonlinear function of state and the disturbance part.
The ESO in its nonlinear version takes the form: With suitable values of  1 ,  2 , . . .,   ,  +1 the observer can track the states so that Note that a third-order observer is required when implementing an ADR controller for a second-order system.With a good estimation of all states of the system and the additional state, the control signal in the ADR controller is defined as where  1 (),  2 (), . . ..,   () are the time-varying gains of the controller.By substituting ( 7) into ( 5), the uncertain and disturbed system (1) can be reduced to a simple disturbancefree system represented as set of linear integrators: The objective in any control is to obtain an accurate estimation and then a cancellation of the total disturbance along with an accurate tracking of the reference.The speed at which the observer estimates the total disturbance (, ) to be cancelled by the controller in real time is crucial.
The procedure to determine the parameters of the ADR controller for the tracking control problem of the chaotic system considered in this study is explained in the next section.

Parameter Optimization of ADR Controller
4.1.Selecting Parameters.Parameter tuning in the ADR controller is usually a manual process.It is based on experience in the control field and is relatively difficult.Moreover, the numerous parameters requiring adjustments by manual methods weaken the control performance.The parameters of the ADR controller are both observer-related and controllerrelated.For the observer, one set of parameters that are required are the gains.With the parameterization technique proposed in [21], the roots of the characteristic polynomial are compared with polynomial thence placing all of the observers poles at − 0 and reducing the characteristic polynomial to Hurwitz-type.The same can be done to determine the controller gains  1 ,  2 , . . .  .The roots of the characteristic polynomial are compared with the following polynomial: so that all the poles are placed at -  , where   is the bandwidth of the feedback controller.A higher bandwidth results in better reference signal tracking and disturbance rejection.Nevertheless, the sensor noise and dynamic uncertainties are two factors that limit the maximum closed-loop bandwidth [19].
If the initial value of the extended observer is different from the initial value of the plant, the phenomenon of peaking appears for very large  0 .This can affect significantly the convergence of the extended observer.A possible way to mitigate peaking is to let the observer and the controller bandwidth vary with time as, for example, [24] where  and  are two positive constants.If   () =   ()/, with constant  that can be chosen from the range [2-10], then, the parameters of the ADR controller, , , and , can be optimized via an appropriate selection of an objective function.

Objective Function.
The selection of the objective function is important for accuracy during parameter setting.To achieve good tracking performance, and at the same time an acceptable control signal, a fitness function  is introduced that is a combination of the integral of the absolute error (IAE) and the integral of the absolute control signal (IAU).
The IAE performance index is adopted to obtain stability and accuracy of the tracking process; the IAU is included to minimize the control signal so that the need to limit the control signal using saturation function is eliminated.In practical applications, there are limiting conditions on the control energy that leads to a failure to achieve the desired performance.The minimization problem is formulated as where  is a positive constant.Setting  large leads to a more accurate tracking and a larger control signal.

Optimization Method.
Genetic algorithms (GAs) [22,23] are optimization methodologies based on Darwinian principles of evolutionary biology.They are often used to solve nonlinear or nondifferentiable optimization problems.As shown in Figure 1, a basic genetic algorithm consists of the following steps: 1. Initial population creation: randomly generate  chromosomes that are possible solutions of the problem.
(2) Calculating the fitness: by running the simulation obtain the fitness of the n chromosomes using the objective function.
(3) Generating the new population: by applying the GA operators until a new population is complete.
(4) Replacing: combine the old with newly generated solutions.
(5) Termination: if the termination criterion is satisfied, stop.
The GA operators in step new population.Selection methods include, for example, the roulette-wheel selection, random selection, and tournament selection.The crossover operation can be single-point, twopoint, or uniform.Mutation prevents the algorithm from being stuck in a local minimum.The mutation can be performed using procedures such as Flipping, Interchanging, and Reversing.The algorithm is terminated when the population converges through a variety of convergence criteria.A couple of popular ones are the fixed number of generations, and the best objective value is no longer changing.

Simulation Results
In this section, two well-known chaotic systems, the Duffing and Genesio-Tesi systems, are considered to demonstrate by numerical experiments the effectiveness of the proposed ADR control scheme.

Duffing System (DS).
A nonlinear oscillator with a cubic stiffness term to describe the hardening spring effect observed in many mechanical problems was introduced by Duffing.Duffing's equation has the form where  is the oscillation displacement,  1 the damping constant,  2 the linear stiffness constant,  3 the cubic stiffness constant,  the excitation amplitude, and  the excitation frequency.To establish the dynamic behavior of the chaotic DS, we select parameter settings:  1 = −1.1, 2 = 0.4,  3 = −1,  = 1.3, and  = 1.8 rad/s.From the DS response for the initial conditions  = 1, ẋ = 2 (Figure 2), the system develops chaotic behavior when no control signal is applied.
To test the performance of the ADR controller in its tracking task, the DS is controlled to follow a reference signal and its derivative given by  1 () = 10 sin ()  2 () = 10 cos () (16) By defining the states of ( 15),  1 =  and  2 = ẋ , it can be rewritten as two first-order ordinary differential equations, 1 () +  cos() as the extended state  3 , ( 17) can be rewritten as The ESO for a second-order system is as follows: The feedback controller is described by To optimize the three parameters , , and , the lower and upper bounds of the parameters are selected as (after a number of experimental simulations)  = [5,200],  = [0.1,1], and  = [2,10].The parameter settings of the GA are as follows: (i) a population size of 20-the tournament selection is applied to all chromosomes as the number of chromosomes in the population is small; (ii) a two-point crossover function and a uniform mutation function; (iii) during the optimization process, the stopping condition is invoked if the maximum number of generations reaches 100; and (iv) to increase the accuracy of tracking, a value of  = 100 is set in objective function.
The optimized controller parameters are then  = 170.0221, = 0.9985, and  = 1.300.The simulation results are presented in Figures 3-6 following the reference signal and its derivative after a short time.A good view of the tracking of the DS of the reference signal is seen in the phase-plane plot (Figure 4).The ESO performs well in estimating the two states of the system and the extended state in less than one second (Figure 5).The control effort required to force the DS follows the reference signal (Figure 6); note that the control signal required to control the DS is limited because the optimization of the ADR controller parameters based on the objective function ( 14) includes the minimization of the control signal as part of it.Therefore, the need to saturate the control signal is eliminated.The ADR controller parameter settings obtained from the minimization process leads to a very fast control of the system's chaotic behavior.

Genesio-Tesi System (GST).
As our second example, the GTS is considered having the form, Mathematical Problems in Engineering   where  1 ,  2 , and  3 are the system states with [ 1 ,  2 ,  3 ]  the full-state vector, , , , and  are unknown positive real constants satisfying  < , ( 1 ,  2 ,  3 , ) is an unknown function depending on the full state and time, () is the control signal to be designed with control input () = 0, ( 1 ,  2 ,  3 , ) = 0, and () = 0.The chaotic response of the GTS (Figure 7) was established using the initial conditions [ 10 = 5.3,  20 = −4, 30 = 2.5] for the variables and a running time of 20 seconds.
To test the performance of the ADRC in the task of tracking, the GTS is controlled to follow the reference signal and its derivatives given by Mathematical Problems in Engineering The uncertain term added is and the external disturbance is To estimate the three states and the total uncertainty, the following fourth-order ESO is required: The feedback controller is described by Using GA optimization, the control parameter settings obtained are  = 24.960, = 0.535,  = 1.985 using the same GA parameters and options adopted in the previous section.The simulation results are presented in Figures 8-12.The time-dependent responses of the states of the GTS (Figure 8) show again that the states exactly follow the reference signal and its derivative after a short time.Figure 9 shows the control  effort required to follow the reference signal.The ESO also performs well in estimating the three states of the system and the extended state (Figure 10). Figure 11 shows the controlled trajectories of GTS in the  1  2 and  1  3 planes.Figure 12 shows the time-varying observer and controller bandwidths.

Conclusion
An ADR controller has been developed to control chaotic systems and applied in two well-known examples to demonstrate its capabilities.The merit of this approach is the use of time-varying parameters for the observer as well as for the feedback controller.The results obtained in simulations clearly show the controller performs well in controlling both chaotic systems with bounded uncertainty and disturbance to any arbitrarily desired reference signal with high accuracy.
The results also show that the ADR method holds promise in providing solutions in controlling chaos.The ADR has been applied to control a second-order DS and a third-order GTS.
A suitable selection of the control parameter settings was observed to influence the control effort and error.Hence, the genetic algorithm was applied to obtain optimal parameter time (s) settings of the ADR controller.In further study, the ADR controller is to be applied to the synchronization of chaos.

Figure 2 :Figure 3 :
Figure 2: Chaotic behavior of DS without the control input in 20 seconds.

Figure 6 :
Figure 6: The time response of the control signal u(t).

Figure 8 :
Figure 8: Time responses of the states.

Figure 10 :Figure 11 :
Figure 10: The performance of the ESO in estimating the states.

Figure 12 :
Figure 12: The time-varying observer and controller bandwidths.