Precision Landing of Unmanned Aerial Vehicle under Wind Disturbance Using Derivative Sliding Mode Nonlinear Disturbance Observer-Based Control Method

: Unmanned aerial vehicles (UAVs) are extensively employed in civilian and military applications because of their excellent maneuverability. Achieving fully autonomous quadrotor flight and precision landing on a wireless charging station in the presence of wind disturbance has become a crucial research topic. This paper presents a composite control technique for UAV altitude and attitude tracking in harsh environments, i.e., wind disturbance. A composite controller was developed based on nonlinear disturbance observer (NDOB) control theory to allow the UAV to land in the presence of random external wind disturbances and ground effects. The NDOB estimated the unknown wind disturbance, and the estimation was fed into the derivative sliding mode nonlinear disturbance observer-based control (DSMNDOBC), allowing the UAV to perform autonomous precision landing. Two loop designs were applied: the inner loop for stabilization and the outer loop for altitude tracking. The quadrotor model dynamics and the proposed controller, DSMNDOBC, were simulated employing MATLAB/Simulink ® , and the results were compared with the one obtained by the proportional derivative (PD) controller and the sliding mode controller (SMC). The simulation results indicated that the DSMNDOBC has superior altitude and attitude control compared to the PD and SMC controllers and better disturbance estimation and attenuation performance.


Introduction 1.Background
Unmanned aerial vehicles (UAVs) are utilized for various purposes, including military, civil engineering, and scientific applications such as delivery services, aerial mapping, surveillance, risk zone inspections, and search and rescue operations [1].Researchers are paying considerable attention to UAVs because of their reliability, maneuverability, relevance in various applications, and affordability.
Quadrotors are an important subclass of UAVs with six degrees of freedom and four control inputs [2], typical of a coupled system with complex dynamics.This multipleinput and multiple-output system is highly nonlinear and underactuated, so excellent maneuverability is required for robust and autonomous mission completion, particularly during precision landing operations under the presence of wind disturbance for automated wireless charging mission.
The design and operation of complex systems involve two separate but connected fields: control systems and communication frameworks.A control system modifies variables to produce the desired performance or behavior.Typically, it comprises a plant, a controller, actuators, and sensors (the system being controlled) to attain desirable set-point values, stabilize the system's behavior, and reduce the impact of disturbances.In contrast, the communication framework concerns the communication infrastructure, such as message formatting, error detection and repair, and data transport.Communication Aerospace 2024, 11, 265 2 of 14 frameworks are generic and can be applied to many systems; however, control systems are intended for specific systems and applications.

Literature Review
Precisely landing a UAV at a wireless charging station (moving or stationary) in a harsh environment, particularly when random wind disturbance is present, is a daunting task.To land precisely on a target, a UAV must accurately follow its trajectory while maintaining stability.When designing precision landing controllers for quadrotors, complications such as ground effects and external disturbances (e.g., wind gusts) can arise [3].To address this control problem, researchers have employed various control algorithms such as the proportional-integral-derivative (PID) controller [4], backstepping method [5], sliding mode control (SMC) [6], model predictive control (MPC) [7], and higher-order SMC [8].However, many control techniques, such as PID control, SMC, and MPC, have primarily been employed to stabilize quadrotors equipped with numerous sensors, inertial measurement units, and cameras.Additional sensor installments on a UAV can result in the quadrotor system gaining additional weight and changing its center of gravity, mass, and inertia, leading to a higher chance of instability.Without eliminating the characteristics of ground effects and exterior disturbances, which are normally mixed with the control inputs, the prelisted controllers become ineffective in properly stabilizing the UAV's flight motion, resulting in unexpected accidents.The disturbances must be measured to ensure the stability and resilience of the system.Therefore, designing high-performance quadrotor controllers remains challenging for researchers [9].
An optical-flow-based approach that adjusts controller gains was proposed for UAV landing [10].A velocity vector field method was proposed in [11].However, these studies were based on UAVs landing on fixed flat platforms.Cabecinhas et al. [12] proposed a robust control method that enabled the landing of a quadrotor on a slope.Ho et al. [13] introduced a method for landing UAVs by addressing the essential selection gain, a problem associated with optical-flow-based UAV landings.To overcome this challenge, Ho et al. [13] utilized a camera to estimate the ground orientation and an adaptive gain selection controller.
Algorithms for landing on stationary platforms have attracted research attention, but they cannot achieve landing in motion.The challenge of landing on mobile platforms was studied in [14], and a vision-based solution was proposed to address the challenge of a quadrotor landing on a vertically moving platform.However, the study was irrelevant to outdoor activities, because a motion capture technique was adopted.In another study [15], global positioning system (GPS) navigation was employed to allow a quadrotor to land employing a backstepping controller.Positioning is the primary factor that significantly influences the effectiveness of UAVs landing.Despite being widely utilized for locating UAVs outdoors [16], GPS suffers significantly in clustering situations, such as densely populated cities, confined valleys, and dense woods.
In addition, in caves and interior spaces with poor GPS reception, ultrasonic rangers, laser rangers, light detection and ranging sensors, and visual cameras are the additional onboard sensors that are most often utilized by researchers [17].An online UAV altitude control technique is proposed to minimize the expected weighted sum age of information for Internet of Things (IoT) devices, with timely data delivery before the data become outdated and lose value [18].Moreover, it leverages the IoT for an efficient and fast transfer of data between the mission computer and the UAV [19,20].By contrast, other methods focus on improving the landing target state estimation.
All the studies discussed above emphasized the role of a robust controller in the precision landing of quadrotors.When designing a robust controller during the precision landing procedure, it is necessary to carefully analyze external random wind disturbances and ground effects, which are the main variables affecting the low-altitude flying performance of quadrotors.

Contribution
This paper presents a derivative-sliding-mode nonlinear disturbance-observer-based control (DSMNDOBC) method for altitude and attitude control of UAVs in the presence of ground effects and random wind disturbances.The main contributions of this study are as follows: • We explored a nonlinear disturbance observer (NDOB) to design a DSMNDOBC.A derivative function is combined with the SMC to achieve dynamic and robust control.
The derivative function can cancel out the overshoot effect, allowing the modeled system to achieve a satisfactory control performance.The SMC component has a fast dynamic response and immunity to changes in plant factors and provides an integration platform for estimated disturbance, such as random wind disturbance.• The integrated NDOB has a low computational demand, with no need for supplemen- tary sensors that can estimate and generate corrective control inputs for the annulment of random external disturbance effects.

•
In numerical simulations, the proposed DSMNDOBC outperformed conventional proportional-derivative (PD) and SMC controllers with regard to altitude and attitude tracking.

Article Organization
The remainder of this paper is organized as follows: Section 2 presents the general quadrotor model.Section 3 presents the problem formulation, adopting the DSMNDOBC design process for altitude and attitude tracking.The simulation results are presented in Section 4. Finally, the conclusions and future research directions are presented in Section 5.

Precision Landing Control
Considering the combined effects of external disturbances such as low visibility, wind gusts, and risky deck motion, it is difficult to satisfy the performance requirements of precision landing.The precise landing of UAVs requires deploying a strong controller that can cancel the effects of inherent and exogenous disturbances while providing stability and a good dynamic response.To build a strong control system for precision landing, external disturbances and ground effects were considered by employing a nonlinear disturbance-observer-based control (NDOBC).The NDOBC can handle uncertainty, underactuation, and disturbances in the system and can be seamlessly integrated with many nonlinear control laws to achieve robustness [22,23].
The NDOBC design approach was as follows: (1) A nonlinear composite controller centered on the SMC was designed to attain stability and satisfy additional performance requirements, with the hypothesis that the disturbance is quantifiable.(2) An NDOB was then built to estimate the disturbance.
(3) The disturbance observer was integrated with a nonlinear controller by substituting the observer's disturbance estimation into the control law.
Considering an unknown external wind disturbance (with unknown amplitude), the equations of motion for the quadrotor along the vertical, roll, pitch, and yaw axes are as follows:

Precision Landing Control
Considering the combined effects of external disturbances such as low visibility, wind gusts, and risky deck motion, it is difficult to satisfy the performance requirements of precision landing.The precise landing of UAVs requires deploying a strong controller that can cancel the effects of inherent and exogenous disturbances while providing stability and a good dynamic response.To build a strong control system for precision landing, external disturbances and ground effects were considered by employing a nonlinear disturbanceobserver-based control (NDOBC).The NDOBC can handle uncertainty, underactuation, and disturbances in the system and can be seamlessly integrated with many nonlinear control laws to achieve robustness [22,23].
The NDOBC design approach was as follows: (1) A nonlinear composite controller centered on the SMC was designed to attain stability and satisfy additional performance requirements, with the hypothesis that the disturbance is quantifiable.(2) An NDOB was then built to estimate the disturbance.
(3) The disturbance observer was integrated with a nonlinear controller by substituting the observer's disturbance estimation into the control law.
Considering an unknown external wind disturbance (with unknown amplitude), the equations of motion for the quadrotor along the vertical, roll, pitch, and yaw axes are as follows: .. .. .. ..
where d z , d ϕ , d θ , and d ψ denote the wind disturbances (m/s and rad/s); r represents the rotor radius (m); d h represents the horizontal distance between two opposite rotor axes in a diagonal direction (m); d v represents the vertical distance from the propeller surface to the ground (m); and δ denotes the ground effects (no unit).The additional term in the denominator represents the aerodynamic interference of other rotors [24].The NDOB was introduced to estimate unknown disturbances, adopting the internal state variable z ndo , as described in the following section.

DSMNDOBC
The detailed mathematical derivations for achieving DSMNDOBC are presented in Figure 2.
the rotor radius (m);  represents the horizontal distance between two opposite rotor axes in a diagonal direction (m);  represents the vertical distance from the propeller surface to the ground (m); and  denotes the ground effects (no unit).The additional term in the denominator represents the aerodynamic interference of other rotors [24].The NDOB was introduced to estimate unknown disturbances, adopting the internal state variable  , as described in the following section.

DSMNDOBC
The detailed mathematical derivations for achieving DSMNDOBC are presented in Figure 2.
where , , , and  denote the system state (no unit), the input vector (kg • m/s ), exogenous wind disturbance (m/s), and the output vector (no unit), respectively.The terms ,  ,  , and ℎ are smooth functions with respect to .Furthermore, an external disturbance  can be generated by the system as follows: The general form of the nonlinear disturbance observer is (1) NDOB Consider a nonlinear system of the form ..
where X, U, d, and Y denote the system state (no unit), the input vector (kg•m/s 2 ), exogenous wind disturbance (m/s), and the output vector (no unit), respectively.The terms f , g 1 , g 2 , and h are smooth functions with respect to x.Furthermore, an external disturbance d can be generated by the system as follows: .
The general form of the nonlinear disturbance observer is .
where z ndo is the internal state variable of the NDOB (m), ξ is the auxiliary variable of the observer (m), p and η denote the nonlinear function and the observer gain to be selected, respectively (no unit), d denotes the disturbance observer estimation (m/s), A is an m × m matrix representing the frequency of the external disturbance (s −1 ) [22], and C is a 1 × m row matrix (no unit).Here, A and C are defined as and the nonlinear observer gain η is estimated as follows: .
where e ndo represents the estimation error of the NDOB (s −1 ).The above equations ensure the global exponential stability of the closed-loop system, regardless of disturbances [25].
Here, the concept of the NDOB is briefly presented.Comprehensive theoretical details can be found in [22,25], where the NDOB is defined for the axes being considered.
z ndo,z can be obtained as .
Similarly, the NDOBs for the attitude dynamics are obtained as follows: . . .
(2) DSMC Derivative sliding mode control (DSMC) is adopted for a quadrotor's altitude and attitude trajectory tracking under the influence of ground effects and random wind disturbances.The tracking error for the motion along the vertical axis is given by For any constant k 1 , the sliding mode surface and its derivative are given as .
The controller inputs U i (i = 1, . . ., 4) are a composite of the sliding mode control law with a derivative component: (39)

PD Controller
Position tracking was achieved by implementing PD controllers as described below.The tracking errors for the x and y positions are defined as follows: .
x and U y = k p,y e y − k d,y . y.
From Equation ( 8), θ d was derived as Similarly, from Equations ( 9) and ( 10), ϕ d was computed as follows: where (x d , y d , ϕ d , θ d ) and (x, y, ϕ, θ) are the desired and current values, respectively, and k p and k d are the proportional and derivative gains to be determined (no unit).

Environmental Setup
The performance of the proposed composite control technique, DSMNDOBC, was evaluated utilizing numerical simulations in a MATLAB/Simulink ® environment.The results of the DSMNDOBC method were compared with those of the SMC and PD control methods.The performance of the controllers was judged by taking the root-mean-square error of the quadrotor trajectory tracking in the z, ϕ, θ, and ψ directions.
A random wind profile was generated, with mean wind speeds of 10 and 20 m/s, utilizing the Dryden continuous wind turbulence model [26].A Dryden wind model block was employed to generate the atmospheric turbulence.White noise was passed through a filter to provide turbulence to the specified velocity spectra.The effects of turbulence on a moving body such as an aircraft were modeled by producing time-varying wind gusts.
In MATLAB/Simulink ® , the Dryden wind model was implemented as a block in the Simulink Library.In this Simulink block, a parameter called "mean wind speed" (representing the resultant wind velocity in the body axis frame) was employed to output time-varying wind gusts (turbulence velocities in x, y, and z and turbulence angular rates in ϕ, θ, and ψ).The model utilized a random number generator to produce a wind disturbance.In our simulations, mean wind speeds of 10 and 20 m/s were utilized, corresponding to maximum wind turbulence velocity amplitudes (in the z-axis) of 5.24 and 10.27 m/s, respectively.
Each DSMNDOBC was subjected to a sample wind profile, as illustrated in Figure 3. Figure 4 presents the disturbance estimates, with the generated disturbances indicated by red lines.The NDOB accurately estimated the disturbances in the system and contributed to the overall performance of the proposed controller.
corresponding to maximum wind turbulence velocity amplitudes (in the z-axis) of 5.24 and 10.27 m/s, respectively.
Each DSMNDOBC was subjected to a sample wind profile, as illustrated in Figure 3. Figure 4 presents the disturbance estimates, with the generated disturbances indicated by red lines.The NDOB accurately estimated the disturbances in the system and contributed to the overall performance of the proposed controller.

Scenario
A total of three scenarios based on two different trajectories (linear and circular) were employed to evaluate the three control methods (DSMNDOBC, PD, and SMC) by introducing three uncertainty media (sensor noise, model uncertainty, and wind turbulence).The quality of the trajectory tracking response of the quadrotor was evaluated based on the root-mean-square error ( ) and maximum error ( ) along the altitude and attitude positions.
(1) Scenario 1 In the first scenario, the performance of the DSMNDOBC was evaluated while the quadrotor tracked a predetermined line trajectory at an altitude of 2 m under wind disturbance conditions of 10 and 20 m/s for 80 s with a sample time of 0.01 s.The  and  results are presented in Tables 1 and 2, respectively, and the simulation results are presented in Figures 5-7.corresponding to maximum wind turbulence velocity amplitudes (in the z-axis) of 5.24 and 10.27 m/s, respectively.Each DSMNDOBC was subjected to a sample wind profile, as illustrated in Figure 3. Figure 4 presents the disturbance estimates, with the generated disturbances indicated by red lines.The NDOB accurately estimated the disturbances in the system and contributed to the overall performance of the proposed controller.

Scenario
A total of three scenarios based on two different trajectories (linear and circular) were employed to evaluate the three control methods (DSMNDOBC, PD, and SMC) by introducing three uncertainty media (sensor noise, model uncertainty, and wind turbulence).The quality of the trajectory tracking response of the quadrotor was evaluated based on the root-mean-square error ( ) and maximum error ( ) along the altitude and attitude positions.
(1) Scenario 1 In the first scenario, the performance of the DSMNDOBC was evaluated while the quadrotor tracked a predetermined line trajectory at an altitude of 2 m under wind disturbance conditions of 10 and 20 m/s for 80 s with a sample time of 0.01 s.The  and  results are presented in Tables 1 and 2, respectively, and the simulation results are presented in Figures 5-7.

Scenario
A total of three scenarios based on two different trajectories (linear and circular) were employed to evaluate the three control methods (DSMNDOBC, PD, and SMC) by introducing three uncertainty media (sensor noise, model uncertainty, and wind turbulence).The quality of the trajectory tracking response of the quadrotor was evaluated based on the root-mean-square error (E RMS ) and maximum error (E MAX ) along the altitude and attitude positions.
(1) Scenario 1 In the first scenario, the performance of the DSMNDOBC was evaluated while the quadrotor tracked a predetermined line trajectory at an altitude of 2 m under wind disturbance conditions of 10 and 20 m/s for 80 s with a sample time of 0.01 s.The E RMS and E MAX results are presented in Tables 1 and 2, respectively, and the simulation results are presented in Figures 5-7.(2) Scenario 2 In the second scenario, a circular trajectory was adopted to assess the performance of the proposed DSMNDOBC under two environmental conditions at the same wind speed.The time-varying feature of a circular trajectory with a sample rate of 100 is mathematically described as ω = 0.01 : π 100 : 2π, [x, y, z] = [20sin ω, 20cos ω, 1.9], and T = 0 : 0.1 : 21.The transient response was tracked along the desired height of 1.9 m with a circular trajectory of a radius of 20 m.The simulation time was 2100 s.
Tables 3 and 4 present the E RMS and E MAX values for the altitude and attitude positions, respectively, over the simulation period.This result reinforces the fact that the DSMNDOBC provides excellent tracking performance compared to the PD and SMC controllers.The simulation results are presented in Figures 8-10.(2) Scenario 2 In the second scenario, a circular trajectory was adopted to assess the performance of the proposed DSMNDOBC under two environmental conditions at the same wind speed.The time-varying feature of a circular trajectory with a sample rate of 100 is mathematically described as  = 0.01: : 2, , ,  = 20 sin  , 20 cos  , 1.9 , and  = 0: 0.1: 21.
The transient response was tracked along the desired height of 1.9 m with a circular trajectory of a radius of 20 m.The simulation time was 2100 s.Tables 3 and 4 present the  and  values for the altitude and attitude positions, respectively, over the simulation period.This result reinforces the fact that the DSMNDOBC provides excellent tracking performance compared to the PD and SMC controllers.The simulation results are presented in Figures 8-10.(3) Scenario 3 In the third scenario, for further robustness analysis of the three control methods, three uncertainty media (sensor noise, model uncertainty, and wind turbulence) were combined and introduced into the simulation with the same line trajectory as in Scenario 1.
First, white Gaussian noise with zero mean was created to simulate sensor noise.For a sampling time of 0.01 s, the variance in the generated noise for the altitude motion was 0.5 and 0.2 for attitude motion tracking, considering that a GPS sensor could have an accuracy of 0.025 m [27].
Second, for the introduction of model uncertainties, a variation of ±50% was assumed for the mass of the quadrotor (), and the inertia ( ,  ,  ) was allowed to vary by ±20%.
Third, wind turbulence with an amplitude of 10.14 / was introduced to analyze the robustness against exogenous wind disturbances.
The trajectory performance analysis results are presented in Tables 5 and 6.   (3) Scenario 3 In the third scenario, for further robustness analysis of the three control methods, three uncertainty media (sensor noise, model uncertainty, and wind turbulence) were combined and introduced into the simulation with the same line trajectory as in Scenario 1.
First, white Gaussian noise with zero mean was created to simulate sensor noise.For a sampling time of 0.01 s, the variance in the generated noise for the altitude motion was 0.5 and 0.2 for attitude motion tracking, considering that a GPS sensor could have an accuracy of 0.025 m [27].
Second, for the introduction of model uncertainties, a variation of ±50% was assumed for the mass of the quadrotor (), and the inertia ( ,  ,  ) was allowed to vary by ±20%.
Third, wind turbulence with an amplitude of 10.14 / was introduced to analyze the robustness against exogenous wind disturbances.
The trajectory performance analysis results are presented in Tables 5 and 6.  (3) Scenario 3 In the third scenario, for further robustness analysis of the three control methods, three uncertainty media (sensor noise, model uncertainty, and wind turbulence) were combined and introduced into the simulation with the same line trajectory as in Scenario 1.
First, white Gaussian noise with zero mean was created to simulate sensor noise.For a sampling time of 0.01 s, the variance in the generated noise for the altitude motion was 0.5 and 0.2 for attitude motion tracking, considering that a GPS sensor could have an accuracy of 0.025 m [27].
Second, for the introduction of model uncertainties, a variation of ±50% was assumed for the mass of the quadrotor (M), and the inertia (I x , I y , I z ) was allowed to vary by ±20%.
Third, wind turbulence with an amplitude of 10.14 m/s was introduced to analyze the robustness against exogenous wind disturbances.
The trajectory performance analysis results are presented in Tables 5 and 6.The proposed controller was further quantified utilizing the mean wind speeds of 22, 24, 26, and 28 m/s to determine the maximum disturbance boundary for the controller.The maximum disturbance bound was 26 m/s.The system diverged at 27 m/s, corresponding to a wind turbulence velocity amplitude of 13.96 m/s.
Although DSMNDOBC seems to result in similar values of both E RMS and E MAX under ±50% of the mass variation, a big difference is observed in both the PD and SMC cases for all investigation samples (z, ϕ, θ, and ψ).This indicates that the DSMNDOBC has superior control, disturbance estimation, and attenuation performance compared with the PD and SMC controllers.

Conclusions
Here, a composite controller was designed utilizing the NDOBC approach to estimate and attenuate the exogenous (wind) and ground effect disturbances that negatively influence quadrotor dynamics.Because a quadrotor is an underactuated system, two loop designs were considered in the algorithm.Stabilization was performed in the inner loop and altitude tracking was performed in the outer loop.
In simulations based on the MATLAB/Simulink ® environment, the proposed composite control technique exhibited a robust performance, as it moved precisely along the specified route when subjected to various wind disturbances.The simulation results indicate that the DSMNDOBC has superior altitude and attitude control and better disturbance estimation and attenuation performance than the PD and SMC controllers.The performance of the DSMNOBC was better, more stable, and more robust.Overall, a UAV with a 2.7 kg mass subjected to a fairly violent turbulence profile provides such an attenuated response, especially considering the sensor noise and model uncertainty in close proximity to the ground.
In the future, the designed controller will be tested on an experimental quadrotor platform for multiple applications, such as autonomous precision landing in harsh environments.Also, we plan to replace the PD controller with the DSMNDOBC to improve the x and y position accuracy by solving the unwanted disturbance estimation errors, bringing improved accuracy to the estimation of roll and pitch angles.
where x, y, and z represent the position of the quadrotor in the inertia frame {E} (m); ϕ, θ, and ψ represent the orientation of the quadrotor in the inertial frame {E} (rad); M represents the vehicle mass (kg); g represents the gravitational acceleration (m/s 2 ); l represents the arm length (m); I x , I y , and I z represent the moment of inertia (kg • m 2 ); and C f m represents the force-to-moment coefficient (no unit).

𝑀
where , , and  represent the position of the quadrotor in the inertia frame {E} (m); , , and  represent the orientation of the quadrotor in the inertial frame {E} (rad);  represents the vehicle mass (kg);  represents the gravitational acceleration (m/s );  represents the arm length (m);  ,  , and  represent the moment of inertia (kg ⋅ m ); and  represents the force-to-moment coefficient (no unit).

Figure 4 .
Figure 4. Wind disturbance estimation of the wind profile generated with a mean wind speed of 20 m/s: (a)  , (b)  , (c)  , and (d)  .

Figure 4 .
Figure 4. Wind disturbance estimation of the wind profile generated with a mean wind speed of 20 m/s: (a)  , (b)  , (c)  , and (d)  .

Figure 4 .
Figure 4. Wind disturbance estimation of the wind profile generated with a mean wind speed of 20 m/s: (a) d z , (b) d ϕ , (c) d θ , and (d) d ψ .

Figure 5 .
Figure 5. Scenario 1 simulation results for altitude and attitude position tracking under mean wind speed of 10 m/s: (a) , (b) , (c) , and (d) .

Figure 6 .
Figure 6.Scenario 1 simulation results for altitude and attitude position tracking under mean wind speed of 20 m/s: (a) , (b) , (c) , and (d) .

Figure 7 .
Figure 7. Three-dimensional trajectory tracking performance for Scenario 1 in the presence of random wind disturbance: (a) 10 m/s and (b) 20 m/s.

Figure 5 .
Figure 5. Scenario 1 simulation results for altitude and attitude position tracking under mean wind speed of 10 m/s: (a) , (b) , (c) , and (d) .

Figure 6 .
Figure 6.Scenario 1 simulation results for altitude and attitude position tracking under mean wind speed of 20 m/s: (a) , (b) , (c) , and (d) .

Figure 7 .
Figure 7. Three-dimensional trajectory tracking performance for Scenario 1 in the presence of random wind disturbance: (a) 10 m/s and (b) 20 m/s.

Figure 5 .
Figure 5. Scenario 1 simulation results for altitude and attitude position tracking under mean wind speed of 10 m/s: (a) , (b) , (c) , and (d) .

Figure 6 .
Figure 6.Scenario 1 simulation results for altitude and attitude position tracking under mean wind speed of 20 m/s: (a) , (b) , (c) , and (d) .

Figure 7 .
Figure 7. Three-dimensional trajectory tracking performance for Scenario 1 in the presence of random wind disturbance: (a) 10 m/s and (b) 20 m/s.

Figure 7 .
Figure 7. Three-dimensional trajectory tracking performance for Scenario 1 in the presence of random wind disturbance: (a) 10 m/s and (b) 20 m/s.

Figure 8 .
Figure 8. Scenario 2 simulation results for altitude and attitude position tracking under mean wind speed of 10 m/s: (a) , (b) , (c) , and (d) .

Figure 9 .
Figure 9. Scenario 2 simulation results for altitude and attitude position tracking under mean wind speed of 20 m/s: (a) , (b) , (c) , and (d) .

Figure 10 .
Figure 10.Three-dimensional trajectory tracking performance for Scenario 2 in the presence of random wind disturbance: (a) 10 m/s and (b) 20 m/s.

Figure 10 .
Figure 10.Three-dimensional trajectory tracking performance for Scenario 2 in the presence of random wind disturbance: (a) 10 m/s and (b) 20 m/s.

Figure 10 .
Figure 10.Three-dimensional trajectory tracking performance for Scenario 2 in the presence of random wind disturbance: (a) 10 m/s and (b) 20 m/s.