Trajectory-Following Control of an Unmanned Aerial–Aquatic Vehicle under Complex Coupling Interferences

: This article explores trajectory-following control for an unmanned aerial–aquatic vehicle (UAAV) navigating complex ocean disturbances and the interplay of air–seawater coupling factors. First, leveraging the backstepping technology, an adaptive algorithm is proposed to tackle the attitude and position following. Additionally, a nonlinear observer is crafted to estimate complex ocean disturbances. The UAAV model, characterized by six degrees of freedom (DOF) and nonlinear properties, experiences significant pose changes when emerging from water, underscoring the critical importance of precise pose control. Finally, stability analysis and numerical simulations are demonstrated to verify the feasibility and validity of the proposed control strategies.


Introduction
Over the past few decades, the challenges associated with designing manned vehicles, characterized by complex layouts, high costs, and low efficiency, have prompted a rapid and distinct development of unmanned vehicles.This evolution is particularly notable in the realms of unmanned aerial vehicles (UAVs) and unmanned underwater vehicles (UUVs) [1][2][3].To enhance the exploration and utilization of marine resources, as well as to meet military requirements, countries worldwide are consistently advancing diverse and innovative unmanned platforms [4].Traditional UAVs typically operate within a single medium, posing challenges in adapting to the progressively complex task requirements.This led to the emergence of an amphibious air-water cross-media vehicle capable of efficiently functioning in two distinct media conditions.
Proposing the military four-in-one safety concept encompassing sea, land, air, and space as the ultimate frontier of security underscores the paramount importance of UAV development.This is particularly significant, with a spotlight on the recently introduced cross-medium vehicle concept [5,6].On the civilian front, traditional underwater robots and drones deployed on ships face limitations, especially in tasks like marine resource surveys.The emergence of transmedia vehicles, however, provides a positive direction in this area.Additionally, the UAAV, as a cross-medium vehicle, offers a distinct advantage over underwater vehicles.It can transition into the air to reduce drag when high-speed travel is necessary [7].Hence, maneuver control of cross-medium vehicle systems has evolved into a captivating subject, given its vast potential applications in missile launches, environmental monitoring, and submarine surfacing [8][9][10].
In terms of motion control, the UAAV performs alternating maneuvers in both air and water, resulting in abrupt changes in propeller forces and rudder moments.In reference [11], an aerodynamic and hydrodynamic model of angular motion was established for a transmedia vehicle featuring a quadrotor layout.Building upon this, a method was introduced for stabilizing water emergence control of a transmedia vehicle in the longitudinal plane.The study delves into the impacts of additional mass and inertial forces on the water emergence process during transmedia motion.Moreover, disturbances in the liquid and intricate collisions between the vehicle and the liquid can negatively impact the motion control of the UAAV, therefore compromising the stability and overall performance of the vehicle.During the execution of a water-exit maneuver, a trajectory tracking controller is established according to [12] to meet the requirements for a large-angle water outlet.
The backstepping control method proves more effective in addressing the nonlinearity issue in motion control [13][14][15].Conversely, the adaptive control algorithm demonstrates a lesser sensitivity to changes in the model parameters of the controlled object.When combined with insights into control-method switching, the application of these methods to cross-medium vehicle control problems becomes profoundly significant for advancing research in cross-medium vehicle control.
The UAAV consistently contends with intricate external disturbances originating from wind, waves, and currents.Particularly at the water-air interface in oceans and lakes, it tends to be more susceptible to external environmental influences.Consequently, the maneuverability of the UAAV in such conditions is notably hampered by the coupling hydrodynamics.Throughout the water emergence process, the UAAV undergoes significant changes in fluid dynamics and is additionally influenced by external factors, notably wind and waves [16].The critical stage for the UAAV is the exit phase at the water surface and air juncture, where the external environment is exceptionally complex, leading to inevitable interference.Therefore, considering these external interferences, ensuring effective maneuver control of the vehicle becomes imperative.
To enhance the control capabilities, a disturbance observer was developed specifically for the disturbances caused by wind and waves [17][18][19].Building upon this, the environmental load associated with wind and waves in the subsequent moment was estimated and integrated into the vehicle's controller, enabling the determination of the vehicle's motion state.The interference observer not only assesses the external environmental load on the vehicle but also provides a level of compensation to enhance the controller's resistance to interference.To align the control research of the vehicle with the real scenario during this stage, adjustments were made to the kinematics and dynamics models of the vehicle in the exit phase.The interference forces and torques acting on the vehicle were meticulously analyzed and modeled [20].This paper concentrates on the vehicle's exit stage, where a trajectory tracking controller is designed.The interference estimation results are then compared with pre-designed interference to evaluate the interference estimator's effectiveness.Finally, simulation verification of the motion control process is conducted.
Currently, there is a scarcity of studies addressing maneuver control issues for UAAVs that account for complex disturbances at the air-seawater interface.The primary challenge in this domain lies in the fact that the unknown and intricate environment can hinder the vehicle's proper and smooth movement.Notably, the UAAV may experience thrill, leading to the unstable movement of the entire system.In response to these challenges, this article introduces an adaptive backstepping control algorithm integrated into a disturbance observer scheme for the UAAV system, considering the presence of coupling disturbances.In comparison to existing works, the innovations of this study can be summarized as follows: (1) We will introduce an adaptive backstepping control scheme designed for the fixedwing UAAV system.Unlike conventional backstepping approaches [14,15], our method dynamically accommodates changing system dynamics, enhancing its adaptability and performance in real-world applications.Moreover, the proposed control scheme exhibits remarkable robustness against complex coupling interferences, which holds significant importance in practical engineering applications.In contrast to PD control [11], backstepping control is a nonlinear technique that systematically addresses complex and nonlinear system dynamics.It employs a recursive design approach, developing the control law step by step, with a focus on canceling out nonlinearities and ensuring system stability.
(2) The kinematic and dynamic models of the UAAV are formulated by altering the order of Euler transformations to account for the motion characteristics during the waterexit process.Subsequently, the considerable nonlinearities and robust coupling inherent in the out-of-water phase, along with time-varying parameters, are addressed.The stability of UAAV control during the out-of-water motion phase was then verified.
(3) The disturbance force and moment resulting from wind and wave loads on the vehicle are analyzed and modeled.Given that the roll angle and bow angle can be approximated as small quantities, the fixed-wing cross-medium vehicle can be simplified as a 4-DOF model.Subsequently, a disturbance observer is designed to estimate wind and wave disturbance loads, which are then used as inputs for the external force of the UAAV dynamic system.Building upon the designed disturbance observer and controller, a mathematical model for wind and wave disturbances is established to facilitate stability analysis.
The structure for the remaining sections of this article is outlined as follows.Section 2 introduces the maneuvering model of the UAAV, encompassing both kinematic and kinetic models.Section 3 establishes models for external environmental disturbances, including winds and waves.The core processes of the motion control algorithm are detailed in Section 4. The effectiveness of the proposed control protocol is demonstrated through numerical simulations in Section 5. Finally, Section 6 provides concluding remarks.

The Maneuvering Model of UAAV
In this section, we will address a foundational concern by first presenting the kinematics and dynamics model of the UAAV.The kinematics and dynamics models of crossmedium vehicles can be categorized into three types based on the medium in which they operate: the underwater movement stage, the air flight stage, and the cross-medium emergence stage.In particular, the cross-medium emergence stage holds significance for the entire maneuver mission, serving as a cornerstone for subsequent research on control methods for cross-medium vehicles.

The Kinematic Model of UAAV
This subsection primarily delves into the motion modeling of the cross-medium vehicle, emphasizing that a well-founded motion model serves as the bedrock for subsequent research on control methods.The UAAV represents a novel aircraft capable of navigating in both air and underwater environments.Moreover, the precision of the motion model establishment directly impacts the accuracy of subsequent issues such as attitude control, speed control, and trajectory tracking.The kinematics and dynamics models of crossmedium vehicles are categorized based on the medium in which they operate, including the underwater movement stage, the air flight stage, and the cross-medium exit stage.
Modeling the UAAV system encompasses understanding both its dynamics and statics.Statics involves the equilibrium of the vehicle under uniform velocity or at a stationary position, while dynamics pertains to the accelerated motion of the vehicle.Kinematics focuses on the geometrical aspects of the vehicle's motion, whereas kinetics addresses the forces responsible for the motion of the UAAV.Consequently, the motion of the UAAV in geodetic and body coordinate systems is illustrated in Figure 1.In the delineation of the aforementioned coordinate system, the six degrees of freedom of the vehicle are described, encompassing both translational and rotational motion.Where O b − X b Y b Z b represents the body coordinate system, O E − X E Y E Z E represents the geodetic coordinate system, O w − X w Y w Z w represents the velocity coordinate system.
In the delineation of the aforementioned coordinate system, the six degrees of freedom of the vehicle are described, encompassing both translational and rotational motion.The underwater gliding of the vehicle and the description of its six degrees of freedom movement can be referenced using the symbol system of the International Towing Tank Conference (ITTC) and the Society of Naval Architecture and Marine Engineering (SNAME) [21].In particular, the motion model can be described as where J(η) stands for the rotation matrix between , which is defined by: T is used to describe the position of the UAAV.
the velocity in the body coordinate system, v 1 = {u v w} T stands for the linear velocity, and v 2 = {p q r} T represents the rotating angular velocity.s, c, t denote the sin, cos, tan, respectively.More specifically, x y z express the linear displacement of surge, sway, and heave direction, respectively.ϕ, θ, φ show the .rotationangle in the above three directions.u, ϑ, ω. express the corresponding angular spin rate.

The Dynamics Model of UAAV
Professor Fossen [22] from Norway proposed a simplified version in vector form for the dynamic equation of six degrees of freedom for rigid bodies, which can be expressed as: where M RB denotes the inertia matrix.C RB (v) represents the Coriolis force and centripetal force matrix.
T is the designed control force and moment, where X Y Z is the longitudinal force, the lateral force, and the vertical force of the vehicle, and K M N is the transverse moment, trim moment and yaw moment of the vehicle.Concretely, the inertia matrix can be expressed as follows: where M RB satisfies the symmetric property and is a constant matrix, i.e., M RB = M T RB > 0, . M RB = 0. Additionally, the Coriolis force matrix also can be shown as follows:

External Environmental Disturbance
The impact of wind and waves on the control of the cross-medium vehicle during movement involves two primary aspects.First, during takeoff at a large angle in the exit strategy, the vehicle is susceptible to the effects of wind and waves.In this process, the vehicle approaches a vertical state, and the interference from waves and wind on the broad wing and shape results in significant disruptions during the water discharge process.Second, the body's center of buoyancy changes with the fluctuation of draft caused by wave action, adversely affecting the vehicle's attitude to stay above the water.

Wind Force Model
The disturbance force and torque induced by the wind load on the vehicle in the wind field can be calculated using the following equation: where v r = v 1 + v 2 × r wp − R T v w is the relative velocity of the vehicle body and the wind field at the center of wind pressure.v 1 is the linear velocity of UAAV, v 2 is the relative angular velocity, R T is the rotation matrix of the vehicle from the inertial frame to the body coordinate system.C d,wind is the wind drag coefficient, representing the degree of influence of wind load on the resistance of the vehicle in the direction of the three axes of the UAAV system.

Wave Force Model
Drawing from the literature on linear wave theory [23], the load resulting from linear waves can be computed using pertinent knowledge from fluid mechanics [24].The wave load on the cylindrical body and wing of an aircraft can be determined by the velocity potential function, wave surface function, and wave surface pressure distribution function corresponding to the characteristics of small amplitude long wave: where ζ a is the wave height relative to the horizontal plane z = 0 in the inertial coordinate system, ω w is the circular frequency unit of a wave is rad•s −1 ; k w = 2π/λ is the wavenumber of a wave, representing the wavenumber in each period with wavelength being λ; The circular frequency ω w and wave number k w of a wave satisfy the dispersion relation k w = ω 2 w /g under wireless water depth.Given the considerably small scale of the UAAV in comparison to the long peak waves of the ocean, the impact of loads resulting from incident waves can be disregarded [25].The wave loads are primarily governed by the Froude-Krylov force, representing the perturbation force on an object under wave loading.This force is obtained through unit element re-integration when the wave pressure field is not influenced by the object in turn.Consequently, the disturbance force and moment due to wave loads on the vehicle can be calculated using the following equation: ∂r × p dyn ndS (11) where dS indicates that the UAAV is in the middle of a wave; S w denotes the wet surface area of the vehicle in the wave; n denotes the outward normal vector of the microelement; ∂r denotes the spatial position of the microelement in the inertial coordinate system.Up to this point, the interference force caused by wind and wave loads on the UAAV can be expressed using Equations ( 6) and (10), while the interference moment is expressed by Equations ( 7) and (11).Subsequently, the interference force and moment resulting from the corresponding assumed wave loads can be obtained using the following equation:

Control Algorithm Design
Given that the UAAV primarily focuses on positional quantities in the longitudinal plane during the water emergence process, the roll angle and bow angle can be approximated as small quantities, i.e., ϕ = ψ = 0 and p = r = 0.The kinematics and dynamics model can be simplified from Equations ( 1) and ( 3) to a mathematical model of longitudinal in-plane plus transverse displacement, i.e., a mathematical model of 4-DOF.The specific form of the state space can be expressed as: where x 1 = [x, y, z, θ] T is the state variable of the system, i.e., linear and angular displace- ments in the longitudinal plane in the geodetic coordinate system; x 2 = [u, v, w, q] T are the first order state variables of the system, i.e., linear and angular velocities in the longitudinal plane of the body coordinate system; y = x 1 is the output of the system that indicates the position and attitude state of the vehicle.u indicates the inputs to the control system, i.e., the inputs to the control forces and moments of the vehicle.d indicates the disturbance of the control system.In this paper, the three phases use a unified disturbance to avoid the critical point in each phase of the control force, which appears to jump to ensure that the response of the control system is continuous.B is the inverse of the quality matrix in the control system; f 1 (x 1 ) is the transformation matrix in the kinematic model equation in the case of 4-DOF, which is the transformation matrix of linear and angular velocities in the airframe coordinate system to the geodetic coordinate system position and attitude.f 2 (x 1 , x 2 ) is the product of the inverse of the mass matrix in the control system and the forces and moments, including static forces: gravity and buoyancy, hydrodynamic forces: inertial and viscous forces [26], in addition to the control forces and moments, as well as the slight disturbance forces designed in this paper, i.e., To facilitate a more effective analysis of the control laws during the out-of-water process, an equivalent approach is employed to conceptualize the navigator as comprising a cylindrical body and approximately rectangular wings.Define the distance from the center of mass to the upper and lower bottom surfaces as l u and l d , respectively.z is the height of the center of mass of the vehicle, defining z = 0 as the interface between water and air; z > 0 denotes air; z < 0 denotes water.The emergent state is partly in the water and partly in the air and is expressed as a range of z ∈ [−l d sin θ, l u sin θ].The craft completely in the water is expressed as z < −l u sin θ, and completely in the air is expressed as z > l d cos θ.Such situations are handled when the dynamics are written in a particular way, and a switched controller is designed [27] to tackle such dynamics.A schematic diagram of UAAV exiting the water across the medium in the presence of wind and wave disturbances is shown in Figure 2. In what follows, the specific steps of the backstepping method of designing a control law will be stated.
Step 1: The error variable and Lyapunov function of the trajectory tracking controller are first defined as (15) where denotes the desired exit trajectory of the navigator, for which the control system needs to guarantee a second-order continuously derivable; denotes the Lyapunov energy function of subsystem 1 used to determine whether the designed controller is asymptotically stable.Subsystem 1 is defined as (16) where is the virtual control law and can be designed as ( )( ) where is the designed control parameter and satisfy .
Define the new error variable , the derivative of the Lyapunov function of the system can be calculated as (18) The asymptotic stability of the controller requires that the Lyapunov function be negatively determined, so the new error variable must be zero, and the design of the controller subsystem 2.
Step 2: Define the new error variable and the corresponding Lyapunov function as (19) On the premise of ignoring the disturbance force, the actual system will be estab- In what follows, the specific steps of the backstepping method of designing a control law will be stated.
Step 1: The error variable and Lyapunov function of the trajectory tracking controller are first defined as where x 1d denotes the desired exit trajectory of the navigator, for which the control system needs to guarantee a second-order continuously derivable; V 1 denotes the Lyapunov energy function of subsystem 1 used to determine whether the designed controller is asymptotically stable.Subsystem 1 is defined as where α 1 is the virtual control law and can be designed as where k 1 is the designed control parameter and satisfy k 1 > 0.
Define the new error variable e 2 = x 2 − α 1 , the derivative of the Lyapunov function of the system can be calculated as The asymptotic stability of the controller requires that the Lyapunov function be negatively determined, so the new error variable must be zero, and the design of the controller subsystem 2.
Step 2: Define the new error variable e 2 and the corresponding Lyapunov function as On the premise of ignoring the disturbance force, the actual system will be established.In the subsequent nonlinear disturbance, the observer will give the disturbance on the UAAV: Thus, the actual control law is designed as where k 2 is the designed control parameter and satisfy k 2 > 0 and the dimensions of the column matrix are related to the system.Substituting the above-proposed controller into the Lyapunov function V 2 , yields .(22) where K = min{2k 1 , 2k 2 }.
To proceed, the Lyapunov function V 2 = 0.5 e 2 1 + e 2 2 of the actual system is satisfied to be negative definite, and the control law designed based on the backstepping method satisfies the asymptotic stabilization requirement of the closed-loop system.
In what follows, the nonlinear disturbance observer will be established to estimate the interference forces and moments near the surface.Designing a nonlinear disturbance observer to estimate the effect of near-surface disturbances on the control system as follows: p(x 1 , x 2 ) = c .
where d denotes the observation of the system disturbance by the nonlinear disturbance observer; f 2 (x 1 , x 2 ) is the nonlinear term; c is a large enough control parameter that can be adjusted by users.x 11 and x 22 satisfy x 11 = x y T and x 22 = z θ T .
Please see [28] for the convergence proof of the above observer.The controller is comprised of an inner loop for attitude control and an outer loop for position control.A nonlinear interferer is designed to estimate interference to the transmedia navigator, and the controller utilizes backstepping to achieve nonlinear control during the transmedia phase.Taking into account the coupling external handover interface disturbance, the final actual control law can be designed as follows: where k 1 and k 2 are the control parameters of the control law.Those are positive userdefined constants.e 1 and e 2 are the defining tracking error; .
x 1d is the derivative of the desired trajectory; d is the estimates of nonlinear interference observers; v 1 and v 2 are the compensated tracking error.The design of the controller and control law for the UAAV in the out-of-water stage has been completed, laying the foundation for the subsequent motion control simulation.

Remark 1.
For the designed control algorithm, the rule in parameter selections can be summarized as follows.According to the above theoretical analysis and the subsequent simulation experience, increasing k 1 and k 2 can improve the control accuracy and speed up the convergence, but also increase the control torque.Therefore, it is necessary to weigh the control performance, power consumption, and other indicators of the control system comprehensively.
Taking the derivative of V 2 and substituting the above-proposed controller into the Lyapunov function V 2 = 1 2 e T 1 e 1 + e T 2 e 2 , yields .
e 1 + e T
where K = min{2k 1 , 2k 2 − 1}, ∆ is the error resulting from the observer.Thus, it can be concluded that the error variable is stable.It is worth mentioning that the error converges to zero as long as KV 2 > ∆.

Numerical Simulation
For the UAAV described in ( 1) and ( 3), the trajectory tracking controller based on the disturbance observer adopts the same target trajectory to ensure that the vehicle can track the up-water trajectory under the disturbance of wind and wave loads, which is represented by the expected position and expected pitch angle.The actual combined wind and wave loads are designed as follows:

Motion Control Simulation
The main content of this section is whether the desired trajectory of the predetermined vehicle can be tracked under the wave load shown in (30).Figure 3 shows the three-dimensional sequence diagram of the vehicle's water emergence process under the disturbance of wind and wave loads.
Next, the expected trajectory is shown in (31) In the above vehicle exit-water trajectory, the initial position of the vehicle is the key factor for the exit.The initial state of the vehicle is selected as x d0 = 0, y d0 = 0, z d0 = −2 m, i.e., η 1d0 = [0 0 − 2] T ; The initial velocity is v(0) = [0.169, 0, 0.1333] T m/s, Initial trim angle is θ d0 = 0.698rad.
After many experiments in the simulation process, a set of parameters with better tracking effects is selected.Under environmental interference, the control parameters and the key parameters of the nonlinear function to be designed in the interference observer are k 1 = k 2 = 10 and c = 50.After the design of the above key parameters is completed, the In the above vehicle exit-water trajectory, the initial position of the vehicle is the key factor for the exit.The initial state of the vehicle is selected as , i.e., ; The initial velocity is , Initial trim angle is .
After many experiments in the simulation process, a set of parameters with better tracking effects is selected.Under environmental interference, the control parameters and the key parameters of the nonlinear function to be designed in the interference observer    In Figure 4, the position curve in the X-axis direction fluctuates and changes, mainly because the disturbed force in the emergence water stage is the superposition of disturbed force in the form of a sinusoidal curve, which is consistent with the real wave disturbance of the actual vehicle.In addition, the UAAV is in a state of large dip angle during the emergence process, which is greatly disturbed by wind and waves, resulting in the negative upward velocity along the X-axis and finally causing fluctuation.In Figure 4, the position curve in the X-axis direction fluctuates and changes, mainly because the disturbed force in the emergence water stage is the superposition of disturbed force in the form of a sinusoidal curve, which is consistent with the real wave disturbance of the actual vehicle.In addition, the UAAV is in a state of large dip angle during the emergence process, which is greatly disturbed by wind and waves, resulting in the negative upward velocity along the X-axis and finally causing fluctuation.The whole process of coming out of the water is divided into three stages: Stage 1 is when the vehicle is completely in the water, Stage 2 is when the vehicle is at the interface In Figure 4, the position curve in the X-axis direction fluctuates and changes, mainly because the disturbed force in the emergence water stage is the superposition of disturbed force in the form of a sinusoidal curve, which is consistent with the real wave disturbance of the actual vehicle.In addition, the UAAV is in a state of large dip angle during the emergence process, which is greatly disturbed by wind and waves, resulting in the negative upward velocity along the X-axis and finally causing fluctuation.
The whole process of coming out of the water is divided into three stages: Stage 1 is when the vehicle is completely in the water, Stage 2 is when the vehicle is at the interface between water and air, and Stage 3 is when the vehicle is in the air, as shown in Figure 5a.At around 0.9 s, the propeller of the vehicle begins to emerge from the surface of the water, and at around 1.8 s, the vehicle is completely detached from the water, ending in two different medium states of water and air.
The simulation results reflect the position and attitude tracking of the vehicle in the process of water emergence.It can be concluded that the vehicle is significantly affected by wind and wave loads in the water emergence stage under the interference of wind and wave loads, and the tracking effect is better when the vehicle is in the single medium stage.The position diagram of the vehicle in three-dimensional space is shown in Figure 5b, which shows the change in the space position of the vehicle in three stages.

Simulation of Interference Estimation
Based on Equation ( 22), an estimation model was established to accurately estimate the disturbance force and torque of the vehicle in the process of water emergence and track the preset expected water emergence trajectory.The motion control simulation experiment of the vehicle under wind and wave interference is carried out to verify the position and attitude control performance of the designed controller under environmental interference.The simulation results of interference force estimation caused by wind and wave load of the vehicle in the process of water discharge are obtained, as shown in Figure 6.The estimated simulation results of the disturbance torque rotating around the Y-axis in the inertial coordinate system are shown in Figure 7.
As shown in Figure 6a,c,e, the simulation results of wind and wave load interference force estimation can reflect that the estimation results obtained by the designed interference observer are compared with the actual input interference force and torque.This reflects that the design of the interference controller can rapidly reduce the estimation error within 0.05 seconds.Figure 6b,d,f can be calculated, showing that after stabilization, the maximum error of the estimated results of the interference observer and the actual input results is less than 3%, indicating that the accuracy of the estimation of the interference observer is reasonable.
periment of the vehicle under wind and wave interference is carried out to verify the position and attitude control performance of the designed controller under environmental interference.The simulation results of interference force estimation caused by wind and wave load of the vehicle in the process of water discharge are obtained, as shown in Figure 6.The estimated simulation results of the disturbance torque rotating around the Y-axis in the inertial coordinate system are shown in Figure 7.As shown in Figure 6a,c,e, the simulation results of wind and wave load interference force estimation can reflect that the estimation results obtained by the designed interference observer are compared with the actual input interference force and torque.This reflects that the design of the interference controller can rapidly reduce the estimation error within 0.05 seconds.Figure 6b,d,f can be calculated, showing that after stabilization, the maximum error of the estimated results of the interference observer and the actual input results is less than 3%, indicating that the accuracy of the estimation of the interference observer is reasonable.The simulation results of the disturbance torque estimation shown in Figure 7 can be analyzed and concluded.According to Figure 7a, it can be seen that the disturbance torque caused by wind and wave load is more accurate and responds faster than the disturbance force estimation.The error chart depicted in Figure 7b shows that compared with the force interference estimation, the initial error is smaller, and the time to reach the stable state of error is shorter.
The simulation results of the control signals shown in Figure 8 can be analyzed and concluded.According to Figure 8, it can be seen that the control signals are all bound and acceptable.Moreover, all the control signals have relatively large values at the beginning, and they converge to a small set with zero as the origin when the system is stable.
ence observer are compared with the actual input interference force and torque.This reflects that the design of the interference controller can rapidly reduce the estimation error within 0.05 seconds.Figure 6b,d,f can be calculated, showing that after stabilization, the maximum error of the estimated results of the interference observer and the actual input results is less than 3%, indicating that the accuracy of the estimation of the interference observer is reasonable.The simulation results of the disturbance torque estimation shown in Figure 7 can be analyzed and concluded.According to Figure 7a, it can be seen that the disturbance torque caused by wind and wave load is more accurate and responds faster than the disturbance force estimation.The error chart depicted in Figure 7b shows that compared with the force interference estimation, the initial error is smaller, and the time to reach the stable state of error is shorter.
The simulation results of the control signals shown in Figure 8 can be analyzed and concluded.According to Figure 8, it can be seen that the control signals are all bound and acceptable.Moreover, all the control signals have relatively large values at the beginning, and they converge to a small set with zero as the origin when the system is stable.In summary, building upon the vehicle exit strategy, this section comprehensively addresses the control of the vehicle during the exit phase under the interference of wind and waves.It takes into account the impact of wind and wave interference estimated by the interference observation state and researches vehicle trajectory tracking.

Conclusions
The key contribution of this paper lies in tracking the water trajectory of a cross-medium vehicle when subjected to disturbances from wind and wave loads and estimating the disturbance using a disturbance observer.The process involves establishing a model for the disturbance force and moment caused by wind and wave loads.Subsequently, a disturbance observer is designed, and the impact of disturbance is incorporated based on the established dynamics model.Building upon this, the paper conducts simulations for the water emergence control of the UAAV under disturbance conditions, verifying the effectiveness of the designed controller in trajectory tracking amid disturbances.Finally, the disturbance force and torque on the vehicle due to wind and waves are estimated and In summary, building upon the vehicle exit strategy, this section comprehensively addresses the control of the vehicle during the exit phase under the interference of wind and waves.It takes into account the impact of wind and wave interference estimated by the interference observation state and researches vehicle trajectory tracking.

Conclusions
The key contribution of this paper lies in tracking the water trajectory of a crossmedium vehicle when subjected to disturbances from wind and wave loads and estimating the disturbance using a disturbance observer.The process involves establishing a model for the disturbance force and moment caused by wind and wave loads.Subsequently,

Figure 1 .
Figure 1.Motion of UAAV in geodetic and body-fixed coordinate systems.

Figure 1 .
Figure 1.Motion of UAAV in geodetic and body-fixed coordinate systems.
of change in the Euler Angle.

Figure 2 .
Figure 2. Schematic diagram of the out-of-water phase.

Figure 2 .
Figure 2. Schematic diagram of the out-of-water phase.

Figure 3 .
Figure 3. Three-dimensional time-series diagram of the water emergence process under wind and wave disturbances for cross-medium vehicles.Next, the expected trajectory is shown in(31) After the design of the above key parameters is com- pleted, the motion control simulation under wind and wave load interference is carried out, and the simulation results are shown in Figures4 and 5 .

Figure 3 .
Figure 3. Three-dimensional time-series diagram of the water emergence process under wind and wave disturbances for cross-medium vehicles.

Figure 4 .
Figure 4. Time response of positional trajectory tracking for out-of-water processes under environmental disturbances in transmedia vehicles.(a) X-axis position tracking results.(b) Y-axis position tracking results.(c) z-axis position tracking results.

Figure 4 .
Figure 4. Time response of positional trajectory tracking for out-of-water processes under environmental disturbances in transmedia vehicles.(a) X-axis position tracking results.(b) Y-axis position tracking results.(c) z-axis position tracking results.

Figure 5 .
Figure 5. Schematic of attitude tracking time response and spatial 3D position of the cross-medium vehicle during water exit.(a) Simulation results of attitude tracking during the water-exit process.(b) Schematic diagram of the three-dimensional position tracking.

Figure 5 .
Figure 5. Schematic of attitude tracking time response and spatial 3D position of the cross-medium vehicle during water exit.(a) Simulation results of attitude tracking during the water-exit process.(b) Schematic diagram of the three-dimensional position tracking.

Figure 6 .
Figure 6.Simulation results of interference force estimation for wind and wave loaded interference observer during the water emergence phase.(a) Estimation of wind and wave load disturbances in the X-axis.(b) Errors of disturbances estimation in the X-axis.(c) Estimation of wind and wave load disturbances in the Y-axis.(d) Errors of disturbances estimation in the Y-axis.(e) Estimation of wind and wave load disturbances in the z-axis.(f) Errors of disturbances estimation in the z-axis.

Figure 6 .
Figure 6.Simulation results of interference force estimation for wind and wave loaded interference observer during the water emergence phase.(a) Estimation of wind and wave load disturbances in the X-axis.(b) Errors of disturbances estimation in the X-axis.(c) Estimation of wind and wave load disturbances in the Y-axis.(d) Errors of disturbances estimation in the Y-axis.(e) Estimation of wind and wave load disturbances in the z-axis.(f) Errors of disturbances estimation in the z-axis.

Figure 7 .
Figure 7. Simulation results of interference moment estimation for wind and wave loaded interference observer during the out-of-water phase.(a) Estimation of wind and wave disturbance moments around the Y-axis direction.(b) Error in estimation around Y-axis direction.

Figure 7 .Figure 8 .
Figure 7. Simulation results of interference moment estimation for wind and wave loaded interference observer during the out-of-water phase.(a) Estimation of wind and wave disturbance moments around the Y-axis direction.(b) Error in estimation around Y-axis direction.J. Mar.Sci.Eng.2023, 11, x FOR PEER REVIEW 15 of 17

Figure 8 .
Figure 8. Simulation results of control signals.(a) Control signal for x channel.(b) Control signal for y channel.(c) Control signal for z channel.(d) Control signal for θ channel.