Open AccessArticle

Adaptive Sliding Mode Attitude Control of Quadrotor UAVs Based on the Delta Operator Framework

CICAEET, School of Automation, Nanjing University of Information Science and Technology, Nanjing 210044, China

Author to whom correspondence should be addressed.

Symmetry 2022, 14(3), 498; https://doi.org/10.3390/sym14030498

Submission received: 27 January 2022 / Revised: 23 February 2022 / Accepted: 25 February 2022 / Published: 28 February 2022

(This article belongs to the Topic Dynamical Systems: Theory and Applications)

Download

Browse Figures

Versions Notes

Abstract

In this paper, a novel adaptive sliding-mode control algorithm is proposed for the attitude control of quadrotor unmanned aerial vehicles (UAVs) under the delta operator framework. First, the delta operator technique is used to discretize the attitude control systems of a quadrotor UAV. Then, based on the linear matrix inequality technique, a linear sliding surface is designed to ensure the asymptotical stability of the quadrotor UAV attitude control system during the sliding motion process. Second, by the estimated external disturbance using a radical basis function (RBF) neural network, an adaptive sliding-mode attitude controller is designed such that the states of the quadrotor UAV attitude systems can be driven towards the desired sliding surface, and thus the attitude control objective of the qudarotor UAV is achieved. Compared with the traditional adaptive sliding-mode control algorithm, the proposed adaptive sliding-mode control algorithm can effectively realize the attitude control of a quadrotor UAV subject to strong disturbances and couplings. Finally, comparisons of the simulation results verify the effectiveness and superiority of the control algorithm proposed in this paper.

Keywords:

quadrotor UAV; attitude control; sliding mode control; delta operator

1. Introduction

Quadrotor unmanned aerial vehicles (UAVs) are multi-rotor aircraft that can take off, land and hover freely. Due to the high flexibility, low cost and strong concealment, they are widely used in various fields, such as rescue and disaster relief, aerial photography mapping, agricultural plant protection, inspection etc. A quadrotor UAV is a typical multi-input multi-output, under-actuated and strongly coupled nonlinear system. In the process of flight, attitude transformation and flight stability are the top priorities; therefore, it is particularly important to study the flight attitude control [1]. The research on the flight control of quadrotor UAVs has become a hotspot in recent years.

Currently, all kinds of control methods, including proportional-integral-differential (PID) control [2], active disturbance rejection control [3,4], sliding mode control [5,6,7,8] and neural network control [9], have been utilized in the attitude control of quadrotor UAVs. The article [2] shows the property of stability and robustness for a nonlinear implicit PID control algorithm of finite-time stabilization of quadrotor UAVs subject to constraints bounded external disturbance. Active disturbance rejection control comes from the improvement of PID control [3].

This is mainly composed of a tracking differentiator, state observer and nonlinear state error feedback control law. In [4], by the combination of adaptive control technique, a linear active disturbance rejection control strategy was developed to address the attitude control problem for the quadrotor UAV systems. Sliding mode control, as a well-known robust control method [5,6], has also been used in quadrotor flight control.

For instance, in [7], a discrete-time sliding-mode control algorithm was proposed for solving the position and attitude tracking control of a small quadrotor UAV. In [8], a sliding-mode attitude tracking control of the quadrotor UAVs with time-varying mass was designed. As an advanced intelligent control algorithm, an RBF neural network has a powerful learning ability, can learn complex uncertain models online and can deal with highly nonlinear control problems. The authors of [9] proposed a fault-tolerant control algorithm for quadrotors, which effectively enhanced the robustness of the system by combining the advantages of non-singular terminal sliding mode control and neural network.

To make full use of the advantages of sliding mode control and adaptive control methods, adaptive sliding mode control design of the quadrotor have also been investigated to improve robustness and adaptability in the literature. In [10], on the basis of Udwadia–Kalaba theory, the adaptive robust tracking control problem is considered for quadrotor UAVs with mismatched uncertainties.

In [11,12], adaptive fuzzy global sliding mode control methods are presented for quandrotor UAVs subject to parameter uncertainties and control chattering problems. A finite-time adaptive integral backstepping sliding mode control design method is introduced in [13], where adaptive control is used to compensate the unknown disturbance upper bound, and semi-global practical stability of the fight attitude control of the quadrotor is achieved.

It is worth noting that most of the results on adaptive sliding mode control of the quadrotor are based on the continuous-time controller, and the controller structure is relative complex. In fact, it was said in [14] that a discrete-time control scheme is much more suitable in practicable situations since the utilization of computer and network in the quadrotor, which is a challenging task has not been well solved so far.

On the other hand, signal discrete sampling in modern control theory is essential due to the wide application of digital controllers. Middleton and Goodwin introduced the delta operator discretization method to the modern control fields [15,16]. Compared with the discretization technique based on traditional shift operator method, the delta operator method is more suitable for high-speed sampling [17,18]. The authors in [19] studied the robust fault-tolerant stabilization of uncertain switched systems under delta operator framework.

In [20,21], the insensitive robust output tracking controller design and

H_{\infty}

filter design for discrete-time systems were studied in a unified delta operator approach framework. The results showed that the delta-domain model was better than the standard shift-domain model in avoiding the inherent numerical ill-conditions at high sampling rates. Recently, using the delta operator approach, the consensus problem of muli-agent systems are well investigated in [22,23]. It is worth noticing that few results have been published associated with the control design of quadrotor UAVs under a delta operator framework, which is a high-speed discrete sampling system.

Inspired by the above research works, this paper is concerned with adaptive sliding-mode attitude control design problem for quadrotor UAVs based on delta operator framework. The main contribution of this paper is summarized as follows. First, the delta operator technique is used to discretize the attitude system of the quadrotor UAV. Then, the designs of linear sliding surface and an adaptive sliding mode reaching control law are shown in the delta domain. The design of the linear sliding surface is on the basis of linear matrix inequality technique, and it ensures the asymptotic stability of the quadrotor UAVs on sliding motion.

By the estimated external disturbance using radial basis function (RBF) neural network, the adaptive sliding mode controller is designed for ensuring all attitudes of the quadrotor UAVs can be driven to the proposed sliding surface and thus the attitude control is achieved. Finally, simulation result comparisons demonstrate the effectiveness and superiority of the proposed attitude control algorithm.

The remainder of this paper is organized as follows. In Section 2, a system model of quadrotor UAVs is presented. In Section 3, a sliding-mode attitude control design algorithm is presented, which includes the design of a sliding surface and the design of an adaptive sliding mode controller. In Section 4, an illustrative simulation example is provided to demonstrate the effectiveness and superiority of the proposed method. Our conclusions are given in Section 5.

2. System Model Description of Quadrotor UAV

The quadrotor UAV has a simple appearance and can take off and land vertically with high control flexibility. It consists of four electrical motors with propellers, which are attached to a rigid cross frame. Supposing the quadrotor UAV has a symmetrical structure, the drag and the thrust forces are proportional to the square of the rotors speed, and ignoring blade flapping and gyro effect, the linearized quadrotor UAV model is obtained from [24,25,26]:

\begin{matrix} \{\begin{matrix} J_{1} \ddot{ϕ} (t) = - K_{1} l \dot{ϕ} (t) + l u_{1} (t) + f_{1} (t), \\ J_{2} \ddot{θ} (t) = - K_{2} l \dot{θ} (t) + l u_{2} (t) + f_{2} (t), \\ J_{3} \ddot{ψ} (t) = - K_{3} \dot{ψ} (t) + c u_{3} (t) + f_{3} (t), \end{matrix} \end{matrix}

(1)

where

ϕ (t)

θ (t)

and

ψ (t)

are the roll angle, pitch angle and yaw angle of the quadrotor UAV, respectively.

J_{1}

J_{2}

and

J_{3}

represent the moments of inertia of each axis of the UAV,

K_{1}

K_{2}

and

K_{3}

are the drag coefficients and are positive constants, l represents the distance from the center of each rotor to the center of the body, c is the proportional coefficient of force and torqu,

u_{i} (t)

represent the control input to the quadrotor UAVs with

u_{1} (t) = F_{1} (t) - F_{2} (t) - F_{3} (t) + F_{4} (t)

u_{2} (t) = F_{1} (t) + F_{2} (t) - F_{3} (t) - F_{4} (t)

and

u_{3} (t) = F_{1} (t) - F_{2} (t) + F_{3} (t) - F_{4} (t)

F_{i} (t)

represents the lift of propeller. In addition,

f_{1} (t)

f_{2} (t)

and

f_{3} (t)

represent unknown external disturbances and coupling terms.

Remark 1.

If the unknown external disturbances and couplings are ignored, the attitude control model of the quadrotor UAVs can be described as follows.

\begin{matrix} \{\begin{matrix} J_{1} \ddot{ϕ} (t) = - K_{1} l \dot{ϕ} (t) + l u_{1} (t), \\ J_{2} \ddot{θ} (t) = - K_{2} l \dot{θ} (t) + l u_{2} (t), \\ J_{3} \ddot{ψ} (t) = - K_{3} \dot{ψ} (t) + c u_{3} (t), \end{matrix} \end{matrix}

(2)

Research works based on model (2) can be seen in [27,28,29,30,31]. Clearly, the dynamic model (1) in this paper is more suitable in practical situations than (2).

Let

x_{1} (t) = [\begin{matrix} ϕ (t) \\ \dot{ϕ} (t) \end{matrix}], x_{2} (t) = [\begin{matrix} θ (t) \\ \dot{θ} (t) \end{matrix}], x_{3} (t) = [\begin{matrix} ψ (t) \\ \dot{ψ} (t) \end{matrix}]

{\bar{A}}_{1} = [\begin{matrix} 0 & 1 \\ 0 & \frac{- K_{1} l}{J_{1}} \end{matrix}]

{\bar{A}}_{2} = [\begin{matrix} 0 & 1 \\ 0 & \frac{- K_{2} l}{J_{2}} \end{matrix}]

{\bar{A}}_{3} = [\begin{matrix} 0 & 1 \\ 0 & \frac{- K_{3}}{J_{3}} \end{matrix}]

{\bar{B}}_{1} = [\begin{matrix} 0 \\ \frac{l}{J_{1}} \end{matrix}]

{\bar{B}}_{2} = [\begin{matrix} 0 \\ \frac{l}{J_{2}} \end{matrix}]

and

{\bar{B}}_{3} = [\begin{matrix} 0 \\ \frac{c}{J_{3}} \end{matrix}]

, the linearized quadrotor UAVs can be rewritten as

\begin{matrix} {\dot{x}}_{i} (t) = {\bar{A}}_{i} x_{i} (t) + {\bar{B}}_{i} u_{i} (t) + {\bar{B}}_{i} f_{i} (t), \end{matrix}

(3)

where

i = 1, 2, 3

x_{i} (t) \in R^{2}

is the state variable of the ith attitude control system,

{\bar{A}}_{i} \in R^{2 \times 2}

and

{\bar{B}}_{i} \in R^{2 \times 1}

are the system matrix and control input matrix, respectively.

Applying the delta operator technique to discretize the system (3), one can find the quadrotor UAVs in the delta domain

\begin{matrix} δ x_{i} (t_{k}) = A_{i} x_{i} (t_{k}) + B_{i} u_{i} (t_{k}) + B_{i} f_{i} (t_{k}), \end{matrix}

(4)

where

δ x_{i} (t_{k})

is the delta operator calculation for the state variable

x_{i} (t_{k})

. It is defined mathematically as

\begin{matrix} δ x_{i} (t_{k}) = \{\begin{matrix} \frac{d}{d t} x_{i} (t), & T = 0, \\ \frac{x_{i} (t_{k} + T) - x_{i} (t_{k})}{T}, & T \neq 0, \end{matrix} \end{matrix}

(5)

where

t_{k} = k T

, T is the sampling period. Clearly, when

T \to 0

, the delta operator

δ x_{i} (t_{k})

is a derivative operation

\dot{x} (t)

and can be used to describe a continuous time dynamic system. When

T = 1

, the delta operator

δ x_{i} (t_{k})

is a difference operation of traditional shift operator, which can be used to describe a traditional discrete time dynamical system. In terms of the definition of the delta operator, one can see that

A_{i} = (e^{{\bar{A}}_{i} T} - I) / T

B_{i} = (\int_{0}^{T} e^{{\bar{A}}_{i} T} {\bar{B}}_{i} d t) / T

i = 1, 2, 3

For the control design of the quadrotor UAVs

(4)

under the delta operator framework, it is assumed that the pair

(A_{i}, B_{i})

is controllable and that column

B_{i}

is full rank,

i = 1, 2, 3

To facilitate the proof, the relevant lemmas are introduced. Lemma 1 shows the definition of the asymptotical stability in the delta domain, Lemma 2 and Lemma 3 will be used to the theoretical proof of Theorem 1, where Lemma 2 is used to the mathematical calculation of delta operator, and Lemma 3 is used as a technique to convert a nonlinear linear matrix inequality into a linear matrix inequality.

Lemma 1

([32]). For the delta operator system

\begin{matrix} δ (x (t_{k})) = g (x (t_{k})), \end{matrix}

(6)

where

g (x (t_{k}))

with

g (0) = 0

, is a linear or nonlinear function, if there exists a positive definite function

V (x (t_{k}))

in the delta domain, such that for any state

x (t_{k})

satisfies

\begin{matrix} δ V (x (t_{k})) = \frac{V (x (t_{k} + T)) - V (x (t_{k}))}{T} < 0, \end{matrix}

(7)

then the system (6) is asymptotical stable in the delta domain.

Lemma 2

([16]). (Differentiation of Product). For any time functions

x (t_{k})

and

y (t_{k})

, we have the following fact.

\begin{matrix} δ [x (t_{k}) y (t_{k})] = (δ x (t_{k})) y (t_{k}) + x (t_{k}) (δ y (t_{k})) + T δ x (t_{k}) δ y (t_{k}) . \end{matrix}

(8)

Lemma 3

([33]). (Schur–Complement Lemma) For a given symmetric matrix

S = [\begin{matrix} S_{11} & S_{12} \\ S_{21} & S_{22} \end{matrix}]

, where

S_{11} = S_{11}^{T}

S_{22} = S_{22}^{T}

S_{12} = S_{21}^{T}

, the inequalities (9)–(11) are equivalent to each other,

\begin{matrix} \begin{matrix} S < 0, \end{matrix} \end{matrix}

(9)

\begin{matrix} \begin{matrix} S_{11} < 0, & S_{22} - S_{21} S_{11}^{- 1} S_{12} < 0, \end{matrix} \end{matrix}

(10)

\begin{matrix} \begin{matrix} S_{22} < 0, & S_{11} - S_{12} S_{22}^{- 1} S_{21} < 0 . \end{matrix} \end{matrix}

(11)

3. Design of the Sliding Mode Attitude Control Algorithm

The design of the sliding-mode attitude control algorithm of the quadrotor UAVs is divided into two steps. The first step is to design a linear sliding surface

σ_{i} (k) = 0

to ensure the stability of the quadrotor UAVs during sliding motion process. The second step is to design an adaptive sliding mode attitude controller such that all attitudes of the quadrotor UAVs can be driven to the designed sliding surface

σ_{i} (k) = 0

3.1. Sliding Surface Design

Let

x_{i} (t_{k}) = [\begin{matrix} x_{i 1} (t_{k}) \\ x_{i 2} (t_{k}) \end{matrix}]

A_{i} = [\begin{matrix} A_{i 11} & A_{i 12} \\ A_{i 21} & A_{i 22} \end{matrix}]

and

B_{i} = [\begin{matrix} 0 \\ B_{i 2} \end{matrix}]

, the quadrotor UAV (4) can be rewritten as:

\begin{matrix} δ x_{i 1} (t_{k}) & = A_{i 11} x_{i 1} (t_{k}) + A_{i 12} x_{i 2} (t_{k}), \end{matrix}

(12)

\begin{matrix} δ x_{i 2} (t_{k}) & = A_{i 21} x_{i 1} (t_{k}) + A_{i 22} x_{i 2} (t_{k}) + B_{i 2} u_{i} (t_{k}) + B_{i 2} f_{i} (t_{k}) . \end{matrix}

(13)

The linear sliding surface is designed as

\begin{matrix} σ_{i} (t_{k}) = S_{i} x_{i} (t_{k}) = 0, \end{matrix}

(14)

where

i = 1, 2, 3

S_{i} = [\begin{matrix} Z_{i} & 1 \end{matrix}]

is the sliding vector with

Z_{i} \in R

to be determined. Substituting

σ_{i} (t_{k}) = Z_{i} x_{i 1} (t_{k}) + x_{i 2} (t_{k})

into (12), the sliding dynamics, i.e., the reduced order quadrotor UAVs system

δ x_{i 1} (t_{k})

is described as

\begin{matrix} δ x_{i 1} (t_{k}) = (A_{i 11} - A_{i 12} Z_{i}) x_{i 1} (t_{k}) . \end{matrix}

(15)

To make the sliding motion on (15) be asymptotical stable, it is necessary to design the appropriate parameter

Z_{i}

of the sliding surface.

Theorem 1.

If there exist positive real parameters

X_{i}

and

W_{i}

and a real number

Y_{i}

, such that the following linear matrix inequality

\begin{matrix} [\begin{matrix} - 2 W_{i} & A_{i 11} X_{i} - A_{i 12} Y_{i} & T W_{i} \\ * & H e (A_{i 11} X_{i} - A_{i 12} Y_{i}) & 0 \\ * & * & - T X_{i} \end{matrix}] < 0 \end{matrix}

(16)

is satisfied for a given sampling period

T > 0

, then the reduced order quadrotor UAVs system (15) is asymptotical stable, where

H e (A_{i 11} X_{i} - A_{i 12} Y_{i}) = A_{i 11} X_{i} - A_{i 12} Y_{i} + {(A_{i 11} X_{i} - A_{i 12} Y_{i})}^{T} .

The symbol ∗ is used to represent a term that is induced by symmetry. Furthermore, the sliding surface is designed as

\begin{matrix} σ_{i} (t_{k}) = [\begin{matrix} Y_{i} X_{i}^{- 1} & 1 \end{matrix}] x_{i} (t_{k}) = 0 . \end{matrix}

(17)

Proof of Theorem 1.

Let

A_{i 1} = A_{i 11} - A_{i 12} Z_{i}

, one can see that the sliding motion obeys

δ x_{i 1} (t_{k}) = A_{i 1} x_{i 1} (t_{k})

. Now, the Lyapunov function candidated is chosen as

\begin{matrix} V_{1} (t_{k}) = x_{i 1}^{T} (t_{k}) P_{i} x_{i 1} (t_{k}), \end{matrix}

(18)

where

P_{i}

is a positive definite symmetric matrix.

In terms of Lemma 2, one has

\begin{matrix} δ V_{1} (t_{k}) = 2 x_{i 1}^{T} (t_{k}) P_{i} δ x_{i 1} (t_{k}) + T {(δ x_{i 1} (t_{k}))}^{T} P_{i} δ x_{i 1} (t_{k}) . \end{matrix}

(19)

Furthermore, let us introduce

\begin{matrix} 0 = - 2 {(δ x_{i 1} (t_{k}))}^{T} {\hat{W}}_{i} [δ x_{i 1} (t_{k}) - A_{i 1} x_{i 1} (t_{k})], \end{matrix}

(20)

where

{\hat{W}}_{i}

is a positive definite matrix, and by the combination of the sliding motion Equations (15) and (20), the state trajectory of (19) along the reduced-order system (15) becomes

\begin{matrix} δ V_{1} (t_{k}) = - {[\begin{matrix} δ x_{i 1} (t_{k}) \\ x_{i 1} (t_{k}) \end{matrix}]}^{T} Θ_{i} [\begin{matrix} δ x_{i 1} (t_{k}) \\ x_{i 1} (t_{k}) \end{matrix}], \end{matrix}

(21)

where

\begin{matrix} Θ_{i} = - [\begin{matrix} T P_{i} - {\hat{W}}_{i} - {\hat{W}}_{i}^{T} & {\hat{W}}_{i} A_{i 1} \\ * & P_{i} A_{i 1} + A_{i 1}^{T} P_{i} \end{matrix}] . \end{matrix}

(22)

It is clear that, when

Θ_{i} > 0

, the reduced-order system of the quadrotor UAV (15) is asymptotically stable by Lemma 1. Now, multiplying the left and right sides of Equation (22) with

diag [{\hat{W}}_{i}^{- 1}, {P_{i}}^{- 1}]

, we have

\begin{matrix} [\begin{matrix} T {\hat{W}}_{i}^{- T} P {\hat{W}}_{i}^{- 1} - 2 {\hat{W}}_{i}^{- 1} & A_{i 1} {P_{i}}^{- 1} \\ * & {P_{i}}^{- 1} A_{i 1}^{T} + A_{i 1} {P_{i}}^{- 1} \end{matrix}] < 0 . \end{matrix}

(23)

Using the Schure–Complement Lemma, one can obtain the condition (16) for asymptotically stable of the sliding motion, where

X_{i} = {P_{i}}^{- 1}

W_{i} = {\hat{W}}_{i}^{- 1}

Y_{i} = Z_{i} X_{i}

. The proof of Theorem 1 is therefore completed. □

3.2. Adaptive Sliding Mode Controller Design

To design the adaptive sliding mode controller, a RBF neural network is first used online to approximate the unknown nonlinear external disturbance

f_{i} (t_{k})

. The RBF neural network usually consists of three parts: Input layer, hidden layer and output layer. How to configure the RBF neural network has been well investigated in the existing works, please see [34] for details. In this paper, the RBF neural network with

2 - 5 - 1

structure is adopted and shown in Figure 1.

The input and output algorithm of a single RBF network is

\begin{matrix} h_{j} (υ) = exp (- \frac{{∥υ - c_{j}∥}^{2}}{2 {b_{j}}^{2}}), \end{matrix}

(24)

\begin{matrix} f (t_{k}) = {(W^{*})}^{T} h_{j} (υ) + η, \end{matrix}

(25)

where

υ

is the input of the neural network, j is the jth node of the hidden layer,

c_{j}

is the center vector of the hidden layer node,

b_{j}

is the width vector of the hidden layer neurons,

W^{*}

is the network approximation ideal weight,

η \leq η_{max}

is the network approximation error.

For the attitude control system of quadrotor UAVs (4), the neural network input is taken as

υ_{i} = {[\begin{matrix} x_{i 1} (t_{k}) & x_{i 2} (t_{k}) \end{matrix}]}^{T}

, and the estimated output of the RBF neural network for the unknown nonlinear function

f_{i} (t_{k})

is:

\begin{matrix} {\overset{⌢}{f}}_{i} (t_{k}) = {\overset{⌢}{W}}_{i}^{T} h_{i j} (υ_{i}), \end{matrix}

(26)

where

{\overset{⌢}{W}}_{i}

is the estimated weight of the network,

i = 1, 2, 3

. The adaptive adjustment law of the estimated weight of the network is designed in the delta domain as follows.

\begin{matrix} δ {\overset{⌢}{W}}_{i} (t_{k}) = γ_{i} σ_{i} (t_{k}) S_{i} B_{i} h_{i j} (υ_{i}) . \end{matrix}

(27)

The designed adaptive sliding mode controller is shown as follows.

\begin{matrix} u_{i} (t_{k}) = - {(S_{i} B_{i})}^{- 1} [S_{i} A_{i} x_{i} (t_{k}) + S_{i} B_{i} {\overset{⌢}{f}}_{i} (t_{k}) + ε_{i} sign (σ_{i} (t_{k})) + β_{i} σ_{i} (t_{k})], \end{matrix}

(28)

where

ε_{i}

and

β_{i}

are adjustable positive parameters,

i = 1, 2, 3

Theorem 2.

Considering the quadrotor UAVs (4) under the delta operator framework, the adaptive sliding-mode attitude controller is designed as shown in (27) and (28), all the attitudes of the quadrotor UAVs can be driven towards the designed sliding surface

σ_{i} (t_{k}) = 0

shown in (14), and the asymptotical stablity is achieved.

Proof of Theorem 2.

Take the Lyapunov function candidated in the delta domain as

\begin{matrix} V_{2} (t_{k}) = \frac{{σ_{i}}^{T} (t_{k}) σ_{i} (t_{k})}{2} + \frac{1}{2 γ_{i}} {\tilde{W}}_{i}^{T} (t_{k}) {\tilde{W}}_{i} (t_{k}), \end{matrix}

(29)

where

γ_{i} > 0

is the network tuning parameter, and

{\tilde{W}}_{i} (t_{k}) = {\overset{⌢}{W}}_{i} - {W_{i}}^{*}

i = 1, 2, 3

. By Lemma 2, one can see that

\begin{matrix} δ V_{2} (t_{k}) = & σ_{i} (t_{k}) δ σ_{i} (t_{k}) + \frac{1}{2 γ_{i}} {(δ {\tilde{W}}_{i} (t_{k}))}^{T} {\tilde{W}}_{i} (t_{k}) + \frac{1}{2 γ_{i}} {\tilde{W}}_{i}^{T} (t_{k}) δ {\tilde{W}}_{i} (t_{k}) + \frac{T}{2} {(∥δ σ_{i} (t_{k})∥)}^{2} \\ + \frac{T}{2 γ} {∥δ {\tilde{W}}_{i} (t_{k}))∥}^{2} \\ = & σ_{i} (t_{k}) δ σ_{i} (t_{k}) + \frac{1}{γ_{i}} {\tilde{W}}_{i}^{T} (t_{k}) (δ {\overset{⌢}{W}}_{i} (t_{k})) + \prod_{i} (t_{k}), \end{matrix}

(30)

where

\prod_{i} (t_{k}) = \frac{T}{2} {(∥δ σ_{i} (t_{k})∥)}^{2} + \frac{T}{2 γ_{i}} {∥δ {\tilde{W}}_{i} (t_{k}))∥}^{2}

i = 1, 2, 3

On the other hand, noticing that

\begin{matrix} δ σ_{i} (t_{k}) = S_{i} [A_{i} x_{i} (t_{k}) + B_{i} u_{i} (t_{k}) + B_{i} f_{i} (t_{k})] . \end{matrix}

(31)

Substituting the designed adaptive sliding mode controller (28) into (31), one can further obtain

\begin{matrix} δ σ_{i} (t_{k}) = S_{i} B_{i} [f_{i} (t_{k}) - {\overset{⌢}{f}}_{i} (t_{k})] - ε_{i} sign (σ_{i} (t_{k})) - β_{i} σ_{i} (t_{k}) . \end{matrix}

(32)

Noticing, from (25) and (26), that

f_{i} (t_{k}) - {\overset{⌢}{f}}_{i} (t_{k}) = - {\tilde{W}}_{i}^{T} (t_{k}) h_{i j} (υ_{i}) + η_{i}

i = 1, 2, 3

, In terms of (30) and (32), it is not difficult to see that

\begin{matrix} δ V_{2} (t_{k}) = & σ_{i} (k) S_{i} B_{i} [- {\tilde{W}}_{i}^{T} (t_{k}) h_{i j} (υ_{i}) + η_{i}] - β_{i} σ_{i} {(t_{k})}^{2} - ε_{i} σ_{i} (t_{k}) s i g n (σ_{i} (t_{k})) \\ + \frac{1}{γ_{i}} {\tilde{W}}_{i}^{T} (t_{k}) δ {\overset{⌢}{W}}_{i} (t_{k}) + \prod_{i} (t_{k}) \\ \leq & (- σ_{i} (k) (S_{i} B_{i})) {\tilde{W}}_{i}^{T} (t_{k}) h_{i j} (υ_{i}) + |σ_{i} (t_{k})| |S_{i} B_{i}| η_{i} - ε_{i} |σ_{i} (t_{k})| \\ - β_{i} {|σ_{i} (t_{k})|}^{2} + \frac{1}{γ_{i}} {\tilde{W}}_{i}^{T} (t_{k}) δ {\overset{⌢}{W}}_{i} (t_{k}) + \prod_{i} (t_{k}) \\ \leq & {\tilde{W}}_{i}^{T} (t_{k}) [- σ_{i} (t_{k}) (S_{i} B_{i}) h_{i j} (υ_{i}) + \frac{1}{γ_{i}} δ {\overset{⌢}{W}}_{i} (t_{k})] + |σ_{i} (t_{k})| |S_{i} B_{i}| η_{i} \\ - β_{i} {|σ_{i} (t_{k})|}^{2} - ε_{i} |σ_{i} (t_{k})| + \prod_{i} (t_{k}) . \end{matrix}

(33)

Substituting the adaptive law (27) into (33), one can see that

\begin{matrix} δ V_{2} (t_{k}) \leq |σ_{i} (t_{k})| |S_{i} B_{i}| η_{i} - ε_{i} |σ_{i} (t_{k})| - β_{i} {|σ_{i} (t_{k})|}^{2} + \prod_{i} (t_{k}) . \end{matrix}

(34)

Since

\prod_{i} (t_{k})

is reasonably bounded [18] and it is also small when the sampling period T is close to 0; then, as long as

ε_{i}

and

β_{i}

are large enough to satisfy

ε_{i} + β_{i} |σ_{i} (t_{k})| ≫ η_{i} |S_{i} B_{i}|

, it is easy to obtain

δ V_{2} (t_{k}) < 0

. Therefore, under the action of the RBF neural-network-based adaptive law (27) and the sliding mode controller (28), the state of the quadrotor UAVs (4) can reach the designed sliding surface

σ_{i} (t_{k}) = 0

, and asymptotical stability can be finally guaranteed. □

The diagram of the RBF neural-network-based adaptive sliding-mode control algorithm is shown in Figure 2, where the RBF network is used to estimate the unknown external disturbances and couplings. The adaptive sliding mode controller was designed to achieve the flight attitude control using the estimated information.

4. Illustrative Example

In this section, simulation results are shown to verified the effectiveness and superiority of the proposed adaptive sliding mode control design method of quadrotor UAVs. In Section 4.1, parameter settings and model of the quadrotor UAVs are given; and in Section 4.2, the design of the linear surfaces is shown. Simulation results are shown in Section 4.3.

4.1. Parameter Settings and Modelling of the Quadrotor UAV

To verify the effectiveness of adaptive sliding-mode attitude control algorithm for quadrotor UAVs proposed in this paper, simulation experiments were performed. The physical parameters of the quadrotor UAVs are shown in Table 1.

The continuous-time models of the roll angle, the pitch angle and the yaw angle are expressed as

\begin{matrix} [\begin{matrix} \dot{ϕ} (t) \\ \ddot{ϕ} (t) \end{matrix}] = [\begin{matrix} 0 & 1 \\ 0 & - 0.724 \end{matrix}] [\begin{matrix} ϕ (t) \\ \dot{ϕ} (t) \end{matrix}] + [\begin{matrix} 0 \\ 4.789 \end{matrix}] u (t) + [\begin{matrix} 0 \\ 4.789 \end{matrix}] f (t) . \end{matrix}

(35)

\begin{matrix} [\begin{matrix} \dot{θ} (t) \\ \ddot{θ} (t) \end{matrix}] = [\begin{matrix} 0 & 1 \\ 0 & - 0.724 \end{matrix}] [\begin{matrix} θ (t) \\ \dot{θ} (t) \end{matrix}] + [\begin{matrix} 0 \\ 4.789 \end{matrix}] u (t) + [\begin{matrix} 0 \\ 4.789 \end{matrix}] f (t) . \end{matrix}

(36)

\begin{matrix} [\begin{matrix} \dot{ψ} (t) \\ \ddot{ψ} (t) \end{matrix}] = [\begin{matrix} 0 & 1 \\ 0 & - 1.364 \end{matrix}] [\begin{matrix} ψ (t) \\ \dot{ψ} (t) \end{matrix}] + [\begin{matrix} 0 \\ 4.789 \end{matrix}] u (t) + [\begin{matrix} 0 \\ 4.789 \end{matrix}] f (t) . \end{matrix}

(37)

Setting the sampling time T = 0.01 s and taking the state

x_{1} (t) = [\begin{matrix} ϕ (t) \\ \dot{ϕ} (t) \end{matrix}]

, now using the delta operator discretization technique in [16], the roll angle system model in the delta operator framework is established as follows.

\begin{matrix} δ x_{1} (t_{k}) = [\begin{matrix} 0 & 0.9964 \\ 0 & - 0.7214 \end{matrix}] x_{1} (t_{k}) + [\begin{matrix} 0.0477 \\ 4.7545 \end{matrix}] u_{1} (t_{k}) + [\begin{matrix} 0.0477 \\ 4.7545 \end{matrix}] f_{1} (t_{k}) . \end{matrix}

(38)

Similarly, let

x_{2} (t) = [\begin{matrix} θ (t) \\ \dot{θ} (t) \end{matrix}]

and

x_{3} (t) = [\begin{matrix} ψ (t) \\ \dot{ψ} (t) \end{matrix}]

, one can also establish the pitch angle system model after the discretization of the delta operator

\begin{matrix} δ x_{2} (t_{k}) = [\begin{matrix} 0 & 0.9964 \\ 0 & - 0.7214 \end{matrix}] x_{2} (t_{k}) + [\begin{matrix} 0.0477 \\ 4.7545 \end{matrix}] u_{2} (t_{k}) + [\begin{matrix} 0.0477 \\ 4.7545 \end{matrix}] f_{2} (t_{k}) . \end{matrix}

(39)

and the yaw angle system model after discretization of the delta operator

\begin{matrix} δ x_{3} (t_{k}) = [\begin{matrix} 0 & 0.9932 \\ 0 & - 1.3547 \end{matrix}] x_{3} (t_{k}) + [\begin{matrix} 0.0476 \\ 4.7241 \end{matrix}] u_{3} (t_{k}) + [\begin{matrix} 0.0476 \\ 4.7241 \end{matrix}] f_{3} (t_{k}) . \end{matrix}

(40)

4.2. Design of Linear Surfaces

For the quadrotor UAVs (38)–(40), by solving the linear matrix inequality (16) with MATLAB LMI control toolbox, we can find

X_{1} = X_{2} = X_{3} = 1.0133

Y_{1} = Y_{2} = 0.3357

and

Y_{3} = 0.3367

. From

Z_{i} = Y_{i} X_{i}^{- 1}

, the sliding vectors are

S_{1} = S_{2} = [\begin{matrix} 0.3313 & 1 \end{matrix}]

and

S_{3} = [\begin{matrix} 0.3323 & 1 \end{matrix}]

. Thus, the sliding surfaces are designed as

σ_{1} (t_{k}) = [\begin{matrix} 0.3313 & 1 \end{matrix}] x_{1} (t_{k}) = 0

σ_{2} (t_{k}) = [\begin{matrix} 0.3313 & 1 \end{matrix}] x_{2} (t_{k}) = 0

and

σ_{3} (t_{k}) = [\begin{matrix} 0.3323 & 1 \end{matrix}] x_{3} (t_{k}) = 0

4.3. Simulation Retsults

For the dynamical model shown in (38)–(40), the effectiveness and superiority of the proposed control algorithm are verified via simulation comparisons. Comparison results of the traditional adaptive sliding-mode control algorithm (AC+SMC) in [25] and the proposed adaptive sliding-mode control algorithm (RBF+SMC) in this paper are given.

In [25],

f_{i} (t_{k})

is supposed as

∥ f_{i} (t_{k}) ∥ \leq η_{1} + η_{2} ∥ x_{i} (t_{k}) ∥

with

η_{1}

and

η_{2}

are unknown positive constants. Clearly, such an assumption is relatively strict, which requires not only linear growth conditions but also ignores the coupling problem. Applying the design method there, adaptive laws and controller used are designed as in (41)–(43).

\begin{matrix} δ {\hat{η}}_{i 1} (t_{k}) = ∥σ_{i} (t_{k})∥ ∥S_{i} B_{i}∥, \end{matrix}

(41)

\begin{matrix} δ {\hat{η}}_{i 2} (t_{k}) = ∥σ_{i} (t_{k})∥ ∥S_{i} B_{i}∥ ∥x_{i} (t_{k})∥, \end{matrix}

(42)

\begin{matrix} u_{i} (t_{k}) = - {(S_{i} B_{i})}^{- 1} [S_{i} A_{i} x_{i} (t_{k}) + (μ_{i} + {\hat{η}}_{i 1} (t_{k}) ∥S_{i} B_{i}∥ + {\hat{η}}_{i 2} (t_{k}) ∥S_{i} B_{i}∥ ∥x_{i} (t_{k})∥) sign (σ_{i} (t_{k}))], \end{matrix}

(43)

where

{\hat{η}}_{i 1} (t_{k})

and

{\hat{η}}_{i 2} (t_{k})

represent the estimated values of unknown interference coefficients

η_{i 1}

and

η_{i 2}

, respectively.

Taking the initial value of the attitude angle of the quadrotor UAVs as

[\begin{matrix} ϕ & θ & ψ \end{matrix}] = [\begin{matrix} 30^{\circ} & - 15^{\circ} & 10^{\circ} \end{matrix}]

, the center vector of each RBF neural network is

c_{j} = [\begin{matrix} - 1 & - 0.5 & 0 \end{matrix}

\begin{matrix} 0.5 & 1 \end{matrix}]

, the width vector is

b_{j} = 3

To show the effectiveness and superiority of the proposed method clearly. Two kinds of disturbances are considered. The first is taken as

\begin{matrix} f_{1} (t_{k}) & = 0.15 + [\begin{matrix} 0.1 sin (0.2 π t_{k}) & 0 \end{matrix}] x_{1} (t_{k}) + 30 cos (x_{11} (t_{k}) x_{12} (t_{k})), \\ f_{2} (t_{k}) & = 0.15 + [\begin{matrix} 0.2 sin (0.4 π t_{k}) & 0 \end{matrix}] x_{2} (t_{k}) + 15 cos (x_{21} (t_{k}) x_{22} (t_{k})), \\ f_{3} (t_{k}) & = 0.15 + [\begin{matrix} 0.4 sin (0.3 π t_{k}) & 0 \end{matrix}] x_{3} (t_{k}) + 40 (1 - x_{31} (t_{k}) x_{32} (t_{k})) . \end{matrix}

The adaptive adjustment gain parameter

μ_{i}

in [25] is taken as 0.5 during simulation. The parameters of the proposed adaptive sliding mode controller are shown in Table 2.

The simulation results of the attitude angle response results of the quadrotor UAVs under the action of the traditional adaptive sliding-mode control algorithm (AC+SMC) and the proposed adaptive sliding-mode control algorithm (RBF+SMC) proposed in this paper are shown in Figure 3, Figure 4 and Figure 5. It can be seen from Figure 3, Figure 4 and Figure 5 that the traditional adaptive sliding-mode control algorithm (AC+SMC) cannot ensure that the attitude angle of the UAVs converges to the origin, the curve of each attitude angle has a distinct mutation and oscillation, and the overshoot is large, which cannot realize the stable control of the attitude of the quadrotor UAV.

The approximation results of the external disturbances

f_{1} (t_{k})

f_{2} (t_{k})

and

f_{3} (t_{k})

by the RBF neural network are shown in Figure 6, Figure 7 and Figure 8. It can be seen that the RBF neural-network-based adaptive scheme can estimate the unknown

f_{i} (t_{k})

quickly and well. In addition, Table 3, shows the performance comparisons for steady-state time and Table 4 shows the approximation time of RBF neural network to

f_{i} (t)

. Clearly, the quadrotor UAVs under proposed control offer a faster convergence speed.

To further illustrate the advantages of the proposed adaptive sliding mode control method, another kind of disturbances is taken as follows.

\begin{matrix} f_{1} (t_{k}) & = 0.1 + x_{22} (t_{k}) x_{32} (t_{k}) + 0.02 x_{22} (t_{k}), \\ f_{2} (t_{k}) & = 0.3 + x_{12} (t_{k}) x_{32} (t_{k}) + 0.03 x_{12} (t_{k}), \\ f_{3} (t_{k}) & = 0.02 + x_{12} (t_{k}) x_{22} (t_{k}) . \end{matrix}

One can see that the cross-coupling is considered here. Now, the design parameters of the traditional adaptive sliding mode control in [25] and the proposed adaptive sliding mode controller are shown in Table 5 and Table 6, respectively.

It can be seen from Figure 9, Figure 10 and Figure 11 that the traditional adaptive sliding-mode control algorithm (AC+SMC) cannot well ensure that the attitude angle of the UAVs converges to the origin. When there is a couplings problem, the curve of each attitude angle has a distinct mutation and oscillation, and the overshoot is large, i.e., it cannot realize the stable control of the attitude of the quadrotor UAV. At the same time, the flight attitude control can still estimate the unknown disturbance well (see Figure 12, Figure 13 and Figure 14) by the proposed design method in this paper (see Figure 9, Figure 10 and Figure 11). Table 7 also shows that the steady-state time of the proposed method is quicker for each attitude angle, and Table 8 shows that the adaptive RBF neural network can estimate

f_{i} (t)

in a quick way.

It can be seen from the above statements that the proposed adaptive sliding-mode control algorithm in the delta operator framework can effectively compensate the unknown nonlinear disturbances and achieve better control performance simultaneously, which is more suitable for use in the flight attitude control of quadrotor UAVs. However, this is based on the assumption that the actuator does not encounter various faults. How to design a sliding mode fault-tolerant attitude controller is an interesting and challenging problem to be solved and will be considered in our future works.

5. Conclusions

In this paper, we proposed an adaptive sliding-mode control algorithm in the delta domain for the attitude control of a quadrotor UAVs subject to external disturbances and couplings. First, the delta operator technique was used to discretize the attitude system of a quadrotor UAV. Then, a linear sliding surface, which can ensure the asymptotic stability of the sliding dynamics, was introduced in terms of the linear matrix inequality technique. Second, an adaptive sliding mode controller, using an RBF neural network to estimate external disturbances and couplings, was designed for the attitude reaching control.

We demonstrated, via Lyapunov stability theory, that the controller guaranteed that all attitudes of the quadrotor UAVs could be driven to the designed sliding surface, and thus attitude control was achieved. Finally, simulation result comparisons verified the effectiveness and superiority of the proposed adaptive sliding-mode attitude control algorithm proposed in this paper. In our future research, how to apply the proposed theoretical method to an actual unmanned system will be our focus and will be studied in depth.

Author Contributions

Conceptualization, B.Z.; formal analysis, Y.W., investigation, H.L.; software, Z.C.; Writing—review and editing, B.Z.; writing—original draft preparation, Y.W. and H.L.; Z.C. provided support in the whole writing of this paper. All authors have read and agreed to the published of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China under grant number 61973169, the Jiangsu Natural Science Foundation under grant number BK20201392, the Qing Lan Project of Jiangsu Province under grant number R2021Q04.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

Noordin, A.; Basri, M.A.M.; Mohamed, Z.; Lazim, I.M. Adaptive PID controller using sliding mode control approaches for quadrotor UAVs attitude and position stabilization. Arab. J. Sci. Eng. 2021, 46, 963–981. [Google Scholar] [CrossRef]
Wang, S.Y.; Polyakov, A.; Zheng, G. Quadrotor stabilization under time and space constraints using implicit PID conttroler. J. Frankl. Institude 2022. [Google Scholar] [CrossRef]
Sun, C.; Liu, M.; Liu, C.; Feng, X.; Wu, H. An industrial quadrotor UAVs control method based on fuzzy adaptive linear active disturbance rejection control. Aerosp. Sci. Technol. 2020, 10, 376. [Google Scholar] [CrossRef]
Wang, Z.J.; Zhao, T. Adaptive-based linear active disturbance rejection attitude control for quadrotor with external disturbances. Trans. Inst. Meas. Control 2022, 44, 286–298. [Google Scholar] [CrossRef]
Zhang, Y.; Nie, Y.; Chen, L. Adaptive fuzzy fault-tolerant control against time-varying faults via a new sliding mode observer method. Symmetry 2021, 13, 1615. [Google Scholar] [CrossRef]
Hong, Q.; Shi, Y.; Chen, Z. Adaptive sliding mode control based on disturbance observer for placement pressure control system. Symmetry 2020, 12, 1057. [Google Scholar] [CrossRef]
Xiong, J.J.; Zhang, G.B. Discreti-time sliding mode control for a quadrotor UAV. Optik 2016, 127, 3718–3722. [Google Scholar] [CrossRef]
Wu, X.W.; Xiao, B.; Qu, Y.H. Modeling and sliding mode-based attitude tracking control of a quadrotor UAVs with time-varying mass. ISA Trans. 2019. [Google Scholar] [CrossRef]
Hou, Z.; Lu, P.; Tu, Z. Nonsingular terminal sliding mode control for a quadrotor UAVs with a total rotor failure. Electronics 2021, 98, 105716. [Google Scholar] [CrossRef]
Yang, S.Y.; Han, J.; Xia, L.; Chen, Y.-H. Adaptive sliding mode attitude controller such that all attitudes of the quadrotor UAVs can be driven to the designed sliding surface. ISA Trans. 2021, 106, 12–30. [Google Scholar] [CrossRef]
Zhang, J.; Ren, Z.; Deng, C.; Wen, B. Adaptive fuzzy global sliding mode control for trajectory tracking of quadrotor UAVs. Nonlinear Dyn. 2019, 97, 609–627. [Google Scholar] [CrossRef]
Nekoukar, V.; Dehkordi, M.M. Robust path tracking of a quadrotor using adaptive fuzzy terminal sliding mode control. Control. Eng. Pract. 2021, 110, 104763. [Google Scholar] [CrossRef]
Eliker, K.; Zhang, W. Finite-time adaptive integral backstepping fast terminal sliding mode control applications on quadrotor UAV. Int. J. Control. Autom. Syst. 2020, 18, 415–430. [Google Scholar] [CrossRef]
Zhang, X.; Wang, Y.; Zhu, G.; Chen, X.; Su, C.-Y. Discrete-time adaptive neural tracking control and its Enperiments for quadrotor unmanned aerial vehicle systems. IEEE/ASME Trans. Mechatronics 2022. [Google Scholar] [CrossRef]
Middleton, R.; Goodwin, G.C. Improved finite word length characteristics in digital control using delta operators. IEEE Trans. Automat. Contr. 1986, 31, 1015–1021. [Google Scholar] [CrossRef]
Middleton, R.H.; Goodwin, G.C. Digital Control and Estimation: A Unified Approach; Prentice Hall: Hoboken, NJ, USA, 1990. [Google Scholar]
Zheng, B.C.; Yu, X.; Xue, Y. Quantized sliding mode control in delta operator framework. Int. J. Robust Nonlinear Control 2018, 28, 519–535. [Google Scholar] [CrossRef]
Yang, H.; Xia, Y.; Shi, P.; Zhao, L. Analysis and Synthesis of Delta Operator Systems; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
Sakthivel, R.; Ramya, L.S.; Ma, Y.K.; Malik, M.; Leelamani, A. Stabilization of uncertain switched discrete-time systems against actuator faults and input saturation. Nonlinear-Anal.-Hybrid Syst. 2020, 35, 519–535. [Google Scholar] [CrossRef]
Guo, X.; Yang, G.-H. Insensitive output feedback H-infinity control of delta operator systems with insensitivity to sampling time jitter. Int. J. Robust Nonlinear Control 2014, 24, 725–743. [Google Scholar] [CrossRef]
Guo, X.; Yang, G.-H. Low-sensitivity H∞ filter design for linear delta operator systems with sampling time jitter. Int. J. Control 2012, 85, 397–408. [Google Scholar] [CrossRef]
Pu, X.; Ren, L.; Liu, Y.; Ji, L.H. Group consensus of multi-agent systems with cooperative-competitive interaction and communication delay in switching topologies networks based on the delta operator method. Neurocomputing 2020, 390, 57–68. [Google Scholar] [CrossRef]
Pu, X.; Ren, L.; Liu, Y.; Pu, R. Couple-group consensus for heterogeneous MASs under switched topologies in cooperative-competitive systems: A hybrid pinning and delta operator skills. Neurocomputing 2021, 441, 335–349. [Google Scholar] [CrossRef]
Chen, G.S. Study on Anti-Jamming and Fault-Tolerant Control of Quadrotor UAV; Nanjing University of Information Science and Technology: Nanjing, China, 2020. [Google Scholar] [CrossRef]
Zhang, Z.C. Sliding Mode Fault-Tolerant Control Algorithm for Delta Operator Time-Delay Systems; Nanjing University of Aeronautics and Astronautics: Nanjing, China, 2016. [Google Scholar]
Zhao, B.; Xian, B.; Zhang, Y.; Zhang, X. Nonlinear robust adaptive tracking control of a quadrotor UAVs via immersion and invariance methodology. IEEE Trans. Ind. Electron. 2015, 62, 2891–2902. [Google Scholar] [CrossRef]
Lendoul, F.; Lara, D.; Fantoni-Coichot, I.; Lozano, R. Real-time nonlinear embedded control for an autonomous quadrotor helicopter. J. Guid. Control. Dyn. 2007, 30, 1049–1060. [Google Scholar] [CrossRef]
Zhao, D.; Sun, C.; Wang, Q.; Yang, W. Neural-network-based PID control for quadrotor aircraft. In Proceedings of the International Conference on Intelligent Science and Big Data Engineering; Springer: Cham, Switzerland, 2015; pp. 287–297. [Google Scholar]
Wang, W.; Ma, H.; Xia, M.; Weng, L.; Ye, X. Attitude adn altitude controller design for quad-rotor type MAVs. Math. Probl. Eng. 2013, 2013, 587098. [Google Scholar]
Wang, S.B.; Han, Y.; Chen, J.; Zhang, Z.C.; LIu, X.Z. Attitude active disturbance rejection control of UAVs attitude based on iterative learning control. Acta Aeronaut. Astronaut. Sin. 2020, 41, 324112. (in Chinese) [Google Scholar] [CrossRef]
Du, H.; Zhu, W.; Wen, G.; Wu, D. Finite time formation control for a group of quadrotor aircraft. Aerosp. Sci. Technol. 2017, 69, 609–616. [Google Scholar] [CrossRef]
Li, H.; Wu, B.; Li, G.; Yang, C. Global Theory of Delta Operator Control and Its Robustness Control; National Defense Inddustry Press: Beijing, China, 2005. [Google Scholar]
Boyd, S.; El Ghaoui, L.; Feron, E.; Balakrishnan, V. Linear matrix inequalities in system and control theory. In Society for Industrial and Applied Mathematics; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1994. [Google Scholar]
Liu, J. Intelligent Control Design and MATLAB Simulation; Tsinghua University Press: Beijing, China; Springer Nature Singapore Pte Ltd.: Singapore, 2018. [Google Scholar]

Figure 1. RBF neural network structure diagram.

Figure 2. Diagram of the RBF neural-network-based adaptive sliding-mode control algorithm.

Figure 3. Response results of the attitude angle

ϕ

Figure 3. Response results of the attitude angle

ϕ

Figure 4. Response results of the attitude angle

θ

Figure 4. Response results of the attitude angle

θ

Figure 5. Response results of the attitude angle

ψ

Figure 5. Response results of the attitude angle

ψ

Figure 6. Neural network approximation results of the external disturbance

f_{1} (t_{k})

Figure 6. Neural network approximation results of the external disturbance

f_{1} (t_{k})

Figure 7. Neural network approximation results of the external disturbance

f_{2} (t_{k})

Figure 7. Neural network approximation results of the external disturbance

f_{2} (t_{k})

Figure 8. Neural network approximation results of the external disturbance

f_{3} (t_{k})

Figure 8. Neural network approximation results of the external disturbance

f_{3} (t_{k})

Figure 9. Response results of the attitude angle

ϕ

Figure 9. Response results of the attitude angle

ϕ

Figure 10. Response results of the attitude angle

θ

Figure 10. Response results of the attitude angle

θ

Figure 11. Response results of the attitude angle

ψ

Figure 11. Response results of the attitude angle

ψ

Figure 12. Neural network approximation results of the external disturbance

f_{1} (t_{k})

Figure 12. Neural network approximation results of the external disturbance

f_{1} (t_{k})

Figure 13. Neural network approximation results of the external disturbance

f_{2} (t_{k})

Figure 13. Neural network approximation results of the external disturbance

f_{2} (t_{k})

Figure 14. Neural network approximation results of the external disturbance

f_{3} (t_{k})

Figure 14. Neural network approximation results of the external disturbance

f_{3} (t_{k})

Table 1. Physical parameters of the quadrotor UAV.

Parameters	Values	Units
m	2.5	kg
l	0.245	m
g	9.8	m/s $^{2}$
c	0.575	m
$K_{1}$	0.148	N s/m
$K_{2}$	0.148	N s/m
$K_{3}$	0.164	N s/m
$J_{1}$	0.05	kg·m $^{2}$
$J_{2}$	0.05	kg·m $^{2}$
$J_{3}$	0.12	kg·m $^{2}$

Table 2. Controller parameters of the quadrotor UAV.

Attitude Angle Systems	$ε_{i}$	$β_{i}$	$γ_{i}$
Roll angle system	100	20	30
Pitch angle system	200	20	30
Yaw angle system	120	20	40

Table 3. Steady-state time of attitude angle under different control algorithms.

Attitude Angle Systems	Roll Angle System	Pitch Angle System	Yaw Angle System
RBF + SMC	16.8 s	15.2 s	13.7 s
AC + SMC	/	/	/

Table 4. Approximation time of RBF neural network to

f_{i} (t)

Table 4. Approximation time of RBF neural network to

f_{i} (t)

Attitude Angle Systems	Roll Angle System	Pitch Angle System	Yaw Angle System
Effective approximation time	4.6 s	7.5 s	1.45 s

Table 5. Controller parameters of the quadrotor UAVs proposed in [25].

Attitude Angle Systems	$μ_{i}$
Roll angle system	4
Pitch angle system	3
Yaw angle system	4

Table 6. Controller parameters of the quadrotor UAVs proposed in this paper.

Attitude Angle Systems	$ε_{i}$	$β_{i}$	$γ_{i}$
Roll angle system	120	20	10
Pitch angle system	200	20	5
Yaw angle system	300	20	5

Table 7. Steady-state time of attitude angle under different control algorithms.

Attitude Angle Systems	Roll Angle System	Pitch Angle System	Yaw Angle System
RBF + SMC	15.7 s	13.7 s	16.4 s
AC + SMC	21.6 s	/	25.8 s

Table 8. Approximation time of RBF neural network to

f_{i} (t)

Table 8. Approximation time of RBF neural network to

f_{i} (t)

Attitude Angle Systems	Roll Angle System	Pitch Angle System	Yaw Angle System
Effective approximation time	1.63 s	1.79 s	2.06 s

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Zheng, B.; Wu, Y.; Li, H.; Chen, Z. Adaptive Sliding Mode Attitude Control of Quadrotor UAVs Based on the Delta Operator Framework. Symmetry 2022, 14, 498. https://doi.org/10.3390/sym14030498

AMA Style

Zheng B, Wu Y, Li H, Chen Z. Adaptive Sliding Mode Attitude Control of Quadrotor UAVs Based on the Delta Operator Framework. Symmetry. 2022; 14(3):498. https://doi.org/10.3390/sym14030498

Chicago/Turabian Style

Zheng, Bochao, Yuewen Wu, Hui Li, and Zhipeng Chen. 2022. "Adaptive Sliding Mode Attitude Control of Quadrotor UAVs Based on the Delta Operator Framework" Symmetry 14, no. 3: 498. https://doi.org/10.3390/sym14030498

APA Style

Zheng, B., Wu, Y., Li, H., & Chen, Z. (2022). Adaptive Sliding Mode Attitude Control of Quadrotor UAVs Based on the Delta Operator Framework. Symmetry, 14(3), 498. https://doi.org/10.3390/sym14030498

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

Article Menu

Adaptive Sliding Mode Attitude Control of Quadrotor UAVs Based on the Delta Operator Framework

Abstract

1. Introduction

2. System Model Description of Quadrotor UAV

3. Design of the Sliding Mode Attitude Control Algorithm

3.1. Sliding Surface Design

3.2. Adaptive Sliding Mode Controller Design

4. Illustrative Example

4.1. Parameter Settings and Modelling of the Quadrotor UAV

4.2. Design of Linear Surfaces

4.3. Simulation Retsults

5. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI