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.
1. 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–8], adaptive control [9–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–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.
The objective of this paper is as follows:
To control chaotic systems with uncertain parameters, in the presence of external disturbances using ADRC method.
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.
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 second-order 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.
2. 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:(1)x˙1=x2x˙2=x3⋮x˙n=fX,t+ΔfX,t+dt+uty=x1twhere X=[x,x˙,…,xn-1]T=[x1,x2,…,xn]T∈Rn is the state vector, fX,t is function of the states and time, not exactly known but assumed to be continuous and bounded, y∈R is the output, and u(t)∈R is the control signal. The functions Δf(X,t),d(t) are the bounded uncertain term and the disturbance term, respectively. In mathematical form,(2)ΔfX,t≤α,dt≤βwhere αandβ are positive constants. The given reference vector Rt is assumed to be bounded and has up to (n−1) bounded derivatives; that is,(3)Rt=r,r˙,…,rn-1T=r1,r2,…,rnTThe aim of a control problem is to design a controller u(t) to force the system output y(t) to follow the given bounded reference signal R(t) and its derivatives, that is, to minimize the tracking error defined as Et=X(t)-R(t),(4)limt→+∞Et=Xt-Rt→0where · is the Euclidian norm (or 2-norm) of a vector.
3. 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 xn+1=f(X,t)+Δf(X,t)+d(t) so that system (1) can be written in state equation form as(5)x˙1=x2x˙2=x3⋮x˙n=xn+1+utx˙n+1=F˙X,ty=x1tHere xn+1=F(X,t) 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:(6)x^˙1=x^2+α1ωoty-x^1x^˙2=x^3+α2ωo2ty-x^1⋮x^˙n=x^n+1+αnωonty-x^1+utx^˙n+1=αn+1ωon+1y-x^1With suitable values of β1,β2,…,βn,βn+1the observer can track the states so that x^1≈x1,x^2≈x2,…,x^n≈xn,x^n+1≈xn+1=F(X,t).
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(7)ut=-x^n+1+k1te^1+k2te^2+⋯+knte^nwhere k1(t),k2(t),….,kn(t) 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 disturbance-free system represented as set of linear integrators: (8)x˙1=x2x˙2=x3⋮x˙n=k1e^1+k2e^2+⋯+kne^nThe 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 F(X,t) 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.
4. Parameter Optimization of ADR Controller4.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 controller-related. 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(9)λs=sn+1+β1sn+β2sn-1+⋯+βn-1s2+βns+βn+1are compared with polynomial(10)Gos=s+ωotn+1thence 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 k1,k2,…kn. The roots of the characteristic polynomial(11)λcs=sn+knsn-1+⋯+k2s+k1are compared with the following polynomial: (12)Gcs=s+ωctnso that all the poles are placed at –ωc, where ωc 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](13)ωot=A1-e-Bt1+e-Btwhere A and B are two positive constants. If ωct=ωo(t)/C, with constant C that can be chosen from the range [2–10], then, the parameters of the ADR controller, A,B,andC, can be optimized via an appropriate selection of an objective function.
4.2. 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 J 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(14)JA,B,C=λ∫0te1tdt+∫0tutdtwhere λ is a positive constant. Setting λ large leads to a more accurate tracking and a larger control signal.
4.3. 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:
Initial population creation: randomly generate n chromosomes that are possible solutions of the problem.
Calculating the fitness: by running the simulation obtain the fitness of the n chromosomes using the objective function.
Generating the new population: by applying the GA operators until a new population is complete.
Replacing: combine the old with newly generated solutions.
Termination: if the termination criterion is satisfied, stop.
Loop: go to Step (2) for fitness evaluation.
Basic genetic algorithm.
The GA operators in step (3) are (a) Selection: based on their fitness, select two parent chromosomes, (b) Crossover: Cross over the parents to form new children, (c) Mutation: mutate new offspring at a randomly selected position in the chromosome, and (d) Acceptance: place new children in the new population. Selection methods include, for example, the roulette-wheel selection, random selection, and tournament selection. The crossover operation can be single-point, two-point, 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.
5. 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.
5.1. 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(15)x¨+a1x˙+a2x+a3x3=ut+qcosωtwhere x is the oscillation displacement, a1 the damping constant, a2 the linear stiffness constant, a3 the cubic stiffness constant, q the excitation amplitude, and ω the excitation frequency. To establish the dynamic behavior of the chaotic DS, we select parameter settings: a1=-1.1,a2=0.4,a3=-1,q=1.3,andω=1.8rad/s. From the DS response for the initial conditions x=1,x˙=2 (Figure 2), the system develops chaotic behavior when no control signal is applied.
Chaotic behavior of DS without the control input in 20 seconds.
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(16)r1t=10sintr2t=10costBy defining the states of (15), x1=xandx2=x˙, it can be rewritten as two first-order ordinary differential equations,(17)x˙1t=x2tx˙2t=a1x1t+a2x2t+a3x13t+qcosωt+utyt=x1tTaking Fx1,x2,ω,t=a1x1t+a2x2t+a3x13t+qcos(ωt) as the extended state x3, (17) can be rewritten as(18)x˙1=x2x˙2=x3+utx˙3=F˙x1,x2,ω,tThe ESO for a second-order system is as follows:(19)x^˙1t=x^2t+3ωotyt-x^1tx^˙2t=x^3t+3ωo2tyt-x^1tx^˙3t=ωo3tyt-x^1tThe feedback controller is described by(20)ut=-x^3+ωc2te^1+2ωce^2To optimize the three parameters A, B, and C, the lower and upper bounds of the parameters are selected as (after a number of experimental simulations) A=5,200,B=0.1,1,andC=[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 A=170.0221,B=0.9985,andC=1.300. The simulation results are presented in Figures 3–6. The time-dependent responses of the states of the DS (Figure 2) show the states exactly 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.
Controlled time response of the states x1,x2.
The controlled trajectory of the DS system.
The performance of the ESO in estimating the states.
The time response of the control signal u(t).
5.2. Genesio–Tesi System (GST)
As our second example, the GTS is considered having the form,(21)x˙1=x2x˙2=x3x˙3=-ax1-bx2-cx3+mx12+Δfx1,x2,x3,t+dt+utwhere x1,x2,andx3 are the system states with x1,x2,x3T the full-state vector, a,b,c,andm are unknown positive real constants satisfying ab<c, ∆fx1,x2,x3,t is an unknown function depending on the full state and time, u(t) is the control signal to be designed with control input ut=0, ∆fx1,x2,x3,t=0, and dt=0. The chaotic response of the GTS (Figure 7) was established using the initial conditions [x10=5.3,x20=-4,x30=2.5] for the variables and a running time of 20 seconds.
Chaotic trajectories of Genesio–Tesi system (GTS).
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(22)r1t=2cost+sintr2t=cost-2sintr3t=-2cost-sintThe uncertain term added is(23)Δfx1,x2,x3,t=0.05sinπx1sin2πx2sin3πx3and the external disturbance is(24)dt=0.002costTo estimate the three states and the total uncertainty, the following fourth-order ESO is required:(25)x^˙1t=x^2t+41ωotx1t-x^1tx^˙2t=x^3t+6ωo2tx1t-x^1tx^˙3t=x^4t+4ωo3tx1t-x^1t+utx^˙4t=α4ωo4tx1t-x^1tThe feedback controller is described by(26)ut=-x^4+ωc3te^1+3ωc2te^2+3ωcte^3Using GA optimization, the control parameter settings obtained are A=24.960,B=0.535,C=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 x1x2andx1x3 planes. Figure 12 shows the time-varying observer and controller bandwidths.
Time responses of the states.
Time response of control input (u(t)).
The performance of the ESO in estimating the states.
The controlled trajectories of the GT system.
The time-varying observer and controller bandwidths.
6. 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 settings of the ADR controller. In further study, the ADR controller is to be applied to the synchronization of chaos.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The author declares that there are no conflicts of interest.
PecoraL. M.CarrollT. L.JohnsonG. A.MarD. J.HeagyJ. F.Fundamentals of synchronization in chaotic systems, concepts, and applications19977452054310.1063/1.166278MR1604666Zbl0933.37030ZhangX.LiuX.ZhuQ.Adaptive chatter free sliding mode control for a class of uncertain chaotic systems201423243143510.1016/j.amc.2014.01.094MR3181281Zbl06893637WangB.ZhongS.Observer-based control on a chaotic system with unknowns and uncertainties201713716717410.1016/j.ijleo.2017.02.079LuoR.SuH.Finite-time control and synchronization of a class of systems via the twisting controller20175562199220710.1016/j.cjph.2017.09.003MR37324492-s2.0-85037042365DadrasS.MomeniH. R.MajdV. J.Sliding mode control for uncertain new chaotic dynamical system2009414185718622-s2.0-6734919222510.1016/j.chaos.2008.07.054Zbl1198.34114SunJ.WangY.WangY.ShenY.Finite-time synchronization between two complex-variable chaotic systems with unknown parameters via nonsingular terminal sliding mode control20168521105111710.1007/s11071-016-2747-8MR3511427Zbl1355.340972-s2.0-84961639091SuX.LiuX.ShiP.YangR.Sliding Mode Control of Discrete-Time Switched Systems with Repeated Scalar Nonlinearities2017629460446102-s2.0-8502849584710.1109/TAC.2016.2626398Zbl1390.93508SunJ.WuY.CuiG.WangY.Finite-time real combination synchronization of three complex-variable chaotic systems with unknown parameters via sliding mode control2017883167716902-s2.0-8501288821810.1007/s11071-017-3338-zZbl1380.34097ChangK.-M.Adaptive control for a class of chaotic systems with nonlinear inputs and disturbances200836246046810.1016/j.chaos.2006.06.090MR2382575WenJ.JiangC.-S.Adaptive fuzzy control for a class of chaotic systems with nonaffine inputs201116147549210.1016/j.cnsns.2010.03.015Zbl1221.931602-s2.0-77955918760HuaC.GuanX.LiX.ShiP.Adaptive observer-based control for a class of chaotic systems200422110311010.1016/j.chaos.2003.12.072Zbl1089.930102-s2.0-1842528955ZhouJ.ErM. J.Adaptive output control of a class of uncertain chaotic systems200756645246010.1016/j.sysconle.2006.12.002MR2319143Zbl1117.930412-s2.0-34047246340JiaN.WangT.Chaos control and hybrid projective synchronization for a class of new chaotic systems201162124783479510.1016/j.camwa.2011.10.069MR2855623DadrasS.MomeniH. R.Control uncertain Genesio-Tesi chaotic system: Adaptive sliding mode approach2009425314031462-s2.0-6765121069210.1016/j.chaos.2009.04.018ZhangD.MeiJ.MiaoP.Global finite-time synchronization of different dimensional chaotic systems20174830331510.1016/j.apm.2017.04.009MR3658856HandaH.SharmaB. B.Novel adaptive feedback synchronization scheme for a class of chaotic systems with and without parametric uncertainty201686506310.1016/j.chaos.2016.02.020MR3479677XingJinZhangFangfangAn Approach of Tracking Control for Chaotic Systems20162016810.1155/2016/97352649735264MR3506127Zbl1347.93142ChangY.-C.A robust tracking control for a class of uncertain chaotic systems20131561752176310.1002/asjc.740MR3130251Zbl1280.930252-s2.0-84886723949YuJ.LeiJ.WangL.Backstepping synchronization of chaotic system based on equivalent transfer function method201713090091310.1016/j.ijleo.2016.11.007LeiJ.WangL.Backstepping synchronous control of chaotic system with reduced number of active inputs201612723113641137310.1016/j.ijleo.2016.09.003HanJ. Q.From PID to active disturbance rejection control200956390090610.1109/TIE.2008.20116212-s2.0-63149185343GaoZ.Scaling and bandwidth-parameterization based controller tuningProceedings of the American Control ConferenceJune 2003Denver, Colo, USA498949962-s2.0-0142199923ZhengQ.GaoZ.Predictive active disturbance rejection control for processes with time delay201453487388110.1016/j.isatra.2013.09.0212-s2.0-84905055201LiuJ.2012Beijing, ChinaTsinghua University Press