Next Article in Journal
BNS: A Detection System to Find Nodes in the Bitcoin Network
Next Article in Special Issue
Coalitional Control Strategy for a Heterogeneous Platoon Application
Previous Article in Journal
RETRACTED: Joshi et al. An Evolutionary Technique for Building Neural Network Models for Predicting Metal Prices. Mathematics 2023, 11, 1675
Previous Article in Special Issue
Higher Order Sliding Mode Control of MIMO Induction Motors: A New Adaptive Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predefined-Time Fractional-Order Tracking Control for UAVs with Perturbation

1
MIET Laboratory, Faculty of Science and Technology, Hassan First University of Settat, Settat 26000, Morocco
2
LIS UMR CNRS 7020, Aix-Marseille University, 13013 Marseille, France
3
Department of Computer & Network Engineering, College of Computer Science and Engineering, University of Jeddah, Jeddah 21959, Saudi Arabia
4
Department of Cybersecurity, College of Computer Science and Engineering, University of Jeddah, Jeddah 23218, Saudi Arabia
*
Authors to whom correspondence should be addressed.
Mathematics 2023, 11(24), 4886; https://doi.org/10.3390/math11244886
Submission received: 8 November 2023 / Revised: 2 December 2023 / Accepted: 4 December 2023 / Published: 6 December 2023
(This article belongs to the Special Issue Control Theory and Applications, 2nd Edition)

Abstract

:
This manuscript describes the design of a controller that assures predefined-time convergence in fractional-order sliding mode control (PTFOSMC) for a quadrotor UAV subjected to matched perturbation. Moreover, predefined-time techniques enable the establishment of a time constraint for convergence as a control parameter, distinguishing them from finite- and fixed-time controllers. The proposed control offers the advantage of sliding mode control, exhibiting rapid response and robust performance for the quadrotor subsystems. Notably, the suggested controller is devoid of terms dependent on the initial conditions of the quadrotor. Additionally, an established switching-type predefined-time controller with fractional-order is introduced to bolster robustness against external disturbances and alleviate the chattering problem associated with the sliding mode technique. The application of the Lyapunov function is employed to analyze the predefined-time stability of the quadrotor utilizing the suggested PTFOSMC. Numerical results are provided to demonstrate the effectiveness of the suggested scheme.

1. Introduction

1.1. Motivations and Background

Quadrotors have grown in popularity in recent years as a result of rapid technical advancement, which has had a considerable impact on current electronic components, such as microcontrollers, sensors, and low-cost development. These advancements work in favor of this platform, propelling it to mainstream status. Many researchers have tackled the topic of trajectory tracking of quadrotor systems in the past decade due to the numerous applications of this vehicle [1,2], including military, highway monitoring [3], photography and videography, geographic mapping of inaccessible terrain and locations, delivery systems [4], surveillance [5], detection, scientific research and development, identification of different objects, search and rescue missions [6], and agriculture. The dynamics of the quadrotor system in real flight involve parameter uncertainty and external disturbances, which negatively impact stabilization and tracking performance. Different control methods have been developed to address these challenges, particularly sliding mode control (SMC) and advanced control systems based on SMC, such as nonsingular fast terminal SMC [7], fixed-time distributed robust SMC [8], adaptive backstepping fast terminal SMC [9], developing PID in sliding mode controllers [10], etc.

1.2. Related Works

One of the most challenging aspects of UAVs is achieving smooth, precise flight. To attain this, optimal, stable controllers must be implemented under external disturbances. A literature review will be offered in this subsection. An investigation of sliding mode controllers in predefined time is suggested by the authors in [11]. They conducted a comparison analysis that emphasized the benefits and drawbacks of each technique used for the flight controller to drive a quadrotor. According to the findings of the comparison analysis, the fractional-order integral terminal sliding-mode control method was chosen as the optimal solution for a quadrotor system. The quadrotor model was obtained using the Lagrange-Euler approach in [12]. To adhere to the position trajectory, a predictive control approach with integral action is formulated, and H control is devised to stabilize the orientation of the quadrotor. The researchers in [13] offer a robust control strategy based on a second-order SMC approach to improve quadrotor position tracking performance. In [14], the external disturbances and unknown states of the quadrotor are calculated using an observer. This is done to estimate the total disturbances and control the quadrotor. The controller and observer are designed in the attitude and position subsystems. An observer is formed to assess outside disturbances and unknown situations. The TSMC is intended to provide finite-time convergence performance through composite learning, neural approximation, and disturbance estimation [15].
The current study provides a FO hybrid finite-time control strategy with fractional-order dynamics to address the trajectory tracking of quadrotors, motivated by previous research and inspired by Reference [16]. The quadrotor system, benefiting from fixed-time stability, can achieve convergence to a stable range within a specified upper-bound convergence time, independent of the initial operational states [17]. As a result, the authors in ref. [18] introduces distributed fixed-time control methods to enhance the performance of power systems. On the flip side, the fixed-time control approach exhibits superior capabilities in surmounting these limitations. Nevertheless, estimating the convergence time of the nonlinear system under the fixed-time control approach is not straightforward; it requires the use of a sophisticated estimation function grounded in tuning parameters [19]. Predefined-time stability can keep the nonlinear system stable within a predetermined upper limited convergence time, which is directly associated with an adjustable variable [20,21]. Although less conservative than the fixed-time control approaches now in use, the system’s upper limited convergence time under predefined-time control methods is still adequate [22,23]. The fixed-time control method and the finite-time control method have drawbacks, which can be more effectively addressed by the predefined-time control strategy. Fractional-order integral/derivative operators have been extensively used in research on the control of engineering systems. A fractional-order error manifold for robot manipulators was developed in [24] to achieve good performances and a more flexible PID controller structure.
In general, sliding-mode control stands out as a potent nonlinear control strategy employing discontinuous control to steer system state trajectories onto a predefined sliding surface. Its straightforward implementation and advantageous resilience characteristics make the SMC approach a popular choice, particularly for its robustness against model parameter uncertainty and external disturbances. Indeed, the effectiveness of this approach, a quality relatively rare among nonlinear design strategies, has spurred the exploration of various new research avenues within the broader field of sliding mode control. The majority of current SMC methods utilize a linear sliding mode surface, leading the system state variables to converge asymptotically or exponentially to the origin within a finite time [14]. To address this problem, terminal SMC was introduced to ensure finite-time stability [7,11]. A hybrid finite-time control technique based on nonsingular terminal SMC and adaptive integral SMC approaches is provided to handle quadrotor path tracking problems with unknown dynamics and disturbances in [25]. The initial conditions of the state variables adversely affect the achieved finite time, with this impact intensifying as the values of these initial states increase. Recent endeavors in this domain have focused on fixed-time stability [17,19]. On the flip side, the settling time holds significance in the trajectory tracking of the quadrotor system. Consequently, a novel concept called predefined-time stability has emerged to enhance state convergence [19,20,22]. Moreover, the aim is to showcase the efficacy and superior performance of the proposed control for the quadrotor system in comparison to existing control methods. Subsequently, communication network and load changes are integrated into simulations using the Matlab/Simulink platform.

1.3. Contributions

Inspired by the works described in the related works subsection and guided by the paper in [26], a novel predefined-time fractional-order sliding mode controller (PTFOSMC) method based on fractional calculus, SMC, and predefined-time stability is suggested to improve the dynamic performance of the quadrotor. The current study examines the disturbances and variations in the drag coefficients of rotational and translational subsystems. The proposed PTFOSMC ensures that the quadrotor system will converge to a stable zone with specified performance within the predetermined time. In comparison with standard SMC [9,10], the PTFOSMC may offer better results, such as faster convergence speed, reduced overshoot, lower error, and lower chattering impact. The contributions of this research can be encapsulated in the following key points:
  • The proposed method completes the stabilization and optimization procedures within a predefined time.
  • The predefined-time fractional-order sliding mode control exhibits minimal overshoot, the smallest error, and faster variable convergence.
  • Utilizing the proposed control, the stability of the quadrotor system within a predefined time is assessed through the application of the Lyapunov function. Simulations are employed to validate the efficacy of the suggested control for the quadrotor system.
  • The simulation results show that, in comparison to current flight controllers, the suggested controller enables the quadrotor to achieve the intended flight trajectory with high accuracy.

2. Preliminary Knowledge

2.1. Predefined-Time Stability

The subsequent section revisits the predefined-time stability of integer-order systems. We examine the characteristics of the integer-order system.
x ˙ = g 1 ( t , x )
where x R n represents the system state. Function g 1 : R + × R n R n is nonlinear, and if we consider the origin as an equilibrium point, then g 1 ( t , 0 ) = 0 . The initial condition is x 0 = x ( 0 ) . The subsequent definitions bear particular significance.
Definition 1 
(Globally finite-time design [27]). Any solution x ( t ) of (1) converges to the origin at some finite time, and the origin of (1) is asymptotically stable; that is, x ( t ) = 0 t T ( x ( 0 ) ) , where T : R n R 0 is the function of the settling time.
Definition 2 
(Globally fixed-time design [28]). If the settling-time function is bounded and the origin of system (1) is finite-time stable, then the basis is fixed-time stable, i.e., T max > 0 : T ( x ( 0 ) ) T max , x ( 0 ) R n .
Definition 3 
(Settling-time design [29]). Let the set of all settling-time function bounds for system (1) be written above:
T = T max R + : T x 0 T max .
Furthermore, the ensuing definitions elucidate the design centered on prescribed-time stability and the design rooted in prescribed-time stability.
Definition 4 
(Prescribed-time design [29]). Take the T set given in (2). In the specific instance for system (1), T m a x can be modified by a specific selection of system parameters, ρ , T max = T max ( ρ ) , referred to the concept of prescribed-time stability. This design is accomplished by selecting T max ( ρ ) T and finding the inverse of the settling-time function, allowing for modifying ρ.
Remark 1. 
It is crucial to emphasize that the precise fixed stabilization time for a system built upon prescribed-time stability remains unknown but is confined within T m a x ( ρ ) . In contrast, a design system employing predefined-time stability features a known stabilization time.
Definition 5 
(Predefined-time design [30]). For a predefined-time constant, T c > 0 , the origin of system (1) is predefined-time stable if it is fixed-time stable and the settling-time function is fulfilled.
T ( x ( 0 ) ) T c , x ( 0 ) R n ,
where  T c  is signified by a predefined time.

2.2. Fractional-Order Calculus

The differentiation and integration of fractional order are studied in fractional calculus. The following operators are defined [31,32].
Definition 6. 
The following definition applies to the function f ( t ) using the Riemann–Liouville fractional derivative order α [31]:
a R L D t α f ( t ) = 1 Γ ( m α ) d m d t m a t f ( τ ) ( t τ ) α m + 1 d τ
where  Γ ( . )  is the Gamma function and  ( m 1 ) < α m .
Definition 7. 
The definition of the fractional-order derivative, as introduced by Caputo for a function f ( t ) , can be expressed as follows [31]:
a C D t α f ( t ) = 1 Γ ( m α ) a t f m ( τ ) ( t τ ) α m + 1 d τ , ( m 1 ) < α < m
with  m N *  and  α R + .
Property 1. 
The Caputo derivative is the equivalent of the following equality [32].
t 0 C D t α t 0 C D t β f ( t ) = t 0 C D t α β f ( t )
where  α > β > 0
Property 2. 
If   1 > α > 0 , we have the ability to write 4 the Caputo derivative [32]
t 0 C D t 1 q t 0 C D t α f ( t ) = t 0 C D t α t 0 C D t 1 q f ( t ) = f ˙ ( t )
To represent the Caputo operator, C D q will be used in place of C D α throughout this paper.
Lemma 1 
([28]). Suppose the following condition can be satisfied:
V ˙ μ V v ς V τ , V ( 0 ) = V 0
where V denotes a positive-definite function; 0 < τ < 1 , v > 1 , μ , and ς are positive constants. The system x ˙ = f ( x , t ) will exhibit fixed-time stability, and the upper limit on the convergence time can be attained through
T T max = 1 μ ( v 1 ) + 1 ζ ( 1 τ )
Lemma 2 
([22]). If V can meet the following predefined-time stabilizing function:
V ˙ 1 μ T s e V μ V 1 μ , 0 < μ < 1
then x ˙ = f ( x , t ) will achieve predefined-time stability within a predefined time T s > 0 . That is, x ˙ = f ( x , t ) will converge to zero for  t T s .
Remark 2. 
As demonstrated in (5), the upper limited convergence time cannot be achieved explicitly because the adjustable parameters affect how the upper limited convergence time of the fixed-time stability is evaluated. The upper bound convergence time of the predefined-time stability is equal to the tuning parameter, as demonstrated in Lemma 2, eliminating the drawbacks of the fixed-time stability.

3. Problem Formulation

The quadrotor employed in this investigation is illustrated in Figure 1. This apparatus features a sturdy body with four rotors. Yaw movement is generated by creating a speed differential between two rotors moving in opposite directions, specifically between rotors (2,4) and (1,3). Vertical displacement of the quadrotor is accomplished by adjusting the total rotor speeds. Forward motion is achieved by modifying the speeds of the (1 and 3) propellers. The lateral motion is created by varying the speed of the propellers (2 and 4). As a result, the body is intrinsically unstable, with a mechanical system that is tightly linked and under-actuated. Furthermore, the quadrotor dynamic model is developed on the following assumptions; the influence of the ground is disregarded, the quadrotor structure is symmetrical, the quadrotor frame and its blades possess rigidity, and the torques and thrust produced by the rotor speeds exhibit a proportionality to the square of the rotor rotation speeds.
The quadrotor provides six outputs [ ϕ , θ , ψ , x , y , z ] T in its translational and rotational motions. Four actual voltage inputs [ ( Π 1 , Π 2 , Π 3 , Π 4 ] T drive these outputs. The earth-fixed frame and the body frames are represented by E = [ x e , y e , z e ] T and B = [ x b , y b , z b ] T , respectively. The quadrotor dynamics can be expressed using the Newton–Euler formalism as follows [7,33].
ϕ ¨ = δ 1 θ ˙ ψ ˙ + δ 2 ω r θ ˙ 2 + δ 3 ϕ ˙ + J x 1 Π 2 + d ϕ θ ¨ = δ 4 ϕ ˙ ψ ˙ + δ 5 ω r ϕ ˙ 2 + δ 6 θ ˙ + J y 1 Π 3 + d θ ψ ¨ = δ 7 ϕ ˙ θ ˙ + δ 8 ψ ˙ + J z 1 Π 4 + d ψ x ¨ = δ 9 x ˙ + Π x + d x y ¨ = δ 10 y ˙ + Π y + d y z ¨ = δ 11 z ˙ + Π z + d z
with
δ 1 = ( J y J z ) / J x , δ 2 = J r / J x , δ 3 = K a x / J x
δ 4 = ( J z J x ) / J y , δ 5 = J r / J y , δ 6 = K a y / J y
δ 7 = ( J x J y ) / J z , δ 8 = K a z / J z , δ 9 = k x / m .
δ 10 = k y / m , δ 11 = k z / m
where m is the total body mass and g denotes gravitational acceleration. Inertia moments of the quadrotor around the x, y, and z axes are denoted by J x , J y , and J z . The torques and thrusts are represented by Π 1 4 , and d x ψ represent the disturbances. The following equation represents the relation between the control signals and the rotor speeds:
Π 1 Π 2 Π 3 Π 4 = 0 0 0 0 ω 1 2 ω 2 2 ω 3 2 ω 4 2
where represents the drag coefficient, ω 1 4 denotes the angular rotor speeds, ♭ and are positive constants. The vehicle must be stabilized and follow the reference trajectory [ ϕ d , θ d , ψ d , X d , Y d , Z d ] T within a finite time. The aim of this study is to formulate a reliable controller capable of supplying thrust magnitude and torque. To streamline the approach to control system design, the following definition relates to the virtual control signals utilized to determine the overall thrust and desired angles.
Π x = ( cos ( ϕ ) sin ( θ ) cos ( ψ ) + sin ( ϕ ) sin ( ψ ) ) Π 1 m Π y = ( cos ( ϕ ) sin ( θ ) sin ( ψ ) sin ( ϕ ) cos ( θ ) ) Π 1 m Π z = cos ( ϕ ) cos ( θ ) Π 1 m g
Based on the position controls, we can calculate the thrust Π 1 and the two attitude angles, ϕ d and θ d :
ϕ d = a t a n [ c ( θ d ) ( Π x s ( ψ d ) Π y c ( ψ d ) ) / ( Π z + g ) ] θ d = a t a n [ ( Π x cos ( ψ d ) + Π y sin ( ψ d ) ) / ( Π z + g ) ] Π 1 = m Π x 2 + Π y 2 + ( Π z + g ) 2

4. Predefined-Time Fractional-Order Controller Design

This section focuses on the design of the flight controller for the quadrotor system. Under the control of the flight controller, the closed-loop stability in the predetermined time is guaranteed, and the trajectories follow their references. The outer and inner loops now use a new predefined-time fractional-order sliding mode control, as seen in Figure 2. In comparison to the backstepping sliding mode and developing PID sliding mode controls, the proposed control enhances the tracking capabilities of the path reference and increases the resilience control of the quadrotor against external disturbances.

4.1. Outer-Loop Control Design

To guarantee the predefined-time stability of the position subsystem and the generation of tilting angles, the main challenge tackled in this section involves the formulation of virtual control laws.
For the elevation subsystem, where the PTFOSMC is employed, the definition of the tracking error for the z ( t ) subsystem and its time derivative can be established as follows:
e x = z z d e ˙ x = z ˙ z ˙ d
Remark 3. 
The predefined-time controller ensures the convergence of the tracking error and its time derivative using the parameter  T f [34], such that
lim t T f e z ( t ) = 0 , lim t T f e ˙ z ( t ) = 0 , a n d e z ( t ) = e ˙ z ( t ) = 0 t > T f
The following predefined-time fractional-order sliding manifold is suggested to stabilize the altitude subsystem.
σ z = I α z { 1 T s z μ z e ( | e z | μ z ) e z 1 μ z } + γ z D 1 α z e z ( t )
where e z is the error tracking, T s z denotes the predefined setting upper bound convergence time, and e z ( t ) μ z 2 = | e z | . sign ( e z ( t ) ) with 0 < μ z 2 < 1 .
When the consensus error e z reaches the fractional-order sliding surface and finds the equivalent control law, the following equations can be achieved:
D α σ z = 0
Assuming that, using (12) and (13), (15) can be written as
1 T s z μ z e ( | e z | μ z ) e z 1 μ z + γ z D 1 e z ( t ) = 0
As a result, it is possible to generate the dynamics of the fractional-order sliding surface in (12), as
e ˙ z ( t ) = 1 γ z T s z μ z e ( | e z | μ z ) e z 1 μ z
The fractional-order sliding surface σ z will converge to zero if the dynamic system in (15) is stable.
Theorem 1. 
Consider the PTFO sliding manifold designed in (12); the convergence to zero in a predefined time T s z of the dynamics system defined in (7) is ensured.
Proof. 
The function is chosen to confirm the system’s stability. Function (16) is chosen to ensure the stability of the z-subsystem:
V 1 z = e z
The time differentiation of V 1 z is achieved as follows:
V ˙ 1 z = sign e z e ˙ z = sign e z D α D 1 α e z = sign e z D α D α z 1 γ z T s z μ z e ( | e z | μ z ) e z 1 μ z = sign e z 1 γ z T s z μ z e ( | e z | μ z ) e z 1 μ z
Due to sign e z × sign e z = 1 , (18) can be achieved as follows:
V ˙ z 1 = sign e i 1 γ z T s z μ z e ( | e z | μ z ) sign e i | e z | 1 μ z = 1 γ z T s z μ z e ( | e z | μ z ) | e z | 1 μ z = 1 γ z T s z μ z e ( | V 1 z | μ z ) | V 1 z | 1 μ z 0
Based on Lemma 2, it is shown that the dynamic system described by (18) can achieve convergence to zero within a predefined time T s z . In order to prove this convergence, the following equation can be rewritten:
V ˙ z 1 = d V z 1 d t = 1 γ z T s z μ z e ( | V 1 z | μ z ) | V 1 z | 1 μ z
Using some calculations, we obtain
d t = γ z T s z μ z e ( | V 1 z | μ z ) | V 1 z | 1 μ z d V 1 z d t = γ z T s z μ z e ( | V 1 z | μ z ) | V 1 z | μ z 1 d V 1 z
By substituting Equation (16), it can be written as  V 1 z = | e z | > 0
d t = γ z T s z μ z V 1 z μ z 1 e V 1 z μ z d V z 1
Substitute u = V 1 z μ z d u = μ z V 1 z μ z 1 d V 1 z
d t = γ z T s z μ z 1 μ z e u d u
By integrating Equation (22) from t 0 to t r z , one can obtain
t r z = t 0 + γ z T s z e u
Undo substitution u = V 1 z μ z .
t r z = t 0 + γ z T s z e V 1 z μ z
The proof from above is now complete.
The following equation can be obtained from properties (3) and (15),
e ¨ z ( t ) = D 1 + α z D 1 α z e z ( t ) = D 1 + α z D α z { 1 γ z T s z μ z e ( | e z | μ z ) e z 1 μ z } = D 1 { 1 γ z T s z μ z e ( | e z | μ z ) e z 1 μ z }
Therefore, the equivalent controller Π z can be obtained as follows:
z ¨ z ¨ d = D 1 { 1 γ z T s z μ z e ( | e z | μ z ) e z 1 μ z }
Π z e q = D 1 { 1 γ z T s z μ z e ( | e z | μ z ) e z 1 μ z } δ 11 z ˙ d z + z ¨ d
After that, system (7) is designed to reach the PTFO sliding manifolds by the switching controller u s .
Π z s = D α z { 1 T c z μ z 2 e ( | σ z | μ z 2 ) σ z 1 μ z 2 + K z sign ( σ z ) }
where T c z > 0 and μ z 2 , K z are positive parameters. Combining (27) and (28), the following law determines the altitude dynamics of the quadrotor:
Π z = Π z e q + Π z s = D 1 { 1 γ z T s z μ z e ( | e z | μ z ) e z 1 μ z } D α z { 1 T c z μ z 2 e ( | σ z | μ z 2 ) σ z 1 μ z 2 + K z sign ( σ z ) } δ 11 z ˙ d z + z ¨ d
Remark 4. 
The proposed controller is divided into two components. The initial segment incorporates two layers of sliding mode variables and one fractional-order operator. The second component of the proposed controller relies on a fractional-order switching controller. Due to the integration of fractional-order operators in its design, the suggested controller exhibits increased robustness to variations in disturbance frequencies.

4.2. Altitude Stability Analysis for a Quadrotor

Theorem 2. 
During a predefined time T s z , the proposed control (29) can regulate the consensus error e 1 z to reach the FOSS.
Proof. 
Generate the Lyapunov function and the first-order time derivative as follows:
V 2 z = σ z V ˙ 2 z = sign σ z σ ˙ z = sign σ z D α D 1 + α σ z
V ˙ 2 z = sign σ z D α e ¨ z ( t ) + D 1 { 1 T s z μ z e ( | e z | μ z ) e z 1 μ z } )
Based on (7) and the PTFOSMC (29), (32) can be yielded as follows:
V ˙ 2 z = sign σ z σ ˙ z = sign σ z D α D 1 + α σ i = sign σ z D α [ δ 11 z ˙ + Π z + d z z ¨ d + D 1 { 1 T s z μ z e ( | e z | μ z ) e z 1 μ z } ) ]
V ˙ 2 z = sign σ z D α [ δ 11 z ˙ D 1 { 1 T s z μ z e ( | e z | μ z ) e z 1 μ z } D α z { 1 T c z μ z 2 e ( | σ z | μ z 2 ) σ z 1 μ z 2 + K z sign ( σ z ) } δ 11 z ˙ d z + z ¨ d + d z z ¨ d + D 1 { 1 T s z μ z e ( | e z | μ z ) e z 1 μ z } ) ]
V ˙ 2 z = sign σ z D α [ D α z { 1 T c z μ z 2 e ( | σ z | μ z 2 ) σ z 1 μ z 2 + K z sign ( σ z ) } ] = sign σ z [ { 1 T c z μ z 2 e ( | σ z | μ z 2 ) σ z 1 μ z 2 + K z sign ( σ z ) } ]
Because of sign ( σ ) sign ( σ ) = 1 , (35) can be achieved as follows:
V ˙ 2 z = sign σ z [ { 1 T c z μ z 2 e ( | σ z | μ z 2 ) | σ z | 1 μ z 2 sign σ z + K z sign ( σ z ) } ] = 1 T c z μ z 2 e ( | σ z | μ z 2 ) | σ z | 1 μ z 2 | σ z | K z | σ z |
V ˙ 2 z = 1 T c z μ z 2 e ( | V 2 z | μ z 2 ) | V 2 z | 1 μ z 2 | σ z | K z | V 2 z | 0
It can be demonstrated using Theorems 1 and 2 that the consensus error e z will reach a stable zone in less time than the upper bound time t T = T s z + T c z .
As a result, the suggested PTFOSMC can stabilize the quadrotor system to the stable area in the upper bound of the given predefined time T c .
Remark 5. 
The position subsystem has three outputs, ( z ( t ) , x ( t ) , y ( t ) ) , but is controlled by only one single control. Hence, the proposed control scheme does not apply to every subsystem. Therefore, to regulate the horizontal position, backstepping with SMC is used.

4.3. Horizontal and Vertical Control Design

In this subsection, the virtual control signals, Π x and Π y , will be created by backstepping in conjunction with SMC to provide the desired angles and thrust. The suggested control strategy for the altitude subsystem can ensure that it quickly converges to the required altitude within the predefined time. The sliding mode surfaces for the x and y subsystems can be defined as:
σ x ( t ) = e ˙ x ( t ) + λ x e x ( t ) σ y ( t ) = e ˙ y ( t ) + λ y e y ( t )
Therefore, the time-derivative of these surfaces is as follows:
σ ˙ x ( t ) = x ¨ x ¨ d + λ x e ˙ x ( t ) σ ˙ y ( t ) = y ¨ y ¨ d + λ y e ˙ y ( t )
Using the controllers created in [33,35], the following laws can be applied to control the horizontal position:
Π x = δ 9 x ˙ + d x + x ¨ d ( t ) e ˙ x ( t ) λ x e x ( t ) K x sign σ x ( t ) Π y = δ 10 y ˙ + d y + y ¨ d ( t ) e ˙ y ( t ) λ y e y ( t ) K y sign σ y ( t )
where λ x ; λ y ; K x ; K y > 0 . The predefined-time FOSMC for the attitude subsystem can be designed using the same procedures as those described for altitude in the following subsection.

4.4. Inner-Loop Control Design

The primary issue investigated in this paragraph is how to provide the rolling, pitching, and yawing control inputs in such a way that the rotating subsystem is predefined-time stable.
The following formulas provide the tracking errors for the ϕ ( t ) , θ ( t ) , and ψ ( t ) subsystems and the time-derivative can be expressed as follows:
e ϕ ( t ) = ϕ ( t ) + ϕ d ( t ) e θ ( t ) = θ ( t ) + θ d ( t ) e ψ ( t ) = ψ ( t ) + ψ d ( t ) , e ˙ ϕ ( t ) = ϕ ˙ ( t ) + ϕ ˙ d ( t ) e ˙ θ ( t ) = θ ˙ ( t ) + θ ˙ d ( t ) e ˙ ψ ( t ) = ψ ˙ ( t ) + ψ ˙ d ( t )
For subsystems ϕ ( t ) , θ ( t ) , and ψ ( t ) , the predefined-time fractional-order sliding mode manifold is given as
σ ϕ = I α ϕ { 1 T s ϕ μ ϕ e ( | e ϕ | μ ϕ ) e ϕ 1 μ ϕ } + γ ϕ D 1 α ϕ e ϕ ( t ) σ θ = I α θ { 1 T s θ μ θ e ( | e θ | μ θ ) e θ 1 μ θ } + γ θ D 1 α θ e θ ( t ) σ z = I α ψ { 1 T s ψ μ ψ e ( | e ψ | μ ψ ) e ψ 1 μ ψ } + γ ψ D 1 α ψ e ψ ( t )
where e Θ is the error tracking with Θ = [ ϕ , θ , ψ ] ; T s Θ denotes the predefined setting upper bound convergence time, and e Θ ( t ) μ Θ 2 = | e Θ | . sign ( e Θ ( t ) ) with 0 < μ Θ 2 < 1 . The following equations can be reached when the consensus error e Θ reaches the fractional-order sliding surface σ Θ = 0 . Considering the use of (41), (42) can be written as follows:
D 1 α ϕ e ϕ ( t ) = D α ϕ { 1 T s ϕ μ ϕ e ( | e ϕ | μ ϕ ) e ϕ 1 μ ϕ } D 1 α θ e θ ( t ) = D α θ { 1 T s θ μ θ e ( | e θ | μ θ ) e θ 1 μ θ } D 1 α ψ e ψ ( t ) = D α ψ { 1 T s ψ μ ψ e ( | e ψ | μ ψ ) e ψ 1 μ ψ }
The following equation can be obtained. The following equation can be obtained from properties (3) and (41):
e ¨ Θ ( t ) = D 1 + α Θ D 1 α Θ e Θ ( t ) = D 1 + α Θ D α Θ { 1 T s Θ μ Θ e ( | e Θ | μ Θ ) e Θ 1 μ Θ } = D 1 { 1 T s Θ μ Θ e ( | e Θ | μ Θ ) e Θ 1 μ Θ }
Therefore, the equivalent controller Π j with j = 2 , 3 , 4 can be obtained as follows:
Θ ¨ Θ ¨ d = D 1 { 1 T s Θ μ Θ e ( | e Θ | μ Θ ) e Θ 1 μ Θ }
Π j e q = D 1 { 1 T s Θ μ Θ e ( | e Θ | μ Θ ) e Θ 1 μ Θ } f ( Θ )
By substituting Equation (45) into Equation (7), the equivalent controller of attitude can be written as
Π 2 e q = J x D 1 { 1 T s ϕ μ ϕ e ( | e ϕ | μ ϕ ) e ϕ 1 μ ϕ } δ 1 θ ˙ ψ ˙ δ 2 ω r θ ˙ 2 δ 3 ϕ ˙ d ϕ Π 3 e q = J y D 1 { 1 T s ϕ μ ϕ e ( | e ϕ | μ ϕ ) e ϕ 1 μ ϕ } δ 4 ϕ ˙ ψ ˙ δ 5 ω r ϕ ˙ 2 δ 6 θ ˙ d θ Π 4 e q = J z D 1 { 1 T s ψ μ ψ e ( | e ϕ | μ ψ ) e ϕ 1 μ ϕ } δ 7 ϕ ˙ θ ˙ δ 8 ψ ˙ d ψ
System (7) is then constructed to use the switching controller, u s , to reach the PTFO sliding manifolds.
Π Θ s = D α Θ { 1 T c Θ μ Θ 2 e ( | σ Θ | μ Θ 2 ) σ Θ 1 μ Θ 2
+ K Θ sign ( σ Θ ) }
where T c Θ > 0 and μ Θ 2 , K Θ are positive parameters with Θ = [ ϕ , θ , ψ ] . Combining (46) and (48), the following law determines the altitude dynamics of the quadrotor:
Π Θ = Π Θ e q + Π Θ s = D 1 { 1 T s Θ μ Θ e ( | e Θ | μ Θ ) e Θ 1 μ Θ } f ( Θ ) D α Θ { 1 T c Θ μ Θ 2 e ( | σ Θ | μ Θ 2 ) σ Θ 1 μ Θ 2 + K Θ sign ( σ Θ ) }
The corresponding PTFOSM controller for the quadrotor attitude is given as
Π 2 = J x D 1 { 1 T s ϕ μ ϕ e ( | e ϕ | μ ϕ ) e ϕ 1 μ ϕ } δ 1 θ ˙ ψ ˙ δ 2 ω r θ ˙ 2 δ 3 ϕ ˙ + K ϕ sign ( σ ϕ ) } d ϕ + ϕ ¨ d D α ϕ { 1 T c ϕ μ ϕ 2 e ( | σ ϕ | μ ϕ 2 ) σ ϕ 1 μ ϕ 2
Π 3 = J y D 1 { 1 T s ϕ μ ϕ e ( | e ϕ | μ ϕ ) e ϕ 1 μ ϕ } δ 4 ϕ ˙ ψ ˙ δ 5 ω r ϕ ˙ 2 δ 6 θ ˙ d θ + K θ sign ( σ θ ) } d θ + ϕ ¨ d D α θ { 1 T c θ μ θ 2 e ( | σ θ | μ θ 2 ) σ θ 1 μ θ 2
Π 4 = J z D 1 { 1 T s ψ μ ψ e ( | e ϕ | μ ψ ) e ϕ 1 μ ϕ } δ 7 ϕ ˙ θ ˙ δ 8 ψ ˙ + K ψ sign ( σ ψ ) } d ψ + ϕ ¨ d D α ψ { 1 T c ψ μ ψ 2 e ( | σ ψ | μ ψ 2 ) σ ψ 1 μ ψ 2
Remark 6. 
The reaching law proposed in this paper has two main advantages: (1) it guarantees the predefined time convergence of the sliding manifold and (2) it copes with complex disturbances.

4.5. Attitude Stability Analysis for a Quadrotor

Theorem 3. 
The errors will achieve convergence to the origin within a predefined time for the fractional-order sliding surface, and stability for the yaw subsystem is attained as the state errors are fulfilled.
Proof. 
The selection of the Lyapunov function for the quadrotor system is as follows:
V 1 ψ = e ψ
The time differentiation of V 1 i is achieved as follows:
V ˙ 1 ψ = sign e ψ e ˙ ψ = sign e ψ D α D 1 α e ψ = sign e ψ D α D α ψ 1 T s ψ μ ψ e ( | e ψ | μ ψ ) e ψ 1 μ ψ = sign e ψ 1 T s ψ μ ψ e ( | e ψ | μ ψ ) e ψ 1 μ ψ
Due to sign e ψ × sign e ψ = 1 , (18) can be achieved as follows:
V ˙ ψ 1 = sign e ψ 1 T s ψ μ ψ e ( | e ψ | μ ψ ) sign e i | e ψ | 1 μ ψ = 1 T s ψ μ ψ e ( | e ψ | μ ψ ) | e ψ | 1 μ ψ = 1 T s ψ μ ψ e ( | V 1 ψ | μ ψ ) | V 1 ψ | 1 μ ψ 0
Based on Lemma 2, it is demonstrated that the dynamic system (system (55)) can converge to zero in a predefined time T s ψ .
The fractional-order sliding surface σ ψ will converge to zero if the dynamic system in (42) is stable.
Theorem 4. 
Consider the PTFO sliding manifold designed in (42); the convergence to zero in a predefined time of the dynamics system defined in (7) is ensured.
Proof. 
To prove Theorem 4, the Lyapunov function candidate of the yaw- subsystem is considered as follows:
V 3 = σ ψ V ˙ 3 = sign σ ψ σ ˙ ψ = sign σ ψ D α D 1 + α σ ψ
V ˙ 2 ψ = sign σ ψ D α e ¨ ψ ( t ) + D 1 { 1 T s ψ μ ψ e ( | e ψ | μ ψ ) e ψ 1 μ ψ }
Based on (7) and the PTFOSMC (49), (32) can be yielded as follows:
V ˙ 2 ψ = sign σ ψ σ ˙ ψ = sign σ ψ D α D 1 + α σ i = sign σ ψ D α [ δ 7 ϕ ˙ θ ˙ + δ 8 ψ ˙ + d ψ + J z 1 Π 4 ψ ¨ d + D 1 { 1 T s ψ μ ψ e ( | e ψ | μ ψ ) e ψ 1 μ ψ } ) ]
V ˙ 2 ψ = sign σ ψ D α [ δ 7 ϕ ˙ θ ˙ + δ 8 ψ ˙ + d ψ ψ ¨ d δ 7 ϕ ˙ θ ˙ δ 8 ψ ˙ d ψ + ψ ¨ d + J z 1 [ J z D 1 { 1 T s ψ μ ψ e ( | e ϕ | μ ψ ) e ϕ 1 μ ϕ } K ψ sign ( σ ψ ) } D α ψ { 1 T c ψ μ ψ 2 e ( | σ ψ | μ ψ 2 ) σ ψ 1 μ ψ 2 ] + D 1 { 1 T s ψ μ ψ e ( | e ψ | μ ψ ) e ψ 1 μ ψ } ) ] = sign σ ψ D α ψ [ K ψ sign ( σ ψ ) + D α ψ { 1 T c ψ μ ψ 2 e ( | σ ψ | μ ψ 2 ) σ ψ 1 μ ψ 2 ] = K ψ | σ ψ | D α ψ { 1 T c ψ μ ψ 2 e ( | σ ψ | μ ψ 2 ) | σ ψ | 1 μ ψ 2 | σ ψ | }
V ˙ 2 ψ = K ψ | V 3 | D α ψ { 1 T c ψ μ ψ 2 e ( | V 3 | μ ψ 2 ) | V 3 | 1 μ ψ 2 | V 3 | } 0 .
It can be demonstrated using Theorems 3 and 4 that the consensus error e ψ will reach a stable zone in less time than the upper bound time t T 2 = T s ψ + T c ψ .
It is clear from the analysis above that the reaching condition of attitude loop stability is ensured.
Example 1. 
The suggested controller is used to stabilize the yaw subsystem using the settings listed in Table 1 and Table 2. The examples in Figure 3 and Figure 4 represent the simulation results for this example. The sliding manifold and the state variable of the yaw subsystem both converge to their origins within the predefined time T f , as can be seen from these results. As illustrated in Figure 3, the smaller the parameters, T c and T s , the less serious the control, indicating that choosing the appropriate T c and T s is important in practical implementations of the recommended control approach, as in Figure 3.

4.6. Stability Global Analysis

Theorem 5. 
Fundamental control law Equations (36) and (59) applied to the investigated system (7) guarantee the overall quadrotor system stability.
Proof. 
The Lyapunov function for the quadrotor system and its time-derivative are given as follows:
V g = V 2 z + V 2 ϕ + V 2 θ + V 2 ψ V ˙ g = V ˙ 2 z + V ˙ 2 ϕ + V ˙ 2 θ + V ˙ 2 ψ
From Equations (36) and (60), we have
V g = 1 T c z μ z 2 e ( | σ z | μ z 2 ) | σ z | 1 μ z 2 | σ z | K z | σ z | D α ϕ { 1 T c ϕ μ ϕ 2 e ( | σ ϕ | μ ϕ 2 ) | σ ϕ | 1 μ ϕ 2 | σ ϕ | } D α θ { 1 T c θ μ θ 2 e ( | σ θ | μ θ 2 ) | σ θ | 1 μ θ 2 | σ θ | } D α ψ { 1 T c ψ μ ψ 2 e ( | σ ψ | μ ψ 2 ) | σ ψ | 1 μ ψ 2 | σ ψ | } K ϕ | σ ϕ | K θ | σ θ | K ψ | σ ψ | 0 .
According to Theorems 1–4, Theorem 5 demonstrates the global stability achieved for both rotational and translational tracking errors through the Lyapunov approach. □
Remark 7. 
The quadrotor system control has two loops, where the inner loop (attitude loop) is faster compared to the outer loop (position loop). The proposed PTFOSMC has been successfully applied to the inner loop to achieve predefined time performances and the z-subsystem has benefited from the proposed controller.
Remark 8. 
The proposed scheme offers several advantages, including the predefined-time convergence of system states, predefined-time convergence in the sliding mode phase, mitigation of the chattering problem, resilience to bounded external disturbances, and heightened robustness.

5. Simulations Results for Quadrotor

In this section, numerical simulations are used to assess the effectiveness of the suggested control approach, which is based on the PTFOSMC method for the path-following problem. The physical characteristics of a quadrotor are listed in Table 1. In Table 2, the PTFOSMC parameters are displayed.
Remark 9. 
The resolution to the chattering problem involves transforming the discontinuous element function into a continuous function, such as the hyperbolic tangent function.

5.1. Scenario 1

In this case, the quadrotor tracks the desired reference without external disturbances. Equation (63) provides more particular information about the desired trajectory.
x d = 0.5 cos π 20 m , y d = 0.5 sin π 20 m z d = 2 2 cos π 2 m ψ d = 0.4 r a d t [ 0 50 ] a n d ψ d = 0.2 r a d t [ 50 80 ]
The tracking abilities of the PTFOSMC without disturbances are displayed in Figure 5, Figure 6 and Figure 7. A good tracking trajectory is produced by the position subsystem’s suggested control strategy, as shown in Figure 5. Moreover, the tracking performance of the yaw, pitch, and roll angles as they converge to their origin values in a predefined time is presented in Figure 6. As seen in this Figure, the yaw, for example, converges to the origin before the predefined time T s ψ = T s ψ 1 + T s ψ 2 = 0.3 s. In addition, one can see that controller u remains continuous and the sliding mode σ ψ = 0 is required for every t T s ψ 2 . Figure 7 shows an excellent tracking trajectory of the desired control without external disturbances.

5.2. Scenario 2

This section is used to assess the stability of a quadrotor when it encounters time-varying disturbances. Quadrotor model equations are expanded to include disturbances. One assumes that the disturbances should be bounded. When the quadrotor flies outside, wind gusts cause these disruptions, which are accelerations. As a result, disturbances are shown in the following ways:
d x = 0.5 sin ( 0.4 t ) m / s 2 10 < t < 30 d y = 0.2 sin ( 0.4 t ) + 0.3 cos ( 0.7 t ) m / s 2 10 < t < 50 d z = 0.4 cos ( 0.7 t ) m / s 2 t > 20 d ϕ = 0.3 cos ( 0.4 t ) rad / s 2 t < 10 d θ = 0.3 cos ( 0.5 t ) rad / s 2 10 < t < 20 d ψ = 0.3 cos ( 0.7 t ) rad / s 2 t > 20
Equations (65) and (66) provide the desired trajectory.
x d = 0.5 cos ( t / 2 ) m t [ 0 , 4 π ] x d = 0.5 m t [ 4 π , 20 ] x d = 0.25 t 4.5 m t [ 20 , 30 ] x d = 3 m t [ 30 , 80 ] , y d = 0.5 sin ( t / 2 ) m t [ 0 , 4 π ] y d = 0.25 t 3.14 m t [ 4 π , 20 ] y d = 5 π m t [ 20 , 30 ] y d = 0.2358 t + 8.94 m t [ 4 π , 20 ] y d = 0.5 m t [ 30 , 80 ]
z d = 0.125 t + 1 m t [ 0 , 4 π ] z d = 0.5 π + 1 m t [ 4 π , 40 ] z d = exp ( 0.2 t + 8.944 ) m t [ 40 , 80 ] , ψ d = 0.4 r a d t [ 0 , 50 ] ψ d = 0.2 r a d t [ 50 , 80 ]
The tracking results obtained in scenario 2 (Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13) demonstrate the effectiveness of the control strategy suggested in this work. The proposed nonlinear control law was found to perform better than the ABSMC and PIDSMC control laws developed in [9,10] for the control of present external disturbances during the flight and following the desired trajectory. The performance attitude and position in terms of tracking performance are shown in Figure 8 and Figure 9. In this regard, the proposed control system provides a shorter settling time than various other controls and we can see that the position and attitude controllers are accurate in tracking the correct values even when the references change quickly. We present the effectiveness of the control inputs of the quadrotor ( Π 1 , Π 2 , Π 3 , Π 4 ) . As can be seen, the smooth convergence of the control signals to their original values ( 7.2598 , 0 , 0 , 0 ) illustrates the efficiency of the PTFOSMC technique. The tracking sliding manifolds made by the system are shown in Figure 11 and Figure 12. These data demonstrate that the tracking surfaces maintain their status of zero. The 3 D flight trajectory is depicted in Figure 13. The results presented above show that, in comparison to previous techniques, the PTFOSMC has successfully tracked the desired trajectory.

5.3. Scenario 3

To further evaluate the effectiveness of the PTFOSMC method, the disturbances mentioned in scenario 2 are added to all state variables. In the current case, the initial conditions are changed; for the yaw subsystem, ψ 0 = 0.2 rad, ψ 0 = 0.4 rad, and ψ 0 = 0.8 rad, for the z-subsystem, z 0 = 0 m, z 0 = 0.2 m, z 0 = 0.5 m, and z 0 = 0.8 m.
The response of the quadrotor system provides evidence that the suggested controller induces a fast response. It should also be noted that there is an upper constraint on the convergence time, which can be predefined as an adjustable control parameter and is independent of the initial conditions, as shown in Figure 14.

5.4. Quantitative Analysis of the Controllers

The preceding curves are used to demonstrate the differences between the controllers because all the strategies produce acceptable results in the ideal scenario. As a result, the performance indices, ISE and ISV, are provided as objective functions to be utilized in the next part.
  • Integral square error (ISE): It is given by ISE = t 0 e t e 2 ( t ) d t
  • Integral of the square value of control signal (ISV): ISV = t 0 t t Π 2 ( t ) d t
Table 3 shows the ISE performance of the three controllers. The proposed method is superior to previous techniques since it can be established through the PTFOSM control that the ISE indices are lower than BSMC and PIDSMC. Moreover, we also notice that the proposed control uses less energy than PIDSMC and the other way around.
The three controllers demonstrate satisfactory performance with moderate energy consumption, regardless of external influences. However, the presence of an additional payload, such as a gust of wind, necessitates extra thrust to ensure optimal performance. As evident from Table 3 and in the second scenario, the PIDSMC technique exhibits the least accuracy ( I S E ψ = 0.01622 ) , followed by the classical MFC technique, ( I S E ψ = 0.02264 ) , whereas PTFOSMC proves to be the most accurate ( I S E ψ = 0.009247 ) . Additionally, it is observed that PIDSMC consumes more energy, while the proposed control consumes less energy.
Table 4 presents the ISV comparative results for the inputs of a quadrotor. It is clear that the performance index values for the proposed control methodology are significantly lower than those for the alternative strategy. The tracking performance of the suggested control approach is then found to be better than that of the method of [9,10] by comparing the simulation results.
Remark 10. 
The proposed control approach was compared to the BSMC and PIDSMC techniques in addressing the quadrotor system tracking issue. While tracking errors in other approaches may only converge to a bounded neighborhood of the origin within a predefined time due to the use of an error as a sliding variable, the suggested control scheme, incorporating a fractional-order manifold, ensures the predefined-time zero-error stability of quadrotor systems, leading to enhanced steady-state performance.
Remark 11. 
In this paper, two nonlinear surfaces are introduced to achieve the finite-time convergence of the quadrotor states.
Remark 12. 
The benefits of the suggested fractional-order predefined-time sliding mode controllers include predefined-time convergence of system states, predefined-time convergence in the reaching phase, alleviation of the chattering phenomenon, and the incorporation of an additional fractional-order differential element that remains insensitive to bounded external disturbances.
The implementation challenges of the proposed control method for the quadrotor’s position and attitude will be examined through four simulations. The complexity of the problem arises from addressing tracking performance, which includes convergence time and steady-state performance for both inner and outer loops, especially in the presence of external disturbances. Moreover, the innovation of the proposed control lies in the formulation of a predefined-time controller based on fractional-order sliding mode control, enhancing the robustness and precision of the quadrotor system.

6. Conclusions

This work suggests a novel control method for tracking quadrotor trajectory. Novel predefined time in fractional-order sliding mode manifolds are developed to stabilize tracking errors to the origin within a predefined time. A systematic evaluation is conducted on the quadrotor system to verify specific conditions on the control parameters, ensuring the stability of zero tracking error. The Newton–Euler method is applied in this research to determine the dynamics of the quadrotor. The SMC can ensure that the quadrotor system will obtain better dynamic performances due to its benefits of quick transient response, high resilience, and simple implementation. The proposed strategy is employed to govern a quadrotor system across various scenarios influenced by escalating disturbance amplitudes and frequencies. Additionally, the PTFOSMC is designed using fractional calculus to reduce the inherent chattering problem of the standard SMC and improve parameter freedom. The simulation findings demonstrate that the quadrotor system using the proposed PTFOSMC might exceed existing control strategies with faster convergence time and higher robustness.

Author Contributions

Conceptualization, A.B., M.L. and S.B.; Methodology, A.B., M.L., S.B. and F.S.A.; Software, A.B. and M.L.; Validation, M.L., S.B., F.S.A. and M.B.; Formal analysis, A.B., M.L. and F.S.A.; Writing—original draft, A.B. and M.L.; Visualization, M.L., S.B. and M.B.; Supervision, M.L., F.S.A. and M.B.; funding acquisition, S.B. and F.S.A. All authors have read and agreed to the published version of the manuscript.

Funding

Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia—project number MoE-IF-UJ-22-04220772-6.

Data Availability Statement

No data were used to support this study.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through project number MoE-IF-UJ-22-04220772-6.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hassanalian, M.; Abdelkefi, A. Classifications, Applications, and Design Challenges of Drones: A Review. Prog. Aerosp. Sci. 2017, 91, 99–131. [Google Scholar] [CrossRef]
  2. Watts, A.C.; Ambrosia, V.G.; Hinkley, E. Unmanned Aircraft Systems in Remote Sensing and Scientific Research: Classification and Considerations of Use. Remote Sens. 2012, 4, 1671–1692. [Google Scholar] [CrossRef]
  3. de Zarzà, I.; de Curtò, J.; Cano, J.C.; Calafate, C.T. Drone-Based Decentralized Truck Platooning with UWB Sensing and Control. Mathematics 2023, 11, 4627. [Google Scholar] [CrossRef]
  4. Eskandaripour, H.; Boldsaikhan, E. Last-Mile Drone Delivery: Past, Present, and Future. Drones 2023, 7, 77. [Google Scholar] [CrossRef]
  5. Gohari, A.; Ahmad, A.; Rahim, R.A.; Supa’at, A.S.M.; Razak, S.A.; Gismalla, M.S.M. Involvement of Surveillance Drones in Smart Cities: A Systematic Review. IEEE Access 2022, 10, 56611–56628. [Google Scholar] [CrossRef]
  6. Mishra, B.; Garg, D.; Narang, P.; Mishra, V.K. Drone-surveillance for Search and Rescue in Natural Disaster. Comput. Commun. 2020, 156, 1–10. [Google Scholar] [CrossRef]
  7. Labbadi, M.; Cherkaoui, M. Robust Adaptive Nonsingular Fast Terminal Sliding-mode Tracking Control for an Uncertain Quadrotor UAV Subjected to Disturbances. ISA Trans. 2020, 99, 290–304. [Google Scholar] [CrossRef] [PubMed]
  8. Mechali, O.; Xu, L.; Xie, X.; Iqbal, J. Theory and Practice for Autonomous Formation Flight of Quadrotors via Distributed Robust Sliding Mode Control Protocol with Fixed-time Stability Guarantee. Control Eng. Pract. 2022, 123, 105150. [Google Scholar] [CrossRef]
  9. Labbadi, M.; Cherkaoui, M. Robust Adaptive Backstepping Fast Terminal Sliding Mode Controller for Uncertain Quadrotor UAV. Aerosp. Sci. Technol. 2019, 93, 105306. [Google Scholar] [CrossRef]
  10. Comert, C.; Kasnakoglu, C. Comparing and Developing PID and Sliding Mode Controllers for Quadrotor. Int. J. Mech. Eng. Robot. Res. 2017, 6, 194–199. [Google Scholar] [CrossRef]
  11. Elyaalaoui, K.; Labbadi, M.; Boubaker, S.; Kamel, S.; Alsubaei, F.S. On Novel Fractional-Order Trajectory Tracking Control of Quadrotors: A Predefined-Time Guarantee Performance Approach. Mathematics 2023, 11, 3582. [Google Scholar] [CrossRef]
  12. Raffo, G.V.; Ortega, M.G.; Rubio, F. An Integral Predictive/nonlinear H∞ Control Structure for a Quadrotor Helicopter. Automatica 2010, 46, 29–39. [Google Scholar] [CrossRef]
  13. Zheng, E.-H.; Xiong, J.-J.; Luo, J.-L. Second Order Sliding Mode Control for a Quadrotor UAV. ISA Trans. 2014, 53, 1350–1356. [Google Scholar] [CrossRef] [PubMed]
  14. Tian, B.; Tian, B.; Lu, H.; Zuo, Z.; Zong, Q.; Zhang, Y. Multivariable Finite-time Output Feedback Trajectory Tracking Control of Quadrotor Helicopters. Int. Robust Nonlinear Control 2018, 28, 281–295. [Google Scholar] [CrossRef]
  15. Xie, W.; Cabecinhas, D.; Cunha, R.; Silvestre, C. Adaptive Backstepping Control of a Quadcopter with Uncertain Vehicle Mass, Moment of Inertia, and Disturbances. IEEE Trans. Ind. Electron. 2022, 69, 549–559. [Google Scholar] [CrossRef]
  16. Guo, Y.; Ma, B. Global Sliding Mode with Fractional Operators and Application to Control Robot Manipulators. Int. J. Control 2019, 92, 1497–1510. [Google Scholar] [CrossRef]
  17. Ni, J.; Liu, L.; Liu, C.; Hu, X.; Li, S. Fast Fixed-time Nonsingular Terminal Sliding Mode Control and Its Application to Chaos Suppression in Power System. Circuits Syst. II Exp. Briefs 2017, 64, 151–155. [Google Scholar] [CrossRef]
  18. Wang, Z.; Wang, J.; La Scala, M. A Novel Distributed-decentralized Fixed-time Optimal Frequency and Excitation Control Framework in a Nonlinear Network-preserving Power System. IEEE Trans. Power Syst. 2021, 36, 1285–1297. [Google Scholar] [CrossRef]
  19. Zeng, T.; Ren, X.; Zhang, Y. Fixed-time Sliding Mode Control and High-gain Nonlinearity Compensation for Dual-motor Driving System. IEEE Trans. Ind. Inform. 2020, 16, 4090–4098. [Google Scholar] [CrossRef]
  20. Muñoz-Vázquez, A.J.; Sanchez-Torres, J.D.; Jiménez-Rodríguez, E.; Loukianov, A.G. Predefined-time Robust Stabilization of Robotic Manipulators. IEEE/ASME Trans. Mechatron. 2019, 24, 2709–2722. [Google Scholar] [CrossRef]
  21. Muñoz-Vázquez, A.J.; Fernández-Anaya, G.; Sanchez-Torres, J.D.; Meléndez-Vázquez, F. Predefined-time Control of Distributed-order Systems. Nonlinear Dyn. 2021, 103, 2689–2700. [Google Scholar] [CrossRef]
  22. Ni, J.; Liu, L.; Tang, Y.; Liu, C. Predefined-time Consensus Tracking of Second-order Multiagent Systems. IEEE Trans. Syst. 2021, 51, 2550–2560. [Google Scholar] [CrossRef]
  23. Benaddy, A.; Labbadi, M.; Elyaalaoui, K.; Bouzi, M. Fixed-Time Fractional-Order Sliding Mode Control for UAVs under External Disturbances. Fractal Fract. 2023, 7, 775. [Google Scholar] [CrossRef]
  24. Muñoz-Vázquez, A.J.; Parra-Vega, V.; Sánchez-Orta, A.; Romero-Galván, G. Fractional PD IλD Error Manifolds for Robust Tracking Control of Robotic Manipulators. J. Dyn. Syst. Meas. Control 2019, 141, 031006-6. [Google Scholar] [CrossRef]
  25. Wang, N.; Deng, Q.; Xie, G.; Pan, X. Hybrid Finite-time Trajectory Tracking Control of a Quadrotor. ISA Trans. 2019, 90, 278–286. [Google Scholar] [CrossRef] [PubMed]
  26. Huang, S.; Wang, J.; Xiong, L.; Liu, J.; Li, P.; Wang, Z. Distributed Predefined-Time Fractional-Order Sliding Mode Control for Power System With Prescribed Tracking Performance. IEEE Trans. Power Syst. 2022, 37, 2233–2246. [Google Scholar] [CrossRef]
  27. Bhat, S.P.; Bernstein, D.S. Finite-time Stability of Continuous Autonomous Systems. SIAM J. Control Optim. 2000, 38, 1430–1443. [Google Scholar] [CrossRef]
  28. Polyakov, A. Nonlinear Feedback Design for Fixed-time Stabilization of Linear Control Systems. IEEE Trans. Autom. Control 2012, 57, 2106–2110. [Google Scholar] [CrossRef]
  29. Fraguela, L.; Angulo, M.T.; Moreno, J.A.; Fridman, L. Design of a prescribed convergence time uniform Robust Exact Observer in the presence of measurement noise. In Proceedings of the IEEE 51st IEEE Conference on Decision and Control, Maui, HI, USA, 10–13 December 2012. [Google Scholar] [CrossRef]
  30. Sanchez-Torres, J.D.; Sanchez-Torres, J.D.; Gómez-Gutiérrez, D.; López, E.; Loukianov, A.G. A Class of Predefined-time Stable Dynamical Systems. IMA J. Math. Control Inf. 2018, 35, i1–i29. [Google Scholar] [CrossRef]
  31. Pourmahmood Aghababa, M. A Novel Terminal Sliding Mode Controller for a Class of Non-autonomous Fractional-order Systems. Nonlinear Dyn. 2013, 73, 679–688. [Google Scholar] [CrossRef]
  32. Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA; Available online: https://shop.elsevier.com/books/fractional-differential-equations/podlubny/978-0-12-558840-9 (accessed on 1 January 1999).
  33. Benaddy, A.; Labbadi, M.; Bouzi, M. Adaptive Nonlinear Controller for the Trajectory Tracking of the Quadrotor with Uncertainties. In Proceedings of the 2020 2nd Global Power, Energy and Communication Conference (GPECOM), Izmir, Turkey, 20–23 October 2020. [Google Scholar] [CrossRef]
  34. Liang, C.-D.; Ge, M.-F.; Ge, M.-F.; Liu, Z.-W.; Ling, G.; Liu, F.; Liu, F. Predefined-time Formation Tracking Control of Networked Marine Surface Vehicles. Control Eng. Pract. 2021, 107, 104682. [Google Scholar] [CrossRef]
  35. Chen, F.; Jiang, R.; Zhang, K.; Jiang, B.; Tao, G. Robust Backstepping Sliding-mode Control and Observer-based Fault Estimation for a Quadrotor UAV. IEEE Trans. Ind. Electron. 2016, 63, 5044–5056. [Google Scholar] [CrossRef]
Figure 1. The four-rotor drone arrangement.
Figure 1. The four-rotor drone arrangement.
Mathematics 11 04886 g001
Figure 2. The proposed improved PTFOSMC method is illustrated in a block diagram.
Figure 2. The proposed improved PTFOSMC method is illustrated in a block diagram.
Mathematics 11 04886 g002
Figure 3. Result e ψ , σ ψ , and Π 4 with T c = 0.05 and T s = 0.2 .
Figure 3. Result e ψ , σ ψ , and Π 4 with T c = 0.05 and T s = 0.2 .
Mathematics 11 04886 g003
Figure 4. Results e ψ , σ ψ , and Π 4 with T c = 0.4 and T s = 0.5 .
Figure 4. Results e ψ , σ ψ , and Π 4 with T c = 0.4 and T s = 0.5 .
Mathematics 11 04886 g004
Figure 5. Quadrotor position without disturbances.
Figure 5. Quadrotor position without disturbances.
Mathematics 11 04886 g005
Figure 6. Quadrotor attitude without disturbances.
Figure 6. Quadrotor attitude without disturbances.
Mathematics 11 04886 g006
Figure 7. The 2D and 3D trajectories of a quadrotor without disturbances.
Figure 7. The 2D and 3D trajectories of a quadrotor without disturbances.
Mathematics 11 04886 g007
Figure 8. Position of a quadrotor.
Figure 8. Position of a quadrotor.
Mathematics 11 04886 g008
Figure 9. Attitude of a quadrotor.
Figure 9. Attitude of a quadrotor.
Mathematics 11 04886 g009
Figure 10. Inputs a quadrotor.
Figure 10. Inputs a quadrotor.
Mathematics 11 04886 g010
Figure 11. Sliding surface for the position subsystem.
Figure 11. Sliding surface for the position subsystem.
Mathematics 11 04886 g011
Figure 12. Sliding surface for the attitude subsystem.
Figure 12. Sliding surface for the attitude subsystem.
Mathematics 11 04886 g012
Figure 13. The 3D trajectories of a quadrotor.
Figure 13. The 3D trajectories of a quadrotor.
Mathematics 11 04886 g013
Figure 14. Initial values of a quadrotor.
Figure 14. Initial values of a quadrotor.
Mathematics 11 04886 g014
Table 1. Quadrotor physical characteristics.
Table 1. Quadrotor physical characteristics.
SymbolValueSymbolValue
m ( kg ) 0.74 k x (m 1 · s · N ) 5.5670 × 10 4
g (s 2 · m) 9.81 k y (m 1 · s · N ) 5.5670 × 10 4
l ( m ) 0.5 k z (m 1 · s · N ) 5.5670 × 10 4
J r ( kg · m 2 )2.03 × 10 5 K a x (rad 1 · s · N ) 5.5670 × 10 4
J x ( kg · m 2 ) 0.004 K a y (rad 1 · s · N ) 5.5670 × 10 4
J y ( kg · m 2 ) 0.004 K a z (rad 1 · s · N ) 5.5670 × 10 4
J z ( kg · m 2 ) 0.0084 ♭ (s 2 · N ) 2.984 × 10 3
Table 2. Control parameters.
Table 2. Control parameters.
ParameterValueParameterValue
T c 1 , c 3 , c 5 0.8 T c z 2
T s ϕ , s θ , s ψ 0.5 T s z 1
k ϕ , θ , ψ 10.51 k z 5
γ ϕ , θ , ψ 20 γ z 4
μ ϕ , θ , ψ 0.2 μ z 0.9
μ ϕ 2 , θ 2 , ψ 2 0.4 μ z 2 0.9
Table 3. Performance index of ISE.
Table 3. Performance index of ISE.
ISE of the Scenario 1
VariableProposed MethodMethod [5]Method [6]
ϕ ( t ) 0.0016980.0033140.03161
θ ( t ) 0.0044480.01050.03038
ψ ( t ) 0.0092450.22620.01622
z ( t ) 0.20360.36460.4262
ISE of the Scenario 2
VariableProposed MethodMethod [5]Method [6]
ϕ ( t ) 0.0002020.0009480.001744
θ ( t ) 0.0063340.0035760.007545
ψ ( t ) 0.0092470.0226400.016220
z ( t ) 4.26 × 10 6 3.76 × 10 6 0.000410
Table 4. Performance index of ISV.
Table 4. Performance index of ISV.
ISV of the Scenario 1
VariableProposed MethodMethod [5]Method [6]
Total thrust7.2597.267.3
Total torques0.0031530.003772−0.003348
ISV of the Scenario 2
VariableProposed MethodMethod [5]Method [6]
Total thrust7.2597.2657.32
Total torques0.00021270.002856−0.0002931
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

Benaddy, A.; Labbadi, M.; Boubaker, S.; Alsubaei, F.S.; Bouzi, M. Predefined-Time Fractional-Order Tracking Control for UAVs with Perturbation. Mathematics 2023, 11, 4886. https://doi.org/10.3390/math11244886

AMA Style

Benaddy A, Labbadi M, Boubaker S, Alsubaei FS, Bouzi M. Predefined-Time Fractional-Order Tracking Control for UAVs with Perturbation. Mathematics. 2023; 11(24):4886. https://doi.org/10.3390/math11244886

Chicago/Turabian Style

Benaddy, Abdellah, Moussa Labbadi, Sahbi Boubaker, Faisal S. Alsubaei, and Mostafa Bouzi. 2023. "Predefined-Time Fractional-Order Tracking Control for UAVs with Perturbation" Mathematics 11, no. 24: 4886. https://doi.org/10.3390/math11244886

APA Style

Benaddy, A., Labbadi, M., Boubaker, S., Alsubaei, F. S., & Bouzi, M. (2023). Predefined-Time Fractional-Order Tracking Control for UAVs with Perturbation. Mathematics, 11(24), 4886. https://doi.org/10.3390/math11244886

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop