Next Article in Journal
Detection of White Leaf Disease in Sugarcane Using Machine Learning Techniques over UAV Multispectral Images
Next Article in Special Issue
UAV Fault Detection Methods, State-of-the-Art
Previous Article in Journal
ARSD: An Adaptive Region Selection Object Detection Framework for UAV Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional-Order Linear Active Disturbance Rejection Control Design and Optimization Based Improved Sparrow Search Algorithm for Quadrotor UAV with System Uncertainties and External Disturbance

School of Astronautics, Beihang University (BUAA), Beijing 100191, China
*
Author to whom correspondence should be addressed.
Drones 2022, 6(9), 229; https://doi.org/10.3390/drones6090229
Submission received: 30 July 2022 / Revised: 27 August 2022 / Accepted: 30 August 2022 / Published: 31 August 2022

Abstract

:
This paper presents a generic robust controller that allows applications in various quadrotor unmanned aerial vehicle (UAV) systems effectively even when facing severe system uncertainties and unknown external disturbances. The fractional-order linear active disturbance rejection control (FOLADRC) scheme has combined the advantages of the fractional-order PID (FOPID) with the linear active disturbance rejection control (LADRC). Firstly, the structure of the FOLADRC-based quadrotor UAV is designed. Then, considering the difficulty of parameter tuning of FOLADRC and the demand for accuracy and rapidity of the controller, the improved sparrow search algorithm is applied. Finally, to illustrate the robustness and effectiveness of FOLADRC, the FOLADRC-based quadrotor UAV is firstly compared with PID and LADRC. The simulation and experiment results show that the FOLADRC method can suppress the influence of system uncertainties and external disturbance effectively, where the superiority compared to PID and LADRC has been demonstrated clearly.

1. Introduction

The quadrotor UAV is widely used in the civil and military field owing to the advantages of hovering, vertical taking off, and landing [1,2]. With the development of the quadrotor UAV, different payloads are carried in it, such as the camera and the mechanical arm, and applied scenarios become more complex with severe system uncertainties and external disturbances [3]. Therefore, a robust and precise quadrotor UAV controller is essential to suppress the effect of complex environments.
Among all control methods, the PID and its variants are the most popular methods for the flight control of the quadrotor UAV, ascribed to implementation simplicity and effectiveness. In some simple application scenarios, the PID is efficient enough to satisfy the task requirements [4,5]. However, some research demonstrates that the performance of PID decreases dramatically in complex scenarios [6,7]. Hence, it is still urgent to design a robust controller to overcome the large payload variations, external disturbances, and system uncertainties.
In recent years, many control methods have been proposed to enhance quadrotor UAV performance, such as sliding mode control (SMC) [8], model predictive control (MPC) [9], and backstepping control [10]. Nevertheless, most methods require relatively accurate system dynamics, which is complicated and time-consuming in a real application. Thus, a simplified model without the consideration of higher-order dynamics is usually established with an estimated value of quadrotor UAV parameters [11]. However, the inaccuracy in quadrotor UAVs significantly influences the anti-disturbance ability of controllers.
Therefore, it is essential to eliminate the bad effect of system inaccuracy. For this purpose, the active disturbance rejection control (ADRC) proposed by Han fully inherits the merits of PID and has strong robustness [12]. The basic idea of ADRC is to estimate and compensate for the total disturbance between the real system and the ideal system. Consequently, the system would become an ideal system simple for controller design if the total disturbance generated by system uncertainties and external disturbance is compensated. The ADRC has been verified to have better performance in disturbance rejection, including the disturbance from parameters perturbation and environment [13,14]. However, the physical meaning of the parameters of ADRC is fuzzy, which is not conducive to practical application. Given this complication, Gao et al. designed the linear active disturbance rejection control (LADRC), which changed the nonlinear controller to PD. The LADRC remarkably decreases the number of tuning parameters while presenting the same control performance compared with ADRC [15]. Recent research suggested that LADRC has a satisfactory performance on quadrotor UAVs [16,17,18].
Although the LADRC could resist the negative influence of system uncertainties and external disturbance, it is relatively sensitive to the estimation of system gain b. Therefore, on the basis of LADRC, FOLADRC is developed by replacing PD with fractional-order PID (FOPID) to solve the problem. The concept of fractional-order PID (FOPID) was first proposed by Podlubny, with the addition of two fractional calculus operations (λ and μ) to the classical PID controller [19]. By changing the parameters of λ and μ, the FOPID can make the robustness of the open-loop gain change in the system zero or close to zero. At this time, the phase angle Bode diagram of the corresponding open-loop transfer function is flat near the cutoff frequency. When the gain b of the system changes, the open-loop cut-off frequency of the system will also change, but the corresponding phase margin of the system only slightly changes, so it shows that FOPID has strong robustness to the change in control gain [20]. Hence, FOLADRC possesses the advantages of both FOPID and LADRC, equipping it with great potential in controller design.
FOLADRC has been widely employed in some fields due to its simple structure and excellent dynamic performance. FOLADRC was applied to the speed control of nonlinear two-mass systems [21]. The results reveal that the controller has great performance with the perturbations. By designing the FOLADRC controller for an active power filter, the current is compensated, and the voltage tracked is achieved [22]. Li et al., investigated FOLADRC schemed for the underdrive commensurate FOS [23]. Moreover, FOLADRC was designed for the heat-flow experiment, which is the integer-order system. The results verify that the controller has great robustness in both FOS and the integer-order system [24]. However, there is no research on the real implementation in quadrotor UAVs. Only a few studies have adopted FOPID to control the attitude of the quadrotor, whose control performance is significantly improved compared with PID [25].
Additionally, the parameter optimization method is vague and only adjusts according to the approximate impact of controller variables on the system. Thus, the main challenge before the FOLADRC is employed in the quadrotor UAV is finding an efficient method of tuning parameters. Gao et al. proposed the performance parameter tuning method based on the frequency domain method by designing the constraint equation of the frequency domain index of the ADRC controller [26]. Unfortunately, it cannot completely guarantee good performance in the time domain. With the purpose of tracking the problem, research has been conducted with evolutionary algorithms to tune the parameters, where the system performance is strongly connected to the objective functions. Ge et al. offered a genetic algorithm (GA) to tune the parameters of ADRC, revealing that the proposed algorithm has better performance in error but slow convergence [27]. He et al. presented a multi-strategy pigeon-inspired optimization (MSPIO) approach to adjust its parameters. Its convergence rate is improved, while the issue of control accuracy is not managed [28]. The sparrow search algorithm (SSA) is also an evolutionary algorithm [29]. SSA has been successfully used in different fields owing to its fast convergence speed [30,31,32]. Nonetheless, it lacks an excellent mutation mechanism and is prone to falling into the local optimum.
Motivated by the above discussions, FOLADRC based on an improved SSA is proposed in this paper. The main contributions of this paper are as follows:
(a)
A realistic and nonlinear model with disturbance and parameter uncertainty of quadrotors is established.
(b)
The FOLADRC is successfully applied to the quadrotors system to overcome the disturbance and parameter uncertainty during flight. The simulation and experiment prove the superior control performance of the method.
(c)
Considering that the traditional SSA optimization algorithm easily falls into the local optimum, this paper develops an improved SSA based on the Levy flight (L-SSA) to address the issue of parameter tuning and optimization in the controller.
The rest of the paper is organized as follows. Section 2 illustrates the dynamic model of a quadrotor UAV. In Section 3, the proposed FOLADRC for attitude control is detailed. Section 4 describes the L-SSA and its implementation on the FOLADRC. Next, the pitch angle tracking simulation and experimental results on the platform are presented in Section 5. Finally, the conclusions of this paper are drawn in Section 6.

2. Dynamic Model of Quadrotor UAV

According to the Newton–Euler equations, the quadrotor UAV is established as six-degree-of-freedom dynamics. The experiment quadrotor UAV is shown in Figure 1.
It is supposed that the quadrotor UAV has a rigid body and is structurally symmetric. For analysis, two coordinate systems, namely the inertial frame (xn, yn, zn) and the body frame (xb, yb, zb), as shown in Figure 2, are considered.
The body coordinate system can coincide with the inertial coordinate system by rotating around the coordinate axis three times, and the angles of the three rotations are the roll angle, pitch angle, and yaw angle ( φ , θ, ψ), respectively. According to the transformation order of X-Y-Z, the transformation matrix R b n from the body coordinate system to the ground coordinate system is obtained:
R b n = [ c o s   θ   c o s   ψ   s i n   θ   c o s   ψ   s i n   φ s i n   ψ   c o s   φ s i n   θ   c o s   ψ   c o s   φ + s i n   ψ   s i n   φ c o s   θ   s i n   ψ   s i n   θ   s i n   ψ   s i n   φ + c o s   ψ   c o s   φ s i n   θ   s i n   ψ   c o s   φ c o s   ψ   s i n   φ s i n   θ   c o s   θ   s i n   φ c o s   θ   c o s   φ ]
The conversion between the body coordinate system and the ground coordinate system satisfies the following equation:
C n = R b n C b
where   C n and C b represent the coordinates under the inertial coordinate system and the body coordinate system, respectively.

2.1. UAV Translation Model

In the body coordinate system, the lift provided by the quadrotor UAV is always along the z-axis. The conversion of lift from the body coordinate system to the inertial coordinate system can be expressed:
F n = R b n F b = R b n [ 0   0   i = 1 4 F i ] T = i = 1 4 F i [ s i n   θ   c o s   ψ   c o s   φ + s i n   ψ   s i n   φ s i n   θ   s i n   ψ   c o s   φ c o s   ψ   s i n   φ c o s   θ   c o s   φ ]
where F i = k i ω i 2 is the lift provided by the ith (i = 1,2,3,4) rotor, ki is the lift coefficient of the quadrotor propeller, and ω i represents the speed of the ith rotor.
In the inertial coordinate system, the air resistance opposite the flight direction suffered by the quadrotor during flight is approximately regarded as proportional to the square of the speed:
f n = [ D x x ˙ 2   D y y ˙ 2   D z z ˙ 2 ] T
where   D x , D y , and   D z   represent the air resistance coefficients in the three coordinate axes of the inertial coordinate system, respectively.
According to Newton’s second law, the translational model of the quadrotor UAV is obtained as follows:
{ x ¨ = 1 m i = 1 4 F i ( s i n   θ   c o s   ψ   c o s   φ + s i n   ψ   s i n   φ ) 1 m D x x ˙ 2 y ¨ = 1 m i = 1 4 F i ( s i n   θ   s i n   ψ   c o s   φ c o s   ψ   s i n   φ ) 1 m D y y ˙ 2 z ¨ = 1 m i = 1 4 F i c o s   θ   c o s   φ g 1 m D z z ˙ 2

2.2. UAV Rotation Model

According to the theorem of momentum, the rotation equation of the quadrotor UAV in the inertial coordinate system satisfies:
M = d d t ( J Ω ) = J Ω ˙ + Ω × J Ω
where M = [ M x b , M y b , M z b ] T , Ω = [ p , q , r ] T , and   J = d i a g [ J x , J y , J z ] . M represents the total torque of the quadrotor, Ω   represents the body rotational angular velocity vector, and J is the moment of inertia matrix of the UAV.
The torque of the UAV is mainly composed of the aerodynamic torque provided by the rotor and the gyro torque generated by the high-speed rotating. The gyro moment can be expressed as:
M gro   = j r [ q p 0 ] ( ω 1 + ω 2 ω 3 + ω 4 ) = [ j r q ( ω 1 + ω 2 ω 3 + ω 4 ) j r p ( ω 1 + ω 2 ω 3 + ω 4 ) 0 ]
where j r represents the moment of inertia of the motor and propeller.
Based on the mechanical formula, the UAV under the ‘X’ shape distribution can obtain its moment component along the three axes:
{ M x b = l ( F 2 + F 3 F 1 F 4 ) + j r q ( ω 1 + ω 2 ω 3 + ω 4 ) M y b = l ( F 3 + F 4 F 1 F 2 ) j r p ( ω 1 + ω 2 ω 3 + ω 4 ) M z b = 2 l C M ( ω 2 2 + ω 4 2 ω 1 2 ω 3 2 )
where l is the distance from the center of gravity of the body to the connecting line between the centers of the two rotors, and C M represents the anti-torque coefficient of the propeller.
The relationship between the attitude angular rate of the quadrotor UAV and the 3-axis angular rate in the body frame coordinate system is:
[ p q r ] = [ φ ˙ ψ ˙ s i n   θ   θ ˙ c o s   φ + ψ ˙ s i n   φ c o s   θ   θ ˙ s i n   φ + ψ ˙ c o s   φ c o s   θ   ]
This paper studies the small change in the attitude angle of the UAV, so the above formula can be simplified as follows:
{ φ ˙ = p θ ˙ = q ψ ˙ = r
Upon substituting Equations (8) and (10) into Equation (6), the rotation mathematical model of quadrotor UAV is finally obtained:
[ φ .. θ .. ψ .. ] = [ ( J y J z J x ) θ ˙ ψ ˙ + j r J x θ ˙ ( ω 1 + ω 2 ω 3 + ω 4 ) + 1 J x l ( F 2 + F 3 F 1 F 4 ) ( J z J x J y ) φ ˙ ψ ˙ j r J y φ ˙ ( ω 1 + ω 2 ω 3 + ω 4 ) + 1 J y l ( F 3 + F 4 F 1 F 2 ) ( J x J y J z ) φ ˙ θ ˙ + 1 J z 2 l C M ( ω 2 2 + ω 4 2 ω 1 2 ω 3 2 ) ]

3. FOLADRC Design for Quadrotor UAV

Accurate attitude control is crucial for the quadrotor UAV to achieve high flight performance. Therefore, FOLADRC is applied to enhance the robustness and rapidity of the quadrotor UAV. We assume that the structure of the quadrotor UAV is symmetrical, therefore the pitch and roll channels have the same control structure. Besides, the yaw controller is designed with the same process as horizontal control. Consequently, only the control scheme of the pitch channel is discussed in this paper. According to Equation (11), the dynamic model of the pitch channel can be simplified as follows:
{ θ ˙ 1 = θ 2 θ ˙ 2 = f ( φ , φ ˙ , θ   , θ ˙ , ψ   , ψ ˙ ) + ω ( t ) + b U y = θ 1
where θ 1 represents the pitch angle of the system, ω ( t ) represents the external disturbance generated by the environment such as gusty wind, and   f ( φ , φ ˙ , θ   , θ ˙ , ψ   , ψ ˙ ) represents the nominal model. U is the pitching moment acting on the quadrotor UAV. b = ( b 0 + b ) represents the controller gain, where b 0 is the estimated value of b and b is the system uncertainty parameter. Since the dynamic of pitch channel is a second-order system, a second-order FOLADRC controller can be applied to control it. FOLADRC consists of the tracking-differentiator (TD), the FOPID controller, and the linear extended state observer (LESO). The structural diagram of the FOLADRC-based system is shown in Figure 3.
(1)
Tracking-differentiator
The input signal, called the desired pitch angle, is supposed to be r, and the equation of TD can be given as:
{ e = θ 1 r f h = fhan ( e , θ 2 , r 0 , h ) θ 1 = θ 1 + h θ 2 θ 2 = θ 2 + h f h
where e represents the pitch error. h and r0 are the filtering factor and speed factor, respectively. The fhan ( e , θ 2 , r 0 , h ) is the optimal control synthetic function and it can be expressed as:
{ d = r 0 h d 0 = r 0 h 2 y = e + h θ 2 a 0 = d 2 + 8 r 0 | y | a = { θ 2 + a 0 d 2 s i g n ( y )   ,   | y | > d 0   θ 2 + y h   ,         | y |   d 0           f h a n = { r 0 s i g n ( a ) , | a | > d r 0 a d ,                                         | a | d
(2)
Error Feedback
Compared with the traditional PID, FOPID introduces two fractional calculus operators, λ and μ, to improve the control performance. Generally, the FOPID controller extends the traditional PID into the fractional-order field, which inherits the advantages of the PID and has stronger robustness. The transfer function of FOPID is described as:
G c ( s ) = K p + K i s λ + K d s μ
where 0 < λ < 2 and 0 < μ < 2. To enhance the dynamic property of the system in the mid-frequency range and simplify parameter tuning, Equation (15) can be rewritten as:
G c ( s ) = K p ( 1 + K i K p s λ ) ( 1 + K d K p s μ )
The structural diagram of FOPID is shown in Figure 4.
(3)
Linear Extended state observer
After estimating the total disturbance by LESO, the disturbance in a system can be compensated. Thus, the system has a strong anti-disturbance ability. The equation based on Equation (12) is established as follows:
{ e 1 = z 1 y z ˙ 1 = z 2 β 1 e 1 z ˙ 2 = z 3 + b 0 u β 2 e 1 z ˙ 3 = β 3 e 1
where z1, z2, and z3 are the estimated values of the pitch angle, the rate of change of the pitch angle, and total disturbance, respectively. β 1 ,   β 2 , and β 3 are the observer gains that need to be tuned, and u is the control input.
Since this paper focuses on the control performance of FOLADRC, taking K p as the reference, the controller bandwidth ω c   = K p   is calculated according to the bandwidth method [15]. Then, the poles of LESO are assigned to (− ω o ,0) according to the observer bandwidth ω o , and ω o = 5 ω c .
In summary, the design of FOLADRC has been completed.

4. Parameter Tuning of FOLADRC Based on L-SSA

4.1. L-SSA

The sparrow population consists of producers, scroungers, and sentinels. Figure 5 shows the relationship between the three types of sparrows. The producers are responsible for finding food for the entire sparrow population and providing foraging directions for all scroungers. The scroungers follow the producers and compete with them for food resources or forage around them. The sentinels will issue an alarm and move to a safe location after discovering the natural enemies.
X i , j N + 1 = { X i , j N + 1 e x p ( i / ( γ N m a x ) ) , R 2 < S t X i , j N + 1 + Q L , R 2 S t
where N represents the current iteration times, X i , j represents the position information of the ith sparrow in the jth dimension, and L is an all-one matrix of dimension 1 × d. Random number R2 ∈ [0, 1], γ ∈ (0, 1), Q ~ N(0, 1), and St is a safe value.
The scroungers search for the producers with the highest energy, and their updated location can be represented as:
X i , j N + 1 = { Q e x p [ ( X w o r s t N X i , j N ) / i 2 ] , ( i > n / 2 ) X b e s t N + 1 + | X i , j N X b e s t N + 1 | A + L , ( i n / 2 )
where X b e s t N + 1 is the optimal position occupied by the producers after the position update, X w o r s t N is the worst position in the current global, and A is a real matrix with values ± 1. When i > n / 2 , the ith hungry scrounger migrates to other regions for foraging to obtain more energy because of having no food.
After discovering the natural enemies, the sentinels at the edge of the population will quickly move to a safe area, and the sentinels in the middle of the population will approach other sparrows by walking randomly. The updated location of sentinels can be described as:
X i , j N + 1 = { X b e s t N + 1 + λ | X i , j N X b e s t N | , ( f i t f , i > f i t f , g ) [ X i , j N + K | X i , j N X w o r s t N | ( f i t f , i f i t f , w + ε ) ] , ( f i t f , i = f i t f , g )
where random number K∈[0, 1], λ ~ N(0, 1), and ε is a small quantity to remove the effect of a possible 0 denominator. f i t f , g and f i t f , w represent the minimum and maximum values of the fitness function, respectively.
As mentioned in the introduction, SSA lacks an excellent mutation mechanism and is prone to fall into local optimum. Consequently, the “Levy flight” strategy of nature is introduced. It uses the random walk strategy to generate large jumps and sharp direction changes in the process of position updates, to avoid the sparrow individuals being bound by local extreme values, and at the same time to open up the foraging search space, which effectively improves the high-dimensional optimization effect of the sparrow algorithm.
The Levi mutation is represented as follows:
X i i t + 1 = { X i i t + ( X i i t X b e s t i t ) L e v y ( b e t a )
where L e v y ( )   is the Levi mutation factor and 1 b e t a 3 .
The flowchart of L-SSA is shown in Figure 6, and the details of the optimization process are as follows:
Step 1. Initialize parameters of L-SSA, such as the sparrow population n, the proportion of producers Pt, iterations Nmax, etc. Set the position of each sparrow randomly and the number of initial iterations as 0.
Step 2. Calculate the fitness of each sparrow to obtain the current fitness limit value and the position of the corresponding sparrow.
Step 3. The n × Pt sparrows with the least fitness are selected as the producers, and the remaining sparrows are scroungers, and the positions are updated according to Equations (18) and (19). Then, sentinels are randomly generated, and their positions are updated according to Equation (20). Calculate the fitness of each sparrow after the location update.
Step 4. Calculate the inertia weight factor η = 1 N N m a x , and adopt the tournament rules to extract sparrows to participate in the tournament. The least fit sparrows are selected for the Levi flight, and the positions are updated according to Equation (21).
Step 5. If the current global fitness minimum value f i t f , i is not greater than the minimum optimal fitness f i t f , g , the optimization is stopped, otherwise step 5 continues.
Step 6. If the iteration time N is equal to the N m a x , the optimization is stopped. If not, reorder and update the population class, and finally loop into step 3.

4.2. Fitness Function and Parameters Tuning

To obtain better control performance, the fitness function based on the time domain index and frequency domain index is set as:
f i t f ( X ) = ( 5 γ 1 0 5 t | e ( X , t ) | d t + 30 γ 2 s i g n ( P m ( X ) P m c m d ) + 10 γ 3 s i g n ( ω f ( X ) ω f c m d ) γ 4 K p )  
where X(X= [ K p , K i , K d , λ , μ ]) represents the optimized parameter variable, t represents the system simulation time, and e(t) represents the pitch angle command tracking the error of the quadrotor UAV. P m and ω f represent the phase angle margin and cut-off frequency of the open-loop control system, respectively. P m _ c m d is the desired system phase angle margin, ω f _ c m d represents the desired cut-off frequency, and γ 1 , γ 2 , γ 3 , γ 4     ( γ 1 + γ 2 + γ 3 + γ 4 = 1 )   are the weight of each item in the fitness function.
According to the control system, let ω f _ c m d = 27 rad/s, P m _ c m d = 60°, and the upper and lower limits of the control parameters [ K p , K i , K d , λ , μ ] are: [1, 0, 0, 0.1, 0.1], [20, 10, 10, 1.9, 1.9]. The design fitness function parameters are shown in Table 1.
The fitness changes in the optimization process of the SSA and L-SSA are shown in Figure 7. As shown in Figure 7, the convergence speed of SSA is very fast. After ten iterations, the fitness value does not change, while for the LSSA, after 80 iterations, the value of fitness function tends to be stable. However, compared with the final optimized fitness function value, L-SSA is significantly better than SSA, so it can be concluded that L-SSA can handle the problem of SSA easily falling into the local optimum.
Eventually, the parameters of FOPID are obtained as kp = 11.9078, ki = 1.6273, kd = 1.3649, λ = 0.92925, and μ = 1.1262. Compared with the parameter tuning method based on frequency domain performance indicators [26], its Bode diagram is shown in Figure 8, the time domain response is shown in Figure 9, the system performance parameters are shown in Table 2, and the time domain indicators are shown in Table 3.
Figure 8 shows that both parameter tuning methods can ensure the phase angle margin and cutoff frequency are near the desired value. As shown in Table 2, the phase angle change rate P ω c at the cutoff frequency of the FOPID open-loop system tuned by the frequency domain optimization algorithm is almost 0, but compared with the FOPID parameters optimized based on the L-SSA, the average phase angle P a v e r near the cutoff frequency is 0. The fitness function index shows that the fitness function of the L-SSA is smaller, and its performance index is better.
It can be seen from Figure 9 and Table 3 that the fractional-order PID controller tuned based on the L-SSA has a smaller adjustment time and overshoot. Consequently, the FOPID open-loop system optimized based on the sparrow algorithm has better time-domain and frequency-domain performance.

5. Simulation and Experimental Analysis

5.1. Simulation Verification

Assume the initial states of the quadrotor are x = z = 0, h = 0 m, V = 0 m/s, θ = φ = ψ = 0°.
To observe the performance of three different controllers, the quadrotor is controlled under three conditions, including undisturbed, external disturbance, and parameter perturbation. The parameters of each controller tuned by L-SSA are shown in Table 4.
  • Case A: Nominal model
The desired pitch angle signal is given as a square wave signal with an amplitude of 10° and a period of 4 s, and the simulation curves of each controller are as in Figure 10.
Figure 10 shows that the control system based on FOLADRC has a smaller overshoot and smaller adjustment time compared with the control system based on PID and LADRC. However, the performance indicators of the three controllers are not very different.
Due to the absence of internal uncertainties and external disturbances, the disturbance values estimated by LADRC and FOLADRC are small, and the LESO does not have large disturbance compensation. Consequently, the performance of the control system based on PID, LADRC, and FOLADRC is similar.
  • Case B: External disturbance
In this condition, the pitching moment disturbance is applied at 3 s and can be equivalent to a sine disturbance, that is
w ( t ) = 0.014 sin   t
The simulation results according to Equation (23) are as follows:
It can be seen from Figure 11 that the quadrotor system under FOLADRC or LADRC has smaller steady-state errors than PID since LESO can observe the external disturbance torque and cancel it. In addition, from Figure 11, it can be known that under external disturbance, the overshoot of FOLADRC has been reduced by 11% compared with LADRC.
  • Case C: Parameter perturbation
Assuming that the actual system deviates from the nominal system, the parameter perturbation Δb0 = 5 is selected for simulation control and the simulation curves are as in Figure 12 and Figure 13:
In this case, the control system based on FOLADRC and LADRC can achieve good control performance, but the overshoot of FOLADRC controller has been reduced by 38%. Furthermore, the FOLADRC and LADRC controllers have more robustness compared with the PID controller.
Generally speaking, under the three conditions, both FOLADRC and LADRC controllers have better control performance, and the FOLADRC controller can track the desired pitch angle faster and more accurately, showing no steady-state errors. Consequently, the quadrotor UAV-based FOLADRC has good control performance and strong robustness.

5.2. Verification of Aircraft Experimental Platform

The flight experiment uses Pixhawk flight control and PX4 firmware. The control period of the flight control is 0.004 s.
To ensure the safety of the flight experiment and the conditions are the same, the quadrotor UAV is placed on the flight test bench. The test bench used for the flight experiment is shown in Figure 14:
It can be seen from Figure 14 that the quadrotor UAV is fixed on the test bench and can be rotated around the center of rotation. Since the quadrotor UAV rotates around the center of rotation of the test bench during the test, disturbance related to the altitude angle ( f t a i = m g d sin   α ) will be generated during the altitude control process due to the influence of the center of gravity. d is the distance from the center of gravity of the quadrotor UAV to the center of rotation of the test bench.
The parameters of controllers are set according to Table 4.
  • Case A: Nominal model
The desired pitch angle signal is given as a step signal with an amplitude of 10°, and the experiment results of each controller are as follows:
It can be seen from Figure 15 that in the initial stage of the quadrotor UAV maneuvering, the system control performance based on LADRC, PID, and LADRC methods overlap, but when the pitch angle is further increased, the tracking performance of the three control methods appears to be relatively different. In the initial stage of command tracking, the disturbance f t a i received by the quadrotor UAV is small due to the small sinα. As the pitch angle continues to increase, the disturbance becomes larger. Because the PID cannot observe the disturbance, it only relies on a small integral term to offset the influence of the disturbance, so there is a large error. The LESO in LADRC can observe the disturbance f t a i in time and cancel it, so the overshoot is small, and the error is quickly eliminated. Compared with LADRC, FOLADRC can also cancel the disturbance error in time, which can perform instruction tracking faster and has stronger robustness.
  • Case B: External perturbation
To further verify the control performance and anti-disturbance characteristics of FOLADRC, an artificial external disturbance was applied to the quadrotor UAV to observe the control performance of the FOLADRC control law. The disturbance experiment is shown in Figure 16.
As shown in Figure 16, a water bottle with a weight of 300 g was placed on one axis of the test bench with an elastic rope when it was in a horizontally stable state. Due to the elasticity of the elastic rope, the water bottle will vibrate up and down and rotate around the axis of the elastic rope, which further causes the water in the bottle to flow, thereby introducing a large nonlinear disturbance into the quadrotor UAV to test the control performance of FOLADRC. The experiment results are shown in Figure 17:
It can be seen from Figure 17 that the control system based on PID produces a large amount of overshoot due to the disturbance of the external load, while the control system based on LADRC and FOLADRC has an increase in overshoot compared to the no-load disturbance. Compared with LADRC, the overshoot of FOLADRC has been reduced by 37% and the adjustment time has been reduced by 60%.
In addition, the anti-disturbance performance of the FOLADRC controller is tested by applying a sudden load disturbance. The experimental results are shown in Figure 18.
It can be seen from Figure 18a that at the moment of applying the disturbance, the control systems based on PID, LADRC, and FOLADRC all deviate from the equilibrium point to a certain extent, among which the PID controller deviates the most from the equilibrium point, while the FOLADRC method has the smallest error.
Figure 18b,c show the changes in the disturbance observations of the control system based on LADRC and FOLADRC after loading the disturbance. When the disturbance load is applied to the quadrotor UAV, the disturbance observations of both controllers increase rapidly, thereby canceling the overshoot caused by the disturbance. The disturbance value has a certain vibration at the beginning. Compared with the LADRC method, the disturbance observation value of the FOLADRC method has smaller vibration and can be stabilized to a certain value in a short time, so that the entire system stabilizes quickly.
Consequently, based on the above simulation experiments and flight tests, it can be seen that the quadrotor UAV control method based on FOLADRC has better control performance and strong robustness in the case of the nominal disturbance-free condition, external load disturbance, and parameter uncertainty.

6. Conclusions

Since the ADRC system is sensitive to the system controller gain, the fractional calculus operator was introduced in this paper to improve the robustness and control performance of the system. This paper studied the design method of the FOLADRC controller for a quadrotor UAV and solved the problem that the attitude control accuracy and stability of the quadrotor UAV are greatly reduced under large parameter uncertainty and strong external disturbance. To verify the feasibility of the designed control system implemented in the quadrotor UAV system, this paper applies the designed FOLADRC controller to the UAV platform in nominal experiments and anti-disturbance experiments. The experimental results show that the designed control system has the advantages of smaller adjusting time, less overshoot, and stronger robustness, which demonstrates that FOLADRC has greater advantages in practical engineering applications.
In addition, a parameter tuning method of the FOLADRC controller based on the L-SSA was proposed, which not only solves the problem that basic SSA easily falls into local optimization, but also settles the issue that the FOLADRC parameter tuning method based on the frequency domain index cannot take into account the rapidity, accuracy, and reliability of the time domain at the same time.

Author Contributions

Conceptualization, J.S. (Jia Song) and Y.H.; methodology, J.S. (Jia Song) and J.S. (Jiangcheng Su); software, J.S. (Jiangcheng Su) and Y.H.; validation, J.S., M.Z. and S.A.; formal analysis, J.S. (Jia Song) and J.S. (Jiangcheng Su); investigation, J.S. (Jia Song) and J.S. (Jiangcheng Su); resources, J.S. (Jia Song), M.Z. and S.A.; data curation, Y.H. and J.S. (Jia Song); writing—original draft preparation, J.S. (Jia Song); writing—review and editing, J.S. (Jia Song) and Y.H.; visualization, J.S. (Jia Song); supervision, J.S. (Jia Song); project administration, J.S. (Jia Song); funding acquisition, J.S. (Jia Song). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Key Laboratory of Hunan Province for Control Technology of Distributed Electric Propulsion Air Vehicle OF FUNDER grant number 62073020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank their colleagues for their constructive suggestions and research assistance throughout this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kim, J.; Gadsden, S.A.; Wilkerson, S.A. A Comprehensive Survey of Control Strategies for Autonomous Quadrotors. Can. J. Electr. Comput. Eng. 2020, 43, 3–16. [Google Scholar] [CrossRef]
  2. Zhang, W.; Gao, F.; Jiang, N.; Zhang, C.; Zhang, Y. High-Temporal-Resolution Forest Growth Monitoring Based on Segmented 3D Canopy Surface from UAV Aerial Photogrammetry. Drones 2022, 6, 158. [Google Scholar] [CrossRef]
  3. Qian, L.; Liu, H. Path Following Control of a Quadrotor UAV with A Cable Suspended Payload Under Wind Disturbances. IEEE Trans. Ind. Electron. 2020, 67, 2019–2021. [Google Scholar] [CrossRef]
  4. Praveen, V.; Pillai, A.S. Modeling and simulation of a quadcopter using PID controller. Int. J. Circuit Theory Appl. 2016, 9, 7151–7158. [Google Scholar]
  5. Noordin, A.; Basri, M.; Mohamed, Z. Sensor Fusion for Attitude Estimation and PID Control of Quadrotor UAV. Int. J. Electr. Electron. Eng. Telecommun. 2018, 7, 183–189. [Google Scholar] [CrossRef]
  6. Zhang, J.; Zhou, L.; Li, C.; Wen, B. Binary observers based control for quadrotor unmanned aerial vehicle with disturbances and measurement delay. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2017, 232, 2807–2820. [Google Scholar] [CrossRef]
  7. Tang, W.; Wang, L.; Gu, J.; Gu, Y. Single Neural Adaptive PID Control for Small UAV Micro-Turbojet Engine. Sensors 2020, 20, 345. [Google Scholar] [CrossRef] [PubMed]
  8. Nguyen, N.P.; Park, D.; Ngoc, D.N.; Xuan-Mung, N.; Huynh, T.T.; Nguyen, T.N.; Hong, S.K. Quadrotor Formation Control via Terminal Sliding Mode Approach: Theory and Experiment Results. Drones 2022, 6, 172. [Google Scholar] [CrossRef]
  9. 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]
  10. Dalwadi, N.; Deb, D.; Ozana, S. Rotor Failure Compensation in a Biplane Quadrotor Based on Virtual Deflection. Drones 2022, 6, 176. [Google Scholar] [CrossRef]
  11. Jiao, R.; Chou, W.; Rong, Y.; Dong, M. Anti-disturbance control for quadrotor UAV manipulator attitude system based on fuzzy adaptive saturation super-twisting sliding mode observer. Appl. Sci. 2020, 10, 3719. [Google Scholar] [CrossRef]
  12. Han, J. From PID to Active Disturbance Rejection Control. IEEE Trans. Ind. Electron. 2009, 56, 900–906. [Google Scholar] [CrossRef]
  13. Castañeda, L.A.; Luviano-Juárez, A.; Ochoa-Ortega, G.; Chairez, I. Tracking control of uncertain time delay systems: An ADRC approach. Control Eng. Pract. 2018, 78, 97–104. [Google Scholar] [CrossRef]
  14. Wu, Z.; Gao, Z.; Li, D.; Chen, Y.; Liu, Y. On transitioning from PID to ADRC in thermal power plants. Control Theory Technol. 2021, 19, 3–18. [Google Scholar] [CrossRef]
  15. Yoo, D.; Yau, S.T.; Gao, Z.Q. Optimal fast tracking observer bandwidth of the linear extended state observer. Int. J. Control 2007, 80, 102–111. [Google Scholar] [CrossRef]
  16. Hu, D.; Pei, Z.; Tang, Z. Single-parameter-tuned attitude control for quadrotor with unknown disturbance. Appl. Sci. 2020, 10, 5564. [Google Scholar] [CrossRef]
  17. Najm, A.A.; Ibraheem, I.K. Altitude and attitude stabilization of UAV quadrotor system using improved active disturbance rejection control. Arab. J. Sci. Eng. 2020, 45, 1985–1999. [Google Scholar] [CrossRef]
  18. Wang, Z.; 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]
  19. Podlubny, I. Fractional-order systems and PIλDμ-controllers. IEEE Trans. Autom. Control 1999, 44, 208–214. [Google Scholar] [CrossRef]
  20. Seyedtabaii, S. New flat phase margin fractional order PID design: Perturbed UAV roll control study. Robot. Auton. Syst. 2017, 96, 58–64. [Google Scholar] [CrossRef]
  21. Erenturk, K. Fractional-Order PI^λD^μ and Active Disturbance Rejection Control of Nonlinear Two-Mass Drive System. IEEE Trans. Ind. Electron. 2013, 60, 3806–3813. [Google Scholar] [CrossRef]
  22. Liu, N.; Cao, S.; Fei, J. Fractional-Order PID Controller for Active Power Filter Using Active Disturbance Rejection Control. Math. Probl. Eng. 2019, 24, 6907570. [Google Scholar] [CrossRef]
  23. Li, Z.; Wei, Y.; Zhou, X.; Wang, J.; Wang, J.; Wang, Y. Differential flatness-based ADRC scheme for underactuated fractional-order systems. Int. J. Robust Nonlinear Control 2020, 30, 2832–2849. [Google Scholar] [CrossRef]
  24. Al-Saggaf, U.M.; Mehedi, I.M.; Mansouri, R.; Bettayeb, M. Fractional order linear ADRC-based controller design for heat-flow experiment. Math. Probl. Eng. 2021, 2021, 7291420. [Google Scholar] [CrossRef]
  25. Chen, Y.; Zhang, T. Quadrotor aircraft control based on fractional order PID. J. Tianjing Polytech. Univ. 2019, 38, 58–63. [Google Scholar]
  26. Gao, K.; Song, J.; Wang, X.; Li, H. Fractional-order proportional-integral-derivative linear active disturbance rejection control design and parameter optimization for hypersonic vehicles with actuator faults. Tsinghua Sci. Technol. 2020, 26, 9–23. [Google Scholar] [CrossRef]
  27. Ge, J.; Guo, X.; Liu, Y. Attitude Control of Quad-Rotor UAV Based on GA-ADRC. Dyn. Syst. Control 2020, 9, 109–120. [Google Scholar] [CrossRef]
  28. Hangxuan, H.; Haibin, D. A multi-strategy pigeon-inspired optimization approach to active disturbance rejection control parameters tuning for vertical take-off and landing fixed-wing UAV. Chin. J. Aeronaut. 2022, 35, 19–30. [Google Scholar]
  29. Xue, J.; Shen, B. A novel swarm intelligence optimization approach: Sparrow search algorithm. Syst. Sci. Control Eng. 2020, 8, 22–34. [Google Scholar] [CrossRef]
  30. Xing, Z.; Yi, C.; Lin, J.; Zhou, Q. Multi-component fault diagnosis of wheelset-bearing using shift-invariant impulsive dictionary matching pursuit and sparrow search algorithm. Measurement 2021, 178, 109375. [Google Scholar] [CrossRef]
  31. Li, X.; Ma, X.; Xiao, F.; Xiao, C.; Wang, F.; Zhang, S. Time-series production forecasting method based on the integration of Bidirectional Gated Recurrent Unit (Bi-GRU) network and Sparrow Search Algorithm (SSA). J. Pet. Sci. Eng. 2022, 208, 109309. [Google Scholar] [CrossRef]
  32. Wang, Y.; Ding, S.; Wang, L.; Du, S. A manifold p-spectral clustering with sparrow search algorithm. Soft Comput. 2022, 2022, 1765–1777. [Google Scholar] [CrossRef]
Figure 1. The quadrotor UAV under consideration.
Figure 1. The quadrotor UAV under consideration.
Drones 06 00229 g001
Figure 2. Body coordinate system and NED coordinate system.
Figure 2. Body coordinate system and NED coordinate system.
Drones 06 00229 g002
Figure 3. FOLADRC-based system structure diagram.
Figure 3. FOLADRC-based system structure diagram.
Drones 06 00229 g003
Figure 4. The structure diagram of FOPID.
Figure 4. The structure diagram of FOPID.
Drones 06 00229 g004
Figure 5. The relationship between the three types of sparrows.
Figure 5. The relationship between the three types of sparrows.
Drones 06 00229 g005
Figure 6. Flowchart of L-SSA.
Figure 6. Flowchart of L-SSA.
Drones 06 00229 g006
Figure 7. The fitness value change graph of the SSA and L-SSA.
Figure 7. The fitness value change graph of the SSA and L-SSA.
Drones 06 00229 g007
Figure 8. Bode diagram of FOPID control system with different parameter tuning methods.
Figure 8. Bode diagram of FOPID control system with different parameter tuning methods.
Drones 06 00229 g008
Figure 9. Time domain response of FOPID control with different parameter tuning methods.
Figure 9. Time domain response of FOPID control with different parameter tuning methods.
Drones 06 00229 g009
Figure 10. Curves of pitch angle under the nominal model.
Figure 10. Curves of pitch angle under the nominal model.
Drones 06 00229 g010
Figure 11. Curves of pitch angle under pitching moment disturbance.
Figure 11. Curves of pitch angle under pitching moment disturbance.
Drones 06 00229 g011
Figure 12. Curves of pitch angle under parameter perturbation.
Figure 12. Curves of pitch angle under parameter perturbation.
Drones 06 00229 g012
Figure 13. Curves of pitch angle under parameter perturbation (partial enlargement).
Figure 13. Curves of pitch angle under parameter perturbation (partial enlargement).
Drones 06 00229 g013
Figure 14. Flying test bench.
Figure 14. Flying test bench.
Drones 06 00229 g014
Figure 15. Quadrotor UAV altitude tracking test pitch angle tracking diagram.
Figure 15. Quadrotor UAV altitude tracking test pitch angle tracking diagram.
Drones 06 00229 g015
Figure 16. Anti-disturbance experiment diagram for Quadrotor UAV.
Figure 16. Anti-disturbance experiment diagram for Quadrotor UAV.
Drones 06 00229 g016
Figure 17. Pitch angle tracking diagram of anti-disturbance flight test for quadrotor UAV.
Figure 17. Pitch angle tracking diagram of anti-disturbance flight test for quadrotor UAV.
Drones 06 00229 g017
Figure 18. Anti-disturbance flight test for quadrotor UAV. (a) Error diagram of an anti-disturbance experiment for quadrotor UAV. (b) z3 observation diagram in LADRC. (c) z3 observation diagram in FOLADRC.
Figure 18. Anti-disturbance flight test for quadrotor UAV. (a) Error diagram of an anti-disturbance experiment for quadrotor UAV. (b) z3 observation diagram in LADRC. (c) z3 observation diagram in FOLADRC.
Drones 06 00229 g018aDrones 06 00229 g018b
Table 1. Fitness function parameters table.
Table 1. Fitness function parameters table.
ParametersValue
γ 1 0.15
γ 2 0.5
γ 3 0.25
γ 4 0.1
ω f _ c m d 27
P m _ c m d 60
Table 2. Frequency domain index and fitness function table of different methods.
Table 2. Frequency domain index and fitness function table of different methods.
P m ω c P ω c   P a v e r   fitf
Frequency domain analysis60.0°26.997 rad/s0.07°/(rad/s)1.30°/(rad/s)−39.81
L-SSA60.9°27.124 rad/s0.51°/(rad/s)0.83°/(rad/s)−50.92
Table 3. System performance in the time-domain of different methods.
Table 3. System performance in the time-domain of different methods.
Performance ParametersFrequency Domain AnalysisL-SSA
Adjustment time (2% error, s)0.5250.503
Overshoot(%)16.910.2
Table 4. Controller parameters.
Table 4. Controller parameters.
kpkikdλμ
FOLADRC11.90781.62731.36490.929251.1262
LADRC11.52681.32991.4928--
PID10.69390.30001.2652--
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Song, J.; Hu, Y.; Su, J.; Zhao, M.; Ai, S. Fractional-Order Linear Active Disturbance Rejection Control Design and Optimization Based Improved Sparrow Search Algorithm for Quadrotor UAV with System Uncertainties and External Disturbance. Drones 2022, 6, 229. https://doi.org/10.3390/drones6090229

AMA Style

Song J, Hu Y, Su J, Zhao M, Ai S. Fractional-Order Linear Active Disturbance Rejection Control Design and Optimization Based Improved Sparrow Search Algorithm for Quadrotor UAV with System Uncertainties and External Disturbance. Drones. 2022; 6(9):229. https://doi.org/10.3390/drones6090229

Chicago/Turabian Style

Song, Jia, Yunlong Hu, Jiangcheng Su, Mingfei Zhao, and Shaojie Ai. 2022. "Fractional-Order Linear Active Disturbance Rejection Control Design and Optimization Based Improved Sparrow Search Algorithm for Quadrotor UAV with System Uncertainties and External Disturbance" Drones 6, no. 9: 229. https://doi.org/10.3390/drones6090229

APA Style

Song, J., Hu, Y., Su, J., Zhao, M., & Ai, S. (2022). Fractional-Order Linear Active Disturbance Rejection Control Design and Optimization Based Improved Sparrow Search Algorithm for Quadrotor UAV with System Uncertainties and External Disturbance. Drones, 6(9), 229. https://doi.org/10.3390/drones6090229

Article Metrics

Back to TopTop