Next Article in Journal
UAV Control Based on Pattern Recognition in Aquaculture Application
Previous Article in Journal
Artificial Intelligence Approach in Aerospace for Error Mitigation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Event-Triggered Prescribed Performance Robust Control for Aggressive Quadrotor Flight

1
School of Aerospace Engineering, Beijing Institute of Technology, Beijing 100081, China
2
Department of Computer Science and Technology, Tsinghua University, Beijing 100084, China
*
Author to whom correspondence should be addressed.
Aerospace 2024, 11(4), 301; https://doi.org/10.3390/aerospace11040301
Submission received: 3 March 2024 / Revised: 2 April 2024 / Accepted: 9 April 2024 / Published: 11 April 2024

Abstract

:
Aggressive flight has become increasingly important for expanding the applications of quadrotors. The typical characteristic of large and rapid changes in commands poses stringent demands on the maneuverability of quadrotors. Ensuring flight stability alone is not enough; dynamic responses must also be selectively constrained, presenting quadcopter flight control with daunting challenges. The prescribed performance control (PPC) method is seen as having the potential to solve this problem by allowing for the constrained control of specified performance, leading to extensive research. However, its practical application still faces challenges, such as the system divergence caused by errors exceeding boundaries due to sudden command mutations. This paper presents a robust dynamic event-triggered PPC (DETPPC) method for an aggressive quadrotor flight. By assessing the direction and proximity of tracking errors approaching constraint boundaries, a dynamic event-triggered compensation mechanism for performance function boundaries is established to mitigate the divergence caused by error surpassing and to preserve preset control over the targeted metrics. Controllers were designed for both the translational and rotational subsystems of the quadrotor, and stability analysis was conducted based on Lyapunov functions. Simulation tests on agile trajectory tracking and abrupt attitude control were carried out, demonstrating the effectiveness of the proposed method.

1. Introduction

In recent years, the application scenarios of quadrotors have expanded into military reconnaissance, strike missions, and commercial aerial performances [1,2,3]. Unlike traditional low-speed or hovering flights, in such applications, commands can change on a large scale, rapidly, or even abruptly [4]. For instance, in military reconnaissance, the tracked target may suddenly evade, necessitating rapid changes in speed and direction by the quadrotor to maintain tracking. These are all typical characteristics of an aggressive flight, imposing stringent demands on the maneuverability of quadrotors and posing challenges to flight control.
During an aggressive flight, quadrotors experience high speeds, rapid attitude changes, and noticeable aerodynamic disturbances, leading to pronounced model nonlinearity and uncertainty [5]. This poses a formidable challenge for the design of quadrotor flight controllers, primarily due to the difficulty of balancing maneuverability and robustness. Existing flight control methods and research, such as backstepping-based control [6], sliding-mode-based control [7], adaptive-based control [8], and neural-network-based control [9], primarily focus on enhancing quadrotors’ robustness against model errors and external disturbances for stable flight [10]. However, they lack attention to and constraints on quadrotors’ transient performance. Combining adaptive nonlinear control with disturbance observers, some studies have been able to achieve finite-time error convergence and ultimately maintain uniformly bounded tracking. However, their convergence speed heavily relies on the initial error. Especially when the initial error is substantial, the convergence time is inevitably longer [11].
To achieve a better transient response aligned with specified constraints, prescribed performance control (PPC) has gradually emerged in quadrotor flight control [12]. It utilizes a prescribed performance function to guide tracking errors to converge in a predetermined residual region within defined constraints, such as maximum overshoot, convergence time, and allowable tracking deviation. Therefore, the design of the specified performance function is crucial, representing a focal point in recent research on PPC. Jia [13] and Veiginis [14] refined the basic performance constraint function based on initial error, allowing for adjustable constraint levels. However, their practicality is limited by the understanding of the initial system states. Such methods are not applicable when the initial error is unknown or the actual situation deviates from the preset initial conditions, such as sudden disturbances or different task instructions. Overly conservative designs can result in slow responses. Bu [15] and Xiao [16] improved control efficiency for nonlinear systems by introducing hyperbolic tangent and exponential functions, creating a performance constraint independent of initial tracking errors. However, the independence from initial conditions may lead to excessive system overshooting in cases of loose initial conditions. In quadrotor applications, Xu et al. [17] integrated PPC with adaptive dynamic programming, enabling finite-time attitude constraint control under external disturbances. Shao et al. [11] established a threshold-triggered expanded state observer based on PPC, enabling adjustments of quadrotor attitude under specific overshoot conditions. Wang et al. [18] integrated PPC with backstepping, accounting for coupling between translational and rotational subsystems, enabling precise position tracking of quadcopter quadrotors with minimal overshoot. However, while the above-mentioned studies quantitatively constrain the targeted indicators, singular problems arise when errors exceed the performance envelope. Therefore, Wang [19] and Hu [20] addressed the issue by introducing terms related to the derivative of the command in the prescribed performance function, allowing constraints to adapt based on changes in task states. However, it is usually difficult to obtain the command form expressed in time. Particularly for abrupt commands like step signals or error steps induced by sudden disturbances, the method becomes inapplicable. But, these situations are extremely common, even typical, in quadrotor aggressive flights. Therefore, to leverage the advantage of PPC in constraining specific performance and improve the control effect of a quadrotor aggressive flight, it is necessary to endow its performance function with the ability to adapt to the out-of-bounds caused by sudden errors.
Drawing inspiration from the above, in this paper, we propose a dynamic event-triggered prescribed performance control (DETPPC) method and design controllers for the translational and rotational subsystem of the quadrotor. Subsequently, the boundedness and stability of the controller based on Lyapunov functions were proved. Moreover, attitude control and aggressive trajectory tracking simulations were conducted under command mutations and external disturbances to demonstrate the effectiveness of the proposed method. The main contributions of this paper are summarized as follows.
  • A dynamic event-triggered performance function boundary adaptive adjustment mechanism is proposed. By assessing the error variation direction and its distance from the performance function boundary, a dynamic compensation factor related to input commands and current errors is established for the performance function, preventing the system divergence caused by errors exceeding the performance function.
  • A robust control framework based on DETPPC is established. Based on this, we designed the controller for the quadrotor translational and rotational subsystems by separately constructing unconstrained error dynamics and performance functions incorporating a dynamic event-triggered compensation factor. Furthermore, stability analysis was conducted based on Lyapunov functions.
  • Attitude control and aggressive trajectory tracking simulations were conducted under command mutations. The results demonstrate that, compared to traditional PPC and cascade PID methods, the proposed DETPPC controller avoids the system divergence caused by errors exceeding the performance function and achieves a faster response and convergence, as well as reduced oscillations.
Based on the above, our method not only preserves the preset constraints on the quadrotor transient performance but also innovatively enhances flight robustness and algorithm practicality by addressing the divergence caused by errors exceeding the performance function during command changes. The following is the content arrangement of this paper. Section 2 gives the problem formulation and preliminaries, including the quadrotor model and the control objective. The design of the dynamic event-triggered prescribed performance function is presented in Section 3. Section 4 and Section 5 construct controllers and stability analysis, respectively. Section 6 presents the simulation and discussion, while the conclusion is given in Section 7.

2. Problem Formulation and Preliminaries

This section primarily discusses the kinematic and dynamic models of quadcopters, serving as the foundation for subsequent research. This paper focuses on a quadcopter with an “X” configuration, as illustrated in Figure 1. Numbers 1–4 represent the propeller number.
Take F b : = [ b x b y b z ] to represent the body-fixed reference frame whose direction is indicated by a yellow arrow in Figure 1. T is the total thrust produced by the propellers and M i , i = x , y , z is the resultant external torque acting on the three axes of the quadrotor. Define the earth reference frame as F e : = [ e x e y e z ] . Figure 2 displays the configuration. l and d represent half the distance between the front and rear rotors and between the center and right rotors, respectively.
Assuming the quadrotor is a rigid body, its kinematics and dynamics model [21] can be expressed as follows:
P ˙ = V V ˙ = g e 3 + 1 m R ( Θ ) F Θ ˙ = Ψ ( Θ ) ω ω ˙ = ω × J ω + M J
where m is the mass of the quadrotor, g is the gravitational acceleration, and e 3 = [ 0 , 0 , 1 ] T is a unitary vector in F e . P = Δ [ x , y , z ] T and V = Δ [ V x , V y , V z ] T are the quadrotor’s position and velocity vector in F e , respectively. Θ = Δ [ ϕ , θ , ψ ] T is the Euler angles to describe the flight attitude of the quadrotor, including the roll, pitch, and yaw angles. ω = Δ [ p , q , r ] T denotes the angular velocity in the F b . R ( Θ ) , as shown in (2), is the rotation matrix from F b to F e . Ψ ( Θ ) is a specific matrix and is defined as (3).
R ( Θ ) = C θ C ψ S ϕ S θ C ψ C ϕ S ψ C ϕ S θ C ψ + S ϕ S ψ C θ S ψ S ϕ S θ S ψ C ϕ C ψ C ϕ S θ S ψ S ϕ C ψ S θ S ϕ C θ C ϕ C θ
where C ( · ) and S ( · ) are abbreviations of cos ( · ) and sin ( · ) , respectively.
Ψ ( Θ ) is a specific matrix and is defined as follows:
Ψ ( Θ ) = 1 S ϕ C θ C ϕ S θ 0 C ϕ S ϕ 0 S ϕ / C θ C ϕ / C θ
F and M in (1) are the total aerodynamic force and moment on the quadrotor, which can be expressed as follows:
F = F p + F e x t M = M p + M e x t
where the subscript p represents the aerodynamics generated by the propellers and e x t represents the external disturbance. The composition of F p and M p is as follows:
F p = F b [ 0 , 0 , T ] T = F b · 0 , 0 , i = 1 4 T p i T M p = F b M x , M y , M z T = F b · diag [ l , d , 1 ] T p 1 + T p 2 + T p 3 T p 4 T p 1 T p 2 + T p 3 T p 4 Q p 1 Q p 2 + Q p 3 + Q p 4
where T p n and Q p n ( n = 1 , 2 , 3 , 4 ) are the thrust and torque of each propeller, respectively.
The goal of quadrotor flight control is to adjust propeller rotation speeds to change T and M p , then alter the quadrotor’s angular velocity, attitude, velocity, and other states to track the desired commands. To fulfil the afore-mentioned objectives, we adopt the following assumptions to facilitate subsequent research and analysis.
Assumption A1. 
The desired state commands are continuous and bounded and the states P , V , Θ , and ω in (1) can be obtained through sensor measurements.
Assumption A2. 
F e x t and M e x t , generated by the external disturbance, are bounded, Lipschitz continuous, and can be regarded as slowly changing.
Furthermore, to implement prescribed performance control on the desired quadrotor states, establishing their relationship with the corresponding control signals is essential. The following section will elaborate on this.

3. Dynamic Event-Triggered Prescribed Performance Function

To achieve prescribed control over the targeted performance, Bechlioulis et al. [22] proposed the PPC method, which aims to drive the tracking error to converge to a designed envelope. Building on this foundation, various PPC-based control methods have been developed. However, existing studies still cannot solve the system divergence due to the errors exceeding the bounds caused by variations in task states after the system reaches a steady state. Therefore, we developed a novel performance function based on dynamic event triggering to facilitate envelope adjustments.

3.1. Prescribed Performance Function Design

For a general control system represented as (6), x ( t ) is the system state, y ( t ) is the system output, and u ( t ) is the control input.
x ˙ ( t ) = f ( x ( t ) ) + g ( x ( t ) ) u ( t ) y ( t ) = x ( t )
The tracking error e ( t ) can be expressed as y ( t ) x r e f ( t ) , where the subscript “ r e f ” represents the reference command. To design prescribed performance control, first define the error manifold as follows:
ε ( t ) = α e ( t ) + e ˙ ( t )
where α > 0 is an adjustable parameter.
To alleviate dependence on the system’s initial state, unlike previous studies, the performance function in this paper is designed in the following form.
β ¯ ( t ) = 1 e a t β + sat e ( t ) Ξ β ̲ ( t ) = β ¯ ( t )
where β represents the maximum allowable boundary for steady-state errors, while a reflects the desired error convergence time. There are designed parameters, β > 0 and α > a 0 . sat e represents the dynamic triggering term, and Ξ is the compensation factor associated with the command-related performance function. Their definition and composition are provided in Section 3.2.

3.2. Dynamic Event-Triggered Mechanism

In this paper, dynamic event triggering serves to detect whether errors exceed bounds, dynamically adjusting the performance function envelope to maintain effective prescribed performance control. The dynamic event-triggered condition is designed as follows:
sat ( e ) = 1 , e ( t ) > k e β & d e ( 1 k e ) β 0 , else
where k e ( 0 , 1 ) is the designed parameter which represents the permissible safety range of error during a steady state. d e = min ( e ( t ) β ̲ ( t ) , e ( t ) + β ¯ ( t ) ) represent the minimum value between the error and the boundary of the performance function. Equation (9) indicates that the triggering of performance function compensation requires the simultaneous satisfaction of two conditions: the absolute error value exceeds the predetermined safety range, and the minimum distance between the error and the upper and lower bounds of the performance function is outside the range.
The dynamic term Ξ in the performance function is designed as follows:
Ξ = x r e f ( t ) e a t + e ( t ) y + ( 1 k e ) β e a t
Remark 1. 
In (8), β ¯ ( t ) , β ̲ ( t ) , β ¯ ˙ ( t ) , and β ̲ ˙ ( t ) are bounded. ε 0 + η constrains the maximum overshoot of the tracking error manifold. β , β represents the maximum range allowed for s ( t ) at a steady state.
Remark 2. 
In (10), the first part introduces x r e f ( t ) to extend the boundaries of β ( t ) when the command changes abruptly. The second part is aimed at enlarging the boundaries to ensure that the deviation between e ( t ) and β ( t ) is not less than ( 1 k e ) β when e ( t ) is uncontrollable. Both are aimed at preventing error overflow and consequent system divergence.
To ensure that ε ( t ) lies within β ¯ ( t ) , β ̲ ( t ) , we define the transformed error before designing the control law. First, the normalized error λ ( t ) is defined as follows:
λ ( t ) = 2 ε ( t ) β ¯ ( t ) + β ̲ ( t ) / β ¯ ( t ) β ̲ ( t )
Then, the transformed error can be designed as follows:
E ( t ) = λ ( t ) 1 λ 2 ( t )
Taking the derivative of E ( t ) with respect to time yields, we have,
E ˙ ( t ) = Λ ε ˙ ( t ) + Φ
where
Λ = 2 1 + λ 2 ( t ) 1 λ 2 ( t ) 2 β ¯ ( t ) β ̲ ( t )
and
Φ = λ ( t ) β ¯ ˙ ( t ) β ̲ ˙ ( t ) β ¯ ˙ ( t ) + β ̲ ˙ ( t ) 2
Theorem 1. 
Supposing that there exists a constant E m > 0 satisfying E ( t ) E m for t 0 , then it follows that λ ( t ) < 1 .
Proof. 
According to β ̲ ( 0 ) = x r e f + β and β ¯ ( 0 ) = x r e f β , then λ ( 0 ) < 1 .
  • For 0 E ( t ) E m for t 0 , there are two cases: (a) 0 λ ( t ) < 1 and (b) λ ( t ) < 1 . The case (b) contradicts that λ ( 0 ) < 1 .
  • For E m E ( t ) < 0 , there have two cases: (a) 1 < λ ( 0 ) < 0 and (b) λ ( t ) > 1 . Case (b) contradicts λ ( 0 ) < 1 .
Therefore, if E ( t ) E m was valid, then λ ( t ) < 1 and β ̲ ( t ) < ε ( t ) < β ¯ ( t ) would exist for t 0 .
Figure 2 illustrates the schematic of the proposed dynamic event-triggered prescribed performance control method. The gray dashed line represents the steady-state boundary before performance function adjustment.
Remark 3. 
It can be observed that, when the command mutation occurs at t = t * , the system error undergoes large-scale variation, exceeding the original steady-state performance function envelope. However, the performance function can dynamically adjust instantly, ensuring envelope error tracking and controlling it within the initially set constraints.
Remark 4. 
Compared with the traditional PPC methods introduced in study [23], the proposed performance functions (8) with the dynamic event-triggered mechanism (9) not only quantitatively constrain the state of the quadcopter but also prevent the system divergence caused by error overflow, thus enhancing practicality and flight safety.

4. Controller Design

The dynamic event-triggered prescribed performance controller scheme of the quadrotor is shown in Figure 3, and the subsequent sections will provide detailed explanations of the design of each component.

4.1. Translational Subsystem Controller Design

Define the tracking error as ˜ = d , where • represents represents any state variable of the quadrotor in (1). Then, considering Assumptions A1 and A2, the error dynamics of the translation subsystem can be represented as follows:
P ˜ ˙ = P ˙ P ˙ d V ˜ ˙ = g e 3 + F p R ( Θ ) m V ˙ d
Considering Θ ˜ = Θ Θ d = ϕ ˜ , θ ˜ , ψ ˜ T , the transformation matrix R ( Θ ) in (16), composed of trigonometric functions of the attitude angles, as shown in (2), can be obtained using the sum and difference identities of trigonometric functions as follows:
R ( Θ ) = C ψ S ψ 0 S ψ C ψ 0 0 0 1 C ϕ ˜ C θ ˜ S θ S ϕ ˜ C ϕ ˜ S θ ˜ S θ d S ϕ ˜ C ϕ ˜ C θ d S ϕ ˜ C ϕ ˜ S θ ˜ C θ S ϕ ˜ C ϕ ˜ C θ ˜ C ϕ d S θ d S ϕ d C ϕ d C θ d
Considering that F p = 0 , 0 , T T , combined with (17), (16) by can be rewritten as (18).
P ˜ ˙ = P ˙ P ˙ d V ˜ ˙ = g e 3 + R ( ψ ) R ˜ F T m V ˙ d
where
R ˜ = C ϕ ˜ C θ ˜ S θ S ϕ ˜ C ϕ ˜ S θ ˜ S θ d S ϕ ˜ C ϕ ˜ C θ d S ϕ ˜ C ϕ ˜ S θ ˜ C θ S ϕ ˜ C ϕ ˜ C θ ˜
F T = F x F y F z = C ϕ d S θ d T S ϕ d T C ϕ d C θ d T
R ( ψ ) = C ψ S ψ 0 S ψ C ψ 0 0 0 1
With reference to (7), define ϵ T = ϵ x , ϵ y , ϵ z T = α P ˜ + P ˜ ˙ as the error manifold of the translational subsystem, where α = diag α x , α y , α z . The performance function of the translational subsystem in this paper is designed as follows:
β ¯ i ( t ) = 1 e a i t β i , + sat e i ( t ) Ξ i β ̲ i ( t ) = β ¯ i ( t )
where α i > a i > 0 , i = x , y , z .
Combining the definition of a normalized error in this paper, the control law for the translational subsystem can be derived as follows:
u T = R ( ψ ) Λ T κ T E T + g e 3
where κ T = diag ( κ x , κ y , κ z ) is a positive definite control gain matrix, E T = diag ( E x , E y , E z ) and Λ T = diag ( Λ x , Λ y , Λ z ) are expressed by (24) and (25).
E i ( t ) = λ i ( t ) 1 λ i 2 ( t )
Λ i = 2 1 + λ i 2 ( t ) 1 λ i 2 ( t ) 2 ( β ¯ i ( t ) β ̲ i ( t ) )
u T = [ u x , u y , u z ] T in (23) essentially represents the linear acceleration command. Generally, based on force equilibrium, the desired attitude angles and total thrust of a quadrotor can be obtained by applying Newton’s second law to the linear acceleration commands. However, this neglects the influence of air resistance, leading to errors in command derivation. Based on an incremental nonlinear dynamic inversion method, the authors of [24] established mapping between quadcopter linear acceleration and horizontal attitude angles, as shown in (26), demonstrating a significant enhancement in the accuracy of attitude angle commands.
m a d , x / T m a d , y / T = m a x / T m a y / T + C ϕ S ψ S ϕ S θ C ψ C ϕ C θ C ψ C ϕ C ψ S ϕ S θ S ψ C ϕ S θ S ψ ϕ d θ d ϕ θ
where T is the total thrust. Variables with subscript d represent desired commands, while the remaining state variables represent measurements from sensors at the current time. Based on (26), considering the definition of u T , then based on the translational subsystem controller, the desired total thrust, roll, and pitch angle of the quadrotor can be obtained as follows:
T d = m u z ϕ d θ d = ϕ θ + R 1 1 m a x u x / T m a y u y / T
where
R 1 = C ϕ S ψ S ϕ S θ C ψ C ϕ C θ C ψ C ϕ C ψ S ϕ S θ S ψ C ϕ S θ S ψ

4.2. Rotational Subsystem Controller Design

Similar to the translational subsystem controller design, the performance function for the attitude control of the quadrotor is selected as (29).
β ¯ j ( t ) = 1 e a j t β j , + sat e j ( t ) Ξ j β ̲ j ( t ) = β ¯ j ( t ) , j = ϕ , θ , ψ
Define the output of the quadrotor attitude controller as u Θ = p d , q d , r d T ; this means the virtual desired angular velocity, which can be expressed in (30).
u Θ = p d q d r d = Ψ ( Θ ) 1 κ Θ E Θ
where κ Θ = diag ( κ ϕ , κ θ , κ ψ ) is a positive definite control gain matrix and E Θ = diag ( E ϕ , E θ , E ψ ) can be expressed by (31).
E j ( t ) = λ j ( t ) 1 λ j 2 ( t ) , λ j ( t ) = Θ ˜ j ( t ) / β ¯ j ( t )
where j = ϕ , θ , ψ .
For angular velocity control, the performance functions are selected as (32) and the control signal u ω can be deduced as (33).
β ¯ k ( t ) = 1 e a k t β k , + sat e k ( t ) Ξ k β ̲ k ( t ) = β ¯ k ( t ) , k = p , q , r
u ω = u p u q u r = κ ω E ω
where κ ω = diag ( κ p , κ q , κ r ) is a positive definite control gain matrix and E ω = diag ( E p , E q , E r ) can be expressed by (34).
E k ( t ) = λ k ( t ) 1 λ k 2 ( t ) , λ k ( t ) = ω ˜ k ( t ) / β ¯ k ( t )
where k = p , q , r .
Therefore, the desired moment of the quadrotor can be obtained as M d = [ M x , d , M y , d , M z , d ] T = J u ω .

5. Stability Analysis

To facilitate derivation, combining (1) and (6), rewrite the quadrotor angular velocity subsystem as follows:
ω ˙ k = f k ( ω k , u k ) + b k , k = p , q , r
where f k ( ω k , u k ) is an unknown non-affine Lipschitz continuous function; b k is the external disturbance.
According to [18], there exist continuous functions, C ̲ k ( ω k , u k ) , C ¯ k ( ω k , u k ) , c ̲ k ( ω k ) , and c ¯ k ( ω k ) , and unknown constants, ν ̲ k and ν ¯ k , k = p , q , r , which make the following inequalities hold.
f k ω , u k C ̲ k ω , u k u k + c ̲ k ( ω ) , u k ν ̲ k 0 f k ω , u k C ¯ k ω , u k u k + c ¯ k ( ω ) , u k ν ¯ k 0
Theorem 2. 
Suppose that both f k ( ω k , u k ) and b k are bounded. By combining the quadrotor translational and rotational subsystem control outputs, (23), (30), and (33), the following properties can be ensured:
  • For t 0 , β ̲ l ( t ) < e l ( t ) < β ¯ l ( t ) exists with l = x , y , z , ϕ , θ , ψ , p , q , r and e l ( t ) = x ˜ , y ˜ , z ˜ , ϕ ˜ , θ ˜ , ψ ˜ , p ˜ , q ˜ , r ˜ .
  • All the state variables in (1) for quadcopters are bounded.
Proof. 
Based on (36) and the extreme value theory, there exist continuous functions, G ̲ k ( u k ) and G ¯ k ( u k ) , such that the following inequalities hold.
f k ( ω k , u k ) G ̲ k , 0 u k ν k f k ( ω k , u k ) G ¯ k , ν k u k 0
where ν k , k = p , q , r are unknown constants.
Thus, according to (36) and (37), it can be obtained that
f k ω k , u k C ̲ k ω k , u k u k + c ̲ k ( ω k ) , u k ν k f k ω k , u k G ̲ k ( ω k ) + C k , m u k ν k , 0 u k ν k f k ω k , u k G ¯ k ( ω k ) + C k , m u k + ν k , ν k u k 0 f k ω k , u k C ¯ k ω k , u k u k + c ¯ k ( ω k ) , u k ν k
Then, rewrite f k ( ω , u k ) as follows:
f k ( ω , u k ) = C k ( ω , u k ) u k + c k ( ω )
with
C k ω k , u k = γ k , 1 C k ω k , u k , u k ν k γ k , 2 C k , m , 0 u k ν k γ k , 3 C k , m , ν k u k 0 γ k , 4 C ¯ k ω k , u k , u k ν k
c k ( ω k ) = c ̲ k ( ω k ) , u k ν k C k , m ν k G ¯ k ( ω k ) , 0 u k ν k C k , m ν k + G ¯ k ( ω k ) , ν k u k 0 c ¯ c ( ω k ) , u k ν k
where C k , m C k ( ω k , u k ) and c k ( ω k ) max c ̲ k ( ω k ) , C k , m ν k G ̲ k ( ω k ) , C k , m ν k + G ¯ k ( ω k ) , c ¯ k ( ω k ) , factors γ k , n [ 1 , + ) , k = p , q , r .
Then, define the normalized error for all states of the quadrotor as λ ( t ) = [ λ x , λ y , λ z , λ ϕ , λ θ , λ ψ , λ p , λ q , λ r ] T , whose derivative λ ˙ l ( t ) is defined as follows:
λ ˙ l ( t ) = λ ¯ ( t , λ l ( t ) )
where l = x , y , z , ϕ , θ , ψ , p , q , r .
Define λ = ( 1 , 1 ) × ( 1 , 1 ) × × ( 1 , 1 ) 9 × , which constrains that λ ( 0 ) λ , and there exists a maximal solution t max that fulfils λ ( t ) λ for t 0 , t max [18]. Thus, set a constant H ¯ l and H ¯ l max β ̲ l ( 0 ) , β ¯ l ( 0 ) , then it can be obtained that e l ( t ) < H ¯ l .
Step 1. 
Differentiate the error manifold of the translational subsystem as follows:
ε ˙ T ( t ) = R ( ψ ) m R ˜ u T + o T
where u T = u T g e 3 and o T = o x , o y , o z T = P ¨ d + α T P ˙ α T P ˙ d .
Define the Lyapunov function candidate V T = E T T K T E T , where E T = E x , E y , E z T and K T = [ K x , K y , K z ] T are a positive definite diagonal matrix. Take the derivative of V T with respect to time as follows:
V ˙ T = E T T K T E ˙ T + E ˙ T T K T E T = E T T K T Λ T R ( ψ ) m R ˜ u T + o T + Φ T + u T T R ˜ T R T ( ψ ) m + o T T + Φ T T Λ T K T E T
where Λ = diag ( Λ x , Λ y , Λ z ) and Φ = [ Φ x , Φ y , Φ z ] T . Choose K T = κ T and submit (23) into (43); we thus have
V ˙ T = E T T κ T Λ T R ( ψ ) m R ˜ R T ( ψ ) Λ T κ T E T + o T + Φ T     + E T T κ T Λ T R ( ψ ) R ˜ T R T ( ψ ) m + o T T + Φ T T Λ T κ T E T = 2 m E T T κ T Λ T R ( ψ ) R ˜ s R T ( ψ ) Λ T κ T E T     + E T T κ T Λ T o T + Φ T + o T T + Φ T T Λ T κ T E T
where R ˜ s = R ˜ + R ˜ T / 2 .
Consider the following boundedness condition:
  • β ¯ i ( t ) , β ¯ ˙ i ( t ) , β ̲ i ( t ) , and β ̲ ˙ i ( t ) , i = x , y , z are bounded;
  • P d , V d , and V ˙ d are bounded;
  • According to extreme value theory, o i and Φ i for t 0 , t max , there is a constant ς T > 0 that satisfies o T + Φ T ς T for t 0 , t max .
Based on Young’s inequality, together with R ( ψ ) R T ( ψ ) = I , it can be deduced that
E T T κ T Λ T o T + Φ T m ς T 2 4 μ o min R ˜ s + μ E T T κ T Λ T R ( ψ ) o min R ˜ s R T ( ψ ) Λ T κ T E T m
where 0 < μ < 1 . Because R ˜ s is positive definite, then we have
2 m μ E T T κ T Λ T R ( ψ ) o min R ˜ s I R ˜ s R T ( ψ ) Λ T κ T E T 0
Use O T to respresent 2 m ( 1 μ ) κ T Λ T R ( ψ ) R ˜ s R T ( ψ ) Λ T κ T , then take (46) and (47) into (45); it can be deduced that
V ˙ T E T T O T E T + m ς T 2 2 μ o min R ˜ s
According to the boundary theory of nonlinear systems based on definiteness [25], it can be obtained that E T E ¯ T max E T ( 0 ) , T for t 0 , t max , where
T ς T o max K T m 2 o min R ˜ s o min O T o min K T μ
Therefore, E i and u T , i , i = x , y , z are bounded and 1 < λ ̲ i λ i ( t ) λ ¯ i < 1 for t 0 , t max .
Step 2. 
Differentiate Θ ˜ with respect to time as follows:
Θ ˜ ˙ = Ψ ( Θ ) ( ω d + ω ˜ ) Θ ˜ ˙ d = Ψ ( Θ ) ω d + o Θ
where o R = Ψ ( Θ ) Θ ˜ Θ ˙ d .
Define the Lyapunov function candidate V Θ = E Θ T K Θ E Θ = V ϕ , V θ , V ψ T , where E Θ = E ϕ , E θ , E ψ T and K Θ = [ K ϕ , K θ , K ψ ] T represent a positive definite diagonal matrix. Take the derivative of V Θ as follows:
V ˙ Θ = E Θ T K Θ E ˙ Θ + E ˙ Θ T K Θ E Θ = E Θ T K Θ Λ Θ Ψ ( Θ ) ω d + o Θ + Φ Θ + ω d T Ψ T ( Θ ) + o Θ T + Φ Θ T Λ Θ K Θ E Θ
where Λ Θ = diag ( Λ ϕ , Λ θ , Λ ψ ) , Φ Θ = [ Φ ϕ , Φ θ , Φ ψ ] T . Choosing K Θ = κ Θ and taking (30) into (51), we have
V ˙ Θ = 2 E Θ T κ Θ Λ Θ κ Θ E Θ + E Θ T κ Θ Λ Θ o Θ + Φ Θ + o Θ T + Φ Θ T Λ Θ κ Θ E Θ
Consider the following boundedness conditions:
  • β j ( t ) and β j ˙ ( t ) , j = ϕ , θ , ψ are bounded;
  • Θ  and  Θ ˙ are bounded by the former step;
  • o j and Φ j are bounded for t 0 , t max .
Because κ Θ and Λ Θ are both positive diagonal matrixes, it can be deduced as follows:
V ˙ j 2 Λ j κ j 2 E j 2 + 2 ς j Λ j κ j E j
Thus, we can conclude that E j E ¯ j max E j ( 0 ) , ¯ j for t 0 , t max with j ς j / κ j . Thus, Θ j , j = ϕ , θ , ψ and ω d , k , k = p , q , r are bounded and 1 λ ̲ j λ j ( t ) λ ¯ j < 1 for t 0 , t max .
Step 3. 
Differentiate ω ˜ k , k = p , q , r respect to time as follows:
ω ˜ ˙ k = f k ( ω k , u k ) ω ˙ d , k = C k ( ω k , u k ) u k + o k
where o k = c k ( ω ) ω ˙ d , k .
Define the Lyapunov function candidate V k = E k 2 / 2 , then its derivative with respect to time can be represented as follows:
V ˙ k = E k E ˙ k = E k Λ k ( C k ( ω k , u k ) u k + o k + Φ k )
Taking (33) into (55), we have
V ˙ k = E k Λ k ( C k ( ω k , u k ) κ k + o k + Φ k )
Consider the following boundedness conditions:
  • β k ( t ) and β ˙ k ( t ) , k = p , q , r are bounded;
  • ω k and ω k ˙ are bounded by the former step;
  • Based on the extreme value theory, o k and Φ k are bounded for t 0 , t max and there exists a constant ς k > 0 that satisfies o k + Φ + k ς k
According to the above conditions, it can be obtained as follows:
V ˙ k C k ( ω k , u k ) Λ k κ k E k 2 + Λ k ς k E k
After that, it can be concluded that E k E ¯ k max E k ( 0 ) , ¯ k for t 0 , t max with k ς k / ( C k , m κ k ) . Therefore, we can obtain that E k and u k , k = p , q , r are bounded and 1 λ ̲ k λ k ( t ) λ ¯ k < 1 for t 0 , t max .
Step 4. 
According to the setting that λ ( t ) λ   for t 0 , t max , where
λ   = i { x , y , z } λ ̲ i , λ ¯ i × j { ϕ , θ , ψ } λ ̲ j , λ ¯ j × k { p , q , r } λ ̲ k , λ ¯ k
Thus, it can be assumed that t max < and λ   λ , and it can be demonstrated that there exists a time instant t t [ 0 , t max ) satisfying λ ( t t ) λ   [25]. It is obvious that such a deduction is contradictory to λ ( t ) λ   .
Hence, it can be concluded that all states of the closed-loop quadrotor based on the proposed dynamic event-triggered PPC remain bounded and t max = . Furthermore, we can obtain the following conclusion:
β ̲ l ( t ) < e l ( t ) < β ¯ l ( t ) , t 0
Henceforth, all proofs are complete.

6. Simulation and Discussion

In order to verify the effectiveness of the proposed method, we provide quadrotor aggressive flight simulations. The parameters of the quadrotor model and DETPPC controllers designed are as shown in Table 1. At the same time, the measurement noise of the sensor is taken into account.
The performance functions of each subsystem are, respectively, designed as (22), (29), and (32), the parameters of which are displayed in Table 1.
This paper selects two controllers for comparison: (1) a traditional prescribed performance controller (TPPC) without dynamic event-triggering, with controller parameters consistent with DETPPC; (2) a Cascade PID controller (P-PID), with parameters set according to Wang [26].
Example 1. 
To verify the improvement in the transient performance of the proposed method for quadrotors facing abrupt commands, this example conducts a simulation of an attitude control under abrupt mutations. The reference attitude commands ( Θ d ), as defined in (60), are employed for attitude control simulation. Θ d in (60) denotes that the roll and pitch attitude angles of the quadrotor were initially subjected to sinusoidal commands. However, at the 3rd and 9th seconds, due to specific flight requirements, their values experienced abrupt changes and were maintained for 3 s each. The initial attitude and angular velocity of the quadrotor are set to Θ 0 = [ 0 , 0 , 0 ] T and ω 0 = [ 0 , 0 , 0 ] T .
Θ d ( deg ) = ϕ d = 20 sin π 10 t θ d = 30 sin π 10 t     t < 3 s & 6 s < t < 9 s & t > 12 s ϕ d = 20 θ d = 30     3 s < t < 6 s ϕ d = 10 θ d = 15     9 s < t < 12 s
Figure 4 and Figure 5 present the tracking results. Due to the absence of performance function compensation based on dynamic event triggering, the TPPC controller diverges when the command abruptly changes at the third second. In comparison with the P-PID method, the DETPPC controller in this study demonstrates faster response and convergence rates, as well as smaller tracking errors. Figure 6 illustrates a comparison of the convergence of Θ ˜ . The purple and green dashed lines represent the proposed DETPPC performance function. The TPPC method encounters an error overshoot at the third second, which is the cause of control divergence. The proposed DETPPC controller demonstrates faster and more stable error convergence than the P-PID controller. In Figure 7, DETPPC’s control gain dynamically adjusts based on Θ ˜ variations. This accounts for DETPPC’s faster convergence and response speed compared to traditional cascade PID methods when handling abrupt commands.
Figure 8 and Figure 9, respectively, depict ω ˜ and the variation of DETPPC control gain with the angular rate tracking errors. Besides the divergence seen in the TPPC method during command mutations, it is evident that, for pitch angular velocity control, the DETPPC controller also shows faster and more stable convergence. Moreover, the tracking error based on the P-PID method occasionally surpasses the expected performance boundaries. For roll angular velocity, DETPPC and P-PID controllers perform similarly, likely due to the small reference command for roll angle.
To sum up, the proposed DETPPC method resolves the divergence problem faced by traditional PPC when confronted with command mutations, while also achieving superior transient performance through tailored constraint design for specified state variables.
Example 2. 
To further validate the afore-mentioned conclusions, this example provides a simulation of agile trajectory tracking for a quadrotor with abrupt command. As errors exceeding performance function boundaries lead to system divergence, the traditional PPC method is not included as a comparison in this simulation. The reference trajectory in this example is a three-dimensional helix, as shown in (61). External disturbances are encountered at the tenth second, persisting for two seconds. Additionally, trajectory command changes occur at the twentieth and thirtieth seconds, respectively. The initial position of the quadcopter is set to P 0 = [ 0 , 0 , 0 ] T .
P d ( m ) = x d = 40 sin π 10 t y d = 40 cos π 10 t 40 z d = t t < 20 s x d = 2 t 40 y = 0 z = 20 20 s < t < 30 s x d = 40 sin π 10 t + 20 y d = 40 cos π 10 t + 40 z = t + 50 30 s < t < 40 s
Figure 10 and Figure 11 present the results of trajectory tracking and position control in three directions, respectively. It can be observed that the DETPPC method enables more accurate trajectory tracking, which is particularly evident during the command mutation at 20 s, where the proposed DETPPC method exhibits reduced overshoot and faster position convergence. Take the root-mean-square (RMS, e p 1 ) values of the tracking errors, as shown in Table 2, to quantitatively evaluate the controller effectiveness. The results demonstrate that the trajectory tracking RMS of the proposed method is only 31% of that of the traditional approach.
The convergence of P ˜ for the quadrotor is illustrated in Figure 12. By imposing performance constraints, DETPPC achieves faster convergence of tracking errors with smaller overshoots through automatic adjustment of control gains, as depicted in Figure 13. In comparison, the position control errors of the P-PID method notably deviate from the desired criteria, especially in the horizontal direction. The position tracking error of the proposed method is only 50% that of the P-PID controller in the horizontal direction, and 75% in the vertical direction. Similarly, the velocity control of the quadcopter in three axes shows the same trend. As depicted in Figure 14, DETPPC facilitates quicker convergence of velocity control, particularly with reduced fluctuations during command changes.
Figure 15, Figure 16 and Figure 17 demonstrate the control comparison of attitude angles and angular rates. For the desired trajectory in (61), the quadrotor’s attitude angles reached 40 degrees. Furthermore, the proposed DETPPC method reaches peak attitude angles nearly 2 s faster than the P-PID controller. Figure 15 demonstrates that, with DERPPC, the quadrotor achieves more agile attitude control. Additionally, in Figure 16, during the 10–12 s period of external disturbance, the proposed method demonstrates rapid attitude convergence, while the P-PID controller exhibits larger and longer fluctuations, especially in roll and pitch angles. ω ˜ in Figure 17 also reflects the same conclusion.

7. Conclusions

Originating from the demands of quadrotor aggressive flight for maneuverability and transient performance, this paper aims to address the control divergence caused by command mutations in existing PPC methods. Through designing a dynamic event-triggered mechanism, and using the direction of error variation and its distance from the envelope of the performance function as triggering criteria, the performance function can be autonomously compensated. Based on this, the performance function was redesigned, and robust prescribed performance controllers for the quadrotor’s translational and rotational subsystems were, respectively, designed. Moreover, a stability analysis of the controller was conducted based on Lyapunov functions. Finally, simulations were conducted to evaluate attitude control and aggressive trajectory tracking under command mutations and external disturbances. The results indicate that the proposed method resolves the divergence of traditional PPC controllers under abrupt command changes, enhancing the algorithm’s robustness. Additionally, compared to the cascaded PID method, the proposed DETPPC controller exhibits faster response and convergence speeds, along with reduced oscillations, in trajectory tracking simulations under sudden command changes. Ultimately, it achieves an RMS trajectory tracking error that is only 31% of the control group. We believe that this work can provide a reference for practical applications, while also indicating areas that require further exploration. It is worth noting that our proposed method has shown that control gains vary with error convergence. Hence, energy consumption during application needs attention, even potentially control saturation. This provides us with further ideas for our future research and for enhancing the practicality of the algorithm.

Author Contributions

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

Funding

This research was funded by the China Postdoctoral Science Foundation (2023M731913) and Civilian Aircraft Research (MJG5-1N21).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fan, B.; Li, Y.; Zhang, R.; Fu, Q. Review on the technological development and application of UAV systems. Chin. J. Electron. 2020, 29, 199–207. [Google Scholar] [CrossRef]
  2. Zhen, Z.; Jiang, J.; Sun, S.; Wang, B. Cooperative Control and Decision of UAV Swarm Operations; National Defense Industry Press: Beijing, China, 2022. [Google Scholar]
  3. Liu, Z.; Wang, X.; Shen, L.; Zhao, S.; Cong, Y.; Li, J.; Yin, D.; Jia, S.; Xiang, X. Mission-oriented miniature fixed-wing UAV swarms: A multilayered and distributed architecture. IEEE Trans. Syst. Man Cybern. Syst. 2020, 52, 1588–1602. [Google Scholar] [CrossRef]
  4. Bry, A.; Richter, C.; Bachrach, A.; Roy, N. Aggressive flight of fixed-wing and quadrotor aircraft in dense indoor environments. Int. J. Robot. Res. 2015, 34, 969–1002. [Google Scholar] [CrossRef]
  5. Jia, J.; Guo, K.; Yu, X.; Guo, L.; Xie, L. Agile flight control under multiple disturbances for quadrotor: Algorithms and evaluation. IEEE Trans. Aerosp. Electron. Syst. 2022, 58, 3049–3062. [Google Scholar] [CrossRef]
  6. Belmouhoub, A.; Medjmadj, S.; Bouzid, Y.; Derrouaoui, S.; Guiatni, M. Enhanced backstepping control for an unconventional quadrotor under external disturbances. Aeronaut. J. 2023, 127, 627–650. [Google Scholar] [CrossRef]
  7. Labbadi, M.; Boudaraia, K.; Elakkary, A.; Djemai, M.; Cherkaoui, M. A continuous nonlinear sliding mode control with fractional operators for quadrotor UAV systems in the presence of disturbances. J. Aerosp. Eng. 2022, 35, 04021122. [Google Scholar] [CrossRef]
  8. 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]
  9. Jabeur, C.B.; Seddik, H. Optimized neural networks-PID controller with wind rejection strategy for a Quad-Rotor. J. Robot. Control JRC 2022, 3, 62–72. [Google Scholar] [CrossRef]
  10. 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]
  11. Shao, X.; Xu, L.; Zhang, W. Quantized control capable of appointed-time performances for quadrotor attitude tracking: Experimental validation. IEEE Trans. Ind. Electron. 2021, 69, 5100–5110. [Google Scholar] [CrossRef]
  12. Cui, G.; Yang, W.; Yu, J.; Li, Z.; Tao, C. Fixed-time prescribed performance adaptive trajectory tracking control for a QUAV. IEEE Trans. Circuits Syst. II Express Briefs 2021, 69, 494–498. [Google Scholar] [CrossRef]
  13. Jia, F.; Wang, X.; Zhou, X. Robust adaptive prescribed performance control for a class of nonlinear pure-feedback systems. Int. J. Robust Nonlinear Control 2019, 29, 3971–3987. [Google Scholar] [CrossRef]
  14. Verginis, C.K.; Bechlioulis, C.P.; Dimarogonas, D.V.; Kyriakopoulos, K.J. Robust distributed control protocols for large vehicular platoons with prescribed transient and steady-state performance. IEEE Trans. Control Syst. Technol. 2017, 26, 299–304. [Google Scholar] [CrossRef]
  15. Bu, X.; Wu, X.; Zhu, F.; Huang, J.; Ma, Z.; Zhang, R. Novel prescribed performance neural control of a flexible air-breathing hypersonic vehicle with unknown initial errors. ISA Trans. 2015, 59, 149–159. [Google Scholar] [CrossRef] [PubMed]
  16. Bu, X.; Xiao, Y. Prescribed performance-based low-computational cost fuzzy control of a hypersonic vehicle using non-affine models. Adv. Mech. Eng. 2018, 10, 1687814018757261. [Google Scholar] [CrossRef]
  17. Xu, J.; Wang, L.; Liu, Y.; Xue, H. Finite-time prescribed performance optimal attitude control for quadrotor UAV. Appl. Math. Model. 2023, 120, 752–768. [Google Scholar] [CrossRef]
  18. Wang, Y.; Hu, J. Robust control for a quadrotor aircraft with small overshoot and high-precision position tracking performance. J. Frankl. Inst. 2020, 357, 13386–13409. [Google Scholar] [CrossRef]
  19. Wang, Y.; Hu, J. Improved prescribed performance control for air-breathing hypersonic vehicles with unknown deadzone input nonlinearity. ISA Trans. 2018, 79, 95–107. [Google Scholar] [CrossRef]
  20. Wang, Y.; Hu, J.; Li, J.; Liu, B. Improved prescribed performance control for nonaffine pure-feedback systems with input saturation. Int. J. Robust Nonlinear Control 2019, 29, 1769–1788. [Google Scholar] [CrossRef]
  21. Padfield, G.D. Helicopter Flight Dynamics: Including a Treatment of Tiltrotor Aircraft; John Wiley & Sons: Hoboken, NJ, USA, 2018. [Google Scholar]
  22. Bechlioulis, C.P.; Rovithakis, G.A. A low-complexity global approximation-free control scheme with prescribed performance for unknown pure feedback systems. Automatica 2014, 50, 1217–1226. [Google Scholar] [CrossRef]
  23. Bu, X. Prescribed performance control approaches, applications and challenges: A comprehensive survey. Asian J. Control 2023, 25, 241–261. [Google Scholar] [CrossRef]
  24. Pfeifle, O.; Fichter, W. Cascaded incremental nonlinear dynamic inversion for three-dimensional spline-tracking with wind compensation. J. Guid. Control Dyn. 2021, 44, 1559–1571. [Google Scholar] [CrossRef]
  25. Sastry, S. Nonlinear Systems: Analysis, Stability, and Control; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013; Volume 10. [Google Scholar]
  26. Wang, L. PID Control System Design and Automatic Tuning Using MATLAB/Simulink; John Wiley & Sons: Hoboken, NJ, USA, 2020. [Google Scholar]
Figure 1. Aerodynamics and configuration of the quadrotor.
Figure 1. Aerodynamics and configuration of the quadrotor.
Aerospace 11 00301 g001
Figure 2. A schematic of the proposed dynamic event-triggered prescribed performance control.
Figure 2. A schematic of the proposed dynamic event-triggered prescribed performance control.
Aerospace 11 00301 g002
Figure 3. Dynamic event-triggered prescribed performance control scheme of the quadrotor.
Figure 3. Dynamic event-triggered prescribed performance control scheme of the quadrotor.
Aerospace 11 00301 g003
Figure 4. The tracking results of ϕ angle.
Figure 4. The tracking results of ϕ angle.
Aerospace 11 00301 g004
Figure 5. The tracking results of θ angle.
Figure 5. The tracking results of θ angle.
Aerospace 11 00301 g005
Figure 6. Comparision of the tracking error of attitude angle.
Figure 6. Comparision of the tracking error of attitude angle.
Aerospace 11 00301 g006
Figure 7. The control gain of DETPPC varies with the tracking error.
Figure 7. The control gain of DETPPC varies with the tracking error.
Aerospace 11 00301 g007
Figure 8. Comparision of the tracking error of angular velocity.
Figure 8. Comparision of the tracking error of angular velocity.
Aerospace 11 00301 g008
Figure 9. The control gain of DETPPC varies with the angular velocity tracking error.
Figure 9. The control gain of DETPPC varies with the angular velocity tracking error.
Aerospace 11 00301 g009
Figure 10. The trajectory tracking results in comparison.
Figure 10. The trajectory tracking results in comparison.
Aerospace 11 00301 g010
Figure 11. The results of position control in three directions.
Figure 11. The results of position control in three directions.
Aerospace 11 00301 g011
Figure 12. The convergence of position tracking errors.
Figure 12. The convergence of position tracking errors.
Aerospace 11 00301 g012
Figure 13. The variation of position control gains with tracking errors.
Figure 13. The variation of position control gains with tracking errors.
Aerospace 11 00301 g013
Figure 14. The convergence of velocity control errors.
Figure 14. The convergence of velocity control errors.
Aerospace 11 00301 g014
Figure 15. The response of the quadrotor’s attitude angles.
Figure 15. The response of the quadrotor’s attitude angles.
Aerospace 11 00301 g015
Figure 16. The convergence of attitude angle control errors.
Figure 16. The convergence of attitude angle control errors.
Aerospace 11 00301 g016
Figure 17. The convergence of angular rate control errors.
Figure 17. The convergence of angular rate control errors.
Aerospace 11 00301 g017
Table 1. The structural parameters and DETPPC controller parameters of the quadcopter.
Table 1. The structural parameters and DETPPC controller parameters of the quadcopter.
ParameterValue
m6.95 kg
l0.65 m
d0.66 m
g9.81 m / s 2
J diag ( 0.256 , 0.266 , 0.437 ) kg m 2
α T diag ( 1.2 , 1.2 , 1.2 )
α Θ diag ( 0.008 , 0.008 , 0.015 )
α ω diag ( 0.015 , 0.015 , 0.032 )
κ T diag ( 0.4 , 0.4 , 0.4 )
κ Θ diag ( 2.8 , 2.8 , 2.8 )
κ ω diag ( 0.015 , 0.015 , 0.015 )
k e , l 0.6 , l = x , y , z , ϕ , θ , ψ , p , q , r
a i 1 , i = x , y , z
a j 4 , j = ϕ , θ , ψ
a k 8 , k = p , q , r
β i , 0.2 , i = x , y , z
β j , 0.05 , j = ϕ , θ , ψ
β k , 0.1 , k = p , q , r
Table 2. RMS of the trajectory tracking errors.
Table 2. RMS of the trajectory tracking errors.
IndexDETPPCP-PID
RMS (m)1.27154.1206
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

Wu, Z.; Ye, J.; Song, T. Dynamic Event-Triggered Prescribed Performance Robust Control for Aggressive Quadrotor Flight. Aerospace 2024, 11, 301. https://doi.org/10.3390/aerospace11040301

AMA Style

Wu Z, Ye J, Song T. Dynamic Event-Triggered Prescribed Performance Robust Control for Aggressive Quadrotor Flight. Aerospace. 2024; 11(4):301. https://doi.org/10.3390/aerospace11040301

Chicago/Turabian Style

Wu, Zeliang, Jianchuan Ye, and Tao Song. 2024. "Dynamic Event-Triggered Prescribed Performance Robust Control for Aggressive Quadrotor Flight" Aerospace 11, no. 4: 301. https://doi.org/10.3390/aerospace11040301

APA Style

Wu, Z., Ye, J., & Song, T. (2024). Dynamic Event-Triggered Prescribed Performance Robust Control for Aggressive Quadrotor Flight. Aerospace, 11(4), 301. https://doi.org/10.3390/aerospace11040301

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