Next Article in Journal
Neural-Network-Assisted Polar Code Decoding Schemes
Previous Article in Journal
Analysis of Telescope Wavefront Aberration and Optical Path Stability in Space Gravitational Wave Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Based on Backpropagation Neural Network and Adaptive Linear Active Disturbance Rejection Control for Attitude of a Quadrotor Carrying a Load

College of Automation and Electronic Engineering, Qingdao University of Science and Technology, Qingdao 266061, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(24), 12698; https://doi.org/10.3390/app122412698
Submission received: 19 November 2022 / Revised: 8 December 2022 / Accepted: 9 December 2022 / Published: 11 December 2022

Abstract

:
A control method combining Backpropagation (BP) neural network and Adaptive Linear Active Disturbance Rejection Control (ALADRC) is proposed for the attitude control problem of quadrotor aircraft. The proposed controller can observe and compensate the total disturbance in the working process of the quadrotor system. At the same time, it has a good ability to suppress the disturbance of the quadrotor load mass change. In addition, adaptive control and BP neural network are used to adjust the controller parameters in real time to solve the problem of difficult parameter tuning. Finally, the stability of the system is proved based on Lyapunov theory. The simulation results show that the quadrotor system can still track the altitude and attitude angle commands stably in the presence of disturbances, and has strong adaptability to load mass changes, so that the quadrotor can still complete the given tasks in the presence of multi-source disturbances.

1. Introduction

Today, quadrotors have become popular in many application scenarios due to their high flexibility and high degrees of freedom, such as power inspection, data collection, surveillance and reconnaissance missions, forest fire rescue, and precision agriculture [1,2,3].
The quadrotor system is a typical underactuated system, which has the characteristics of nonlinearity, strong coupling, and susceptibility to interference, so it is difficult to establish an accurate mathematical model. A way to design a controller with high stability and strong robustness is a difficult research point. Due to the advantages of simple structure and easy implementation, linear controllers are widely used in quadrotor control systems. Research by Shakeel et al. [4] presented that the Linear Quadratic Regulator (LQR) controller has a good control effect by analyzing and evaluating several different control schemes of the quadrotor. However, many scholars still try to use modern control methods to improve the control performance of quadrotor controllers, such as predictive control [5,6], sliding film control [7], LQR [8,9] and so on. Based on backstepping control and Lyapunov theory, research by Rodríguez et al. [10] presented a robust control algorithm in order to increase the endurance time of the quadrotor and reduce unnecessary energy dissipation. Research by Nigro et al. [11] presented a new type of quadrotor aircraft by adding a driving gimbal mechanism to provide multi-directional thrust. At the same time, a double closed-loop motion control scheme was proposed to enhance the robustness of the system. However, the above studies are all based on the nominal conditions of the quadrotor, that is, the influence of disturbance is ignored. This leads to the poor control effect of the controller designed in some scenes with large disturbances.
The quadrotors are usually affected by external unknown wind disturbances and unmodeled dynamics within the system, which will lead to poor stability of the quadrotor system. Designing a robust controller to overcome the influence of disturbances to ensure the stability of the quadrotor system has always been a key topic for some scholars. For example, considering the uncertainty of system parameters, research by Fan et al. [12] presented a control scheme combining neural network and adaptive control, and provided the simulation results. Based on the backstepping control theory, research by Chen et al. [13] presented a nonlinear adaptive integral backstepping control algorithm. The simulation results prove that the controller still has strong robustness in the presence of wind disturbance. However, these studies did not address specific mathematical models of wind disturbance, nor did they explain how wind disturbance acts on the system.
Nowadays, in order to solve the problem of unknown disturbance of the quadrotor system, many scholars have carried out in-depth research on the quadrotor control system based on the Active Disturbance Rejection Control (ADRC) theory. The ADRC does not require an accurate mathematical model of the controlled object, and its structure is simple, so it can be used in the controller design of a quadrotor system [14]. Aiming at the problem that the ADRC controller has many parameters, it is difficult to tune these parameters. Research by Zhang et al. [15] applied the fractional fuzzy particle swarm optimization (FOFPSO) algorithm to optimize the ADRC controller parameters, and it was compared with several other optimization algorithms to optimize the ADRC parameters. Finally, the authors propose a novel wind field model; this model is employed to prove the effectiveness of the controller under the action of wind disturbance. Research by Yuan and Zhao et al. [16,17], aiming at the problem of grid voltage fluctuation in power grid systems such as wind power grid-connected inverters, proposed an improved LADRC controller to improve the disturbance observation accuracy of the system by reducing the phase lag of the observer. The results showed that the improved method has better rapidity and immunity.
Linear Active Disturbance Rejection (LADRC) is widely used because of the small number of controller parameters and convenient adjustment [18,19,20,21]. Although the number of controller parameters that need to be adjusted by LADRC is much smaller than that of ADRC, it is still a difficult task to readjust the controller parameters for different systems or the same system in different working scenarios. At the same time, Linear Extension State Observer (LESO) in LADRC often has the problem of limited bandwidth, which may cause poor control effect in practical applications [22].
For the problem of quadrotor attitude control, this paper designs a quadrotor attitude control strategy based on LADRC. The control strategy not only considers the internal and external disturbances of the system, but also has a strong adaptability to load quality changes. In order to solve the problem of LADRC controller having too many parameters and being difficult to adjust manually, BP neural network and adaptive control algorithm are proposed in this paper. The BP neural network is used to dynamically optimize and adjust the bandwidth value of LESO in real time, so that the controller can still maintain good control effect in the complex and changeable working environment. By adjusting the values of k p and k d in the PD controller in real time, adaptive control reduces the parameter error caused by manual parameter adjustment, and compensates the observation error caused by bandwidth selection, which improves the control accuracy of LADRC and the dynamic performance of the system. Finally, in order to verify the effectiveness of the control scheme, a simulation experiment is designed to verify the attitude tracking ability of the controller.

2. Dynamic Model of Quadrotor Aircrafts

The coordinate system and structure of the quadrotor are shown in Figure 1. We use six variables to describe the motion of the quadrotor, namely the quadrotor position x , y , z , the attitude angle ϕ , θ and ψ [23]. The centroid of quadrotor moves upward continuously from below the origin of the fuselage mount when the load is gradually reduced. Referring to the body coordinate system B and the ground fixed coordinate system E, the dynamic model of the variable mass quadrotor aircraft is established [24].
As shown in Figure 1, the rotation matrix from coordinate system B to E is R B E [25,26].
The positional motion of the quadrotor satisfies Newton’s second law,
{ X ˙ E = V E ( M Δ m ) V ˙ E = L E + G M E + G Δ m E + f 1 E ,
where X E = [ x y z ] T ; V E is the linear velocity; M is the mass of the quadrotor; L E is the total lift; G Δ m E = [ 0 0 Δ m g ] T ; G M E = [ 0 0 M g ] T ; f 1 E = [ k x ˙ k y ˙ k z ˙ ] T is the matrix of air resistance.
The lift force of the quadrotor is expressed as
L B = [ L x L y L z ] T = [ 0 0 L ] T ,
where L = i = 1 4 L i = L 1 + L 2 + L 3 + L 4 ; from R B E , we can obtain L E :
L E = R B E L B = L [ sin ψ sin ϕ + cos ϕ cos ψ sin θ sin ϕ cos ψ + cos ϕ sin ψ sin θ cos ϕ cos θ ] .
In summary, we can obtain the positional dynamics equation of the quadrotor.
{ x ¨ = k f x ˙ ( M     Δ m ) + L ( sin ψ sin ϕ   +   cos ψ cos ϕ sin θ ) ( M     Δ m ) y ¨ = k f y ˙ ( M     Δ m ) + L ( cos ψ sin ϕ   +   sin ψ cos ϕ sin θ ) ( M     Δ m ) z ¨ = k f z ˙ ( M     Δ m ) + L cos ϕ cos θ ( M     Δ m ) g .
As the load mass decreases, the centroid of quadrotor gradually moves upward, thereby affecting the attitude angle of the quadrotor. Because the quadrotor is loaded by the mounting platform, the two can be considered as a whole. The initial coordinates of the centroid of quadrotor are ( x 1 , y 1 , z 1 ) , the coordinates of the centroid of the part of the load that is reduced during operation are ( x 2 , y 2 , z 2 ) . When the load mass decreases, the change in the coordinates of the centroid of quadrotor is shown in the following formula:
{ x a = ( x 1 M x 2 Δ m ) / ( M Δ m ) y a = ( y 1 M y 2 Δ m ) / ( M Δ m ) z a = ( z 1 M z 2 Δ m ) / ( M Δ m ) ,
where the coordinates of the centroid of the quadrotor are ( x a , y a , z a ) . According to Figure 1, initial coordinate of the centroid of the quadrotor is ( 0 , 0 , 0 ) , the coordinate of the centroid of the reduced load mass during the working process is ( 0 , 0 , z c ) . Substituting it into Formula (6), the coordinate of the centroid of the quadrotor can be obtained as ( 0 , 0 , z a ) .
z a = z c m / ( M Δ m ) .
According to the Newton–Euler equation, the attitude dynamics equation can be expressed as [27].
The attitude angle motion of the quadrotor satisfies the Newton–Euler equation, and its attitude dynamic model can be expressed as [28,29]
{ I 3 × 3 τ ˙ B + τ B × I 3 × 3 τ B = ς B ς B = Q B + f 2 B + P B ,
where I 3 × 3 = d i a g ( I x , I y , I z ) is the inertia matrix; τ B = [ ϕ ˙ θ ˙ ψ ˙ ] T is the angular velocity; ς B is the resultant moment; P B and Q B are the pulling torque and reaction torque of the four-wheel drive motor, respectively; f 2 B = [ l k f ϕ ˙ l k f θ ˙ k f ψ ˙ ] T is the matrix of air resistance.
P B = [ b 1 l 2 z a 2 ( L 1 + L 2 + L 3 L 4 ) b 1 l 2 z a 2 ( L 1 L 2 + L 3 + L 4 ) 0 ] ,
Q B = [ 0 0 b 2 ( L 1 L 2 + L 3 L 4 ) ] .
In summary, we can obtain attitude dynamics equation of the quadrotor.
{ ϕ ¨ = l k f ϕ ˙ / I x + b 1 l 2 z a 2 ( L 1 + L 2 + L 3 L 4 ) / I x + ( I y I z ) ψ ˙ θ ˙ / I x ψ ¨ = k f ψ ˙ / I z + b 2 ( L 1 L 2 + L 3 L 4 ) / I z + ( I x I y ) ϕ ˙ θ ˙ / I z θ ¨ = l k f θ ˙ / I y + b 1 l 2 z a 2 ( L 1 L 2 + L 3 + L 4 ) / I y + ( I z I x ) ϕ ˙ ψ ˙ / I y ,
where ψ , ϕ , and θ are the yaw angle, roll angle and pitch angle, respectively.
To simplify the quadrotor dynamic equation, ( U 1 , U 2 , U 3 , U 4 ) is used to represent the relevant parameters.
[ U 1 U 2 U 3 U 4 ] = [ 1 ( M Δ M ) 1 ( M Δ M ) 1 ( M Δ M ) 1 ( M Δ M ) b 1 I x b 1 I x b 1 I x b 1 I x b 1 I y b 1 I y b 1 I y b 1 I y b 2 I z b 2 I z b 2 I z b 2 I z ] [ L 1 L 2 L 3 L 4 ] .
In summary, we can obtain the dynamic equation of the quadrotor.
{ x ¨ = k f x ˙ / ( M Δ m ) + U 1 ( sin ϕ sin ψ + cos ψ cos ϕ sin θ ) y ¨ = k f y ˙ / ( M Δ m ) + U 1 ( sin ϕ cos ψ + sin ψ cos ϕ sin θ ) z ¨ = k f z ˙ / ( M Δ m ) g + U 1 cos ϕ cos θ θ ¨ = U 3 l 2 z a 2 + ( I z I x ) ϕ ˙ ψ ˙ / I y l k f θ ˙ / I y ϕ ¨ = U 2 l 2 z a 2 + ( I y I z ) θ ˙ ψ ˙ / I x l k f ϕ ˙ / I x ψ ¨ = U 4 + ( I x I y ) ϕ ˙ θ ˙ / I z k f ψ ˙ / I z .
Equation (13) can also be written as [30]
{ x ˙ 1 = x 2 x ˙ 2 = A u ( n ) + B + c ,
where x 1 = [ z ϕ θ ψ ] T , u ( n ) = [ U 1 U 2 U 3 U 4 ] T ;
A = [ cos θ cos ϕ 0 0 0 0 l 2 z a 2 0 0 0 0 l 2 z a 2 0 0 0 0 1 ] B = [ g k f z ˙ / ( M Δ M ) l k f ϕ ˙ / I x + ( I y I z ) ψ ˙ θ ˙ / I x l k f θ ˙ / I y + ( I z I x ) ψ ˙ ϕ ˙ / I y k f ψ ˙ / I z + ( I x I y ) ϕ ˙ θ ˙ / I z ] .

3. Control Scheme Design

Aiming at the problems of strong coupling, nonlinearity and mass change of the quadrotor system, we proposed a control scheme combining BP neural network, adaptive controller and LADRC. The BP-ALDRC control system of quadrotor under the condition of continuous time-varying wind field disturbance and changing load quality is established. Figure 2 is a block diagram of the control scheme.
In this Figure, x d is the given signal, [ z 1 z 2 z 3 ] is the observation vector of LESO to the system state variable, [ β 1 β 2 β 3 ] is the bandwidth gain vector of LESO, [ γ 1 γ 2 ] is the approximate value of the given signal and its differential, u p is a nonlinear state error feedback control law with PD controller output, u is the control law of the system, σ is the error combination, k ^ p , k ^ d are the estimated values of PD controller parameters k p , k d , respectively, M ^ c is the estimated value of load mass change M c , x is the state variable of the system, b 0 is the compensation factor.
The advantage of using LADRC to control the quadrotor is that there is no need to decouple the channels. LADRC directly estimates system internal disturbances such as system coupling and environmental external disturbances in the working process as the total disturbance for observation and compensation. Aiming at the difficulty of setting LADRC parameters and avoiding errors caused by manual parameter adjustment, we use adaptive control to adjust the PD controller parameters in real time.

4. Control Methodology

4.1. Design of LADRC Controller

The dynamic Equation (12) can be transformed into
{ ϑ ˙ 1 = ϑ 2 ϑ ˙ 2 = q ( t ) + d ( t ) + G u ( x ) y = ϑ 1 ,
where ϑ = [ ϑ 1 ϑ 2 ] T are measurable state variables; y is the output signal; f ( t ) is an unknown nonlinear function of the system; d ( t ) is the external disturbance. G a is partially known, and the known part is denoted as G 0 .
All uncertainties in the system and all disturbances received by the system are regarded as the total disturbance, which is defined as
F ( ϑ , u ( ϑ ) ) = f ( t ) + [ g ( t ) b 0 ] u ( ϑ ) + d ( t ) .
To observe the total disturbance, the extended state space is defined ζ = [ ζ 1 ζ 2 ζ 3 ] , where ζ 1 = ϑ 1 , ζ 2 = ϑ ˙ 1 and ζ 3 = F ( ϑ , u ( ϑ ) ) .
The nonlinear system (12) can be transformed into
{ ζ ˙ = M ζ + N u + T h y = U ζ ,
where M = [ 0 1 0 0 0 1 0 0 0 ] , N = [ 0 b 0 0 ] , T = [ 0 0 1 ] , D = [ 1 0 0 ] .
According to (16), LESO is designed as follows:
{ z ˙ = M z ^ + N u + γ ( y y ^ ) y ^ = U z ,
where z = [ z 1 z 2 z 3 ] is the state vector of the observer, which is used to observe the state of system ζ = [ ζ 1 ζ 2 ζ 3 ] ; γ = [ γ 1 γ 2 γ 3 ] T is the observation gain vector of LESO.
The gain value of LESO is obtained from its characteristic equation:
λ ( s ) = s 3 + ω 1 s 2 + ω 2 s + ω 3 = ( s + ω 0 ) 3 ,
where ω 0 is the bandwidth of LESO, denoted as ω = [ 3 ω 0 3 ω 0 2 ω 0 3 ] T .
Remark 1.
The high bandwidth gain of LESO helps to improve the control accuracy of LADRC, but it will cause the output of the controller to increase, and the actuator will consume more energy, which is not conducive to the long-term flight of the quadrotor. This study introduces BP neural network to improve the overall performance of the controller to solve this problem.
The extended state space observation value of LESO is defined as z ^ 1 = x ^ 1 , z ^ 2 = x ^ 2 , z ^ 3 = F ^ ( x , u ( n ) ) . The control law can be written in the following form:
u = ( u p z ^ 3 ) / b 0 ,
where u p = k p ( x d z ^ 1 ) + k d ( x ˙ d z ^ 2 ) + x ¨ d is the output of the PD controller.

4.2. Parameter Tuning Adaptive Law Design

By applying APC to dynamically adjust the k p and k d values in the PD controller, the observation error caused by bandwidth selection is compensated. It can not only simplify the controller parameter tuning process, but also improve the control precision of LADRC and the dynamic performance of the system.
Remark 2.
In order to simplify the adjustment of parameter k p , k d and avoid errors caused by manual parameter adjustment, we introduce an adaptive law which can adjust k p , k d timely. Formulas (36) and (37) provide the adaptive law of k p , k d .
Assumption 1: Assume that a given signal x d ( t ) is smooth and bounded such that x d ( t ) X d is true. Here X d is a known constant.
Define tracking error:
e 1 = x d z 1 .
Derivation can be obtained as follows:
e ˙ 1 = x ˙ d z ˙ 1 = x ˙ d z 2 .
Define the filter tracking error:
σ = [ ξ T   1 ] [ e 1 e ˙ 1 ] ,
where ξ T = a 1 is a suitable gain vector such that when σ 0 , e 1 0 is satisfied.
Derivation of σ can be obtained:
σ ˙ = x ¨ d z ˙ 2 + [ 0 ξ T ] [ e 1 e ˙ 1 ] .
From the output of LESO, we can obtain
x ˙ 2 = q ( t ) + N u ( x ) + d ( t ) = F ( x , u ( x ) ) .
The input of the control is as follows:
δ = μ σ + u p [ 0 ξ T ] [ e 1 e ˙ 1 ] ,
where μ is any positive number.
In summary, we can obtain
σ ˙ = μ σ + k p ( x d z ^ 1 ) + k d ( x ˙ d z ^ 2 ) .

4.3. BP Neural Network Design

The continuous time-varying wind field disturbance and the continuous reduction of the load mass will cause the control effect of the controller to deteriorate. LESO observes the system status and the total disturbance. However, in the face of continuous time-varying wind field disturbance and changing load quality, the fixed observer bandwidth will lead to larger observation errors and worse control effects. To solve this problem, this research employs BP neural network (BP) to adjust the bandwidth gain value of LESO online, enhance the observation performance of LESO, and compensate for adverse effects of wind disturbances and changes in load mass.
Remark 3.
LESO can observe each state variable and total disturbance of the quadrotor system according to the input and output of the system. However, due to the discretization of the observer, the bandwidth value of the LESO cannot exceed the sampling period of the system, otherwise it will lead to the divergence of the LESO observations. If the bandwidth value is small, the LESO observation accuracy will be low. Therefore, we introduce APC to compensate the observation error caused by bandwidth selection by dynamically adjusting the k p and k d values in the PD controller, so as to improve the control accuracy of LADRC and the dynamic performance of the system. This approach will make LADRC more adaptable and extensive.
The structure of BP is shown in Figure 3.
The output error is defined as
e 2 = x d x 1 ,
where x d is a given signal, x 1 is the state variable of the system; derivation of e 2 can be obtained in the following way:
e ˙ 2 = x ˙ d x ˙ 1 = x ˙ d x 2 .
The input of BP neural network can be expressed as
I N = [ e 2 e ˙ 2 ] T .
The input and output of the hidden layer are
{ H a i n = a = 1 6 j a i n I N H a o u t = tan s i g ( H a i n ) a = 1 , 2 , 3 , 4 , 5 , 6 .
The hidden unit is 6, and the excitation function adopts the tangent sigmoid transfer function:
tan s i g ( k ) = 2 e 2 k + 1 1 .
Remark 4.
BP neural network adjusts the bandwidth gain vector  [ β 1 β 2 β 3 ] of LESO in real time according to the error information of the system. When the system produces output error, the BP neural network will output and adjust  [ β 1 β 2 β 3 ] according to the system error information to reduce the observation error of LESO, thereby reducing the system output error and improving the control accuracy of the system.
The BP neural network that adjusts the LESO bandwidth gain value is divided into three layers of network structure; the input signal of the BP is output error e 2 , the differential signal e ˙ 2 of the output error and the bias value k . The number of neurons in the hidden layer is 6 to ensure the strong ability of BP neural network to approximate the error. The output layer contains three neurons, corresponding to the gain vector of LESO. The output of BP neural network can be expressed by the following Formula (34):
y ( k ) = f ( i = 1 l ω i k g ( m = 1 M ω m i v m ) ) .
Since the BP neural network will normalize the input signal [ e 2 e ˙ 2 ] , its output signal will also be mapped to the range ( 1 1 ) , so the output of the BP cannot be directly used as the bandwidth gain value of LESO. To make the BP−ALADRC controller converge, the output is denormalized to meet the requirement of the gain vector [ β 1 β 2 β 3 ] of LESO in the BP−ALADRC controller. The gain vector of LESO can be obtained as
β k = β k + f ( i = 1 l ω i k f ( m = 1 M ω m i v m ) ) · k .

4.4. Mass Change Adaptive Law Design

The mass of the load carried by the quadrotor may change continuously during the actual working process, which is not conducive to the stability of the system. Therefore, we establish a mass adaptive law based on the height channel to reduce the adverse effects of mass changes on system stability.
According to Formula (12), we can obtain
χ ˙ z = U 1 cos θ cos ϕ g k f z ˙ / ( M Δ m ) z ¨ d + [ 0 ξ T ] [ e 1 e ˙ 1 ] .
According to the equivalent deterministic principle, we can obtain
U 1 = α χ z [ 0 ξ T ] [ e 1 e ˙ 1 ] + z ¨ d + g + k f z ˙ / ( M Δ m ^ ) cos ϕ cos θ ,
where α > 0 is an adjustable parameter. Bringing (37) into (36) can obtain
χ ˙ z = k f z ˙ M Δ M ^ α χ z k f z ˙ M Δ M .
Then, M ^ g = 1 M     Δ m ^ , M g = 1 M     Δ m are defined, then the estimation error can be expressed as
M ˜ g = M ^ g M g .
Substituting Formula (39) into Formula (38) can obtain
χ ˙ z = α χ z + k f z ˙ M ^ g k f z ˙ M g .

5. Stability Analysis

To ensure that the system tracking error converges to the zero neighborhood, we design the following adaptive law:
k ^ ˙ p = [ δ k p ( x d z ^ 1 ) φ ] / k ˜ p ,
k ^ ˙ d = [ δ k d ( x ˙ d z ^ 2 ) ζ ] / k ˜ d ,
M ^ ˙ g = κ χ ˙ z k f z ˙ .
Remark 5.
To avoid k ˜ p and k ˜ d being zero, k ^ ˙ p and k ^ ˙ d are integrated. The above formula can be rewritten as
k ^ p = 0 t [ δ k p ( x d z ^ 1 ) φ ] / k ˜ p d τ + k p 0 ,
k ^ d = 0 t [ δ k d ( x ˙ d z ^ 2 ) ζ ] / k ˜ d d τ + k d 0 .
Appropriate constants k p 0 and k d 0 are chosen to make k ˜ p = k p k ^ p and k ˜ d = k d k ^ d non-zero.
Proof of Theorem 1.
A positive definite Lyapunov function is designed as follows.
V = 1 2 σ 2 + 1 2 χ z 2 + 1 2 k ˜ p φ k ˜ p + 1 2 k ˜ d ζ k ˜ d + 1 2 M ˜ g κ M ˜ g .
Taking the derivative of Equation (46),
V ˙ = σ σ ˙ + χ z χ ˙ z + k ˜ p φ 1 k ˜ ˙ p + k ˜ d ζ 1 k ˜ ˙ d + M ˜ g κ 1 M ˜ ˙ g .
Bringing (31) and (40) into (47) can obtain
V ˙ = σ [ μ σ + k ^ p ( ν d z ^ 1 ) + k ^ d ( ν ˙ d z ^ 2 ) ] + χ z [ α χ z + k f z ˙ M ^ g k f z ˙ M g ] + k ˜ p φ 1 k ˜ ˙ p + k ˜ d ζ 1 k ˜ ˙ d + M ˜ g κ 1 M ˜ ˙ g = μ σ 2 + k ^ p ( ν d z ^ 1 ) σ α χ z 2 + k f z ˙ M ^ g χ z + k ^ d ( ν ˙ d z ^ 2 ) σ k f z ˙ M g χ z + k ˜ p φ 1 k ˜ ˙ p + k ˜ d ζ 1 k ˜ ˙ d + M ˜ g κ 1 M ˜ ˙ g μ σ 2 + k ^ p ( ν d z ^ 1 ) σ α χ z 2 + k f z ˙ M ^ g χ z + k ^ d ( ν ˙ d z ^ 2 ) σ k f z ˙ M g χ z k ˜ p φ 1 k ˜ ˙ p + k ˜ d ζ 1 k ˜ ˙ d + M ˜ g κ 1 M ˜ ˙ g
Bringing (41), (42) and (43) into (48) can obtain
V ˙ μ σ 2 + k ^ p ( ν d z ^ 1 ) σ α χ z 2 + k f z ˙ M ^ g χ z + k ^ d ( ν ˙ d z ^ 2 ) σ k f z ˙ M g χ z k ˜ p φ 1 [ δ k ^ p ( ν d z ^ 1 ) φ ] / k ˜ p + k ˜ d ζ 1 [ δ k d ( ν ˙ d z ^ 2 ) ζ ] / k ˜ d + M ˜ g κ 1 κ χ ˙ z k f z ˙
Formula (49) is simplified to obtain
V ˙ μ σ 2 α χ z 2 ,
where μ and α are positive parameters. Simplifying Equation (49) can obtain
V ˙ 0 .
Through the above analysis, according to the Lyapunov stability theory, it can be proved that the system is uniformly stable.
According to Barbalat theorem, we can obtain
lim t χ z ( t ) = 0 , lim t δ ( t ) = 0 .
Therefore, when t , δ ( t ) 0 ; this can prove that the quadrotor system is asymptotically stable.

6. Simulation Results and Discussion

The designed controller is tested through numerical simulation in MATLAB. The following three control methods are employed to carry out comparative simulation experiments on the loaded quadrotor. The initial state selected during simulation is [ 0 0 0 ] rad, [ 0 0 0 ] m. The quadrotor model parameters are shown in Table 1.
Case Study 1.
Quadrotor mass M = 5   k g , load mass m = 3   k g ; given input signal: Z d = 6   m ,   θ d = 35   ,   ϕ d = 45   ,   ψ d = 55   . The simulation results of BP−ALADRC are compared with LADRC and DSADRC. Figure 4 is the output curves, and Figure 5 is the output curves and error curves. In the altitude channel, BP−ALADRC has the fastest tracking speed, BP−ALADRC and DSADRC can reach stability within 1 s and 1.2 s, respectively, while LADRC needs 2.3 s to achieve stable tracking. In the channels of pitch and row angle, while maintaining fast track speed, BP−ALADRC has significantly less overshoot than DSADRC and LADRC. In the yaw angle channel, BP−ALADRC was significantly faster to respond. In summary, the BP−ALADRC controller has the fastest tracking speed and smallest overshoot compared to DSADRC and LADRC, and can quickly and stably reach system stabilization, which proves the effectiveness of the BP−ALADRC controller.
Case Study 2.
In this test, we use BP−LADRC without PD adaptive law as a comparison. To simulate the external interference received by the quadrotor during the working process, we use Gaussian white noise with a variance of 2.5 and a mean value of 0 to represent the external interference. As shown in Figure 6.
According to Figure 7 and Figure 8, when the quadrotor is disturbed, the tracking curve fluctuation of BP−ALADRC to the reference input is the smallest; BP−ALADRC has the fastest tracking speed, while LADRC has the slowest tracking speed; the overshoot of BP−ALADRC is significantly smaller than that of BP−LADRC without PD adaptive law; to sum up, the introduction of BP can reduce the overshoot and speed up the system response speed. At the same time, the addition of the PD adaptive law further increases the stability of the system, eliminates the adverse effects of parameter deviations on the system, further reduces the system overshoot. Meanwhile, system response speed is further accelerated in the height channel, and the control effect of the controller is improved.
From Figure 9, when the quadrotor is affected by external disturbances, the parameters k p and k d of the PD controller will be adjusted timely on the basis of the system state, which can eliminate the parameter deviation.
Case Study 3.
The load mass carried by the quadrotor will change in real time; in this test, we make the load mass decrease continuously, and the system input signal remains unchanged. The BP−ALADRC controller is compared with mass adaptive control law and the BP−ALADRC controller without mass adaptive control law. Figure 10 is a load mass change diagram.
Figure 11 and Figure 12 show that the BP−ALADRC with the mass adaptive law can effectively reduce the error caused by the load mass change, and the BP−ALADRC without the mass adaptive law has significantly larger overshoots in the pitch angle and roll angle channels. At the same time, the response speed of the system is not as good as that of the BP−ALADRC controller with the quality adaptive law.
In Case Study 1, the control scheme designed in this paper and the control methods of LADRC [30] and DSADRC are applied to the attitude control of the quadrotor system. It can be seen that the controller designed in this paper has less overshoot and faster response speed in the dynamic response process. In Case Study 2, we conducted interference experiments. By adding white noise interference to the system, we verified that the controller has a strong ability to suppress interference. Finally, Case Study 3 shows that the controller has strong adaptability to the change of load mass. All the above case studies show that the designed controller can track the given signal quickly and stably, effectively reducing the adverse effects of interference on the quadrotor system. It is proved that the controller designed in this paper has more significant signal tracking ability, interference suppression ability and better system adjustment ability. In addition, the introduction of adaptive control effectively compensates the observation error of LESO, making the controller more robust. Finally, when the load mass of the quadrotor changes, the designed mass adaptive law can effectively overcome the adverse effects of mass changes on the system and enhance the stability of the system.

7. Conclusions

In order to improve the attitude control performance of the quadrotor, this paper considers the change of the load mass of the quadrotor, establishes a dynamic model of the quadrotor with variable mass, and designs a composite control scheme based on LADRC. For the problems of many parameters of LADRC controller, difficult adjustment and large error of manual parameter adjustment, this paper proposes BP neural network and adaptive control to optimize and adjust LADRC parameters. The BP neural network is used to dynamically optimize and adjust the bandwidth value of LESO in real time, which improves the disturbance observation accuracy of the system. At the same time, the introduction of adaptive control reduces the observation error of LESO and the parameter deviation of PD controller. These ensure the stability of the quadrotor system in the complex and changeable working environment. Finally, all the case studies show that for the same given attitude angle and altitude, the designed control system has the advantages of small overshoot, fast system response and strong robustness compared with LADRC and DSADRC controllers. In addition, the designed mass adaptive law can also effectively reduce the adverse effects of load mass changes on the system and further enhance the stability of the system. At present, due to the limitations of experimental conditions, this paper only carried out theoretical research. In the future work, we will verify the effectiveness of the control scheme in practical applications.

Author Contributions

Y.G.: contributed to the conception of the study, performed the experiment, performed the data analyses and wrote the manuscript; G.Z.: contributed to the conception of the study; T.Z.: helped perform the analysis with constructive discussions. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

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

  1. Comba, L.; Biglia, A.; Aimonino, D.R.; Gay, P. Unsupervised detection of vineyards by 3D point-cloud UAV photogramme-try for precision agriculture. Comput. Electron. Agric. 2018, 155, 84–95. [Google Scholar] [CrossRef]
  2. Tokekar, P.; Vander Hook, J.; Mulla, D.; Isler, V. Sensor planning for a symbiotic UAV and UGV system for precision agriculture. IEEE Trans. Robot. 2016, 32, 1498–1511. [Google Scholar] [CrossRef]
  3. Roberge, V.; Tarbouchi, M.; Labonté, G. Fast Genetic Algorithm Path Planner for Fixed-Wing Military UAV Using GPU. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 2105–2117. [Google Scholar] [CrossRef]
  4. Shakeel, T.; Arshad, J.; Jaffery, M.H.; Rehman, A.U.; Eldin, E.T.; Ghamry, N.A.; Shafiq, M. A Comparative Study of Control Methods for X3D Quadrotor Feedback Trajectory Control. Appl. Sci. 2022, 12, 9254. [Google Scholar] [CrossRef]
  5. Islam, M.; Okasha, M.; Sulaeman, E. A model predictive control (mpc) approach on unit quaternion orientation based quadrotor for trajectory tracking. Int. J. Control. Autom. Syst. 2019, 17, 2819–2832. [Google Scholar] [CrossRef]
  6. Eskandarpour, A.; Sharf, I. A constrained error-based MPC for path following of quadrotor with stability analysis. Nonlinear Dyn. 2020, 99, 899–918. [Google Scholar] [CrossRef]
  7. Falcon, R.; Rios, H.; Dzul, A. Comparative analysis of continuous sliding-modes control strategies for quad-rotor robust tracking. Control. Eng. Pract. 2019, 90, 241–256. [Google Scholar] [CrossRef]
  8. Okyere, E.; Bousbaine, A.; Poyi, G.T.; Joseph, A.K.; Andrade, J.M. Lqr controller design for quad-rotor helicopters. J. Eng. 2019, 17, 4003–4007. [Google Scholar] [CrossRef]
  9. Jian, P.; Changlong, L. UAV attitude control with LQR controller based on extended state observer. Syst. Simul. 2018, 30, 753–759. [Google Scholar]
  10. Rodríguez-Guerrero, L.; Benítez-Morales, A.; Santos-Sánchez, O.-J.; García-Pérez, O.; Romero-Trejo, H.; Ordaz-Oliver, M.-O.; Ordaz-Oliver, J.-P. Robust Backstepping Control Applied to UAVs for Pest Recognition in Maize Crops. Appl. Sci. 2022, 12, 9075. [Google Scholar]
  11. Nigro, M.; Pierri, F.; Caccavale, F. Control of an Omnidirectional UAV for Transportation and Manipulation Tasks. Appl. Sci. 2021, 11, 10991. [Google Scholar] [CrossRef]
  12. Fan, Y.; Guo, H.; Han, X.; Chen, X. Research and Verification of Trajectory Tracking Control of a Quadrotor Carrying a Load. Appl. Sci. 2022, 12, 1036. [Google Scholar] [CrossRef]
  13. Chen, X.; Zhao, Y.; Fan, Y. Adaptive Integral Backstepping Control for a Quadrotor with Suspended Flight. In Proceedings of the 2020 5th International Conference on Automation. Control and Robotics Engineering (CACRE), Dalian, China, 19–20 September 2020; pp. 226–234. [Google Scholar]
  14. Han, J.Q. From PID to Active Disturbance Rejection Control. IEEE Trans. Ind. Electron. 2009, 56, 900–906. [Google Scholar] [CrossRef]
  15. Zhang, Q.; Wei, Y.; Li, X. Quadrotor Attitude Control by Fractional-Order Fuzzy Particle Swarm Optimization-Based Active Disturbance Rejection Control. Appl. Sci. 2021, 11, 11583. [Google Scholar] [CrossRef]
  16. Yuan, C.; Zhou, X.; Ma, Y. DC Bus Voltage Control of Wind Power Inverter Based on First-Order LADRC. IEEE Access 2021, 10, 3263–3274. [Google Scholar] [CrossRef]
  17. Zhao, F.; Wang, X.; Zhu, T. Low-Frequency Passivity-Based Analysis and Damping of Power-Synchronization Controlled Grid-Forming Inverter. IEEE J. Emerg. Sel. Top. Power Electron. 2022. [Google Scholar] [CrossRef]
  18. Hua, L.; Zhang, J.; Li, D.; Xi, X. Fault-Tolerant Active Disturbance Rejection Control of Plant Protection of Unmanned Aerial Vehicles Based on a Spatio-Temporal RBF Neural Network. Appl. Sci. 2021, 11, 4084. [Google Scholar] [CrossRef]
  19. Guo, B.; Bacha, S.; Alamir, M.; Mohamed, A.; Boudinet, C. LADRC applied to variable speed micro-hydro plants: Experimental validation. Control. Eng. Pract. 2019, 85, 290–298. [Google Scholar] [CrossRef]
  20. Zhou, Y.; Luo, D.; Wu, B.; Cheng, B.; Lin, Q. Active vibration isolation system based on the LADRC algorithm for atom interferometry. Appl. Opt. 2020, 59, 3487–3493. [Google Scholar] [CrossRef]
  21. Geng, X.; Zhu, Q.; Liu, T.; Na, J. U-model based predictive control for nonlinear processes with input delay. J. Process Control. 2019, 75, 156–170. [Google Scholar] [CrossRef] [Green Version]
  22. Wang, H.; Huang, H. Property and applications of extended state observer. Control. Decis. 2013, 28, 1078–1082. [Google Scholar]
  23. Zhang, Y.; Chen, Z.; Sun, M. Trajectory tracking control for a quadrotor unmanned aerial vehicle based on dynamic surface active disturbance rejection control. Trans. Inst. Meas. Control. 2020, 42, 2198–2205. [Google Scholar] [CrossRef]
  24. Wu, J.; Peng, H.; Chen, Q.; Penge, X. Modeling and control approach to a distinctive quadrotor helicopter. ISA Trans. 2014, 53, 173–185. [Google Scholar] [CrossRef] [PubMed]
  25. Cai, Z.; Lou, J.; Zhao, J.; Wu, K.; Liu, N.; Wang, Y.X. Quadrotor trajectory tracking and obstacle avoidance by chaotic grey wolf optimization-based active disturbance rejection control. Mech. Syst. Signal Process. 2019, 128, 636–654. [Google Scholar] [CrossRef]
  26. 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]
  27. Tang, P.; Lin, D.; Zheng, D.; Fan, S.; Ye, J. Observer based finite-time fault tolerant quadrotor attitude control with actuator faults. Aerosp. Sci. Technol. 2020, 104, 105968. [Google Scholar] [CrossRef]
  28. Nekoukar, V.; Dehkordi, N.M. Robust path tracking of a quadrotor using adaptive fuzzy terminal sliding mode control. Control. Eng. Pract. 2021, 110, 104763. [Google Scholar] [CrossRef]
  29. Wang, B.; Yu, X.; Mu, L.; Zhang, Y. Disturbance observer-based adaptive fault-tolerant control for a quadrotor helicopter subject to parametric uncertainties and external disturbances. Mech. Syst. Signal Process. 2019, 120, 727–743. [Google Scholar] [CrossRef]
  30. Liang, X.; Li, J.; Zhao, F. Attitude Control of Quadrotor UAV Based on LADRC Method. In Proceedings of the 2019 Chinese Control and Decision Conference (CCDC), Nanchang, China, 3–5 June 2019; pp. 1924–1929. [Google Scholar]
Figure 1. Structure diagram of a quadrotor.
Figure 1. Structure diagram of a quadrotor.
Applsci 12 12698 g001
Figure 2. Control scheme structure diagram.
Figure 2. Control scheme structure diagram.
Applsci 12 12698 g002
Figure 3. The structure of BP.
Figure 3. The structure of BP.
Applsci 12 12698 g003
Figure 4. System output: (a) Altitude and Roll angle output; (b) Pitch and Yaw angle output.
Figure 4. System output: (a) Altitude and Roll angle output; (b) Pitch and Yaw angle output.
Applsci 12 12698 g004
Figure 5. System output errors: (a) Altitude and Roll angle output errors; (b) Pitch and Yaw angle output errors.
Figure 5. System output errors: (a) Altitude and Roll angle output errors; (b) Pitch and Yaw angle output errors.
Applsci 12 12698 g005
Figure 6. Gaussian noise simulation diagram.
Figure 6. Gaussian noise simulation diagram.
Applsci 12 12698 g006
Figure 7. System output: (a) Altitude and Roll angle output; (b) Pitch and Yaw angle output.
Figure 7. System output: (a) Altitude and Roll angle output; (b) Pitch and Yaw angle output.
Applsci 12 12698 g007aApplsci 12 12698 g007b
Figure 8. System output errors: (a) Altitude and Roll angle output errors; (b) Pitch and Yaw angle output errors.
Figure 8. System output errors: (a) Altitude and Roll angle output errors; (b) Pitch and Yaw angle output errors.
Applsci 12 12698 g008
Figure 9. k p and k d adaptive values: (a) k p adaptive value; (b) k d adaptive value.
Figure 9. k p and k d adaptive values: (a) k p adaptive value; (b) k d adaptive value.
Applsci 12 12698 g009
Figure 10. Load mass variation diagram.
Figure 10. Load mass variation diagram.
Applsci 12 12698 g010
Figure 11. System output: (a) Altitude and Roll angle output; (b) Pitch and Yaw angle output.
Figure 11. System output: (a) Altitude and Roll angle output; (b) Pitch and Yaw angle output.
Applsci 12 12698 g011
Figure 12. System output errors: (a) Altitude and Roll angle output errors; (b) Pitch and Yaw angle output errors.
Figure 12. System output errors: (a) Altitude and Roll angle output errors; (b) Pitch and Yaw angle output errors.
Applsci 12 12698 g012
Table 1. Parameters of the quadrotor model.
Table 1. Parameters of the quadrotor model.
SymbolValue
k 0.21   N s 2 / r a d 2
M 5   k g
I x 0.8   k g m 2
I y 0.8   k g m 2
I z 1.5   k g m 2
g 9.8   m / s 2
l 0.15   m
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/app122412698

AMA Style

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. Applied Sciences. 2022; 12(24):12698. https://doi.org/10.3390/app122412698

Chicago/Turabian Style

Gao, Yunlong, Guixin Zhu, and Tong Zhao. 2022. "Based on Backpropagation Neural Network and Adaptive Linear Active Disturbance Rejection Control for Attitude of a Quadrotor Carrying a Load" Applied Sciences 12, no. 24: 12698. https://doi.org/10.3390/app122412698

APA Style

Gao, Y., Zhu, G., & Zhao, T. (2022). Based on Backpropagation Neural Network and Adaptive Linear Active Disturbance Rejection Control for Attitude of a Quadrotor Carrying a Load. Applied Sciences, 12(24), 12698. https://doi.org/10.3390/app122412698

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