Auto-Landing of Moving-Mass Actuated Unmanned Aerial Vehicles Based on Linear Active Disturbance Rejection Control

: Unlike the roll motion of the unmanned aerial vehicle (UAV) controlled by the ailerons, the moving-mass actuated unmanned aerial vehicle (MAUAV) uses the motion of the mass block inside the wing to generate the roll moment. The light weight and severe coupling of lateral and longitudinal motion of this type of small UAV make its landing control a challenging task. Considering the above problems, the dynamic models of MAUAV are ﬁrst established. Then, forward velocity, altitude, attitude, and moving-mass position controllers are designed separately to make the MAUAV track a given path during the landing process. Linear active disturbance rejection control (LADRC) is introduced in the design process of all four controllers, compensating for unknown disturbances in the system. Simulation results show that the proposed control scheme can achieve fast and accurate tracking of forward velocity and ﬂight trajectory commands with good robustness to model uncertainties.


Introduction
In recent years, thanks to the advantages of low cost, simple structure, strong stealth, and high mobility, unmanned aerial vehicles (UAVs) have been widely used in civil and military fields, such as meteorological monitoring, aerial surveying and mapping, unmanned cargo system, and battlefield surveillance [1][2][3][4][5].The application of UAVs in these fields can reduce the consumption of human resources, improve work efficiency, and at the same time, ensure the safety of personnel to reduce the casualty rate.According to the classification of wing structure, UAVs are mainly divided into two types: fixed-wing and rotary-wing.They each have their advantages and disadvantages, but in general, fixed-wing UAVs are more widely used than rotary-wing UAVs and capable of long-distance flight missions.Traditionally, fixed-wing UAVs realize their attitude movement through the control of ailerons, rudders, and elevators.With the development of moving-mass control techniques in many different fields (including satellites [6,7], re-entry vehicles [8][9][10], and underwater vehicles [11], etc.), researchers have also started to focus on applying this control scheme to the flight control of UAVs.As a reliable alternative to the traditional scheme, the movingmass control scheme allows the controlled mass blocks to be mounted inside the aircraft, and the attitude motion of UAVs can be controlled by driving the movement of the mass blocks.Therefore, this control scheme can optimize the aerodynamic shape, simplify the wing structure to reduce the manufacturing difficulty of the wing shape, and improve the integrity of the shape to further enhance the stealthiness of the UAVs.
The landing phase is one of the most critical and accident-prone phases of the entire mission execution process for fixed-wing UAVs.A typical UAV landing consists of three phases: approach, glide slope, and flare.In order to reduce accidents during the landing process of UAVs, it is necessary to design an auto-landing control strategy with good control performance and high precision in order that UAVs can closely track the desired path in the three landing stages and finally land safely.The application of superior autolanding control strategies will also greatly improve the recovery rate of UAVs, and thus save costs.The more common methods used for UAV control include PID control, H ∞ , sliding mode control, model predictive control, backstepping, adaptive control, linearquadratic Gaussian approach, robust adaptive control, etc. [12][13][14][15][16][17][18].They are applied in all phases of UAV flight (including the landing phase, which has the highest accident rate) to improve flight stability and safety.In particular, an attitude controller based on the PID method was designed to overcome the loss of stability in the flying aircraft [12].Rao and Go [13] proposed an automatic landing controller for UAVs using sliding mode technology, and avoided the problems of chatter and high feedback gain by selecting suitable reaching laws.Mihai [19] designed backstepping and dynamic inversion-based control strategies to assist the automatic landing of UAVs under the superposition of various wind disturbances, which are measured by the estimator and compensated in the design process of the controller.With the development of intelligent systems, more and more intelligent control algorithms are being applied to the flight control of UAVs.A convolutional neural network is trained to predict the heading angle of the aircraft by using image inputs in the case of beacon failure [20].Authors in [21] proposed a deep learning-based vision transformer particle region-based convolutional neural network for navigation when the fixed-wing UAV was autonomous landing.
It is not difficult to find that the above control methods are all applied to fixed-wing UAVs driven by ailerons.For moving-mass actuated fixed-wing UAVs, many related studies [22][23][24][25] have been carried out by Erturk et al.Based on dynamics analysis and simulation tests of MAUAVs under different flight conditions, they pointed out that the moving-mass scheme can reduce fuel consumption to some extent since this method is able to eliminate the drag generated by aerodynamic surface deflection.In [16], authors concluded that the controllability of moving-mass control is similar or higher than that with aerodynamic surface solutions when the UAV is flying at low speeds.As mentioned earlier, the development of artificial intelligence in recent years has led to its increasingly widespread use in fixed-wing UAVs.At the same time, researchers have also used AI technology to control MAUAVs.Qiu et al. [26] proposed a fixed-wing UAV structure with a combination of rudder, elevator, and movable slider to control the attitude.The offset of the internal slider is able to generate a rolling moment to control the roll motion of the UAV, thus the parameters of the slider are important factors affecting the dynamic performance of the vehicle.The maximum rolling moment provided by this solution depends to a large extent on the mass and displacement of the slider.In addition, this moment is related to the attitude of the UAV, which decreases with the increasing roll (and pitch) angle within a certain range when the mass and displacement of the slider remain unchanged.Based on this structure, they analyzed its dynamic characteristics and designed an adaptive sliding mode attitude controller incorporating the fuzzy system and RBF neural network for it.Based on deep reinforcement learning, an end-to-end attitude controller of MAUAVs was designed and the robustness of this controller was verified by random initialization and parameter perturbation during simulation [27].However, these studies related to MAUAVs did not design detailed and complete automatic landing control strategies for the landing phase.Both the coupling of different channels and model uncertainties in MAUAV systems bring difficulties to the design of auto-landing control strategies.
The linear active disturbances rejection controller (LADRC) was first proposed in the literature [28].This control method works by expanding the unknown disturbances into additional states of the system, and then estimating the values of all states using a linear state estimator.Finally, the control input of the system is designed based on PI control and the estimated values.Compared with the above-mentioned control methods, this control method not only has a simple design form, the selection of control parameters follows certain rules (reducing the possibility of artificial blind selection of parameter values), but also can obtain excellent tracking performance.The simplicity of the structure, the regularity of the parameters, and the accuracy of the control results are more conducive to the practical application in the system with unknown disturbance.A practical application of LADRC to the sensorless control of internal permanent-magnet synchronous motor was carried out by authors in [29].They used linear extended state observers (LESO) to achieve high-performance position estimation and good robustness of the system over full-speed range.The LADRC has also been extensively studied for use in the field of aerial vehicles [30][31][32].In [31], a fuzzy controller and LADRC are combined to simultaneously obtain their adaptive and anti-interference capabilities, and subsequently increase the response speed of a quadrotor UAV and improve the dynamic performance of the system.Similarly, Wang et al. [32] designed a mass adaptive control method combined with sliding mode control (SMC) and LADRC for a quadrotor load UAV.This is due to the fact that the bandwidth limitation of LADRC can be overcome by enhancing the robustness of the system through SMC.And the LESO is able to estimate the unmodelled dynamics caused by SMC chattering and external unknown disturbances.In this paper, in the view of achievements of the LADRC method in the field of control, four controllers are designed with the help of this method and work together to realize the autonomous landing of MAUAV.The main contributions of this paper are summarized in the following four aspects: 1.
The characteristics and dynamic models of the moving-mass actuated unmanned aerial vehicle are presented and an automatic landing control strategy is designed, which has not been studied yet for this type of UAV.Both lateral and longitudinal movements of the UAV are taken into account during landing and the tracking error in both planes is kept within a small range.This is in contrast to many papers [33,34] where only separate lateral and longitudinal dynamics are studied in auto-landings.

2.
The control strategy of the MAUAV is divided into an altitude-attitude slider system and a forward velocity tracking system.The former includes three controllers of altitude, attitude, and slider to provide the expected driving force for slider and expected deflections for rudder and elevator; the latter establishes a throttle command to maintain a steady forward velocity, due to the fact that accurate control of forward velocity is essential for an efficient and safe landing.In addition, compared with the attitude control of the same type of UAV [26], this paper considers the influence of wind disturbance during the whole landing process, which is more in line with the actual situation.

3.
The design of all four controllers in this paper is based on the LADRC, which simplifies the difficulty of the controller design and the complexity of the algorithm.The designed controllers are not based on accurate dynamic models and are able to estimate and compensate for external disturbances, offering better applicability than existing controllers used in the field of mass-actuated UAVs [24,35].4.
The LuGre friction model has been introduced into the moving-mass position model to provide a more detailed description of the friction behavior of the slider and to make the simulation environment more realistic.
The remainder of this paper is arranged as follows: In Section 2, by analyzing the motion characteristics of the MAUAV, its dynamic models, such as attitude, position, and slider displacement models, are established.Based on these dynamic models, the controller design models of forward velocity, altitude, attitude, and moving-mass position are sorted out in Section 3, and four corresponding controllers with disturbance compensation functions are designed.In addition, by drawing the actual flight trajectory of the MAUAV and comparing it with the expected trajectory in the simulation test, the specific landing process of the MAUAV is demonstrated, and the effectiveness of the auto-landing method composed of four controllers is verified in Section 4. Finally, Section 5 includes some conclusions of the paper.

Problem Formulation and Dynamic Models 2.1. Configuration of the MAUAV
The schematic diagram of the studied moving-mass actuated UAV is shown in Figure 1.Unlike conventional fixed-wing UAVs fitted with ailerons, the roll motion of the MAUAV is achieved by mounting a movable slider in the lateral direction.The longitudinal movement and forward power mode remain unchanged and are still controlled by the elevator and propeller, respectively.Compared to the dual-slider scheme in literature [24], the rudder was retained in the single-slide scheme.This kind of scheme not only reduces the manufacturing difficulty of the airframe, but also alleviates the space constraint problem, which in turn improves the control capability of the vehicle by installing a larger slider [26].Considering the special characteristics of the MAUAV, the definition of the body frame is given first.The origin O is the center of mass when the airframe is not fitted with the slider.The x b -axis points out the nose along the airframe's longitudinal axis; the y b -axis points out the right wing; the z b -axis points out the belly and completes a right-handed rectangular coordinate.
Drones 2023, 7, x FOR PEER REVIEW 4 of 24 method composed of four controllers is verified in Section 4. Finally, Section 5 includes some conclusions of the paper.

Configuration of the MAUAV
The schematic diagram of the studied moving-mass actuated UAV is shown in Figure 1.Unlike conventional fixed-wing UAVs fitted with ailerons, the roll motion of the MAUAV is achieved by mounting a movable slider in the lateral direction.The longitudinal movement and forward power mode remain unchanged and are still controlled by the elevator and propeller, respectively.Compared to the dual-slider scheme in literature [24], the rudder was retained in the single-slide scheme.This kind of scheme not only reduces the manufacturing difficulty of the airframe, but also alleviates the space constraint problem, which in turn improves the control capability of the vehicle by installing a larger slider [26].Considering the special characteristics of the MAUAV, the definition of the body frame is given first.The origin O is the center of mass when the airframe is not fitted with the slider.The x b -axis points out the nose along the airframe's longitudinal axis; the y b -axis points out the right wing; the z b -axis points out the belly and completes a righthanded rectangular coordinate.

Wind Description
The wind acts as a disturbance by altering the attitude and velocity of the UAV.Define g V as the velocity of the airframe relative to the ground, a V as the velocity of the airframe relative to the surrounding air mass, and w V as the velocity of the air mass rel- ative to the ground.These velocities can be expressed in the body frame in terms of components along the x b , y b , and z b axes as . The relationship between ground velocity, air velocity, and wind velocity is given by = + .
For the purposes of simulation, we will assume that the total wind vector consists of the steady ambient wind s w V and random winds g w
The steady wind can be represented in the

Wind Description
The wind acts as a disturbance by altering the attitude and velocity of the UAV.Define V g as the velocity of the airframe relative to the ground, V a as the velocity of the airframe relative to the surrounding air mass, and V w as the velocity of the air mass relative to the ground.These velocities can be expressed in the body frame in terms of components along the x b , y b , and z b axes as The relationship between ground velocity, air velocity, and wind velocity is given by For the purposes of simulation, we will assume that the total wind vector consists of the steady ambient wind V w s and random winds V w g (including wind gusts and atmospheric disturbances), i.e., V w = V w s + V w g .The steady wind can be represented in the inertial frame as V i w s = [w n s , w e s , w d s ] T , where w n s , w e s , w d s are the speeds of the steady wind in the north, east, and down directions, respectively.The stochastic component of the wind is typically expressed in the aircraft body frame as V b w g = u w g , v w g , w w g T .The two most widely used models of continuous gusts are the Dryden and von Karman models [36].The Dryden turbulence model is mathematically simpler and the forming filters in this model Drones 2023, 7, 539 5 of 23 are derived from the spectral square roots of the spectrum equations.The definitions of the longitudinal, lateral, and vertical component spectra functions are described as (s+ Va Lv ) 2 ; where σ u , σ v , and σ w are the intensities of the turbulence; L u , L v , and L w are spatial wavelengths; r is the airspeed of the vehicle.Then, the total winds in the body frame can be expressed as where R b i is the rotation matrix from inertial frame to body frame given as where ϕ, θ, ψ are roll, pitch, and yaw angles of the MAUAV, and their units are radians in the following dynamic models.
The derivative of the wind is T , which is calculated as [19] .(5)

Dynamic Model of the MAUAV
The translational dynamic equations describing the components of the velocity relative to the ground are shown in body frame as follows [26]: . .
where p, q, r are the angular rates of roll, pitch, and yaw angles, respectively; q = ρV 2 a /2 is the UAV dynamic pressure depending on the nominal flight velocity (V a = u 2 r + v 2 r + w 2 r ) and ρ is the atmospheric density; S is the wing area of the MAUAV; m S is the total mass of the whole system, which is the sum of the slider mass m 1 and airframe mass m B , i.e., m S = m 1 + m B ; µ 1 = m 1 /m S is the mass ratio of the slider; g is the gravitational acceleration; T is slider's position vector in the body frame (since the slide is fixedly mounted on the airframe and is placed along the y b -axis of the body frame, it is possible to . F P is the thrust force, which can be expressed as [26] where S prop is the area swept by the propeller; C prop is the propeller's efficiency factor; k motor is a constant depending on the properties of the selected engine; δ t is the throttle control size of the engine, which varies from 0 to 1; the expressions of the dimensionless coefficients C X , C Y , C Z are described as [1] Drones 2023, 7, 539 6 of 23 where δ r is the rudder deflection, which is measured in radian in this model; b is the wing span; are coefficients of the lateral force; α and β represent the attack and the sideslip angles, their expressions are [1] the lift and the drag coefficients C L and C D are described as [1] where c is the average aerodynamic chord length of the wing; δ e is the elevator deflection, and its unit is radian in the dynamic models; e is the Oswald efficiency factor, which ranges between 0.8 and 1.0; AR is the wing aspect ratio, which is equal to b 2 /S; are coefficients of the lift force and drag force.
Define ω = [p, q, r] T as the angular rates vector of MAUAV, then the roll, pitch, and yaw angular rates relative to the body frame can be described in the rotational dynamic equations as follows [26]: where M is the external moment in the body frame, which is described as [26] where Ω = k Ω δ t is the propeller speed (rotation angular velocity) and k T p is a constant coefficient depending on the selected propeller; C l , C m , C n are the aerodynamic coefficients with the expressions [26] where denote related aerodynamic coefficients.And J is the inertia moment, which is described as [26] where where J x B , J y B , J z B , J xy B , J yz B , J xz B represent the inertia moments of the airframe when it is not fitted with the slider.
Drones 2023, 7, 539 7 of 23 Define Φ = [ϕ, θ, ψ] T as an angular vector, then the rotational kinematic equations of the MAUAV can be obtained as [26] .Φ = Aω, (16) where A is a transformation matrix described as Newton's second law applied to the moving mass can be stated as where is the time derivative in the inertial frame, V 1 is the velocity vector of the moving mass with respect to the inertial frame, and F 1 is the sum of all external forces acting on the moving mass.
The derivative of velocity taken in the inertial frame can be written in terms of the derivative in the body frame and the angular velocity as where is the time derivative in the body frame and ω b/i is the angular velocity of the MAUAV with respect to the inertial frame.
Combining ( 18) and ( 19), we can easily apply Newton's second law by expressing the forces and velocities in the body frame as where Expanding Equation (20) to take the component representation along the y b -direction of the body frame, then the dynamic motion of the moving mass in the body frame can be expressed as where F b 1y = F u + m 1 g sin ϕ cos θ − F f , F u represents the active control force exerted by the servo mechanism on the movable mass and F f is the friction that hinders the motion of the mass block.To precisely model the nonlinear friction, the LuGre model [37] is extensively utilized, which can be described by where σ 0 , σ 1 , σ 2 are friction force parameters, which can be physically interpreted as the stiffness of the bristles between two contact surfaces, damping coefficient of the bristles, and viscous coefficient, respectively; the unmeasurable internal friction state z physically stands for the average deflection of the bristles between two contact surfaces, and its characteristic is governed by We assume that the inertial coordinates are defined as being aligned with the runway.The origin of the inertial system is the desired touchdown point of MAUAV and X-axis of inertial coordinates coincides with the runway centerline, which is directed north.Then, the translational kinematic equations can be obtained as [19,38] .
where X is the horizontal displacement; Y is the lateral deviation relative to the runway; and H is the flight altitude.

Design of the Controllers
This section describes the specific design process of the four controllers in the automatic landing system (Figure 2).These controllers enable the automatic landing of the MAUAV by providing the desired values of the four controls: the driving force of the slider-F u , elevator deflection-δ e , rudder deflection-δ r , and throttle command-δ t .According to the functional classification, the control strategy of this UAV is divided into two main parts: an altitude-attitude slider system and a forward velocity tracking system.The former consists of three controllers: an altitude controller, an attitude controller, and a movingmass position controller.The altitude controller generates the desired pitch angle signal based on the desired altitude trajectory.The attitude controller generates the desired slider displacement, elevator deflection, and rudder deflection based on the desired pitch, roll, and yaw signals.The moving-mass position controller generates the desired driving force of the slider based on the desired slider displacement.The latter develops a throttle input to ensure that the MAUAV tracks desired forward velocity throughout the landing process.All controllers are designed using LADRC, with the help of which unknown disturbances are compensated and the form of controllers is not complicated.The controller design results are slightly different depending on the system order.Figure 3 illustrates the process of LADRC-based controller design for first-order and second-order systems.We assume that the inertial coordinates are defined as being aligned with the runway.The origin of the inertial system is the desired touchdown point of MAUAV and X-axis of inertial coordinates coincides with the runway centerline, which is directed north.Then, the translational kinematic equations can be obtained as [19,38] ) where X is the horizontal displacement; Y is the lateral deviation relative to the runway; and H is the flight altitude.

Design of the Controllers
This section describes the specific design process of the four controllers in the automatic landing system (Figure 2).These controllers enable the automatic landing of the МАUAV by providing the desired values of the four controls: the driving force of the slider-F u , elevator deflection-δ e , rudder deflection-δ r , and throttle command-δ t .According to the functional classification, the control strategy of this UAV is divided into two main parts: an altitude-attitude slider system and a forward velocity tracking system.The former consists of three controllers: an altitude controller, an attitude controller, and a moving-mass position controller.The altitude controller generates the desired pitch angle signal based on the desired altitude trajectory.The attitude controller generates the desired slider displacement, elevator deflection, and rudder deflection based on the desired pitch, roll, and yaw signals.The moving-mass position controller generates the desired driving force of the slider based on the desired slider displacement.The latter develops a throttle input to ensure that the MAUAV tracks desired forward velocity throughout the landing process.All controllers are designed using LADRC, with the help of which unknown disturbances are compensated and the form of controllers is not complicated.The controller design results are slightly different depending on the system order.Figure 3 illustrates the process of LADRC-based controller design for first-order and second-order systems.

Forward Velocity Controller
Equation (25) shows that the position of the MAUAV relative to the ground runway during the approach, glide slope, and flare is affected by the speed of vehicle.Considering that the ground speed of the MAUAV during landing with low wind disturbances mainly depends on the forward velocity u, it is necessary to control the forward velocity to make the MAUAV land efficiently and safely.Based on Equation ( 6), we can derive the forward velocity dynamic model as ( ) From Equations ( 7), ( 8), ( 10) and ( 26), the design model can be written as where ( ) According to the design model (27), the expansion state space can be expressed as where u f is designed as an extended state.
To estimate the two states of the new system, we construct the LESO as where û and ˆu f represent the estimate of u and u f ; 1 u β and 2 u β are the tuning parame- ters and are determined by the value of the bandwidth of the LESO, which is denoted as ω .Perfect estimations of the true value of the states in the system can be achieved when adjusting these two gains by the following relationships [39]: ; , Based on the accurate estimation of the system state, the final throttle command has the following form:

Forward Velocity Controller
Equation (25) shows that the position of the MAUAV relative to the ground runway during the approach, glide slope, and flare is affected by the speed of vehicle.Considering that the ground speed of the MAUAV during landing with low wind disturbances mainly depends on the forward velocity u, it is necessary to control the forward velocity to make the MAUAV land efficiently and safely.Based on Equation ( 6), we can derive the forward velocity dynamic model as From Equations ( 7), ( 8), ( 10) and ( 26), the design model can be written as where According to the design model (27), the expansion state space can be expressed as .u = f u + g u u δ t ; .
where f u is designed as an extended state.
To estimate the two states of the new system, we construct the LESO as where û and fu represent the estimate of u and f u ; β u 1 and β u 2 are the tuning parameters and are determined by the value of the bandwidth of the LESO, which is denoted as ω u 0 .Perfect estimations of the true value of the states in the system can be achieved when adjusting these two gains by the following relationships [39]: Based on the accurate estimation of the system state, the final throttle command has the following form: where k u 1 = ω u c represents a positive feedback gain, ω u c is the desired closed-loop natural frequencies of the controller; u c is the desired forward velocity.The specific design idea of the forward velocity controller can be referred to Figure 3a.

Altitude Controller
A combination of lateral and longitudinal motion control is required during the landing of the MAUAV.In the longitudinal plane, the pitch angle is controlled to ensure that the MAUAV can track the altitude trajectory for the three processes of approach, glide slope, and flare.From Equation (25), we can know that the altitude of the MAUAV relative to the ground can be expressed as where Based on the altitude design model (33), the expansion state space can be expressed as . .
In a similar way to the design of the forward velocity control, the LESO and the desired pitch angle are designed as follows: .
where H c is the desired altitude; the values of the adjustment parameters β H 1 , β H 2 , and k H 1 are taken according to the following rules [39]:

Attitude Controller
In the previous section on the design of the altitude controller, we have obtained the control command for the pitch angle, making it as one of the desired output commands in the attitude controller.In addition, the goal of lateral motion control is to minimize the lateral error between the command trajectory and the position of the MAUAV.The lateral motion guidance method in this paper is derived from [1], which constructs a commanded course angle χ c at every spatial point relative to the desired straight-line path that results in the MAUAV moving toward the path.The commanded course angle is described by the following equation [1]: where r 0 = r x , r y , r z is the origin of the desired path; a = [a n , a e , a d ] is a unit vector whose direction indicates the desired direction of travel, where a n , a e , a d expresses the north, east, and down components of the unit direction vector.The straight-line path in this paper coincides with the centerline of the runway, and it is assumed that its origin r 0 = [−3016, 0, 0] and its direction along the inertial north direction, i.e., a = [1, 0, 0]; χ ∞ ∈ (0, π/2] is an approach angle for large distance from straight-line path, k path is a positive constant that influences the rate of the transition from χ ∞ to zero,e ay is the cross-track error.
X and the coordinated turn condition is described by where is the velocity of the airframe relative to the inertial frame.Then, the desired tracking control of course angle is achieved by having first-order error dynamics in χ as shown .
where error e χ = χ − χ c , k χ is positive constant and an expression for desired ϕ c is obtained as where ζ = χ − ψ is the crab angle.The crab angle is small when the MAUAV is flying in low wind condition, we can assume cos ζ ≈ 1 and ( 40) can be simplified as [40] Then, we have the desired change rate of yaw angle ψ c as [1] .
According to the rotational kinematic equations of the MAUAV above ( 16), the desired tracking control of attitude angles is achieved by having first-order error dynamics in Φ as shown .
where error Rearranging the above equations and substituting equation, desired rotation rates ω c = [p c , q c , r c ] T are designed as follows: According to the rotational dynamic equations of the MAUAV above (11), we get . .
where u ω = y b 1c , δ e , δ r T , f ω and g ω are given by where Then, the LESO and control input u ω are designed as follows: .
To facilitate the handling of the control problem of MAUAV system, an assumption is described as follows: .

Convergence Analysis of LESO
The four controllers designed above all employ LESO to estimate the states and disturbances of the system.Considering the similarity of the principles, an example analysis of the convergence of the LESO in the moving-mass position controller is given below, the convergence of the other three observers can be argued by the same method [41].50) and ( 52), the observer estimation error can be shown as . . where Theorem 1. Assuming .f y is bounded, there exists a constant σ i > 0 and a finite T 1 > 0 in order that | x i (t)| ≤ σ i , i = 1, 2, 3, ∀t ≥ T 1 > 0, and ω , for some positive integer k.
Proof of Theorem 1. Solving (55), it follows that Let f y is bounded, namely, .

Stabitily Analysis of LADRC
Theorem 2. Assuming .f y is bounded, there exists a constant ρ i > 0 and a finite T 3 > 0 in order for some positive integer j.
Proof of Theorem 2. From (52) one has According to (68) and Theorem 1, one has where e A e (t−τ) A x x(τ)dτ.It follows that and Since A e is Hurwitz, there exists a finite time T 2 > 0 in order that for all t ≥ T 3 , i = 1, 2, and for all t ≥ T 3 .From (71), (72), and (75), we obtain It can be derived by the same method that the tracking errors between desired and actual outputs in subsystems ( 29), (34), and (45) are also bounded after a period of time.

Numerical Simulation Setup
The task of this section is to verify the effectiveness of the proposed control scheme for the MAUAV landing process using numerical simulations.A sampling period of 0. Our control target is to enable the MAUAV to fly along the designed path and eventually land on the designated location on the ground runway (inertial coordinates of the touchdown point are [X td , Y td , H td ] = [0, 0, 0]) under the coordinated work of forward velocity, attitude, altitude, and moving-mass position controllers.In order to achieve this goal, it is necessary to give a reasonable expected path, expected forward velocity, and to select the appropriate controller parameters.Clearly, the desired lateral displacement command is Y c = 0, and the desired height trajectory is designed as follows: Specifically, in the approach stage, the height of the MAUAV is kept constant at 100 m; in the glide slope stage, the MAUAV declines along a straight line with an inclination angle of about −2.6 • ; in the flare stage, the MAUAV declines along an exponential curve, and the leveling decision height is about 15 m.To improve the tracking accuracy in the inertial X-axis direction of the MAUAV when landing on the runway, the commanded value of forward velocity is designed as u c = −X(t)/(75.4− t) when t < 75.4 s.When t = 75.4s, the commanded forward velocity is kept the same as the value taken at 75.38 s to avoid singularity.In addition, the desired three attitude angles are determined based on the desired flight trajectory.The desired roll angle ϕ c is obtained based on Equation (40), and the desired yaw angle ψ c is updated by Equation (41).The desired pitch angle θ c is determined by the height trajectory and its value is selected according to Equation (34).
In order to obtain an excellent control effect, it is necessary to select the appropriate regulation parameters for the designed controllers.The bandwidth and the desired closed-loop natural frequencies for the forward velocity, altitude, attitude, and moving-mass position controllers are chosen as:

Simulation Results
Simulations were performed using the physical and control parameters described above, and the final simulation results are shown in Figures 4-9. Figure 4 shows the forward velocity of the MAUAV and the throttle control size of the engine.The propeller throttle opening size is adjusted according to the designed forward velocity controller (31) to quickly control the forward velocity to the desired value.In Figure 5, we can find the slider offset, elevator offset, and rudder offset designed by the attitude controller.As can be seen from the figure, the signals are relatively smooth and continuous, except for some slight oscillations in the local area.The driving force designed by the moving-mass position controller can be found in Figure 6.As shown, the driving force is mainly activated in the first 10 s, since the lateral deflection of MAUAV is eliminated by adjusting the roll angle and a certain amount of force is needed to quickly move the slider, and thus change the roll angle.After the lateral deflection is close to 0 m, the drive maintains the wing level through a small range of outputs under the influence of wind disturbance.Figure 7a-c represent the roll angle, pitch angle, and yaw angle, respectively.During the whole landing process, the three attitude angles are smooth and continuous.The three-dimensional trajectory of the MAUAV during landing is shown in Figure 8.In order to demonstrate the accuracy of the landing process more intuitively, we show the path of MAUAV in X-H planes and lateral deviation in Figure 9.It is not difficult to find that the actual path is consistent with the desired path throughout the landing process except for the approach phase, which is due to the relatively large airspeed fluctuations in this phase.In addition, the tracking accuracy in the X and H directions of the MAUAV eventually reach errors of about 0.01 m and 0.2 m, respectively.In addition, the lateral deviation is rapidly reduced during the approach stage and finally converged to the runway centerline at about 15 s.When the MAUAV lands on the runway, the lateral deviation of it reaches about 0.02 m.Overall, the effectiveness of the proposed four controllers can be demonstrated by the above simulation results.These controllers work together to assist the MAUAV to perform landing tasks according to the desired path and ultimately land on the runway precisely.
during the approach stage and finally converged to the runway centerline at abo When the MAUAV lands on the runway, the lateral deviation of it reaches about Overall, the effectiveness of the proposed four controllers can be demonstrated above simulation results.These controllers work together to assist the MAUAV to p landing tasks according to the desired path and ultimately land on the runway pr   during the approach stage and finally converged to the runway centerline at abo When the MAUAV lands on the runway, the lateral deviation of it reaches about Overall, the effectiveness of the proposed four controllers can be demonstrated above simulation results.These controllers work together to assist the MAUAV to p landing tasks according to the desired path and ultimately land on the runway pr

Conclusions
In this paper, the automatic landing control of a fixed-wing UAV with movin rather than ailerons has been studied.Considering the specificity of the structure studied UAV, its control design models were established based on the analysis of dy models, and then four controllers were designed based on these design models LADRC.Among them, the forward velocity controller is designed to make the M

Conclusions
In this paper, the automatic landing control of a fixed-wing UAV with movin rather than ailerons has been studied.Considering the specificity of the structur studied UAV, its control design models were established based on the analysis of d models, and then four controllers were designed based on these design model LADRC.Among them, the forward velocity controller is designed to make the M

Conclusions
In this paper, the automatic landing control of a fixed-wing UAV with moving mass rather than ailerons has been studied.Considering the specificity of the structure of the studied UAV, its control design models were established based on the analysis of dynamic models, and then four controllers were designed based on these design models using LADRC.Among them, the forward velocity controller is designed to make the MAUAV reach the desired speed quickly from the initial speed by controlling the engine oil port opening.The designed altitude controller is able to generate the corresponding pitch angle signal based on the desired altitude trajectory and serve as the pitch angle tracking command in the attitude controller.The designed attitude controller can ensure that the MAUAV moves at the desired attitude angles, and its generated mass block movement signal serves as a tracking command for the moving-mass position controller.Then, the movement of the slider is achieved by controlling the magnitude of the driving force, which is taken by the designed moving-mass position controller.The simulation results demonstrate the accuracy and effectiveness of the proposed control strategy.Specifically, the synergy of the four controllers allows the MAUAV to land on the runway according to a given trajectory.Not only the actual longitudinal trajectory matches the desired altitude, but also the lateral deviation of the MAUAV when landing on the runway is small.The design of the controllers in this paper is based on LADRC, and thus the control algorithm itself lacks novelty.In view of this shortcoming, we will endeavor to design new algorithms, which takes into account the characteristics of MAUAVs, to improve the safety and accuracy of the landing process in our future researches.

Figure 1 .
Figure 1.Schematic diagram of the moving-mass actuated UAV.

Figure 1 .
Figure 1.Schematic diagram of the moving-mass actuated UAV.
nonlinear static function is used to describe different friction effects and can be parameterized to characterize the Stribeck effect.Typically, the function is modeled by g .+ ( f s − f c )e −( .Stribeck velocity; f c and f s are the levels of the normalized Coulomb friction and stiction force, respectively; and sgn is a standard signum function.
the Stribeck velocity; f c and f s are the levels of the normalized Coulomb fric- tion and stiction force, respectively; and sgn is a standard signum function.

Figure 2 .
Figure 2. Simplified block diagram of the automatic landing system for MAUAV.

Figure 2 .Figure 3 .
Figure 2. Simplified block diagram of the automatic landing system for MAUAV.

Figure 3 .
Figure 3. Block diagrams of the controller design based on LADRC: (a) first-order; (b) second-order.

Figure 4 .
Figure 4. (a) Forward velocity; (b) throttle control size of the engine.

Figure 4 .
Figure 4. (a) Forward velocity; (b) throttle control size of the engine.

Figure 4 .
Figure 4. (a) Forward velocity; (b) throttle control size of the engine.

Figure 8 .Figure 8 .
Figure 8.The trajectory of the MAUAV during the whole landing process.