Next Article in Journal
Pursuing Optimization Using Multimodal Transportation System: A Strategic Approach to Minimizing Costs and CO2 Emissions
Next Article in Special Issue
Applications of Voronoi Diagrams in Multi-Robot Coverage: A Review
Previous Article in Journal
Development of Hydroacoustic Localization Algorithms for AUV Based on the Error-Corrected WMChan-Taylor Algorithm
Previous Article in Special Issue
Integral Sliding Mode Output Feedback Control for Unmanned Marine Vehicles Using T–S Fuzzy Model with Unknown Premise Variables and Actuator Faults
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Distributed Heterogeneous Formation Control for UAV-USVs with Input Quantization

1
College of Navigation, Dalian Maritime University, Dalian 116026, China
2
School of Instrumentation Science and Engineering, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2024, 12(6), 975; https://doi.org/10.3390/jmse12060975
Submission received: 25 April 2024 / Revised: 1 June 2024 / Accepted: 7 June 2024 / Published: 11 June 2024
(This article belongs to the Special Issue Unmanned Marine Vehicles: Navigation, Control and Sensing)

Abstract

:
This paper investigates the cooperative formation trajectory tracking problem for heterogeneous unmanned aerial vehicle (UAV) and multiple unmanned surface vessel (USV) systems with input quantization performance. Firstly, at the kinematic level, a distributed guidance law based on an extended state observer (ESO) is designed to compensate for the unknown speed of neighbor agents for expected trajectory tracking, and subsequently at the dynamic level, an ESO is utilized to estimate model uncertainties and environmental disturbances. Following that, a linear analytic model is employed to depict the input quantization process, and the corresponding adaptive quantization controller is designed without necessitating prior information on quantization parameters. Based on the input-to-state stability, the stability of the proposed control structure is proved, and all the signals in the closed-loop system are ultimately bounded. Finally, a simulation study is provided to show the efficacy of the proposed strategy.

1. Introduction

In recent times, owing to escalating societal demands, there has been a significant advancement in the development of unmanned systems [1,2,3]. These systems offer benefits such as compactness, agility, affordability, and the elimination of casualty risks [4,5,6]. They can operate in hazardous zones or areas unsuitable for manned vehicles, broaden operational possibilities, and are cost effective [7,8]. To further enhance the operational efficiency of these unmanned systems, the focus on formation strategies in unmanned systems has gained prominence in practical engineering applications [9,10].
The application of unmanned surface vehicles (USVs) in ocean engineering has expanded significantly, which can be attributed to their advantages such as an enhanced load capacity, greater operational convenience, and reduced cost of missions [11]. Nevertheless, the limited observational range of USVs poses challenges in locating targets during maritime search and rescue operations [12]. To address this limitation, integrating unmanned aerial vehicles (UAVs) with USV systems to form heterogeneous multi-agent systems has proven beneficial [13]. A USV supplemented by the agility of a UAV and its superior target search capabilities effectively broadens the communication scope at sea among USVs, thereby expanding the operational area of these heterogeneous systems [14,15]. Therefore, exploring the formation dynamics within these heterogeneous systems holds substantial importance in engineering.
Significant advancements in the cooperative tracking control of UAVs and USVs have been made, with notable methodologies including the leader–follower method [16], virtual structure method [17], behavior-based approach [18], and model prediction method [19]. Kai Xue [20] proposes a distributed collaborative tracking method based on leader–follower for a heterogeneous multi-agent system composed of UAVs and multiple USVs. Haitao Liu [21] proposes a distributed formation control protocol for UAV-USV heterogeneous systems to ensure that the target trajectory is tracked within a fixed time and a predetermined formation configuration is formed under various uncertainties. Shaoshi Li [22] proposes a formation control scheme based on velocity estimation, which only requires the relative azimuth vector between adjacent aircraft to achieve the desired formation. Dapeng Huang [19] proposes a distributed model predictive formation control algorithm for heterogeneous multi-agent systems based on directed topology, which realizes the coordination of heterogeneous formation between multiple UAVs and USVs. In literature [23], a global fixed-time adaptive neural network non-singular fast terminal sliding formation controller is designed to ensure that the trajectory of the target is tracked within a fixed time and a predetermined formation configuration is formed. Huiyan Li [24] designed a linear time-varying UAV-USV system based on distributed control, proposed the necessary conditions to maintain the consistency of the system, and expounded the effects of inertia matrix parameters and the heterogeneity between UAVs and USVs on the controller performance. The collective research efforts of these scholars have paved the way for enhanced understanding and practical implementation of heterogeneous multi-agent formation and control within the realm of UAVs and USVs.
Nonetheless, it is worth noting that the aforementioned scholarly works have yet to address the intricate challenge posed by the constrained maritime communication bandwidth. For navigational tasks, the transfer of data among sensor components is dependent on the use of communication channels. Consequently, this requires the information to be quantized and encoded before its transmission. Moreover, in light of the channel bandwidth constraints, maintaining a seamless system functioning within the defined bandwidth limits necessitates incorporating quantization technologies. The application of these quantization methods enables a strategic decrease in communication frequency and effectively manages the number of controllers within a specific bandwidth range.
At present, some studies have considered the quantitative problems existing in the practice of physical systems [25,26,27,28,29,30,31,32,33]. However, the input quantization of heterogeneous UAV and USV formation control systems has not been considered in the literature.
Based on the above problems, this paper designs the adaptive distributed heterogeneous formation control algorithm with input quantization. The research aims to advance cooperative tracking control strategies for UAV and USV formations, focusing on input quantization. This paper proposes an adaptive distributed heterogeneous formation control algorithm that incorporates input quantization to manage the constraints posed by limited communication bandwidth in maritime environments.
In contrast to the existing body of research, this paper stands out for the following salient attributes:
(1) Novelty in problem scope: Unlike prior studies [14,19,20,21,22,23,24,34,35,36,37,38], these research pioneers the exploration of the intricate domain concerning heterogeneous UAV and USV formation control coupled with input quantization. The integration of these two crucial factors is unprecedented in the field.
(2) Advanced estimation technique: Distinguishing itself from works such as [17,20,21,39,40,41], this paper leverages an extended state observer (ESO) methodology to assess the influence of external environmental factors deftly. Furthermore, it presents a strategic resolution to the challenges posed by heterogeneity by employing a unified model approach.
(3) Enhanced quantization methodology: In comparison to studies like [5,6,8,15,25,26,27,30,31,32,33], this paper introduces a unique approach to the quantization process. It elucidates this process via a linear time-varying model, which offers notable advantages. Most notably, it eliminates the need for explicit quantization parameter information, thus streamlining the implementation process.
The organization of this paper is as follows. Section 2 introduces some necessary preliminaries and states the problem formulation. Section 3 expounds on the controller design. Section 4 provides the closed-loop system stability analysis. Section 5 demonstrates the simulation studies. Section 6 concludes the paper.

2. Problem Formulation

Consider a heterogeneous multi-agent system (HMAS) consisting of one leader UAV and N follower USVs. According to the results, in [42] the dynamic model of a quadrotor UAV can be described as
p ¨ a x = cos ϕ sin θ cos ψ + sin ϕ sin ψ u p m a d x p ˙ a x m a + Δ a x , p ¨ a y = cos ϕ sin θ sin ψ + sin ϕ cos ψ u p m a d y p ˙ a y m a + Δ a y , p ¨ a z = cos θ cos ϕ u p m a d z p ˙ a z m a g + Δ a z ,
ϕ ¨ = θ ˙ ψ ˙ J a y J a z J a x J a r J a x θ ˙ d ¯ + τ ϕ J a x d ϕ ϕ ˙ J a x , θ ¨ = ϕ ˙ ψ ˙ J a z J a x J a y J a r J a y ϕ ˙ d ¯ + τ θ J a y d θ θ ˙ J a y , ψ ¨ = ϕ ˙ θ ˙ J a x J a y J a z + τ ψ J a z d ψ ψ ˙ J a z ,
where p a x , p a y , p a z T indicates the spatial location information of the UAV; ϕ , θ , ψ T indicates the attitude information of the UAV; u p is defined as the control thrust of the UAV; τ ϕ , τ θ , τ ψ T indicates the three control torques of the UAV; m a is the mass of the UAV; g indicates the gravitational acceleration; d ¯ indicates the overall residual rotor angle; d x , d y , d z , d ϕ , d θ , d ψ denote the translational drag coefficients; J a x , J a y , J a z are the torques of the inertia; J a r is the torque of the rotor’s inertia; and Δ a x , Δ a y , Δ a z indicate the external interference encountered by the UAV.
Inspired by these results, considering external disturbances and parametric uncertainties, the three-degrees-of-freedom model of a UAV (1) can be redefined as follows:
p ¨ a = g a u a + f a + Δ a .
where p a = p a x , p a y , p a z T denotes the position of the UAV, f a = d x p ˙ a x / m a , d y p ˙ a y / m a , d z p ˙ a z / m a g T , g a = diag 1 / m a , 1 / m a , 1 / m a , Δ a = Δ a x , Δ a y , Δ a z T , and u a = u a x , u a y , u a z T is defined as
u a x = cos ϕ sin θ cos ψ + sin ϕ sin ψ u p , u a y = cos ϕ sin θ sin ψ + sin ϕ cos ψ u p , u a z = cos θ cos ϕ u p .
The kinematic and dynamic equation of the i-th USV in the XY plane is defined as [43]
x ˙ s i = μ s i cos ψ s i v s i sin ψ s i , y ˙ s i = μ s i sin ψ s i + v s i cos ψ s i , ψ ˙ s i = r s i .
μ ˙ s i = f μ s i α i + 1 m μ s i τ μ s i f + w μ s i , v ˙ s i = f v s i α i + 1 m v s i w v s i , r ˙ s i = f r s i α i + 1 m r s i τ r s i f + w r s i .
and
f μ s i ( α i ) = 1 m μ s i ( m v s i v s i r s i d μ s i μ s i d μ s i 1 | μ s i | μ s i ) , f v s i ( α i ) = 1 m μ s i ( m μ s i μ s i r s i d v s i v s i d v s i 1 | v s i | v s i ) , f r s i ( α i ) = 1 m r s i ( ( m μ s i m v s i ) μ s i v s i d r s i r s i d r s i 1 | r s i | r s i ) ,
where x s i , y s i is defined as the position information of the i-th USV; ψ s i is defined as the yaw angle of the i-th USV; α i = μ s i , v s i , r s i T are the surge, sway, and yaw velocity, respectively; m μ s i , m v s i , m r s i are the inertial mass; f μ s i α i , f v s i α i , f r s i α i are the nonlinear unknown functions consisting of the unmodeled hydrodynamics and Coriolis forces; and τ μ s i f and τ r s i f are the surge force and the yaw moment. w μ s i , w v s i , w r s i are the bounded disturbances caused by waves, wind, and ocean currents.
The hand position method is a control technique employed to manage the movement and positioning of USVs or similar underactuated systems. Underactuated systems have fewer control inputs than degrees of freedom, making traditional control strategies less effective. The hand position method addresses this challenge by defining a virtual control point (hand point) on the vehicle, allowing for more effective control. To deal with the underactuated USV motion model described by (4) and (5), the idea of the hand position method is considered. Define the front point p s i x , p s i y of the USVs as the hand point, which can be formulated as
p s i x = x s i + L s i cos ψ s i , p s i y = y s i + L s i sin ψ s i .
where L s i is the distance between the actual position x s i , y s i and the new defined hand point p s i x , p s i y , which is shown in Figure 1.
Taking the second derivative of (7) shows that
p ¨ s i x = μ ˙ s i cos ψ s i v ˙ s i + L s i r ˙ s i sin ψ s i μ s i r s i sin ψ s i v s i r s i + L s i r s i 2 cos ψ s i , p ¨ s i y = μ ˙ s i sin ψ s i v ˙ s i + L s i r ˙ s i cos ψ s i + μ s i r s i cos ψ s i v s i r s i + L s i r s i 2 sin ψ s i ,
substituting (6) into (8) yields that
p ¨ s i x = f s i x ( β ) + cos ψ s i m μ s i τ μ f L s i sin ψ s i m r s i τ r f + w d i x , p ¨ s i y = f s i y ( β ) + sin ψ s i m μ s i τ μ f + L s i cos ψ s i m r s i τ r f + w d i y .
where
f s i x ( β ) = f μ ( α ) cos ψ s i f v ( α ) + L s i f r ( α ) sin ψ s i μ s i r s i sin ψ s i v s i r s i + L s i r s i 2 cos ψ s i , f s i y ( β ) = f μ ( α ) sin ψ s i + f v ( α ) + L s i f r ( α ) cos ψ s i + μ s i r s i cos ψ s i v s i r s i + L s i r s i 2 sin ψ s i ,
w d i x ( β ) = w μ s i m μ s i cos ψ s i w v s i m μ s i + L s i w r s i m r s i sin ψ s i , w d i y ( β ) = w μ s i m μ s i sin ψ s i + w v s i m μ s i + L s i w r s i m r s i cos ψ s i ,
with β = μ s i , v s i , r s i , ψ s i .
Based on (9), the dynamics model of the i-th USV can be described as
p ¨ s i = f s i x y + Ω s i ψ s i ω s i u s i + w d i x y .
where p s i = p s i x , p s i y T is the position of the i-th i Π 2 USV, and f s i x y = f s i x , f s i y T , Ω s i ψ s i = cos ψ s i , sin ψ s i ; sin ψ s i , cos ψ s i , u s i = τ μ , τ r T , ω s i = diag 1 / m μ s i , L s i / m r s i , w d i x y = w d i x , w d i y T .
To ensure the unity of the USV model and the UAV model, this formula is expanded and further converted to
p ¨ s i x p ¨ s i y = f s i x L s i sin ψ s i m r s i τ r i f s i y + sin ψ s i m μ s i τ μ i + cos ψ s i m μ s i 0 0 L s i cos ψ s i m r s i τ μ i τ r i + w d i x w d i y .
Combining (3) and (11), the unified model of the HMAS in the XY plane can be obtained:
x ˙ i 1 = x i 2 x ˙ i 2 = F x i + G x i u x i + Δ x i .
When Equation (12) represents the UAV model. Where x i 1 = x a 1 R 2 , x i 2 = x a 2 R 2 is defined as the position and speed information of the UAV in the X Y plane, F x i = f a x y is defined as unknown functions, G x i = g a x y is defined as coefficient matrix, Δ x i = Δ a x y indicates the external interference encountered by the UAV, u x i = u a x y is defined as coefficient matrix control input signal. When Equation (12) represents the USV model. where x i 1 = x s i 1 R 2 , x i 2 = x s i 2 R 2 are the position and velocity of the i-th USV, F x i = F s i = f s i x L s i sin ψ s i m r s i τ r i , f s i y + sin ψ s i m μ s i τ μ i T is defined as unknown functions, G x i = G s i is defined as a coefficient matrix, Δ x i = w d i x y is the bounded disturbances, u x i = u s i is defined as a coefficient matrix control input signal.
Similarly, consider the UAV model in the Z-axis direction as
p ˙ a z = v a z , v ˙ a z = f a z + 1 m a u a z + Δ a z = F a z + G a z u a z + Δ a z .
where p a z is the altitude of the UAV; v a z is defined as the speed of the UAV in the Z-axis, F a z = d z p ˙ a z d z p ˙ a z m a m a g , and G a z = 1 m a .
The quantification technique can convert the initial continuous signal into a signal that is segmented and constant; the HMAS model with input quantization is represented as
U ˙ s i V ˙ s i = f s i x L s i sin ψ s i m r s i Q ( τ r i ) f s i y + sin ψ s i m μ s i Q ( τ μ i ) + cos ψ s i m μ s i 0 0 L s i cos ψ s i m r s i Q ( τ μ i ) Q ( τ r i ) + w d i x w d i y .
In a similar way,
U ˙ a V ˙ a = d x p ˙ a x m a d y p ˙ a y m a + m a 1 0 0 m a 1 Q ( u a x ) Q ( u a y ) + Δ a x Δ a y .
Combining (14) and (15), the HMAS model with input quantization is represented as
x ˙ i 1 = x i 2 , x ˙ i 2 = F x i + G x i Q ( u x i ) + Δ x i .
Based on (16), the unified dynamic model of an HMAS in the XY plane can be obtained as
U i V i = F i + G i Q τ i + Δ i ,
where F i = F i U , F i V T , Q τ i = Q ( τ i U ) , Q ( τ i V ) T , G i = d i a g [ G i U , G i V ] , Δ i = Δ i U , Δ i V T , and τ i = u x i .
The uniform quantization has the characteristic of easy realization; it can be expressed as follows:
Q ( τ i U ) = k i U r o u n d ( τ i U k i U ) ,
Q ( τ i V ) = k i V r o u n d ( τ i V k i V ) ,
where k i U and k i V denote the quantized indexes. The round function can return a value rounded to a specified number of decimal places. During the quantization process, quantization errors will be generated, and the control system’s performance could be reduced. A continuous and differentiable time-varying reference trajectory is defined as
p 0 ( t ) = [ x 0 ( t ) , y 0 ( t ) ] T R 2 ,
G = { V , E } represents the USV and UAV communication topology diagram, Where V = { n 0 , n 1 , , n N } is the set of all nodes combined, denoting N agents in formation; E = { ( i , j ) V × V } is a set of all edges, and A = a i j R N × N is the adjacency matrix, when ( n i , n j ) E , a i j = 1 ; if not, a i j = 0 .
Assumption 1.
The aerodynamic drag coefficients d i x , d i y , and d i z are bounded and unknown. The quadrotor UAV experiences external perturbations Δ x , Δ y , Δ z , which are confined within certain bounds, fulfilling the conditions Δ a x Δ ¯ a x , Δ a y Δ ¯ a y , Δ a z Δ ¯ a z . Here, Δ ¯ a x , Δ ¯ a y , and Δ ¯ a z represent unknown positive constants. The USV experiences external disturbances w μ s i , w ν s i , w r s i , which are confined within certain limits, fulfilling the conditions w μ s i w ¯ μ s i , w v s i w ¯ v s i , w r s i w ¯ r s i . In this context, w ¯ μ s i , w ¯ v s i , and w ¯ r s i signify unknown positive constants.
Remark 1.
The complexity of operating environments introduces unpredictable aerodynamic forces, complicating the acquisition of precise system parameters. Moreover, considering that UAVs and USVs are not operated under extreme weather conditions, it renders Assumption 1 both justified and pragmatic.
In this paper, a distributed controller is considered for each unmanned agent to ensure that each USV and UAV can track the expected time-varying trajectory and ensure
lim t p i ( t ) p i d ( t ) p 0 ( t ) l i 0 ,
where l i 0 is a positive constant, and p i d ( t ) R 2 is the position deviation relative to the parameterized path. p i ( t ) = [ x i ( t ) , y i ( t ) ] T R 2 .

3. Controller Design

This section presents a formation control strategy that employs the ESO, designed to direct the HMAS to adhere to a specified trajectory within the XY plane. Following this, a decentralized formation control mechanism is formulated, focusing on the altitude tracking of the UAV in the Z-axis. The formation control flowchart of the HMAS is shown in Figure 2.

3.1. Kinematic Controller Design

The formation error can be defined as
z i = j = 1 N a i j ( p i p j p i j d ) + a i 0 ( p i p 0 p i d ) ,
where p i j d = p i d p j d R 2 is the relative position deviation between HMAS.
Then, the time derivative of (20) can be obtained as follows:
z ˙ i = j = 1 N a i j { U i , V i T [ U j , V j ] T p ˙ i j d } + a i 0 U i , V i T a i 0 p ˙ 0 a i 0 p ˙ i d ,
which can be rewritten as
z ˙ i = d i U i , V i T a i 0 p ˙ 0 + σ i ,
where d i = j = 0 N a i j , and σ i is defined as σ i = j = 1 N a i j [ U j , V j ] T j = 0 N a i j p i j d .
Assumption 2.
Since the velocity and acceleration of HMAS are both upper bounds, the unknown term σ i is assumed to satisfy σ ˙ i σ * , where σ * is a positive constant. The extended state observer can be used to observe the unknown term σ i .
The ESO is designed as follows:
z ^ ˙ i = ϑ i + σ ^ i k i 1 z ˜ i , σ ^ ˙ i = k i 2 z ˜ i .
where z ˜ i = z ^ i z i is defined as the ESO estimation error and ϑ i = d i U i , V i T a i 0 p ˙ 0 . Based on (22) and (23), it can be found that
z ˜ ˙ i = σ ˜ i k i 1 z ˜ i , σ ˜ ˙ i = k i 2 z ˜ i σ ˙ i .
Equation (24) can be rewritten as
E ˜ ˙ i 1 = A i E ˜ i 1 + B i ,
where E ˜ i 1 = [ z ˜ i , σ ˜ i ] T , and
A i = k i 1 1 k i 2 0 , B i = 0 2 σ ˙ i .
At the kinematic level, a distributed control law considering ESO is designed:
U i c V i c = d 1 [ χ i z i + a i 0 p ˙ 0 σ ^ i ] ,
where χ i = d i a g [ k i x , k i y ] R 2 is positive constant.

3.2. Dynamic Controller Design

The control objective of the dynamic subsystem is
lim t e i U = U i U i d δ 1 , lim t e i V = V i V i d δ 2 ,
where δ 1 , δ 2 is positive constant, Let Q ( τ i U ) = q 1 i U ( t ) τ i U + q 2 i U ( t ) , Q ( τ i V ) = q 1 i V ( t ) τ i V + q 2 i V ( t ) , and
q 1 i U ( t ) = Q τ i U ( t ) τ i U ( t ) τ i U ( t ) a 1 τ i U ( t ) < a , q 1 i V ( t ) = Q τ i V ( t ) τ i V ( t ) τ i V ( t ) a 1 τ i V ( t ) < a , q 2 i U ( t ) = 0 τ i U ( t ) a Q τ i U ( t ) τ i U ( t ) τ i U ( t ) < a , q 2 i V ( t ) = 0 τ i V ( t ) a Q τ i V ( t ) τ i V ( t ) τ i V ( t ) < a ,
where q 1 i U ( t ) , q 1 i V ( t ) are unknown. Since the positive and negative values of the signal are not changed by the quantization process, q 1 i U ( t ) , q 1 i V ( t ) are obtained. If τ i U ( t ) < a , τ i V ( t ) < a , then Q ( τ i U ( t ) ) , Q ( τ i V ( t ) ) are bounded since q 2 i U ( t ) , q 2 i v ( t ) are bounded. q 2 i U ( t ) q ¯ 2 , q 2 i V ( t ) q ¯ 2 .
Consider the following two sliding mode functions:
s i U = c i U 0 t e i U d t + e i U ( t ) , s i V = c i V 0 t e i V d t + e i V ( t ) ,
where s i = [ s i U , s i V ] T , c i U > 0 and c i V > 0 . The time derivative of (30) shows that
s ˙ i U = c i U e i U + e ˙ i U = G i U q 1 i U ( t ) τ i U + G i U q 2 i U ( t ) + Δ i U + F i U U ˙ i d + c i U e i U .
s ˙ i V = c i V e i V + e ˙ i V = G i V q 1 i V ( t ) τ i V + G i V q 2 i V ( t ) + Δ i V + F i V V ˙ i d + c i V e i V .
Let d i U ( X ) = Δ i U + F i U , d i V ( X ) = Δ i V + F i V , it follows that
s i U s ˙ i U s i U l i U s i U η i U sgn s i U + l i U s i U + η i U sgn s i U + 1 2 G i U s i U + d i U ( X ) + c i U e i U U ˙ i d + s i U G i U q 1 i U τ i U + 1 2 G i U q ¯ 2 i U 2 .
s i V s ˙ i V s i V l i V s i V η i V sgn s i V + l i V s i V + η i V sgn s i V + 1 2 G i V s i V + d i V ( X ) + c i V e i V V ˙ i d + s i V G i V q 1 i V τ i V + 1 2 G i V q ¯ 2 i V 2 .
where the ‘sgn(x)’ function is a mathematical function that extracts the sign of a real number. It is defined as follows: If the input x is negative, the output of ‘sgn(x)’ function is −1. If the input x is zero, the output of ‘sgn(x)’ function is 0. If the input x is positive, the output of ‘sgn(x)’ function is 1.
Since the specific function value of d i ( X ) = [ d i U , d i V ] T is unknown, the extended state observer can be used to observe it:
p ^ ¨ i = G i Q τ i U , Q τ i V T + d ^ i ( X ) b i 1 p ˜ ˙ i , d ^ ˙ i ( X ) = b i 2 p ˜ ˙ i ,
where p ˜ ˙ i = p ^ ˙ i p ˙ i is the ESO estimation error; we have
p ˜ ¨ i = d ˜ i ( X ) b i 1 p ˜ ˙ i , d ˜ ˙ i ( X ) = b i 2 p ˜ ˙ i d ˙ i ( X ) ,
Which can be reshaped into
E ˜ ˙ i 2 = C i E ˜ i 2 + D i ,
where E ˜ i 2 = p ˜ ˙ i , d ˜ i ( X ) T , and
C i = b i 1 1 b i 2 0 , D i = 0 2 d ˙ i ( X ) .
Define
κ i U = l i U s i U + η i U sgn s i U + 1 2 G i U s i U + d ^ i U ( X ) + c i U e i U U ˙ i d ,
where l i U and η i U are positive constants. It follows that
s i U s ˙ i U l i U s i U 2 η i U s i U + s i U κ i U + s i U d ˜ i U ( X ) + s i U G i U q 1 i U τ i U + 1 2 G i U q ¯ 2 i U 2 ,
where d ˜ i U ( X ) = d i U ( X ) d ^ i U ( X ) .
Similarly, let κ i V = l i V s i V + η i V sgn s i V + 1 2 G i V s i V + d ^ i V ( X ) + c i V e i V V ˙ i d ; it follows that
s i V s ˙ i V l i V s i V 2 η i V s i V + s i V κ i V + s i V d ˜ i V ( X ) + s i V G i V q 1 i V τ i V + 1 2 G i V q ¯ 2 i V 2 .
Remark 2.
Because q 1 i U ( t ) and q 1 i V ( t ) are time-varying and unknown, adaptive estimation is adopted. The lower bounds of q 1 i U ( t ) and q 1 i V ( t ) are used for estimation to prevent singular problems that arise when the estimate is zero. Time-varying gains μ i U and μ i V are defined as 1 / q 1 i U min and 1 / q 1 i V min , respectively, where q 1 i U min and q 1 i V min represent the lower bounds of q 1 i U ( t ) and q 1 i V ( t ) .
The dynamic controllers and adaptive laws are designed as follows:
τ i U = G i U 1 s i U μ ^ i U 2 κ i U 2 s i U μ ^ i U κ i U + ρ i U , τ i V = G i V 1 s i V μ ^ i V 2 κ i V 2 s i V μ ^ i V κ i V + ρ i V ,
μ ^ ˙ i U = γ 1 s i U κ i U γ 1 ς i U μ ^ i U , μ ^ ˙ i V = γ 2 s i V κ i V γ 2 ς i V μ ^ i V ,
where ρ i U , ρ i V , γ 1 , γ 2 , ς i U , ς i V are positive constants.
The error subsystem generated by z i , s i U , s i V , μ ˜ i U , μ ˜ i V is as follows:
z ˙ i = χ i z i σ ˜ i , s ˙ i U = c i U e i U + e ˙ i U , s ˙ i V = c i V e i V + e ˙ i V , μ ˜ ˙ i U = μ ^ ˙ i U μ ˙ i U = γ 1 s i U κ i U γ 1 ς i U μ ^ i U , μ ˜ ˙ i V = μ ^ ˙ i V μ ˙ i V = γ 2 s i V κ i V γ 2 ς i V μ ^ i V .

3.3. Altitude Controller Design of UAV

The height error system is defined as
e a z p = p a z p a z d , e a z v = v a z v a z d ,
where p a z d is defined as the desired position information, and v a z d is defined as the desired velocity information.
Consider the following input and adaptive laws:
u a z = G ^ a z 1 σ z e ζ e ζ κ ^ 11 h T h 2 ς 1 2 e ζ κ ^ 12 2 ς 2 2 + v ˙ a z d k ζ e a z v , κ ^ ˙ 11 = l 11 k 11 κ ^ 11 + e ζ T e ζ h T h 2 ζ 1 2 , κ ^ ˙ 12 = l 12 k 12 κ ^ 12 + e ζ T e ζ 2 ς 2 2 , G ^ ˙ a z = Pr o j G ^ a z , G ¯ a z { S } = 0 if G ^ a z = G ¯ a z , S 0 or G ^ a z = G ^ a z , S 0 , S otherwise ,
where S = ι 13 ( k 13 G ^ a z + e ζ u a z ) , and G ¯ a z , G a z is the lower and maximum bound of G ^ a z . k11, k12, k13, σz, l11, l12, l13, kζ, ζ1, ζ2 are positive parameters to be designed.
Remark 3.
An adaptive law (43) is used to estimate quantization parameters to ensure that the designed dynamic control law (42) does not depend on quantization information. Quantization parameters can be adjusted according to the needs of system performance.

4. Stability Analysis

Theorem 1.
System (25) consists of the said input σ ˙ i and states z ˜ i , σ ˜ i . The system is input-state stable, so that for every bounded input over an infinite time horizon, the state of the system also remains bounded.
Proof. 
The following Lyapunov function is designed:
V ˙ 1 = i = 1 N E ˜ i 1 T P i E ˜ i 1 ,
where P i represents a positive definite matrix which complies with A i T P i + P i T A i I . Then, the time derivative of (25) can be obtained as follows:
V ˙ 1 i = 1 N { E i 1 T E i 1 + E i 1 T P i B i σ ˙ i } .
Due to
E i 1 P i B i σ ˙ i θ ¯ i 1 ,
we obtain
V ˙ 1 i = 1 N ( 1 θ ¯ i 1 ) E i 1 2 ,
where 0 < θ ¯ i 1 < 1 . For input σ ˙ i , the subsystem (25) is stable from input to state, and
E i 1 ( t ) λ max P i λ min P i max E i 1 t 0 e ϱ i 1 t t 0 , P i B i σ * θ ¯ i 1 , t t 0 ,
where ϱ i 1 = 1 θ ¯ i 1 / λ max P i .
Theorem 2.
System (37) consists of the said input d ˙ i ( X ) and states p ˜ ˙ i , d ˜ i ( X ) . The system is input-state stable so that for every bounded input over an infinite time horizon, the state of the system also remains bounded.
Proof. 
The following Lyapunov function is designed:
V ˙ 2 = i = 1 N E ˜ i 2 T Q i E ˜ i 2 ,
where Q i represents a positive definite matrix which complies with
C i T Q i + Q i T C i I .
Then, the time derivative of (37) can be obtained as follows:
V ˙ 2 i = 1 N E i 2 T E i 2 + E i 2 T Q i D i d ˙ i ( X ) ,
since
E i 2 Q i D i d ˙ i ( X ) θ ¯ i 2 ,
we obtain V ˙ 2 i = 1 N ( 1 θ ¯ i 2 ) E i 2 2 , where 0 < θ ¯ i 2 < 1 . For input d ˙ i ( X ) , the subsystem (37) is stable from input to state, and
E i 2 ( t ) λ max Q i λ min Q i max E i 2 t 0 e ϱ i 2 t t 0 , Q i D i d ( X ) * θ ¯ i 2 , t t 0 ,
where ϱ i 2 = 1 θ ¯ i 2 / λ max Q i . □
Theorem 3.
The dynamic error subsystem (44), viewed as the system with the states being z i , s i U , s i V , μ ˜ i U , μ ˜ i V and the inputs being q 2 i U , q 2 i V , μ i U , μ i V , σ ˜ i , is stable from input to state.
Proof. 
Define the following Lyapunov function:
V 3 = i = 1 N 1 2 s i U 2 + 1 2 γ 1 μ i U μ ˜ i U 2 + 1 2 s i V 2 + 1 2 γ 2 μ i V μ ˜ i V 2 + 1 2 z i T z i
where μ i U > 0 and μ i V > 0 .
Then, the time derivative of (57) can be obtained as follows:
V ˙ 3 = i = 1 N s i U s ˙ i U + 1 γ 1 μ i U μ ˜ i U μ ^ ˙ i U + s i V s ˙ i V + 1 γ 2 μ i V μ ˜ i V μ ^ ˙ i V + z i T z ˙ i .
Substituting Equations (40)–(43) into (58), it can be seen that
V ˙ 3 i = 1 N l i U 1 2 s i U 2 η i U s i U + s i U κ i U + 1 2 G i U q ¯ 2 i U 2 q 1 i U s i U 2 μ ^ i U 2 κ i U 2 s i U μ ^ i U κ i U + ρ i U + 1 γ 1 μ i U μ ˜ i U γ 1 s i U κ i U γ 1 ς i U μ ^ i U l i V 1 2 s i V 2 η i V s i V + s i V κ i V + 1 2 G i V q ¯ 2 i V 2 q 1 i V s i V 2 μ ^ i V 2 κ i V 2 s i V μ ^ i V κ i V + ρ i V + 1 γ 2 μ i V μ ˜ i V γ 2 s i V κ i V γ 2 ς i V μ ^ i V + z i T z ˙ i + 1 2 d ˜ i U 2 ( X ) + 1 2 d ˜ i V 2 ( X ) .
From a a 2 ρ + a < ρ and considering q 1 i U q 1 i U min = 1 / μ i U > 0 and q 1 i V q 1 i V min = 1 / μ i V > 0 , we have
V ˙ 3 i = 1 N l i U 1 2 s i U 2 η i U s i U + 1 2 G i U q ¯ 2 i U 2 + 1 μ i U ρ i U 1 μ i U μ ˜ i U ς i U μ ^ i U l i V 1 2 s i V 2 η i V s i V + 1 2 G i V q ¯ 2 i V 2 + 1 μ i V ρ i V 1 μ i V μ ˜ i V ς i V μ ^ i V + z i T z ˙ i + 1 2 d ˜ i U 2 ( X ) + 1 2 d ˜ i V 2 ( X ) .
Noting that
μ ˜ i U μ ^ i U 1 2 μ ˜ i U 2 + 1 2 μ i U 2 , μ ˜ i V μ ^ i V 1 2 μ ˜ i V 2 + 1 2 μ i V 2 ,
it follows that (60) can be put into
V ˙ 3 i = 1 N l i U 1 2 s i U 2 η i U s i U + 1 2 G i U q ¯ 2 i U 2 + 1 μ i U ρ i U 1 2 μ i U ς i U μ ˜ i U 2 + 1 2 ς i U μ i U l i V 1 2 s i V 2 η i V s i V + 1 2 G i V q ¯ 2 i V 2 + 1 μ i V ρ i V 1 2 μ i V ς i V μ ˜ i V 2 + 1 2 ς i V μ i V + z i T z ˙ i + 1 2 d ˜ i U 2 ( X ) + 1 2 d ˜ i V 2 ( X ) .
where D i U = G i U q ¯ 2 i U 2 / 2 + ρ i U / μ i U + ς i U μ i U / 2 , D i V = G i V q ¯ 2 i V 2 / 2 + ρ i V / μ i V + ς i V μ i V / 2 .
Equation (62) can be reshaped into
V ˙ 3 i = 1 N l i U 1 2 s i U 2 η i U s i U 1 2 μ i U ς i U μ ˜ i U 2 + D i U l i V 1 2 s i V 2 η i V s i V 1 2 μ i V ς i V μ ˜ i V 2 + D i V + z i T z ˙ i + 1 2 d ˜ i U 2 ( X ) + 1 2 d ˜ i V 2 ( X ) .
Substituting (44) into (63), it can be seen that
V ˙ 3 i = 1 N l i U 1 2 s i U 2 l i V 1 2 s i V 2 1 2 μ i U ς i U μ ˜ i U 2 1 2 μ i V ς i V μ ˜ i V 2 z i T χ i z i + D i U + D i V z i T σ ˜ i + 1 2 d ˜ i U 2 ( X ) + 1 2 d ˜ i V 2 ( X ) ,
it follows that
V ˙ 3 i = 1 N l i U 1 2 s i U 2 l i V 1 2 s i V 2 1 2 μ i U ς i U μ ˜ i U 2 1 2 μ i V ς i V μ ˜ i V 2 λ min χ i 1 2 z i 2 + D i U + D i V + 1 2 σ ˜ i 2 + 1 2 d ˜ i U 2 ( X ) + 1 2 d ˜ i V 2 ( X ) .
where z = [ z 1 T , z 2 T , , z N T ] T , s U = [ s 1 U , s 2 U , , s N U ] T , s V = [ s 1 V , s 2 V , , s N V ] T , μ ˜ U = [ μ ˜ 1 U , μ ˜ 2 U , , μ ˜ N U ] T , μ ˜ V = [ μ ˜ 1 V , μ ˜ 2 V , , μ ˜ N V ] T , ħ 3 = min i = 1 , , N [ l i U 1 2 , l i V 1 2 , ς i U / 2 μ i U , ς i V / 2 μ i V , λ min ( χ i ) 1 2 ] , E 3 = [ s U , s V , μ ˜ U , μ ˜ V , z ] .
This indicates that
V ˙ 3 ( 1 θ ¯ 3 ) ħ 3 E 3 2 θ ¯ 3 ħ 3 E 3 2 + i = 1 N ( D i U + D i V + 1 2 E i 1 2 + 1 2 E i 2 2 ) ,
noting that
E 3 i = 1 N 2 D i U + 2 D i V + E i 1 2 + E i 2 2 2 θ ¯ 3 ħ 3 ,
renders
V ˙ 3 1 θ ¯ 3 ħ 3 E 3 2 ,
where 0 < θ ¯ 2 < 1 ; it is concluded that the subsystem (44) is input-state stable, and
E 3 ( t ) λ max P c λ min P c max E 3 t 0 e ϱ 3 t t 0 i = 1 N 2 D i U + 2 D i V + E i 1 2 + E i 2 2 2 θ ¯ 3 ħ 3 , t t 0 ,
where ϱ 3 = 2 ħ 3 1 θ ¯ 3 / λ max P c and P c = diag 1 , 1 2 γ 1 μ 1 U , , 1 2 γ 1 μ N U , 1 2 γ 2 μ 1 V , , 1 2 γ 2 μ N V .
The stability of the closed-loop control system is analyzed as follows. □
Theorem 4.
Consider (16), (27), (42) and (43), the cascade system is input-state stable, and all signals are uniformly ultimately bounded.
Proof. 
Considering the stability theorem for cascaded systems proposed by Sontag and Wang in [44] and combining Theorems 1 and 2, it follows that the closed-loop system which is composed of (25) and (44) is input-state stable. From (51) and (69), E 3 ( t ) is eventually bounded by
E 3 ( t ) λ max P c λ min P c i = 1 N 2 D i U + 2 D i V + E i 2 2 2 θ ¯ 3 ħ 3 + P i B i 2 σ * 2 2 θ ¯ 3 ħ 3 θ ¯ i 1 2 .
This proof is complete. □
Theorem 5.
Suppose Assumptions 1 and 2 are valid and the condition v a z d p ˙ a z d = 0 is met. Under these circumstances, if the control scheme for formation tracking along with the adaptive law is implemented as outlined in (46), it leads to the UAV successfully tracking the desired height signal. Additionally, the tracking error e a z p is ensured to be uniformly ultimately bounded.
Proof. 
The time derivative of (45) can be obtained as follows:
e ˙ a z p = e a z v , e ˙ a z v = F a z + G a z u a z + Δ a z v ˙ a z d .
In system (71), regarding e a z v as the virtual control input. Design the virtual control signal ζ z = k ζ e a z p to ensure the stability of system (71).
Design the following Lyapunov function:
V a z p = 1 2 e a z p T e a z p ,
then, the time derivative of (72) can be obtained as follows:
V ˙ a z p = k ζ e a z p T e a z p .
Define a new error as
e ζ = e a z v ζ z ,
taking the time derivative of (74) and substituting (46) into it, it follows that
e ˙ ζ = σ ζ e ζ + F a z + Δ a z + G ˜ a z u a z e ζ κ ^ 11 h ( X ) T h ( X ) 2 ς 1 2 e ζ κ ^ 12 2 ς 2 2 .
Define the following Lyapunov function:
V z = 1 2 e a z p T e a z p + 1 2 e ζ T e ζ + κ ˜ i 11 2 2 ι 11 + κ ˜ i 12 2 2 l 12 + Tr G ˜ a z T G ˜ a z 2 l 13 .
where κ ˜ 11 = κ 11 κ ^ 11 , κ ˜ 12 = κ 12 κ ^ 12 , and G ˜ a z = G a z G ^ a z . The time derivative of (76) can be obtained as follows:
V ˙ z k ζ e a z p T e a z p σ z e ζ T e ζ + e ζ T F a z + e ζ T Δ a z + e ζ T G ˜ a z u a z e ζ T e ζ κ ^ 11 h ( X ) T h ( X ) 2 ς i 1 2 e ζ T e ζ κ ^ 12 2 ς i 2 2 + κ ˜ 11 κ ˜ ˙ 11 ι 11 + κ ˜ 12 κ ˜ ˙ 12 ι 12 + Tr G ˜ a z T G ˜ ˙ a z l 13 .
Since F a z is an unknown term, the RBFNN system is used to approximate it, F a z = W a * h ( X ) + ε a , where W a * represents the ideal weight, ε a represents the estimation error, ε a ε ¯ a , and h ( X ) denotes Gaussian functions.
Noting that
e ζ T W a * h ( X ) e ζ T e ζ κ 11 h ( X ) T h ( X ) 2 ς i 1 2 + ς i 1 2 2 , e ζ T ε a + Δ a z e ζ T e ζ κ i 12 2 ς i 2 2 + ς i 2 2 2 ,
where κ 11 = W a * T W a * and κ 12 = ε ¯ a + Δ ¯ a z T ε ¯ a + Δ ¯ a z .
Substituting (78) into (77), it follows that
V ˙ z k ζ e a z p T e a x p σ z e ζ T e ζ + e ζ T e ζ κ 11 h ( X ) T h ( X ) 2 ς 1 2 + e ζ T e ζ κ 12 2 ς 2 2 e ζ T e ζ κ ^ 11 h T ( X ) h ( X ) 2 ς 1 2 e ζ T e ζ κ ^ i 2 2 ς 2 2 κ ˜ 12 k 12 κ ^ 12 + e ζ T e ζ 2 ς 2 2 κ ˜ 11 k 11 κ ^ 11 + e ζ T e ζ h T ( X ) h ( X ) 2 ς 1 2 + e ζ T G ˜ a z u a z Tr G ˜ a z T k 13 G ^ a z + e ζ u a z + ς 1 2 2 + ς 2 2 2 k ζ e a z p T e a z p σ z e ζ T e ζ k 11 2 κ ˜ 11 2 k 12 2 κ ˜ 12 2 k 13 2 G ˜ a z F 2 + k 11 2 κ 11 2 + k 12 2 κ 12 2 + k 13 2 G a z F 2 + ς 1 2 2 + ς 2 2 2 ϑ z V z + v z .
where ϑ z = min 2 k ζ , 2 σ z , ι 11 k 11 , ι 12 k 12 , ι 13 k 13 > 0 , v z = k 11 2 κ 11 2 + k 12 2 κ 12 2 + k 13 2 G a z F 2 + ς 1 2 2 + ς 2 2 2 .
According to (79), we can obtain
V ˙ z ϑ z V z + v z .
The altitude tracking errors and velocity errors of the UAV are uniformly ultimately bounded. □

5. Illustrative Examples

In this section, we demonstrate the effectiveness of the controller presented in this paper by constructing a formation system comprising one UAV and four USVs, as outlined below. The model parameters of the HMAS are given in Table 1.
If B = d i a g { 0 , 0 , 0 , 0 , 1 } is the adjacency weight matrix between agents and the virtual leader, then the adjacency matrix A and the Laplacian matrix L are as follows:
A = 0 1 1 0 1 1 0 1 1 0 1 1 0 0 0 0 1 0 0 0 1 0 0 0 0 L = 3 1 1 0 1 1 3 1 1 0 1 1 2 0 0 0 1 0 1 0 1 0 0 0 1
Case 1. Let leader UAV navigate along the time-varying parametric trajectory p 0 = [ 1.2 t , 0.2 t ] T . The initial states of five HMAS are set as p 1 = [ 0 , 0 , 0 ] T , p 2 = [ 7 , 7 , 0 ] T , p 3 = [ 7 , 7 , 8 ] T , p 4 = [ 14 , 14 , 0 ] T , p 5 = [ 14 , 14 , 0 ] T . In the XY plane, the desired formation pattern is defined as p 1 d = [ 0 , 0 ] T , p 2 d = [ 7 , 7 ] T , p 3 d = [ 7 , 7 ] T , p 4 d = [ 14 , 14 ] T , p 5 d = [ 14 , 14 ] T .
According to Equations (30), (39), (40), (42) and (43), the controller parameters are selected sequentially as c i U = c i V = 10 , l i U = l i V = 30 , η i U = η i V = 2 , ρ i U = ρ i V = 0.02 , ς i U = ς i V = 0.2 , γ 1 = γ 2 = 2 , k i U = k i V = 0.2 . The external disturbances are given as Δ a x z = 0.3 cos ( 0.2 t ) , Δ a x y = [ 0.2 cos ( 0.5 t ) , 0.8 cos ( t ) ] T , ω d i x y = [ 1.1 cos ( 0.5 t ) , 0.2 sin ( 2 t ) ] T , p a z = 8 . The total duration of the simulation run is T z = 400 s .
Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11 show the simulation results for the underactuated HMAS by using the proposed controller with input quantization. The trajectories of the HMAS in the 3D environment are depicted in Figure 3. It is clear from the figure that the HMAS set off from a given initial position, maintain the desired distance from each other, and maintain a special formation under the influence of external interference and system internal error. Figure 4 and Figure 5 depict the longitudinal and lateral velocity of four USVs and one UAV in the XY plane. As can be seen from the figure, the speed synchronization between the USV and the UAV is achieved in around 30s. Figure 6 and Figure 7 describe the formation tracking errors of four USVs and one UAV in the XY plane in both the lateral and longitudinal directions. It is clear that the HMAS can accomplish the formation mission under a small error. Figure 8 describes the longitudinal control input τ i U and Q ( τ i U ) with and without input quantization. Figure 9 describes the lateral control input τ i V and Q ( τ i V ) with and without input quantization, respectively. Under comparison, the quantization process is able to guarantee the control signal can maintain the quality of the control scheme without losing control performance. The sliding mode function of USVs and UAV is depicted in Figure 11. The performance of ESO in estimating d i ( X ) is depicted in Figure 11, which clearly demonstrates the effectiveness of the proposed solution for the accurate compensation of uncertainties and external disturbances within heterogeneous systems.
Case 2. To further discuss the effectiveness of the proposed controller, we compare the proposed quantization controller with the strategy proposed in [27,45,46]. Let leader UAV navigate along the reference trajectory p 0 = [ 0.2 t , 2 sin 0.1 t 3 ] T . The other initial conditions and control design parameters are chosen the same as their counterparts in Case 1.
The simulation results are presented in Figure 12, Figure 13, Figure 14, Figure 15 and Figure 16. Figure 12 shows the formation trajectories of the five HMAS. As can be seen from the figure, in the XY plane the five HMAS set off from a given initial position, maintain the desired distance from each other, and maintain a special formation under the influence of external interference and system internal error. Figure 13 and Figure 14 depict the comparison curve of longitudinal control input Q ( τ U ( t ) ) and lateral control input Q ( τ V ( t ) ) in different methods, respectively. Figure 15 and Figure 16 present the longitudinal and lateral formation tracking error comparisons between the two control strategies, respectively. The analysis of these figures reveals that the implementation of a linear time-varying model for representing the quantizers leads to reduced tracking errors in the proposed control strategy, thereby enhancing its tracking accuracy.

6. Conclusions

This paper outlines a distributed time-varying formation control framework tailored for USVs and UAVs. The process commences with the unification of the heterogeneous agent models for UAVs and USVs, Subsequently, a trajectory-guidance law incorporating ESO considerations is devised, enabling the pursuit of time-varying trajectories on the XY plane. Within the dynamic subsystems, ESO mechanisms are harnessed to identify unknown system functions, facilitating the development of linear time-varying models that aptly characterize quantizers. Remarkably, this approach obviates the necessity for quantizer parameter information. The closed-loop system’s input-state stability is rigorously demonstrated, affirming the system’s robust performance. In culmination, the efficacy of the proposed adaptive distributed formation tracking controller is rigorously verified through comprehensive simulation results.
The scope of this paper is currently limited to the heterogeneous formation consisting of a single UAV and multiple USVs. Future research should extend this work to address the cooperative dynamics of multiple UAVs and USVs in more complex scenarios. To enhance the effectiveness of the proposed algorithm, future investigations will aim to integrate transformations between formations in the presence of both dynamic and static obstacles. This integration will not only improve the algorithm’s robustness but also validate its performance in terms of network lifetime and communication delay, which are crucial for real-world applications. Additionally, our future objectives include the development of methodologies to reduce communication overload and conserve communication resources within a quantized environment. By incorporating performance parameters such as network lifetime and communication delay, we aim to provide a more comprehensive evaluation of the system’s effectiveness and practicality.

Author Contributions

Conceptualization, W.L. and J.N.; methodology, Y.H., Z.L. and J.N.; software, Y.H. and Z.L.; validation, Y.H.; writing—original draft preparation, Y.H. and Z.L.; writing—review and editing, Y.H., W.L., J.N., Z.L. and X.Y.; supervision, W.L. and J.N.; project administration, W.L. and Z.L.; funding acquisition, X.Y. and J.N. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grants 51939001, 61976033, and in part by the Liaoning Revitalization Talents Program under Grant XLYC 1908018, XLYC 1807046, and part by the Fundamental Research Funds for the Central Universities under Grant 3132023151.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Yang, Z.; Li, S.; Yu, D.; Chen, C.P. BLS-based formation control for nonlinear multi-agent systems with actuator fault and input saturation. Nonlinear Dyn. 2022, 109, 2657–2673. [Google Scholar] [CrossRef]
  2. Ma, L.; Wang, Y.L.; Fei, M.R.; Pan, Q.K. Cross-dimensional formation control of second-order heterogeneous multi-agent systems. ISA Trans. 2022, 127, 188–196. [Google Scholar] [CrossRef] [PubMed]
  3. Ennong, T.; Ye, L.; Teng, M.; Yulei, L.; Yueming, L.; Jian, C. Design and experiment of a sea-air heterogeneous unmanned collaborative system for rapid inspection tasks at sea. Appl. Ocean Res. 2024, 143, 103856. [Google Scholar] [CrossRef]
  4. Hou, H.Q.; Liu, Y.J.; Liu, L.; Lan, J. Adaptive fuzzy formation control for heterogeneous multi-agent systems using time-varying IBLFs. Nonlinear Dyn. 2023, 111, 16077–16091. [Google Scholar] [CrossRef]
  5. Mazare, M.; Asharioun, H.; Davoudi, E.; Mokhtari, M. Distributed finite-time neural network observer-based consensus tracking control of heterogeneous underwater vehicles. Ocean Eng. 2023, 272, 113882. [Google Scholar] [CrossRef]
  6. Ren, H.; Ma, H.; Li, H.; Wang, Z. Adaptive Fixed-Time Control of Nonlinear MASs With Actuator Faults. IEEE/CAA J. Autom. Sin. 2023, 10, 1252–1262. [Google Scholar] [CrossRef]
  7. Liu, Z.; Nojavanzadeh, D.; Saberi, A.; Stoorvogel, A.A. Scale-Free Collaborative Protocol Design for Output Synchronization of Heterogeneous Multi-Agent Systems With Nonuniform Communication Delays. IEEE Trans. Netw. Sci. Eng. 2022, 9, 2882–2894. [Google Scholar] [CrossRef]
  8. Ren, H.; Cheng, Z.; Qin, J.; Lu, R. Deception attacks on event-triggered distributed consensus estimation for nonlinear systems. Automatica 2023, 154, 111100. [Google Scholar] [CrossRef]
  9. Khoshnevisan, L.; Liu, X. Resilient neural network-based control of nonlinear heterogeneous multi-agent systems: A cyber-physical system approach. Nonlinear Dyn. 2023, 111, 19171–19185. [Google Scholar] [CrossRef]
  10. Sun, F.; Wu, X.; Kurths, J.; Zhu, W. Group consensus of heterogeneous multi-agent systems with packet loss and unknown speed of second-order agents in cooperative–competitive networks. Nonlinear Dyn. 2022, 110, 3447–3461. [Google Scholar] [CrossRef]
  11. Liu, Y.; Wang, J.; Shi, Y.; He, Z.; Liu, F.; Kong, W.; He, Y. Unmanned airboat technology and applications in environment and agriculture. Comput. Electron. Agric. 2022, 197, 106920. [Google Scholar] [CrossRef]
  12. Cui, Y.; Peng, L.; Li, H. Filtered probabilistic model predictive control-based reinforcement learning for unmanned surface vehicles. IEEE Trans. Ind. Inform. 2022, 18, 6950–6961. [Google Scholar] [CrossRef]
  13. Ning, J.; Ma, Y.; Liu, L.; Li, T. Fuzzy course tracking control of an unmanned surface vehicle with event-triggered mechanism and input quantization. Ocean Eng. 2023, 288, 115704. [Google Scholar] [CrossRef]
  14. Mazare, M. Distributed adaptive fault tolerant formation control for multiple underwater vehicles: Free-will arbitrary time approach. Ocean Eng. 2023, 286, 115601. [Google Scholar] [CrossRef]
  15. Huang, Y.; Li, W.; Ning, J.; Li, Z. Formation Control for UAV-USVs Heterogeneous System with Collision Avoidance Performance. J. Mar. Sci. Eng. 2023, 11, 2332. [Google Scholar] [CrossRef]
  16. Huang, C.; Zhang, X.; Zhang, G.; Deng, Y. Robust practical fixed-time leader–follower formation control for underactuated autonomous surface vessels using event-triggered mechanism. Ocean Eng. 2021, 233, 109026. [Google Scholar] [CrossRef]
  17. Li, J.; Zhang, G.; Li, B. Robust adaptive neural cooperative control for the USV-UAV based on the LVS-LVA guidance principle. J. Mar. Sci. Eng. 2022, 10, 51. [Google Scholar] [CrossRef]
  18. Tan, G.; Zhuang, J.; Zou, J.; Wan, L. Coordination control for multiple unmanned surface vehicles using hybrid behavior-based method. Ocean Eng. 2021, 232, 109147. [Google Scholar] [CrossRef]
  19. Huang, D.; Li, H.; Li, X. Formation of Generic UAVs-USVs System Under Distributed Model Predictive Control Scheme. IEEE Trans. Circuits Syst. II Express Briefs 2020, 67, 3123–3127. [Google Scholar] [CrossRef]
  20. Xue, K.; Wu, T. Distributed consensus of USVs under heterogeneous UAV-USV multi-agent systems cooperative control scheme. J. Mar. Sci. Eng. 2021, 9, 1314. [Google Scholar] [CrossRef]
  21. Liu, H.; Weng, P.; Tian, X.; Mai, Q. Distributed adaptive fixed-time formation control for UAV-USV heterogeneous multi-agent systems. Ocean Eng. 2023, 267, 113240. [Google Scholar] [CrossRef]
  22. Li, S.; Wang, X.; Wang, S.; Zhang, Y. Distributed Bearing-Only Formation Control for UAV-UWSV Heterogeneous System. Drones 2023, 7, 124. [Google Scholar] [CrossRef]
  23. Liu, S.; Jiang, B.; Mao, Z.; Ma, Y. Adaptive Fault-Tolerant Formation Control of Heterogeneous Multi-Agent Systems under Directed Communication Topology. Sensors 2022, 22, 6212. [Google Scholar] [CrossRef] [PubMed]
  24. Li, H.; Li, X. Distributed consensus of heterogeneous linear time-varying systems on UAVs–USVs coordination. IEEE Trans. Circuits Syst. II Express Briefs 2019, 67, 1264–1268. [Google Scholar] [CrossRef]
  25. Hao, L.Y.; Yang, G.H. Robust fault tolerant control based on sliding mode method for uncertain linear systems with quantization. ISA Trans. 2013, 52, 600–610. [Google Scholar] [CrossRef] [PubMed]
  26. Hao, L.Y.; Zhang, H.; Li, H.; Li, T.S. Sliding mode fault-tolerant control for unmanned marine vehicles with signal quantization and time-delay. Ocean Eng. 2020, 215, 107882. [Google Scholar] [CrossRef]
  27. Hao, L.Y.; Zhang, H.; Guo, G.; Li, H. Quantized sliding mode control of unmanned marine vehicles: Various thruster faults tolerated with a unified model. IEEE Trans. Syst. Man Cybern. Syst. 2019, 51, 2012–2026. [Google Scholar] [CrossRef]
  28. Hayakawa, T.; Ishii, H.; Tsumura, K. Adaptive quantized control for linear uncertain discrete-time systems. Automatica 2009, 45, 692–700. [Google Scholar] [CrossRef]
  29. Wang, C.; Wen, C.; Lin, Y.; Wang, W. Decentralized adaptive tracking control for a class of interconnected nonlinear systems with input quantization. Automatica 2017, 81, 359–368. [Google Scholar] [CrossRef]
  30. Xiong, S.; Wu, Q.; Wang, Y. Distributed Coordination of Heterogeneous Multi-agent Systems with Dynamic Quantization and L 2- L Control. Int. J. Control Autom. Syst. 2020, 18, 2468–2481. [Google Scholar] [CrossRef]
  31. Zheng, B.C.; Yang, G.H. Quantized output feedback stabilization of uncertain systems with input nonlinearities via sliding mode control. Int. J. Robust Nonlinear Control 2014, 24, 228–246. [Google Scholar] [CrossRef]
  32. Choi, Y.H.; Yoo, S.J. Distributed quantized feedback design strategy for adaptive consensus tracking of uncertain strict-feedback nonlinear multiagent systems with state quantizers. IEEE Trans. Cybern. 2021, 52, 7069–7083. [Google Scholar] [CrossRef] [PubMed]
  33. Fang, X.; Zhou, J.; Wen, G. Location game of multiple unmanned surface vessels with quantized communications. IEEE Trans. Circuits Syst. II Express Briefs 2021, 69, 1322–1326. [Google Scholar] [CrossRef]
  34. Salam, A.; Javaid, Q.; Ahmad, M.; Wahid, I.; Arafat, M.Y. Cluster-Based Data Aggregation in Flying Sensor Networks Enabled Internet of Things. Future Internet 2023, 15, 279. [Google Scholar] [CrossRef]
  35. Liu, H.; Meng, Q.; Peng, F.; Lewis, F.L. Heterogeneous formation control of multiple UAVs with limited-input leader via reinforcement learning. Neurocomputing 2020, 412, 63–71. [Google Scholar] [CrossRef]
  36. Alam, M.M.; Arafat, M.Y.; Moh, S.; Shen, J. Topology control algorithms in multi-unmanned aerial vehicle networks: An extensive survey. J. Netw. Comput. Appl. 2022, 207, 103495. [Google Scholar] [CrossRef]
  37. Cheng, W.; Zhang, K.; Jiang, B.; Simani, S. Neural Network Observer-Based Prescribed-Time Fault-Tolerant Tracking Control for Heterogeneous Multiagent Systems With a Leader of Unknown Disturbances. IEEE Trans. Aerosp. Electron. Syst. 2023, 59, 9042–9053. [Google Scholar] [CrossRef]
  38. Lv, C.; Chen, J.; Yu, H.; Chi, J.; Yang, Z. Adaptive NN state error PCH trajectory tracking control for unmanned surface vessel with uncertainties and input saturation. Asian J. Control 2023, 25, 3903–3919. [Google Scholar] [CrossRef]
  39. Huang, T.; Chen, Z.; Gao, W.; Xue, Z.; Liu, Y. A USV-UAV Cooperative Trajectory Planning Algorithm with Hull Dynamic Constraints. Sensors 2023, 23, 1845. [Google Scholar] [CrossRef]
  40. Wang, N.; Ahn, C.K. Coordinated trajectory-tracking control of a marine aerial-surface heterogeneous system. IEEE/ASME Trans. Mechatron. 2021, 26, 3198–3210. [Google Scholar] [CrossRef]
  41. Li, J.; Zhang, G.; Shan, Q.; Zhang, W. A novel cooperative design for USV-UAV systems: 3D mapping guidance and adaptive fuzzy control. IEEE Trans. Control Netw. Syst. 2023, 10, 564–574. [Google Scholar] [CrossRef]
  42. Chen, F.; Jiang, R.; Zhang, K.; Jiang, B.; Tao, G. Robust backstepping sliding-mode control and observer-based fault estimation for a quadrotor UAV. IEEE Trans. Ind. Electron. 2016, 63, 5044–5056. [Google Scholar] [CrossRef]
  43. Park, B.S.; Kwon, J.W.; Kim, H. Neural network-based output feedback control for reference tracking of underactuated surface vessels. Automatica 2017, 77, 353–359. [Google Scholar] [CrossRef]
  44. Sontag, E.D.; Wang, Y. On characterizations of the input-to-state stability property. Syst. Control Lett. 1995, 24, 351–359. [Google Scholar] [CrossRef]
  45. Hao, L.Y.; Zhao, Z.H. Quantized Sliding Mode Fault-Tolerant Control for Unmanned Marine Vehicles with Thruster Saturation. J. Mar. Sci. Eng. 2023, 11, 309. [Google Scholar] [CrossRef]
  46. Hao, L.Y.; Yu, Y.; Li, T.S.; Li, H. Quantized output-feedback control for unmanned marine vehicles with thruster faults via sliding-mode technique. IEEE Trans. Cybern. 2021, 52, 9363–9376. [Google Scholar] [CrossRef]
Figure 1. Motion model of USV.
Figure 1. Motion model of USV.
Jmse 12 00975 g001
Figure 2. Formation control flowchart of HMAS.
Figure 2. Formation control flowchart of HMAS.
Jmse 12 00975 g002
Figure 3. Motion model of USVs and UAV in case 1.
Figure 3. Motion model of USVs and UAV in case 1.
Jmse 12 00975 g003
Figure 4. Longitudinal velocity of USVs and UAV.
Figure 4. Longitudinal velocity of USVs and UAV.
Jmse 12 00975 g004
Figure 5. Lateral velocity of USVs and UAV.
Figure 5. Lateral velocity of USVs and UAV.
Jmse 12 00975 g005
Figure 6. Longitudinal formation tracking errors.
Figure 6. Longitudinal formation tracking errors.
Jmse 12 00975 g006
Figure 7. Lateral formation tracking errors.
Figure 7. Lateral formation tracking errors.
Jmse 12 00975 g007
Figure 8. Comparison of longitudinal control inputs.
Figure 8. Comparison of longitudinal control inputs.
Jmse 12 00975 g008
Figure 9. Comparison of lateral control inputs.
Figure 9. Comparison of lateral control inputs.
Jmse 12 00975 g009
Figure 10. Sliding mode function of USVs and UAV.
Figure 10. Sliding mode function of USVs and UAV.
Jmse 12 00975 g010
Figure 11. Estimated performance of ESO.
Figure 11. Estimated performance of ESO.
Jmse 12 00975 g011
Figure 12. Motion model of USVs and UAV in case 2.
Figure 12. Motion model of USVs and UAV in case 2.
Jmse 12 00975 g012
Figure 13. Comparison of longitudinal control inputs.
Figure 13. Comparison of longitudinal control inputs.
Jmse 12 00975 g013
Figure 14. Comparison of lateral control inputs.
Figure 14. Comparison of lateral control inputs.
Jmse 12 00975 g014
Figure 15. Comparison of longitudinal formation tracking errors.
Figure 15. Comparison of longitudinal formation tracking errors.
Jmse 12 00975 g015
Figure 16. Comparison of lateral formation tracking errors.
Figure 16. Comparison of lateral formation tracking errors.
Jmse 12 00975 g016
Table 1. Model Parameters.
Table 1. Model Parameters.
ParameterValue
m a 2
g a 9.8
J a x , J a y , J a z 1.5
d a x , d a y , d a z 0.012
m μ s i 25.8
m v s i 33.8
m r s i 2.76
d μ s i 0.725
d v s i 0.89
d r s i 1.9
d μ b i 1 1.33
d v b i 1 36.47
d r b i 1 0.75
d r b i 1.9
f i v 36.5 v v 0.8896 v 0.805 v r m i μ μ r
f i r 0.75 r r 1.90 r + 0.08 v r
+ ( m i u m i v ) μ v 1.0948 μ r
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ning, J.; Huang, Y.; Liu, Z.; Li, W.; Yue, X. Adaptive Distributed Heterogeneous Formation Control for UAV-USVs with Input Quantization. J. Mar. Sci. Eng. 2024, 12, 975. https://doi.org/10.3390/jmse12060975

AMA Style

Ning J, Huang Y, Liu Z, Li W, Yue X. Adaptive Distributed Heterogeneous Formation Control for UAV-USVs with Input Quantization. Journal of Marine Science and Engineering. 2024; 12(6):975. https://doi.org/10.3390/jmse12060975

Chicago/Turabian Style

Ning, Jun, Yuyang Huang, Zihan Liu, Wei Li, and Xingwang Yue. 2024. "Adaptive Distributed Heterogeneous Formation Control for UAV-USVs with Input Quantization" Journal of Marine Science and Engineering 12, no. 6: 975. https://doi.org/10.3390/jmse12060975

APA Style

Ning, J., Huang, Y., Liu, Z., Li, W., & Yue, X. (2024). Adaptive Distributed Heterogeneous Formation Control for UAV-USVs with Input Quantization. Journal of Marine Science and Engineering, 12(6), 975. https://doi.org/10.3390/jmse12060975

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop