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 and 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
,
,
, 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
[
25,
26].
The positional motion of the quadrotor satisfies Newton’s second law,
where
;
is the linear velocity;
is the mass of the quadrotor;
is the total lift;
;
;
is the matrix of air resistance.
The lift force of the quadrotor is expressed as
where
; from
, we can obtain
:
In summary, we can obtain the positional dynamics equation of the quadrotor.
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
, the coordinates of the centroid of the part of the load that is reduced during operation are
. When the load mass decreases, the change in the coordinates of the centroid of quadrotor is shown in the following formula:
where the coordinates of the centroid of the quadrotor are
. According to
Figure 1, initial coordinate of the centroid of the quadrotor is
, the coordinate of the centroid of the reduced load mass during the working process is
. Substituting it into Formula (6), the coordinate of the centroid of the quadrotor can be obtained as
.
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]
where
is the inertia matrix;
is the angular velocity;
is the resultant moment;
and
are the pulling torque and reaction torque of the four-wheel drive motor, respectively;
is the matrix of air resistance.
In summary, we can obtain attitude dynamics equation of the quadrotor.
where
,
, and
are the yaw angle, roll angle and pitch angle, respectively.
To simplify the quadrotor dynamic equation, (
,
,
,
) is used to represent the relevant parameters.
In summary, we can obtain the dynamic equation of the quadrotor.
Equation (13) can also be written as [
30]
where
,
;
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, is the given signal, is the observation vector of LESO to the system state variable, is the bandwidth gain vector of LESO, is the approximate value of the given signal and its differential, is a nonlinear state error feedback control law with PD controller output, is the control law of the system, is the error combination, are the estimated values of PD controller parameters , respectively, is the estimated value of load mass change , is the state variable of the system, 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
where
are measurable state variables;
is the output signal;
is an unknown nonlinear function of the system;
is the external disturbance.
a is partially known, and the known part is denoted as
.
All uncertainties in the system and all disturbances received by the system are regarded as the total disturbance, which is defined as
To observe the total disturbance, the extended state space is defined , where , and .
The nonlinear system (12) can be transformed into
where
,
,
,
.
According to (16), LESO is designed as follows:
where
is the state vector of the observer, which is used to observe the state of system
;
is the observation gain vector of LESO.
The gain value of LESO is obtained from its characteristic equation:
where
is the bandwidth of LESO, denoted as
.
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
,
,
. The control law can be written in the following form:
where
is the output of the PD controller.
4.2. Parameter Tuning Adaptive Law Design
By applying APC to dynamically adjust the and 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, and avoid errors caused by manual parameter adjustment, we introduce an adaptive law which can adjust, timely. Formulas (36) and (37) provide the adaptive law of, .
Assumption 1: Assume that a given signal is smooth and bounded such that is true. Here is a known constant.
Derivation can be obtained as follows:
Define the filter tracking error:
where
is a suitable gain vector such that when
,
is satisfied.
Derivation of
can be obtained:
From the output of LESO, we can obtain
The input of the control is as follows:
where
is any positive number.
In summary, we can obtain
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 theandvalues 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
where
is a given signal,
is the state variable of the system; derivation of
can be obtained in the following way:
The input of BP neural network can be expressed as
The input and output of the hidden layer are
The hidden unit is 6, and the excitation function adopts the tangent sigmoid transfer function:
Remark 4. BP neural network adjusts the bandwidth gain vector 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 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
, the differential signal
of the output error and the bias value
. 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):
Since the BP neural network will normalize the input signal
, its output signal will also be mapped to the range
, 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
of LESO in the BP−ALADRC controller. The gain vector of LESO can be obtained as
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
According to the equivalent deterministic principle, we can obtain
where
is an adjustable parameter. Bringing (37) into (36) can obtain
Then,
,
are defined, then the estimation error can be expressed as
Substituting Formula (39) into Formula (38) can obtain
5. Stability Analysis
To ensure that the system tracking error converges to the zero neighborhood, we design the following adaptive law:
Remark 5. To avoidandbeing zero,
andare integrated. The above formula can be rewritten as Appropriate constants and are chosen to make and non-zero.
Proof of Theorem 1. A positive definite Lyapunov function is designed as follows.
Taking the derivative of Equation (46),
Bringing (31) and (40) into (47) can obtain
Bringing (41), (42) and (43) into (48) can obtain
Formula (49) is simplified to obtain
where
and
are positive parameters. Simplifying Equation (49) can obtain
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
Therefore, when , ; 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
rad,
m. The quadrotor model parameters are shown in
Table 1.
Case Study 1. Quadrotor mass
, load mass
; given input signal:
. 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
and
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.