This paper proposes an application of linear flatness control along with active disturbance rejection control (ADRC) for the local stabilization and trajectory tracking problems in the underactuated ball and rigid triangle system. To this end, an observer-based linear controller of the ADRC type is designed based on the flat tangent linearization of the system around its corresponding unstable equilibrium rest position. It was accomplished through two decoupled linear extended observers and a single linear output feedback controller, with disturbance cancelation features. The controller guarantees locally exponentially asymptotic stability for the stabilization problem and practical local stability in the solution of the tracking error. An advantage of combining the flatness and the ADRC methods is that it possible to perform online estimates and cancels the undesirable effects of the higher-order nonlinearities discarded by the linearization approximation. Simulation indicates that the proposed controller behaves remarkably well, having an acceptable domain of attraction.
Instituto Politécnico Nacional20160268201616371. Introduction
In the last 3 decades, there has been increasing interest in the control of underactuated mechanical systems. These systems are characterized by having fewer actuators than degrees of freedom [1]. This class of mechanical systems is becoming popular in many control applications, such as spacecraft, aerial robotic systems, underwater vehicles, and flexible robotics, to mention a few [2–4]. Possible advantages associated with these systems are cost reduction, lighter structures, and smaller dimensions, among others (see [5–7] for a comprehensive treatment of this class of systems). In this context, we are interested in the output feedback asymptotic stabilization and the output feedback trajectory control problems of the uncertain ball and rigid triangle (BRT) system. This system consists of a rigid isosceles triangle, OAB, hung from the vertex, O, and able to rotate around this vertex by the action of the torque, τ (see Figure 1). The identically rigid beams, OA¯ and OB¯, of mass m0 and length R, are attached to the left and right ends, respectively, of a rigid homogeneous beam, AB¯, of length 2L and mass M. Over the AB¯ beam, there is a ball of mass m, which can move freely forwards and backwards on a direction parallel to it. As is evident, this ball moves when the triangle rotates and the AB¯ beam bends either side. As the figure shows, the system is quite nonlinear due to the gravitational forces, the Coriolis and centripetal forces, and the acceleration couplings. On the other hand, by simple comparison, it is easy to see that this system is quite similar to but more complex than the traditional ball and beam system [6], because both dynamics of this system are completely coupled. Therefore, this system is not feedback-linearizable. In addition, it does not have a well defined relative degree [8].
The ball and rigid triangle system.
Fortunately, the linearized tangent model of this system is locally controllable around the unstable equilibrium point, meaning that it is locally flat. Hence, the stabilization and tracking problems both can be solved locally, from a combined perspective of flatness approach and active disturbance rejection control (ADRC). To accomplish this, we can propose two decoupled extended linear observers, assuming that only the position variables of this system are available. These observers allow us to simultaneously estimate the time derivatives of the nonavailable flat output and recover the uncertain underlying nonlinear dynamics. Subsequently, these estimations, together with the ADRC approach, allow us to propose a control scheme to solve the aforementioned control problems. The main difference between flatness controllers [9–12] and the ADRC scheme for flat systems is the fact that traditional flatness-based controls need perfect knowledge of the plant, while ADRC schemes for flat systems may largely ignore unknown nonlinearities and exogenous additive perturbation inputs in the input-to-flat output dynamics [13, 14]. We underscore that the brilliant idea of designing a controller with the capacity of lumped compensation of endogenous and exogenous disturbances, by means of an observer-based control, was proposed by Han [15]. He was the first to introduce the concept of the ADRC. ADRC controllers have led to new paradigmatic traditional nonlinear control problems in which disturbances, internal and external, are actively estimated and rejected [16, 17]. This concept has been used lately for controlling some kinds of underactuated systems in conjunction with the differential flatness approach with promising results (see [14, 18–22]). Recall that we can use this approach for systems in which tangent linear approximation is locally controllable.
In this article, an ADRC scheme with a flatness-based approach is proposed to practically solve the output feedback stabilization and the output feedback trajectory tracking problem for the BRT system. Our control approach assumes a lack of knowledge of the system parameters, the nonlinearities, and exogenous disturbance signals. The scheme not only estimates the unknown dynamics and the unknown state variables, but also reduces the tracking control problem to that defined on a chain of integrators after online active disturbance cancelations. Finally, the control algorithm is tested in several numerical simulations, showing excellent results for the output stabilization problem and output tracking problem. It is worth mentioning that the basin idea of this methodology is sustained by the use of high-gain observers, which are used as identifiers of the uncertain dynamics, assuming that in some operational region this uncertainty can be considered as being smooth and bounded. This kind of methodology has been used in [23–26].
The rest of this study is organized as follows. In Section 2, we obtain the nonlinear model of the BRT system by using the Euler Lagrange method. We normalize the obtained nonlinear model; then we linearize it around to its unstable equilibrium, and we show that the incremental linear model is flat with respect to a suitable output. In Section 3, we briefly describe the ADRC approach and introduce a useful theorem related to its stability. In Section 4, we propose the two decoupled extended observers and design a suitable ADRC controller. The results of the numerical simulations showing the effectiveness of our proposal are presented in Section 5, while Appendix is devoted to final remarks and future work.
Notation. In the forthcoming developments, we use the following configuration:(1)x˙=Aix+Eiu;y=CiTx,where x∈Ri denotes the system state, u∈R is the system single input, and matrices Ai∈Ri×i,Ei∈Ri, and Ci∈Ri are defined according to the well-known Brunovsky form. Finally, λi{X} stands for the eigenvalues of matrix X.
2. The Ball and Rigid Triangle (BRT) Model
Obtaining the Mathematical Model. We first consider a fixed reference frame, located at the fixed point O, with the horizontal and vertical unidirectional vectors, given by i→ and j→, respectively (see Figure 1). The generalized coordinates are denoted by the vector q=[q1,q2]. The coordinate q1<L is the cart position, measured from itself to the gravity center of the beam. The coordinate q2<π/2 is the angle formed by the vertical line passing through O and the imaginary line passing through O and the beam gravity center. To simplify the dynamical model and capture its nonlinear nature, the masses of OA¯ and OB¯ beams are considered equal to zero. Evidently, the distance between the fixed point O and the gravity center of the AB¯ beam is LR=R2-L2. After some physical considerations, it is easy to show that the kinetic energy of the triangle structure is defined as(2)KMq˙=12Ibq˙22+12MLR2q˙22,where Ib is the total moment of inertia of the AB¯ beam. Similarly, the kinetic energy related to the mass, m, is(3)Kmq,q˙=12mLR2q˙22+q12q˙22+2LRq˙1q˙2+q˙12.Therefore, the total kinetic energy of the system is(4)Kq,q˙=Kmq,q˙+KMq˙.The potential energies associated with the beam and the moving mass are, respectively, given by(5)Pmq=Mg1-cosq2;Pmq=mgLR1-cosq2+q1sinq2.In addition, the total potential energy of the system is defined as(6)Pq,q˙=Pmq+PMq.From (4) and (6), we can write the Lagrangian function as(7)Lq,q˙=Kq,q˙-Pq.Therefore, the Euler Lagrange equations are(8)mq¨1+mR0q¨2-mq1q˙22+mgsinq2=0;mR0q¨1+Iqq¨2+2mq1q˙1q˙2+M+mgR0sinq2+mgq1cosq2=τ,where(9)Iq=Ib+mLR2+mq12.To simplify the algebraic manipulations in the forthcoming developments, the above equations are normalized using the following transformations:(10)x1=q1LR;x2=q2;u=τmg;ϵ=tgLR,leading to the following simpler system:(11)x¨1+x¨2-x1x˙22+sinx2=0;x¨1+1+x12+Nx¨2+2x1x˙1x˙2+1+μsinx2+x1cosx2=u,with N=Ib/(mLR2) and μ=M/m, where “·” stands for differentiation with respect to the dimensionless ϵ.
Evidently, the system above admits the following Euler-Lagrange representation:(12)Mnxx¨+Cnx,x˙x˙+Gnx=U,where x=(x1,x2) and(13)Mnx=1111+x12+N;Gnx=∂∂xVnx;Cnx,x˙=0x1x˙2x1x˙2x1x˙1;U=0u,and Vnx=-1+μcosx2+x1sinx2. Finally, it is easy to see that system (8) verifies the following properties: (i) Mn(x)>0; (ii) the matrix S=Mn(x)-2Cn(x,x˙) is skew symmetric; and (iii) the operator u→x˙2 is passive; that is, the energy function time derivative E˙n=x˙2u, where En=x˙TMnxx˙/2+Vnx (notice that the sign of En is not well defined).
Problem Statement. Given the system model described in (8), the objective of this work consists of two goals: (1) designing a linear stabilizing controller to simultaneously bring the rigid triangle and the cart to the zero position x1=x2=0, assuming that all variables are initialized within a small vicinity of the origin, and (2) solving locally the output feedback trajectory tracking problem.
Motivation. It is easy to see that system (11) is unstable around the zero equilibrium point. In addition, this system cannot be feedback-linearizable using either static or dynamic state feedback. This fact is easy to verify by simply comparing the structure of the BRT system with the structure of the complete model of the ball and beam system [6]. On the other hand, we strongly believe that the problem of stabilizing the system (11) has not been solved, until now, using the energy shape based control method (see [27–29]). Additionally, there has been no analytical expression for the closed loop energy function. Here, we provide an approximated control solution by considering a linear version of the original model of the BRT system, in conjunction with the ADRC approach. In our opinion, this mechanical system can be applied to balance loads. For instance, the model can be extended to a three-dimensional case, which represents a container inside which a heavy object can slide freely due to load and unload maneuvers, where we want to prevent the object from hitting and damaging. Finally, we mention that this model can be considered as a simplification of the ball and plate system or the ball and beam system if the radius ball is considered in the model [30, 31].
2.1. Flatness of the Linearized Model
Consider the tangent linearization of system (11) around the following unstable equilibrium points:(14)x1=0,x˙1=0,x2=0,x˙2=0,u=0.This leads to the following linearized system:(15)x¨1δ+x¨2δ+x2δ=0;x¨1δ+1+Nx¨2δ+1+μx2δ+x1δ=uδ,where x1δ=x1 and x2δ=x2 are the incremental states of the linearized system. Then, the incremental flat output, Fδ, is(16)Fdelta=x1δ+x2δ.Indeed, all variables can be parameterized in terms of Fδ and a finite number of its time derivatives, as follows:(17)x1δ=Fδ+Fδ2;x2δ=-Fδ2;x˙1δ=Fδ1+Fδ3;x˙2δ=-Fδ3;uδ=-Fδ4N+x1+μx2.
2.2. A Flatness-Based Pole Placement Approach for Stabilization
The state-dependent input-coordinate transformation,(18)uδ=-vδN+x1+μx2,shows that system (15) is equivalent to the following chain of integrators:(19)Fδ4=vδ,where vδ is the new stabilizing feedback controller, fixed as(20)vδ=F∗4-∑i=04kiFδi-F∗i,where the set of coefficients k1,k2,k3,k4 is selected, such that the closed loop characteristic polynomial,(21)ps=s4+k4s3+k3s2+k2s+k1is Hurwitz, and F∗=F∗t is the desired reference signal, with bounded time derivatives up to order 4. For a detailed treatment of the flatness approach, we suggest [9, 11].
Remark 1.
Recall that the time derivatives Fδ(i),i={0,1,2,3} are given in the original coordinates as(22)Fδ=x1δ+x2δ;Fδ1=x˙1δ+x˙2δ;Fδ2=-x2δ;Fδ3=-x˙2δ.In the following section, we develop a control strategy for the BRT system, based on the ADRC.
Numerical Simulation. To show the closed loop response of the nonlinear normalized system (11) to the derived flatness controller (18), two simulations were carried out. To this end, the physical parameters were fixed as(23)m=0.1Kg;M=0.5Kg;Lr=0.5m;L=0.5m;Ib=0.13Kgm2,while the closed loop characteristic polynomial was chosen to be(24)ps=s2+2ζωns+ωn22.
3. ADRC to Locally Control the Ball and Rigid Triangle (BRT) System
In this section, we solve once again the previously mentioned control problem, assuming that only the position is available for measurements and that some dynamics of the original system are unknown. Before developing our control strategy, we give a summary on the ADRC subject.
3.1. Summary of ADRC
Let us consider the uncertain plant of nth-order, described by the following set of nonlinear differential equations:(25)x˙=Anx+Enbu+fx,d;y=CnTx,where x=x1,x2,…,xnT∈Rn is the state; u and y∈R are, respectively, the input and output; b is a known constant; d is the external unknown disturbance; and f(·) stands for the unknown nonlinear time-varying dynamics. Next, consider the following extended observer (EOS) for system (25):(26)x^˙e=An+1x^e+Bn+1u+ly-y^;y^e=Cn+1Tx^e,where BnT=[0,0,…,1,0]∈Rn, and the estimation vector x^e∈Rn+1 is defined as(27)x^e≔x^1,x^2,…,x^n,x^n+1T=x^,f^T,and l is the gains vector of the observer, set as(28)l=w0α1,w02α2,…,w0n+1αn+1T,where w0>0, with constants αi defined as(29)αi≔n+1!i!n+1-i!;i=1,2,…,n+1.According to this fact, it is clear that the matrix,(30)A_=-α110⋯0-α201⋯0⋮⋮⋮⋯0⋮⋮⋯1-αn+10⋯0,is Hurwitz.
Theorem 2.
Consider system (25), in closed loop with(31)u=-1bKTx^-yr+x^n+1-yrn,where x^=[x^1,…,x^n]T and x^n+1 are computed via (26) and (27), with K∈Rn selected, such that(32)Hn=An-EnKTis a Hurwitz matrix, and yrT=yr,y˙r,…,yrn-1 is the reference vector, with yrj≤Γj, for j={0,…,n+1}, provided that
∇xf and fd=∂f/∂d are bounded within the domain of interest, with d,d˙∈L∞, and f(x,d) being continuous;
λi{Hn}≠λj{w0A_}, for all i and j, where i={1,2,…,n}, and j={1,2,…,n+1};
w0>max1,HnA_-1. Then, the output feedback dynamic, defined by (26) and (31), is practically stable. Moreover, x~e=xe-x^e is bounded with the ultimate bound given by O(1/w0), where xe=x,fT. Finally, if Δt=0, where,(33)Δt=fdd˙+∂f∂xnyrn+∇xfAnyr,
then x~e exponentially converges to zero. The proof of this theorem can be found in Appendix.
4. A GPI Observer-Based Active Disturbance Rejection Controller
From the linearized control model (15), we use the following simplified perturbed model for the unknown nonlinear system (11), defined as(34)Fδ4=-uδN+ηt,where ηt acts in place of all the following: the state dependent expressions, all the higher-order terms (h.o.t.) neglected by the linearization, possibly of nonmodeled dynamics, and, finally, the external unknown disturbances affecting the system. All of these uncertain terms are grouped together as a single time-varying function, which in our case has the following form:(35)ηt=1NFδ+Fδ2-μNFδ2+h.o.t.Assuming that η(t) satisfies Assumption (A1) found in Theorem 2, it is possible to solve the output feedback trajectory control problem of the uncertain system (34). That is, we can propose an admissible reference trajectory, F∗=F∗(t), in order to carry out some controlled maneuver tasks. Evidently, this trajectory has to be chosen, such that it is very close to the rest position (x1δ=0,x2δ=0). Recall that the incremental linear control model (34) is valid, as long as the trajectories of the original system are in a vicinity close to the origin. In order to solve this problem, we defined the trajectory tracking control errors as(36)ei=Fδi-1-F∗i-1t,withi=1,2,3,4.Based on this definition, we can see that the errors evolve to(37)e˙1=e2;e˙2=e3;e˙3=e4;e˙4=-uδN+ηt.Since variables x1δ and x2δ are available, then, from Remark 1, the flat outputs Fδ=x1δ+x2δ and Fδ(2)=-x2δ are also available, implying that the tracking errors e1=Fδ-F∗ and e3=Fδ(2)-F∗(2) can be measured simultaneously. Hence, a set of decoupled high-gain extended linear Luenberger observers for the estimation of the unknown variables e2,e4,η is proposed as follows:(38)e^˙1=e^2+2w0e1-e^1;e^˙2=e3+w02e1-e^1;e^˙3=e^4+3w0e3-e^3;e^˙4=-uδN+η^+3w02e3-e^3;η^˙=w03e3-e^3,where w0>1. Notice that the above linear observers are in agreement with the EOS proposed in (26) and (28). Having proposed the high-gain observers for the estimation of the unknown variables, we proceed to design the linear controller based on the ADRC approach. That is, the input uδ will be synthesized with an active disturbance canceling strategy for the uncertain η, in terms of its estimated value η^, and the use of variables e^2 and e^4 instead of the actual variables e2 and e4, respectively. Therefore, according to (18) and (31), we can propose uδ as(39)uδ=Nk4e^4+k3e3+k2e^2+k1e1-F∗4+η^,where KT=k1,k2,k3,k4 must be selected, such that the matrix(40)H4=010000100001-k1-k2-k3-k4is Hurwitz (see definition in (32)).
Remark 3.
Matrix A_ is defined by(41)A_=-310-301-100(see (30)), with λjA_=-1 and j=1,2,3. On the other hand, the first and second equations of (38) are neglected, because the observation error e~1=e1-e^1 satisfies the following equation:(42)e~¨1+2w0e~˙1+e~1=0.This implies that e~2 exponentially converges to zero.
From the discussion above, the following can be concluded.
Fact 1.
η(t), defined in (35), satisfies the conditions in Assumption (A1).
Fact 2.
Selecting, for instance, KT, such that λjH4=-1 for i=1,2,3,4 and w0≫1, we can always fulfill the conditions in Assumptions (A2) and (A3).
From Facts 1 and 2 and according to Theorem 2, we can assure that the system locally asymptotically and exponentially converges to zero, if F∗=0. Besides, if F∗t≠0 and its time derivatives up to order 4 are bounded, then we can always assure that the tracking error is ultimately bounded, where the confined region of the tracking error can be as small as desired.
Remark 4.
In order to avoid the large initial peaking phenomena found in the response of observer variables, we suggest using a clutch function to smooth these transient peaking responses in all the observer variables used in the controller. The “clutch” is defined as a time function smoothly increasing from 0 to 1, during a small time interval [0,α). We fix the smoothing function as(43)sft=1fort>α;sinrπt2αfor0≤t<α,where r is a suitably large positive even integer (see [19]). Therefore, the “smoothing” of the observer variables may be implemented as(44)z^f=z^sft,where z=e2,e4,η.
Remark 5.
As already mentioned, an advantage of combining the flatness and the ADRC methods is that it allows online estimates and cancels the undesirable effects of the higher-order nonlinearities discarded by the linearization approximation. Because our result is based on this advantage, it is important to provide some arguments to validate it. The ADRC is based on restriction to (1) flat systems and their corresponding input-to-flat outputs dynamics (multivariable or monovariable) and (2) a brute force exact linearization of the underlying input output dynamics (monovariable) or the set of statically or dynamically decoupled nonlinear dynamics, written in Isidori’s canonical form. This is achieved by online estimating and online feedback cancelling of absolutely everything that perturbs the forced linear dynamics from its desired nominal behaviour (trajectory tracking or stabilization). This uncertainty cancelling includes poorly known additive expressions containing state-dependent nonlinearities and the unpredictable effects of unmodeled dynamics and of external unknown disturbances. Very many examples and applications regarding the effectiveness of this technique as well as the theoretical results backing the methodology actually constitute sufficient proof of the assertion made in Abstract of the paper. References [13–19] are vivid proof of the range of applications and effectiveness of the method, not to mention the underlying simplicity. The following articles also contain sufficient results and applications which validate the technique in a generous manner [32–36]. Reference [37] offers complete survey up to 2014.
Numerical Simulation. To test the performance of the proposed ADRC in conjunction with the flatness controller, we carried out two numerical simulations, with the following setup:(45)m=0.1Kg;M=0.5Kg;Lr=0.5m;L=0.5m;Ib=0.13Kgm2,and the closed loop characteristic polynomial chosen as(46)ps=s2+2ζωns+ωn22,with ζ=0.707 and ωn=1.2. Notice that Ib was computed via the formula Ib=M(Lr2+L2/12). The aim of the first simulation was to solve the stabilization of the BRT system at the rest equilibrium position. The other simulation solved the output feedback trajectory control. To this end, we use the control equation, (39), with the decoupled corresponding observer (38). The observer parameter was fixed as w0=10, and the clutch parameters were set as α=0.1[s] and r=3.
The Output Feedback Stabilizing Controller. The closed loop stabilization to the GPI observer based on the ADRC of the system (8) was carried out. To this end, the reference trajectory, F∗, for the flat output, Fδ(t), was set to zero, as required by the stabilization task. The pair of initial conditions, for equal numbers of executions, was selected as(47)p1=x1=0,x˙1=0.5m/s,x2=-0.2rad,x˙2=-0.1rad/s,p2=x1=0.1m,x˙1=0.3m/s,x2=0,x˙2=-0.1.Figure 2 shows the closed loop responses of the ADRC for the stabilization of system (8) at the rest position for each initial condition. From this figure we can see that, before 4 seconds elapse, the position variables and the control action go to zero.
Closed loop response of the ADRC for two different initial conditions: p1 and p2.
The Output Feedback Trajectory Tracking Problem. A more challenging control task was selected for the second experiment, where the sine signal F∗t=0.4sin(t/2) was selected as the reference trajectory. To this end, we fixed the initial condition at the origin and used the same setup as in the previous simulation.
Figure 3 shows a window from 5 to 30 seconds of the evolution of the tracking errors ei=Fδ(i-1)-F∗(i-1)(t), with i=1,2,3,4. The tracking errors e1, e2, e3, and e4, are, respectively, in the ranges of 10-6, 10-5, 10-4, and 10-3. From these results we can claim that the system is able to effectively track the proposed reference after 4 seconds has elapsed. Additionally, the closed loop response showed a considerably large domain of attraction.
A window time, from 5 to 30 seconds, of the tracking trajectory error when the reference is a sinusoidal signal.
Finally, Figure 4 shows the phase space of the position and angular variables. There, we can see that those variables have a limit cycle in their corresponding phase spaces.
Phase space of the position and angular variables.
5. Conclusions
The differential flatness approach, in conjunction with ADRC, allows systematic solutions to a number of interesting nonlinear control problems. In this instance, we have exploited the local flatness property for the efficient stabilization and tracking of the underactuated BRT. It is important to emphasize that this system is not feedback-linearizable and that its relative degree is not well defined. Moreover, as far as we know, the stabilization of the BRT remained unsolved by using either the shaping energy approach or the IDA-PBC method. Nevertheless, this problem can be partially solved using linear control theory. The fact that the tangent linearization of this system, around the unstable equilibrium point, is locally controllable implies that the system is also locally flat. This allows us to use the robust ADRC in the efficient online estimation of the locally neglected nonlinearities and their active feedback cancelation. The solution, which is quite robust with respect to unmodeled disturbances and neglected nonlinearities, is, in fact, a linear controller with an online compensator. It is based on a set of decoupled linear extended observers and a single linear output feedback controller, with disturbance cancelation features. The proposed controller guarantees locally exponentially asymptotic stability for the stabilization problem and practical and local stability in the solution of the tracking error. To assess the effectiveness of the proposed methodology, numerical simulations were carried out. From the simulation results, we can demonstrate that the proposed controller behaves remarkably well, having an acceptable domain of attraction.
AppendixProof of Theorem 2
The following proof is based on the previous works [38–40]. First of all, we define the actual extended state xe=x1,x2,…,xn,fT∈Rn+1. Hence, system (25) can be rewritten as(A.1)x˙e=An+1xe+Bn+1u+En+1η;ye=Cn+1Txe,with η:=f˙. Evidently, system (26) is a Luenberger observer for the system above. Defining the observation error vector as x~e=xe-x^e, we can easily obtain, from (A.1) and (26), the following equation:(A.2)x~˙e=An+1-lCn+1Tx~e+En+1η,which is equivalent to(A.3)ξ˙=w0A_ξ+w0-n+1En+1η,where we use the following transformation:(A.4)x~e≔diagw0,w02,…,w0n+1ξ=Wn+1ξ.Using the definition of η=f˙ and (33), we obtain(A.5)ξ˙=w0A_ξ+w0-n+1∂f∂xnEn+1∇xfAn-KTe+K¯TWn+1ξ+w0-n+1En+1Δt,where the definition of A_ was defined in (30).
Stability Analysis of the Trajectory Tracking Error e=x-yr. Substituting the value of (31) in (25) leads to(A.6)x˙=Anx+En-KTx^-yr+w0n+1ξn+1+yrn.Expressing the last equation in terms of the tracking error, we obtain(A.7)e˙=Hne+EnK¯TWn+1ξ.where K¯T=[KT1] and Hn=An-EnKT. Now, we analyse the convergence of (A.7) and (A.5). Then, both equations are easily rewritten as(A.8)e˙ξ˙=HnEnK¯TWn+10w0A_eξ+0n×1En+1w0n+1Φ·,where (A.9)Φx,e,ξ,t=F1xe+F2xξ+Δt;F1·=∂f∂xn∇xfAn-KT;F2·=∂f∂xnK¯TWn+1.Hence, Φ(·) is continuous and bounded, according to Assumption (A1). Introduce the following change of coordinates:(A.10)ζ=e+Nξ,where N is a constant matrix, which is a solution of the following Sylvester equation:(A.11)EnK¯TWn+1=HnN-w0NA_.Assumption (A2) assures existence and uniqueness of N [41]. Consequently, (A.10) can be rewritten as (A.12)ζ˙ξ˙=Hn00w0A_ζξ+0n×1En+1w0n+1Φ·.Because Hn is selected to be Hurwitz, we conclude that ζ converges exponentially to zero. Finally, we need to prove boundedness and convergence of the system:(A.13)ξ˙=w0A_ξ+w0-n+1En+1Φx,0,ξ,t.To this end, we need to introduce the following auxiliary lemma to estimate the norm of the matrix N, without the need to compute it.
Lemma A.1.
Under Assumptions (A2) and (A3), the matrix N of (A.11) is bounded by(A.14)N≤w0nkN,kN≔w0K¯A_-1w0-HnA_-1>0.
Besides, for K¯ that is sufficiently large, then lim∥K∥→∞N≤w0n+1; because the proof of this fact is not difficult, it is omitted.
Resuming the stability proof, because A_ is Hurwitz, there exists P=PT>0, such that PA_+A_TP=-2In+1. Defining the Lyapunov function V=1/2ξTPξ, it is easy to see that its time derivative, along the trajectories of (A.13), can be upper bounded by(A.15)V˙≤-w0ξ2+λmaxPkF1kNw0+kF2w0n+1ξ2+λmaxPkΔξw0n+1,where Δ≤kΔ, F1≤kF1, and F2≤kF2. From the above, we need to analyse two cases:
(kΔ=0) Selecting w0>1, such that(A.16)λmaxPkF1kNw0+kF2w0n+1=kξ<w0,then V˙<0, concluding that ξ exponentially converges to zero. Notice that (A.16) can be always satisfied selecting the free observer parameter w0≫1.
(kΔ≠0) If (A.16) holds, then from (A.15) the ultimate boundedness of ξ with ultimate bound is(A.17)ξ≤λmaxPkΔw0n+2w0-kξ.Notice that w0-kξ≫1, by selecting w0 that is sufficiently large. Hence, ξ~O(1/w0(n+2)). On the other hand, from (A.10), since ζ→0, the error e~O(1/w02). That is, the conditions of Theorem 2 are satisfied.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This research was supported by the Centro de Investigación en Computación of the Instituto Politécnico Nacional (CIC-IPN) and by the Secretaría de Investigación y Posgrado of the Instituto Politécnico Nacional (SIP-IPN), under Research Grants 20160268 and 20161637. This research was done while Dr. Carlos Aguilar-Ibanez was on sabbatical leave from the Departamento de Mecatrónica del CINVESTAV.
FantoniI.LozanoR.2002SpringerSarrasI.AcostaJ. Á.OrtegaR.MahindrakarA. D.Constructive immersion and invariance stabilization for a class of underactuated mechanical systems20134951442144810.1016/j.automatica.2013.01.059MR30440262-s2.0-84876665463González-VázquezS.Moreno-ValenzuelaJ.Motion control of a quadrotor aircraft via singular perturbations201310, article 368Aguilar-AvelarC.Moreno-ValenzuelaJ.A composite controller for trajectory tracking applied to the Furuta pendulum20155728629410.1016/j.isatra.2015.02.009SpongM. W.Underactuated mechanical systems1998Berlin, GermanySpringer135150Olfati-SaberR.2000Massachusetts Institute of TechnologyGaoB.ChenH.ZhangX.Control design for underactuated mechanical systems: a survey2006105541Sira-RamirezH.On the control of the ‘ball and beam’ system: a trajectory planning approachProceedings of the IEEE Conference on Decision and Control2000Sydney, AustraliaRouchonP.FliessM.LévineJ.MartinP.Flatness, motion planning and trailer systemsProceedings of the 32nd IEEE Conference on Decision and ControlDecember 1993San Antonio, Tex, USAIEEE270027052-s2.0-0027758507FliessM.LévineJ.MartinP.RouchonP.Design of trajectory stabilizing feedback for driftless at systemsProceedings of the Third ECC1995Rome, Italy18821887LévineJ.2009Berlin, GermanySpringer Science & Business Media10.1007/978-3-642-00839-9MR2723930FliessM.LévineJ.MartinP.RouchonP.Flatness and defect of non-linear systems: introductory theory and examples19956161327136110.1080/00207179508921959MR16135572-s2.0-0029327262Ramírez-NeriaM.Sira-RamírezH.Garrido-MoctezumaR.Luviano-JuárezA.Linear active disturbance rejection control of underactuated systems: the case of the Furuta pendulum201453492092810.1016/j.isatra.2013.09.0232-s2.0-84905121508Sira-RamírezH.López-UribeC.Velasco-VillaM.Linear observer-based active disturbance rejection control of the omnidirectional mobile robot2013151516310.1002/asjc.523MR30157582-s2.0-84872296738HanJ.From PID to active disturbance rejection control200956390090610.1109/tie.2008.20116212-s2.0-63149185343GaoZ.HuangY.HanJ.An alternative paradigm for control system design5Proceedings of the 40th IEEE Conference on Decision and Control (CDC '01)December 2001Orlando, Fla, USA457845852-s2.0-0035705733GaoZ.HuS.JiangF.A novel motion control design approach based on active disturbance rejection5Proceedings of the 40th IEEE Conference on Decision and Control (CDC '01)December 2001Orlando, Fla, USA4877488210.1109/.2001.9809802-s2.0-0035707031LiR.LiT.ZhengQ.LiQ.Ship tracking control based on linear active disturbance rejection controlProceedings of the 2012 3rd International Conference on Intelligent Control and Information Processing (ICICIP '12)July 2012Dalian, ChinaIEEE20120510.1109/icicip.2012.63914532-s2.0-84872305023MoralesR.Sira-RamírezH.SomolinosJ. A.Linear active disturbance rejection control of the hovercraft vessel model20159610010810.1016/j.oceaneng.2014.12.0312-s2.0-84921038617Coral-EnriquezH.Cortés-RomeroJ.RamosG. A.Robust active disturbance rejection control approach to maximize energy capture in variable-speed wind turbines201320131239674010.1155/2013/396740MR3080125Cortés-RomeroJ.Rojas-CubidesH.Coral-EnriquezH.Sira-RamírezH.Luviano-JuárezA.Active disturbance rejection approach for robust fault-tolerant control via observer assisted sliding mode control201320131260952310.1155/2013/609523MR3142657LiR.LiT.BuR.ZhengQ.ChenC. L. P.Active disturbance rejection with sliding mode control based course and path following for underactuated ships20132013974371610.1155/2013/7437162-s2.0-84890028348KhalilH. K.High-gain observers in nonlinear feedback controlProceedings of the International Conference on Control, Automation and Systems (ICCAS '08)October 2008Seoul, South KoreaIEEE475710.1109/iccas.2008.46947052-s2.0-58149086895PugaS.BonillaM.MalabreM.LozanoR.Singularly perturbed implicit control law for linear time varying SISO systems201424101530154910.1002/rnc.2946MR32080452-s2.0-84900023996PugaS.BonillaM.MoogC. H.MalabreM.LozanoR.Singularly perturbed feedback linearization for SISO nonlinear systems with measurement of the stateProceedings of the European Control Conference (ECC '15)July 2015Linz, Austria3250325510.1109/ecc.2015.7331035RubioJ. D. J.OchoaG.BalcazarR.PachecoJ.Uniform stable observer for the disturbance estimation in two renewable energy systems20155815516410.1016/j.isatra.2015.05.0032-s2.0-84943584179BlochA. M.LeonardN. E.MarsdenJ. E.Matching and stabilization by the method of controlled lagrangians2Proceedings of the 37th IEEE Conference on Decision and Control1998Piscataway, NJ, USAIEEE14461451OrtegaR.SpongM. W.Gómez-EsternF.BlankensteinG.Stabilization of a class of underactuated mechanical systems via interconnection and damping assignment20024781218123310.1109/tac.2002.800770MR19174332-s2.0-0036687767OrtegaR.Garcia-CansecoE.Interconnection and damping assignment passivity-based control: a survey200410543245010.3166/ejc.10.432-450MR21343242-s2.0-27844531885GaoQ.HouY.LiK.SunZ.WangC.HouR.Neural network based active disturbance rejection control of a novel electrohydraulic servo system for simultaneously balancing and positioning by isoactuation configuration201620169492109510.1155/2016/4921095Galvan-ColmenaresS.Moreno-Armendáriz,M. A.de Jesús RubioJ.Ortíz-RodriguezF.YuW.Aguilar-IbáñezC. F.Dual PD control regulation with nonlinear compensation for a ball and plate system201420141089420910.1155/2014/894209MR31985702-s2.0-84899989200ZhengQ.GaoL. Q.GaoZ.On validation of extended state observer through analysis and experimentation2012134202450510.1115/1.40053642-s2.0-84855913911GuoB.-Z.ZhaoZ.-L.Weak convergence of nonlinear high-gain tracking differentiator20135841074108010.1109/tac.2012.2218153MR30665752-s2.0-84894532856GuoB.-Z.ZhaoZ.-L.On the convergence of an extended state observer for nonlinear systems with uncertainty201160642043010.1016/j.sysconle.2011.03.0082-s2.0-79955794731GuoB.-Z.ZhaoZ.-L.On convergence of the nonlinear active disturbance rejection control for MIMO systems20135121727175710.1137/110856824MR30456482-s2.0-84879627920ChenW.-H.YangJ.GuoL.LiS.Disturbance-observer-based control and related methods—an overview20166321083109510.1109/tie.2015.2478397MadońskiR.HermanP.Survey on methods of increasing the efficiency of extended state disturbance observers201556182710.1016/j.isatra.2014.11.0082-s2.0-84929271631AtassiA. N.KhalilH. K.A separation principle for the stabilization of a class of nonlinear systems19994491672168710.1109/9.788534MR17098632-s2.0-0032594209ZhouW.ShaoS.GaoZ.A stability study of the active disturbance rejection control problem by a singular perturbation approach200939–12491508MR25011342-s2.0-65549098107KhalilH. K.PralyL.High-gain observers in nonlinear feedback control2014246993101510.1002/rnc.3051MR31811732-s2.0-84896543901BhatiaR.RosenthalP.How and why to solve the operator equation AX-XB=Y199729112110.1112/s0024609396001828MR1416400