Quadcopter UAVs Extended States/Disturbance Observer-Based Nonlinear Robust Backstepping Control

A trajectory tracking control for quadcopter unmanned aerial vehicle (UAV) based on a nonlinear robust backstepping algorithm and extended state/disturbance observer (ESDO) is presented in this paper. To obtain robust attitude stabilization and superior performance of three-dimension position tracking control, the construction of the proposed algorithm can be separated into three parts. First, a mathematical model of UAV negatively influenced by exogenous disturbances is established. Following, an extended state/disturbance observer using a general second-order model is designed to approximate undesirable influences of perturbations on the UAVs dynamics. Finally, a nonlinear robust controller is constructed by an integration of the nominal backstepping technique with ESDO to enhance the performance of attitude and position control mode. Robust stability of the closed-loop disturbed system is obtained and guaranteed through the Lyapunov theorem without precise knowledge of the upper bound condition of perturbations. Lastly, a numerical simulation is carried out and compared with other previous controllers to demonstrate the great advantage and effectiveness of the proposed control method.


Introduction
The unmanned aerial vehicles also known as drones are a great achievement of science and technology. They are widely used in various fields of interest from academic research to industrial applications. Thus, there are an increasing number of scientific studies to expand the range of activities and enhance the operating capacity of these vehicles to reach their limit. The majority of this research has significantly been achieved in several major areas, such as a reliable controller construction [1], innovative airframe optimization design [2], efficient communication [3], and energy management for UAVs [4]. Nevertheless, there always exist many further challenges to construct an efficient flight controller to guarantee rigorous stability and rapid adaptation of the vehicle's behavior in various flight environments due to different causes, such as: (i) The quadcopter UAV suffers from the intrinsic characteristic of an under-actuated system, meaning that the vehicle has less actuator control input than the number of degrees of freedom (DOF), which must be controlled. In addition, the various sources of parametric uncertainties and actuator failure are also considered as an essential characteristic needing to be handled; (ii) UAVs are also known as complex systems with an extremely nonlinear characteristic, open-loop instability, and strong coupled system with multiple input and multiple output (MIMO) variables in which the interaction between the rotational and translational dynamic model is highly complicated; (iii) The operational principle of multi-copter UAVs is based on the aerodynamic force generated by the spinning of multi-rotors and propellers. Thus, the vehicle is extremely sensitive to any interference of external disturbances, particularly gusts of wind.
To overcome the control challenge, a variety of works to improve the control performance of quadcopter UAVs have been conducted in recent decades. Many traditional linear control techniques have been introduced by using the Proportional Integral Derivative (PID) control law [5] and the Linear Quadratic Regulator (LQR) algorithm [6]. Although the PID and LQR can manage the hovering behavior of the vehicle and handle some simple actions, strict stability of the control system cannot be guaranteed in the presence of uncertainties and/or external perturbations. Furthermore, to design these linear controllers, a linearization process needs to be carried out surrounding single or multiple operating points. Thus, these methods still face the disadvantage of breaking down performance once the vehicle diverges from the considered points.
To attenuate the undesirable effects of nonlinearity parts and uncertainties/disturbances generated from various sources, plentiful control techniques have been introduced, such as a nonlinear control method presented in [7,8], whereby the controller guarantees a null tracking error and robust stability of the vehicle in the appearance of perturbations. In [9], an adaptive controller for trajectory tracking of UAVs is introduced by using an accurate parameterization algorithm. In [10], a nonlinear controller for the attitude and position tracking control of quadcopter UAVs is introduced by integration between a commandfiltered backstepping algorithm and the gain scheduled method. In [11,12], a combination of model predictive algorithm and robust control (H ∞ ) is conducted for trajectory tracking and attitude stabilization of a quadcopter UAV to obtain a null steady-state error under the effect of disturbances. However, in all these controller methods, the nominal control performance needs to be sacrificed to handle the stabilization of the control system in the presence of multiple uncertainties and/or disturbance sources.
With a simple procedure of construction and strict robustness to counteract uncertainties/disturbances, the sliding mode control (SMC) and extended SMC have been largely developed in the UAV control system, such as second-order sliding mode control [13][14][15][16], backstepping and integral SMC [17,18], and adaptive SMC [19,20]. It is well-known as a robust control in the SMC community, the nonsingular terminal SMC (NSTSMC) is an interesting technique [21], because the algorithm ensures both finite-time stability of sliding surface and the convergence of system states to desired trajectories. Motivated by the nonsingular terminal SMC and a nonlinear disturbance observer-based control [22,23], a composite control algorithm [24] is presented by an integration of the single hidden layer feedforward network and an NSTSMC to ensure the zero convergence of tracking error in a finite-time. In [25], another research based on a finite-time disturbance observer and NSTSMC is also developed to enhance the tracking control. However, these controller algorithms are only suitable for simple desired trajectories because several assumptions of these methods may not be practical in some cases. In [26], the concept of a virtual control variable is used to restructure the underactuated UAV system to obtain a fully-actuated model. Afterward, an advanced SMC is designed to enhance the control performance of the attitude and position subsystem. However, to reduce an inherent chattering effect in the SMC technique, the high gain parameters of the switching law must be avoided. Thus, this controller cannot strictly guarantee both precise tracking control and a strong anti-chattering effect.
The backstepping control technique and its extended methods are known as efficient recursive nonlinear control algorithms for a particular class of engineering system in which the controller is constructed from plentiful subsystem models. The process of designing a backstepping algorithm is to start from the origin of a complex system based on the form of state feedback, and substitute to the next step of the procedure. In this process, auxiliary variables and virtual control inputs are proposed to achieve the stable performance of each subsystem gradually. The process of backing out and substituting is terminated once an ultimate law is obtained [27][28][29]. Due to the advantage of the recursive structure, the backstepping method is renowned for its efficiency in the group of robust control techniques for complex engineering models, especially for unmanned aerial vehicles. Although there are many previous works that have solved the tracking control problem in both attitude stabilization and position control [30][31][32], the control performance is also ruined in the various working situations because it is not easy to estimate the variation in the uncertainties, exogenous perturbations, and actuator faults. To deal with this problem, the other extended backstepping controllers using disturbance observer and fault-estimation are employed [33][34][35][36]. In [37], an extended state observer is designed for online approximation in both the unknown velocity states and the perturbations affecting subsystem dynamics. Following, the dynamic surface controller is presented by using the estimated values to ensure the convergence of the tracking error to zero. However, in this research, the disturbance observer is constructed by using UAV position states computed by an integration of the velocity feedback. Obviously, these position states are not accurate in long-term operation because the accumulated error is always progressive during the integration process. Furthermore, the upper bound knowledge of uncertainties and/or exogenous perturbations is required for the observer design step. In [38], an active disturbance rejection control (ADRC) is introduced by using an auxiliary variable and cascade control method to ameliorate the robustness and anti-perturbation of the UAV control system. However, the drawback of ADRC is to cover some indeterminate coefficients consisting of tracking differentiators. Furthermore, the conflict between these parameters may occur in the controller tuning process.
The inspiration of this study is to deal with the aforementioned drawbacks of the existing research. A nonlinear robust backstepping control based on an extended state/disturbance observer (RBCESDO) for attitude stabilization and trajectory tracking control is presented to improve the performance of quadcopter UAVs under the various influences of uncertainties and/or external disturbances. In this method, the auxiliary variables, virtual controls, and ultimate control law are recursively constructed by using the approximated values of the system states and disturbances given by the ESDO. The contribution of this study can be stated as follows: (i) The proposed algorithm overcomes the drawbacks of previous methods in the requirement of full state measurement. The ESDO is able to estimate the velocity state of the vehicle once this parameter cannot be directly measured. Thus, the implementation cost for data acquisition may be reduced and the influence of high measurement noise generated from the velocity sensor is also alleviated. (ii) The unmeasured velocity states and lumped perturbations are estimated by the presented ESDO integrating with advantages of the recursive structure of the backstepping technique, the convergence of tracking errors is always guaranteed. (iii) The numerical simulation is fully executed in both attitude and position control mode; these performance results are compared with other existing control methods to confirm the strict stability and efficiency of the presented control scheme in both the convergence error and anti-disturbance capacity. (iv) Finally, unlike the existing methods, the upper bounds of the uncertainties and/or external disturbances are not demanded during the steps of designing the proposed control scheme.
The remainder of this article is organized as follows. The mathematical model of a quadcopter UAV and problem description is presented in Section 2. In Section 3, the steps of constructing a nonlinear robust backstepping controller based on ESDO and stability analysis are derived. A numerical simulation conducted and compared with other existing methods is given in Section 4. Final statements of this study are given in Section 5 through some conclusions.

Mathematical Model and Problem Description
In this section, the mathematical model of a quadcopter UAV is constructed in the presence of external disturbances: There are some assumptions are given as follows: Hypothesis 1. The structure of a quadcopter UAV is considered as a rigid body.

Hypothesis 2.
The coefficients of four rotors and parameters of propellers are identical.
The quadcopter UAV can be considered as a 6-DOFs rigid body, which is an idealized expression of invariable volume. Thus, the behavior of the quadcopter (i.e., rotational motion and translational motion) can be described through the movement of a particle positioned in the UAV center of mass. Obviously, a dynamic model of the vehicle entirely depends on aerodynamic force and torque generated by the speed of propellers mounted on four rotors. Thus, it is necessary to consider that the properties of four rotors and propellers are identical to facilitate the formulation of a mathematical model. Through two essential coordinate frames, i.e., Earth frame {E}, and body frame {B} (as shown in Figure 1), the dynamical model of a quadcopter UAV can be described from two subsystems of translational motion (Q) and rotational motion (Θ) by using the Newton-Euler formula as follows: where Q = [x, y, z] T ∈ R 3 is the absolute position of the quadrotor UAV with respect to the earth frame {E}; Θ = [φ, θ, ψ] T ∈ R 3 , φ, θ ∈ (−π/2, π/2), represent a vector of Euler angles (i.e., the roll, pitch, and yaw); G = [0, 0, mg] T ∈ R 3 with m ∈ R + and g = 9.81 m/s 2 denote the total mass of the vehicle and gravity acceleration constant, respectively; ∈ R 3 is angular velocity; I = diag I x , I y , I z ∈ R 3×3 denotes the matrix of the inertial moment along x, y, and z axes; F aero = [0, 0, f z ] T ∈ R 3 and M aero = τ φ , τ θ , τ ψ T ∈ R 3 are the vectors of aerodynamic force and torque generated by multi-rotors and propellers of the vehicle, in which f z [N] is the thrust force; τ φ , τ θ , τ ψ [N.m] represent the roll, pitch, and yaw torques, respectively. These values can be computed by: . . , 4 represent the aerodynamic force and torque produced by the rotor and propeller i-th; b, d ∈ R + are thrust and drag coefficients, respectively; Ω i is the angular velocity of motor i-th; l is the arm length of the vehicle frame.
and ( ) is a unit vector along the z-axis.
In addition, in order to design an efficient control law for a quadcopter UAV, it cannot ignore the undesirable influence of disturbances and/or uncertainties on the dynamical system. Hence, are considered as the effects of perturbations/uncertainties on the translational and rotational subsystems, respectively. In addition, T E B (Θ) and R E B (Θ) represent the translational and rotational transformation matrices from the body frame {B} to earth frame {E}, respectively. These values can be computed by and where S (•) , C (•) , and T (•) represent the notations of sin(•), cos(•), and tan(•), respectively. It is proposed to consider that the variations in the three Euler angles are unprogressive maneuvers, meaning that the rotational behaviors of the roll, pitch, and yaw angles are surrounding the origin. Thus, R E B (Θ) is closely equivalent to the identity matrix, i.e., R E B (Θ) = I 3×3 ; M gyros denotes the gyroscopic torque generated by the total residual angular velocity of four rotors, i.e., Ω = . . , 4, effecting on a rigid body of the vehicle, its value is given by: whereê 3 = [0, 0, 1] T is a unit vector along the z-axis.
In addition, in order to design an efficient control law for a quadcopter UAV, it cannot ignore the undesirable influence of disturbances and/or uncertainties on the dynamical system. Hence, F dis = ξ x , ξ y , ξ z T and M dis = ξ φ , ξ θ , ξ ψ T are considered as the effects of perturbations/uncertainties on the translational and rotational subsystems, respectively. From Equations (1)-(5), the mathematical model of a quadcopter UAV can be clearly described by: In order to construct a robust backstepping controller-based ESDO for improving the tracking performance of a quadcopter UAV, an equivalent MIMO control system under the effects of disturbances and/or uncertainties can be described by a state-space model as follows, ..
where Z ∈ R 6 , u ∈ R 6 , and ξ ∈ R 6 represent the vectors of system states, control input, and lumped perturbations and/or uncertainties, respectively. Its detailed descriptions can be given as follows, T and the terms of g(Z, t) ∈ R 6 and Υ ∈ R 6×6 are given as follows: The objective of controller design is to derive a robust backstepping control law guaranteeing a good convergence of the output state, Z, to the desired trajectory Z d ∈ R 6 , in the appearance of uncertain parameters and/or external perturbations. It can be seen that the quadcopter UAV is an underactuated mechanical system because the vehicle has four control inputs f z , τ φ , τ θ , τ ψ but there are six degrees of freedom to be controlled (x, y, z, φ, θ, ψ). Hence, in order to design a robust controller for a full dynamic system, the attitude control for Euler angles (φ, θ, ψ), and position control for the three dimensions (x, y, z) are simultaneously constructed by using two separated subsystems (i.e., rotational and translational dynamics). The attitude subsystem comprising roll, pitch, and yaw rotation is controlled by three torques τ φ , τ θ , τ ψ (i.e., u 4 , u 5 , u 6 ); the vertical position (z) is controlled by a control input u z obtained from Equations (6)-(9) as follows: Therefore, the aerodynamic force f z can be derived from the controller u z as follows, The remaining states of the control system (i.e., horizontal positions x and y), are indirectly controlled through desired attitude states (φ d , θ d , ψ d ) and vertical control force f z achieved from Equations (6)-(9) as follows: where u x , u y , u z (i.e., u 1 , u 2 , u 3 ) are designed controllers for translational behaviors along x, y, and z axes, respectively; the heading angle, ψ d , is known parameters in advance.

Robust Backstepping Control-Based ESDO
A robust backstepping controller based on an extended state/disturbance observer for a quadcopter UAV is presented in this section in the following steps: Step 1: ESDO is designed by using the second-order system to observe the unmeasured system states consisting of the effectiveness of disturbances/uncertainties on the whole control system.
Step 2: A full robust backstepping controller is derived by an integration of the approximated values from the proposed ESDO to reimburse the negative effects of uncertainties and/or external disturbances on a quadcopter UAV. In addition, a stability analysis is also fulfilled to demonstrate a good convergence of the output states to desired trajectories. During this procedure, the precise knowledge of the upper bounds of perturbations is not required.

Extended State/Disturbance Observer (ESDO)
This subsection investigates a method to observe the unmeasured states of the vehicle system including the lumped uncertainties and/or external disturbances. To derive an ESDO, a MIMO nonlinear system is obtained from Equation (7) as follows, It can be rewritten as the shorter form: .
Z, ξ T ∈ R 12 denotes a new vector of output state and disturbance in- The time-varying extended state/disturbance observer for a quadcopter UAV is proposed as follows: .
The stability of the proposed ESDO can be achieved through the observer error model obtained by Equations (15) and (16) as follows: .
Obviously, it is always feasible to select constant values of α 1i , β 1i ∈ R + , i = 1, . . . , 6 such that the matrix Φ 1 meets the Hurwitz stability theorem meaning that the eigenvalues of matrix Φ 1 , (i.e., λ = −λ j < 0, j = 1, 2, . . . , 12), are arbitrarily located in the left half side of the complex plane (LHP). Thus, Theorem 1. Let assumes that there always exist positive constants of δ, σ, ϕ ∈ R + in such a way that the following expression always satisfies: where λ max {Φ 1 } = −ϕ < 0 is a maximum eigenvalue of the matrix Φ 1 . The observation error, χ(t) given by Equation (17), is strongly convergent to a small ball containing the origin zero if the presented ESDO given in Equation (16) is selected to observe the system states and lumped perturbations of a quadcopter UAV.
Proof of Theorem 1. The solution, χ(t), of the differential Equation (18) can be easily achieved by [39] as follows: The inequality, , is always satisfied with the assumption given in Equation (20). In addition, the term of σδ E 1 e −ϕ(t−τ) is a positive scalar: thus, the bound of approximate error, χ(t), can be obtained as follows: Due to e −ϕ(t−t 0 ) → 0, t 0 ≤ ∀t → ∞ , thus χ(t) is bounded by: The proof of Theorem 1 is completed. (23), we can conclude that the observation error, χ(t) , is ultimately bounded by a constant µ, if the observation gains α 1i , β 1i ∈ R + , i = 1, . . . , 6 are properly chosen such that matrix Φ 1 satisfies the Hurwitz criteria. Therefore, .

Z(t) and ξ(t)
are also convergent to a small region containing the origin zero, i.e., .

Robust Backstepping Controller Design
In this subsection, a robust backstepping controller is constructed by a combination of the approximated value from the proposed ESDO above and backstepping control technique to reimburse the negative effects of uncertainties and/or external disturbances on a quadcopter UAV. The presented controller guarantees the convergence of the output state, Z = [x, y, z, φ, θ, ψ] T , to the desired trajectory Z d = [x d , y d , z d , φ d , θ d , ψ d ] T . A general scheme of the proposed controller for a quadcopter UAV is shown in Figure 2. Let z 1 = [ z 11 , z 1i , . . . , z 16 ] T ∈ R 6 , i = 1, 2, . . . , 6 defines a vector of real tracking error In order to derive a backstepping controller, an auxiliary desired state, Z aux ∈ R 6 , used for a virtual control system, is presented as follows: where k 1 = diag(k 11 , k 1i , . . . , k 16 ), k 2 = diag(k 21 , k 2i , . . . , k 26 ), k 1i , k 2i ∈ R + , i = 1, . . . , 6 are diagonal matrices of controller gains. Time derivative the auxiliary state is computed from Equation (26) as follows: .  In order to derive a backstepping controller, an auxiliary desired state, 6 aux  Z , used for a virtual control system, is presented as follows: From Equations (25) and (28), it can be seen that where 66    is positive diagonal matrix given by, Let z 2 = [ z 21 , z 2i , . . . , z 26 ] T ∈ R 6 denote a vector of virtual tracking error computed by the following equation From Equations (25) and (28), it can be seen that where Ψ ∈ R 6×6 is positive diagonal matrix given by, Substituting Equations (28) and (29), the virtual tracking error z 2 is computed by: where P ∈ R 6×6 is positive diagonal matrix given by: The time derivative of z 2 is obtained from Equations (7) and (31) as follows, In order to stabilize the attitude control and also guarantee a good trajectory tracking of position control for a quadcopter UAV, the proposed robust backstepping controller-based ESDO, u(t), can be designed as the following equation, The stability analysis of the control system by using the controller u(t) is presented in the next section.

Stability Analysis of the Proposed Robust Backstepping Controller
The stability analysis of the control system is presented in this subsection as follows, implies that for any finite constant ε.
The solution, W(t), of differential Equation (37) can be achieved by [39] as follows The proof of Lemma 1 is completed.

Theorem 2.
Consider the disturbed dynamic system (7). The proposed controller presented in Equation (34) can guarantee the ultimate convergence of the system output states of a quadcopter UAV to desired trajectories, i.e., the trajectory tracking errors z 1 , z 2 are bounded by a small region neighborhood of zero.

Proof of Theorem 2. A Lyapunov function candidate for the trajectory tracking error is chosen by,
It is easy to achieve the first-time derivative of the function V( z 1 , z 2 ) presented in Equation (39) from Equations (29) and (33) as follows, Substituting the proposed controller, u(t), given in Equation (34), and matrix P in Equation (32) into Equation (40) and doing some mathematical conversion steps, the value of . V( z 1 , z 2 ) can be obtained as follows, where M, N ∈ R 6×6 are positive diagonal constant matrices, its values can be computed from Equations (30) and (32)  , i = 1, 2, . . . , 6 (42) In addition, the values of z T 1 M z 1 , z T 2 N z 2 , and z T 2 P ξ from Equation (41) can be computed from Equations (32), (42), and (43) as follows, Using Young's inequality for | z 2i | and ξ i , the following expression is always satisfied, From Equations (44)-(47), the derivative function .
V( z 1 , z 2 ) given in Equation (41) can be obtained as follows, Let ρ, η ∈ R + be positive constants, which always satisfy the following conditions, Hence, from Equations (24), (39) and (47)-(49), the below inequality always satisfies, According to Lemma 1, the solution V(t) for inequality (50) can be obtained as follows, From Equation (51), it is obvious, when t → +∞ , the function V(t) is constantly convergent by Theorem 2 is completely proven. (52), we can conclude that the function V( z 1 , z 2 ) is bounded when t → +∞ , and its value converges to ηµ 2 /(4ρ). It implies that the trajectory tracking errors z 1 , z 2 are convergent to a small region containing the origin zero.

Simulation Results and Discussions
In this Section, a numerical simulation is carried out to verify the effectiveness of the proposed Robust Backstepping Control based on Extended State/Disturbance Observer for a quadcopter UAV through a scenario of trajectory tracking test, where different types of external perturbations and parametric uncertainties are considered. A comparative flight simulation between the proposed controller and other existent algorithms such as Sliding Mode Control [26] and Active Disturbance Rejection Control (ADRC) [38] is conducted on the same vehicle model with an identical working environment to emphasize the superior performance of the proposed method. The physical system parameters of a quadcopter UAV and the controller coefficients to perform a numerical simulation are provided in Tables 1 and 2, respectively, in which the controller parameters given in Table 2 are selected to ensure the rapid response of the UAV system and minor trajectory tracking error. A scenario of numerical simulation can be described as follows: the quadcopter UAV starts to takeoff from an initial state Z 0 = [x 0 , y 0 , z 0 , φ 0 , θ 0 , ψ 0 ] T = [8, 0, 0, 0, 0, 0.1] T . A control mission of the vehicle is to drive the UAV's position closely tracking a typical desired trajectory given by Equation (53). During the working process, the flight of the quadcopter is influenced by randomly lumped perturbations.
   x d = 5 + 5 cos t/ 1 + sin 2 t y d = 5 sin t cos t/ 1 + sin 2 t z d = 2 (53) From the simulation result shown in Figure 3a,b, the tracking performance of the proposed RBCESDO is significantly better than that of the Sliding Mode Control [26] and Active Disturbance Rejection Control [38] in both two-dimensional (2D) and threedimensional (3D) trajectories in general. According to Figure 3a, the RBCESDO algorithm exhibits a rapid altitude achievement (z = 2 m). In addition, the horizontal position tracking (x and y) also responds faster and more accurately than the other considered methods. As a result, the vehicle position rapidly tracks the desired trajectory at x = 9 m, y = 1.5 m, while the corresponding position is slowly convergent by using SMC and ADRC because of the powerfully approximate capacity of the proposed extended state/disturbance observer. Following the result shown in Figure 4, the undesirable lumped perturbations' influence on the attitude dynamic ξ φ , ξ θ , ξ ψ and on the position dynamic ξ x , ξ y , ξ z are closely estimated by ξ φ ,ξ θ ,ξ ψ and ξ x ,ξ y ,ξ z , respectively. These approximated values are integrated with the proposed backstepping controller given in Equation (34) to compensate for the effects of disturbances/uncertainties. Thus, the performance of attitude states (roll, pitch, and yaw angle) and position states (x, y, and z) of the proposed RBCESDO is more precisely convergent to the reference trajectories compared with other methods as shown in Figures 5a-c and 6.    , respectively. These approximated values are integrated with the proposed backstepping controller given in Equation (34) to compensate for the effects of disturbances/uncertainties. Thus, the performance of attitude states (roll, pitch, and yaw angle) and position states (x, y, and z) of the proposed RBCESDO is more precisely convergent to the reference trajectories compared with other methods as shown in Figures 5a-c and 6.    (34) to compensate for the effects of disturbances/uncertainties. Thus, the performance of attitude states (roll, pitch, and yaw angle) and position states (x, y, and z) of the proposed RBCESDO is more precisely convergent to the reference trajectories compared with other methods as shown in Figures 5a-c and 6.   On the other hand, a general advantage of all three comparative methods (i.e., SMC, ADRC, and RBCESDO) is to eliminate the chattering effect of the controller, as shown in Figure 7. However, to remove the chattering phenomenon in the SMC technique, the high gain parameters of the switching law must be avoided. Thus, the SMC controller cannot be strictly guaranteed in both a tracking precision of the quadcopter UAV to a reference path and a strong anti-chattering effect. In this simulation, the moderate switching controller gains of the SMC method are selected instead of the large ones to ensure the anti-chattering capacity, resulting in the larger tracking errors of roll, pitch, and yaw angles, and larger error of positions are generated in comparison with the proposed RBCESDO method (as shown in Figure 5d-f for attitude control and Figure 6 for position control), while these fluctuation errors are more amplified with ADRC control algorithm as exhibited in Figures 5g-i and 6.
To emphasize a significantly superior performance of the proposed controller compared to the other algorithms, a convergence to zero of the attitude tracking error (φ error , θ error , ψ error ) and position tracking error (x error , y error , z error ) of the three considered methods (RBCESDO, SMC, and ADRC) is also portrayed in Figures 8 and 9, which exhibit that the proposed controller can obtain the best performance with the smallest tracking error and fastest response. More specifically, the maximum of the attitude tracking error of the RBCESDO reaches about Max{φ error } ≈ 0.00051 [rad] for roll, Max{θ error } ≈ 0.00079 [rad] for pitch, and Max{ψ error } ≈ 0.0013 [rad] for yaw (in Figure 9a) for yaw. In addition, the position control of UAV is also highly improved as exhibited in Figure 9b,  In summary, the proposed algorithm (RBCESDO) seems to provide a much better tracking performance for attitude and position control compared to the considered methods in both rapid response and minor tracking error.  On the other hand, a general advantage of all three comparative methods (i.e., SMC, ADRC, and RBCESDO) is to eliminate the chattering effect of the controller, as shown in Figure 7. However, to remove the chattering phenomenon in the SMC technique, the high gain parameters of the switching law must be avoided. Thus, the SMC controller cannot  On the other hand, a general advantage of all three comparative methods (i.e., SMC, ADRC, and RBCESDO) is to eliminate the chattering effect of the controller, as shown in Figure 7. However, to remove the chattering phenomenon in the SMC technique, the high gain parameters of the switching law must be avoided. Thus, the SMC controller cannot be strictly guaranteed in both a tracking precision of the quadcopter UAV to a reference   To emphasize a significantly superior performance of the proposed controller compared to the other algorithms, a convergence to zero of the attitude tracking error ( ) In summary, the proposed algorithm (RBCESDO) seems to provide a much better tracking performance for attitude and position control compared to the considered methods in both rapid response and minor tracking error.

Conclusions and Future Works
In this study, a robust backstepping control based on an extended state/disturbance observer for trajectory tracking is developed to improve the performance of the quadcopter UAV under the various influences of uncertainties and/or external perturbations. In the proposed algorithm, the ultimate control law is recursively constructed by using the

Conclusions and Future Works
In this study, a robust backstepping control based on an extended state/disturbance observer for trajectory tracking is developed to improve the performance of the quadcopter UAV under the various influences of uncertainties and/or external perturbations. In the proposed algorithm, the ultimate control law is recursively constructed by using the

Conclusions and Future Works
In this study, a robust backstepping control based on an extended state/disturbance observer for trajectory tracking is developed to improve the performance of the quadcopter UAV under the various influences of uncertainties and/or external perturbations. In the proposed algorithm, the ultimate control law is recursively constructed by using the ESDO method and defined auxiliary variables. The Lyapunov theorem is used to analyze and prove the robust stability and powerful convergence of the system states to reference paths. Finally, the numerical simulation is fully executed in both attitude and position control mode; these simulation results are also compared with the Sliding Mode Control and Active Disturbance Rejection Control methods to strongly confirm the superior performance of the presented control scheme in both desired trajectory convergence and anti-disturbance capacity. In addition, precise knowledge of the upper bounds of uncertainties and/or external disturbances is not demanded during the steps of designing the proposed controller.
Another significant contribution of the proposed algorithm comes from the real aspect. The ESDO is able to estimate both the velocity state of the UAV's system and wind speed once these parameters cannot be directly measured due to sensor limitations. Thus, the real implementation cost for data acquisition may be reduced, and the influences of sensor noise are also removed. Future works are going to discover a finite-time convergence of the system states to a desired trajectory.