Next Article in Journal
Parametric Study of Structured UTM Separation Recommendations with Physics-Based Monte Carlo Distribution for Collision Risk Model
Previous Article in Journal
Learning to Propose and Refine for Accurate and Robust Tracking via an Alignment Convolution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unknown Input Observer-Based Fixed-Time Trajectory Tracking Control for QUAV with Actuator Saturation and Faults

School of Electrical Engineering, Hebei University of Science and Technology, Shijiazhuang 050018, China
*
Author to whom correspondence should be addressed.
Drones 2023, 7(6), 344; https://doi.org/10.3390/drones7060344
Submission received: 14 April 2023 / Revised: 19 May 2023 / Accepted: 23 May 2023 / Published: 25 May 2023

Abstract

:
The trajectory tracking control problem of a quadrotor unmanned aerial vehicle (QUAV) subject to external disturbances, inertia uncertainties, actuator faults, and input saturation is addressed in this paper. In contrast with previous works, input saturation herein refers to rotor speed saturation rather than thrust and torque saturation. First, the control system is decoupled into translational and rotational subsystems. Then, for both subsystems, two novel fixed-time unknown input observers (UIO) based on disturbance filtering are developed to estimate the lumped disturbance rapidly and precisely without awareness of the boundary of disturbances. Furthermore, fixed-time tracking controllers for translational and rotational subsystems are proposed based on the estimation values provided by UIO to stabilize tracking errors into a small region in fixed time regardless of the initial values. The theoretical analysis based on the Lyapunov method is presented to demonstrate the stability. Finally, the simulation results show that the proposed control method is effective. The comparison simulation is carried out to validate superiority of the proposed observer and its advantage can be summed up as: (1) the upper bound of the disturbance or its derivative is not needed; (2) the estimation results are smoother and the observation precision is higher due to the absence of sign function; (3) the mutant disturbance can be also estimated quickly and precisely.

1. Introduction

In recent years, research on quadrotor unmanned aerial vehicles (QUAVs) has been widely conducted. The interest is sparked by the broad range of tasks for which it can be applied, including aerial photography, inspection monitoring, and agricultural plant protection, all while enjoying the benefits of small size, compact structure, and flexible maneuverability [1,2,3]. However, the QUAV is a multivariable, strongly coupled, nonlinear under-actuated system with four inputs but six degrees of freedom, which makes the control design and stability analysis more complicated [4,5]. As a result, trajectory tracking control of QUAV has become a complicated matter that has been extensively studied.
Traditional linear control strategies, such as PID [6] and LQR [7], are only applicable to linearization models, and control performance for nonlinear systems is hard to guarantee. As a result, numerous nonlinear control schemes, including but not limited to feedback linearization control [8], backstepping control [9,10], and sliding mode control (SMC) [11,12,13], have been developed for the QUAV in order to mitigate the influence of nonlinear characteristics. Moreover, the speed of convergence is regarded as an important parameter for the QUAV’s control systems. Therefore, control systems based on finite-time control have gained significant attention [2,14,15,16,17]. For example, in [14], the authors proposed a fast terminal sliding mode controller (FTSMC) to achieve trajectory tracking control of the QUAV, and the tracking error for desired altitude and attitude was guaranteed to converge to zero in finite time. A finite-time nonlinear disturbance observer (NDO) was also presented to estimate external disturbance acting on altitude and attitude channels of the QUAV. An adaptive command-filtered backstepping sliding mode control scheme was proposed to achieve finite-time tracking control of the QUAV system under modeling uncertainties and external disturbances in [16]. However, the convergence time of finite-time control depends on the initial values of the system. Consequently, fixed-time control is proposed by researchers, and the error can converge to zero in fixed time independent of initial values. Xia and Son [18] addressed the fixed-time control problem of autonomous ship landing operations of vertical take-off and landing UAVs subject to external disturbances. An anti-saturation adaptive fault-tolerant control scheme with fixed-time prescribed performance for the longitudinal model of fixed wing UAV under actuator faults, control input saturation, angle of attack asymmetric constraint, and uncertainties was presented in [19]. Ref. [20] proposed a fixed-time adaptive fast super-twisting disturbance observer to estimate the unknown external disturbance, and a fixed-time controller was designed using a universal barrier Lyapunov function to satisfy asymmetric tracking error constraints.
The external disturbances, inertia uncertainties, and actuator faults often need to be considered in the trajectory tracking control for QUAV [4,11,12,16]. The most efficient strategy for dealing with external disturbances, inertia uncertainties, and actuator faults is to compensate for their negative impacts on the system via exact estimation. As a result, numerous disturbance observer-based control approaches, such as those based on the extended state observer (ESO) [21,22], the super-twisting disturbance observer (STDO) [20,23,24], and the sliding mode disturbance observer (SMDO) [25,26], have been developed. Nonetheless, the observers mentioned above all require that disturbance is continuously differentiable and its derivative is bounded or known beforehand. When the disturbance changes suddenly, its derivative is infinite, hence the observation performance cannot be guaranteed in this case. Moreover, the sign function is adopted in ESO and SMDO, so the estimation of disturbance may be not smooth. To overcome this constraint, by introducing a reference auxiliary system, a nonlinear observer was developed to accomplish a high-accuracy estimation for the lumped disturbance in fixed time without being aware of the boundary of disturbances in [27]. However, the design process of this observer is complicated. The control problem of nonlinear robotic systems suffering from actuator failure was investigated in [28]. Based on the state transformation, a finite-time unknown input observer (UIO) was designed to estimate the lumped unknown input. Nevertheless, the observer in [28] can only estimate the lumped unknown input in finite time.
Furthermore, the output capability of the actuator is often limited; hence, input saturation should be frequently considered in trajectory tracking control. In [29], the finite-time trajectory tracking problem for a novel 12-rotor UAV with input saturation was investigated. The authors of [30] addressed the finite-time tracking control problem of the QUAV subject to external disturbances, parametric uncertainties, actuator faults, and input saturation. Cao and Lynch [31] studied the position control of a QUAV with state and input constraints using an inner-outer loop control structure. However, all the works of literature cited above simply supply the saturation values of the control force and torque, which is inexhaustive. The control input signals of the QUAV are determined by the rotor speed of the four motors, and the rotor speed of each motor is different and has an upper limit. In fact, as long as one motor’s rotor speed is saturated, the control input signals related to that motor’s rotor speed are constrained. Therefore, rotor speed saturation is considered in this paper.
Motivated by the preceding discussions, this work focuses on developing a fixed-time control method for the QUAV trajectory tracking control when external disturbances, inertia uncertainties, actuator faults, and input saturation are considered. In comparison to the current investigation results, the major contributions of this work can be summarized as follows:
  • The fixed-time UIOs are designed to estimate the lumped disturbance for translational and rotational subsystems, with the estimation error able to converge to zero in fixed time. The proposed observer has three advantages and can overcome the limitations described in [20,21,22,23,24,25,26]. Firstly, it does not require prior knowledge of the upper bound of disturbance or its derivative. Additionally, the sign function is absent in UIO, which implies that the observation of disturbance is smoother and the estimation precision is higher. Furthermore, it can estimate the mutant disturbance accurately and rapidly.
  • Based on the estimation values provided by fixed-time UIO, fixed-time tracking controllers are proposed for both subsystems to stabilize the tracking error into a small region in fixed time. In contrast with the existing controllers with finite-time stability [2,14,15,16,17], the proposed fixed-time controllers’ convergence time is independent of initial values.
  • External disturbances, inertia uncertainties, actuator faults, and input saturation are all considered in this paper. In [4,11,12,16,28], input saturation is not considered. Although input saturation is considered in [29,30,31], they simply give the saturation values of control input signals. In this paper, input saturation is constraining the rotor speed rather than the control input signals, which is more reasonable.
The rest of this paper can be summarized as follows. In Section 2, the model of the QUAV, the control objectives, and some useful lemmas are given. The control schemes for translational and rotational subsystems are proposed in Section 3. The simulation results show the effectiveness of the proposed control methods in Section 4. Finally, the conclusions are drawn in Section 5.

2. Mathematical Model and Problem Formulation

In this section, the mathematical model of QUAV is established first. Then, for the convenience of controller design, the model is simplified. Some useful lemmas are also provided.

2.1. Mathematical Model of QUAV

The QUAV framework with a symmetrical structure is investigated as illustrated in Figure 1. The flight motions can be subdivided into translational motion and rotational motion. To describe the kinematics and dynamics of QUAV, the earth-fixed frame E 0 : o e x e y e z e and the body-fixed frame B 0 : o b x b y b z b are defined as shown in Figure 1. For the translational motion, p = [ x , y , z ] T with respect to E 0 and v = [ u , v , w ] T with respect to B 0 stand for the position and the velocity of the QUAV, respectively. For the rotational motion, Θ = [ ϕ , θ , ψ ] T in E 0 and ω = [ p , q , r ] T in B 0 represent the Euler angle and the angular velocity of the QUAV, respectively, where the three channels represent roll, pitch, and yaw channels in sequence.
Based on the previous descriptions and the Newton–Euler equation, the mathematical model of translational and rotational motion can be expressed as
{ p ˙ = v v ˙ = g e z + T m R e z + Δ 1
{ Θ ˙ = W ω J ω ˙ = ω × J ω + τ + Δ 2
where g is the gravitational acceleration and e z = [ 0 , 0 , 1 ] T , m is the total mass, Δ 1 and Δ 2 are external disturbances, T and τ are the total thrust and attitude control torque vector, respectively. J = J 0 + Δ J denotes the inertia matrix, J 0 is the nominal part of the inertia matrix and Δ J is the uncertainty part. From [32], the inverse of the inertia matrix J can be expressed as
J 1 = ( J 0 + Δ J ) 1 = ( J 0 ) 1 ( J 0 + J 0 Δ J 1 J 0 ) 1 = J 0 1 + Δ J ˜
R is the rotation matrix from body-fixed frame to earth-fixed frame
R = [ C θ C ψ C ψ S θ S ϕ S ψ C ϕ C ψ S θ C ϕ + S ψ S ϕ C θ S ψ S ψ S θ S ϕ + C ψ C ϕ S ψ S θ C ϕ C ψ S ϕ S θ C θ S ϕ C θ C ϕ ]
where S represents sin ( ) and C denotes cos ( ) . The skew-symmetric matrix ω × and matrix W are
ω × = [ 0 r q r 0 p q p 0 ]
W = [ 1 tan θ sin ϕ tan θ cos ϕ 0 cos ϕ sin ϕ 0 sin ϕ / cos θ cos ϕ / cos θ ]

2.2. Problem Formulation

The mathematical model of translational and rotational subsystems with external disturbances is presented in Equations (1) and (2). Moreover, inertia uncertainties, actuator faults, and input saturation are also considered in this paper.
Considering uncertainties of inertia matrix, M denotes the derivative of the Euler angle, then the mathematical model of rotational subsystem can be arranged as
{ Θ ˙ = Μ M ˙ = W J 0 1 τ + f
where f = W ˙ ω + W J 1 ( ω × J ω + Δ 2 ) + W Δ J ˜ τ .
The total thrust and attitude control torque are determined by the combination of four rotors’ speed. The relationship between the total thrust T and attitude control torque τ = [ τ x τ y τ z ] T with the rotor speed can be obtained as
T = c 1 ν ϖ = [ c T c T c T c T ] [ ϖ 1 2 ϖ 2 2 ϖ 3 2 ϖ 4 2 ] T
τ = c 2 ν ϖ = [ 0 d c T 0 d c T d c T 0 d c T 0 c M c M c M c M ] [ ϖ 1 2 ϖ 2 2 ϖ 3 2 ϖ 4 2 ]
where c T and c M are constants, d denotes the distance from rotor endpoint to QUAV center of gravity, and ϖ i is the rotor speed, i = 1 , 2 , 3 , 4 .
The actuators of QUAV are four motors, and the multiplicative fault for actuators is considered here. The actuator faults can be described as
ϖ d i = ρ i ϖ c i i = 1 , 2 , 3 , 4
where ρ i [ 0 , 1 ] , ϖ c i represents the commanded rotor speed by controller and ϖ d i is the desired rotor speed. Then, the total thrust and attitude control torque can be rewritten as
T d = c 1 ν ϖ d = c 1 ρ ν ϖ c
τ d = c 2 ν ϖ d = c 2 ρ ν ϖ c
where ρ = diag ( ρ 1 2 , ρ 2 2 , ρ 3 2 , ρ 4 2 ) stands for the health condition of actuators, ν ϖ d = [ ϖ d 1 2 , ϖ d 2 2 , ϖ d 3 2 , ϖ d 4 2 ] T , ν ϖ c = [ ϖ c 1 2 , ϖ c 2 2 , ϖ c 3 2 , ϖ c 4 2 ] T , T d , and τ d are the desired total thrust and attitude control torque, respectively.
However, the actual total thrust and attitude control torque may not reach the desired value due to the upper limit of rotor speed. Therefore, input saturation is considered in this paper. Let ϖ max i denote the saturation limit of ϖ i . The ϖ i is limited in the following form
ϖ i = { ϖ d i , ϖ d i ϖ max i ϖ max i , ϖ d i > ϖ max i i = 1 , 2 , 3 , 4
where ϖ max i is a constant normally. However, when actuator faults occur, the saturation limit of rotor speed will also be influenced. Therefore, the saturation limit is different based on whether it occurs with actuator faults or without actuator faults in this paper.
Combining Equations (10) and (13), the relationship between the actual rotor speed ϖ i and the commanded rotor speed by controller ϖ c i can be obtained as
ϖ i = { ρ i ϖ c i , ρ i ϖ c i ϖ max i ϖ max i , ρ i ϖ c i > ϖ max i i = 1 , 2 , 3 , 4
For the translational subsystem, we define the virtual control vector as u = [ u x , u y , u z ] T = g e z + T m R e z . According to Equation (8), the virtual control vector can be expressed as
u = g e z + T m R e z = g e z + 1 m c 1 ν ϖ R e z = g e z + 1 m c 1 ν ϖ c R e z + 1 m c 1 ( ν ϖ ν ϖ c ) R e z = g e z + T c m R e z + 1 m c 1 ( ν ϖ ν ϖ c ) R e z = u c + 1 m c 1 ( ν ϖ ν ϖ c ) R e z
where T c = c 1 ν ϖ c denotes the commanded total thrust by controller. Then, we can obtain a simplified mathematical model
{ p ˙ = v v ˙ = u c + d
where u c = [ u c x , u c y , u c z ] T is the controller to be designed for the translational subsystem, and d = [ d x , d y , d z ] T is the lumped disturbance containing external disturbances, actuator faults, and input saturation and can be expressed as
d = 1 m c 1 ( ν ϖ ν ϖ c ) R e z + Δ 1
According to Equation (9), we can obtain
τ = c 2 ν ϖ = c 2 ν ϖ c + c 2 ( ν ϖ ν ϖ c ) = τ c + c 2 ( ν ϖ ν ϖ c )
where τ c = [ τ c x , τ c y , τ c z ] T = c 2 ν ϖ c represents the attitude control torque to be designed. Substituting Equation (18) into Equation (7), the mathematical model of the rotational subsystem can be rewritten as
{ Θ ˙ = M M ˙ = W J 0 1 τ c + D
where D = [ D x , D y , D z ] T is the lumped disturbance of the rotational subsystem containing external disturbances, inertia uncertainties, actuator faults, and input saturation and can be expressed as
D = W J 0 1 c 2 ( ν ϖ ν ϖ c ) + f
The control objective of this work is to design fixed-time UIOs to estimate the lumped disturbance d and D quickly and accurately without awareness of the boundary of disturbances or the derivatives. Then, the fixed-time tracking controllers u c and τ c are designed to stabilize tracking error to zero in fixed time regardless of initial values.

2.3. Notation and Preliminaries

Consider the nonlinear dynamical system
x ˙ = f ( x ) , x ( 0 ) = x 0
where x = [ x 1 , x 2 , , x n ] T R n is the state vector and f ( x ) is the nonlinear function.
Definition 1
[33]. The origin of system (21) is globally fixed time stable if it is globally uniformly finite time stable and the settling time function  T is globally bounded. There exists a finite constant T max R + such that T T max and x ( t ) = 0 for t T and x 0 R n . Let r = [ r 1 , , r n ] T R n , r i > 0 be the weight vector. For any λ > 0 , x R n , and the dilation mapping can be defined as Λ r ( x ) = [ λ r 1 x 1 , λ r 2 x 2 , , λ r n x n ] T R n .
Definition 2
[34]. A vector field f ( x ) R n is said to be homogeneous with degree  k with respect to weight vector  r = [ r 1 , , r n ] T R n if we have  f i ( Λ r ( x ) ) = λ k + r i f i ( x ) for any  λ > 0 , x R n .
Definition 3
[34]. A vector field f ( x ) R n is said to be homogeneous in the 0-limit with triple  ( r 0 , k 0 , f 0 ( x ) ) , where  r 0 R n is the weight vector,  k 0 is the homogeneous degree, and  f 0 ( x ) is the approximating function for  lim λ 0 max x R n \ 0 λ k 0 Λ r 0 1 f ( Λ r 0 ( x ) ) f 0 ( x ) = 0 .
Definition 4
[34]. A vector field f ( x ) R n is said to be homogeneous in the  -limit with triple  ( r , k , f ( x ) ) , where  r R n is the weight vector,  k is the homogeneous degree, and  f ( x ) is the approximating function for  lim λ max x R n \ 0 λ k Λ r 1 f ( Λ r ( x ) ) f ( x ) = 0 .
Definition 5
[35]. A vector field  f ( x ) R n is said to be homogeneous in the bi-limit if it is homogeneous in the 0-limit and  -limit simultaneously.
Lemma 1
[35]. For nonlinear system (21), we suppose the vector field  f ( x ) R n is homogeneous in the 0-limit and  -limit simultaneously with triples  ( r 0 , k 0 , f 0 ( x ) ) and  ( r , k , f ( x ) ) . If the origin of system  x ˙ = f ( x ) x ˙ = f 0 ( x ) x ˙ = f ( x ) is globally asymptotically stable, then the origin of (21) is fixed time stable when  k > 0 > k 0 holds.
Lemma 2
[33]. Suppose  Ψ = R n and there exists a Lyapunov function  V ( x ( t ) ) that satisfies  V ˙ ( x ( t ) ) ( l 1 V ϑ 1 ( x ( t ) ) + l 2 V ϑ 2 ( x ( t ) ) ) . Parameters  l 1 , l 2 > 0 , 0 < ϑ 1 < 1 and ϑ 2 > 1 are positive constants, then the origin of system (21) is fixed time stable and T ( x ( 0 ) ) 1 / [ l 1 ( 1 ϑ 1 ) ] + 1 / [ l 2 ( ϑ 2 1 ) ] for any given initial conditions x ( 0 ) R n .
Notation 1.
x = x T x denotes the Euclidean norm, where x = [ x 1 , x 2 , , x n ] T R n . 2) For any  x = [ x 1 , x 2 , , x n ] T R n s i g n ( x ) = x / x s i g p ( x ) = x p s i g n ( x ) sgn ( x ) = [ s i g n ( x 1 ) , , s i g n ( x n ) ] T x a = [ | x 1 | a s i g n ( x 1 ) , , | x n | a s i g n ( x n ) ] T .

3. Control Scheme

The control block diagram of the proposed control method is presented in Figure 2. The translational subsystem is used to achieve position control with the goal reference p d = [ x d , y d , z d ] T . The rotational subsystem is used to realize attitude control with the command Θ d = [ ϕ d , θ d , ψ d ] T , which is specifically generated via the attitude command generator. The attitude command generator can be obtained as
{ ϕ d = arcsin ( u c x sin ψ d u c y cos ψ d u c x 2 + u c y 2 + ( u c z + g ) 2 ) θ d = arctan ( u c x cos ψ d + u c y sin ψ d u c z + g )
T c = m ( u c x ( sin θ cos ψ cos ϕ + sin ψ sin ϕ ) + u c y ( sin θ sin ψ cos ϕ cos ψ sin ϕ ) + ( u c z + g ) cos θ cos ϕ )
Remark 1.
Unlike the previous articles, the desired command x d , y d , z d , ψ d comes from a real trajectory, which considers various constraints. The pseudospectral method is used to generate a discrete trajectory, and then the time-continuous trajectory, which is the reference trajectory in this paper, can be obtained by polynomial curve fitting.

3.1. Design of Translational Subsystem

In this subsection, the fixed-time UIO will be designed for translational subsystem to estimate the lumped disturbance d with a unique feature that it does not require the upper bound of disturbance or its derivative. Then, a fixed-time tracking controller u c will be presented based on the simplified mathematical model in Equation (16) to stabilize tracking error into a small region in fixed time.
We design the fixed-time UIO first. To begin with, a new state η = [ η 1 , η 2 , η 3 ] T is introduced as follows
η ( t ) = L 1 ( u c + L 2 η ) d t + L 1 v
where both L 1 = diag [ l 11 , l 12 , l 13 ] , L 2 = diag [ l 21 , l 22 , l 23 ] are positive definitions of the diagonal matrix. Differentiating η with respect to time, combining Equation (16), we can obtain the relationship between the new state η and the lumped disturbance d
η ˙ = L 1 L 2 η + L 1 d
Theorem 1.
Consider the translational subsystem with external disturbances, actuator faults, and input saturation described by Equation (16), define d ^ = [ d ^ x , d ^ y , d ^ z ] T as the estimation of the lumped disturbance d , if the observer is designed as
{ η ˙ = L 1 ( u c + L 2 η ) + L 1 v ˙ η ^ ˙ = η ˙ + H 1 s i g m 1 ( η η ^ ) + H 2 s i g n 1 ( η η ^ ) d ^ = L 2 η ^ + L 1 1 η ˙
where η ^ = [ η ^ 1 , η ^ 2 , η ^ 3 ] T denotes the estimation of η , both H 1 = diag [ h 11 , h 12 , h 13 ] , H 2 = diag [ h 21 , h 22 , h 23 ] represent the positively defined matrix and 0 < m 1 < 1 , n 1 > 1 are positive constants. Then, the lumped disturbance d can be approximated in fixed time through d ^ .
Proof of Theorem 1.
Define η ˜ = η η ^ = [ η ˜ 1 , η ˜ 2 , η ˜ 3 ] T , choose the Lyapunov function as V 0 = 1 2 η ˜ T η ˜ . Differentiating V 0 with respect to time, we can obtain
V ˙ 0 = η ˜ T η ˜ ˙     = η ˜ T ( η ˙ ( η ˙ + H 1 s i g m 1 ( η ˜ ) + H 2 s i g n 1 ( η ˜ ) ) )     = η ˜ m 1 1 η ˜ T H 1 η ˜ T η ˜ n 1 1 η ˜ T H 2 η ˜ T     = η ˜ m 1 1 i = 1 3 h 1 i η ˜ i 2 η ˜ n 1 1 i = 1 3 h 2 i η ˜ i 2     h 1 i , min η ˜ m 1 + 1 h 2 i , min η ˜ n 1 + 1     2 ( m 1 + 1 ) / 2 h 1 i , min V 0 ( m 1 + 1 ) / 2 2 ( n 1 + 1 ) / 2 h 2 i , min V 0 ( n 1 + 1 ) / 2
where h 1 i , min = min { h 1 i } , h 2 i , min = min { h 2 i } , i = 1 , 2 , 3 . From Lemma 2, it can be concluded that η ˜ will converge to zero in fixed time, and the convergence time can be estimated by
T 0 = 2 ( 1 m 1 ) / 2 h 1 i , min ( 1 m 1 ) + 2 ( 1 n 1 ) / 2 h 2 i , min ( n 1 1 )
Define d ˜ = d d ^ = [ d ˜ x , d ˜ y , d ˜ z ] T as the estimation error between d and d ^ . Combining Equation (25), we can obtain
d ˜ = d d ^ = L 2 η + L 1 1 η ˙ ( L 2 η ^ + L 1 1 η ˙ )     = L 2 η ˜
According to Equation (28), we can conclude that the observation error d ˜ of fixed-time UIO will converge to origin within fixed time.
The proof is completed. □
Remark 2.
Compared with ESO in [21,22] and STDO in [25,26], the proposed fixed-time UIO in Equation (26) does not include the sign function and the prior knowledge of the disturbance, so the observer is simple to design and there is no chattering when estimating the disturbance.
Then, the fixed-time tracking controller is proposed based on the accurate estimation of the lumped disturbance d . Firstly, the tracking error of translational subsystem can be defined as
{ e p = p p d e v = v v d
where p d = [ x d , y d , z d ] T and v d = [ u d , v d , w d ] T are desired commands of position and speed, respectively, and satisfy p ˙ d = v d .
Based on the definition of tracking error in Equation (30), the fixed-time tracking controller for the translational subsystem is designed as
u c = v ˙ d k p ( s i g α 1 ( e p ) + s i g β 1 ( e p ) ) k v ( s i g α 2 ( e v ) + s i g β 2 ( e v ) ) d ^
where α 1 = α / ( 2 α ) , α 2 = α , α ( 0 , 1 ) , β 1 = β 2 / ( 2 β 2 ) , β 2 = ( 4 3 α ) / ( 3 2 α ) , k p , k v are positive constants.
Theorem 2.
Consider the translational subsystem described in Equation (16). If the lumped disturbancedcan be estimated accurately by fixed-time UIO presented in Equation (26), then the tracking error could converge to zero in fixed-time with the proposed controller in Equation (31).
Proof of Theorem 2.
According to the fixed-time UIO, the estimation error d ˜ will converge to zero for t T 0 . Then, we can obtain the error equation from Equation (16) and Equation (30)
{ e ˙ p = e v e ˙ v = k p ( s i g α 1 ( e p ) + s i g β 1 ( e p ) ) k v ( s i g α 2 ( e v ) + s i g β 2 ( e v ) )
For the convenience of bi-limit homogeneity analysis, based on Equation (32), a vector field is defined as
f ( e p , e v ) = [ e v k p ( s i g α 1 ( e p ) + s i g β 1 ( e p ) ) k v ( s i g α 2 ( e v ) + s i g β 2 ( e v ) ) ]
Two auxiliary vector fields are designed as
f 0 = [ e v k p s i g α 1 ( e p ) k v s i g α 2 ( e v ) ]
f = [ e v k p s i g β 1 ( e p ) k v s i g β 2 ( e v ) ]
To achieve 0-limit homogeneity with triple ( r 0 , k 0 , f 0 ) of Equation (33), we set k 0 = 1 and define the dilation mapping as Λ r 0 ( e p , e v ) = [ λ ( 2 α ) / ( 1 α ) e p , λ 1 / ( 1 α ) e v ] T .
The detailed 0-limit homogeneity approximation is presented as follows
λ k 0 Λ r 0 1 f ( Λ r 0 ( e p , e v ) ) = λ 1 [ λ 2 α 1 α λ 1 1 α e v λ 1 1 α ( λ 2 α 1 α α 1 k p s i g α 1 ( e p ) + λ 2 α 1 α β 1 k p s i g β 1 ( e p ) + λ 1 1 α α 2 k v s i g α 2 ( e v ) + λ 1 1 α β 2 k v s i g β 2 ( e v ) ) ] = [ e v k p s i g α 1 ( e p ) k v s i g α 2 ( e v ) λ α 1 α ( λ 2 α 1 α β 1 k p s i g β 1 ( e p ) + λ 1 1 α β 2 k v s i g β 2 ( e v ) ) ] = [ e v k p s i g α 1 ( e p ) k v s i g α 2 ( e v ) ( λ 4 k p s i g β 1 ( e p ) + λ 4 2 α 3 2 α k v s i g β 2 ( e v ) ) ]
Then, we can obtain
lim λ 0 max ( e p , e v ) R n \ 0 λ k 0 Λ r 0 1 f ( Λ r 0 ( e p , e v ) ) f 0 = lim λ 0 max ( e p , e v ) R n \ 0 0 λ 4 k p s i g β 1 ( e p ) + λ 4 2 α 3 2 α k v s i g β 2 ( e v ) = 0
According to Definition 3, the Equation (33) is homogeneous in 0-limit.
In a similar way, the homogeneity in -limit with triple ( r , k , f ) of Equation (33) is also analyzed. We set homogeneous degree k = 1 , and define the dilation mapping as Λ r ( e p , e v ) = [ λ ( 2 β 2 ) / ( β 2 1 ) e p , λ 1 / ( β 2 1 ) e v ] T . Then, we have
λ k Λ r 1 f ( Λ r ( e p , e v ) ) = λ 1 [ λ 2 β 2 β 2 1 λ 1 β 2 1 e v λ 1 β 2 1 ( λ 2 β 2 β 2 1 α 1 k p s i g α 1 ( e p ) + λ 2 β 2 β 2 1 β 1 k p s i g β 1 ( e p ) + λ 1 β 2 1 α 2 k v s i g α 2 ( e v ) + λ 1 β 2 1 β 2 k v s i g β 2 ( e v ) ) ] = [ e v k p s i g β 1 ( e p ) k v s i g β 2 ( e v ) λ β 2 β 2 1 ( λ 2 β 2 β 2 1 α 1 k p s i g α 1 ( e p ) + λ 1 β 2 1 α 2 k v s i g α 2 ( e v ) ) ] = [ e v k p s i g β 1 ( e p ) k v s i g β 2 ( e v ) ( λ 4 k p s i g α 1 ( e p ) + λ 2 ( 2 α ) k v s i g α 2 ( e v ) ) ]
Then, we can obtain
lim λ max ( e p , e v ) R n \ 0 λ k Λ r 1 f ( Λ r ( e p , e v ) ) f = lim λ max ( e p , e v ) R n \ 0 0 λ 4 k p s i g α 1 ( e p ) + λ 2 ( 2 α ) k v s i g α 2 ( e v ) = 0
According to Definition 4, the Equation (33) is homogeneous in -limit. Therefore, the vector field described in Equation (33) is homogeneous in bi-limit.
Then, we prove that Equations (33)–(35) are globally asymptotically stable.
For Equation (33), we define the Lyapunov function as
V 1 = 2 k p ( α 1 + 1 ) e p β 1 + 1 + 2 k p ( β 1 + 1 ) e p α 1 + 1 + ( α 1 + 1 ) ( β 1 + 1 ) e v T e v
Differentiating V 1 with respect to time, then we can obtain
V ˙ 1 = 2 k p ( α 1 + 1 ) ( β 1 + 1 ) e p β 1 1 e p T e ˙ p + 2 k p ( α 1 + 1 ) ( β 1 + 1 ) e p α 1 1 e p T e ˙ p   + 2 ( α 1 + 1 ) ( β 1 + 1 ) e v T e ˙ v     = 2 k p ( α 1 + 1 ) ( β 1 + 1 ) ( e p α 1 1 + e p β 1 1 ) e p T e ˙ p   + 2 ( α 1 + 1 ) ( β 1 + 1 ) e v T ( k p ( s i g α 1 ( e p ) + s i g β 1 ( e p ) ) k v ( s i g α 2 ( e v ) + s i g β 2 ( e v ) ) )     = 2 k p ( α 1 + 1 ) ( β 1 + 1 ) ( e p α 1 1 + e p β 1 1 ) e p T e v   2 k p ( α 1 + 1 ) ( β 1 + 1 ) ( e p α 1 1 + e p β 1 1 ) e v T e p   2 k v ( α 1 + 1 ) ( β 1 + 1 ) ( e v α 2 + 1 + e v β 2 + 1 )     = 2 k v ( α 1 + 1 ) ( β 1 + 1 ) ( e v α 2 + 1 + e v β 2 + 1 ) 0
Let V · = 0 ; thus, we can obtain e v = 0 . Then, substituting it into Equation (33), we have e ˙ p = 0 , e ˙ v = k p ( s i g α 1 ( e p ) + s i g β 1 ( e p ) ) , which is the set of all points where V ˙ 1 = 0 . Supposing { e p 0 , e v = 0 } , then there always exists e ˙ v 0 , namely e v not always equal to 0 . Therefore, { e p 0 , e v = 0 } is not an invariant set. When { e p = 0 , e v = 0 } , there would always exist e ˙ p = 0 , e ˙ v = 0 , namely e p = 0 , e v = 0 would always hold. Consequently, the only invariant set is { e p = 0 , e v = 0 } . According to LaSalle invariant theorem, { e p = 0 , e v = 0 } is the global equilibrium point and global asymptotical stability can be guaranteed.
For Equation (34), the Lyapunov function can be designed as
V 2 = 2 k p e p α 1 + 1 + ( α 1 + 1 ) e v T e v
The time derivative of V 2 is
V ˙ 2 = 2 k p ( α 1 + 1 ) e p α 1 1 e p T e v + 2 ( α 1 + 1 ) e v T ( k p s i g α 1 ( e p ) k v s i g α 2 ( e v ) )   = 2 k v ( α 1 + 1 ) e v α 2 + 1   0
According to the LaSalle invariant theorem, Equation (34) is globally asymptotically stable.
For Equation (35), the global asymptotical stability is also proved with the following Lyapunov function
V 3 = 2 k p e p β 1 + 1 + ( β 1 + 1 ) e v T e v
According to Lemma 1, the origin of Equation (32) is fixed-time stable.
The proof is completed. □

3.2. Design of Rotational Subsystem

In this subsection, the fixed-time UIO for rotational subsystem is presented to estimate accurately the lumped disturbance D . Then, a fixed-time tracking controller τ c is proposed based on the simplified mathematical model of rotational subsystem in Equation (19).
Like the translational subsystem, firstly, a new state ς = [ ς 1 , ς 2 , ς 3 ] T is designed as
ς ( t ) = L 3 ( W J 0 1 τ c + L 4 ς ) d t + L 3 M
where both L 3 = diag [ l 31 , l 32 , l 33 ] , L 4 = diag [ l 41 , l 42 , l 43 ] are positive define diagonal matrices. Differentiating ς ( t ) with respect to time, combining Equation (19), then we can obtain
ς ˙ = L 3 L 4 ς + L 3 D
Theorem 3.
Consider the rotational subsystem with external disturbances, inertia uncertainties, actuator faults, and input saturation described by Equation (19), then define D ^ = [ D ^ x , D ^ y , D ^ z ] T , ς ^ = [ ς ^ 1 , ς ^ 2 , ς ^ 3 ] T as the estimation of D and ς respectively, if the observer is designed as
{ ς ˙ = L 3 ( W J 0 1 τ c + L 4 ς ) + L 3 M ˙ ς ^ ˙ = ς ˙ + H 3 s i g m 2 ( ς ς ^ ) + H 4 s i g n 2 ( ς ς ^ ) D ^ = L 4 ς ^ + L 3 1 ς ˙
where both H 3 = diag [ h 31 , h 32 , h 33 ] , H 4 = diag [ h 41 , h 42 , h 43 ] are positive define matrices and 0 < m 2 < 1 , n 2 > 1 are positive constants. Then, the lumped disturbance D can be approximated in fixed time through D ^ .
Proof of Theorem 3.
Define ς ˜ = ς ς ^ = [ ς ˜ 1 , ς ˜ 2 , ς ˜ 3 ] T , then construct the Lyapunov function as V 4 = 1 2 ς ˜ T ς ˜ . Differentiating V 4 with respect to time, we can obtain
V ˙ 4 = ς ˜ T ς ˜ ˙         = ς ˜ T ( ς ˙ ( ς ˙ + H 3 s i g m 2 ( ς ˜ ) + H 4 s i g n 2 ( ς ˜ ) ) )         = ς ˜ m 2 1 ς ˜ T H 3 ς ˜ ς ˜ n 2 1 ς ˜ T H 4 ς ˜         = ς ˜ m 2 1 i = 1 3 h 3 i ς ˜ i 2 ς ˜ n 2 1 i = 1 3 h 4 i ς ˜ i 2         h 3 i , min ς ˜ m 2 + 1 h 4 i , min ς ˜ n 2 + 1         2 ( m 2 + 1 ) / 2 h 3 i , min V 4 ( m 2 + 1 ) / 2 2 ( n 2 + 1 ) / 2 h 4 i , min V 4 ( n 2 + 1 ) / 2
where h 3 i , min = min { h 3 i } , h 4 i , min = min { h 4 i } , i = 1 , 2 , 3 . From Lemma 2, we can conclude that ς ˜ will converge to zero in fixed time, and the convergence time can be estimated by
T 1 = 2 ( 1 m 2 ) / 2 h 3 i , min ( 1 m 2 ) + 2 ( 1 n 2 ) / 2 h 4 i , min ( n 2 1 )
Define D ˜ = D D ^ = [ D ˜ x , D ˜ y , D ˜ z ] T as the estimation error between D and D ^ . Combining Equation (46), we can obtain
D ˜ = D D ^ = L 4 ς + L 3 1 ς ˙ ( L 4 ς ^ + L 3 1 ς ˙ )     = L 4 ς ˜
According to Equation (49), we can conclude that the observation error D ˜ will converge to origin within fixed time.
The proof is completed. □
Then, the controller of rotational subsystem is designed.
The tracking error can be defined as
{ e 1 = Θ Θ d e 2 = M M d
where Θ d = [ ϕ d , θ d , ψ d ] T and M d = [ ϕ ˙ d , θ ˙ d , ψ ˙ d ] T stand for the desired Euler angle and its derivative, respectively, and satisfy Θ ˙ d = M d .
Based on the precise estimation of the lumped disturbance D , the fixed-time tracking controller τ c is designed as
{ τ c = J 0 W 1 ( M ˙ d μ D ^ ) μ = k 1 ( s i g a 1 ( e 1 ) + s i g b 1 ( e 1 ) ) + k 2 ( s i g a 2 ( e 2 ) + s i g b 2 ( e 2 ) )
where a 1 = a / ( 2 a ) , a 2 = a , a ( 0 , 1 ) , b 1 = b 2 / ( 2 b 2 ) , b 2 = ( 4 3 a ) / ( 3 2 a ) , k 1 , k 2 are positive constants.
Theorem 4.
Consider the rotational subsystem described by Equation (19). If the lumped disturbanceDcan be estimated precisely by the fixed-time UIO presented in Equation (47), then the tracking error could converge to zero in fixed-time with the designed controller in Equation (52).
Proof of Theorem 4.
According to Equation (47), the estimation error D ˜ will converge to zero for t T 1 . Then, we can obtain
{ e ˙ 1 = e 2 e ˙ 2 = k 1 ( s i g a 1 ( e 1 ) + s i g b 1 ( e 1 ) ) k 2 ( s i g a 2 ( e 2 ) + s i g b 2 ( e 2 ) )
According to the proof process of Theorem 2, the origin of Equation (53) is fixed-time stable.
The proof is completed. □

4. Simulation and Discussion

The simulation results of the proposed control scheme are presented in this section. Then, comparative simulations are carried out to demonstrate the superiority of the fixed-time UIO proposed in this paper.

4.1. Simulation Results

The parameters of the QUAV are given as m = 0.5   kg , g = 9.8   m / s 2 , c T = 2.98 × 10 6 , c M = 1.14 × 10 7 , d = 0.225   m , J = 10 3 × diag [ 4.9 4.9 8.8 ]   kg m 2 , J 0 = 10 3 × diag [ 4.5 4.5 8 ]   kg m 2 . The health condition of actuators is represented as
{ ρ 1 = ρ 2 = ρ 3 = ρ 4 = 1 0 t < 10 ρ 1 = ρ 2 = ρ 3 = ρ 4 = 0.9 , 10 t 15 ρ 1 = ρ 2 = ρ 3 = ρ 4 = 1 , t > 15
In fact, the saturation limit of rotor speed is influenced in the case of actuator faults, so it can be set as
ϖ max i = { δ ϖ ¯ max i , with actuator faults ϖ ¯ max i , without actuator faults
where δ is determined by the actual situation and is set to δ = 0.9 , and ϖ ¯ max i represents the saturation limit without actuator faults, ϖ ¯ max 1 = ϖ ¯ max 2 = ϖ ¯ max 3 = ϖ ¯ max 4 = 900   rpm .
The external disturbances of the translational and rotational subsystem are supposed to be
Δ 1 = 0.1 × [ 1 + sin ( 0.3 t ) + cos ( 0.8 t ) 1 + sin ( 0.6 t ) cos ( 0.4 t ) 1 + sin ( 0.5 t ) + cos ( 0.2 t ) ] N / kg
Δ 2 = 0.05 × [ 1 + sin ( 0.3 t ) + cos ( 0.8 t ) 1 + sin ( 0.6 t ) cos ( 0.4 t ) 1 + sin ( 0.5 t + π / 2 ) + cos ( 0.2 t + π / 2 ) ] Nm
The parameters of fixed-time tracking controller and fixed-time UIO for the translational subsystem are set as k p = 0.3 , k v = 0.9 , α = 0.8 , m 1 = 0.8 , n 1 = 1.2 , L 1 = diag [ 50 50 50 ] , L 2 = H 1 = H 2 = diag [ 1 1 1 ] .
The parameters of the fixed-time tracking controller and fixed-time UIO for rotational subsystem are set as k 1 = 80 , k 2 = 10 , a = 0.8 , m 2 = 0.8 , n 2 = 1.2 , L 3 = diag [ 40 40 40 ] , L 4 = H 3 = H 4 = diag [ 1 1 1 ] .
The tracking results of position and speed are presented in Figure 3, where V = u 2 + v 2 + w 2 . Figure 4 shows the trajectory tracking results from 3D view and X-Y view. Obviously, the proposed fixed-time tracking controller for the translational subsystem can track the desired command accurately. Figure 5 presents the simulation results of the Euler angle and its tracking error, and we can conclude that the tracking performance of the designed fixed-time tracking controller for the rotational subsystem is satisfactory, and the tracking error could rapidly converge to zero. The rotor speed commanded by the controller and the actual rotor speed are presented in Figure 6, and we can see that the commanded rotor speed increases in order to track the desired command in the event of actuator faults. The saturation limit of rotor speed decreases during actuator faults. Furthermore, even though the controller orders a high rotor speed, the actual rotor speed does not surpass the saturation limit we set. Figure 7 depicts the controller’s commanded thrust and torque as well as the actual thrust and torque. The thrust and torque commanded, as well as the actual thrust and torque, are different due to the difference between the commanded rotor speed and the actual rotor speed at the start and in the case of actuator faults. The estimation results of lumped disturbances for both translational and rotational subsystem are drawn in Figure 8, and the observation performance of the proposed fixed-time UIO is acceptable. Though the lumped disturbances are very large at first, which is caused by the fact that the actual rotor speed does not match the commanded rotor speed in the starting stage, the proposed fixed-time UIO can estimate the lumped disturbances quickly and precisely.

4.2. Comparison

To verify the performance of the proposed fixed-time UIO (FUIO), simulations using the fixed-time extend state observer (FESO) in [22] and the fixed-time integral sliding mode disturbance observer (FISMDO) in [25] are carried out for comparison.
The FESO for rotational subsystem is designed as
{ Θ ^ ˙ = M ^ + p 1 s i g m 3 ( ξ 1 ) + q 1 s i g n 3 ( ξ 1 ) M ^ ˙ = W J 0 1 τ c + D ^ + p 2 s i g m 4 ( ξ 1 ) + q 2 s i g n 4 ( ξ 1 ) D ^ ˙ = p 3 s i g m 5 ( ξ 1 ) + q 3 s i g n 5 ( ξ 1 ) + γ 1 s i g n ( ξ 1 )
where ξ 1 = Θ Θ ^ , Θ ^ , M ^ , D ^ are the estimations of Θ , M , D , respectively. The parameters are set as m 3 = 0.7 , m 4 = 2 m 3 1 , m 5 = 3 m 3 2 , n 3 = 1.2 , n 4 = 2 n 3 1 , n 5 = 3 n 3 2 , p 1 = q 1 = 10 , p 2 = q 2 = 27 , p 3 = q 3 = 74 , γ 1 = 3 > D ˙ .
The FESO for translational subsystem is designed as
{ p ^ ˙ = v ^ + p 4 s i g m 6 ( ξ 2 ) + q 4 s i g n 6 ( ξ 2 ) v ^ ˙ = u c + d ^ + p 5 s i g m 7 ( ξ 2 ) + q 5 s i g n 7 ( ξ 2 ) d ^ ˙ = p 6 s i g m 8 ( ξ 2 ) + q 6 s i g n 8 ( ξ 2 ) + γ 2 s i g n ( ξ 2 )
where ξ 2 = p p ^ , p ^ , v ^ , d ^ are the estimations of p , v , d , respectively. The parameters are given as m 6 = 0.8 , m 7 = 2 m 6 1 , m 8 = 3 m 6 2 , n 6 = 1.2 , n 7 = 2 n 6 1 , n 8 = 3 n 6 2 , p 4 = q 4 = 10 , p 5 = q 5 = 27 , p 6 = q 6 = 15 , γ 2 = 1 > d ˙ .
The FISMDO for the rotational subsystem is designed as
{ s ˙ 0 = D D ^ s 1 = s ˙ 0 + [ λ 1 s 0 a 3 + λ 2 s 0 b 3 + λ 3 s ˙ 0 a 4 + λ 4 s ˙ 0 b 4 ] d s D ^ ˙ = λ 1 s 0 a 3 + λ 2 s 0 b 3 + λ 3 s ˙ 0 a 4 + λ 4 s ˙ 0 b 4 + λ 5 s 1 a 5 + λ 6 s 1 b 5 + γ 3 sgn ( s 1 )
where λ 1 = λ 3 = 2 , λ 2 = λ 4 = 4 , λ 5 = λ 6 = 0.5 , a 3 = a 4 = a 5 = 0.8 , b 3 = b 4 = b 5 = 1.2 , γ 3 = 3 > D ˙ .
The FISMDO for the translational subsystem is designed as
{ s ˙ 2 = d d ^ s 3 = s ˙ 2 + [ κ 1 s 2 a 6 + κ 2 s 2 b 6 + κ 3 s ˙ 2 a 7 + κ 4 s ˙ 2 b 7 ] d s d ^ ˙ = κ 1 s 2 a 6 + κ 2 s 2 b 6 + κ 3 s ˙ 2 a 7 + κ 4 s ˙ 2 b 7 + κ 5 s 3 a 8 + κ 6 s 3 b 8 + γ 4 sgn ( s 3 )
where κ 1 = κ 3 = 2 , κ 2 = κ 4 = 4 , κ 5 = κ 6 = 0.5 , a 6 = a 7 = a 8 = 0.8 , b 6 = b 7 = b 8 = 1.2 , γ 4 = 1 > d ˙ .
Figure 9 depicts the estimation error of lumped disturbances, employing three observers for translational and rotational subsystems. Obviously, the estimation error of the three observers could converge to zero. However, the estimation performance of the three observers is also different. Compared with FESO, the estimation error converges faster when using the FUIO and the FISMDO to estimate the lumped disturbances. Moreover, according to Figure 9c, we can find that the FUIO converges fastest when the lumped disturbance changes abruptly owing to actuator faults. Additionally, since the FUIO does not include the sign function, its estimation error converges more smoothly and accurately than the FESO and FISMDO. The tracking error of the position and the Euler angle when using different disturbance observers are shown in Figure 10, and we find that the convergence speed of tracking error is faster in the case of actuator faults and the tracking precision is higher when using the FUIO because it has better estimation performance.
Remark 3.
When the rotor speed saturation is considered, the lumped disturbance is very large. For ease of analysis, the rotor speed saturation is not considered in the comparison simulations.
Remark 4.
As the estimation precision of the disturbance observer can affect the lumped disturbance, the lumped disturbance is different when using different observers to estimate the lumped disturbance. Hence, we compare the estimation performance of the three observers using the estimation error.
Remark 5.
From Equations (58)–(61), we know that the FESO and FISMDO are designed based on the boundary of the derivative of the lumped disturbance. Therefore, the estimation performance of the FESO and FISMDO cannot be guaranteed in theory when the lumped disturbance changes sharply.
Remark 6.
According to the estimation error of different observers presented in Figure 9, we can conclude the advantages of the FUIO. First, the sign function is absent in FUIO, so there is no chattering and the estimation precision is higher. Additionally, the design process of FUIO does not need prior knowledge of the boundary of the disturbance or its derivative. Therefore, the design process becomes easier and the application is expanded. Moreover, the FUIO can observe the lumped disturbance accurately and quickly when the lumped disturbance changes suddenly due to actuator faults.

4.3. Discussion

4.3.1. Differences between Rotor Speed Saturation and Thrust or Torque Saturation

In contrast with [29,30,31], the rotor speed saturation is considered in this paper instead of thrust and torque saturation. According to Equation (8), we find that the total thrust is related to the rotor speed of all four motors. Therefore, as soon as the rotor speed of any of the motors reaches saturation, the total thrust will be affected. When only thrust saturation is considered, the role of each rotor speed cannot be reflected. Hence, rotor speed saturation is considered in this paper, which is more reasonable.

4.3.2. The Feature of Proposed Observer

In the process of observer design, a new state η is introduced in Equation (24), and its time derivative in Equation (25) shows that the new state η is the first-order filtering of lumped disturbance d . Using the new state, the observer is designed in Equation (26). Instead of estimating the lumped disturbance directly, the observer first observes the new state and then uses the filtering relationship between the disturbance and the new state to estimate the disturbance. The advantage of this is that the upper bound of the derivative of lumped disturbance is not required and the lumped disturbance can be estimated accurately when it changes suddenly.

4.3.3. The Future Work

In this section, simulation results that consider the external disturbance, inertia uncertainties, actuator faults, and input saturation are presented to validate the effectiveness and superiority of the proposed control scheme. Although the simulation results are satisfactory, the flight experiment verification is also indispensable, which we will explore in future work.

5. Conclusions

This paper investigates the trajectory tracking control problem for a QUAV in the presence of external disturbances, inertia uncertainties, actuator faults, and input saturation. The simulation results show that the high-accuracy estimation of the lumped disturbance can be obtained in fixed time, employing the proposed fixed-time UIO. Compared with other observers, such as the observers in [20,21,22,23,24,25,26], the advantages of proposed fixed-time UIO can be summarized as follows: (a) there is no chattering in the estimation of lumped disturbance and the estimation precision is better; (b) prior knowledge of the boundary of the disturbance or its derivative is not required; (c) the lumped disturbance can be observed accurately and quickly when it changes suddenly due to actuator faults. The fixed-time tracking controllers for translational and rotational subsystems could drive the tracking error into a small region in fixed time regardless of the initial values with the help of the proposed fixed-time UIO. Rigorous convergence analysis and numerical simulation results have demonstrated the effectiveness and superiority of the proposed control scheme.

Author Contributions

Conceptualization, S.S. and Y.Z.; methodology, S.S. and Y.Z.; software, S.X.; validation, S.X., S.S. and Y.Z.; formal analysis, X.W.; investigation, S.X. and S.S.; resources, S.X., S.S. and Y.Z.; data curation, S.X.; writing—original draft preparation, S.X.; writing—review and editing, S.S. and Y.Z.; visualization, S.X.; supervision, Y.Z.; project administration, Y.Z.; funding acquisition, S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant Number 61903122), the Natural Science Foundation of Hebei Province (Grant Number F2021208015), the Science and Technology Project of Hebei Education Department (Grant Number BJ2021003), and the National Natural Science Foundation of China (Grant Number U20A20198 and 62003129).

Data Availability Statement

No applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Serrano, M.E.; Gandolfo, D.C.; Scaglia, G.J.E. Trajectory tracking controller for unmanned helicopter under environmental disturbances. ISA Trans. 2020, 106, 171–180. [Google Scholar] [CrossRef] [PubMed]
  2. Cao, C.Y.; Wei, C.S.; Liao, Y.X.; Zhang, Y.C.; Li, J. On novel trajectory tracking control of quadrotor UAV: A finite-time guaranteed performance approach. J. Franklin Inst. 2022, 359, 8454–8483. [Google Scholar] [CrossRef]
  3. Li, B.; Song, C.; Bai, S.X.; Huang, J.Y.; Ma, R.; Wan, K.F.; Neretin, E. Multi-UAV trajectory planning during cooperative tracking based on a Fusion Algorithm integrating MPC and standoff. Drones 2023, 7, 196. [Google Scholar] [CrossRef]
  4. Kidambi, K.B.; Fermuller, C.; Aloimonos, Y.; Xu, H. Robust nonlinear control-based trajectory tracking for quadrotors under uncertainty. IEEE Control Syst. Lett. 2021, 5, 2042–2047. [Google Scholar] [CrossRef]
  5. Blas, L.A.; Davila, J.; Salazar, S.; Bonilla, M. Robust trajectory tracking for an uncertain UAV based on active disturbance rejection. IEEE Control Syst. Lett. 2022, 6, 1466–1471. [Google Scholar] [CrossRef]
  6. Singhal, K.; Kumar, V. Robust trajectory tracking control of non-holonomic wheeled mobile robots using an adaptive fractional order parallel fuzzy PID controller. J. Franklin Inst. 2022, 359, 4160–4215. [Google Scholar] [CrossRef]
  7. Mathiyalagan, K.; Sangeetha, G. Finite-time stabilization of nonlinear time delay systems using LQR based sliding mode control. J. Franklin Inst. 2019, 356, 3948–3964. [Google Scholar] [CrossRef]
  8. Lee, D.; Kim, H.J.; Sastry, S. Feedback linearization vs. adaptive sliding mode control for a quadrotor helicopter. Int. J. Control Autom. Syst. 2009, 7, 419–428. [Google Scholar] [CrossRef]
  9. Koksal, N.; An, H.; Fidan, B. Backstepping-based adaptive control of a quadrotor UAV with guaranteed tracking performance. ISA Trans. 2020, 105, 98–110. [Google Scholar] [CrossRef] [PubMed]
  10. Liu, W.Q.; Cheng, X.H.; Zhang, J.J. Command filter-based adaptive fuzzy integral backstepping control for quadrotor UAV with input saturation. J. Franklin Inst. 2023, 360, 484–507. [Google Scholar] [CrossRef]
  11. Labbadi, M.; Boukal, Y.; Cherkaoui, M.; Djemai, M. Fractional-order global sliding mode controller for an uncertain quadrotor UAVs subjected to external disturbances. J. Franklin Inst. 2021, 358, 4822–4847. [Google Scholar] [CrossRef]
  12. Labbadi, M.; Cherkaoui, M. Adaptive fractional-order nonsingular fast terminal sliding mode based robust tracking control of quadrotor UAV with Gaussian random disturbances and uncertainties. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 2265–2277. [Google Scholar] [CrossRef]
  13. Wang, Q.; Wang, W.; Suzuki, S.; Namiki, A.; Liu, H.X.; Li, Z.R. Design and implementation of UAV velocity controller based on reference model sliding mode control. Drones 2023, 7, 130. [Google Scholar] [CrossRef]
  14. Huang, D.Q.; Huang, T.P.; Qin, N.; Li, Y.A.; Yang, Y. Finite-time control for a UAV system based on finite-time disturbance observer. Aerosp. Sci. Technol. 2022, 129, 107825. [Google Scholar] [CrossRef]
  15. Cheng, W.L.; Jiang, B.; Zhang, K.; Ding, S.X. Robust finite-time cooperative formation control of UGV-UAV with model uncertainties and actuator faults. J. Franklin Inst. 2021, 358, 8811–8837. [Google Scholar] [CrossRef]
  16. Wang, J.H.; Alattas, K.A.; Bouteraa, Y.; Mofid, O.; Mobayen, S. Adaptive finite-time backstepping control tracker for quadrotor UAV with model uncertainty and external disturbance. Aerosp. Sci. Technol. 2023, 133, 108088. [Google Scholar] [CrossRef]
  17. Liu, B.J.; Li, A.J.; Guo, Y.; Wang, C.Q. Adaptive distributed finite-time formation control for multi-UAVs under input saturation without collisions. Aerosp. Sci. Technol. 2022, 120, 107252. [Google Scholar] [CrossRef]
  18. Xia, K.W.; Son, H.S. Adaptive fixed-time control of autonomous VTOL UAVs for ship landing operations. J. Franklin Inst. 2020, 357, 6175–6196. [Google Scholar] [CrossRef]
  19. Tan, J.; Dong, Y.F.; Shao, P.Y.; Qu, G.M. Anti-saturation adaptive fault-tolerant control with fixed-time prescribed performance for UAV under AOA asymmetric constraint. Aerosp. Sci. Technol. 2022, 120, 107264. [Google Scholar] [CrossRef]
  20. Cui, L.; Hou, X.Y.; Zuo, Z.Q.; Yang, H.J. An adaptive fast super-twisting disturbance observer-based dual closed-loop attitude control with fixed-time convergence for UAV. J. Franklin Inst. 2022, 359, 2514–2540. [Google Scholar] [CrossRef]
  21. Chen, L.L.; Liu, Z.B.; Dang, Q.Q.; Zhao, W.; Wang, G.D. Robust trajectory tracking control for a quadrotor using recursive sliding mode control and nonlinear extended state observer. Aerosp. Sci. Technol. 2022, 128, 107749. [Google Scholar] [CrossRef]
  22. Shao, S.S.; Wang, S.; Zhao, Y.J. Fixed time output feedback control for quadrotor unmanned aerial vehicle under disturbances. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2022, 236, 3554–3566. [Google Scholar] [CrossRef]
  23. Cui, L.; Zhang, R.Z.; Yang, H.J.; Zuo, Z.Q. Adaptive super-twisting trajectory tracking control for an unmanned aerial vehicle under gust winds. Aerosp. Sci. Technol. 2021, 115, 106833. [Google Scholar] [CrossRef]
  24. Gonzalez, J.A.C.; Pena, O.S.; Morales, J.D.L. Observer-based super twisting design: A comparative study on quadrotor altitude control. ISA Trans. 2021, 109, 307–314. [Google Scholar] [CrossRef] [PubMed]
  25. Su, B.; Wang, H.B.; Wang, Y.L. Dynamic event-triggered formation control for AUVs with fixed-time integral sliding mode disturbance observer. Ocean Eng. 2021, 240, 109893. [Google Scholar] [CrossRef]
  26. Dong, Z.; Liu, L.; Wang, S.P. Sliding mode disturbance observer-based adaptive dynamic inversion fault-tolerant control for Fixed-Wing UAV. Drones 2022, 6, 295. [Google Scholar] [CrossRef]
  27. Xuan-Mung, N.; Golestani, M. Energy-efficient disturbance observer-based attitude tracking control with fixed-time convergence for spacecraft. IEEE. Trans. Aero. Elec. Syst. 2022, 1–10. [Google Scholar] [CrossRef]
  28. Xiong, C.G.; Yang, L.; Zhou, B.; Chen, Y. Finite-time fault-tolerant control of robotic systems with uncertain dynamics. Int. J. Control Autom. Syst. 2022, 20, 2681–2690. [Google Scholar] [CrossRef]
  29. Fu, C.Y.; Tian, Y.T.; Huang, H.Y.; Zhang, L.; Peng, C. Finite-time trajectory tracking control for a 12-rotor unmanned aerial vehicle with input saturation. ISA Trans. 2018, 81, 52–62. [Google Scholar] [CrossRef]
  30. Liu, K.; Wang, R.J.; Wang, X.D.; Wang, X.X. Anti-saturation adaptive finite-time neural network based fault-tolerant tracking control for a quadrotor UAV with external disturbances. Aerosp. Sci. Technol. 2021, 115, 106790. [Google Scholar] [CrossRef]
  31. Cao, N.; Lynch, A.F. Inner-outer loop control for quadrotor UAVs with input and state constraints. IEEE Trans. Contr. Syst. T 2015, 24, 1797–1804. [Google Scholar] [CrossRef]
  32. Lee, D. Fault-tolerant finite-time controller for attitude tracking of rigid spacecraft using intermediate quaternion. IEEE Trans. Aero Elec Sys. 2021, 57, 540–553. [Google Scholar] [CrossRef]
  33. Polyakov, A. Nonlinear feedback design for fixed-time stabilization of linear control system. IEEE Trans. Automat Contr. 2012, 57, 2106–2110. [Google Scholar] [CrossRef]
  34. Bernuau, E.; Efimov, D.; Perruquetti, W. On homogeneity and its application in sliding mode control. J. Franklin Inst. 2014, 351, 1866–1901. [Google Scholar] [CrossRef]
  35. Andrieu, V.; Praly, L.; Astolfi, A. Homogeneous approximation, recursive observer design, and output feedback. SIAM J. Control Optim. 2009, 47, 1814–1850. [Google Scholar] [CrossRef]
Figure 1. The framework of QUAV.
Figure 1. The framework of QUAV.
Drones 07 00344 g001
Figure 2. Control block diagram of the proposed control method.
Figure 2. Control block diagram of the proposed control method.
Drones 07 00344 g002
Figure 3. Position and speed of QUAV: (a) position x ; (b) position y ; (c) position z ; (d) speed.
Figure 3. Position and speed of QUAV: (a) position x ; (b) position y ; (c) position z ; (d) speed.
Drones 07 00344 g003
Figure 4. Three-dimensional trajectory of QUAV: (a) 3D view; (b) X-Y view.
Figure 4. Three-dimensional trajectory of QUAV: (a) 3D view; (b) X-Y view.
Drones 07 00344 g004
Figure 5. Euler angle and its tracking error: (a) Roll angle; (b) pitch angle; (c) yaw angle; (d) error.
Figure 5. Euler angle and its tracking error: (a) Roll angle; (b) pitch angle; (c) yaw angle; (d) error.
Drones 07 00344 g005aDrones 07 00344 g005b
Figure 6. Rotor speed: (a) Rotor 1; (b) rotor 2; (c) rotor 3; (d) rotor 4.
Figure 6. Rotor speed: (a) Rotor 1; (b) rotor 2; (c) rotor 3; (d) rotor 4.
Drones 07 00344 g006
Figure 7. Thrust and torque: (a) thrust; (b) torque-X; (c) torque-Y; (d) torque-Z.
Figure 7. Thrust and torque: (a) thrust; (b) torque-X; (c) torque-Y; (d) torque-Z.
Drones 07 00344 g007aDrones 07 00344 g007b
Figure 8. Observation results of lumped disturbances: (a) disturbance d x ; (b) disturbance d y ; (c) disturbance d z ; (d) disturbance D x ; (e) disturbance D y ; (f) disturbance D z .
Figure 8. Observation results of lumped disturbances: (a) disturbance d x ; (b) disturbance d y ; (c) disturbance d z ; (d) disturbance D x ; (e) disturbance D y ; (f) disturbance D z .
Drones 07 00344 g008
Figure 9. Estimation error of lumped disturbances: (a) error of dx; (b) error of d y ; (c) error of d z ; (d) error of D x ; (e) error of D y ; (f) error of D z .
Figure 9. Estimation error of lumped disturbances: (a) error of dx; (b) error of d y ; (c) error of d z ; (d) error of D x ; (e) error of D y ; (f) error of D z .
Drones 07 00344 g009
Figure 10. Tracking error of position and Euler angle: (a) error of x ; (b) error of y ; (c) error of z ; (d) error of ϕ ; (e) error of θ ; (f) error of ψ .
Figure 10. Tracking error of position and Euler angle: (a) error of x ; (b) error of y ; (c) error of z ; (d) error of ϕ ; (e) error of θ ; (f) error of ψ .
Drones 07 00344 g010
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

Shao, S.; Xu, S.; Zhao, Y.; Wu, X. Unknown Input Observer-Based Fixed-Time Trajectory Tracking Control for QUAV with Actuator Saturation and Faults. Drones 2023, 7, 344. https://doi.org/10.3390/drones7060344

AMA Style

Shao S, Xu S, Zhao Y, Wu X. Unknown Input Observer-Based Fixed-Time Trajectory Tracking Control for QUAV with Actuator Saturation and Faults. Drones. 2023; 7(6):344. https://doi.org/10.3390/drones7060344

Chicago/Turabian Style

Shao, Shikai, Shuangyin Xu, Yuanjie Zhao, and Xiaojing Wu. 2023. "Unknown Input Observer-Based Fixed-Time Trajectory Tracking Control for QUAV with Actuator Saturation and Faults" Drones 7, no. 6: 344. https://doi.org/10.3390/drones7060344

APA Style

Shao, S., Xu, S., Zhao, Y., & Wu, X. (2023). Unknown Input Observer-Based Fixed-Time Trajectory Tracking Control for QUAV with Actuator Saturation and Faults. Drones, 7(6), 344. https://doi.org/10.3390/drones7060344

Article Metrics

Back to TopTop