Next Article in Journal
A High-Reliability Photoelectric Detection System for Mars Sample Return’s Orbiting Sample
Next Article in Special Issue
Disturbance Observer-Based Backstepping Terminal Sliding Mode Aeroelastic Control of Airfoils
Previous Article in Journal
Large Debris Removal: Using Features of Attitude Motion for Load Factor Regulation during Re-Entry
Previous Article in Special Issue
Adaptive Incremental Nonlinear Dynamic Inversion Control for Aerial Manipulators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design and Experimental Validation of an Adaptive Multi-Layer Neural Network Observer-Based Fast Terminal Sliding Mode Control for Quadrotor System

by
Zainab Akhtar
1,
Syed Abbas Zilqurnain Naqvi
1,
Yasir Ali Khan
2,
Mirza Tariq Hamayun
3 and
Salman Ijaz
2,*
1
Mechatronics and Control Engineering Department, University of Engineering and Technology, Lahore 54000, Pakistan
2
Member of Control System Laboratory, University of Nottingham, Ningbo 315104, China
3
Department of Electrical and Computer Engineering, COMSATS University Islamabad (CUI), Lahore 54000, Pakistan
*
Author to whom correspondence should be addressed.
Aerospace 2024, 11(10), 788; https://doi.org/10.3390/aerospace11100788
Submission received: 19 July 2024 / Revised: 12 September 2024 / Accepted: 13 September 2024 / Published: 24 September 2024
(This article belongs to the Special Issue Challenges and Innovations in Aircraft Flight Control)

Abstract

:
This study considers the numerical design and practical implementation of a new multi-layer neural network observer-based control design technique for unmanned aerial vehicles systems. Initially, an adaptive multi-layer neural network-based Luenberger observer is designed for state estimation by employing a modified back-propagation algorithm. The proposed observer’s adaptive nature aids in mitigating the impact of noise, disturbance, and parameter variations, which are usually not considered by conventional observers. Based on the observed states, a nonlinear dynamic inversion-based fast terminal sliding mode controller is designed to attain the desired attitude and position tracking control. This is done by employing a two-loop control structure. Numerical simulations are conducted to demonstrate the effectiveness of the proposed scheme in the presence of disturbance, parameter uncertainty, and noise. The numerical results are compared with current approaches, demonstrating the superiority of the proposed method. In order to assess the practical effectiveness of the proposed method, hardware-in-loop simulations are conducted by utilizing a Pixhawk 6X flight controller that interfaces with the mission planner software. Finally, experiments are conducted on a real F450 quadrotor in a secured laboratory environment, demonstrating stability and good tracking performance of the proposed MLNN observer-based SMC control scheme.

1. Introduction

Recently, unmanned aerial vehicles (UAVs) gained considerable attention among engineers and researchers, with applications in several industries such as environmental monitoring, disaster response, precision agriculture, infrastructure inspection, and numerous other domestic applications [1,2]. Recent technological advancements in UAVs have enabled cost-effective autonomous navigation [3]. Controlling UAVs poses a greater challenge due to the presence of adverse forces like aerodynamic factors (such as wind gusts) and sensor noise, which can potentially cause instability or vehicle crash. Hence, addressing these issues is crucial for ensuring safe navigation.
In the existing work, a variety of control algorithms and architecture have been presented to improve the tracking accuracy of UAV systems. The strategies encompass adaptive control [4], reinforcement learning [5], data-driven based fault-tolerant control [6], fuzzy fractional PID controller [7], and adaptive neural network based control [8]. Furthermore, studies have explored ways to enhance control strategies for accurately following quadrotor flight paths under the influence of disturbances. Some recent control techniques include disturbance observer based H∞ control proposed in [9], disturbance observer-based control for intense maneuvering [10], a non-singular finite-time tracking controller using disturbance observers [11], robust backstepping sliding mode control [12], adaptive finite-time backstepping control which estimates disturbance [13], linear quadratic controllers with integral action [14], Fast Terminal Sliding Mode Controller (FTSMC) [15], and adaptive fuzzy control integrated with finite-time control methods [16].
To develop an efficient control system for UAVs, it is crucial to have accurate understanding of the plant states in order to improve control and navigation. While it is possible to use suitable sensors to collect this data, their expensive price and susceptibility to sensor noise and external influences make this method unfeasible in some situations. To solve this problem, various observer- schemes are proposed for UAV systems, such as Kalman filters based state estimation [17,18,19], sliding mode observers [20,21], and high-gain observers [22,23,24]. Nevertheless, they come with certain limitations. The classical Kalman filter requires the complete and accurate dynamics of the system, whereas the Extended Kalman Filter (EKF) and Unscented Kalman Filter (UKF) extends the Kalman filter to handle non-linearity—through linearization or unscented transforms. While classical Kalman filters rely on the accurate system model to provide precise state estimates, they can exhibit high susceptibility to noise and parameter variations. On the other hand, sliding mode observers (SMOs) and high-gain observers (HGOs) are robust to model uncertainties and require a general model structure. However, SMO can suffer from chattering, which introduces high-frequency oscillations, while HGOs are sensitive to noise and require precise gain tuning to prevent instability or degraded performance.
Neural network (NN)-based methods have emerged as potential solutions to overcome the constraints of traditional observers. These methods have become more important in UAV applications [25,26], which demand extremely reliable control and state estimation strategies. In recent years, several NN-observer-based methods have been proposed. For example, the authors in [27] designed fuzzy model predictive control utilizing an NN observer for a constrained hypersonic vehicle. The work in [28] presented a disturbance observer employing a Radial Basic Neural Network (RBFNN) for estimating unknown disturbances for tracking control of multiple UAVs, ensuring finite-time convergence and distributed coordination. Similarly, [29], introduces a disturbance observer based on RBFNN for controlling the attitude and altitude of a UAV amidst disturbances and model uncertainties. In [30], an SMO-based NN method is introduced for a quadcopter to mitigate the impacts of bounded disturbances. In [31], an intelligent sensor scheme based on NN is proposed for attitude estimation within a Global Positioning System (GPS) navigation model, which enhances navigation accuracy by compensating for delay-dependent errors. In [32], a Single Hidden Layer Neural Network (SHLNN)-based observer is suggested for a quadrotor UAV, incorporating a sigmoid activation function to enhance robustness against measurement noise. To address disturbances and uncertainty, a controller using an NN observer is proposed for a missile proposed in [33] holding worthwhile results. Additionally, [34] introduces a multi-layer neural network-based observer for estimating the states of an induction motor. In [35] a disturbance observer is combined with a multi-layer exponential sliding mode controller to address disturbances and parameter uncertainties encountered in a UAV. In [36], a deep learning-based Kalman filter for state estimation in multi-rotor UAVs is proposed and evaluated under the influence of noise. In [37], a single-layer neural network-based observer for a UAV using a back-propagation weight update law is used to eliminate disturbances effects. The work in [38] introduces a back-propagation-based attitude control of a quadrotor, with the study focusing solely on one varying parameter, namely mass variation, in their simulation results. The study in [39] integrates an adaptive neural network and an extended state observer to handle system uncertainties and external disturbances. The proposed observer scheme employs a SHLNN to estimate the disturbances.
Although many control techniques have been designed and tested, most of these approaches are primarily validated through numerical simulations rather than actual real-world applications. This gap emphasizes the necessity for more actual testing and validation in order to connect theoretical modes and operational performance, ensuring the reliability and effectiveness of quadrotor systems in diverse applications. Some recent work on real-time implementation of control algorithms on UAV based on full state feedback includes the implementation of second-order SMC in [40], integral backstepping-based SMC in [41], super-twisting SMC [42,43], adaptive PID control in [44], and model predictive control in [45]. The study in [46] investigates three nonlinear control methods: an existing strategy, Nonlinear Dynamic Inversion (NDI), and Incremental NDI (INDI) for stabilizing a tail-sitter UAV during hover. These were simulated and validated using Hardware-in-the-Loop (HITL) and experimental flight test results. Despite the mentioned scheme has been successfully deployed. However, there is still room for investigation specifically in the presence of unmeasured system states, measurement noise, and perturbation.
This research introduces an adaptive multi-layer neural network (MLNN) observer-based fast terminal SMC scheme for quadrotor UAVs assuming that only the positions of the plant are measured. Initially, the adaptive MLNN observer is employed for estimating the system’s states accurately. Subsequently, the estimated state information is utilized by the nonlinear dynamic inversion (NDI) based fast terminal SMC law to enhance system robustness against unmodeled dynamics and perturbation. The Lyapunov stability proof of the proposed MLNN observer is also provided, which is an extension of the Lyapunov proof for a single-layer NN proposed in [34]. To show the practical applicability of the suggested technique, hardware-in-loop simulations are performed using Pixhawk 6X flight controller and mission planner software. Finally, tests are performed on real quadrotor F450 system to show the efficacy of the suggested scheme. The key contributions of the proposed method are detailed below:
  • In comparison to the state observers cited in [9,17,22], the designed MLNN observer can estimate the system state information without requiring precise knowledge of the system dynamics and ensures finite-time convergence.
  • The NDI-based fast terminal SMC method offers significant benefits in terms of robustness, convergence speed, handling nonlinear dynamics, computational efficiency, and reduced chattering.
  • The hardware-in-loop simulation (HILS) results using Pixhawk 6X flight controller demonstrate good performance of the proposed scheme. In addition, flight tests performed on a real quadrotor demonstrate commendable tracking accuracy.

2. Dynamics of Quadrotor UAV

A quadrotor possesses six DOF but is equipped with only four propellers, rendering it an under-actuated system. The assumptions outlined for the quadrotor in this study are taken from [47]. A quadrotor UAV encompasses twelve states, namely, [ x ,   y ,   z ,   x ˙ ,   y ˙ ,   z ˙ ,   ϕ ,   θ ,   ψ ,   p ,   q ,   r ] , where x ,   y , and z represent the translational positions, while ϕ ,   θ and ψ , denote the roll, ptich and yaw angles, and p ,   q and r, are their angular velocities Figure 1 shows the structure of an F450 quadrotor.
Its nonlinear equations consists of a rigid body kinematic, dynamic, control effectiveness and propulsor model as follows.

2.1. Kinematic Model

Lets p e = x y z T be the position vector and v e = x ˙ y ˙ z ˙ T denotes the velocity vector in an inertial frame:
p ˙ e = v e
Similarly, Θ b = ϕ θ ψ T denotes the attitude vector and ω b = p q r T represents the velocity vector in body fixed frame. The relationship between attitude rate and angular velocities is given as
Θ ˙ b = R W ω b
where
R w = 1 S ϕ T θ C ϕ T θ 0 C ϕ S ϕ 0 S ϕ C θ C ϕ C θ
where ϕ ,   θ and ψ denote the roll, pitch and yaw angle, and p ,   q and r represents their angular rates.

2.2. Dynamics Model of Quadrotor System

The position dynamics is represented as:
v ˙ e = m G e + R e ( T e + F d e )
where G e = [ 0 , 0 , g ] T ,   T e = [ 0 , 0 , T 1 ] T , T 1 represents the total thrust developed by the propellers and is unidirectional, F d e is the vector representation of aerodynamic forces in the Earth frame and R e denotes the rotational matrix that transform vector from body frame to Earth frame and is defined as follows:
R e = 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 θ
The attitude dynamics of the quadrotor system in the body-axes frame are defined as follows:
I · ω b = ω b × ( I · ω b ) + τ g + τ b + M d
where τ b = [ τ 2 , τ 3 , τ 4 ] T denotes the moments generated by propeller, and I = d i a g ( I p x , I p y , I p z ) denotes the moment of inertia matrix, M d defines the aerodynamics torque, and τ g denotes the gyroscopic torque generated by rotating propellers can be expressed as
τ g = I r ( Ω r × ω b )
The term Ω r denotes the residual speed of rotors, defined as
Ω r = Ω 1 + Ω 2 Ω 3 + Ω 4
where, Ω 1 2 ,   Ω 2 2 ,   Ω 3 2 , and Ω 4 2 speed of rotors 1, 2, 3 and 4 respectively.

2.3. Thrust and Moment Model

The relationship between input torque and angular speed of individual propellers are described as under:
T τ 2 τ 3 τ 4 = b b b b b l b l b l b l b l b l b l b l d d d d Ω 1 2 Ω 2 2 Ω 3 2 Ω 4 2
To facilitate the controller design, a small angle approximation is considered. By redefining the state variables
x 1 = x , x 2 = x ˙ , x 3 = y , x 4 = y ˙ x 5 = z , x 6 = z ˙ , x 7 = ϕ , x 8 = p , x 9 = θ , x 10 = q , x 11 = ψ , x 12 = r ,
The state space form of UAV becomes:
x ˙ 1 = x 2 x ˙ 2 = ( C x 7 C x 9 C x 11 + S x 7 S x 11 ) T 1 m x ˙ 3 = x 4 x ˙ 4 = ( C x 7 S x 9 S x 11 S x 7 C x 11 ) T 1 m x ˙ 5 = x 6 x ˙ 6 = g + ( C x 7 C x 9 ) T 1 m
x ˙ 7 = x 8 x ˙ 8 = x 10 x 11 I p y I p z I p x I p r I p x x 10 Ω r + τ 2 I p x x ˙ 9 = x 10 x ˙ 10 = x 8 x 12 I p z I p x I p y I p r I p y x 8 Ω r + τ 3 I p y x ˙ 11 = x 12 x ˙ 12 = x 8 x 10 I p x I p y I p z + τ 4 I p z
The parameters of the quadrotor F450 are given in Table 1 and are taken from [48].
Remark 1.
In this paper, all the parameters given in Table 1 are known, and the measurements of the linear positions and Euler angles are assumed to be available. The goal is to develop such a control law that ensures that linear coordinates and yaw angle accurately track the desired positions.

2.4. A Multi-Layer Neural Network Architecture

In Figure 2, a common structure of an MLNN is depicted consisting of L layers where o 0 = o 0 1 o 0 2 o 0 n [ 0 ] is input layer, o L = o L 1 o L 2 o L n [ L ] is output layer, and o i = o i 1 o i 2 o i n [ 1 ] is the hidden layer. n [ i ] are the units in the i t h layer, for i = 1 , 2 , L . Each layer o 1 , L is expressed in the following form:
o i = λ W p i o i 1 + b p i
where W p i , and b p i are the weights and bias values for the MLNN layers, and λ ( ) represents the sigmoid activation function:
λ ( ) = 1 1 + e x p ( )

3. Proposed Method

The overall control strategy is illustrated in Figure 3. In this work, we utilize a dual-loop control architecture to achieve the target x-y translation and altitude control of UAV system. The outer-loop tracks the x and y position by providing the desired roll ϕ d and pitch θ d commands to the inner-loop, where the inner-loop tracks the altitude, ϕ d ,   θ d , and ψ d inputs by generating control inputs T 1 ,   τ 2 ,   τ 3 and τ 4 .

3.1. Multi-Layer Neural Network-Based Luenberger Observer Design

This section proposes an adaptive MLNN-based Luenberger observer for the quadrotor system to estimate the unknown system states p ,   q , and r in the attitude dynamic system. The attitude dynamics of the quadrotor system, given in (11), can be written to the form
x ˙ p ( t ) = A p x p ( t ) + G p ( x p , u p , d p ) y p ( t ) = C p x p ( t )
where x p ( t ) = c o l ( x 7 , x 8 , x 9 , x 10 , x 11 , x 12 ) R 6 , u p ( t ) = ( τ 2 , τ 3 , τ 4 ) T R 3 , d p ( t ) = Ω r R 1 and y p ( t ) R 3 are the states, input, disturbance, and output vectors respectively. In (14), A p = 0 3 I 3 0 3 0 3 denotes system matrix, G p ( x p , u p , d p ) = c o l x 10 x 11 I p y I p z I p x I p r I p x x 10 Ω r + τ ϕ I p x , x 8 x 12 I p z I p x I p y I p r I p y x 8 Ω r + τ θ I p y , x 8 x 10 I p x I p y I p z + τ ψ I p z represents a nonlinear function contain unknown states which is to be estimated by the MLNN observer and C p = I 3 0 3 × 3 denotes output matrix respectively. By following the results of Kolmogorov theorem, a non-linear function G p ( x , u , d ) is described in terms of activation function and ideal weights [37]:
G p ( x p , u p , d p ) = W L λ 0 ( W p L 1 , W p 1 , b p L 1 , b p 1 , o 0 ) + b p L + ( x p )
where ( x p ) is an neural network approximation error satisfying ( x p ) ) N , where N is a positive number. An assumptions considered for designing the observer-based controller is that since the pair ( A p , C p ) is observable (rank(obsv( A p , C p )) = 6), the proposed NN-based Luenberger observer structure for attitude dynamic system (11) is given as [37]:
x ^ ˙ p ( t ) = A p x ^ p ( t ) + L M y p ( t ) y ^ p ( t ) + G ^ p ( x ^ p , u p , d p ) G ^ p ( x ^ p , u p , d p ) = W ^ p L λ 0 W ^ p L 1 , W ^ p 1 , b ^ p L 1 , , b ^ p 1 , o ^ 0 + b ^ p L y ^ p ( t ) = C p x ^ p ( t )
where x ^ p ( t ) represents the estimate of x p ( t ) , y ^ p ( t ) denote the estimated output, L M R 6 × 3 represents the observer gain and G ^ p ( x ^ p , u p , d p ) is the unknown function to be estimated by the MLNN observer, whereas W ^ p i and b ^ p i define the NN-based learned weight matrices and bias vectors respectively shown in Figure 2. The error is e p ( t ) = x p ( t ) x ^ p ( t ) , and its derivative e ˙ p ( t ) after substituting (11) and (16) and adding and subtracting the term W ^ p L λ 0 W ^ p L 1 , W ^ p 1 , b ^ p L 1 , , b ^ p 1 , o ^ 0 is obtained as
e ˙ p ( t ) = H p e p ( t ) + e p b L + Υ ( t ) + e w p L λ 0 W ^ p L 1 , W ^ p 1 , b ^ p L 1 , b ^ 1 , o ^ 0
where e W p L = W p L W ^ p L , e p b L = b p L b ^ p L , H p = A p L M C p and Υ ( t ) denotes a bounded disturbance which ensures Υ Υ ¯ , for a positive constant Υ ¯ , given as:
Υ ( t ) = W p L λ 0 W p L 1 , , W p 1 , b p L 1 , b p 1 , o 0 λ 0 W ^ p L 1 , W ^ p 1 , b ^ p L 1 , b ^ p 1 , o ^ 0 + ( x p ) . This study utilizes the following cost function.
χ p = 1 2 e y T e y
By employing the chain rule on (18) we get:
𝜕 χ p 𝜕 W ^ p i = d i a g ( 𝜕 C p 𝜕 o ^ i ) T 𝜕 o ^ i 𝜕 o ^ i
where 𝜕 χ p 𝜕 o ^ i = 𝜕 χ p 𝜕 e y 𝜕 e y 𝜕 e p 𝜕 e p 𝜕 o ^ L 𝜕 o ^ L 𝜕 o ^ L 1 𝜕 o ^ L 1 𝜕 o ^ L 2 𝜕 o ^ i + 1 𝜕 o ^ i , 𝜕 χ p 𝜕 e y = e y T , 𝜕 e y 𝜕 e p = C p and 𝜕 e p 𝜕 o ^ L = H p 1 . For e y ( t ) = C p e p ( t ) , after putting 𝜕 χ p 𝜕 e y , 𝜕 e y 𝜕 e p , 𝜕 e p 𝜕 o ^ L in (19), we obtain:
𝜕 χ p 𝜕 W ^ p i = d i a g e p T C p T C p H p 1 𝜕 o ^ L 𝜕 o ^ L 1 𝜕 o ^ L 1 𝜕 o ^ L 2 𝜕 o ^ i + 1 𝜕 o ^ i T 𝜕 o ^ i 𝜕 W ^ p i
where
𝜕 o ^ L 𝜕 o ^ L 1 = W ^ p L 𝜕 o ^ i 𝜕 W ^ p i = λ W ^ p i o ^ i 1 + b ^ p i o ^ i 1 T

3.2. Multi-Layer Neural Network Weight Update Law

This work employs a weight-tuning strategy that combines a modified back-propagation algorithm with gradient scaling, as detailed in references [34,37].
W ^ ˙ p i = ϑ i j = i + 1 L 1 W ^ p j + ς 𝜕 χ p 𝜕 W ^ p i Υ e p W ^ p i b ^ ˙ p i = μ i b ^ p L j = i + 1 L W ^ p j + ς 𝜕 χ p 𝜕 b ^ p i Υ e b ^ j i
where ϑ i L , and μ i L represent the learning rates, while Υ is a positive constant. Inserting (20) in (21):
W ^ ˙ p i = ϑ i j = i + 1 L 1 W ^ p j + ς d i a g ( 𝜕 χ p 𝜕 o ^ i T ) λ ´ ( W p i o i 1 + b p i ) o i 1 T Υ e p W ^ p i b ^ ˙ p i = μ i b ^ p L j = i + 1 L W ^ p j + ς d i a g ( 𝜕 χ p 𝜕 o ^ p i T ) λ ´ ( W p L 1 o L 2 + b p L 1 ) Υ e p b ^ p i
where ς is chosen to be smaller. Consequently, the output layer’s adaptive law is derived as
W ^ ˙ p L = d i a g ( e p T T L ) T o L 1 T Υ e p W ^ p L
Theorem 1.
Suppose P o is a positive definite, while Q o is another matrix, satisfying the equality:
H p T P o + P o H p = Q o
then the error (17) remains within bounds if the subsequent condition is fulfilled:
e p > 2 P λ + i = 1 L 1 Υ L 1 Υ i 2 Υ 2 i 2 + i = 1 L 1 Υ 1 Υ 3 i 2 + i = 1 L 1 Υ L 1 ϱ 2 i 2 i = 1 L 1 Υ 1 ϱ 3 i 2 ζ m i n ( Q )
Proof. 
To establish the proposed MLNN-observer’s stability, we extend the single-layer neural network architecture [34,37] to a multi-layer with the modified back-propagation weight update rule. In Appendix A, the proof of this theorem is provided. □
Remark 2.
This study develops an MLNN-based observer (16) to predict the internal states of the quadrotor, such as angular rates. The observer error dynamics (17) is stable if the conditions given in (24)–(25) are satisfied.

3.3. NDI Based Fast Terminal SMC Law

This section proposes an NDI-based fast terminal SMC law to attain the position and attitude-tracking control of a UAV. Since the attitude dynamic system depends on the estimated states obtained from the MLNN observer, we first consider the detailed formulation of the attitude dynamic control law. Meanwhile, the other control laws are designed along similar lines. The proposed control law is designed in two steps. First, the NDI controller is developed to cancel the nonlinearities without considering the disturbance effect. Then, the adaptive fast terminal SMC law is integrated, which provides robustness against disturbance/uncertainty. Therefore, NDI-based fast terminal SMC control law can be written in the following form:
u p ( t ) = u p n ( t ) + u p l ( t )
where u p l ( t ) is the NDI component of the control law without considering disturbances, and u p n ( t ) denotes the fast terminal SMC component, which is the adaptive part to deal with disturbances. According to (11), the attitude dynamics of the quadrotor are affine in input. Since, the attitude angle states are available at the output, the output dynamics of the attitude subsystem (14) can also be written into the form
y ˙ p ( t ) = C p A p x p + C p B p u p + C p G p ( x p , t ) + C p D p d p ( t )
where B p = [ 0 3 × 3 d i a g ( 1 / I p x , 1 / I p y , 1 / I p z ) ] T , C p = [ I 3 × 3 0 3 × 3 ] , D p = c o l ( 0 3 × 3 , I p r I p x x 10 , I p r I p y x 8 , 0 ) and d p ( t ) = Ω r .

3.3.1. NDI Control Law

The attitude dynamics of the quadrotor (12), after neglecting the effect of disturbance/uncertainty, can be written as
y ˙ p ( t ) = C p A p x p + C p B p u p + C p G p ( x p , t )
Define the desired attitude states y p d [ ϕ d , θ d , ψ d ] T and its derivatives y ˙ d . The control input is then designed as
u p n = ( C p B p ) 1 ( C p A p x p + C p G p ( x ^ p , t ) + y ˙ d )
where x ^ p denotes the estimated state acquired from the MLNN observer and y ˙ d denotes the desired dynamics defined as
y ˙ d = k p e y + k i e y d t
where e y = y d y denotes tracking error, k p = d i a g ( k p ϕ , k p θ , k p ψ ) and k i = d i a g ( k i ϕ , k i θ , k i ψ ) are NDI controller parameter that determine the desired attitude dynamics.

3.3.2. Fast Terminal SMC Law

The next step is integrating fast terminal SMC law to enhance system robustness. Therefore, the fast terminal SMC law u p S ( t ) is defined as
u p s ( t ) = K f s i g n ( s p ) η f s p β f
where K f = d i a g ( κ f ϕ , κ f θ , κ f ψ ) , η f = d i a g ( η f ϕ , η f θ , η f ψ ) are the sliding gains and 0 < β f < 1 is fast terminal SMC design parameters, s p ( t ) 3 × 1 denotes the sliding surface defined as
s p ( t ) = e y ( t ) + λ f e y α f ( t )
where λ f > 0 and 0 < α f < 1 are design parameters.
Remark 3.
The parameters λ f and η f of fast terminal SMC law play crucial roles in determining the effectiveness and robustness of the controller design. The parameter λ f influences on convergence rate and transient response while η f impacts robustness and steady-state error. Achieving an optimal balance between control effort, robustness and swift response in the control system requires precise tuning of these parameters.

3.4. The Stability Analysis

Substitute the control law (26)–(31) into the system dynamics (27), we get
e ˙ y ( t ) = K f s i g n ( s p ) η f s p β f + y ˙ d + C p ( G p ( x p , t ) G ^ p ( x ^ p , t ) ) + C p D p d p ( t )
The derivative of switching function (32) after substituting the (33) is obtained as
s ˙ p = K f s i g n ( s p ) η f s p β f + y ˙ d + C p ( G p ( x p , t ) G ^ p ( x ^ p , t ) ) + C p D p d p ( t ) ( 1 + λ f α f e y α f 1 )
To prove system stability, define the Lyapunov function as
V p = 1 2 s p T ( t ) s p ( t )
Taking the time derivative of Lyapunov function and substitute (34), we get
V ˙ p = s p T K f s i g n ( s p ) η f s p β f + y ˙ d + C p ( G p ( x p , t ) G ^ p ( x ^ p , t ) ) + C p D p d p ( t ) ( 1 + λ f α f e y α f 1 )
Assuming the nonlinear function G p ( x p , t ) satisfies the Lipschitz condition than we can write [49,50]
C p ( G p ( x p , t ) G ^ p ( x ^ p , t ) ) L f e y ( t )
where L f is Lipschitz constant. Utilizing (37), we get
V ˙ p s p T K f η f s p β f + y ˙ d + L f e y ( t ) + C p D p d p ( t ) ( 1 + λ f α f e y α f 1 )
To ensure V ˙ p < 0 , the following conditions on K f and η f must be satisfied
K f y ˙ d + L f e y ( t ) + C p D p d p ( t ) 1 + λ f α f e y α f 1 , η f 0
Hence, V ˙ p < 0 . This ensures the convergence of the sliding surface s p ( t ) to zero, thereby stabilizing the system despite disturbances.
Remark 4.
The final control law (26) is derived by incorporating the NDI control law (29) and fast terminal SMC law (31), provided that the sliding gains are chosen according to (39), and the estimated state information is obtained from (16).
The outer loop position control is achieved by generating the desired roll ϕ d and pitch commands θ d commands to the inner loop subsystem. The control laws for outer loop is designed in similar way as
u x = m T 1 ( x ˙ d k x s i g n ( s x ) η x s x β x ) u y = m T 1 ( y ˙ d k y s i g n ( s y ) η y s y β y ) T = m C x 7 C x 9 ( z ˙ d k z s i g n ( s z ) η z s z β z )
where x d ,   y d and z d are the desired x ,   y and z positions, k x ,   k y ,   k z and η x ,   η y ,   η z are the sliding gains, and s x ,   s y ,   s z are the sliding surfaces defined as:
s x = e x + λ f x e x α f x ( t ) s y = e y + λ f y e y α f y ( t ) s z = e z + λ f z e z α f z ( t )
here e x = x d x ,   e y = y d y ,   e z = z d z . ϕ d and θ d are the desired roll and pitch commands calculated as
ϕ d = s i n 1 u x s i n ψ d u y c o s ψ d θ d = s i n 1 u x s i n ψ d u y c o s ψ d c o s ϕ d

4. Parameter Setting and Results

This section discusses the efficacy of the scheme designed in the previous section by testing it on a nonlinear model of the quadrotor system. To validate the simulation results, HIL simulations are performed using Pixhawk 6X flight controller. Then, the experiments are conducted on real-time hardware setup of F450 in a constrained indoor environment.

4.1. Parameter Settings

The MLNN observer is trained using data generated from the quadrotor UAV model. The inputs to the MLNN are ϕ ,   θ ,   ψ , while the outputs are p ^ ,   q ^ , and r ^ , as depicted in Figure 4. The neural network undergoes training using input-output data to approximate states closely resembling the actual ones. The neural network is structured with four hidden layers, each comprising 10, 12, 12, and 20 neurons, respectively, all employing a sigmoid activation function. The network’s initial weights and bias terms are set with small random values, and these weights are subsequently refined and adjusted using the modified back-propagation algorithm. During training, the input data is processed through the network’s layers, and the difference between predicted and actual outputs is calculated. This error is then propagated backward through the network to compute gradients, which guide the adjustment of weights and biases via a modified backpropagation algorithm. The network undergoes multiple iterations or epochs of this process, refining its parameters to minimize the error and improve predictive accuracy so that the estimated states closely resemble the actual ones, thus improving the network’s ability to make accurate predictions. After training, the trained MLNN model is placed in run-time simulations to estimate states and evaluate the observer’s performance. The controller and observer gains are given in Table 2.

4.2. Results

This subsection initially evaluates the performance of the suggested scheme with numerical simulations on the nonlinear model (1)–(9). Then, HIL simulations are performed using Pixhawk 6X autopilot. Finally, real tests are performed on the F450 quadrotor system mounted on the frame for safe indoor environment.

4.3. Numerical Simulations

In this subsection, the effectiveness of the proposed method is tested on a nonlinear model of aircraft given in (1)–(6). A sine and cosine function of amplitude 6 m is applied to desired x and y positions. The desired z position ramps up from 0 to 10 m over 50 s, reaching a height of 5 m at 60 s, and remains constant thereafter. Our trajectory represents a gradually ascending circular path, which is commonly used in realistic UAV applications to generate circular ascending flight patterns [51,52,53]. This data is used to train the neural network. The proposed method is tested by considering two different conditions; the first one is nominal scenarios, and the other considers the effect of disturbance and parameter variations. To evaluate the performance of the proposed scheme subject to parameter variation, a performance index is defined as
e M S E i = 1 t f t i t i t f e i ( t ) 2 d t
where e M S E i represents the mean square error and e i is the error between i t h desired and state estimated by the proposed scheme.

4.3.1. Nominal Condition

In this subsection, the simulations are conducted without incorporating noise or wind gusts into the plant dynamics. The results depicted in Figure 5 illustrate the state tracking results of the proposed observer-based controller under nominal conditions. The comparison is made with the adaptive neural network extended state observer (ANNESO) employing a single hidden layer neural network architecture given in [39]. The simulation results showcase outstanding tracking performance for the x ,   y , and z positions and the attitude angles for the proposed scheme. The plots of control input are depicted in Figure 6, and the mean square error calculation in the nominal condition is given in Table 3.

4.3.2. Results with Noise, Disturbance and Parameter Variations

To validate the performance of the proposed MLNN observer-based control scheme, we presented three different cases for parameter variation ( 5 % , 10 % , and 15 % ) and wind-gust condition (mean velocity of 4 m/s, 6 m/s and 8 m/s), For each case, the noise with a variance of 0.02 and zero mean values is added to the p ,   q ,   r and z states. Wind gusts effect using the Dryden wind gust model [54] are added into the non-linear plant equations.
Figure 7 and Figure 8 showcases the tracking of states and control input results respectively under noise, wind gust, and parameter variation across cases A, B, and C, respectively. Noise and external aerodynamic effects during flight can lead to variations in the parameters of UAVs, including mass, moments of inertia along the x ,   y , and z axes, drag forces, and coefficient of thrust. The adaptive gains of the MLNN observer mitigate the impact of such factors, facilitating the estimation of trajectories that closely resemble the reference trajectories. Despite parameter variations in Case B and C causing slight deviations from the desired trajectory, the proposed adaptive MLNN effectively learns and rejects the nonlinearities. Additionally, the NDI-based Fast Terminal SMC incorporates robustness to further enhance performance. The comparison of the proposed scheme is made with an existing work presented in [44]. The results are compared considering only worst case scenario of parameters variations and disturbance presented in Case C. These simulation plots indicate that the proposed scheme exhibits less deviation and noise impact due to its adaptive nature, compared to [44]. The MSE values for the Cases A, B, and C are provided in Table 4. Additionally, the MSE comparison between the proposed scheme and [44] for Case C is also provided which clearly shows that the MSE values of the proposed scheme are lower than those of the scheme suggested in [44] which ensure better tracking accuracy.

4.4. Hardware-in-the-Loop Simulations

In this section, we conducted HIL simulations utilizing the Pixhawk 6X flight controller to showcase the algorithm’s viability in practical scenarios. These simulations offer concrete proof of the effectiveness and applicability of our control algorithm during real-time UAV operations. The Pixhawk flight controller board is interfaced with the Mission Planner platform. These simulations aimed to validate the controller design implemented on the Pixhawk hardware board, incorporating UAV and MLNN observer dynamics within Simulink. Flight missions were programmed in QGroundControl (QGC), guiding the flight controller algorithm to follow the designated missions. The connection setup for HIL simulations is given in Figure 9, which illustrates the communication setup between MATLAB, Pixhawk, and QGround Control. The communication relies on the Mavlink protocol, recognized as the standard communication protocol for UAVs. In HIL architecture shown in Figure 10, the control algorithm operates on the Pixhawk 6X, while the computational dynamics of the UAV and observer dynamics are performed on a host computer. Simulink’s Monitor and Tune simulation tool is utilized for real-time logging and tuning of controller settings. The host computer computes dynamics, transmitting the data to the PX4 Autopilot via the USB port. The PX4 autopilot calculates thrust and torque based on the received data, establishes a hardware loop, and sends them back to the host. An additional serial port on the Pixhawk 6X is employed to establish communication between the PX4 Autopilot and Simulink for utilizing the monitor and tune tool. The altitude and attitude tracking results of the proposed scheme are showcased in Figure 11, while position tracking results are given in Figure 12. The results clearly demonstrate that the proposed MLNN observer-based NDI fast terminal SMC scheme closely follows the desired trajectories. The control inputs displayed in Figure 13 are bounded within a specific range.

4.5. Experimental Results on Quadrotor Test-Rig

The real-time hardware setup comprises a F450 X-configuration quadrotor frame, Pixhawk 6X Autopilot, RC Transmitter and receiver (FS-i6) given in Figure 14. The propeller system includes 4 permanent magnet synchronous motors (PMSM) with KV-950 and 20A Esc, powered by a 5000 mAh Li-po battery with a low-battery tester. The experimental setup is mounted on a safety stand that is only capable of performing attitude tracking. The Pixhawk 6X is equipped with a comprehensive sensor suite for precise navigation and control. It includes three ICM-42688 IMUs, featuring balanced gyro technology for measuring linear acceleration and angular velocity which are essential for flight stabilization and smooth movement. The system also incorporates two barometers, the ICP20100 and BMP388, for accurate altitude measurement, and a BMM150 magnetometer to detect orientation relative to Earth’s magnetic field. These sensors work together to provide reliable flight data, enabling autonomous operations and real-time flight adjustments. The quadcopter is securely mounted on a safety aluminum stand with a spherical ball joint to ensure the safety of the environment and personnel during testing, allowing observation of three primary movements: roll, pitch, and yaw. The MLNN-based observer is simulated using a high-processing laptop, which communicates with the Pixhawk via a wired RS232 protocol. The control algorithm operates on the Pixhawk 6x, with all necessary hardware devices interfaced accordingly. The desired roll, pitch, and yaw commands are inputted to the quadrotor using the RC transmitter, while the RC receiver receives Pulse Position Modulation (PPM) commands from the transmitter and forwards them as reference signals to the autopilot. The Pixhawk autopilot executes the desired attitude commands based on the state inputs from the MLNN observer and sensors, generating thrust and torque commands in the form of Pulse Width Modulation (PWM) signals for the motors. The safety stand allowed us to capture real-time flight data without compromising the accuracy of the results. Thus, the data presented reflects practical flight behavior, despite the UAV being physically secured for safety reason. The attitude tracking results with wind disturbance are depicted in Figure 15, demonstrating commendable performance when implemented on a real F450 quadrotor.

5. Conclusions

This work proposes an adaptive MLNN-observer-based NDI fast terminal SMC scheme for quadrotor UAV. The adaptive MLNN-observer is based on a modified back-propagation algorithm that alleviates the impacts of wind disturbances, parameter uncertainties, and noise by learning the non-linearities in the data. A rigorous stability analysis is performed to guarantee the convergence of observer error dynamics. Then, NDI-based fast terminal SMC laws are designed with both inner and outer loops, which provide attitude and position tracking control while contributing towards additional robustness. To test the efficacy of the proposed approach, simulations are first performed in different operating conditions. To validate the numerical simulations, HIL simulations are conducted using Pixhawk 6X flight controller. Finally, the hardware test are performed in constraint indoor environment hardware-in-loop simulations using a real quadrotor F450 platform to validate the effectiveness of the proposed method. In future, the MLNN based observer will be extended to estimate both sensor and actuators fault and external disturbance and closed loop performance will be analyzed in the presence of sensor and actuator malfunctioning.

Author Contributions

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

Funding

This work is supported by Natural Science Foundation of China with project code 62250410367.

Data Availability Statement

No data available for this study.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

We define the Lyapunov candidate function V p as follows:
V p = 1 2 e p T P e p + i = 1 L 1 2 t r ( e W p i T e W p i ) + i = 1 L 1 2 e b p i T e b p i
The derivative of (A1) results in:
V ˙ p = 1 2 e ˙ p T P e + 1 2 e p T P e ˙ p + i = 1 L 1 2 t r ( e p W i T e ˙ W p i ) + i = 1 L 1 2 e p b p i T e ˙ b p i
where e p is described in (17), e ˙ W p i = W ˙ p i W ˙ ^ p i and e ˙ b p i = b ˙ p i b ˙ ^ i . After replacing e ˙ W p i , and e ˙ b p i in (A2), the following equation is obtained:
V ˙ p = 1 2 e p T Q e + e p T P [ e W p L λ 0 ( W p L 1 , , W p 1 , b p L 1 , b p 1 , o 0 ) + λ ( t ) + e b p L ] + i = 1 L 2 t r [ e W p i T ϑ i j = i + 1 L 1 W ^ p j + ς d i a g ( 𝜕 χ p 𝜕 o ^ i T ) λ ´ ( W p i o i 1 + b p i ) o i 1 T Υ e p W ^ p i ] + t r [ e w p L 1 T d i a g ( 𝜕 χ p 𝜕 o ^ L 1 T ) λ ´ ( W p L 1 o L 2 + b p L 1 ) o L 2 T Υ e p W ^ p L 1 ]   +   t r [ e w p L T d i a g ( e T T L ) T o L 1 T + Υ e p W ^ p L ] + i = 1 L 2 e b p i T [ μ i b ^ p L j = i + 1 L W ^ p j + ξ d i a g ( 𝜕 χ p 𝜕 o ^ i T ) λ ´ ( W p L 1 o L 2 + b p L 1 ) Υ e p b ^ p i ] + e b p L T [ M p L e p + Υ e p b ^ p L ]
Let λ λ 0 M , b p i b p M i and W p i W p M i for all i = 1 L 1 , we get:
V ˙ p 1 2 ζ m i n Q e p 2 + L 1 L 1 e P ( e W p L λ 0 M + λ ) + L 1 L 1 e p P e b p L + i = 1 L 1 e p e W p i 2 Υ i W ^ p L + Υ e p i = 1 L 1 W p M i e p w i e p w i 2 + L 1 L 1 e p × N p L e W p L   +   L 1 L 1 Υ e p ( W M p L e W p L e W p L 2 ) + i = 1 L 1 e p e b p i 2 ϱ i b ^ p L   +   Υ e p i = 1 L 1 b M p i e b p i e b p i 2 + L 1 L 1 e p e b p L   +   L 1 L 1 Υ e p ( b M p L e b p L e b p L 2 )
Let N p i = ϑ i B T H p T C p T C p , M p i = μ i B T H p T C p T C p for i = 1 L 1 , Υ 2 i = Υ W M p L + N p L + P λ 0 M 2 Υ Υ i 2 L 1 and Υ 3 i = Υ i W M p L + Υ W M p i 2 Υ 1 . Using the terms that include e W p i and e W p L in (A4). Let Υ i = N p i j = 1 L 1 λ M i , and adding and subtracting Υ 2 i 2 e p and Υ 3 i 2 e p | , we obtain:
e p ( 1 L 1 P λ + Υ L 1 Υ i 2 Υ 2 i 2 + ( Υ 1 ) ) Υ 3 i 2 Υ L 1 Υ i 2 Υ 2 i e W p L 2 Υ 1 Υ 3 i e W p i 2 Υ i e W p L e W p i 2 )
(a1)
We assume Υ L 1 > Υ i 2 and Υ > 1 Υ > ( L 1 ) Υ i 2 and Υ > 1 for i = 1 L 1 .
Let ϱ 2 i = Υ b M p L + 1 + P 2 Υ ϱ i 2 L 1 ) ) and ϱ 3 i = 2 ϱ i b M L + Υ b M p i 2 ( Υ 1 ) .
Using the terms that include e b p i and e b L in (A4). Let 2 ϱ i = M p i j = i L 1 λ m j , add and subtract the terms ϱ 2 i 2 e p and ϱ 3 i 2 e p , we get
e p ( e b p L 2 Υ L 1 ϱ i 2 + e b p L ( Υ b M p L L 1 + 1 L 1 + P L 1 ) + e b p i Υ b M p i + 2 ϱ i b M p L ( ϱ i e b p L e b p i 2 ) e b p i 2 Υ 1 ) = e p ( Υ L 1 ϱ i 2 ) ϱ 2 i 2 + Υ 1 ϱ 3 i 2 Υ L 1 ϱ i 2 ϱ 2 i e b p L 2 Υ 1 ϱ 3 i e b p i 2 ϱ i e b p L e b p i 2
(a2)
For Υ L 1 > ϱ i 2 and Υ > 1 Υ > ( L 1 ) ϱ i 2 and Υ > 1
If the conditions (a1) and (a2) hold:
V ˙ p 1 2 λ m i n Q e p 2 + e p P λ + i = 1 L 1 Υ L 1 Υ i 2 × Υ 2 i 2 + i = 1 L 1 Υ 1 Υ 3 i 2 + i = 1 L 1 Υ L 1 ϱ i 2 ϱ 2 i 2 i = 1 L 1 Υ 1 ϱ 3 i 2
V ˙ p is negative definite if:
e p > 2 P η + i = 1 L 1 Υ L 1 Υ i 2 Υ 2 i 2 + i = 1 L 1 Υ 1 Υ 3 i 2 + i = 1 L 1 Υ L 1 ϱ 2 i 2 i = 1 L 1 Υ 1 ϱ 3 i 2 γ m i n ( Q )
The condition in (A8) must be met for ensuring the negative definiteness of V p . So, according to the standard Lyapunov theorem, the error is bounded, confirming its convergence and thereby validating the observer’s stability.

References

  1. Wang, E.; Sun, J.; Liang, Y.; Zhou, B.; Jiang, F.; Zhu, Y. Modeling, Guidance, and Robust Cooperative Control of Two Quadrotors Carrying a “Y”-Shaped-Cable-Suspended Payload. Drones 2024, 8, 103. [Google Scholar] [CrossRef]
  2. Hanif, A.S.; Han, X.; Yu, S.-H. Independent Control Spraying System for UAV-Based Precise Variable Sprayer: A Review. Drones 2022, 6, 383. [Google Scholar] [CrossRef]
  3. Yang, T.S.; Yang, F.; Li, D. A New Autonomous Method of Drone Path Planning Based on Multiple Strategies for Avoiding Obstacles with High Speed and High Density. Drones 2024, 8, 205. [Google Scholar] [CrossRef]
  4. Hu, Z.; Jin, X. Adaptive formation control architectures for a team of quadrotors with multiple performance and safety constraints. Int. J. Robust Nonlinear Control 2023, 33, 8183–8204. [Google Scholar] [CrossRef]
  5. Zhao, W.; Liu, H.; Wang, B. Model-free attitude synchronization for multiple heterogeneous quadrotors via reinforcement learning. Int. J. Intell. Syst. 2021, 36, 2528–2547. [Google Scholar] [CrossRef]
  6. Zhao, W.; Liu, H.; Lewis, F.L. Data-driven fault-tolerant control for attitude synchronization of nonlinear quadrotors. IEEE Trans. Autom. Control 2021, 66, 5584–5591. [Google Scholar] [CrossRef]
  7. Xian, W.; Qi, Q.; Liu, W.; Liu, Y.; Li, D.; Wang, Y. Control of quadrotor robot via optimized nonlinear type-2 fuzzy fractional PID with fractional filter: Theory and experiment. Aerosp. Sci. Technol. 2024, 151, 109286. [Google Scholar] [CrossRef]
  8. Zhang, X.; Li, H.; Zhu, G.; Zhang, Y.; Wang, C.; Wang, Y.; Su, C.-Y. Finite-Time Adaptive Quantized Control for Quadrotor Aerial Vehicle with Full States Constraints and Validation on QDrone Experimental Platform. Drones 2024, 8, 264. [Google Scholar] [CrossRef]
  9. Lyu, X.; Zhou, J.; Gu, H.; Li, Z.; Shen, S.; Zhang, F. Disturbance Observer Based Hovering Control of Quadrotor Tail-Sitter VTOL UAVs Using H Synthesis. IEEE Robot. Autom. Lett. 2018, 3, 2910–2917. [Google Scholar] [CrossRef]
  10. Castillo, A.; Sanz, R.; Garcia, P.; Qiu, W.; Wang, H.; Xu, C. Disturbance observer-based quadrotor attitude tracking control for aggressive maneuvers. Control Eng. Pract. 2019, 82, 14–23. [Google Scholar] [CrossRef]
  11. El-Sousy, F.F.; Alattas, K.A.; Mofid, O.; Mobayen, S.; Asad, J.H.; Skruch, P.; Assawinchaichote, W. Non-singular finite time tracking control approach based on disturbance observers for perturbed quadrotor unmanned aerial vehicles. Sensors 2022, 22, 2785. [Google Scholar] [CrossRef]
  12. 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]
  13. Martins, L.; Cardeira, C.; Oliveira, P. Inner-outer feedback linearization for quadrotor control: Two-step design and validation. Nonlinear Dyn. 2022, 110, 479–495. [Google Scholar] [CrossRef]
  14. Liu, K.; Yang, P.; Wang, R.; Jiao, L.; Li, T.; Zhang, J. Observer-based adaptive fuzzy finite-time attitude control for quadrotor UAVs. IEEE Trans. Aerosp. Electron. Syst. 2023, 59, 8637–8654. [Google Scholar] [CrossRef]
  15. Derrouaoui, S.H.; Bouzid, Y.; Guiatni, M. Nonlinear robust control of a new reconfigurable unmanned aerial vehicle. Robotics 2021, 2, 76. [Google Scholar] [CrossRef]
  16. Koo, S.M.; Sands, T. Bilinear interpolation of three–dimensional gain–scheduled autopilots. Sensors 2023, 24, 13. [Google Scholar] [CrossRef]
  17. Xiong, J.-J.; Zheng, E.-H. Optimal kalman filter for state estimation of a quadrotor UAV. Optik 2015, 126, 2862–2868. [Google Scholar] [CrossRef]
  18. Zhang, T.; Liao, Y. Attitude measure system based on extended Kalman filter for multi-rotors. Comput. Electron. Agric. 2017, 134, 19–26. [Google Scholar] [CrossRef]
  19. Khamseh, H.B.; Ghorbani, S.; Janabi-Sharifi, F. Unscented Kalman filter state estimation for manipulating unmanned aerial vehicles. Aerosp. Sci. Technol. 2019, 92, 446–463. [Google Scholar] [CrossRef]
  20. Cai, W.; She, J.; Wu, M.; Ohyama, Y. Disturbance suppression for quadrotor UAV using sliding-mode-observer-based equivalent-input-disturbance approach. ISA Trans. 2019, 92, 286–297. [Google Scholar] [CrossRef]
  21. Zhao, Z.; Cao, D.; Yang, J.; Wang, H. High-order sliding mode observer-based trajectory tracking control for a quadrotor UAV with uncertain dynamics. Nonlinear Dyn. 2020, 102, 2583–2596. [Google Scholar] [CrossRef]
  22. Du, J.; Hu, X.; Liu, H.; Chen, C.P. Adaptive robust output feedback control for a marine dynamic positioning system based on a high-gain observer. IEEE Trans. Neural Netw. Learn. Syst. 2015, 26, 2775–2786. [Google Scholar] [CrossRef] [PubMed]
  23. Xu, H.; Wang, J. Distributed observer-based control law with better dynamic performance based on distributed high-gain observer. Int. J. Syst. Sci. 2020, 51, 631–642. [Google Scholar] [CrossRef]
  24. Nicosia, S.; Tornambè, A. High-gain observers in the state and parameter estimation of robots having elastic joints. Syst. Control Lett. 1989, 13, 331–337. [Google Scholar] [CrossRef]
  25. Muliadi, J.; Kusumoputro, B. Neural network control system of UAV altitude dynamics and its comparison with the PID control system. J. Adv. Transp. 2018, 2018, 3823201. [Google Scholar] [CrossRef]
  26. Gu, W.; Valavanis, K.P.; Rutherford, M.J.; Rizzo, A. UAV Model-based Flight Control with Artificial Neural Networks: A Survey. J. Intell. Robot. Syst. 2020, 100, 1469–1491. [Google Scholar] [CrossRef]
  27. Ma, Y.; Cai, Y. A fuzzy model predictive control based upon adaptive neural network disturbance observer for a constrained hypersonic vehicle. IEEE Access 2017, 6, 5927–5938. [Google Scholar] [CrossRef]
  28. Wang, D.; Zong, Q.; Tian, B.; Shao, S.; Zhang, X.; Zhao, X. Neural network disturbance observer-based distributed finite-time formation tracking control for multiple unmanned helicopters. ISA Trans. 2018, 73, 208–226. [Google Scholar] [CrossRef]
  29. Matassini, T.; Shin, H.-S.; Tsourdos, A.; Innocenti, M. Adaptive control with neural networks-based disturbance observer for a spherical UAV. IFAC-PapersOnLine 2016, 49, 308–313. [Google Scholar] [CrossRef]
  30. Stebler, S.; Campobasso, M.; Kidambi, K.; MacKunis, W.; Reyhanoglu, M. Dynamic neural network-based sliding mode estimation of quadrotor systems. In Proceedings of the IEEE 2017 American Control Conference (ACC), Seattle, WA, USA, 24–26 May 2017; pp. 2600–2605. [Google Scholar]
  31. Jaradat, M.A.K.; Abdel-Hafez, M.F. Non-linear autoregressive delay-dependent INS/GPS navigation system using neural networks. IEEE Sens. J. 2016, 17, 1105–1115. [Google Scholar] [CrossRef]
  32. Boudjedir, H.; Bouhali, O.; Rizoug, N. Adaptive neural network control based on neural observer for quadrotor unmanned aerial vehicle. Adv. Robot. 2014, 28, 1151–1164. [Google Scholar] [CrossRef]
  33. Zhao, B.; Xu, S.; Guo, J.; Jiang, R.; Zhou, J. Integrated strapdown missile guidance and control based on neural network disturbance observer. Aerosp. Sci. Technol. 2019, 84, 170–181. [Google Scholar] [CrossRef]
  34. Lakhal, A.; Tlili, A.; Braiek, N.B. Neural network observer for nonlinear systems application to induction motors. Int. J. Control Autom. 2010, 3, 1–16. [Google Scholar]
  35. Ullah, M.; Zhao, C.; Maqsood, H.; Humayun, M.; Hassan, M.U. Exponential Sliding Mode Control Based on a Neural Network and Finite-Time Disturbance Observer for an Autonomous Aerial Vehicle Exposed to Environmental Disturbances and Parametric Uncertainties. J. Control. Autom. Electr. Syst. 2022, 33, 1659–1670. [Google Scholar] [CrossRef]
  36. Al-Sharman, M.K.; Zweiri, Y.; Jaradat, M.A.K.; Al-Husari, R.; Gan, D.; Seneviratne, L.D. Deep-learning-based neural network training for state estimation enhancement: Application to attitude estimation. IEEE Trans. Instrum. Meas. 2019, 69, 24–34. [Google Scholar] [CrossRef]
  37. Yuning, J.; Rasool, M.A.U.; Qian, B.; Farid, G.; Chaudary, S.T. An Adaptive Neural Network State Estimator for Quadrotor Unmanned Air Vehicle. Int. J. Adv. Comput. Sci. Appl. 2019, 10, 316–321. [Google Scholar] [CrossRef]
  38. Gao, Y.; Zhu, G.; Zhao, T. Based on backpropagation neural network and adaptive linear active disturbance rejection control for attitude of a quadrotor carrying a load. Appl. Sci. 2022, 12, 12698. [Google Scholar] [CrossRef]
  39. Shen, S.; Xu, J.; Chen, P.; Xia, Q. Adaptive Neural Network Extended State Observer-Based Finite-Time Convergent Sliding Mode Control for a Quad Tiltrotor UAV. IEEE Trans. Aerosp. Electron. Syst. 2023, 59, 6360–6373. [Google Scholar] [CrossRef]
  40. Munoz, F.; Espinoza, E.S.; Salazar, S.; Gonzalez, I.; Lozano, R. Second Order Sliding Mode Controllers for Altitude Control of a Quadrotor UAS: Real-Time Implementation in Outdoor Environments. Neurocomputing 2017, 233, 61–71. [Google Scholar] [CrossRef]
  41. Antonio-Toledo, M.E.; Sanchez, E.N.; Alanis, A.Y.; Flórez, J.A.; Perez-Cisneros, M.A. Real-Time Integral Backstepping with Sliding Mode Control for a Quadrotor UAV. IFAC-PapersOnLine 2018, 51, 549–554. [Google Scholar] [CrossRef]
  42. Kahouadji, M.; Mokhtari, M.R.; Choukchou-Braham, A.; Cherki, B. Real-time attitude control of 3 DOF quadrotor UAV using modified Super Twisting algorithm. J. Frankl. Inst. 2020, 357, 2681–2695. [Google Scholar] [CrossRef]
  43. Firdaus, A.R.; Tokhi, M.O. Real-time embedded system of super twisting-based integral sliding mode control for quadcopter UAV. In Proceedings of the 2019 2nd International Conference on Applied Engineering (ICAE), Batam, Indonesia, 2–3 October 2019; IEEE: Piscataway, NJ, USA, 2019; pp. 1–7. [Google Scholar]
  44. Noordin, A.; Mohd Basri, M.A.; Mohamed, Z. Real-Time Implementation of an Adaptive PID Controller for the Quadrotor MAV Embedded Flight Control System. Aerospace 2023, 10, 59. [Google Scholar] [CrossRef]
  45. Borbolla-Burillo, P.; Sotelo, D.; Frye, M.; Garza-Castañón, L.E.; Juárez-Moreno, L.; Sotelo, C. Design and Real-Time Implementation of a Cascaded Model Predictive Control Architecture for Unmanned Aerial Vehicles. Mathematics 2024, 12, 739. [Google Scholar] [CrossRef]
  46. Athayde, A.; Moutinho, A.; Azinheira, J.R. Experimental Nonlinear and Incremental Control Stabilization of a Tail-Sitter UAV with Hardware-in-the-Loop Validation. Robotics 2024, 13, 51. [Google Scholar] [CrossRef]
  47. Quan, Q. Introduction to Multicopter Design and Control, 1st ed.; Springer: Singapore, 2017. [Google Scholar] [CrossRef]
  48. Dai, X.; Ke, C.; Quan, Q.; Cai, K.Y. RFlySim: Automatic test platform for UAV autopilot systems with FPGA-based hardware-in-the-loop simulations. Aerosp. Sci. Technol. 2021, 114, 106727. [Google Scholar] [CrossRef]
  49. Hamayun, M.T.; Ijaz, S.; Bajodah, A.H. Output integral sliding mode fault tolerant control scheme for LPV plants by incorporating control allocation. IET Control Theory Appl. 2017, 11, 1959–1967. [Google Scholar] [CrossRef]
  50. Khalil, H. Nonlinear Systems, 3rd ed.; Michigan State University: East Lansing, MI, USA, 2002. [Google Scholar]
  51. Rahaman-Noronha, M.; William, B.; Bramesfeld, G. Circular Flight-Path Optimization for a Solar-Powered UAV Flying in Horizontal Winds. In Proceedings of the AIAA SCITECH 2023 Forum, National Harbor, MD, USA, 23–27 January 2023. [Google Scholar]
  52. Giagkos, A.; Tuci, E.; Wilson, M.S.; Charlesworth, P.B. UAV flight coordination for communication networks: Genetic algorithms versus game theory. Soft Comput. 2021, 25, 9483–9503. [Google Scholar] [CrossRef]
  53. Baca, T.; Stepan, P.; Spurny, V.; Hert, D.; Penicka, R.; Saska, M.; Thomas, J.; Loianno, G.; Kumar, V. Autonomous landing on a moving vehicle with an unmanned aerial vehicle. J. Field Robot. 2019, 36, 874–891. [Google Scholar] [CrossRef]
  54. Hoblit, F.M. Gust Loads on Aircraft: Concepts and Applications; American Institute of Aeronautics and Astronautics: Washington, DC, USA, 1988. [Google Scholar]
Figure 1. Structure of F450 quadrotor UAV.
Figure 1. Structure of F450 quadrotor UAV.
Aerospace 11 00788 g001
Figure 2. An MLNN structure.
Figure 2. An MLNN structure.
Aerospace 11 00788 g002
Figure 3. NDI-based fast terminal sliding mode controller design.
Figure 3. NDI-based fast terminal sliding mode controller design.
Aerospace 11 00788 g003
Figure 4. Overall observer-based controller framework of the proposed scheme.
Figure 4. Overall observer-based controller framework of the proposed scheme.
Aerospace 11 00788 g004
Figure 5. Tracking response under nominal condition without disturbance and noise.
Figure 5. Tracking response under nominal condition without disturbance and noise.
Aerospace 11 00788 g005
Figure 6. Control inputs under nominal condition.
Figure 6. Control inputs under nominal condition.
Aerospace 11 00788 g006
Figure 7. States tracking response with parameter variations, disturbance and noise.
Figure 7. States tracking response with parameter variations, disturbance and noise.
Aerospace 11 00788 g007
Figure 8. Control inputs with parameter variations, disturbance and noise.
Figure 8. Control inputs with parameter variations, disturbance and noise.
Aerospace 11 00788 g008
Figure 9. Hardware-in-Loop architecture.
Figure 9. Hardware-in-Loop architecture.
Aerospace 11 00788 g009
Figure 10. HIL setup.
Figure 10. HIL setup.
Aerospace 11 00788 g010
Figure 11. HIL simulation altitude and attitude tracking.
Figure 11. HIL simulation altitude and attitude tracking.
Aerospace 11 00788 g011
Figure 12. HIL simulations of position tracking.
Figure 12. HIL simulations of position tracking.
Aerospace 11 00788 g012
Figure 13. Control inputs HILS.
Figure 13. Control inputs HILS.
Aerospace 11 00788 g013
Figure 14. F450 quadrotor for real-time implementation.
Figure 14. F450 quadrotor for real-time implementation.
Aerospace 11 00788 g014
Figure 15. Attitude angle tracking on real quadrotor system with wind disturbance effect.
Figure 15. Attitude angle tracking on real quadrotor system with wind disturbance effect.
Aerospace 11 00788 g015
Table 1. Parameter values.
Table 1. Parameter values.
ParametersMagnitudesDetails
I p x 2.11 × 10 2 kgm2Inertia along X-axis
I p y 2.19 × 10 2 kgm2Inertia along Y-axis
I p z 3.66 × 10 2 kgm2Inertia along Z-axis
I p r 1.28 × 10 4 kgm2Inertia of Rotor
b1.105 × 10 5 Ns2Thrust coefficient
d7.3 × 10 2 ms2Coefficient of drag
m1.4 kgMass
l2.25 × 10 1 mArm length
Table 2. Parameter values of controller/observer.
Table 2. Parameter values of controller/observer.
Gain k p ϕ k p θ k p ψ k p x k p y k p z
Value552.510105
Gain k i ϕ k i θ k i ψ k i x k i y k i z
Value3.753.752.53.53.255
Gain κ f ϕ κ f θ κ f ψ κ f x κ f y κ f z
Value331.52.52.54
Gain η f ϕ η f θ η f ψ η f x η f y κ f z
Value0.50.510.250.250.01
Gain γ β Υ λ f x λ f y λ f z
Value2382050.0080.080.080.05
Table 3. Comaprison of mean square error values for nominal case.
Table 3. Comaprison of mean square error values for nominal case.
ParametersProposed Scheme[39]
e M S E x 3.3 × 10 6 6.2 × 10 6
e M S E y 2.5 × 10 6 4.3 × 10 6
e M S E z 2.01 × 10 6 3.3 × 10 6
e M S E ϕ 1.3 × 10 5 3.9 × 10 5
e M S E θ 2.13 × 10 5 4.4 × 10 5
e M S E ψ 3.4 × 10 5 5.1 × 10 5
Table 4. MSE values under noise, disturbance, and parameter variations.
Table 4. MSE values under noise, disturbance, and parameter variations.
ParametersCase ACase BCase CCase C [44]
e M S E x 2.5 × 10 5 6.1 × 10 5 3.1 × 10 4 4.9 × 10 3
e M S E y 1.63 × 10 5 5.7 × 10 5 2.4 × 10 4 3.8 × 10 3
e M S E z 1.1 × 10 5 2.9 × 10 5 2.1 × 10 4 8.3 × 10 3
e M S E ϕ 8.3 × 10 5 2.3 × 10 4 5.2 × 10 4 7.3 × 10 3
e M S E θ 9.13 × 10 5 3.5 × 10 4 9.8 × 10 3 4.13 × 10 2
e M S E ψ 2.35 × 10 4 4.2 × 10 4 6.7 × 10 4 3.4 × 10 3
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

Akhtar, Z.; Naqvi, S.A.Z.; Khan, Y.A.; Hamayun, M.T.; Ijaz, S. Design and Experimental Validation of an Adaptive Multi-Layer Neural Network Observer-Based Fast Terminal Sliding Mode Control for Quadrotor System. Aerospace 2024, 11, 788. https://doi.org/10.3390/aerospace11100788

AMA Style

Akhtar Z, Naqvi SAZ, Khan YA, Hamayun MT, Ijaz S. Design and Experimental Validation of an Adaptive Multi-Layer Neural Network Observer-Based Fast Terminal Sliding Mode Control for Quadrotor System. Aerospace. 2024; 11(10):788. https://doi.org/10.3390/aerospace11100788

Chicago/Turabian Style

Akhtar, Zainab, Syed Abbas Zilqurnain Naqvi, Yasir Ali Khan, Mirza Tariq Hamayun, and Salman Ijaz. 2024. "Design and Experimental Validation of an Adaptive Multi-Layer Neural Network Observer-Based Fast Terminal Sliding Mode Control for Quadrotor System" Aerospace 11, no. 10: 788. https://doi.org/10.3390/aerospace11100788

APA Style

Akhtar, Z., Naqvi, S. A. Z., Khan, Y. A., Hamayun, M. T., & Ijaz, S. (2024). Design and Experimental Validation of an Adaptive Multi-Layer Neural Network Observer-Based Fast Terminal Sliding Mode Control for Quadrotor System. Aerospace, 11(10), 788. https://doi.org/10.3390/aerospace11100788

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