Next Article in Journal
Air-to-Ground Path Loss Model at 3.6 GHz under Agricultural Scenarios Based on Measurements and Artificial Neural Networks
Next Article in Special Issue
On the Dynamics of Flexible Wings for Designing a Flapping-Wing UAV
Previous Article in Journal
An Attention-Based Odometry Framework for Multisensory Unmanned Ground Vehicles (UGVs)
Previous Article in Special Issue
A Hybrid Global/Reactive Algorithm for Collision-Free UAV Navigation in 3D Environments with Steady and Moving Obstacles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fixed-Time Extended Observer-Based Adaptive Sliding Mode Control for a Quadrotor UAV under Severe Turbulent Wind

by
Armando Miranda-Moya
1,
Herman Castañeda
1,* and
Hesheng Wang
2
1
Tecnologico de Monterrey, School of Sciences and Engineering, Monterrey 64849, Mexico
2
Department of Automation, Key Laboratory of System Control and Information Processing of Ministry of Education, Key Laboratory of Marine Intelligent Equipment and System of Ministry of Education, Shanghai Engineering Research Center of Intelligent Control and Management, Shanghai Jiao Tong University, Shanghai 200240, China
*
Author to whom correspondence should be addressed.
Drones 2023, 7(12), 700; https://doi.org/10.3390/drones7120700
Submission received: 22 October 2023 / Revised: 27 November 2023 / Accepted: 28 November 2023 / Published: 9 December 2023
(This article belongs to the Special Issue Conceptual Design, Modeling, and Control Strategies of Drones-II)

Abstract

:
This paper presents a fixed-time extended state observer-based adaptive sliding mode controller evaluated in a quadrotor unmanned aerial vehicle subject to severe turbulent wind while executing a desired trajectory. Since both the state and model of the system are assumed to be partially known, the observer, whose convergence is independent from the initial states of the system, estimates the full state, model uncertainties, and the effects of turbulent wind in fixed time. Such information is then compensated via feedback control conducted by a class of adaptive sliding mode controller, which is robust to perturbations and reduces the chattering effect by non-overestimating its adaptive gain. Furthermore, the stability of the closed-loop system is analyzed by means of the Lyapunov theory. Finally, simulation results validate the feasibility and advantages of the proposed strategy, where the observer enhances performance. For further demonstration, a comparison with an existent approach is provided.

1. Introduction

External disturbances due to wind produce adverse effects on quadrotor unmanned aerial vehicles (UAVs) when flying outdoors. Depending on its intensity, the rotorcraft can suffer from small erratic movements to complete system instability. To avoid the latter, advanced robust control strategies such as [1,2,3,4,5] have been proposed to handle the nonlinearities of the UAV in perturbed conditions. However, their rejection against strong enough disturbances tends to be slow since it relies on a feedback regulation defined by the error between the states of the quadrotor and their respective setpoints [6]. Moreover, it has been proven in [7,8] that the vehicle can maintain its flight stability but cannot execute proper trajectory tracking. Thus, approaches such as [9,10,11,12] employ a disturbance rejection method where the total disturbance affecting the system is estimated through a disturbance observer and later compensated on the control inputs of the robot.
Frequently, external perturbations are established as step or periodic functions [13,14,15,16,17]. However, such functions do not reflect the actual wind dynamics. Although wind evolution is hard to describe mathematically due to its random conduct, the US Military standard MIL-STD-1797A [18] endorses the Dryden and Von Karman wind models to analyze the performance of aerial vehicles in environments that involve turbulent air behavior at different altitude levels. These numerical approximations are implemented in Rodriguez-Mata et al. [19] where a nonlinear high-gain disturbance observer-based PD attitude controller is designed for a quadrotor UAV. The effects of the wind disturbances represented by the Dryden model are analyzed through the evolution of the torque of the motors. In Zhang et al. [20], disturbance observer-based backstepping and sliding mode controllers are developed to manage the linear and angular dynamics of a rotorcraft, respectively. Furthermore, the approaches from Aboudonia et al. [21,22] present a disturbance observer-linear feedback control scheme for a quadrotor UAV susceptible to mismatched perturbations. The low-altitude turbulent wind effects are analyzed in the mathematical model of the rotorcraft as drag forces that actuate over the vehicle. Moreover, the same authors propose a disturbance observer-sliding mode strategy for a quadrotor UAV considering partial motor failure. The wind turbulent field given by the Von Karman model is also analyzed in the form of drag forces [23].
Previous designs require the complete knowledge of the full-state measurements of the system which may not be available due to a lack of instruments or sensor noise conditions. Then, in [24], Han proposes an active disturbance rejection controller where an extended state observer (ESO) provides estimations of both the full-state and slow time-varying perturbations and uncertainties by adding an extra state to the mathematical model of the plant. Following the prior approach, the authors of [25] propose a backstepping-sliding mode strategy to control a quadrotor using the measurement and disturbance estimations from a linear ESO. Additionally, Shao et al. [26] introduce an ESO-based dynamic surface controller that improves the performance of a backstepping scheme. However, its robustness against external disturbances can still be enhanced by employing a sliding mode controller instead. Then, in ref. [27], the authors present a nonlinear extended state observer-based fuzzy sliding mode technique to drive the pose of a rotorcraft. Considering a first-order sliding method as the auxiliary controller, a fuzzy logic approach is implemented to attenuate the chattering of the control inputs of the vehicle. Nonetheless, the effectiveness of this strategy in reducing the chattering effect strongly depends on an appropriate number of fuzzy rules and heuristics which may not provide results as good as a defined adjustment function.
On the other hand, to guarantee the complete stability of the quadrotor UAV under an observer-based control scheme, the convergence of the observer must be much faster than the convergence of the controller. Extended state observers with asymptotic and finite-time equilibrium achievement are proposed in [28,29], respectively. Even though finite-time strategies offer a better performance regarding estimation and error reduction, their convergence time increases as the initial conditions of the error of the system rise. In addition, these values could be unknown during practice. As a solution, in [30], Polyakov introduces a fixed-time stable system whose convergence is independent from its initial conditions. This approach has been adopted for the design of ESO-based controllers in [31,32,33,34]. In that sense, Mechali et al. [35] make a complete analysis of the feasibility of fixed-time ESO-based terminal sliding mode schemes for the position and attitude driving of a quadrotor UAVs subject to external disturbances described by the Dryden wind model. Similarly, Li et al. [36] introduce a fixed-time ESO-based dynamic surface control whose design follows a fast terminal sliding mode approach. Nevertheless, different from real conditions, the authors consider external disturbances in the form of periodic functions only. Finally, Liu et al. [37] present the design of an ESO-based algorithm for a fixed-wing UAV. The observer is employed to estimate the uncertainties that affect the vehicle when performing a high angle of attack. Nonetheless, prior approaches rely on fixed-gain control strategies with heuristically chosen parameters. This may lead to the implementation of larger gains than necessary, resulting in excessive energy delivery and possible overfitting of the controller. To avoid the latter, adaptive strategies can be designed to provide the proper amount of control effort to the system without overestimating the control parameters.
Then, the main contribution of this paper relies on the design of a fixed-time extended state observer (FxtESO)-based adaptive sliding mode controller (ASMC) for a quadrotor UAV subject to severe wind turbulent gusts described by the Von Karman model. Through the assumption of partial state measurements, the FxtESO provides estimations regarding the full state of the UAV, the complex random behavior of the turbulent wind gust, and model uncertainties in fixed time, independently of the initial state and without complete knowledge of the system. The flight control is performed by the ASMC, which offers properties such as robustness against bounded external perturbations and uncertainties, chattering attenuation due to the adjustment of its adaptive gain, and practical finite-time state convergence. Furthermore, the closed-loop stability of the system is guaranteed through an analysis based on the Lyapunov theory. Simulation results demonstrate the improvements in the performance of the rotorcraft due to the implementation of the observer. Moreover, the effectiveness and advantages such as robustness, estimation of the full state and perturbations and its adaptive control gain, which modulates the control input as needed of the proposed approach are compared with an existing observed-based control methodology.
The paper is divided as follows: Section 2 addresses the dynamics of the quadrotor UAV and the Von Karman turbulent wind model. The design of the FxtESO-based ASMC scheme is described in Section 3. Simulation results are revealed and discussed in Section 4. Finally, the conclusions are presented in Section 5.

2. Preliminaries

The definitions and lemmas that sustain the FxtESO-based ASMC scheme are addressed. First, we assume a scalar state x R and a nonlinear function sig r ( x ) = | x | r sign ( x ) , with r ( 0 , 1 ) , | x | as the absolute value of x, and
sign ( x ) = 1 if x < 0 0 if x = 0 1 if x > 0 .
Now, we consider the following single-input single-output (SISO) nonlinear system with control input U R :
x ˙ = f ( x , U ) , x 0 = x ( t 0 ) .
Definition 1
([38]). The equilibrium point x = 0 is asymptotically stable if there is a continuously differentiable function V that fulfills V ( 0 ) = 0 , and V > 0 , V ˙ < 0 for all x { 0 } that exists in the domain of solutions of (2).
Definition 2
([39]). The origin of (2) is finite-time stable if it is asymptotically stable and any solution of the system achieves the equilibrium point x = 0 after a time T m ( x 0 ) . In other words, x ( t , x 0 ) = 0 , t T m ( x 0 ) .
Definition 3
([40]). The solution of system (2) is said to be practical finite-time stable if, with x 0 , there is value Γ > 0 and time T k ( Γ , x 0 ) < that satisfy | x ( t ) | < Γ , t t 0 + T k ( Γ , x 0 )
Definition 4
([41,42]). The origin of (2) is fixed-time stable if it is globally finite-time stable and there is a bounded settling time function T f ( x ) such that 0 < T f T m a x .
Lemma 1
([43]). System
ϑ ˙ 1 = ϑ 2 ϰ 1 ϑ 1 ϖ 1 sign ( ϑ 1 ) ϑ ˙ 2 = ϑ 3 ϰ 2 ϑ 1 ϖ 2 sign ( ϑ 1 ) ϑ ˙ n 1 = ϑ n ϰ n 1 ϑ 1 ϖ n 1 sign ( ϑ 1 ) ϑ ˙ n = ϰ n ϑ 1 ϖ n sign ( ϑ 1 ) ϖ 1 = ϖ ¯
is homogeneous of degree ϖ ¯ 1 with respect to weights { ( i 1 ) ϖ ¯ ( i 2 ) } for 1 i n and ϖ i = i ϖ ¯ ( i 1 ) for 2 i n , if ϖ ¯ > 1 1 n 1 .
Lemma 2
([44]). We assume that vector field ν f is homogeneous of degree k with respect to vector field ν g . Then, the origin is a finite-time-stable equilibrium under ν f if and only if the origin is an asymptotically stable equilibrium under ν f and k < 0 .
Lemma 3
([30]). We assume that there is a continuous unbounded function V. We set Θ s R n ; it is globally fixed-time attractive if V ( q ) = 0 q Θ s and any q ( t ) satisfies V ˙ ( q ) ( α V a ( q ( t ) ) + β V b ( q ( t ) ) ) c for α , β , a , b , c > 0 , with a c < 1 , b c > 1 , and settling time function
T f 1 α c ( 1 a c ) + 1 β c ( b c 1 ) .
Lemma 4
([40,45]). We consider system (2) and suppose that there is a continuous function V ( x ) that fulfills V ( 0 ) = 0 and V ( x ) R + , x 0 . Then, the origin of (2) is practical finite-time stable if
V ˙ ( x ) h V ( x ) L + W ,
where h , W > 0 and 0 < L < 1 . Moreover, the trajectories of x are bounded in finite-time as
lim Y Y 0 x V L ( x ) W ( 1 L ) h ,
with 0 < Y 1 , 0 < Y 0 < 1 and settling time function
T k V 1 L ( x 0 ) h Y 0 ( 1 L ) .

3. UAV Dynamics and the Von Karman Wind Model

In this section, the mathematical expressions that describe both rotorcraft and external disturbances in the form of turbulent wind are introduced. The Newton–Euler dynamics of the quadrotor contemplate an inertial reference frame ( X 0 , Y 0 , Z 0 ) with origin at a point in the Earth and a body reference frame ( X 1 , Y 1 , Z 1 ) with origin in the center of mass of the UAV (see Figure 1). Moreover, the angular dynamics is considered to be the same in both frames since we assume small roll and pitch rotations [46]. Then, the equations of motion of the aerial vehicle are defined as
d ˙ = Rv ,
m v ˙ = f m ω × v ,
f = T h e 3 + R 1 m g e 3 + R 1 δ ,
Θ ˙ = ω ,
ω ˙ = J 1 ( τ ω × J ω ) ,
where d = [ x , y , z ] T represents the inertial frame linear position of the UAV, v = [ u , v , w ] T is the body frame linear velocity, Θ = [ ϕ , θ , ψ ] T is the attitude of the UAV, while ω is its angular velocity; R is the body to inertial frame rotation matrix, vector f involves the linear forces that influence the quadrotor, δ = [ δ x , δ y , δ z ] T describes the external disturbances that affect the quadrotor, T h is the total thrust of the four motors, i.e., T h = T h 1 + T h 2 + T h 3 + T h 4 , m is the mass of the aerial vehicle, g is the gravity acceleration, and e 3 = [ 0 , 0 , 1 ] T is a direction vector. Furthermore, assuming mechanical symmetry, the inertia matrix J is established as
J = J x x 0 0 0 J y y 0 0 0 J z z .
Finally, the torques generated by the motors are defined in vector τ = [ τ ϕ , τ θ , τ ψ ] T , where
τ ϕ τ θ τ ψ = l sin ( π / 4 ) ( T h 1 T h 2 T h 3 + T h 4 ) l sin ( π / 4 ) ( T h 1 + T h 2 T h 3 T h 4 ) τ 1 + τ 2 τ 3 + τ 4 ,
with l as the distance from the center of each motor to the center of mass of the UAV, sin ( π / 4 ) as the angle between them, and T h i , τ i for i = 1 , 2 , 3 , 4 as the thrust and torque from each motor, respectively, which can be expressed as
T h i = C T Ω i 2 ,
τ i = C τ Ω i 2 .
with C T as the thrust constant, C τ as the moment constant, and Ω i as the angular velocity of motor i. Now, assuming the existence of model uncertainties regarding the angular dynamics and the mass of the rotorcraft, the inertial frame equations of motion of the quadrotor are rewritten in the scalar form as
ϕ ¨ = τ ϕ J x x ,
θ ¨ = τ θ J y y ,
ψ ¨ = τ ψ J z z ,
x ¨ = T h m u ( s ϕ s ψ + c ϕ c ψ s θ ) + δ x ( t ) m u ,
y ¨ = T h m u ( c ϕ s ψ s θ + c ψ s ϕ ) + δ y ( t ) m u ,
z ¨ = T h m u ( c ϕ c θ ) + g + δ z ( t ) m u ,
where m u is a nominal value of m such that m u m while s and c are the sine and cosine functions, respectively. External disturbance terms δ x , δ y , and δ z are analyzed as drag forces due to turbulent wind through the following equation [23]:
δ * = 1 2 ρ B * A * ( d ˙ * v w ) 2 sign ( d ˙ * v w ) ,
where = { x , y , z } , ρ is the air density, B is the drag coefficient in the inertial frame axis; A is the area of the quadrotor in each axis plane, and v w is the inertial frame wind velocity given by
v w x = v ¯ w x + ζ u , v w y = v ¯ w y + ζ v , v w z = v ¯ w z + ζ w ,
where v ¯ w is the inertial frame mean wind velocity while ζ u , ζ v and ζ w are wind deviations described by the Von Karman wind turbulence model [47]:
ζ u ( S ) n u ( S ) = η u 2 L u π V ( 1 + 0.25 L u V S ) 1 + 1.357 L u V S + 0.1987 ( L u V ) 2 S 2 ,
ζ v ( S ) n v ( S ) = η v L v π V ( 1 + 2.7478 L v V S + 0.3398 ( L v V ) 2 S 2 ) 1 + 2.9958 L v V S + 1.9754 ( L v V ) 2 S 2 + 0.1539 ( L v V ) 3 S 3 ,
ζ w ( S ) n w ( S ) = η w L w π V ( 1 + 2.7478 L w V S + 0.3398 ( L w V ) 2 S 2 ) 1 + 2.9958 L w V S + 1.9754 ( L w V ) 2 S 2 + 0.1539 ( L w V ) 3 S 3 ,
with white noise signals n u , n v , n w and low-altitude parameters
L w = z ,
L u = L v = z ( 0.177 + 0.000823 z ) 1.2 ,
and
η w = 0.1 W 20 ,
η u = η v = η w ( 0.177 + 0.000823 z ) 0.4 ,
where L u , L v , and L w are turbulence scale lengths, and
V = ( d ˙ x v w x ) 2 + ( d ˙ y v w y ) 2 + ( d ˙ z v w z ) 2 .
Terms η u , η v , and η w depend on W 20 , whose value denotes the wind velocity at an altitude of 6 m according to regulation MIL-HDBK-1797 [18] and specifies whether the turbulence intensity is light ( W 20 = 15 knots), moderate ( W 20 = 30 knots), or severe ( W 20 = 45 knots).

4. FxtESO-Based ASMC Design

This section describes the design of a fixed-time extended state observer-based adaptive sliding mode strategy. The control objective is to drive a quadrotor UAV while tracking a desired trajectory in the presence of turbulent wind gusts. Assuming that the position measurements of the quadrotor are the only available states, the FxtESO supplies full-state and disturbance estimations within fixed-time convergence, whereas the ASMC drives the pose of the quadrotor while providing practical finite-time convergence of the state, robustness against bounded external disturbances and model uncertainties, and chattering attenuation. Then, the flight control scheme is composed of an inner loop that deals with the fast dynamics of the vehicle and an outer loop that handles the linear position of the UAV and the external disturbances rejection; see Figure 2.

4.1. Fixed-Time Extended State Observer

We consider the following SISO nonlinear system:
ξ ˙ 1 = ξ 2 ,
ξ ˙ 2 = a ( ξ ) + b ( ξ ) u r ,
c = ξ 1 ,
where ξ = [ ξ 1 , ξ 2 ] T denotes the state variables, u r is the control input, c is the output of the system, and a ( ξ ) , b ( ξ ) are continuously differentiable functions. Now, extending to a third-order system, we have
ξ ˙ 1 = ξ 2 ,
ξ ˙ 2 = ξ 3 + b ^ ( ξ ) u r ,
ξ ˙ 3 = γ ( ξ ) ,
c = ξ 1 ,
where | b ^ ( ξ ) | > 0 is a bounded nominal model of b ( ξ ) , and ξ 3 = Δ δ + a ^ ( ξ ) is a continuously differentiable and globally bounded augmented state that describes the total disturbance composed of those external perturbations and model uncertainties affecting the system. Finally, | γ ( ξ ) | 0 denotes the dynamics of the extended state. Moreover, estimations of state variables ξ 1 , ξ 2 and ξ 3 are given by FxtESO [48]:
ξ ^ ˙ 1 = ξ ^ 2 + κ 1 sig α 1 ( ξ ˜ ) + ε 1 sig β 1 ( ξ ˜ ) , ξ ^ ˙ 2 = ξ ^ 3 + κ 2 sig α 2 ( ξ ˜ ) + ε 2 sig β 2 ( ξ ˜ ) + b ^ ( ξ ) u r , ξ ^ ˙ 3 = κ 3 sig α 3 ( ξ ˜ ) + ε 3 sig β 3 ( ξ ˜ ) + χ sign ( ξ ˜ ) ,
with estimation error ξ ˜ = ξ 1 ξ ^ 1 , parameters κ 1 , κ 2 , κ 3 , ε 1 , ε 2 , ε 3 R + , and χ > γ ( ξ ) . Considering a small enough scalar φ 1 , parameters α 1 , 2 , 3 are defined as α 1 ( 1 φ 1 , 1 ) , α 2 = 2 α 1 1 , and α 3 = 3 α 1 2 . In that sense, by choosing φ 2 to be sufficiently small, terms β 1 , 2 , 3 are established as β 1 ( 1 , 1 + φ 2 ) , β 2 = 2 β 1 1 , and β 3 = 3 β 1 2 . Furthermore, by taking Equation (36), the control law of the augmented system yields
u r = ξ ^ 3 + ν a b ^ ( ξ ^ ) ,
with ν a as an auxiliary control strategy designed to control the positioning and orientation of the quadrotor UAV.

4.2. Attitude Controller

In order to design a controller for the rotational motion, the angular dynamics equations of the quadrotor are rewritten as
ϖ ˙ = a ( ϖ ) + b ( ϖ ) τ ,
c a = ς a ( ϖ ) ϖ ,
with states vector ϖ = [ ϕ , ϕ ˙ , θ , θ ˙ , ψ , ψ ˙ ] T , control input τ = [ τ ϕ , τ θ , τ ψ ] T , ς a ( ϖ ) R 3 × 6 , and c a = [ ϕ , θ , ψ ] T as the available state. Introducing the desired states ϖ d = [ ϕ d , ϕ ˙ d , θ d , θ ˙ d , ψ d , ψ ˙ d ] T and the error vector e a = [ ϕ ϕ d , θ θ d , ψ ψ d ] T , a sliding surface vector is designed as
σ ϖ = λ a e a + e ˙ a ,
with λ a = diag ( λ ϕ , λ θ , λ ψ ) and λ ϕ , θ , ψ > 0 . Then, the changing rate of (43) yields
σ ˙ ϖ = λ a e ˙ a + e ¨ a .
Equations (16)–(18) are substituted in (44) to perform a feedback linearization in the form of (40) as follows [49]:
τ ϕ τ θ τ ψ = J x x ( λ ϕ ( ϕ ˙ ϕ ˙ d ) + ϕ ¨ d ξ ^ 3 , ϕ + ν ϕ J y y ( λ θ ( θ ˙ θ ˙ d ) + θ ¨ d ξ ^ 3 , θ + ν θ J z z ( λ ψ ( ψ ˙ ψ ˙ d ) + ψ ¨ d ξ ^ 3 , ψ + ν ψ ,
where ξ ^ 3 , ϕ , ξ ^ 3 , θ , and ξ ^ 3 , ψ are the total disturbance estimates for the attitude subsystem while ν ϕ , ν θ and ν ψ are auxiliary controllers described by the adaptive sliding mode strategy
ν = K 1 ( t ) sig 1 / 2 ( σ ) K 2 σ ,
with = { ϕ , θ , ψ } , K 2 > 0 , and K 1 as the adaptive gain with adaptability law
K ˙ 1 ( t ) = k sign ( σ μ ) K 1 > k min k min K 1 k min ,
where k > 0 is the adaptation rate regulator, μ > 0 is the adaptiveness threshold, and k min R + ensures the delivery of the lowest control effort to the system. Then, K 1 increases its magnitude when condition | σ | > μ is met. On the other hand, once | σ | μ is true, the adaptive gain starts decreasing until achieving the value of k min . Moreover, the chattering attenuation property of the ASMC can be attributed to two factors: first, structure sig 1 / 2 ( · ) guarantees smooth and continuous responses from (46), and second the non-overestimation of gain K 1 due to its designed dynamics (47).

4.3. Positioning Controller

We consider the translation dynamics of the quadrotor as
ϱ ˙ = a p ( ϱ ) + b p ( ϱ ) T h + o ( t ) ,
c p = ς p ( ϱ ) ϱ ,
with state vector ϱ = [ x , x ˙ , y , y ˙ , z , z ˙ ] T , control input T h , which is equivalent to the total propulsion force of the motors, o ( t ) = 1 / m [ 0 , δ x ( t ) , 0 , δ y ( t ) , 0 , δ z ( t ) ] T as the external disturbances vector, ς p ( ϖ ) R 3 × 6 , and c p = [ x , y , z ] T as the available state. Moving on, the desired linear state and the position error vectors are, respectively, defined as ϱ d = [ x d , x ˙ d , y d , y ˙ d , z d , z ˙ d ] T and e p = [ x x d , y y d , z z d ] T . Hence, the sliding surface is established as
σ ρ = λ p e p + e ˙ p ,
where λ p = diag ( λ x , λ y , λ z ) and λ x , y , z > 0 . Thus,
σ ˙ ρ = λ p e ˙ p + e ¨ p .
Since the translation dynamics of the quadrotor are underactuated, feedback linearization cannot be performed directly after substituting (19)–(21) in (51). Therefore, the following virtual auxiliary control strategy is adopted:
ν v x ν v y ν v z = ( T h / m u ) ( s ϕ ^ s ψ ^ + c ϕ ^ c ψ ^ s θ ^ ) ( T h / m u ) ( c ϕ ^ s ψ ^ s θ ^ + c ψ ^ s ϕ ^ ) ( T h / m u ) ( c ϕ ^ c θ ^ ) .
Recalling that = { x , y , z } , virtual controller ν v is described by the following ASMC with disturbance compensation:
ν v = K 1 ( t ) sig 1 / 2 ( σ ) K 2 σ ξ ^ 3 , ,
where the adaptability of K 1 was previously defined in (47), and ξ ^ 3 , is the total disturbance estimation provided by FxtESO. Additionally, to fulfill the desired x and y trajectories, the desired roll and pitch rotations are determined by
θ d = asin ( m u / T h ) ν v x s ψ d s ϕ d ) c ψ d c ϕ d ,
ϕ d = asin ( m u / T h ) ( s ψ d ν v x c ψ d ν v y ) .
Next, the z-dynamics are stabilized by control input T h as
T h = m u c ϕ ^ c θ ^ ν v z ξ ^ 3 , z + z ¨ d λ z ( z ^ ˙ z ˙ d ) .

4.4. Stability Analysis

The stability analysis of the extended state observer-based adaptive controller is addressed. Then, the following Theorem is stated:
Theorem 1.
The sliding mode controller (46) with adaptive gain dynamics (47) using the fixed-time estimations from (39) in a closed loop with a nonlinear system (7)–(11) offers a practical finite-time convergence of the state error at bounded time t T k ( σ 0 ) .
Proof. 
In order to provide the stability, we have consider three stages:
FxtESO Convergence Analysis
We consider estimation errors
ξ ˜ i = ξ i ξ ^ i , for i = { 1 , 2 , 3 }
with time-derivative
ξ ˜ ˙ 1 = ξ ˜ 2 κ 1 sig α 1 ( ξ ˜ 1 ) ε 1 sig β 1 ( ξ ˜ 1 ) , ξ ˜ ˙ 2 = ξ ˜ 3 κ 2 sig α 2 ( ξ ˜ 1 ) ε 2 sig β 2 ( ξ ˜ 1 ) , ξ ˜ ˙ 3 = γ ( ξ ) κ 3 sig α 3 ( ξ ˜ 1 ) ε 3 sig β 3 ( ξ ˜ 1 ) χ sign ( ξ ˜ 1 ) .
Then, to prove the fixed-time convergence of the observer, the analysis of the estimation error dynamics (58) is divided into two sections [50]. First, we consider subsystem
ξ ˜ ˙ 1 = ξ ˜ 2 κ 1 sig α 1 ( ξ ˜ 1 ) , ξ ˜ ˙ 2 = ξ ˜ 3 κ 2 sig α 2 ( ξ ˜ 1 ) , ξ ˜ ˙ 3 = κ 3 sig α 3 ( ξ ˜ 1 ) ,
and terms
K = κ 1 1 0 κ 2 0 1 κ 3 0 0 ; ξ ˜ α = [ ξ ˜ 1 1 ϱ 1 , ξ ˜ 2 1 ϱ 2 , ξ ˜ 3 1 ϱ 3 ] T ,
with ϱ i = ( i 1 ) α ¯ ( i 2 ) for i = 1 , 2 , 3 . Referring to the work from Perruquetti et al. [43], system ξ ˜ ˙ α = K ξ ˜ α is achievable by choosing α ¯ = 1 . Since K is Hurwitz, ξ ˜ ˙ α is asymptotically stable. Then, there is a Lyapunov function V ( ξ ˜ α ) = ξ ˜ α T P ξ ˜ α that fulfills the following:
  • Matrix P is the solution of Lyapunov equation K T P + P K = I with identity matrix I R 3 × 3 .
  • V ˙ ( x ˜ α ) = ξ ˜ ˙ α T P ξ ˜ α + ξ ˜ α T P ξ ˜ ˙ α = ξ ˜ α T K T P + ξ ˜ α T P K ξ ˜ α = ξ ˜ α T ( K T P + P K ) ξ ˜ α = ξ ˜ α T I ξ ˜ α < 0 .
Moreover, Equation (61) is also valid for α ¯ ( 1 φ 1 , 1 ) , 0 < φ 1 < 1 [51]. Thus, V ( ξ ˜ α ) is a Lyapunov candidate function for error Subsystem (59). In addition, according to Lemma 1, (59) is a homogeneous system of degree α ¯ 1 < 0 . Consequently, Lemma 2 holds for this case and
V ˙ ( ξ ˜ α ) = ς α V ι 1 ( ξ ˜ α ) ,
where ς α > 0 and 0 < ι 1 < 1 . Now, we consider the second half of (58):
ξ ˜ ˙ 1 = ξ ˜ 2 ε 1 sig β 1 ( ξ ˜ 1 ) , ξ ˜ ˙ 2 = ξ ˜ 3 ε 2 sig β 2 ( ξ ˜ 1 ) , ξ ˜ ˙ 3 = ε 3 sig β 3 ( ξ ˜ 1 ) ,
and we let
E = ε 1 1 0 ε 2 0 1 ε 3 0 0 ; ξ ˜ β = [ ξ ˜ 1 1 υ 1 , ξ ˜ 2 1 υ 2 , ξ ˜ 3 1 υ 3 ] T ,
with υ j = ( j 1 ) β ¯ ( j 2 ) for j = 1 , 2 , 3 . Similarly to the previous analysis, if β ¯ = 1 , system ξ ˜ ˙ β = E ξ ˜ β is asymptotically stable. Thus, there is a Lyapunov function V a ( ξ ˜ β ) = ξ ˜ β T P a ξ ˜ β with derivative V ˙ a ( ξ ˜ β ) = ξ ˜ β T I ξ ˜ β < 0 , which also holds when β ¯ ( 1 , 1 + φ 2 ) with small enough φ 2 . By means of Lemma 1, Subsystem (63) is homogeneous of degree β ¯ 1 > 0 . Furthermore, referring to Basin et al. [51], the subsystem can be written in the form
V ˙ ( ξ ˜ β ) = ς β V a ι 2 ( ξ ˜ β ) ,
with ς β > 0 and ι 2 > 1 . Hence, there is Lyapunov function V b ( ξ ˜ ) for error System (58) with derivative V ˙ b ( ξ ˜ ) = ( ς α V b ι 1 ( ξ ˜ ) + ς β V b ι 2 ( ξ ˜ ) ) . Then, according to Lemma 3, observer error ξ ˜ = [ ξ ˜ 1 , ξ ˜ 2 , ξ ˜ 3 ] T converges to zero in a fixed time defined by [52]
T f λ max 1 α ¯ ( P ) ς α ( 1 α ¯ ) + 1 ς β ( β 1 ) λ min ( P a ) β 1 ,
where λ max ( P ) is the maximum eigenvalue of P , and λ min ( P a ) is the minimum eigenvalue of P a . Moreover, since ξ ˜ 1 , ξ ˜ ˙ 1 0 as t T f , ξ ˜ 2 and ξ ˜ 3 also tend to zero. Then, ξ ˜ = 0 when t > T f , which means that the equivalent control identity given by
ξ ˜ ˙ 3 = γ ( ξ ) χ s i g n ( ξ ˜ 1 ) = 0
holds [53]. We notice that ξ ˜ 3 = 0 cannot be accomplished in practice due to several factors such as noise conditions, sampling frequency, and the disturbance dynamics [50]. However, assuming bounded perturbations and uncertainties, the fixed-time attractiveness of (39) allows the reaching of a small vicinity ϰ > 0 such that ξ ˜ 3 ϰ , after t > T f .
Adaptive Sliding Mode Control Stability
By taking Equations (40), (44), and (51), the closed-loop FxtESO-based ASMC system, in a single-input single-output form, is expressed as
σ ˙ = K 1 ( t ) sig 1 / 2 ( σ ) K 2 σ + ξ ˜ 3 ,
recalling that ξ ˜ 3 = ξ 3 ξ ^ 3 represents the estimation error for the total disturbance. Then, the Lyapunov candidate function for the stability analysis of the system is selected as [54]
V ( σ ) = 1 2 σ 2 + 1 2 ( K 1 ( t ) K M ) 2 ,
where K M is a value of adaptive gain K 1 that satisfies inequality ( K 1 ( t ) K M ) < 0 and [49]
K M ( t ) > 1 | σ | 1 2 ( δ m a x K 2 | σ | ) ,
0 < | ξ 3 | δ m a x .
While condition K 1 ( t ) > K m i n holds, the time differentiation of (69) is
V ˙ ( σ ) = σ σ ˙ + ( K 1 ( t ) K M ) ( k sign ( | σ | μ ) ) .
By substituting (68) and (71) in (72),
V ˙ ( σ ) | σ | ( K 1 ( t ) | σ | 1 2 sign ( σ ) K 2 σ + δ m a x ) + ( K 1 ( t ) K M ) ( k sign ( | σ | μ ) ) .
Now, terms ± K M | σ | 3 2 are introduced and applied as follows:
V ˙ ( σ ) | σ | ( K M | σ | 1 2 K 2 | σ | + δ m a x ) + ( K 1 ( t ) K M ) ( k sign ( | σ | μ ) | σ | 3 2 ) .
For the sake of simplicity, we let h p = K M | σ | 1 2 + K 2 | σ | δ m a x R + . Then, (74) yields
V ˙ ( σ ) h p | σ | + ( K 1 ( t ) K M ) ( k sign ( | σ | μ ) | σ | 3 2 ) .
Then, the new terms ± h q | K 1 ( t ) K M | with | K 1 ( t ) K M | = ( K 1 ( t ) K M ) and h q > 0 are added up in (75), which results in
V ˙ ( σ ) h p | σ | h q | K 1 ( t ) K M | + | K 1 ( t ) K M | ( h q k sign ( | σ | μ ) + | σ | 3 2 ) h p 2 2 | σ | h q 2 2 | K 1 ( t ) K M | + | K 1 ( t ) K M | ( h q k sign ( | σ | μ ) + | σ | 3 2 ) min h p 2 , h q 2 | σ | 2 + | K 1 ( t ) K M | 2 + | K 1 ( t ) K M | ( h q k sign ( | σ | μ ) + | σ | 3 2 ) .
The previous equation can be rewritten in the form of (4) from Lemma 4 with h = min { h p 2 , h q 2 } , L = 0.5 and W = | K 1 ( t ) K M | ( h q k sign ( | σ | μ ) + | σ | 3 2 ) > 0 . Therefore, referring to Definition 3, there is a small neighborhood of zero described by Γ > 0 such that the achievement of the sliding mode in a practical finite-time can be divided into two phases: (1) the reaching phase when | σ | > Γ and (2) the controller convergence when | σ | Γ at time
T k V 1 2 ( σ 0 ) 0.5 min h p 2 , h q 2 Y 0 .
FxtESO-Based ASMC Stability
The previously proven fixed-time convergence of FxtESO and the practical finite-time convergence of the ASMC permit the separation principle to be fulfilled. Hence, both algorithms, controller and observer, can be designed independently [55]. Then, to avoid ambiguities, the gains of the extended observer, and thus its dynamics, are chosen such that the estimated states converge to the actual state at time T f sufficiently faster than controller convergence T k , i.e., T f T k . The aforementioned is guaranteed as long as the total disturbance affecting the system remains in bounds | ξ 3 | δ max . Thus, strong enough uncertainties and external disturbances, or long delays in the input or output, still can make the system unstable. The proof is complete. □

5. Results

The performance of the proposed observer-based control strategy was verified through computer simulations in the Matlab/Simulink environment. The 60 s virtual experiment was run at a frequency of 100 Hz and consisted of a quadrotor UAV executing a stationary flight followed by a spiral trajectory defined by
x d ( t ) = 0 if 0 t 5 sin 0.3 ( t 5 ) d t if 5 < t 60 , y d ( t ) = 0 if 0 t 5 cos 0.3 ( t 5 ) d t if 5 < t 60 , z d ( t ) = 2 if 0 t 10 2 0.1 ( t 10 ) if 10 < t 60 , ψ d ( t ) = 0 if 0 t 5 0.3 ( t 5 ) if 5 < t 60 .
The rotorcraft was completely exposed to severe turbulent wind fields described by the Von Karman model with parameters v ¯ w x = 8 m / s , v ¯ w y = 10 m / s , v ¯ w z = 6 m / s , W 20 = 45 knots = 23.15 m / s , and ρ = 1.2 Kg / m 3 , where the resultant disturbances can be appreciated in Figure 3. The quadrotor was programmed with initial conditions d ( 0 ) = [ 0 , 0 , 0 ] T m, Θ ( 0 ) = [ 0 , 0 , 0 ] T rad, and specifications J x x = 0.0411 Kgm 2 , J y y = 0.0478 Kgm 2 , J z z = 0.0599 Kgm 2 , l = 0.225 m , K T = 1.356 × 10 5 N / ( rad 2 / s 2 ) , K τ = 2.036 × 10 7 Nm / ( rad 2 / s 2 ) , B x = B y = 0.3 , B z = 0.5 , A x = A y = 0.1 m 2 , and A z = 0.2 m 2 . Moreover, the vehicle was set with a mass of m = 2 Kg. However, an approximate mass of m u = 1.75 Kg was included in the observer-based controller to show its effectiveness under model uncertainties. Additionally, white Gaussian noise signals of a 100 Hz frequency were applied to the outputs of the quadrotor. These signals are illustrated in Figure 4 and were established with mean 0.005 and standard deviation 0.001 for position measurements and mean 0.001 and standard deviation 0.001 for attitude measurements. Finally, the maximum motor thrust input was bounded to 30 T h < 0 N for design purposes. The chosen tuning parameters of the ASMC were λ x = λ y = 1 , λ z = 2 , K 2 x = K 2 y = 1 , K 2 z = 2 , k x = 0.2 , k y = 0.3 ,   k z = 0.5 , μ x = μ y = μ z = 0.05 , k m i n x = k m i n y = k m i n z = 0.01 for the position subsystem, and λ ϕ = λ θ = 8 , λ ψ = 2 , K 2 ϕ = K 2 θ = 6 , K 2 ψ = 0.1 , k ϕ = k θ = 0.5 , k ψ = 0.3 , μ ϕ = μ θ = μ ψ = 0.05 , k m i n ϕ = k m i n θ = 0.1 , and k m i n ψ = 0.01 for the attitude subsystem. To demonstrate that FxtESO with gains κ 1 x , y , z , = κ 1 ϕ , θ , ψ = ε 1 x , y , z = ε 1 ϕ , θ , ψ = 10 , κ 2 x , y , z , = κ 2 ϕ , θ , ψ = ε 2 x , y , z = ε 2 ϕ , θ , ψ = 25 , κ 3 x , y , z , = κ 3 ϕ , θ , ψ = ε 3 x , y , z = ε 3 ϕ , θ , ψ = 25 , α 1 x , y , z = α 1 ϕ , θ , ψ = 0.75 , α 2 x , y , z = α 2 ϕ , θ , ψ = 0.4 , α 3 x , y , z = α 3 ϕ , θ , ψ = 0.1 , β 1 x , y , z = β 1 ϕ , θ , ψ = 1.2 , β 2 x , y , z = β 2 ϕ , θ , ψ = 1.4 , β 3 x , y , z = β 3 ϕ , θ , ψ = 1.6 , χ x , y , z = χ ϕ , θ , ψ = 0.1 recovers the target system, its initial conditions were set as d ^ ( 0 ) = [ 0 , 0 , 0.1 ] T m, and Θ ^ ( 0 ) = [ 0.05 , 0.05 , 0.1 ] T rad.

5.1. Disturbance Compensation

A comparison involving the FxtESO-based ASMC with (w/DC) and without (w/o DC) disturbance compensation was given to demonstrate the advantages of its implementation on the performance of the UAV. The applied turbulent wind gusts reached peak velocities of max ( | v w x | ) = 15.83 m / s , max ( | v w y | ) = 16 m / s , and max ( | v w z | ) = 11.10 m / s which implied maximum drag forces of max ( | δ x | ) = 4.58 N , max ( | δ y | ) = 5.03 N , and max ( | δ z | ) = 7.51 N . The latter is equivalent to 23.36 % , 25.65 % , and 38.30 % of the total weight of the aerial vehicle. The execution of the trajectory of the rotorcraft, from a top view, is portrayed in Figure 5a. It is evident that the approach including the disturbance compensation action offered a superior tracking performance with slighter oscillations and deviations from the desired position. In addition, Figure 5b,c demonstrate two important differences among executions in the x and y axes. First, the disturbance rejection-based control system presented a more robust reaction to perturbations within the first seconds of the simulation, during the hover flight execution. Second, we note that it took longer for the feedback regulation-based sliding mode strategy to become close to the desired state due to the existence of both external disruptions and model uncertainties. Moving on, the vertical performance of the UAV is given in Figure 5d. The vehicle managed through the ASMC without disturbance compensation struggled in achieving the reference altitude at t 10 since the turbulent wind fields were continuously pushing the quadrotor down. After t 10 , the discrepancies between both controllers were less perceptible. However, the zoomed-in section of the graph denotes slight oscillations on the obtained trajectory without disturbance compensation.
The heading performance of the UAV is displayed in Figure 6a. In this particular case, both controllers delivered a similar execution throughout the simulation except at interval 40 t 50 where the disturbance compensation-free ASMC response presented slight oscillations. The executed roll and pitch rotations of the rotorcraft are shown in Figure 6b. During the time interval 0 < t < 20 , it is visible that the outputs from the proposed approach were smoother than those from the controller without disturbance compensation. This is expected since the disturbance compensation action demands faster responses from the system. Furthermore, we notice that once t 20 , unlike FxtESO-ASMC, the classic ASMC strategy began to occasionally provide high-frequency rotations. The latter suggests that the vehicle was exposed to strong disturbances whose magnitude may be close to the maximum boundary it can handle without reaching an unstable state.
The evolution of the adaptive gains for the positioning subsystem is portrayed in Figure 6c. Naturally, the feedback regulation-based controller requires higher magnitudes from the adaptive gains. Moreover, we note that K 1 x and K 1 y exhibited an increasing tendency, whereas the needed effort from K 1 z was significantly larger than the perturbation rejection-based approach. As a consequence, the control input related to the total motor thrust leaned towards presenting chattering, as seen in Figure 6d. Therefore, the implementation of a disturbance compensation action in the control input of the system improves the trajectory execution of the UAV and aids in keeping a lower value from the adaptive gains of the controller while attenuating the chattering phenomenon in exchange for a regularly larger control input magnitude.

5.2. Comparison Study

The operation of the proposed FxtESO-based control strategy was compared with the nonlinear extended state observer-based adaptive sliding mode controller (NESO-ASMC) from [56] described by NESO:
ξ ^ ˙ 1 = ξ ^ 2 Q 1 f ( ξ ˜ ) ,
ξ ^ ˙ 2 = ξ ^ 3 Q 2 f ( ξ ˜ ) ,
ξ ^ ˙ 3 = Q 3 f ( ξ ˜ ) ,
with
f ( ξ ˜ ) = | ξ 1 ξ ^ 1 | ι sign ( ξ 1 ξ ^ 1 ) if | ξ 1 ξ ^ 1 | > Λ ξ 1 ξ ^ 1 Λ 1 ι if | ξ 1 ξ ^ 1 | Λ
and the auxiliary control strategy
ν = ς 1 sig 1 / 2 ( σ ) ς 2 2 sign ( σ ) ,
where ς 2 = ς 3 ς 1 , ς 3 > 0 , and the adaptability law is
ς ˙ 1 = ς 4 sign ( | σ | μ ς ) ς 1 > ς min ς min ς 1 ς min .
For comparative analysis, the NESO-based ASMC parameters were chosen as Q 1 q = 20 , Q 2 q = 247 , Q 3 q = 1500 , ι q = 0.7 , and Λ q = 0.05 for q = { x , y , z , ϕ , θ , ψ } , ς 3 x = ς 3 y = 0.01 , ς 3 z = 0.2 , ς 4 x = ς 4 y = 0.5 , ς 4 z = 3 , ς min x = ς min y = 0.01 , ς min z = 0.1 , μ ς x = μ ς y = μ ς z = 0.05 , ς 3 ϕ = ς 3 θ = 0.3 , ς 3 ψ = 1 , ς 4 ϕ = ς 4 θ = ς 4 ψ = 1 , ς min ϕ = ς min θ = ς min ψ = 0.1 , and μ ς ϕ = μ ς θ = μ ς ψ = 0.05 .
The evolution of the disturbance estimation errors from both observer-based methods is displayed in Figure 7. Although their estimation performance under sensor noise conditions is accurate, we can notice that the NESO provided a larger peaking phenomenon than the fixed-time observer (see ξ ˜ 3 z ) due to the initial estimation error ξ ˜ z = 0.1 and the UAV mass uncertainty. As a consequence, Figure 8a shows that the FxtESO-based ASMC granted a more accurate trajectory tracking response with lower-amplitude oscillations at some points of the operation. Additionally, the trajectory execution in the z-axis from Figure 8b clearly portrays the effects of the peaking phenomenon on the NESO-based ASMC since the quadrotor was not able to take off as fast as the vehicle driven by the FxtESO-ASMC. Furthermore, the displayed trajectories along the x and y from Figure 8c,d demonstrate that the proposed method offered stronger robustness specifically during the hover flight stage of the trajectory where the NESO-ASMC presented deviations between 1 and 3 m. The effect of prior errors in the attitude subsystem is given in Figure 8e. We notice that, unlike in the proposed approach, the influence of the total disturbance in the system caused the quadrotor managed under the NESO-ASMC execute high-frequency rotations within the first 10 s of the simulation. On the other hand, Figure 8f illustrates that both observer-based controllers provided an effective tracking of the desired heading.
The state estimation performance of the employed observers for the yaw output is given in Figure 9a. Despite the presence of sensor noise and model uncertainties, both methods achieved the real state within 0.5 s. However, the proposed FxtESO offered smoother and more straightforward state tracking since it did not present any overshoot. The evolution of the positioning subsystem adaptive gains is illustrated in Figure 9b. At the start of the simulation, the FxtESO-ASMC approach required approximately 400 % less magnitudes of K 1 x , K 1 y , and K 1 z than the NESO-ASMC scheme. Furthermore, at time t > 20 , the adjustments of parameters K 1 x and K 1 y from both observer-based controllers gradually settled down and stayed at a range of values denoted by [0.5,2). Regarding adaptive gain K 1 z at t > 20 , it is easy to note that it had a more active evolution under the NESO-ASMC algorithm, whereas the same parameter from the FxtESO-ASMC kept within the same magnitude interval as K 1 x and K 1 y . The applied control inputs for the quadrotor are portrayed in Figure 9c,d, respectively. The thrust provided by the NESO-ASMC at time 0 < t 5 shows the effect of the peaking phenomenon on the propulsion of the system. While the signal from the FxtESO-ASMC showed no peaking. Moreover, it allowed operation of the rotorcraft with an average thrust of 22 N with few peaks between 26 and 28 N, which is suitable for a 2 Kg vehicle. Finally, the advantage of the proposed approach with respect to the executed torques is evident. The fixed-time observer-based controller supplied rotational control inputs within the range of ± 0.1 Nm, which is available considering the determined characteristics of the UAV. On the other hand, the torques given by the NESO-based algorithm remained within the range of ± 0.3 Nm, including a few short-time spikes that reached ± 0.5 Nm.
A quantitative analysis for prior results is presented in Table 1. Performance indicators such as the mean squared error,
M S E = 1 N i = 1 N e i 2 , for i = { 1 , 2 , , N } ,
the root mean squared error,
R M S E = 1 N i = 1 N e i 2 , for i = { 1 , 2 , , N } ,
the integral absolute error,
I A E = 0 T | e ( t ) | d t ,
the integral squared error,
I S E = 0 T e ( t ) 2 d t ,
and the integral time absolute error,
I T A E = 0 T | e ( t ) | t d t ,
were employed to evaluate and corroborate the advantages of the proposed method. Statistically, it can be concluded that both observer-based strategies provided the same position estimation performance since their discrepancies in this aspect were not significant. However, the computed data show the superiority of the FxtESO-ASMC approach in the disturbance estimation and position error areas. Furthermore, the information related to the norm of control suggests that the adaptive fixed-time observer-based controller needed less control effort than the similar, NESO-based controller.

6. Conclusions

A fixed-time extended state observer-based adaptive sliding mode control was designed for driving a quadrotor UAV subject to severe-intensity turbulent wind gusts while tracking a spiral trajectory. The external disturbances, denoted by the Von Karman wind model, and the effects of model uncertainties regarding the angular dynamics and the mass of the rotorcraft were rejected via fixed-time augmented state estimation and feedback control loop compensation, regardless of the initial conditions of the system. Furthermore, the adaptive gain from the controller allowed the usage of the necessary control input, preserving robustness. Moreover, the closed-loop system stability was guaranteed by means of Lyapunov theory. Finally, simulations in Simulink/Matlab demonstrated the advantages of disturbance compensation in improving the performance of the UAV and the effectiveness of the proposed control scheme by making comparisons against a nonlinear extended state observer-based adaptive controller.

Author Contributions

Conceptualization, A.M.-M., H.C. and H.W.; methodology, A.M.-M. and H.C.; software, A.M.-M.; validation, H.C. and H.W.; writing—original draft preparation, A.M.-M.; writing—review and editing, A.M.-M., H.C. and H.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data is contained within the article.

Acknowledgments

We thank Multi-Robot System Laboratory of the Tecnologico de Monterrey, and the CONACyT scholarship for support in the development of this project.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ASMCAdaptive sliding mode control
ESOExtended state observer
FxtESOFixed-time extended state observer
NESONonlinear extended state observer
UAVUnamnned Aerial Vehicle

References

  1. Najm, A.A.; Ibraheem, I.K. Nonlinear PID controller design for a 6-DOF UAV quadrotor system. Eng. Sci. Technol. Int. J. 2019, 22, 1087–1097. [Google Scholar] [CrossRef]
  2. Zhou, L.; Zhang, J.; She, H.; Jin, H. Quadrotor UAV flight control via a novel saturation integral backstepping controller. Automatika 2019, 60, 193–206. [Google Scholar] [CrossRef]
  3. Labbadi, M.; Cherkaoui, M. Robust adaptive backstepping fast terminal sliding mode controller for uncertain quadrotor UAV. Aerosp. Sci. Technol. 2019, 93, 105306. [Google Scholar] [CrossRef]
  4. Labbadi, M.; Cherkaoui, M.; El houm, Y.; Guisser, M. Modeling and robust integral sliding mode control for a quadrotor unmanned aerial vehicle. In Proceedings of the 2018 IEEE 6th International Renewable and Sustainable Energy Conference (IRSEC), Rabat, Morocco, 5–8 December 2018; 2018; pp. 1–6. [Google Scholar]
  5. Nekoukar, V.; Dehkordi, N.M. Robust path tracking of a quadrotor using adaptive fuzzy terminal sliding mode control. Control. Eng. Pract. 2021, 110, 104763. [Google Scholar] [CrossRef]
  6. Li, S.; Yang, J.; Chen, W.H.; Chen, X. Disturbance Observer-Based Control: Methods and Applications; CRC Press: Boca Raton, FL, USA, 2014. [Google Scholar]
  7. Miranda-Moya, A.; Castañeda, H.; Wang, H. Turbulent Wind Gusts Estimation and Compensation via High-Gain Extended Observer-based Adaptive Sliding Mode for a Quadrotor UAV. In Proceedings of the IEEE 2021 International Conference on Unmanned Aircraft Systems (ICUAS), Athens, Greece, 15–18 June 2021; pp. 476–481. [Google Scholar]
  8. Davoudi, B.; Taheri, E.; Duraisamy, K.; Jayaraman, B.; Kolmanovsky, I. Quad-rotor flight simulation in realistic atmospheric conditions. AIAA J. 2020, 58, 1992–2004. [Google Scholar] [CrossRef]
  9. Gao, Z.; Guo, G. Command-filtered fixed-time trajectory tracking control of surface vehicles based on a disturbance observer. Int. J. Robust Nonlinear Control. 2019, 29, 4348–4365. [Google Scholar] [CrossRef]
  10. Lyu, X.; Zhou, J.; Gu, H.; Li, Z.; Shen, S.; Zhang, F. Disturbance Observer Based Hovering Control of Quadrotor Tail-sitter VTOL UAVs Using H Synthesis. IEEE Robot. Autom. Lett. 2018, 3, 2910–2917. [Google Scholar] [CrossRef]
  11. Wu, K.; Fan, B.; Zhang, X. Trajectory following control of UAVs with wind disturbance. In Proceedings of the 2017 36th Chinese Control Conference (CCC), Dalian, China, 26–28 June 2017; pp. 4993–4997. [Google Scholar]
  12. Yu, G.; Chen, Y.; Chen, Z.; Wu, H.; Cheng, L. Design of terminal sliding mode controller for a quadrotor UAV with disturbance observer. In Proceedings of the 2020 IEEE 39th Chinese Control Conference (CCC), Shenyang, China, 27–29 July 2020; pp. 2072–2077. [Google Scholar]
  13. Qin, L.; Zhou, W.; Li, L.; Jiang, W. Active disturbance rejection control system design for quadrotor. In Proceedings of the 2017 36th Chinese Control Conference (CCC), Dalian, China, 26–28 June 2017; pp. 6530–6534. [Google Scholar]
  14. Wu, C.; Yan, J.; Lin, H.; Wu, X.; Xiao, B. Fixed-time disturbance observer-based chattering-free sliding mode attitude tracking control of aircraft with sensor noises. Aerosp. Sci. Technol. 2021, 111, 106565. [Google Scholar] [CrossRef]
  15. Wu, S.; Chen, L.; Zhang, D.; Chen, J.; Shao, X. Disturbance observer based fixed time sliding mode control for spacecraft proximity operations with coupled dynamics. Adv. Space Res. 2020, 66, 2179–2193. [Google Scholar] [CrossRef]
  16. Wang, X.; Guo, J.; Tang, S.; Qi, S. Fixed-time disturbance observer based fixed-time back-stepping control for an air-breathing hypersonic vehicle. ISA Trans. 2019, 88, 233–245. [Google Scholar] [CrossRef]
  17. Smith, J.; Su, J.; Liu, C.; Chen, W.H. Disturbance observer based control with anti-windup applied to a small fixed wing UAV for disturbance rejection. J. Intell. Robot. Syst. 2017, 88, 329–346. [Google Scholar] [CrossRef]
  18. MIL-STD-1797A; Department of Defense Handbook: Flying Qualities of Piloted Aircraft. University Reprints 2017; Department of Defense: Washington, DC, USA, 1995.
  19. Rodriguez-Mata, A.E.; Gonzalez-Hernandez, I.; Rangel-Peraza, J.G.; Salazar, S.; Leal, R.L. Wind-gust compensation algorithm based on high-gain residual observer to control a quadrotor aircraft: Real-time verification task at fixed point. Int. J. Control. Autom. Syst. 2018, 16, 856–866. [Google Scholar] [CrossRef]
  20. Zhang, Z.; Wang, F.; Guo, Y.; Hua, C. Multivariable sliding mode backstepping controller design for quadrotor UAV based on disturbance observer. Sci. China Inf. Sci. 2018, 61, 112207. [Google Scholar] [CrossRef]
  21. Aboudonia, A.; Rashad, R.; El-Badawy, A. Time domain disturbance observer based control of a quadrotor unmanned aerial vehicle. In Proceedings of the 2015 XXV International Conference on Information, Communication and Automation Technologies (ICAT), Sarajevo, Bosnia and Herzegovina, 29–31 October 2015; pp. 1–6. [Google Scholar]
  22. Aboudonia, A.; El-Badawy, A.; Rashad, R. Disturbance observer-based feedback linearization control of an unmanned quadrotor helicopter. Proc. Inst. Mech. Eng. Part J. Syst. Control. Eng. 2016, 230, 877–891. [Google Scholar] [CrossRef]
  23. Aboudonia, A.; Rashad, R.; El-Badawy, A. Composite hierarchical anti-disturbance control of a quadrotor UAV in the presence of matched and mismatched disturbances. J. Intell. Robot. Syst. 2018, 90, 201–216. [Google Scholar] [CrossRef]
  24. Han, J. From PID to active disturbance rejection control. IEEE Trans. Ind. Electron. 2009, 56, 900–906. [Google Scholar] [CrossRef]
  25. Wang, H.; Li, N.; Wang, Y.; Su, B. Backstepping sliding mode trajectory tracking via extended state observer for quadrotors with wind disturbance. Int. J. Control. Autom. Syst. 2021, 19, 3273–3284. [Google Scholar] [CrossRef]
  26. Shao, X.; Liu, J.; Cao, H.; Shen, C.; Wang, H. Robust dynamic surface trajectory tracking control for a quadrotor UAV via extended state observer. Int. J. Robust Nonlinear Control. 2018, 28, 2700–2719. [Google Scholar] [CrossRef]
  27. Zhiyuan, C.; Yanyang, L.; Yanhua, S.; Hongyu, C.; Bin, W.; Mingqi, H.; Rao, Y. Fuzzy sliding mode control for rotorcraft aerial manipulator with extended state observer. In Proceedings of the 2020 IEEE: Chinese Automation Congress (CAC), Shanghai, China, 6–8 November 2020; pp. 1710–1714. [Google Scholar]
  28. Khalil, H.K. Extended high-gain observers as disturbance estimators. Sice J. Control. Meas. Syst. Integr. 2017, 10, 125–134. [Google Scholar] [CrossRef]
  29. Liu, J.; Sun, M.; Chen, Z.; Sun, Q. Super-twisting sliding mode control for aircraft at high angle of attack based on finite-time extended state observer. Nonlinear Dyn. 2020, 99, 2785–2799. [Google Scholar] [CrossRef]
  30. Polyakov, A. Nonlinear feedback design for fixed-time stabilization of linear control systems. IEEE Trans. Autom. Control. 2011, 57, 2106–2110. [Google Scholar] [CrossRef]
  31. Biao, T.; Xingling, S.; Wei, Y.; Wendong, Z. Fixed time output feedback containment for uncertain nonlinear multiagent systems with switching communication topologies. ISA Trans. 2021, 111, 82–95. [Google Scholar] [CrossRef] [PubMed]
  32. Hu, D.; Zhang, S.; Zou, A.M. Velocity-free fixed-time attitude cooperative control for spacecraft formations under directed graphs. Int. J. Robust Nonlinear Control. 2021, 31, 2905–2927. [Google Scholar] [CrossRef]
  33. Cui, L.; Jin, N.; Chang, S.; Zuo, Z.; Zhao, Z. Fixed-time ESO based fixed-time integral terminal sliding mode controller design for a missile. ISA Trans. 2021, 125, 237–251. [Google Scholar] [CrossRef]
  34. Jin, R.; Rocco, P.; Geng, Y. Observer-based fixed-time tracking control for space robots in task space. Acta Astronaut. 2021, 184, 35–45. [Google Scholar] [CrossRef]
  35. Mechali, O.; Xu, L.; Huang, Y.; Shi, M.; Xie, X. Observer-based fixed-time continuous nonsingular terminal sliding mode control of quadrotor aircraft under uncertainties and disturbances for robust trajectory tracking: Theory and experiment. Control. Eng. Pract. 2021, 111, 104806. [Google Scholar] [CrossRef]
  36. Li, B.; Ban, H.; Gong, W.; Xiao, B. Extended state observer-based finite-time dynamic surface control for trajectory tracking of a quadrotor unmanned aerial vehicle. Trans. Inst. Meas. Control 2020, 42, 2956–2968. [Google Scholar] [CrossRef]
  37. Liu, J.; Sun, M.; Chen, Z.; Sun, Q. Output feedback control for aircraft at high angle of attack based upon fixed-time extended state observer. Aerosp. Sci. Technol. 2019, 95, 105468. [Google Scholar] [CrossRef]
  38. Khalil, H.K. Nonlinear Systems, 3rd ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2002; Volume 115. [Google Scholar]
  39. Fridman, L.; Barbot, J.P.; Plestan, F. Recent trends in sliding mode control. In IET Control, Robotics 552 and Sensors Series Volume 102; The Institution of Engineering and Technology: London, UK, 2016. [Google Scholar]
  40. Zhu, Z.; Xia, Y.; Fu, M. Attitude stabilization of rigid spacecraft with finite-time convergence. Int. J. Robust Nonlinear Control. 2011, 21, 686–702. [Google Scholar] [CrossRef]
  41. Zhang, J.; Yu, S.; Yan, Y. Fixed-time output feedback trajectory tracking control of marine surface vessels subject to unknown external disturbances and uncertainties. ISA Trans. 2019, 93, 145–155. [Google Scholar] [CrossRef]
  42. Polyakov, A.; Fridman, L. Stability notions and Lyapunov functions for sliding mode control systems. J. Frankl. Inst. 2014, 351, 1831–1865. [Google Scholar] [CrossRef]
  43. Perruquetti, W.; Floquet, T.; Moulay, E. Finite-time observers: Application to secure communication. IEEE Trans. Autom. Control. 2008, 53, 356–360. [Google Scholar] [CrossRef]
  44. Bhat, S.P.; Bernstein, D.S. Geometric homogeneity with applications to finite-time stability. Math. Control. Signals Syst. 2005, 17, 101–127. [Google Scholar] [CrossRef]
  45. Wang, H.; Chen, B.; Lin, C.; Sun, Y.; Wang, F. Adaptive finite-time control for a class of uncertain high-order non-linear systems based on fuzzy approximation. IET Control. Theory Appl. 2017, 11, 677–684. [Google Scholar] [CrossRef]
  46. Castañeda, H.; Gordillo, J.L. Spatial Modeling and Robust Flight Control Based on Adaptive Sliding Mode Approach for a Quadrotor MAV. J. Intell. Robot. Syst. 2019, 93, 101–111. [Google Scholar] [CrossRef]
  47. Gage, S. Creating a unified graphical wind turbulence model from multiple specifications. In Proceedings of the AIAA Modeling and Simulation Technologies Conference and Exhibit, Austin, TX, USA, 1–14 August 2003; p. 5529. [Google Scholar]
  48. Zhang, J.; Yu, S.; Yan, Y.; Wu, D. Fixed-time output feedback sliding mode tracking control of marine surface vessels under actuator faults with disturbance cancellation. Appl. Ocean. Res. 2020, 104, 102378. [Google Scholar] [CrossRef]
  49. Castañeda, H.; Rodriguez, J.; Gordillo, J.L. Continuous and smooth differentiator based on adaptive sliding mode control for a quad-rotor MAV. Asian J. Control. 2019, 23, 661–672. [Google Scholar] [CrossRef]
  50. Zhang, L.; Wei, C.; Wu, R.; Cui, N. Fixed-time extended state observer based non-singular fast terminal sliding mode control for a VTVL reusable launch vehicle. Aerosp. Sci. Technol. 2018, 82, 70–79. [Google Scholar] [CrossRef]
  51. Basin, M.; Yu, P.; Shtessel, Y. Finite-and fixed-time differentiators utilising HOSM techniques. IET Control. Theory Appl. 2017, 11, 1144–1152. [Google Scholar] [CrossRef]
  52. Zhang, J.; Yu, S.; Yan, Y. Fixed-time extended state observer-based trajectory tracking and point stabilization control for marine surface vessels with uncertainties and disturbances. Ocean. Eng. 2019, 186, 106109. [Google Scholar] [CrossRef]
  53. Tian, B.; Zuo, Z.; Yan, X.; Wang, H. A fixed-time output feedback control scheme for double integrator systems. Automatica 2017, 80, 17–24. [Google Scholar] [CrossRef]
  54. Plestan, F.; Shtessel, Y.; Bregeault, V.; Poznyak, A. New methodologies for adaptive sliding mode control. Int. J. Control. 2010, 83, 1907–1919. [Google Scholar] [CrossRef]
  55. Khalil, H.K. High-Gain Observers in Nonlinear Feedback Control; SIAM: Philadelphia, PA, USA, 2017. [Google Scholar]
  56. Castañeda, H.; Salas-Peña, O.S.; de León-Morales, J. Extended observer based on adaptive second order sliding mode control for a fixed wing UAV. ISA Trans. 2017, 66, 226–232. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Quadrotor UAV. Reference frames.
Figure 1. Quadrotor UAV. Reference frames.
Drones 07 00700 g001
Figure 2. FxtESO–based ASMC scheme.
Figure 2. FxtESO–based ASMC scheme.
Drones 07 00700 g002
Figure 3. Resultant external disturbances described by the Von Karman wind model.
Figure 3. Resultant external disturbances described by the Von Karman wind model.
Drones 07 00700 g003
Figure 4. Applied sensor noise to the position and attitude measurements.
Figure 4. Applied sensor noise to the position and attitude measurements.
Drones 07 00700 g004
Figure 5. Disturbance compensation analysis. Trajectory tracking performance of (a) Top view, (b) y response, (c) x coordinate and (d) z axis evolution.
Figure 5. Disturbance compensation analysis. Trajectory tracking performance of (a) Top view, (b) y response, (c) x coordinate and (d) z axis evolution.
Drones 07 00700 g005
Figure 6. Disturbance compensation analysis. Attitude performance and control input evolution. (a) Yaw angle (b) Roll and pitch responses, (c) Adaptive gains, (d) Thrust.
Figure 6. Disturbance compensation analysis. Attitude performance and control input evolution. (a) Yaw angle (b) Roll and pitch responses, (c) Adaptive gains, (d) Thrust.
Drones 07 00700 g006
Figure 7. ESO–based controller comparison. Disturbance estimation error.
Figure 7. ESO–based controller comparison. Disturbance estimation error.
Drones 07 00700 g007
Figure 8. FxtESO and NESO–based adaptive controller comparison. Trajectory tracking: (a) top view, (b) z axis, (c) x response, (d) y behavior, (e) Roll and pitch estimation, (f) yaw estimation.
Figure 8. FxtESO and NESO–based adaptive controller comparison. Trajectory tracking: (a) top view, (b) z axis, (c) x response, (d) y behavior, (e) Roll and pitch estimation, (f) yaw estimation.
Drones 07 00700 g008
Figure 9. Quadrotor UAV trajectory. FxtESO and NESO–based adaptive controller comparison. (a) Yaw convergence, (b) Adaptive gains, (c) Required thrust, (d) Generated torques.
Figure 9. Quadrotor UAV trajectory. FxtESO and NESO–based adaptive controller comparison. (a) Yaw convergence, (b) Adaptive gains, (c) Required thrust, (d) Generated torques.
Drones 07 00700 g009
Table 1. Quantitative analysis of the results.
Table 1. Quantitative analysis of the results.
Disturbance Estimation
MSERMSEIAEISEITAE
x-axis
FxtESO-ASMC0.04690.21669.93902.8160302.7
NESO-ASMC0.05860.242111.003.5190305.2
y-axis
FxtESO-ASMC0.08750.295813.515.25377.2
NESO-ASMC0.10590.325414.346.3550364.6
z-axis
FxtESO-ASMC0.46980.685429.2528.19729.3
NESO-ASMC1.09881.048233.3765.94699.5
Position Estimation
MSERMSEIAEISEITAE
x-axis
FxtESO-ASMC1.41 × 10 6 1.19 × 10 3 0.05698.50 × 10 5 1.7110
NESO-ASMC1.38 × 10 6 1.17 × 10 3 0.05628.30 × 10 5 1.6484
y-axis
FxtESO-ASMC1.43 × 10 6 1.19 × 10 3 0.05728.62 × 10 5 1.7180
NESO-ASMC1.40 × 10 6 1.18 × 10 3 0.05658.43 × 10 5 1.69
z-axis
FxtESO-ASMC5.43 × 10 6 2.33 × 10 3 0.06553.26 × 10 4 1.799
NESO-ASMC6.95 × 10 6 2.63 × 10 3 0.06974.17 × 10 4 1.768
MSERMSEIAEISEITAE
x-axis
FxtESO-ASMC0.00750.08664.19900.4472104.3
NESO-ASMC0.05660.23796.08703.3990114.7
FxtESO-ASMC w/o DC0.04130.20329.1862.481234.9
y-axis
FxtESO-ASMC0.01110.10543.59100.664077.02
NESO-ASMC0.29250.540811.6517.55117.30
FxtESO-ASMC w/o DC0.05360.23158.79703.214153.8
z-axis
FxtESO-ASMC0.04650.21563.10202.791034.13
NESO-ASMC0.14990.38726.08708.994039.97
FxtESO-ASMC w/o DC0.09170.30288.04205.5077.95
Quadrotor UAV. Norms of Control
T h τ
FxtESO-ASMC1717.70395.9815
NESO-ASMC1721.801041.7
FxtESO-ASMC w/o DC1716.80785.2556
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Miranda-Moya, A.; Castañeda, H.; Wang, H. Fixed-Time Extended Observer-Based Adaptive Sliding Mode Control for a Quadrotor UAV under Severe Turbulent Wind. Drones 2023, 7, 700. https://doi.org/10.3390/drones7120700

AMA Style

Miranda-Moya A, Castañeda H, Wang H. Fixed-Time Extended Observer-Based Adaptive Sliding Mode Control for a Quadrotor UAV under Severe Turbulent Wind. Drones. 2023; 7(12):700. https://doi.org/10.3390/drones7120700

Chicago/Turabian Style

Miranda-Moya, Armando, Herman Castañeda, and Hesheng Wang. 2023. "Fixed-Time Extended Observer-Based Adaptive Sliding Mode Control for a Quadrotor UAV under Severe Turbulent Wind" Drones 7, no. 12: 700. https://doi.org/10.3390/drones7120700

APA Style

Miranda-Moya, A., Castañeda, H., & Wang, H. (2023). Fixed-Time Extended Observer-Based Adaptive Sliding Mode Control for a Quadrotor UAV under Severe Turbulent Wind. Drones, 7(12), 700. https://doi.org/10.3390/drones7120700

Article Metrics

Back to TopTop