Next Article in Journal
Applications of GANs to Aid Target Detection in SAR Operations: A Systematic Literature Review
Next Article in Special Issue
Differentiator- and Observer-Based Feedback Linearized Advanced Nonlinear Control Strategies for an Unmanned Aerial Vehicle System
Previous Article in Journal
Aerodynamic Interaction Minimization in Coaxial Multirotors via Optimized Control Allocation
Previous Article in Special Issue
Performance Estimation of Fixed-Wing UAV Propulsion Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Nonlinear Control with Estimation of Disturbances and Parameter Uncertainties for UAVs and Integrated Brushless DC Motors

by
Claudia Verónica Vera Vaca
1,2,
Stefano Di Gennaro
3,4,
Claudia Carolina Vaca García
1 and
Cuauhtémoc Acosta Lúa
1,4,*
1
Department of Technological Sciences, University of Guadalajara, La Cienega University Center, Av. Universidad 1115, Ocotlan 47820, Jalisco, Mexico
2
Department of Industrial Engineering, The National Technological Institute of Mexico Campus Ocotlan, Av. Instituto Tecnologico S/N, Ocotlan 47829, Jalisco, Mexico
3
Department of Information Engineering, Computer Science and Mathematics, Via Vetoio, Loc. Coppito, 67100 L’Aquila, Italy
4
Center of Excellence DEWS, University of L’Aquila, Via Vetoio, Loc. Coppito, 67100 L’Aquila, Italy
*
Author to whom correspondence should be addressed.
Drones 2024, 8(9), 447; https://doi.org/10.3390/drones8090447
Submission received: 4 July 2024 / Revised: 7 August 2024 / Accepted: 19 August 2024 / Published: 30 August 2024

Abstract

:
Unmanned Aerial Vehicles (UAVs) have become increasingly prevalent in various applications, ranging from surveillance to package delivery. Achieving precise control of UAV position while enhancing robustness against uncertainties and disturbances remains a critical challenge. In this study, we propose a robust nonlinear control system for a UAV and its actuators, focusing on accurately controlling the position reference vector and improving robustness against parameter uncertainties and external disturbances. The control strategy employs two control loops: an outer loop for the UAV frame and an inner loop for the UAV actuators. The outer loop generates the required angular velocities for the actuators to follow the reference position vector using the UAV’s output and the inner loop ensures that the actuators track these angular velocity references. Both control loops utilize PI-like controllers for simplicity. The proposed system incorporates nonlinear control techniques and estimation strategies for disturbances and parameter variations, enabling dynamic adaptation to changing environmental conditions. Numerical simulations were performed using both Simulink® and the simulated PX4 Autopilot environment, showing the effectiveness of the proposed control system in achieving precise position control and robust performance for both the UAV and its actuators in the presence of uncertainties and disturbances. These results underscore the potential applicability of the control system in other UAV operational scenarios.

1. Introduction

The extensive utilization of Unmanned Aerial Vehicles (UAVs) in modern times has had a significant influence on several domains, including surveillance, exploration, package delivery, and aiding in rescue missions [1,2]. Within this particular framework, the proficient management of UAVs becomes a pivotal field of study, aiming to enhance the precision of position tracking and the resilience of the controller against unpredictable parameters and external disturbances. UAVs have transformed various industries by providing innovative solutions to traditional challenges, thus increasing operational efficiency and effectiveness [3,4,5,6].
Traditionally, linear control methods like PID and LQR have been widely adopted for their simplicity and effectiveness in basic operations. However, these methods often fall short in handling the complex dynamics and external disturbances that UAVs encounter in different scenarios [7,8,9]. Consequently, various nonlinear control strategies, including sliding mode control [10,11,12,13], backstepping [14,15,16,17,18], and adaptive control [19,20,21,22], have been explored to enhance UAV robustness and adaptability in uncertain environments.
Robust control methods that incorporate uncertainty and disturbance estimators have shown promise in maintaining control performance under varying conditions [23,24,25,26,27]. High-order sliding mode (HOSM) estimators [28,29] and adaptive dynamic programming are examples of approaches used to estimate and compensate for external disturbances and parameter uncertainties, ensuring UAVs can maintain stable flight and accurate trajectory tracking despite significant environmental and dynamic uncertainties. The necessity for such advanced control systems is underscored by their practical applications. For instance, the European project COMP4DRONES highlights the importance of robust flight control methodologies for precision agriculture, where UAVs monitor plant health and efficiently apply fertilizers [30]. Similarly, NASA’s UAS Traffic Management (UTM) project focuses on developing a traffic management system for autonomous UAV operations to ensure safe and efficient integration into the national airspace, demonstrating the need for robust and adaptive control strategies in diverse operational environments [31]. Additionally, studies have recently further examined various robust control approaches, such as robust dynamic sliding mode control [32,33,34], neural networks [35,36,37], and deep reinforcement learning [38,39], to enhance UAV performance in dynamic and uncertain environments.
Now, shifting the focus to the integration of UAV control with motor dynamics, few studies address this combined approach although these do not consider robust control. An example is the development of a velocity sensor for real-time backstepping control of a multirotor, considering actuator dynamics to enhance performance [40]. Similarly, a trajectory tracking and attitude control method for quadrotors that includes actuator dynamics, improving overall system performance, has been proposed [41]. Low-level control strategies incorporating actuator dynamics for multirotor UAVs have also demonstrated improved control accuracy and stability [42]. Further, the integration of UAV control with motor dynamics is also addressed in [43,44,45] to enhance trajectory tracking performance.
In contrast to existing research, our work differentiates by the following key contributions:
  • Detailed mathematical models of the UAV and its brushless DC motors are presented, highlighting the relationship between them and providing a robust foundation for control design.
  • A dual-loop control system is developed, with an inner loop managing the UAV’s actuators and an outer loop controlling the UAV’s frame to ensure precise position tracking through motor control.
  • A robust position control design is implemented, integrating control for the UAV and its motors, and incorporating estimation and compensation of external disturbances and parameter variations.
  • The proposed control strategies are validated through numerical simulations in both Simulink® and PX4 Autopilot environments, demonstrating effectiveness across different UAV platforms and operational scenarios.
This study aims to fill these gaps by providing a comprehensive and integrated control solution for UAVs.
The paper is organized as follows. In Section 2, the mathematical model of a UAV and their actuators is recalled. In Section 3, the problem is formulated and a controller is designed to solve this problem. In Section 4 the controller is adapted to ensure its effectiveness even with the presence of parameter variations and external disturbances, whereas in Section 5 some simulation results are given, showing the effectiveness of the proposed solution. Finally, some conclusions are given, along with some future activities in Section 6.

2. Mathematical Model

2.1. UAV Dynamic Model

A quadrotor is a four-rotor aerial vehicle that employs two coordinate systems for describing its motion: the Earth fixed frame {E} ( O E , X E , Y E , Z E ) which is affixed to the Earth’s surface and is employed to depict motion concerning the Earth’s surface, and the Body fixed frame {B} ( O B , X B , Y B , Z B ) , which is attached to the center of mass of the quadcopter’s body and describes its motion relative to itself. As can be seen in Figure 1, four identical brushless DC motors are securely mounted on the robust, cross-shaped frame of the quadcopter. Motors 1 and 2 rotate counterclockwise, while motors 3 and 4 rotate clockwise with respect to the positive Z B -axis direction at an angular velocity ω i 0 that generates forces F i 0 ( i = 1 , 2 , 3 , 4 ) [46].
The variables states used to describe the position, velocity and attitude of the quadrotor are: p v η ω , where p = [ x , y , z ] T denotes the translational position vector, and η = [ ϕ , θ , ψ ] T denotes the UAV’s attitude with respect to {E}, where ϕ [ π / 2 , π / 2 ] is the roll angle about the X E –axis, θ [ π / 2 , π / 2 ] is the pitch angle about the Y E –axis, and ψ [ π , π ] is the yaw angle about the Z E –axis. Moreover v = [ u , v , w ] T and ω = [ p , q , r ] T are the linear and angular velocities, represented in {E} and {B}, respectively [8,47,48,49].
To transform vectors expressed in {B} into vectors expressed in {E}, the rotation matrix R ( η ) and the transformation matrix W ( η ) are used, respectively
R ( η ) = cos ψ cos θ cos ψ sin θ sin ϕ sin ψ cos ϕ cos ψ sin θ cos ϕ + sin ψ sin ϕ sin ψ cos θ sin ψ sin θ sin ϕ + cos ψ cos ϕ sin ψ sin θ cos ϕ cos ψ sin ϕ sin θ cos θ sin ϕ cos θ cos ϕ
W ( η ) = 1 tan θ sin ϕ tan θ cos ϕ 0 cos ϕ sin ϕ 0 sec θ sin ϕ sec θ cos ϕ
The dynamic model is developed using the Newton-Euler method. The translational dynamics are expressed in {E}, while the rotation dynamics are expressed in {B}. Thus, the dynamic equations of motion can be represented as follows [25,50,51]
p ˙ = v v ˙ = 1 m R ( η ) F p r o p + 1 m F g r a v + 1 m F a e r o η ˙ = W ( η ) ω ω ˙ = I 1 S ω I ω + T p r o p + T g y r o + T a e r o
with m and I = d i a g ( I x x , I y y , I z z ) which represent the mass and the inertia matrix of the UAV. S ω ( ω ) is the so-called dyadic representation of ω , which can be expressed as a skew-symmetric matrix so-called dyadic representation of ω . F p r o p and T p r o p are the forces and moments generated by the propellers. F g r a v represents the force due to the gravity, F a e r o and T a e r o denote aerodynamic forces and torques acting on the UAV due to the external disturbances, and T g y r o defines the gyroscopic effects caused by propeller rotations. These variables are defined as
F p r o p = 0 0 U 1 = 0 0 i = 1 4 F i , T p r o p = U 2 U 3 U 4 = 2 2 l ( F 1 F 2 F 3 + F 4 ) 2 2 l ( F 1 F 2 + F 3 F 4 ) κ τ ( F 1 + F 2 F 3 F 4 ) , F g r a v = 0 0 m g , F a e r o = F x F y F z , T a e r o = T ϕ T θ T ψ , T g y r o = i = 1 4 ( 1 ) i I r Ω r S ω Z B
where the force produced by each motor is determined by F i = κ b ω i 2 ( i = 1 , 2 , 3 , 4 ) , with κ b representing the thrust factor. Additionally, l denotes the distance from the center of mass of the UAV to the motor shaft, κ τ stands for the drag factor, I r corresponds to the motor inertia, Ω r = ω 1 ω 2 + ω 3 + ω 4 which represents a cross configuration, and Z B = [ 0 0 1 ] T . Moreover, F j ( j = x , y , z ) denotes the aerodynamic force for each applicable vehicle orientation, computed as F j = 1 2 ρ C j V a w 2 A , where C j represents the aerodynamic coefficients, ρ stands for air density, A is the projected area of UAV element normal to the incident flow and, V a w = V V w i n d indicates the aircraft’s velocity relative to the wind [52]. Furthermore, T k ( k = ϕ , θ , ψ ) corresponds to the aerodynamic torques induces by the aerodynamic forces estimated by T a e r o = l c × F a e r o , where l c corresponds to the position vector from the center of mass to the point of application of the force.
When the quadrotor Euler angles ( ϕ , θ , ψ ) are considered small ( ϕ , θ , ψ 0 ), W ( η ) simplifies to the identity matrix I 3 x 3 . This assumption is valid because the quadcopter’s movements are gradual and gentle [53]. Therefore, the body rates equations can be approximated as [ ϕ ˙ , θ ˙ , ψ ˙ ] [ p , q , r ] . This approach is detailed in [47,48].
Consequently, under the small angle assumption, and considering (1) and (2), the quadrotor’s dynamics can be summarized as
m x ¨ = sin ψ sin ϕ + cos ψ sin θ cos ϕ U 1 + F x m y ¨ = cos ψ sin ϕ + sin ψ sin θ cos ϕ U 1 + F y m z ¨ = m g + cos θ cos ϕ U 1 + F z I x x ϕ ¨ = θ ˙ ψ ˙ ( I y y I z z ) I r θ ˙ Ω r + l U 2 + T ϕ I y y θ ¨ = ϕ ˙ ψ ˙ ( I z z I x x ) + I r ϕ ˙ Ω r + l U 3 + T θ I z z ψ ¨ = ϕ ˙ θ ˙ ( I x x I y y ) + U 4 + T ψ
where the control inputs U = [ U 1 , U 2 , U 3 , U 4 ] T are given by U = M Ω r 2 , where M is the matrix of forces and moments generated by the motors which in this UAV configuration corresponds to
M = κ b κ b κ b κ b κ p κ p κ p κ p κ p κ p κ b κ p κ τ κ τ κ τ κ τ
with a thrust factor κ b = C τ ρ D 4 and a drag factor κ τ = C D ρ A / 2 , where C τ is the propeller’s nondimensional thrust coefficient, C D is the drag coefficient, ρ is the air density, D is the propeller diameter, and A is reference area of the UAV [54]. κ p is obtained by κ b sin ( 2 / 2 ) since the rotors are placed at 45 degrees from axis X B and Y B as can be seen in Figure 1.
Remark 1.
The main advantage of this approach, which involves handling the control inputs separately from the angular velocities, is related to the practicality of the controller design. Specifically, by obtaining the control inputs and applying U = M Ω r 2 to relate them to the angular velocities, the quadrotor model remains unchanged. Only the matrix M would need to be modified to obtain the corresponding angular velocities. This method ensures flexibility and ease of implementation across various UAV configurations, allowing any designed controller to be adapted for use in other multi-rotor vehicles (with 6, 8, or different numbers of rotors) by simply changing the respective equation that links the model inputs with the rotor signals [45,55].

2.2. Actuators Dynamic Model

The motors used by the quadrotor are Brushless DC (BLDC) type. While drones typically use three-phase BLDC motors for better performance and efficiency, in this case, single-phase motors are assumed to simplify the design and reduce system complexity. This simplification allows for more straightforward modeling and control implementation, which is particularly beneficial in the initial stages of system development. The single-phase BLDC motor model reduces the number of variables and equations, making it easier to analyze and understand the fundamental dynamics of the UAV system [25,56]. These motors consist of a series circuit comprising a voltage source, resistor, inductor, and voltage drop across the motor as shown in Figure 2.
The mathematical model of a single-phase Brushless DC motor is based on electrical and mechanical equations that describe its behavior. The electrical equation used is based on Kirchhoff’s voltage law, which relates the applied voltage V, internal resistance R , inductance L, and the electromotive force E m generated by the motor, which depends on the angular velocity ω m and is related to the electromotive constant of the motor K b . On the other hand, the mechanical equation relates the external load τ l , the angular velocity of the motor ω m , and the produced torque T m which can be obtained from the motor torque constant K m multiplied by the current i. Therefore the mathematical model of the motor is summarized by
V = R i + L d dt i + K b ω m J m d dt ω m = K m i B m ω m τ l
where J m refers to motor moment of inertia and B m is the damping coefficient.

3. Design of a Nonlinear Control with Perfectly Known Parameters

The control strategy considers two control loops, an inner loop is used for the UAV’s actuators and an outer loop for the UAV frame. The overall control goal is to follow the reference vector p r e f using the vehicle’s output p. To do this, a UAV controller (outer control loop) that generates the required angular velocities for the actuators ω m , r e f is created. After then, the actuators controller (inner control loop) is made to track the references ω m , r e f using the actual angular velocities of the actuators ω m . In this case, both are PI–like controllers to simplify their operation.

3.1. UAV Controller

Applying a PI-like controller for the trajectory of a UAV involves designing a control system that allows the drone to follow a desired route in a precise and stable manner. In our control approach, we have divided the task into two fundamental stages: position and attitude control.
The position controller is defined considering the error variable e p = p p r e f , where the control problem consists of designing a controller such that lim t e p = 0 . This controller is broken down into three main components. We will call the first subsystem S 1 , which is responsible for controlling the altitude of the UAV. The second component, known as subsystem S 2 , is used for controlling the displacement on the X E -axis by imposing θ r e f to the attitude controller. Finally the third, subsystem S 3 , controls linear displacement in the Y E -axis by imposing ϕ r e f to the attitude controller.
On the other hand, the attitude controller is responsible for controlling the roll ( ϕ ), pitch ( θ ), and yaw ( ψ ) movements by considering the error variable e η = η η r e f where the controller aim is that lim t e η = 0 . From this controller we can consider a fourth subsystem called S 4 , used to exclusively control the yaw movement.
Taking the aforementioned into account and introducing some variable modifications in (3), where x = x 1 , x ˙ 1 = x 2 , y = y 1 , y ˙ 1 = y 2 , z = z 1 , z ˙ 1 = z 2 , ϕ = ϕ 1 , ϕ ˙ 1 = ϕ 2 , θ = θ 1 , θ ˙ 1 = θ 2 , ψ = ψ 1 , ψ ˙ 1 = ψ 2 , we obtain
S 1 = z ˙ 1 = z 2 z ˙ 2 = g + cos θ 1 cos ϕ 1 U 1 m + F z m S 2 = x ˙ 1 = x 2 x ˙ 2 = sin ψ 1 sin ϕ 1 + cos ψ 1 sin θ 1 cos ϕ 1 U 1 m + F x m θ ˙ 1 = θ 2 θ ˙ 2 = ϕ 2 ψ 2 ( I z z I x x ) + I r ϕ 2 Ω r + l U 3 + T θ / I y y S 3 = y ˙ 1 = y 2 y ˙ 2 = cos ψ 1 sin ϕ 1 + sin ψ 1 sin θ 1 cos ϕ 1 U 1 m + F y m ϕ ˙ 1 = ϕ 2 ϕ ˙ 2 = θ 2 ψ 2 ( I y y I z z ) I r θ 2 Ω r + l U 2 + T ϕ / I x x S 4 = ψ ˙ 1 = ψ 2 ψ ˙ 2 = ϕ 2 θ 2 ( I x x I y y ) + U 4 + T ψ / I z z
Assumption 1.
We assume that all variables for the control of the UAV are measured accurately.
This assumption is justified based on the availability and reliability of modern sensor technologies used in UAVs. GPS systems provide accurate global positioning data for outdoor environments, while motion capture and ultra-wideband (UWB) localization systems offer high-precision position measurements for indoor applications. Attitude data, including roll, pitch, and yaw angles, are obtained from Inertial Measurement Units (IMUs) that combine accelerometers and gyroscopes, with magnetometer readings often integrated through sensor fusion algorithms like the Kalman filter to enhance accuracy. Linear velocity can be derived by differentiating position data or integrating acceleration data from IMUs, and angular velocity is directly measured by gyroscopes within the IMU. The aerodynamic forces ( F x , F y , F z ) and moments ( T ϕ , T θ , T ψ ) are known from measurements obtained using pitot tubes. Pitot tubes measure the airspeed relative to the UAV, which allows the calculation of different aerodynamic forces and, subsequently, the corresponding torques due to the wind. These advanced measurement technologies ensure that the UAV’s state variables are reliably known, thereby supporting effective and stable flight control.
By using the error variable e p = p p r e f of the position controller and considering (5), we can obtain the tracking error for subsystems S 1 , S 2 , and S 3 as e z 1 = z 1 z 1 , r e f , e x 1 = x 1 x 1 , r e f and e y 1 = y 1 y 1 , r e f respectively. The dynamics of these errors are given by
e ˙ z 1 = z 2 z ˙ 1 , r e f
e ˙ x 1 = x 2 x ˙ 1 , r e f
e ˙ y 1 = y 2 y ˙ 1 , r e f
The given control problem will be solved under the following assumption.
Assumption 2.
The reference signals z 1 , r e f , x 1 , r e f and y 1 , r e f and their derivatives z ˙ 1 , r e f , x ˙ 1 , r e f and y ˙ 1 , r e f are bounded.
In the Assumption 2, the signals z 1 , r e f , x 1 , r e f and y 1 , r e f are assumed to be constant and are considered bounded since they have fixed values that do not change over time. Consequently, their derivatives z ˙ 1 , r e f , x ˙ 1 , r e f and y ˙ 1 , r e f are also considered bounded.
Without loss of generality and following the same order of subsystems in (5), the virtual controls k z 1 e z 1 k 0 , z 1 I e z 1 = z 2 , d z ˙ 1 , r e f ; k x 1 e x 1 k 0 , x 1 I e x 1 = x 2 , d x ˙ 1 , r e f and k y 1 e y 1 k 0 , y 1 I e y 1 = y 2 , d y ˙ 1 , r e f are used, where I e j is an integrative term defined by I e ˙ j = e j , k j 1 corresponds to the proportional gain and k 0 , j 1 is the integral gain with j = z , x , y . In this manner it is feasible to get z 2 , d , x 2 , d and y 2 , d as follows
z 2 , d = z ˙ 1 , r e f k z 1 e z 1 k 0 , z 1 I e z 1
x 2 , d = x ˙ 1 , r e f k x 1 e x 1 k 0 , x 1 I e x 1
y 2 , d = y ˙ 1 , r e f k y 1 e y 1 k 0 , y 1 I e y 1
Now we can define the errors e z 2 = z 2 z 2 , d , e x 2 = x 2 x 2 , d and e y 2 = y 2 y 2 , d and upon derivation, we obtain e ˙ z 2 = z ˙ 2 z ˙ 2 , d , e ˙ x 2 = x ˙ 2 x ˙ 2 , d and e ˙ y 2 = y ˙ 2 y ˙ 2 , d . Therefore, the dynamics of the tracking error of S 1 , S 2 , and S 3 can be expressed as
e ˙ z 2 = g + 1 m cos θ 1 cos ϕ 1 U 1 + F z m z ˙ 2 , d
e ˙ x 2 = 1 m ( sin ψ 1 sin ϕ 1 + cos ψ 1 sin θ 1 cos ϕ 1 ) U 1 + F x m x ˙ 2 , d
e ˙ y 2 = 1 m ( cos ψ 1 sin ϕ 1 + sin ψ 1 sin θ 1 cos ϕ 1 ) U 1 + F y m y ˙ 2 , d
where z ˙ 2 , d = z ¨ 1 , r e f k z 1 e ˙ z 1 k 0 , z 1 I ˙ e z 1 , x ˙ 2 , d = x ¨ 1 , r e f k x 1 e ˙ x 1 k 0 , x 1 I ˙ e x 1 and y ˙ 2 , d = y ¨ 1 , r e f k y 1 e ˙ y 1 k 0 , y 1 I ˙ e y 1 .
However, in order to achieve lim t e x = 0 and lim t e y = 0 , which ensures accurate tracking of translation along the x and y axes, it is essential to control the inclination angles pitch ( θ 1 ) and roll ( ϕ 1 ) respectively. In this manner, by the stabilization of e z 1 , e x 1 and e y 1 , the reference values for sin θ 1 and sin ϕ 1 can be fixed as follows
sin θ 1 , d = m cos ϕ 1 cos ψ 1 U 1 ( x ˙ 2 , d 1 m ( sin ϕ 1 sin ψ 1 ) U 1 F x m k x 2 e x 2 k 0 , x 2 I e x 2 )
sin ϕ 1 , d = m cos ψ 1 U 1 ( y ˙ 2 , d 1 m ( cos ϕ 1 sin θ 1 sin ψ 1 ) U 1 F y m k y 2 e y 2 k 0 , y 2 I e y 2 )
In this way, to impose the references, θ 1 and ϕ 1 are now treated as tracking errors as depicted below
e ϕ 1 = sin ϕ 1 sin ϕ 1 , d
e θ 1 = sin θ 1 sin θ 1 , d
Here, sin ϕ 1 , d and sin θ 1 , d represent the desired values of the sine of the roll and pitch angles, respectively. Hence ϕ 1 , d = arcsin ϕ 1 , d and θ 1 , d = arcsin θ 1 , d .
Through derivation of (17) and (18), and by applying the chain rule of differentiation, we obtain the dynamics of the errors which are given by
e ˙ ϕ 1 = cos ϕ 1 ϕ 2 d dt sin ϕ 1 , d
e ˙ θ 1 = cos θ 1 θ 2 d dt sin θ 1 , d
By the utilization of the virtual controls k ϕ 1 e ϕ 1 k 0 , ϕ 1 I e ϕ 1 = cos ϕ 1 ϕ 2 , d d dt sin ϕ 1 , d for (17), and k θ 1 e θ 1 k 0 , θ 1 I e θ 1 = cos θ 1 θ 2 , d d dt sin θ 1 , d for (18), where I e k is an integrative term defined by I e ˙ k = e k , k k 1 corresponds to the proportional gain and k 0 , k 1 is the integral gain with k = ϕ 1 , θ 1 , it is feasible to get ϕ 2 , d and θ 2 , d as follows
ϕ 2 , d = 1 cos ϕ 1 d dt sin ϕ 1 , d k ϕ 1 e ϕ 1 k 0 , ϕ 1 I e ϕ 1
θ 2 , d = 1 cos θ 1 d dt sin θ 1 , d k θ 1 e θ 1 k 0 , θ 1 I e θ 1
Now, by defining the tracking errors e ϕ 2 = ϕ 2 ϕ 2 , d , e θ 2 = θ 2 θ 2 , d and their respective derivation e ˙ ϕ 2 = ϕ ˙ 2 ϕ ˙ 2 , d , e ˙ θ 2 = θ ˙ 2 θ ˙ 2 , d , the error dynamics are obtained
e ˙ ϕ 2 = I y y I z z I x x θ 2 ψ 2 I r I x x θ 2 Ω r + l I x x U 2 + T ϕ I x x ϕ ˙ 2 , d
e ˙ θ 2 = I z z I x x I y y ϕ 2 ψ 2 + I r I y y ϕ 2 Ω r + l I y y U 3 + T θ I y y θ ˙ 2 , d
On the other hand, through the utilization of the error variable e η = η η r e f of the attitude controller and considering (5), we can derive the tracking error for the subsystem S 4 as e ψ 1 = ψ 1 ψ 1 , r e f . The dynamics of this error is given by e ˙ ψ 1 = ψ 2 ψ ˙ 1 , r e f . Hence, through the utilization of the virtual control k ψ 1 e ψ 1 k 0 , ψ 1 I e ψ 1 = ψ 2 , d ψ ˙ 1 , r e f , the dynamics of the tracking error in the subsystem S 4 can be seen as follows
e ˙ ψ 2 = I x x I y y I z z ϕ 2 θ 2 + U 4 I z z + T ψ I z z ψ ˙ 2 , d
where ψ ˙ 2 , d = ψ ¨ 1 , r e f k ψ 1 e ˙ ψ 1 k 0 , ψ 1 I e ˙ ψ 1 .
Therefore, based on (12), (23), (24), and (25), the proposed control inputs are
U 1 = m c o s θ 1 c o s ϕ 1 z ˙ 2 , d + g F z m k z 2 e z 2 k 0 , z 2 I e z 2 U 2 = I x x l ϕ ˙ 2 , d θ 2 ψ 2 ( I y y I z z ) I x x + I r θ 2 Ω r I x x T ϕ I x x k ϕ 2 e ϕ 2 k 0 , ϕ 2 I e ϕ 2 U 3 = I y y l θ ˙ 2 , d ϕ 2 ψ 2 ( I z z I x x ) I y y I r ϕ 2 Ω r I y y T θ I y y k θ 2 e θ 2 k 0 , θ 2 I e θ 2 U 4 = I z z ψ ˙ 2 , d ϕ 2 θ 2 ( I x x I y y ) I z z T ψ I z z k ψ 2 e ψ 2 k 0 , ψ 2 I e ψ 2
where k i , k 0 , i > 0 , i = z 2 , ϕ 2 , θ 2 , ψ 2 .
To establish a strong connection between the UAV controller and the actuators, it is essential to convert the control signals produced by (26) into reference angular velocities that the actuators can interpret and execute. This transformation is accomplished through the resolution of Ω r from U = M Ω r 2 .

3.2. Actuators Controller

As mentioned in Section 2.1, the quadcopter is equipped with four identical motors, operating according to the dynamics defined by Equation (4).
Assumption 3.
We assume that all variables for the control of the BLDC motor are measured accurately.
This assumption is justified based on the advanced measurement technologies and control systems available for BLDC motors. The rotor speed ( ω m , i ) and load torque ( τ l , i ) are typically measured using rotary encoders and torque sensors, respectively. The current ( i i ) is measured using current sensors such as Hall effect sensors. The inertia ( J m , i ), damping coefficient ( B m , i ), motor constant ( K m , i ), back-EMF constant ( K b , i ), resistance ( R i ), and inductance ( L i ) are known from the motor specifications provided by the manufacturer. The voltage ( V i ) is measured using voltage sensors. These technologies and methods ensure that all relevant variables are accurately known, supporting precise and effective control of the BLDC motor.
In this framework, a PI–like controller for the actuators is formulated to guarantee that the angular velocities of the motors, ω m , i , i = 1 , , 4 , closely track the desired references. Therefore, lim t e ω , m , i = 0 , considering the error variable e ω , m , i = ω m , i ω m , i , r e f .
Therefore, from (4) it is possible to obtain the error dynamics as follows
e ˙ ω , m , i = K m , i J m , i i i B m , i J m , i ω m , i τ l , i J m , i d dt ω m , i , r e f
where it is necessary to impose the current by
i i , r e f = J m , i K m , i ( B m , i ω m , i J m , i + τ l , i J m , i + ω ˙ i , r e f k ω , i e ω , i k 0 , ω , i I e ω , i )
with k ω , i , k 0 , ω , i > 0 , where the PI–like controller have been applied. Finally, it is possible to get the current error dynamics by considering the current error e i , i = i i i i , r e f , such that lim t e i , i = 0
e ˙ i , i = 1 L i ( R i i i K b , i ω m , i L i + V i ) d dt i i , r e f
from where it is possible to derive the voltage input which has to be applied to the motor
V i = L i ( R i i i L i + K b , i ω m , i L i + d dt i i , r e f k i , i e i , i k 0 , i , i I e i , i )
with k i , i , k 0 , i , i > 0 .

4. Design of a Robust Nonlinear Control with Estimation of Disturbances and Parameter Uncertainties

In this section, a robust control is designed to address uncertainties in system parameters and the presence of external disturbances. The proposed control strategy aims to enhance the stability and performance of the UAV despite variations in system dynamics and external factors that could potentially disrupt its operation.
The robust control framework integrates nonlinear control techniques and disturbance estimation strategies to dynamically adapt to changing environmental conditions. By incorporating PI-like controllers within both the outer and inner control loops, the system effectively mitigates the impact of uncertainties and disturbances on the UAV’s performance.
A general schematic of the proposed robust control system is shown in Figure 3. This diagram illustrates the interaction between the various components of the control architecture, including the measurement of state variables, the generation of control inputs, and the implementation of feedback mechanisms to ensure precise trajectory tracking and stable flight performance.

4.1. UAV Robust Control with Disturbance Estimation

To achieve robustness of the closed–loop system with respect to parameter variations and external disturbances, the dynamics of the tracking errors of the system (12)–(14), and (23)–(25) are rewritten
e ˙ z 2 = g + 1 m ( cos θ 1 cos ϕ 1 ) U 1 + F z m z ˙ 2 , d + δ z 1 e ˙ x 2 = 1 m ( sin ψ 1 sin ϕ 1 + cos ψ 1 sin θ 1 cos ϕ 1 ) U 1 + F x m x ˙ 2 , d + δ x 1 e ˙ y 2 = 1 m ( sin ψ 1 sin θ 1 cos ϕ 1 cos ψ 1 sin ϕ 1 ) U 1 + F y m y ˙ 2 , d + δ y 1 e ˙ ϕ 2 = I y y I z z I x x θ 2 ψ 2 I r I x x θ 2 Ω r + l I x x U 2 + T ϕ I x x ϕ ˙ 2 , d + δ ϕ 1 e ˙ θ 2 = I z z I x x I y y ϕ 2 ψ 2 + I r I y y ϕ 2 Ω r + l I y y U 3 + T θ I y y θ ˙ 2 , d + δ θ 1 e ˙ ψ 2 = I x x I y y I z z ϕ 2 θ 2 + U 4 I z z + T ψ I z z ψ ˙ 2 , d + δ ψ 1
with
δ z 1 = 1 m 1 m cos θ 1 cos ϕ 1 U 1 + F z m F z m δ x 1 = 1 m 1 m sin ψ 1 sin ϕ 1 + cos ψ 1 sin θ 1 cos ϕ 1 U 1 + F x m F x m δ y 1 = 1 m 1 m sin ψ 1 sin θ 1 cos ϕ 1 cos ψ 1 sin ϕ 1 U 1 + F y m F y m δ ϕ 1 = I y y I z z I x x I y y I z z I x x θ 2 ψ 2 I r I x x I r I x x θ 2 Ω r + l I x x l I x x U 2 + T ϕ I x x T ϕ I x x δ θ 1 = I z z I x x I y y I z z I x x I y y ϕ 2 ψ 2 + I r I y y I r I y y ϕ 2 Ω r + l I x x l I y y U 3 + T θ I y y T θ I y y δ ψ 1 = I x x I y y I z z I x x I y y I z z ϕ 2 θ 2 + 1 I z z 1 I z z U 4 + T ψ I z z T ψ I z z
Remark 2.
Within the vectors δ j and δ k there are the parametric variations and the external disturbances.
The following statement solves the posed control problem with estimation the variation of external disturbances, where · 1 / 2 = | · | 1 / 2 sign ( · ) , · 0 = sign ( · )
Theorem 1.
Under Assumptions 1 and 2, the robust control
ξ ˙ z 1 = λ z 1 ξ z 1 e z 1 1 / 2 + ξ z 2 + 1 m cos θ 1 cos ϕ 1 U 1 + F z m g z ˙ 2 , d ξ ˙ z 2 = λ z 2 ξ z 1 e z 1 0 ξ ˙ x 1 = λ x 1 ξ x 1 e x 1 1 / 2 + ξ x 2 + 1 m sin ψ 1 sin ϕ 1 + cos ψ 1 sin θ 1 cos ϕ 1 U 1 + F x m x ˙ 2 , d ξ ˙ x 2 = λ x 2 ξ x 1 e x 1 0 ξ ˙ y 1 = λ y 1 ξ y 1 e y 1 1 / 2 + ξ y 2 + 1 m sin ψ 1 sin θ 1 cos ϕ 1 cos ψ 1 sin ϕ 1 U 1 + F y m y ˙ 2 , d ξ ˙ y 2 = λ y 2 ξ y 1 e y 1 0 ξ ˙ ϕ 1 = λ ϕ 1 ξ ϕ 1 e ϕ 1 1 / 2 + ξ ϕ 2 + I y y I z z I x x θ 2 ψ 2 I r I x x θ 2 Ω r + l I x x U 2 + T ϕ I x x ϕ ˙ 2 , d ξ ˙ ϕ 2 = λ ϕ 2 ξ ϕ 1 e ϕ 1 0 ξ ˙ θ 1 = λ θ 1 ξ θ 1 e θ 1 1 / 2 + ξ θ 2 + I z z I x x I y y ϕ 2 ψ 2 + I r I y y ϕ 2 Ω r + l I y y U 3 + T θ I y y θ ˙ 2 , d ξ ˙ θ 2 = λ θ 2 ξ θ 1 e θ 1 0 ξ ˙ ψ 1 = λ ψ 1 ξ ψ 1 e ψ 1 1 / 2 + ξ ψ 2 + I x x I y y I z z ϕ 2 θ 2 + 1 I z z U 4 + T ψ I z z ψ ˙ 2 , d ξ ˙ ψ 2 = λ ψ 1 ξ ψ 1 e ψ 2 0 U 1 = m c o s θ 1 c o s ϕ 1 z ˙ 2 , d + g F z m k z 2 e z 2 k 0 , z 2 I e z 2 ξ z 2 U 2 = I x x l ( ϕ ˙ 2 , d θ 2 ψ 2 ( I y y I z z ) I x x + I r θ 2 Ω r I x x T ϕ I x x k ϕ 2 e ϕ 2 k 0 , ϕ 2 I e ϕ 2 ξ ϕ 2 ) U 3 = I y y l ( θ ˙ 2 , d ϕ 2 ψ 2 ( I z z I x x ) I y y I r ϕ 2 Ω r I y y T θ I y y k θ 2 e θ 2 k 0 , θ 2 I e θ 2 ξ θ 2 ) U 4 = I z z ψ ˙ 2 , d ϕ 2 θ 2 ( I x x I y y ) I z z T ψ I z z k ψ 2 e ψ 2 k 0 , ψ 2 I e ψ 2 ξ ψ 2 I ˙ e z 2 = e z 2 I ˙ e ϕ 2 = e ϕ 2 I ˙ e θ 2 = e θ 2 I ˙ e ψ 2 = e ψ 2
ensures that the tracking error e z 2 , e ϕ 2 , e θ 2 , e ψ 2 asymptotically converge to zero and the estimation errors ϵ x 2 = ξ x 2 δ x 1 , ϵ y 2 = ξ y 2 δ y 1 , ϵ z 2 = ξ z 2 δ z 1 , ϵ ϕ 2 = ξ ϕ 2 δ ϕ 1 , ϵ θ 2 = ξ θ 2 δ θ 1 and ϵ ψ 2 = ξ ψ 2 δ ψ 1 converge to zero in finite–time and the quantities δ j and δ k are bounded by finite positive constants Δ j , m a x and Δ k , m a x for
λ j , 1 > 0 , λ j , 2 > 3 Δ j , m a x + 2 Δ j , m a x 2 / λ j , 1 2 λ k , 1 > 0 , λ k , 2 > 3 Δ k , m a x + 2 Δ k , m a x 2 / λ k , 1 2 k j 2 > 0 , k 0 , j 2 > 0 , k k 2 > 0 , k 0 , k 2 > 0
for j = x , y , z and k = ϕ , θ , ψ , and with Δ j , m a x , Δ k , m a x R appropriate bounds.
Proof. 
Considering the decoupled subsystem
e ˙ z 1 = λ z 1 e z 1 1 / 2 + e z 2 e ˙ z 2 = λ z 2 e z 1 0 δ ˙ z
Let us consider the following Lyapunov candidate
V z = 1 2 E z T P z E z
with
P z = λ z 1 2 + 4 λ z 2 λ z 1 λ z 1 2 = P z T > 0 , E z = e z 1 1 / 2 e z 2
E ˙ z = 1 2 | e z 1 | 1 / 2 λ z 1 E z δ ˙ z 2 | e z 1 | 1 / 2 e z 1 0 λ z 0 E z , λ z 1 = λ z 1 1 2 λ z 2 0 , λ z 0 = 0 0 2 0
Therefore,
V ˙ z = λ z 1 2 | e z 1 | 1 / 2 E z T Q z 1 E z δ ˙ z 2 | e z 1 | 1 / 2 e z 1 0 E z T Q z 0 E z
where
Q z 1 = P 1 λ z 1 + λ z 0 T P z 2 λ z 1 = λ z 1 2 + 2 λ z 2 λ z 1 λ z 1 1 Q z 0 = P z λ z 0 = 2 λ z 1 2 2 0
One can observe that δ ˙ z in (35) is bounded. It is readily checked that, since the derivatives of the external perturbations are assumed bounded, δ ˙ z is bounded for any finite time interval. Let Δ z , max be such that | δ ˙ z | Δ z , max for any finite time interval. To take into account the perturbation δ ˙ z affecting the system (35), V ˙ z can be bounded as follows
V ˙ z = λ z 1 2 | e z 1 | 1 / 2 E z T Q z 1 E z δ ˙ z 2 | e z 1 | 1 / 2 e z 1 0 E z T Q z 0 E z λ z 1 2 | e z 1 | 1 / 2 E z T Q z 1 E z Δ z , max 2 | e z 1 | 1 / 2 E z T Q z 0 E z = λ z , 1 2 | e z 1 | 1 / 2 E z T Q z E z λ z 1 2 | e z 1 | 1 / 2 λ min Q z E z 2 2 γ z V z 1 / 2 , γ z = λ z 1 λ min P z 2 λ max P z λ min Q z E z = | e z 1 1 / 2 | | e z 2 | Q z = λ z 1 2 + 2 λ z 2 2 Δ z , max λ z 1 + 2 Δ z , max λ z 1 λ z 1 + 2 Δ z , max λ z 1 1
with λ z 1 > 0 , λ z 2 > 0 , Q z > 0 , and E z 2 2 = E z 2 2 = | e z 1 1 / 2 | 2 + | e z 2 | 2 = | e z 1 | + e z 2 2 . Hence, the boundedness of the signals involved for any finite time interval, and the fact that V ˙ z γ z V z 1 / 2 , ensure that one gets the finite–time convergence of E z to the origin.
Analogously, the following subsystems
e ˙ x 1 = λ x 1 e x 1 1 / 2 + e x 2 e ˙ x 2 = λ x 2 e x 1 0 δ ˙ x e ˙ y 1 = λ y 1 e y 1 1 / 2 + e y 2 e ˙ y 2 = λ y 2 e y 1 0 δ ˙ y e ˙ ϕ 1 = λ ϕ 1 e ϕ 1 1 / 2 + e ϕ 2 e ˙ ϕ 2 = λ ϕ 2 e ϕ 1 0 δ ˙ ϕ e ˙ θ 1 = λ θ 1 e θ 1 1 / 2 + e θ 2 e ˙ θ 2 = λ θ 2 e θ 1 0 δ ˙ θ e ˙ ψ 1 = λ ψ 1 e ψ 1 1 / 2 + e ψ 2 e ˙ ψ 2 = λ ψ 2 e ψ 1 0 δ ˙ ψ
can be studied in a similar way, using the following Lyapunov candidates
V x = 1 2 E x T P x E x , P x = λ x 1 2 + 4 λ x 2 λ x 1 λ x 1 2 = P x T > 0 , E x = e x 1 1 / 2 e x 2 V y = 1 2 E y T P y E y , P y = λ y 1 2 + 4 λ y 2 λ y 1 λ y 1 2 = P y T > 0 , E y = e y 1 1 / 2 e y 2 V ϕ = 1 2 E ϕ T P ϕ E ϕ , P ϕ = λ ϕ 1 2 + 4 λ ϕ 2 λ ϕ 1 λ ϕ 1 2 = P ϕ T > 0 , E ϕ = e ϕ 1 1 / 2 e ϕ 2 V θ = 1 2 E θ T P θ E θ , P θ = λ θ 1 2 + 4 λ θ 2 λ θ 1 λ θ 1 2 = P θ T > 0 , E θ = e θ 1 1 / 2 e θ 2 V ψ = 1 2 E ψ T P ψ E ψ , P ψ = λ ψ 1 2 + 4 λ ψ 2 λ ψ 1 λ ψ 1 2 = P ψ T > 0 , E ψ = e ψ 1 1 / 2 e ψ 2
with Δ x , max , Δ y , max , Δ ϕ , max , Δ θ , max and Δ ψ , max the bound for | δ ˙ x | , | δ ˙ y | , | δ ˙ ϕ | , | δ ˙ θ | and | δ ˙ ψ | for any finite time interval. □

4.2. Closed Loop Stability in UAV

In this section the stability analysis of the proposed nonlinear robust control scheme is performed. However, according to Theorem 1, it can be demonstrated that the variables ξ x 2 , ξ y 2 , ξ z 2 , ξ ϕ 2 , ξ θ 2 , ξ ψ 2 tends to δ x 1 , δ y 1 , δ z 1 , δ ϕ 1 , δ θ 1 , δ ψ 1 in finite–time.

4.2.1. Subsystems S1 and S4

To perform the closed-loop analysis for subsystem S1, which controls the altitude of the UAV, we consider Equation (6) and the dynamics of the tracking errors with parameter variations and external disturbances e ˙ z 2 from (31) along with the control input U 1 in (33). The closed-loop system can be represented as
( I ˙ e z 1 e ˙ z 1 I ˙ e z 2 e ˙ z 2 ) = ( 0 1 0 0 k 0 , z 1 k z 1 0 1 0 0 0 1 0 0 k 0 , z 2 k z 2 ) ( I e z 1 e z 1 I e z 2 e z 2 )
Analyzing (36), we can obtain the eigenvalues using Maple, a symbolic and numeric computing software. The eigenvalues are
ρ z 1 = k z 1 / 1 ± k z 1 2 4 k 0 , z 1 / 2 ρ z 2 = k z 2 / 2 ± k z 2 2 4 k 0 , z 2 / 2
which shows that all eigenvalues of the matrix (36) that satisfy Re { ρ z 1 , ρ z 2 } < 0 for any values of k 0 , z 1 > 0 , k z 1 > 0 , k 0 , z 2 > 0 , k z 2 > 0 , i.e., that the matrix is Hurwitz [57]. Therefore, the tracking error e z 1 , e z 2 are asymptotically stable.
The stability analysis of Subsystem S4, which is used to control the yaw movement, is performed similarly.

4.2.2. Subsystems S2 and S3

Similarly to Section 4.2.1, the analysis for the translational motion in X E -axis described in S 2 of (5) is performed as follows
( I ˙ e x 1 e ˙ x 1 I ˙ e x 2 e ˙ x 2 I ˙ e θ 1 e ˙ θ 1 I ˙ e θ 2 e ˙ θ 2 ) = ( 0 1 0 0 0 0 0 0 k 0 , x 1 k x 1 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 k 0 , x 2 k x 2 0 cos ϕ 1 cos ψ 1 U 1 m 0 0 0 0 0 0 0 1 0 0 0 0 0 0 k 0 , θ 1 k θ 1 0 cos θ 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 k 0 , θ 2 k θ 2 ) ( I e x 1 e ˙ x 1 I e x 2 e ˙ x 2 I e θ 1 e ˙ θ 1 I e θ 2 e ˙ θ 2 )
Analyzing the matrix (38), we can obtain the eigenvalues using Maple, a symbolic and numeric computing software. The eigenvalues are
ρ x 1 = k x 1 / 2 ± k x 1 2 4 k 0 , x 1 / 2 ρ x 2 = k x 2 / 2 ± k x 2 2 4 k 0 , x 2 / 2 ρ θ 1 = k θ 1 / 2 ± k θ 1 2 4 k 0 , θ 1 / 2 ρ θ 2 = k θ 2 / 2 ± k θ 2 2 4 k 0 , θ 2 / 2
Analogously of Section 4.2.1, all eigenvalues of the matrix (38) satisfy Re { ρ x 1 , ρ x 2 , ρ θ 1 , ρ θ 2 } < 0 for any values of k x 1 > 0 , k x 2 > 0 , k 0 , x 1 > 0 , k 0 , x 2 > 0 , k θ 1 > 0 , k θ 2 > 0 , k 0 , θ 1 > 0 , k 0 , θ 2 > 0 , i.e., that the matrix is Hurwitz [57]. Therefore, the tracking error e x 1 , e x 2 , e θ 1 , e θ 2 are asymptotically stable.
The stability analysis of Subsystem S3, which is responsible for controlling the linear displacement along the Y E -axis, is performed similarly to Subsystem S2.

4.3. Actuators Robust Control with Disturbance Estimation

In the same way as the system described at Section 4.1, to achieve robustness of the closed–loop system (4) with respect to parameter variations and external disturbances, the dynamics of the tracking errors of the actuator system, (27) and (29) are rewritten
e ˙ ω , i = K m , i J m , i i i B m , i J m , i ω m , i τ l , i J m , i d dt ω i , r e f + δ ω , i e ˙ i , i = 1 L i R i i i K b , i ω m , i + V i d dt i i , r e f + δ i , i
considering
δ ω , i = K m , i J m , i K m , i J m , i i i B m , i J m , i B m , i J m , i ω m , i τ l , i J m , i τ l , i J m , i δ i , i = R i L i R i L i i i K b , i L i K b , i L i ω m , i + 1 L i 1 L i V i
To estimate the variation of parameters in this system, we have that
Theorem 2.
Under Assumption 3, the robust control is proposed as
ξ ˙ ω , 1 , i = λ ω , 1 , i ξ ω , 1 , i e ω , i 1 / 2 + ξ ω , 2 , i + K m , i J m , i i i B m , i J m , i ω m , i τ l , i J m , i d dt ω i , r e f ξ ˙ ω , 2 , i = λ ω , 2 , i ξ ω , 1 , i e ω i 0 ξ ˙ i , 1 , i = λ i , 1 , i ξ i , 1 , i e i , i 1 / 2 + ξ i , 2 , i R i L i i i K b , i L i ω m , i + V i L i d dt i i , r e f ξ ˙ i , 2 , i = λ i , 2 , i ξ i , 1 , i e i , i 0 i i , r e f = J m , i K m , i ( B m , i J m , i ω m , i + τ l , i J m , i + d dt ω i , r e f k ω , i e ω , i k 0 , ω , i I e ω , i ξ ω , 1 , i ) V i = L i R i L i i i + K b , i L i ω m , i + d dt i i , r e f k i , i e i , i k 0 , i , i I e i , i ξ i , 2 , i I ˙ e ω , i = e ω , i I ˙ e i , i = e i , i
ensures that the tracking error e ω , i , e i , i asymptotically converge to zero and that the estimation errors ξ ω , 2 , i δ ω , i and ξ i , 2 , i δ i , i converge to zero in finite–time and the quantities δ ω , i and δ i , i are bounded by finite positive constants Δ ω , i , m a x and Δ i , i , m a x for
λ ω , 1 , i > 0 , λ ω , 2 , i > 3 Δ ω , i , m a x + 2 Δ ω , i , m a x 2 / λ ω , 1 , i 2 λ i , 1 , i > 0 , λ i , 2 , i > 3 Δ i , i , m a x + 2 Δ i , i , m a x 2 / λ i , 1 , i 2 k ω , i > 0 , k 0 , ω , i > 0 , k i , i > 0 , k 0 , i , i > 0 .
for i = 1 , 2 , 3 , 4 and with Δ ω , i , m a x , Δ i , i , m a x R appropriate bounds.
Proof. 
Let consider the following decoupled subsystems
e ˙ ω , 1 , i = λ ω , 1 , i e ω , 1 , i 1 / 2 + e ω , 2 , i e ˙ ω , 2 , i = λ ω , 2 , i e ω , 1 , i 0 δ ˙ ω , i e ˙ i , 1 , i = λ i , 1 , i e i , 1 , i 1 / 2 + e i , 2 , i e ˙ i , 2 , i = λ i , 2 , i e i , 1 , i 0 δ ˙ i , i
which can be studied, using the following Lyapunov candidates
V ω , i = 1 2 E ω , i T P ω , i E ω , i P ω , i = λ ω , 1 , i 2 + 4 λ ω , 2 , i λ ω , 1 , i λ ω , 1 , i 2 = P ω , i T > 0 , E ω , i = e ω , 1 , i 1 / 2 e ω , 2 , i V i , i = 1 2 E i , i T P i , i E i , i P i , i = λ i , 1 , i 2 + 4 λ i , 2 , i λ i , 1 , i λ i , 1 , i 2 = P i , i T > 0 , E i , i = e i , 1 , i 1 / 2 e i , 2 , i
Therefore, as previous proof, one can observe that δ ˙ ω and δ ˙ i in (44) are bounded. It can be readily verified that, given the derivatives of the external perturbations are also assumed to be bounded, δ ˙ ω and δ ˙ i remain bounded over any finite time interval. □

4.4. Closed Loop Stability in Actuators

In this section the stability analysis of the proposed nonlinear robust control scheme is performed. However, according to Theorem 2, it can be demonstrated that the variables ξ ω , 2 , i , and ξ i , 2 , i , tends to δ ω , i and δ i , i in finite–time.
To perform the closed-loop analysis for the control of the actuators, we consider Equations (27) and (29) and the dynamics of the tracking errors with parameter variations and external disturbances in (32) along with the control input V i in (42). The closed-loop system can be represented as
( I ˙ e ω , i e ˙ ω , i I ˙ e i , i e ˙ i , i ) = ( 0 1 0 0 k 0 , ω , i k ω , i 0 K m , i J m , i 0 0 0 1 0 0 k 0 , i , i k i , i ) ( I e ω , i e ω , i I e i , i e i , i )
Analyzing (45), we can obtain the eigenvalues using Maple, a symbolic and numeric computing software. The eigenvalues are
ρ ω , i = k ω , i / 2 ± k ω , i 2 4 k 0 , ω , i / 2 ρ i , i = k i , i / 2 ± k i , i 2 4 k 0 , i , i / 2
which shows that all eigenvalues of the matrix (45) that satisfy Re { ρ ω , i , ρ i , i } < 0 for any values of k 0 , ω , i > 0 , k ω , i > 0 , k 0 , i , i > 0 , k i , i > 0 , i.e., that the matrix is Hurwitz [57]. Therefore, the tracking error e z 1 , e z 2 are asymptotically stable.

5. Simulations Results

This section presents the simulation results of the Robust Controller applied to a UAV with integrated Brushless DC motors, in the presence of parametric uncertainties and external disturbances. Numerical simulations were performed with Simulink® and the simulated PX4 Autopilot environment [58].
Simulink®, developed by MathWorks, is a simulation and model-based design environment for dynamic and embedded systems. It provides an interactive graphical interface and customizable set of block libraries that allow for the modeling, simulation, and analysis of multidomain dynamic systems. Simulink is particularly well-suited for designing control systems due to its ability to model complex systems, run simulations to predict system behavior, and analyze performance. The use of Simulink in this study allows for the comprehensive testing of the control algorithms in a simulated environment before implementation, ensuring that the algorithms perform as expected under various conditions and disturbances.
The algorithms were then tested in the simulated PX4 Autopilot environment via numerical simulations to confirm their performance under real-world settings before being put into use. PX4 is an open-source flight control software designed for controlling autonomous vehicles, including drones. It is widely used in both research and commercial applications due to its flexibility, robustness, and compatibility with various hardware platforms. PX4 supports numerous flight modes and features, making it an ideal platform for testing and validating flight control algorithms. By using the simulated PX4 Autopilot environment, we can evaluate how the control algorithms could perform on actual UAV hardware, providing a realistic assessment of their effectiveness and robustness in handling parametric uncertainties and external disturbances.
Also, the combination of Simulink® for initial design and simulation, followed by validation in the simulated PX4 Autopilot environment, ensures a thorough evaluation of the control algorithms. This approach allows for the detection and resolution of potential issues in the simulation phase, reducing the risk of failure in real-world applications. The dual-phase testing strategy enhances the reliability and robustness of the control system, making it suitable for deployment in various UAV applications. The specifications of the UAV used are: a mass (m) of 1.3 kg, an arm length (l) of 0.175 m, a thrust coefficient ( C τ ) of 0.0048, a drag coefficient ( C D ) of 0.0002351, an air density ( ρ ) of 1.225 kg/m3, a propeller diameter (D) of 0.24 m, moments of inertia about the X B -axis ( I x x ) and Y B -axis ( I y y ) both equal to 0.081 kg·m2, and a moment of inertia about the vertical axis ( I z z ) of 0.142 kg·m².
Concerning the motors, their parameters include: a back electromotive force coefficient ( K b ) of 0.01038 Vs/rad, a torque coefficient ( K m ) of 0.04 N·m, an electrical resistance (R) of 0.2 Ω , a damping coefficient ( B m ) of 0.0002 N·m·s/rad, a motor inertia ( J m ) of 0.0000049 kg·m2, and an inductance (L) of 0.7 mH.

5.1. Environmental Disturbances Acting on the Quadrotor

The principal environmental disturbance acting on the quadrotor is the one applied by the wind. Considering its behavior, wind can be classified into two principal categories: constant wind and turbulent wind. Unfortunately, in many scenarios, the induced speed within a particular wind field is not uniform. For this reason, in our model, variable wind gusts are modeled as lateral wind blasts. These gusts are characterized by their velocity v W = ( v W , X , v W , Y ) T in the ground-fixed ( X E , Y E ) frame. The wind velocity components are modeled as
v W , X = V w cos ( β 0 + Δ β 0 ) + 0.025 N v W , Y = V w sin ( β 0 + Δ β 0 ) + 0.025 N
where V w = V w , 0 + Δ V w sin ( ω w ) , with V w , 0 as the nominal wind magnitude, β 0 is the nominal angle between the X E -axis and the wind vector, and N is a uniform distribution. Additionally, Δ β 0 = 0.2 β 0 and Δ V w = 0.1 V w , 0 .
In the vehicle-fixed ( X B , Y B ) frame, the wind velocity components are
v w , x = v W , X cos α z + v W , Y sin α z v w , y = v W , X sin α z + v W , Y cos α z
The resulting wind velocity v w is the sum of the apparent wind velocity v x due to the vehicle’s forward motion minus v w , x and the apparent wind velocity v y due to the vehicle’s lateral motion minus v w , y , namely v a w , x = v x v w , x and v a w , y = v y v w , y .
Wind gusts induce longitudinal and lateral forces F x and F y , respectively, as well as roll and yaw moments. The forces and torques due to wind are expressed as
F x = 1 2 A x ρ C d v a w , x 2 F y = 1 2 A y ρ C d v a w , y 2 T ϕ = l c , z F y T θ = l c , z F x T ψ = l c , x F y l c , y F x
where A x and A y represent the vehicle’s front and lateral surfaces, ρ is the air density, C d is aerodynamic coefficient, and l c , x , l c , y are the distances between the center of mass and the center of pressure of each side of the vehicle.
Assumption 4.
For the purposes of this simulation, only the wind components in the x and y directions are considered. This simplification is justified by the following:
 1. 
Lateral wind components (x and y) have a more significant impact on UAV behavior compared to vertical wind components (z).
 2. 
Reducing the model’s complexity facilitates the design and implementation of the controller.
 3. 
Vertical wind components are typically smaller and less perturbative, making this a reasonable and accurate assumption for most UAV operational scenarios.

5.2. Simulation Results of Simulink®

The effectiveness of the suggested control strategy was tested by doing thorough simulations. These simulations were essential for assessing the control system’s performance in different situations and disturbances, guaranteeing its resilience and accuracy in trajectory tracking and stability control.
The performance of the controller has been tested making use of the following references x 1 , r e f = 1 m, y 1 , r e f = 1 m and z 1 , r e f = 1 m.
The system’s initial conditions for position were initialized as follows: starting at x 1 ( 0 ) = 0 m, y 1 ( 0 ) = 0 m, and z 1 ( 0 ) = 0 m. Velocity had an initial state of x 2 ( 0 ) = 0 m/s, y 2 ( 0 ) = 0 m/s, z 2 ( 0 ) = 0 m/s. As for attitude, the system commenced with ϕ 1 ( 0 ) = 0 rad, θ 1 ( 0 ) = 0 rad, and ψ 1 ( 0 ) = 0 rad and the rotational velocities in the body frame were initialized at ϕ 2 ( 0 ) = 0 rad/s, θ 2 ( 0 ) = 0 rad/s, and ψ 2 ( 0 ) = 0 rad/s.
Furthermore, regarding the motors, the initial condition for the current vector was set as i i ( 0 ) = 0 .
To evaluate the robustness of the system, the following parameter variations had been considered: m = 1.05 m , I x x = 1.05 I x x , I y y = 1.05 I y y , I z z = 1.05 I z z , K m , i = 1.05 K m , i , B m , i = 1.05 B m , i , J m , i = 1.05 K m , i and L i = 1.1 L i .
The control gains for the active control are k x 1 = 2 , k y 1 = 2 , k z 1 = 2 , k x 2 = 5 , k y 2 = 5 , k z 2 = 11 , k d , x 2 = 0.5 , k d , y 2 = 0.5 , k ϕ 1 = 35 , k 0 , ϕ 1 = 25 , k θ 1 = 35 , k 0 , θ 1 = 25 , k ψ 1 = 5 , k d , ψ 1 = 0.5 , k ϕ 2 = 60 , k 0 , ϕ 2 = 25 , k θ 2 = 60 , k 0 , θ 2 = 25 , k ψ 2 = 5 , k 0 , ψ 2 = 0.5 , λ x 1 = 9 , λ x 2 = 2 , λ y 1 = 9 , λ y 2 = 2 , λ z 1 = 15 , λ z 2 = 2 , λ ϕ 1 = 10 , λ ϕ 2 = 1 , λ θ 1 = 10 , λ θ 2 = 1 , λ ψ 1 = 5 , λ ψ 2 = 1 . Concerning the motors, the control gains are: k ω , i = 30 and k i , i = 1 , λ ω , 1 , i = 10 , λ ω , 2 , i = 5 , λ i , 1 , i = 10 , λ i , 2 , i = 15 .
Moreover, the wind behavior is given by a nominal velocity of V w , 0 = 20 km/h, a nominal angle of β 0 = 7 π / 4 rad, and a wind oscillation frequency of ω w = 0.1 rad/s. The wind forces exerted on the UAV, along with the torques induced by the wind force, are depicted in Figure 4.
The results are summarized in Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14. Figure 5 shows the quadrotor’s precise positioning illustrated in the three dimensions x 1 , y 1 , and z 1 . The figure marks the initial position of the quadrotor with a green circle at coordinates ( 0 , 0 , 0 ) and the final position with a red circle at coordinates ( 1 , 1 , 1 ) . This visualization highlights the path taken by the UAV from its starting point to its target position, providing a clear representation of its movement through space relative to the specified position references x 1 , r e f , y 1 , r e f , and z 1 , r e f . On the other hand, Figure 6a–c portray the quadrotor’s precise positioning in the three dimensions x 1 , y 1 , and z 1 with respect to their corresponding position references x 1 , r e f , y 1 , r e f , and z 1 , r e f . Simultaneously, Figure 6d–f display the associated tracking errors e x , e y , and e z , providing a visual representation of the outer loop’s performance in the UAV frame. It is evident from these figures that the proposed controller ensures remarkably minimal tracking errors at this outer loop level.
Figure 7 shows the linear velocities x 2 , y 2 , and z 2 of the UAV over time. Specifically, Figure 7a illustrates the linear velocity x 2 , Figure 7b depicts the linear velocity y 2 , and Figure 7c presents the linear velocity z 2 . These subfigures highlight how the UAV’s linear velocities evolve, providing insights into the system’s dynamic response. Additionally, Figure 7d–f display the tracking errors for these linear velocities, e x 2 , e y 2 , and e z 2 , respectively. These errors represent the difference between the desired and actual velocities, showcasing the controller’s effectiveness in minimizing deviations.
Similarly, Figure 8 presents the orientation dynamics of the UAV. Figure 8a–c show the roll angle ϕ 1 , pitch angle θ 1 , and yaw angle ψ 1 , respectively, as functions of time. These plots illustrate the UAV’s rotational behavior. Corresponding tracking errors for these angles, e ϕ 1 , e θ 1 , and e ψ 1 , are depicted in Figure 8d–f. These errors highlight the precision of the attitude control achieved by the proposed system.
Figure 9 shows the angular velocities ϕ 2 , θ 2 , and ψ 2 of the UAV. Specifically, Figure 9a illustrates the angular velocity ϕ 2 in the x-axis, Figure 9b depicts the angular velocity θ 2 in the y-axis, and Figure 9c presents the angular velocity ψ 2 in the z-axis. These subfigures highlight the rotational dynamics and the system’s response. Figure 9d–f display the tracking errors for these angular velocities, e ϕ 2 , e θ 2 , and e ψ 2 , respectively, showcasing the effectiveness of the inner loop control.
Conversely, Figure 10a–d represent the angular velocities of the quadrotor’s four actuators ω 1 , ω 2 , ω 3 , and ω 4 in tandem with their corresponding reference values ω 1 , r e f , ω 2 , r e f , ω 3 , r e f , and ω 4 , r e f . Simultaneously, Figure 10e–h illustrate the respective tracking errors for the inner loop of the UAV actuators. The clear results underscore the effectiveness of the proposed controller in achieving minimal tracking errors within this inner loop.
Transitioning to Figure 11a–d, these figures present the current profiles of the quadrotor’s four actuators i 1 , i 2 , i 3 , and i 4 alongside their reference values i 1 , r e f , i 2 , r e f , i 3 , r e f , and i 4 , r e f . Complementary to this, Figure 11e–h exhibit the corresponding tracking errors.
On the other hand, Figure 12 presents a comprehensive analysis of the estimation process and error dynamics. The subfigures a, b, c, delineate the performance of estimators ξ x 1 , ξ y 1 , and ξ z 1 in tracking the corresponding errors e x , e y , and e z , providing insights into the estimation precision. Additionally, the variations δ x 1 , δ y 1 , and δ z 1 , tracked by estimators ξ x 2 , ξ y 2 , and ξ z 2 , are portrayed in subfigures d, e, f. This visual representation underscores the effectiveness of the estimation process and its ability to minimize errors while accurately capturing state variations.
Finally, Figure 13 and Figure 14 delve into the inner loop dynamics of the UAV’s actuators, shedding light on the estimation and error characteristics. Figure 13a–d illustrates the performance of estimators ξ ω , 1 , 1 , ξ ω , 1 , 2 , ξ ω , 1 , 3 , and ξ ω , 1 , 4 in tracking the corresponding errors e ω , 1 , e ω , 2 , e ω , 3 , and e ω , 4 , while Figure 14a–d showcases the estimators ξ i , 1 , 1 , ξ i , 1 , 2 , ξ i , 1 , 3 , and ξ i , 1 , 4 monitoring the current associated errors e i , 1 , e i , 2 , e i , 3 , and e i , 4 , offering a comprehensive view of the estimation performance. Furthermore, Figure 13e–h and Figure 14e–h depict the variations δ ω , 1 , δ ω , 2 , δ ω , 3 , δ ω , 4 , δ i , 1 , δ i , 2 , δ i , 3 , and δ i , 4 tracked by estimators. The visual representation collectively underscores the effectiveness of the proposed controller in minimizing errors and accurately estimating the angular velocities and currents within the inner loop.

5.3. Simulation Results of PX4 Autopilot Environment

The “PX4 Autopilot Support from UAV Toolbox” is a comprehensive software tool that facilitates the design, simulation, implementation and validation of flight control algorithms for unmanned aerial vehicles (UAVs) using the PX4 autopilot. Its use focused on verifying and evaluating controller performance in various environments, allowing for extensive validation prior to deployment in real flight conditions.
Therefore, for a direct comparison of the controller’s performance, identical flight references, initial conditions, UAV parameters, parametric variations and external disturbances were maintained, consistent with those utilized in the Simulink® simulations.
The control gains in this case are: k x 1 = 0.5 , k y 1 = 0.5 , k z 1 = 315 , k 0 , z 1 = 0.014 , k x 2 = 222 , k 0 , z 2 = 0.2 , k y 2 = 222 , k 0 , z 2 = 0.2 , k z 2 = 138 , k 0 , z 2 = 2.15 , k d , z 2 = 0.45 , k ϕ 1 = 30.6 , k θ 1 = 30.6 , k ψ 1 = 115 .
The results are summarized in Figure 15 where, Figure 15a–c portray the quadrotor’s precise positioning in the three dimensions x 1 , y 1 , and z 1 with respect to their corresponding position references x 1 , r e f , y 1 , r e f , and z 1 , r e f . Simultaneously, Figure 15d–f display the associated tracking errors e x , e y , and e z , providing a visual representation of the outer loop’s performance in the UAV frame. It is evident from these figures that the proposed controller ensures remarkably minimal tracking errors at this outer loop level.

6. Conclusions

Throughout this inquiry, we have designed a complete integrated control system tailored for a UAV equipped with multiple actuators. The main goal was to achieve accurate control over the position reference vector while simultaneously enhancing the controller’s robustness against external disturbances and parameter uncertainties. Numerical simulations were performed using both Simulink® and the simulated PX4 Autopilot environment. The results demonstrated the effectiveness of the proposed control system in achieving precise position control and robust performance for both the UAV and its actuators under various conditions. The control strategy employed two control loops: an outer loop for the UAV frame and an inner loop for the UAV actuators. The outer loop generated the required angular velocities for the actuators to follow the reference position vector, using the UAV’s output. The inner loop ensured that the actuators tracked these angular velocity references. Both control loops utilized PI-like controllers for simplicity. The incorporation of nonlinear control techniques and estimation strategies for disturbances and parameter variations enabled dynamic adaptation to changing environmental conditions. These results underscore the potential applicability of the control system in other UAV operational scenarios.
As part of our future endeavors, we aim to explore the optimization of energy consumption within the system. This includes considering changes to the UAV’s design and architecture, particularly in the context of reconfiguration following the failure of one or more actuators. Investigating the effect of such changes on energy efficiency is consistent with our overarching goal of expanding UAV technologies to not only improve control precision but also address environmental concerns and maximize resource consumption. Additionally, the computational intensity and controller runtime, power consumption, and hardware implementation aspects will be explored in further studies, which are beyond the scope of the current dense article.

Author Contributions

Conceptualization, C.A.L.; methodology, C.A.L., S.D.G. and C.V.V.V.; software, C.V.V.V.; investigation, C.V.V.V.; resources, C.V.V.V.; writing—original draft preparation, C.V.V.V.; writing—review and editing, C.A.L. and C.C.V.G.; supervision, S.D.G.; project administration, C.A.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data that support the findings of this study are available on request from the corresponding author.

Acknowledgments

Vera Vaca C. V. gratefully acknowledges the scholarship number 827938 provided by CONACYT to carry out her postgraduate studies from which this research was derived.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Nex, F.; Armenakis, C.; Cramer, M.; Cucci, D.A.; Gerke, M.; Honkavaara, E.; Kukko, A.; Persello, C.; Skaloud, J. UAV in the advent of the twenties: Where we stand and what is next. ISPRS J. Photogramm. Remote Sens. 2022, 184, 215–242. [Google Scholar] [CrossRef]
  2. Sivakumar, M.; Tyj, N.M. A literature survey of unmanned aerial vehicle usage for civil applications. J. Aerosp. Technol. Manag. 2021, 13, e4021. [Google Scholar] [CrossRef]
  3. Syed, F.; Gupta, S.K.; Hamood Alsamhi, S.; Rashid, M.; Liu, X. A survey on recent optimal techniques for securing unmanned aerial vehicles applications. Trans. Emerg. Telecommun. Technol. 2021, 32, e4133. [Google Scholar] [CrossRef]
  4. Alvares, P.; Silva, L.; Magaia, N. Blockchain-based solutions for UAV-assisted connected vehicle networks in smart cities: A review, open issues, and future perspectives. Telecom 2021, 2, 108–140. [Google Scholar] [CrossRef]
  5. Srivastava, A.; Prakash, J. Techniques, answers, and real-world UAV implementations for precision farming. Wirel. Pers. Commun. 2023, 131, 2715–2746. [Google Scholar] [CrossRef]
  6. Budiyono, A.; Higashino, S.-I. A review of the latest innovations in UAV technology. J. Instrum. Autom. Syst. 2023, 10, 7–16. [Google Scholar] [CrossRef]
  7. Lopez-Sanchez, I.; Moreno-Valenzuela, J. PID control of quadrotor UAVs: A survey. Annu. Rev. Control 2023, 56, 100900. [Google Scholar] [CrossRef]
  8. Rinaldi, M.; Primatesta, S.; Guglieri, G. A comparative study for control of quadrotor UAVs. Appl. Sci. 2023, 13, 3464. [Google Scholar] [CrossRef]
  9. Maaruf, M.; Mahmoud, M.S.; Ma’arif, A. A survey of control methods for quadrotor UAV. Int. J. Robot. Control Syst. 2022, 2, 652–665. [Google Scholar] [CrossRef]
  10. Nguyen, N.P.; Mung, N.X.; Thanh, H.L.N.N.; Huynh, T.T.; Lam, N.T.; Hong, S.K. Adaptive sliding mode control for attitude and altitude system of a quadcopter UAV via neural network. IEEE Access 2021, 9, 40076–40085. [Google Scholar] [CrossRef]
  11. Baek, J.; Kang, M. A synthesized sliding-mode control for attitude trajectory tracking of quadrotor UAV systems. IEEE/ASME Trans. Mechatronics 2023, 28, 2189–2199. [Google Scholar] [CrossRef]
  12. Shao, X.; Sun, G.; Yao, W.; Liu, J.; Wu, L. Adaptive sliding mode control for quadrotor UAVs with input saturation. IEEE/ASME Trans. Mechatronics 2021, 27, 1498–1509. [Google Scholar] [CrossRef]
  13. Noordin, A.; Mohd Basri, M.A.; Mohamed, Z.; Mat Lazim, I. Adaptive PID controller using sliding mode control approaches for quadrotor UAV attitude and position stabilization. Arab. J. Sci. Eng. 2021, 46, 963–981. [Google Scholar] [CrossRef]
  14. 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]
  15. Wang, J.; 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]
  16. Liu, K.; Wang, R. Antisaturation command filtered backstepping control-based disturbance rejection for a quadarotor UAV. IEEE Trans. Circuits Syst. II Express Briefs 2021, 68, 3577–3581. [Google Scholar] [CrossRef]
  17. Bao, C.; Guo, Y.; Luo, L.; Su, G. Design of a fixed-wing UAV controller based on adaptive backstepping sliding mode control method. IEEE Access 2021, 9, 157825–157841. [Google Scholar] [CrossRef]
  18. Xu, L.X.; Ma, H.J.; Guo, D.; Xie, A.H.; Song, D.L. Backstepping sliding-mode and cascade active disturbance rejection control for a quadrotor UAV. IEEE/ASME Trans. Mechatronics 2020, 25, 2743–2753. [Google Scholar] [CrossRef]
  19. Kazim, M.; Azar, A.T.; Koubaa, A.; Zaidi, A. Disturbance-rejection-based optimized robust adaptive controllers for UAVs. IEEE Syst. J. 2021, 15, 3097–3108. [Google Scholar] [CrossRef]
  20. Zhou, L.; Xu, S.; Jin, H.; Jian, H. A hybrid robust adaptive control for a quadrotor UAV via mass observer and robust controller. Adv. Mech. Eng. 2021, 13, 16878140211002723. [Google Scholar] [CrossRef]
  21. Wang, J.; Zhu, B.; Zheng, Z. Robust adaptive control for a quadrotor UAV with uncertain aerodynamic parameters. IEEE Trans. Aerosp. Electron. Syst. 2023, 59, 8313–8326. [Google Scholar] [CrossRef]
  22. Ristevski, S.; Koru, A.T.; Yucelen, T.; Dogan, K.M.; Muse, J.A. Experimental results of a quadrotor UAV with a model reference adaptive controller in the presence of unmodeled dynamic. In Proceedings of the AIAA Scitech 2022 Forum, San Diego, CA, USA, 3–7 January 2022; p. 1381. [Google Scholar] [CrossRef]
  23. García, O.; Ordaz, P.; Santos-Sánchez, O.J.; Salazar, S.; Lozano, R. Backstepping and robust control for a quadrotor in outdoors environments: An experimental approach. IEEE Access 2019, 7, 40636–40648. [Google Scholar] [CrossRef]
  24. Islam, S.; Liu, P.X.; El Saddik, A. Robust control of four-rotor unmanned aerial vehicle with disturbance uncertainty. IEEE Trans. Ind. Electron. 2014, 62, 1563–1571. [Google Scholar] [CrossRef]
  25. Bianchi, D.; Di Gennaro, S.; Di Ferdinando, M.; Acosta Lúa, C. Robust control of UAV with disturbances and uncertainty estimation. Machines 2023, 11, 352. [Google Scholar] [CrossRef]
  26. Xia, K.; Shin, M.; Chung, W.; Kim, M.; Lee, S.; Son, H. Landing a quadrotor UAV on a moving platform with sway motion using robust control. Control Eng. Pract. 2022, 128, 105288. [Google Scholar] [CrossRef]
  27. Din, A.F.U.; Mir, I.; Gul, F.; Al Nasar, M.R.; Abualigah, L. Reinforced learning-based robust control design for unmanned aerial vehicle. Arab. J. Sci. Eng. 2023, 48, 1221–1236. [Google Scholar] [CrossRef]
  28. Zhao, Z.; Cao, D.; Yang, J.; Wang, H. High-order sliding mode observer-based trajectory tracking control for a quadrotor UAV with uncertain dynamics. Nonlinear Dyn. 2020, 102, 2583–2596. [Google Scholar] [CrossRef]
  29. Chandra, A.; Lal, P.P.S. Higher order sliding mode controller for a quadrotor UAV with a suspended load. IFAC-PapersOnLine 2022, 55, 610–615. [Google Scholar] [CrossRef]
  30. COMP4DRONES. Agriculture Use Case—COMP4DRONES Project. Available online: https://www.comp4drones.eu/project-info/use-cases/agriculture/ (accessed on 24 June 2024).
  31. NASA. UAS Traffic Management (UTM) Project. Available online: https://www.nasa.gov/directorates/armd/past-armd-projects/uas-traffic-management-utm-project/ (accessed on 24 June 2024).
  32. Fethalla, N.; Saad, M.; Michalska, H.; Ghommam, J. Robust observer-based dynamic sliding mode controller for a quadrotor UAV. IEEE Access 2018, 6, 45846–45859. [Google Scholar] [CrossRef]
  33. Sanwale, J.; Dahiya, S.; Trivedi, P.; Kothari, M. Robust fault-tolerant adaptive integral dynamic sliding mode control using finite-time disturbance observer for coaxial octorotor UAVs. Control Eng. Pract. 2023, 135, 105495. [Google Scholar] [CrossRef]
  34. Ha, L.N.N.T.; Hong, S.K. Robust dynamic sliding mode control-based PID–super twisting algorithm and disturbance observer for UAVs. Electronics 2019, 8, 760. [Google Scholar] [CrossRef]
  35. Wei, M.; Zheng, L.; Li, H.; Cheng, H. Adaptive Neural Network-based Model Path-Following Contouring Control for Quadrotor Under Diversely Uncertain Disturbances. IEEE Robot. Autom. Lett. 2024. [Google Scholar] [CrossRef]
  36. Madebo, M.M.; Abdissa, C.M.; Lemma, L.N.; Negash, D.S. Robust Tracking Control for Quadrotor UAV with External Disturbances and Uncertainties Using Neural Network Based MRAC. IEEE Access 2024, 12, 36183–36201. [Google Scholar] [CrossRef]
  37. Verberne, J.; Moncayo, H. Robust control architecture for wind rejection in quadrotors. In Proceedings of the International Conference on Unmanned Aircraft Systems (ICUAS), Atlanta, GA, USA, 11–14 June 2019; pp. 152–161. [Google Scholar] [CrossRef]
  38. Wan, K.; Gao, X.; Hu, Z.; Wu, G. Robust motion control for UAV in dynamic uncertain environments using deep reinforcement learning. Remote Sens. 2020, 12, 640. [Google Scholar] [CrossRef]
  39. Wan, K.; Li, B.; Gao, X.; Hu, Z.; Yang, Z. A learning-based flexible autonomous motion control method for UAV in dynamic unknown environments. J. Syst. Eng. Electron. 2021, 32, 1490–1508. [Google Scholar] [CrossRef]
  40. Mayorga-Macías, W.A.; González-Jiménez, L.E.; Meza-Aguilar, M.A.; Luque-Vega, L.F. Velocity Sensor for Real-Time Backstepping Control of a Multirotor Considering Actuator Dynamics. Sensors 2020, 20, 4229. [Google Scholar] [CrossRef]
  41. Wang, B.; Ghamry, K.A.; Zhang, Y. Trajectory tracking and attitude control of an unmanned quadrotor helicopter considering actuator dynamics. In Proceedings of the 2016 35th Chinese Control Conference (CCC), Chengdu, China, 27–29 July 2016; IEEE: New York, NY, USA, 2016; pp. 10795–10800. [Google Scholar] [CrossRef]
  42. Di, J.; Kang, Y.; Ji, H.; Wang, X.; Chen, S.; Liao, F.; Li, K. Low-level control with actuator dynamics for multirotor UAVs. Robot. Intell. Autom. 2023, 43, 290–300. [Google Scholar] [CrossRef]
  43. Ferry, N. Quadcopter Plant Model and Control System Development with MATLAB/Simulink Implementation. Ph.D. Thesis, Rochester Institute of Technology, Rochester, NY, USA, 2017. Available online: http://www.ritravvenlab.com/uploads/1/1/8/4/118484574/ferry.pdf (accessed on 24 June 2024).
  44. Bouabdallah, S. Design and Control of Quadrotors with Application to Autonomous Flying. Ph.D. Thesis, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland, 2007. Available online: https://infoscience.epfl.ch/record/95939?v=pdf (accessed on 24 June 2024).
  45. Mayorga-Macías, W.A. Embedded Control System for a Multi-Rotor Considering Motor Dynamics. Ph.D. Thesis, ITESO—Universidad Jesuita de Guadalajara, San Pedro Tlaquepaque, Mexico, 2020. Available online: https://rei.iteso.mx/items/6cd58f47-b751-4033-ae7f-1d3657d6c258 (accessed on 24 June 2024).
  46. ArduPilot Dev Team. Connect ESCs and Motors. Available online: https://ardupilot.org/copter/docs/connect-escs-and-motors.html (accessed on 6 March 2024).
  47. Cárdenas R, C.A.; Grisales, V.H.; Collazos Morales, C.A.; Cerón-Muñoz, H.D.; Ariza-Colpas, P.; Caputo-Llanos, R. Quadrotor modeling and a pid control approach. In Proceedings of the Intelligent Human Computer Interaction: 11th International Conference, IHCI 2019, Allahabad, India, 12–14 December 2019; Springer: Berlin/Heidelberg, Germany, 2020; pp. 281–291. [Google Scholar] [CrossRef]
  48. Beard, R.W.; McLain, T.W. Small Unmanned Aircraft: Theory and Practice; Princeton University Press: Princeton, NJ, USA, 2024. [Google Scholar]
  49. Choi, J.; Cheon, D.; Lee, J. Robust landing control of a quadcopter on a slanted surface. Int. J. Precis. Eng. Manuf. 2021, 22, 1147–1156. [Google Scholar] [CrossRef]
  50. Guillén-Bonilla, J.T.; Vaca García, C.C.; Di Gennaro, S.; Sánchez Morales, M.E.; Acosta Lúa, C. Vision-Based Nonlinear Control of Quadrotors Using the Photogrammetric Technique. Math. Probl. Eng. 2020, 2020, 5146291. [Google Scholar] [CrossRef]
  51. Acosta Lúa, C.; Vaca García, C.C.; Di Gennaro, S.; Castillo-Toledo, B.; Sánchez Morales, M.E. Real-Time Hovering Control of Unmanned Aerial Vehicles. Math. Probl. Eng. 2020, 2020, 2314356. [Google Scholar] [CrossRef]
  52. Schneider, E.E.; Dale, T. Guidelines for Manned Space Flight Experiments; NASA: Washington, DC, USA, 1971. Available online: https://ntrs.nasa.gov/api/citations/19710016459/downloads/19710016459.pdf (accessed on 6 March 2024).
  53. Nagaty, A.; Saeedi, S.; Thibault, C.; Seto, M.; Li, H. Control and navigation framework for quadrotor helicopters. J. Intell. Robot. Syst. 2013, 70, 1–12. [Google Scholar] [CrossRef]
  54. NASA Glenn Research Center. Drag Equation. Available online: https://www1.grc.nasa.gov/beginners-guide-to-aeronautics/drag-equation/ (accessed on 6 March 2024).
  55. Gonzalez-Sierra, J.; Dzul, A.; Martinez, E. Formation control of distance and orientation based-model of an omnidirectional robot and a quadrotor UAV. Robot. Auton. Syst. 2022, 147, 103921. [Google Scholar] [CrossRef]
  56. Morbidi, F.; Cano, R.; Lara, D. Minimum-energy path generation for a quadrotor UAV. In Proceedings of the 2016 IEEE International Conference on Robotics and Automation (ICRA), Stockholm, Sweden, 16–21 May 2016; IEEE: New York, NY, USA, 2016; pp. 1492–1498. [Google Scholar]
  57. Khalil, H.K. Nonlinear Systems; Prentice Hall: Saddle River, NJ, USA, 2006. [Google Scholar]
  58. MathWorks. UAV Toolbox Support Package for PX4 Autopilots. Available online: https://la.mathworks.com/help/uav/px4-spkg.html?s_tid=CRUX_topnav (accessed on 24 June 2024).
Figure 1. Quadrotor schematic.
Figure 1. Quadrotor schematic.
Drones 08 00447 g001
Figure 2. Brushless DC model.
Figure 2. Brushless DC model.
Drones 08 00447 g002
Figure 3. Robust Control Scheme.
Figure 3. Robust Control Scheme.
Drones 08 00447 g003
Figure 4. (a) External forces applied by the wind. F x component (red), F y component (blue); (b) External torques applied by the force of the wind. T ϕ (red), T θ (blue), and T ψ (green).
Figure 4. (a) External forces applied by the wind. F x component (red), F y component (blue); (b) External torques applied by the force of the wind. T ϕ (red), T θ (blue), and T ψ (green).
Drones 08 00447 g004
Figure 5. Three-dimensional trajectory of the UAV showing the initial point (green circle) at ( 0 , 0 , 0 ) and the final point (red circle) at ( 1 , 1 , 1 ) .
Figure 5. Three-dimensional trajectory of the UAV showing the initial point (green circle) at ( 0 , 0 , 0 ) and the final point (red circle) at ( 1 , 1 , 1 ) .
Drones 08 00447 g005
Figure 6. (a) Position x 1 (black) and reference x 1 , r e f (red dash–dot line); (b) position y 1 (black) and reference y 1 , r e f (red dash–dot line); (c) position z 1 (black) and reference y 1 , r e f (red dash–dot line); (d) tracking error e x = x 1 x 1 , r e f ; (e) tracking error e y = y 1 y 1 , r e f ; (f) tracking error e z = z 1 z 1 , r e f .
Figure 6. (a) Position x 1 (black) and reference x 1 , r e f (red dash–dot line); (b) position y 1 (black) and reference y 1 , r e f (red dash–dot line); (c) position z 1 (black) and reference y 1 , r e f (red dash–dot line); (d) tracking error e x = x 1 x 1 , r e f ; (e) tracking error e y = y 1 y 1 , r e f ; (f) tracking error e z = z 1 z 1 , r e f .
Drones 08 00447 g006
Figure 7. (a) Linear velocity x 2 ; (b) linear velocity y 2 ; (c) linear velocity z 2 ; (d) tracking error e x 2 = x 2 x 2 , d ; (e) tracking error e y 2 = y 2 y 2 , d ; (f) tracking error e z 2 = z 2 z 2 , d .
Figure 7. (a) Linear velocity x 2 ; (b) linear velocity y 2 ; (c) linear velocity z 2 ; (d) tracking error e x 2 = x 2 x 2 , d ; (e) tracking error e y 2 = y 2 y 2 , d ; (f) tracking error e z 2 = z 2 z 2 , d .
Drones 08 00447 g007
Figure 8. (a) Roll ϕ 1 ; (b) pitch θ 1 ; (c) yaw ψ 1 ; (d) tracking error e ϕ 1 = sin ϕ 1 sin ϕ 1 , d ; (e) tracking error e θ 1 = sin θ 1 sin θ 1 , d ; (f) tracking error e ψ 1 = ψ 1 ψ 1 , r e f .
Figure 8. (a) Roll ϕ 1 ; (b) pitch θ 1 ; (c) yaw ψ 1 ; (d) tracking error e ϕ 1 = sin ϕ 1 sin ϕ 1 , d ; (e) tracking error e θ 1 = sin θ 1 sin θ 1 , d ; (f) tracking error e ψ 1 = ψ 1 ψ 1 , r e f .
Drones 08 00447 g008
Figure 9. (a) Angular velocity in x, ϕ 2 ; (b) angular velocity in y, θ 2 ; (c) angular velocity in z, ψ 2 ; (d) tracking error e ϕ 2 = ϕ 2 ϕ 2 , d ; (e) tracking error e θ 2 = θ 2 θ 2 , d ; (f) tracking error e ψ 2 = ψ 2 ψ 2 , d .
Figure 9. (a) Angular velocity in x, ϕ 2 ; (b) angular velocity in y, θ 2 ; (c) angular velocity in z, ψ 2 ; (d) tracking error e ϕ 2 = ϕ 2 ϕ 2 , d ; (e) tracking error e θ 2 = θ 2 θ 2 , d ; (f) tracking error e ψ 2 = ψ 2 ψ 2 , d .
Drones 08 00447 g009
Figure 10. (a) Angular velocity ω 1 (black) and reference ω 1 , r e f (red dash–dot line); (b) angular velocity ω 2 (black) and reference ω 2 , r e f (red dash–dot line); (c) angular velocity ω 3 (black) and reference ω 3 , r e f (red dash–dot line); (d) angular velocity ω 4 (black) and reference ω 4 , r e f (red dash–dot line); (e) tracking error e ω 1 = ω 1 ω 1 , r e f ; (f) tracking error e ω 2 = ω 2 ω 2 , r e f ; (g) tracking error e ω 3 = ω 3 ω 3 , r e f ; (h) tracking error e ω 4 = ω 4 ω 4 , r e f .
Figure 10. (a) Angular velocity ω 1 (black) and reference ω 1 , r e f (red dash–dot line); (b) angular velocity ω 2 (black) and reference ω 2 , r e f (red dash–dot line); (c) angular velocity ω 3 (black) and reference ω 3 , r e f (red dash–dot line); (d) angular velocity ω 4 (black) and reference ω 4 , r e f (red dash–dot line); (e) tracking error e ω 1 = ω 1 ω 1 , r e f ; (f) tracking error e ω 2 = ω 2 ω 2 , r e f ; (g) tracking error e ω 3 = ω 3 ω 3 , r e f ; (h) tracking error e ω 4 = ω 4 ω 4 , r e f .
Drones 08 00447 g010
Figure 11. (a) Current i 1 (black) and reference i 1 , r e f (red dash–dot line); (b) current i 2 (black) and reference i 2 , r e f (red dash–dot line); (c) current i 3 (black) and reference i 3 , r e f (red dash–dot line); (d) current i 4 (black) and reference i 4 , r e f (red dash–dot line); (e) tracking error e i 1 = i 1 i 1 , r e f ; (f) tracking error e i 2 = i 2 i 2 , r e f ; (g) tracking error e i 3 = i 3 i 3 , r e f ; (h) tracking error e i 4 = i 4 i 4 , r e f .
Figure 11. (a) Current i 1 (black) and reference i 1 , r e f (red dash–dot line); (b) current i 2 (black) and reference i 2 , r e f (red dash–dot line); (c) current i 3 (black) and reference i 3 , r e f (red dash–dot line); (d) current i 4 (black) and reference i 4 , r e f (red dash–dot line); (e) tracking error e i 1 = i 1 i 1 , r e f ; (f) tracking error e i 2 = i 2 i 2 , r e f ; (g) tracking error e i 3 = i 3 i 3 , r e f ; (h) tracking error e i 4 = i 4 i 4 , r e f .
Drones 08 00447 g011
Figure 12. (a) Estimator ξ x 1 (black) and error e x (red dash–dot line); (b) estimator ξ y 1 (black) and error e y (red dash–dot line); (c) estimator ξ z 1 (black) and error e z (red dash–dot line); (d) estimator ξ x 2 (black) and variation δ x 1 (red dash–dot line); (e) estimator ξ y 2 (black) and variation δ y 1 (red dash–dot line); (f) estimator ξ z 2 (black) and variation δ z 1 (red dash–dot line).
Figure 12. (a) Estimator ξ x 1 (black) and error e x (red dash–dot line); (b) estimator ξ y 1 (black) and error e y (red dash–dot line); (c) estimator ξ z 1 (black) and error e z (red dash–dot line); (d) estimator ξ x 2 (black) and variation δ x 1 (red dash–dot line); (e) estimator ξ y 2 (black) and variation δ y 1 (red dash–dot line); (f) estimator ξ z 2 (black) and variation δ z 1 (red dash–dot line).
Drones 08 00447 g012
Figure 13. (a) Estimator ξ ω , 1 , 1 (black) and error e ω , 1 (red dash–dot line); (b) estimator ξ ω , 1 , 2 (black) and error e ω , 2 (red dash–dot line); (c) estimator ξ ω , 1 , 3 (black) and error e ω , 3 (red dash–dot line); (d) estimator ξ ω , 1 , 4 (black) and error e ω , 4 (red dash–dot line); (e) estimator ξ ω , 2 , 1 (black) and variation δ ω , 1 (red dash–dot line); (f) estimator ξ ω , 2 , 2 (black) and variation δ ω , 2 (red dash–dot line); (g) estimator ξ ω , 2 , 3 (black) and variation δ ω , 3 (red dash–dot line); (h) estimator ξ ω , 2 , 4 (black) and variation δ ω , 4 (red dash–dot line).
Figure 13. (a) Estimator ξ ω , 1 , 1 (black) and error e ω , 1 (red dash–dot line); (b) estimator ξ ω , 1 , 2 (black) and error e ω , 2 (red dash–dot line); (c) estimator ξ ω , 1 , 3 (black) and error e ω , 3 (red dash–dot line); (d) estimator ξ ω , 1 , 4 (black) and error e ω , 4 (red dash–dot line); (e) estimator ξ ω , 2 , 1 (black) and variation δ ω , 1 (red dash–dot line); (f) estimator ξ ω , 2 , 2 (black) and variation δ ω , 2 (red dash–dot line); (g) estimator ξ ω , 2 , 3 (black) and variation δ ω , 3 (red dash–dot line); (h) estimator ξ ω , 2 , 4 (black) and variation δ ω , 4 (red dash–dot line).
Drones 08 00447 g013
Figure 14. (a) Estimator ξ i , 1 , 1 (black) and error e i , 1 (red dash–dot line); (b) estimator ξ i , 1 , 2 (black) and error e i , 2 (red dash–dot line); (c) estimator ξ i , 1 , 3 (black) and error e i , 3 (red dash–dot line); (d) estimator ξ i , 1 , 4 (black) and error e ω , 4 (red dash–dot line); (e) estimator ξ i , 2 , 1 (black) and variation δ i , 1 (red dash–dot line); (f) estimator ξ i , 2 , 2 (black) and variation δ i , 2 (red dash–dot line); (g) estimator ξ i , 2 , 3 (black) and variation δ i , 3 (red dash–dot line); (h) estimator ξ i , 2 , 4 (black) and variation δ i , 4 (red dash–dot line).
Figure 14. (a) Estimator ξ i , 1 , 1 (black) and error e i , 1 (red dash–dot line); (b) estimator ξ i , 1 , 2 (black) and error e i , 2 (red dash–dot line); (c) estimator ξ i , 1 , 3 (black) and error e i , 3 (red dash–dot line); (d) estimator ξ i , 1 , 4 (black) and error e ω , 4 (red dash–dot line); (e) estimator ξ i , 2 , 1 (black) and variation δ i , 1 (red dash–dot line); (f) estimator ξ i , 2 , 2 (black) and variation δ i , 2 (red dash–dot line); (g) estimator ξ i , 2 , 3 (black) and variation δ i , 3 (red dash–dot line); (h) estimator ξ i , 2 , 4 (black) and variation δ i , 4 (red dash–dot line).
Drones 08 00447 g014
Figure 15. (a) Position x 1 (black) and reference x 1 , r e f (red dash–dot line); (b) position y 1 (black) and reference y 1 , r e f (red dash–dot line); (c) position z 1 (black) and reference y 1 , r e f (red dash–dot line); (d) tracking error e x = x 1 x 1 , r e f ; (e) tracking error e y = y 1 y 1 , r e f ; (f) tracking error e z = z 1 z 1 , r e f .
Figure 15. (a) Position x 1 (black) and reference x 1 , r e f (red dash–dot line); (b) position y 1 (black) and reference y 1 , r e f (red dash–dot line); (c) position z 1 (black) and reference y 1 , r e f (red dash–dot line); (d) tracking error e x = x 1 x 1 , r e f ; (e) tracking error e y = y 1 y 1 , r e f ; (f) tracking error e z = z 1 z 1 , r e f .
Drones 08 00447 g015
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

Vera Vaca, C.V.; Di Gennaro, S.; Vaca García, C.C.; Acosta Lúa, C. Robust Nonlinear Control with Estimation of Disturbances and Parameter Uncertainties for UAVs and Integrated Brushless DC Motors. Drones 2024, 8, 447. https://doi.org/10.3390/drones8090447

AMA Style

Vera Vaca CV, Di Gennaro S, Vaca García CC, Acosta Lúa C. Robust Nonlinear Control with Estimation of Disturbances and Parameter Uncertainties for UAVs and Integrated Brushless DC Motors. Drones. 2024; 8(9):447. https://doi.org/10.3390/drones8090447

Chicago/Turabian Style

Vera Vaca, Claudia Verónica, Stefano Di Gennaro, Claudia Carolina Vaca García, and Cuauhtémoc Acosta Lúa. 2024. "Robust Nonlinear Control with Estimation of Disturbances and Parameter Uncertainties for UAVs and Integrated Brushless DC Motors" Drones 8, no. 9: 447. https://doi.org/10.3390/drones8090447

APA Style

Vera Vaca, C. V., Di Gennaro, S., Vaca García, C. C., & Acosta Lúa, C. (2024). Robust Nonlinear Control with Estimation of Disturbances and Parameter Uncertainties for UAVs and Integrated Brushless DC Motors. Drones, 8(9), 447. https://doi.org/10.3390/drones8090447

Article Metrics

Back to TopTop