Next Article in Journal
UAV Relay Energy Consumption Minimization in an MEC-Assisted Marine Data Collection System
Next Article in Special Issue
Collision-Free Formation-Containment Tracking of Multi-USV Systems with Constrained Velocity and Driving Force
Previous Article in Journal
Adaptive Fuzzy Quantized Control for a Cooperative USV-UAV System Based on Asynchronous Separate Guidance
Previous Article in Special Issue
A Multi-Ship Collision Avoidance Algorithm Using Data-Driven Multi-Agent Deep Reinforcement Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Formation Control for UAV-USVs Heterogeneous System with Collision Avoidance Performance

1
College of Navigation, Dalian Maritime University, Dalian 116026, China
2
College of Information Technology, Jilin Normal University, Siping 136000, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2023, 11(12), 2332; https://doi.org/10.3390/jmse11122332
Submission received: 9 October 2023 / Revised: 24 November 2023 / Accepted: 7 December 2023 / Published: 10 December 2023
(This article belongs to the Special Issue Maritime Autonomous Surface Ships)

Abstract

:
This paper investigates the cooperative formation trajectory tracking problem for heterogeneous unmanned aerial vehicle (UAV) and multiple unmanned surface vessel (USV) systems with collision avoidance performance. Firstly, a formation control protocol based on extended state observer (ESO) is proposed to ensure that the UAV and the USVs track the target trajectory simultaneously in the XY plane. Then, the collision avoidance control strategy of USV formation based on artificial potential field (APF) theory is designed. Specifically, the APF method is improved by reconstructing the repulsive potential field to make the collision avoidance action of USVs more in line with the requirements of International Regulations for Preventing Collisions at Sea (COLREGs). Following that, an altitude controller for the UAV is proposed to maintain the cooperative formation of the heterogeneous systems. Based on the input-to-state stability, the stability of the proposed control structure is proven, 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 the past few years, unmanned systems have gained significant prominence in the realm of industrial advancement, attributed to their diverse usage in areas such as reconnaissance [1,2], marine pollution tracking [3,4], and advanced traffic management [5,6]. These systems principally encompass unmanned aerial vehicles (UAVs), unmanned ground vehicles (UGVs), unmanned surface vehicles (USVs), and autonomous underwater vehicles (AUVs) [7,8]. Although each unmanned system can handle tasks, its capacity limits its ability to deal with more complex tasks. Therefore, a homogeneous or heterogeneous system is more effective than a single unmanned system in completing intricate tasks [9,10,11].
The development of USV has been widely employed in ocean engineering due to its benefits, which include increased loading capacity, more convenience, and cheaper mission costs [12,13]. However, finding the rescue target in rescue missions and maritime searches is difficult owing to the restricted observation range of USVs [14,15]. Fortunately, incorporating UAVs into USV systems to construct heterogeneous multi-agent systems can compensate for this shortcoming [16,17]. A USV usually boosts its capacity to operate at sea by using the UAV’s flexibility and utilizing the powerful target search ability while simultaneously widening the communication range of the sea between USVs, therefore, the heterogeneous systems extend the effective working area [18,19]. As a result, it is critical to investigate the formation challenges of a heterogeneous system, which has primarily engineering importance.
For the cooperative tracking control problem of UAVs and USVs, many research results have been addressed, among which typical control methods include the leader following method [20], virtual structure method [21], behavior-based approach [22], and model prediction method [23]. In [24], a series of coordinate transformations have been developed to convert the tracking error dynamics of a marine aerial-surface heterogeneous system into translation–rotation cascade forms. This eliminates variability between the UAV and the USV in terms of both translation and rotation dynamics. In [25], the research proposes an adaptive fault-tolerant time-varying formation control scheme. This scheme is particularly tailored to address challenges posed by actuator failures, parameter uncertainties, and external disturbances within the framework of a directed communication topology.
In the area of cooperative control of heterogeneous formations of UAVs and USVs, a globally applicable fixed-duration adaptive neural network is designed in [18] that employs a nonsingular rapid terminal sliding approach for formation control. This protocol is engineered to accurately follow a specified trajectory and achieve a predefined formation arrangement within a set timeframe, effectively addressing a range of uncertainties. This is accomplished by harnessing the combined strengths of adaptive methodologies and Radial Basis Function (RBF) neural networks. In [26], the research proposes a velocity estimation-based control strategy that comprises a distributed observer for estimating each vehicle’s reference velocity. In [23], a distributed model predictive control algorithm tailored for heterogeneous systems is presented, characterized by a directed topological structure.
It is worth noting that the above research results tend to focus on the establishment of air–sea cooperation and the realization of trajectory tracking without considering the problem of collision avoidance [27,28]. During the mission, a minimum safe distance has to be maintained between any two unmanned systems to prevent collisions and to avoid damage to the overall performance of the search and rescue cooperative mission. Therefore, another key issue for the heterogeneous formation system is collision avoidance. Collision avoidance includes not only avoiding collisions between unmanned systems or between unmanned systems and other individuals but also avoiding collisions between unmanned systems and obstacles. Collision avoidance can be further divided into static collision avoidance based on sensor information and dynamic collision avoidance based on an unknown environment. Collision avoidance between aerial vehicles and obstacles can be divided into two ways: overall collision avoidance and changing formation collision avoidance.
There are relatively few research results on the cooperative formation control of heterogeneous systems with collision avoidance and obstacle avoidance. To avoid collision between ships, Refs. [29,30] propose an unmanned ship formation control method based on guaranteed performance, but the collision avoidance of obstacles is not addressed. For the collision avoidance problem between UAVs with static and dynamic obstacles, a cooperative controller for multiple unmanned ships based on artificial potential function is proposed in [12,31], while in the area of cooperative collision avoidance and obstacle avoidance for heterogeneous systems, research results are more limited. In [32], a distributed heterogeneous cooperative tracking control approach is proposed based on the leader-following method, and the artificial potential field (APF) algorithm is used to construct a control strategy with a collision avoidance mechanism. However, the research results between vehicles ignore the nautical engineering practice and poorly consider the COLREGs.
Based on the preceding analysis, this paper investigates the cooperative formation trajectory tracking problem for UAV-USV heterogeneous systems with model uncertainty and external disturbances. To solve this problem, a formation control protocol based on extended state observer (ESO) is proposed to ensure that UAV and USV track the target trajectory simultaneously. The collision avoidance control strategy of USV formation based on improved APF theory is designed. The collision avoidance problem between multiple USV formations formed under UAV coordination is solved by innovatively introducing the ship encounter situation and danger evaluation index into the artificial potential field. The key contributions can be summarized as follows:
(1)
Compared with the existing results in [12,29,30,31], which only study the cooperative trajectory tracking problem of UAV and USV heterogeneous systems, this paper explores the collision avoidance protocol for USV formation under UAV cooperation with navigation practice.
(2)
Compared with the existing results in [18,25,26,32], which estimate that the system’s indeterminate terms rely on RBF neural networks and fuzzy logic, etc., this paper employs ESO to realize the accurate compensation of uncertainties and external disturbances in heterogeneous systems.
(3)
Compared with the existing results in [31,32,33,34,35,36,37], this paper innovatively introduces the ship encounter situation and danger evaluation index into the APF approach, and the improved APF method for heterogeneous cooperative control collision avoidance decision is more in line with the navigation practice.
The organization of this paper is structured as follows. Section 2 formulates the problem. Section 3 expounds on the controller design and the closed-loop system stability analysis. Section 4 demonstrates the simulation illustrations. Section 5 concludes this paper.

2. Preliminaries and Problem Statement

2.1. Problem Formulation

Consider the heterogeneous systems consisting of one UAV and N USVs. First, the dynamic models of the UAV and N USVs are presented. They are used to illustrate a unified dynamic model for the heterogeneous system. For ease of use, let Π = { 1 , 2 , , N } . According to the results in [38], the dynamic model of the quadrotor UAV can be described as
p ¨ a x = sin ϕ a sin ψ a + cos ϕ a sin θ a cos ψ a u a p m a d a x p ˙ a x m a + Δ a x , p ¨ a y = sin ϕ a cos ψ a + cos ϕ a sin θ a sin ψ a u a p m a d a y p ˙ a y m a + Δ a y , p ¨ a z = cos θ a cos ϕ a u a p m a d a z p ˙ a z m a g a + Δ a z ,
ϕ a ¨ = θ a ˙ ψ a ˙ J a y J a z J a x J a r J a x θ a ˙ d ¯ a + τ a ϕ J a x d a ϕ ϕ a ˙ J a x , θ a ¨ = ϕ a ˙ ψ a ˙ J a z J a x J a y J a r J a y ϕ a ˙ d ¯ a + τ a θ J a y d a θ θ a ˙ J a y , ψ a ¨ = ϕ a ˙ θ a ˙ J a x J a y J a z + τ a ψ J a z d a ψ ψ a ˙ J a z ,
where ϕ a , θ a , ψ a T is the attitude state, p a x , p a y , p a z T is the position state, τ a ϕ , τ a θ , τ a ψ are the three control torques, u a p is the control thrust, g a is the gravitational acceleration, m a denotes the mass, d ¯ a is the overall residual rotor angle, J a x , J a y , J a z represent the moments of the inertia, d a x , d a y , d a z , d a ϕ , d a θ , d a ψ represent the translational drag coefficients, Δ a x , Δ a y , Δ a z are the external disturbances encountered, and J a r denotes the moment of rotor’s inertia.
Inspired by the results in [25], the UAV’s rotational dynamics can be stabilized separately. Therefore, in light of external disturbances and parametric uncertainties, the UAV model (1) is redefined in the following manner
p ¨ a = g a u a + f a + Δ a
where p a = p a x , p a y , p a z T is the position, f a = [ d a x p ˙ a x / m a , d a y p ˙ a y / m a , d a z p ˙ a z / m a g a ] T , g a = diag 1 / m a , 1 / m a , 1 / m a , Δ a = Δ a x , Δ a y , Δ a z T , u a = u a x , u a y , u a z T denotes the new control signal which is given as [39]
u a x = sin ϕ a sin ψ a + cos ϕ a sin θ a cos ψ a u p u ay = sin ϕ a cos ψ a + cos ϕ a sin θ a sin ψ a u p u a z = cos θ a cos ϕ a u p
In the horizontal plane, the i-th i Π USV’s kinematic and dynamic equations are described as [40]
x ˙ b i = μ b i cos ψ b i v b i sin ψ b i , y ˙ b i = μ b i sin ψ b i + v b i cos ψ b i , ψ ˙ b i = r b i
μ ˙ b i = f μ b i α i + 1 m μ b i τ μ b i f + w μ b i , v ˙ b i = f v b i α i + 1 m v b i w v b i r ˙ b i = f r b i α i + 1 m r b i τ r b i f + w r b i
and
f μ b i α i = 1 m μ b i m v b i v b i r b i d μ b i μ b i d μ b i 1 μ b i μ b i f v b i α i = 1 m v b i m μ b i μ b i r b i d v b i v b i d v b i 1 v b i v b i , f r b i α i = 1 m r b i m μ b i m v b i μ b i v b i d r b i r b i d r b i 1 r b i r b i ,
where ψ b i denotes the yaw angle; x b i , y b i denotes the position; α i = μ b i , v b i , r b i T denote the surge, sway, and yaw velocity, respectively; m μ b i ,   m v b i ,   m r b i represent the inertial mass; f μ b i α i , f v b i α i , f r b i α i denote the nonlinear dynamics consisting of Coriolis forces and the unmodeled hydrodynamics; w μ b i ,   w v b i ,   w r b i represent the bounded disturbances; and τ μ b i f and τ r b i f represent the surge force and the yaw moment.
Because the USVs’ motion model given in (4) and (5) is underactuated, a hand position technique is employed to compensate. The USVs’ front point p b i x , p b i y is defined as the hand point, which can be expressed as
p b i x = x b i + L b i cos ψ b i p b i y = y b i + L b i sin ψ b i
where L b i denotes the distance between the new defined hand point p b i x , p b i y and the actual position x b i , y b i , which is depicted in Figure 1.
The second derivative of (7) yields the following result
p ¨ b i x = μ ˙ b i cos ψ b i v ˙ b i + L b i r ˙ b i sin ψ b i μ s i r b i sin ψ b i v b i r b i + L b i r b i 2 cos ψ b i p ¨ b i y = μ ˙ b i sin ψ b i v ˙ b i + L b i r ˙ b i cos ψ b i + μ s i r b i cos ψ b i v b i r b i + L b i r b i 2 sin ψ b i
Substituting (6) into (8) gives that
p ¨ b i x = f b i x ( β ) + cos ψ b i m μ b i τ μ f L b i sin ψ b i m r b i τ r f + w b i x p ¨ b i y = f b i y ( β ) + sin ψ b i m μ b i τ μ f + L b i cos ψ b i m r b i τ r f + w b i y
where
f b i x ( β ) = f μ ( α ) cos ψ b i f v ( α ) + L b i f r ( α ) sin ψ b i μ b i r b i sin ψ b i v b i r b i + L b i r b i 2 cos ψ b i f b i y ( β ) = f μ ( α ) sin ψ b i + f v ( α ) + L b i f r ( α ) cos ψ b i + μ b i r b i cos ψ b i v b i r b i + L b i r b i 2 sin ψ b i
w b i x ( β ) = w μ b i m μ b i cos ψ b i w v b i m μ b i + L b i w r b i m r b i sin ψ b i w b i y ( β ) = w μ b i m μ b i sin ψ b i + w v b i m μ b i + L b i w r b i m r b i cos ψ b i
with β = μ b i , v b i , r b i , ψ b i T .
Based on (9), the i-th USV position dynamics can be described as
p ¨ b i = f b i x y + Ω b i ψ b i ω b i u b i + w b i x y
where p b i = p b i x , p b i y T is the i-th USV position, f b i x y = f b i x , f b i y T , Ω b i ψ b i = cos ψ b i , sin ψ b i ; sin ψ b i , cos ψ b i , u b i = τ μ , τ r T , ω b i = diag 1 / m μ b i , L b i / m r b i , w b i x y = w b i x , w b i y T .
Combining (3) and (10), the unified model of the heterogeneous systems can be described as
x ˙ i 1 = x i 2 x ˙ i 2 = F x i + G x i u x i + Δ x i
When (11) represents the UAV model. Where x i 1 = p a x , p a y T = x a 1 R 2 , x i 2 = [ p ˙ a x , p ˙ a y ] T = x a 2 R 2 , F x i = F a x y = [ d x p ˙ a x / m a , d y p ˙ a y / m a ] T , G x i = G a x y = diag { 1 / m a , 1 / m a } , Δ x i = Δ a x y = [ Δ a x , Δ a y ] T , u x i = u a x y = u a x , u a y T .
When (11) represents the USV model. Where x i 1 = p b i x , p b i y T = x b i 1 R 2 , x i 2 = p ˙ b i x , p ˙ b i y T = x b i 2 R 2 , F x i = f b i x y = f b i x , f b i y T , G x i = Ω b i ψ b i ω b i , Δ x i = w b i x y = w b i x , w b i y T ,   u x i = u b i = τ μ , τ r T .
The virtual leader is defined and its motion model is described as follows
x ˙ l 1 = x l 2 x ˙ l 2 = F l
where x l 1 R 2 ,   x l 2 R 2 denote the position and velocity state vectors, respectively. F l R 2 is a smooth unknown nonlinear function.
Define the system’s error variable as
x ¯ i 1 = x i 1 x l 1 δ i x ¯ i 2 = x i 2 x l 2
where δ i R 2 represents the desired relative position vector between i-th of the agent of the heterogeneous systems and the virtual leader.
Assumption A1. 
The aerodynamic drag coefficients d i x , d i y , and d i z are bounded and unknown. The unknown nonlinear function F l is bounded. Then, there exists a positive constant α l such that F l < α l .
Assumption A2. 
The quadrotor UAV experiences external perturbations Δ a x , Δ a y , Δ a 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 μ b i , w ν b i , w r b i which are confined within certain limits, fulfilling the conditions w μ b i w ¯ μ b i , w v b i w ¯ v b i , w r b i w ¯ r b i . In this context, w ¯ μ b i , w ¯ v b i , and w ¯ r b i signify unknown positive constants.
Assumption A3. 
The velocity and acceleration of the agent and the virtual leader are bounded, and they satisfy v l 0 , v i 0 .
Assumption A4. 
G x i is a symmetric matrix and its eigenvalues λ 1 , λ 2 , , λ m are satisfied such that 0 < ϑ i < λ 1 < λ 2 < λ m < , where ϑ i is a positive constant.
Theorem 1
([41]). For any point in time, the function V ( t ) is continuous and positive, with its initial state being within limits. Given that the inequality V ˙ ( t ) γ V ( t ) + μ is satisfied, where γ > 0 and μ > 0 , it can be deduced that
V ( t ) μ γ 1 e γ t + V ( 0 ) e γ t .
Theorem 2
([42]). At all moments, the function S ( t ) remains positive and continuous, with its initial condition S ( 0 ) being constrained. Should the condition S ˙ ( t ) > q S ( t ) be valid for t t 0 0 with q > 0 , it leads to the inference that
S ( t ) > e q t t 0 S t 0 .

2.2. Algebraic Graph Theory

G = ( V g , ε g , A g ) is an undirected graph in this paper, with V g = v 1 , v 2 , , v n being the node-set, ε g V g × V g being the edge set, and A g = a i j being the adjacency matrix, of which all the elements are non-negative. The adjacency matrix A g depicts the effectiveness of communication from agent j to agent i, where a i j denotes the communication weight corresponding to the edge, and diagonal element a i i = 0 , a i j > 0 . As an undirected graph, the adjacency element a i j = a j i must be guaranteed to be satisfied. If any pair of unique nodes can be connected by an undirected path, then the graph is connected. The Laplacian matrix L g = l i j R n × n is defined as follows
L g = D g A g
where D g = diag d 1 , d 2 , , d n ,   d i = j = 1 n a i j and i = 1 , 2 , , n . Defining the leader adjacency weight matrix as B g = diag b 1 , b 2 , , b n , where b i > 0 , indicates an exchange of information between agent i and the leader, otherwise, b i = 0 , i = 1 , 2 , , n . It is presupposed that there is always at least one agent linked to the leader and b 1 + b 2 + + b n > 0 .
Theorem 3
([43]). G must be irreducible for the graph to be linked.
Theorem 4
([43]). L g ˜ = L g + B g = l 11 + b 1 l 1 n l n 1 l n n + b n are positive, where B g is the leader adjacency weight matrix if the Laplacian matrix L g of the undirected graph G is irreducible.

2.3. Improved Artificial Potential Field and Virtual Repulsion

The essential concept of the APF technique is that each agent is viewed as a high-potential field for the control objective of collision avoidance. Any agent that is close to the others will be repelled from their potential fields by the repulsive force. Specifically, the APF method has fewer design parameters and a significantly simpler controller architecture. Some results regard the agent as a particle, which makes collision avoidance issues impractical. In actuality, the domain of agents should be considered during an encounter situation. As illustrated in Figure 2, the domain of USV can be defined as three concentric circles with unequal radii.
As is shown in Figure 2, the repulsive force field is separated into the following four areas. The area outside the first circular area is defined as the Safe Zone, where d i j is greater than the detection distance r 3 of the potential field, there is no repulsive force in this area. The area inside the first circular area and outside the second circular is defined as the Negotiation Collision Avoidance Zone, where the distance d i j is larger than r 2 and less than or equal to r 3 . If the incoming ship is in zones A, B, or C in the area, the i-th USV is a giving-way vessel or has the same avoidance responsibility, and the i-th USV should take collision avoidance actions to avoid the collision. If the incoming ship is in zone D, the i-th USV is a stand-on vessel. The area inside the second circular area and outside the third circular area is defined as the Emergency Collision Avoidance Zone, where the distance d i j is larger than r 1 and less than or equal to r 2 . In the area where the distance between the other ship and the i-th USV is larger than and less than r 2 , the i-th USV should perform emergency collision avoidance actions that can momentarily violate the COLREGs. The area within the third circular area is defined as the Prohibited Zone, which is shown by the red circular limit with radius r 1 in Figure 2. All other USVs are prohibited from entering this area.
Where d i j = x i x j represents the relative position variable between the i-th USV and the incoming j-th USV, d i j represents the distance variable, and r 3 is the collision hazard detection distance, we conclude that the j-th USV is defined as a collision avoidance neighbor Π c of the i-th USV. Once its collision avoidance neighbor Π c holds d i j = x i x j 2 r 1 , it can be defined that the systems have been collided.
When the judgment condition 355 T r < 360 or 0 T r < 67.5 , r 2 < d i j r 3 are met, the incoming j-th USV is in Zone A and Zone B. Meanwhile, when the judgment condition 67.5 < T r 112.5 , r 2 < d i j r 3 are met, the incoming j-th USV is in Zone C. The i-th USV is a give-way ship. According to the COLREGs, it should turn right to give way. The repulsive force field is shown as follows [34]:
φ i j ( p , v ) = η d R j e θ m θ 1 1 d i j r 2 1 ρ 0
when the judgment condition r 1 < d i j r 2 is met, the i-th USV needs emergency collision avoidance. The repulsive force field is shown as follows:
φ i j ( p , v ) = η d R j 1 d i j r 1 1 ρ 0 2 + v o t cos θ 2
When d i j > r 2 or the incoming j-th USV is in Zone D, the repulsion field is not defined, that is, the repulsion is zero.
Where p and v are the position and speed, T r denotes the relative position of the incoming j-th USV and the i-th USV, η d is the positive scaling factor for the USVs, R j is the radius of the puffing circle of the other ship, ρ 0 denotes the repulsive potential field’s influence radius of other ships, θ m is the angle between the maximum relative position line, θ is the angle between p o t , v o t , p o t denotes the relative position line of the i-th USV and other ships, and v o t denotes the relative velocity.
The virtual repulsive force f i j c a is defined as the negative gradient of the repulsive potential function φ i j ( p , v ) as follows:
f i j c a = ( p , v ) φ i j ( p , v )
where ( p , v ) represents a negative gradient along p and v.
The total virtual repulsive force term of the whole systems is induced from (14)–(16) as follows:
u i c a = ω i ϑ i j Π c ( p , v ) φ i j ( p , v )
where ω i denotes the positive gain parameter and ϑ i is a positive constant.

3. Main Results

This section introduces a formation control protocol utilizing an extended state observer, aimed at guiding heterogeneous systems to follow a predefined trajectory in the XY plane, accounting for model uncertainties and external disturbances. Subsequently, a decentralized formation controller is developed to manage the height-tracking control of the UAV along the Z axis.

3.1. Controller Design Based on ESO

Rewrite (11) as
x ˙ i 1 = x i 2 x ˙ i 2 = F i + G x i u x i
where F i = F T i , F L i T = F x i + Δ x i , F T i and F L i denote the transverse and longitudinal orientation of F i and F 1 denotes the sum of the nonlinear unknown term and the external disturbance encountered of the UAV. F 2 , F 3 , F 4 , F 5 denote the sum of the nonlinear unknown term and the external disturbance encountered by different USVs, respectively.
For the purpose of approximating F i , the design of an extended state observer has been formulated as follows
x ^ ˙ i 1 = x ^ i 2 + α i 1 ε i k x i x ^ ˙ i 2 = G x i u x i + F ^ i + α i 2 ε i 2 k x i F ^ ˙ i = α i 3 ε i 3 k x i
where x ^ i 1 ,   x ^ i 2 ,   F ^ i represent the observer state, ε i > 0 , α i 1 R ,   α i 2 R ,   α i 3 R are positive constants, k x i = x i 1 x ^ i 1 . According to (19), we have
x ˜ ˙ i 1 = x ˜ i 2 + α i 1 ε i k x i x ˜ ˙ i 2 = F ˜ i + α i 2 ε i 2 k x i F ˜ ˙ i = α i 3 ε i 3 k x i F ˙ i
where x ˜ i 1 = x ^ i 1 x i 1 is the ESO estimation error.
Define η i as
η i = η i 1 η i 2 η i 3 T
where
η i 1 = x i 1 x ^ i 1 ε i 2 , η i 2 = x i 2 x ^ i 2 ε i , η i 3 = F i F ^ i
The observation error equation of state can be written as
ε i η i ˙ = A ¯ i η i + ε i B ¯ i F ˙ i
where
A ¯ i = α i 1 1 0 α i 2 0 1 α i 3 0 0 , B ¯ i = 0 0 1
The error in position and velocity for heterogeneous systems within the XY plane is defined as follows
e i x 1 = j Π a i j x 1 i x 1 j δ i δ j + b i x 1 i x 1 l δ i
e i x 2 = j Π a i j x 2 i x 2 j + b i x 2 i x 2 l
The formation control term is designed as follows
u i f = 1 ϑ i F ^ i k i e i x 1 + e i x 2
Furthermore, the distributed formation controller with collision avoidance performance can be defined as follows
u i = u i f + u i c a
Remark 1. 
In the formation controller (25), the distributed formation control term u i f is used by tracking the time-varying velocity variable of the leader. The virtual repulsive force term u i c a is used to achieve the collision avoidance objective among the agents. By designing the action functions (14) and (15) of the APF, they are non-conflicting items.
Using (17) and (24), distributed formation controller (25) can be given as
u i = 1 ϑ i F ^ i k i e i x 1 + e i x 2 ω i ϑ i j Π c ( p , v ) φ i j ( p , v )
Take the derivative of (13) and substitute (11) and (12) into it. Define the derivatives of the error variables and write them in vector form as follows
Z ¯ ˙ = M I m Z ¯ + 0 n m F x i + 0 n m U + 0 n m Δ x i 0 n m F l
where Z ¯ ˙ = x ¯ ˙ 11 T , , x ¯ ˙ n 1 T , x ¯ ˙ 12 T , , x ¯ ˙ n 2 T T , U = G x i u x i , M = 0 n I n 0 n 0 n , ⊗ denotes the Kronecker product.

3.2. Altitude Controller Design for UAV

The error system of altitude is defined as
e z p = p a z c p
e z v = v a z c v
where the desired position signal is c p and the desired velocity signal is c v .
The adaptive laws and control input of formation tracking are constructed as
u a z = H ^ a z 1 σ e ζ e ζ κ ^ 11 ω 1 T ω 1 2 ς 1 2 e ζ κ ^ 12 2 ς 2 2 + c ˙ v k ζ e z v
κ ^ ˙ 11 = l 11 k 11 κ ^ 11 + e ζ T e ζ ω 1 T ω 1 2 ς 1 2
κ ^ ˙ 12 = l 12 k 12 κ ^ 12 + e ζ T e ζ 2 ς 2 2
H ^ ˙ a z = Proj H ̲ a z , H ¯ a z { F } = 0 , if H ^ a z = H ¯ a z and F 0 or H ^ a z = H ̲ a z and F 0 F , otherwise
where F = l 13 k 13 H ^ a z + e ζ u a z ,   H ̲ a z , and H ¯ a z are the parameter H a z ’s lower bound and upper bound of, respectively, where σ , k 11 ,   k 12 ,   k 13 ,   k ζ ,   l 11 ,   l 12 ,   l 13 ,   ς 1 ,   ς 2 denote positive parameters.

3.3. Stability Analysis

  • Part A. Proof of the stability of the extended state observer
Consider the following Lyapunov function
V 0 = i = 1 N + 1 V 0 i
where V 0 i is given as
V 0 i = ε i η i T P i η i
V ˙ 0 i = ε i η i ˙ T P i η i + ε i η i T P i η i ˙ = ( A ¯ i η i + ε i B ¯ i F ˙ i ) T P i η i + η i T P i ( A ¯ i η i + ε i B ¯ i F ˙ n i ) = η i T A ¯ i T P i η i + ε i ( B ¯ i F ˙ n i ) T P η i + η i T P i A ¯ i η i + ε η i T P i B ¯ i F ˙ n i = η i T A ¯ i T P i + P i A ¯ i η i + 2 ε i η i T P i B ¯ i F ˙ n i η i T Q i η i + 2 ε i P i B ¯ i · η i · | F ˙ i |
and
V ˙ 0 i λ min ( Q i ) η i 2 + 2 ε i | F ˙ i | P i B ¯ i η i
Therefore, the convergence condition of the observer satisfies
η i 2 ε i | F ˙ i | P B ¯ i λ min ( Q i )
where Q i denotes any given symmetric positive definite matrix, λ min ( Q i ) is the minimum eigenvalue of Q i , and there is a three-by-three symmetric positive definite matrix P i that satisfies the following equation
A ¯ i T P i + P i A ¯ i + Q i = 0
Remark 2. 
Based on Assumptions 2 and 3, both the velocity and acceleration of the agent, as well as the external disturbances, are bounded. From F i = F x i + Δ x i , it is evident that F i encompasses the agent’s velocity and acceleration, as well as external disturbances. So, | F ˙ i | is bounded.
  • Part B. Proof of the stability of the system
Define the Lyapunov quadratic scalar function as follows
V 1 = 1 2 Z ¯ T P g I m Z ¯
where P = 2 L ˜ g L ˜ g L ˜ g L ˜ g , L ˜ g = L g + B g .
According to Theorem 4, the Lyapunov quadratic scalar function V 1 is a positive definite function.
Substitute (27) into V 1 ˙ and expand the partial matrix as follows
V 1 ˙ = Z ¯ T P g I m Z ¯ ˙ = 1 2 Z ¯ T M T P g + P g M I m Z ¯ + Z ¯ T P g I m 0 n m F x i + U + Δ x i F l = Z ¯ T 0 n L ˜ g L ˜ g L ˜ g I m Z ¯ + i = 1 n e i x 1 + e i x 2 T F x i + G x i u x i + Δ x i F l
Substituting the distributed formation controller (24) into (41), we can obtain
V 1 ˙ = Z ¯ T 0 n L ˜ g L ˜ g L ˜ g I m Z ¯ + i = 1 n e i x 1 + e i x 2 T F x i i = 1 n e i x 1 + e i x 2 T F l + i = 1 n e i x 1 + e i x 2 T × F ^ i k i e i x 1 + e i x 2 ω i j Π c ( p , v ) φ i j ( p , v ) = Z ¯ T 0 n L ˜ g L ˜ g L ˜ g I m Z ¯ i = 1 n e i x 1 + e i x 2 T F l i = 1 n k i e i x 1 + e i x 2 2 i = 1 n e i x 1 + e i x 2 T × ω i j Π c ( p , v ) φ i j ( p , v )
In the case of an undirected graph, the Laplacian matrix exhibits symmetry, and the virtual repulsive force term u i c a becomes negligible once collision avoidance is achieved through the repulsive potential function φ i j ( p , v ) .
Therefore, the aforementioned term V = i = 1 n e i x 1 + e i x 2 T ω i j Π c ( p , v ) φ i j ( p , v ) = 0 , based on Assumption 1 and Young’s inequality as follows
e i x 1 + e i x 2 T F l ξ i e i x 1 + e i x 2 2 + α l 2 4 ξ i
where ξ i is a positive constant.
And then the results (42) can be derived as follows
V 1 ˙ Z ¯ T 0 n L ˜ g L ˜ g L ˜ g I m Z ¯ + i = 1 n ξ i e i x 1 + e i x 2 2 + i = 1 n α l 4 ξ i i = 1 n k i e i x 1 + e i x 2 2 = Z ¯ T 0 n L ˜ g L ˜ g L ˜ g I m Z ¯ + i = 1 n ξ i k i e i x 1 + e i x 2 2 + i = 1 n α l 2 4 ξ i
Let ρ i satisfy the limit of ρ i k i ξ i and ρ i > 0 , the inequality is transformed into
V 1 ˙ Z ¯ T 0 n L ˜ g L ˜ g L ˜ g I m Z ¯ i = 1 n ρ i e i x 1 + e i x 2 2 + i = 1 n α l 2 4 ξ i = Z ¯ T ρ i L ˜ g 2 L ˜ g 2 L ˜ g 2 L ˜ g 2 0 n L ˜ g L ˜ g L ˜ g I m Z ¯ + i = 1 n α l 2 4 ξ i = Z ¯ T ρ i L ˜ g 2 ρ i L ˜ g 2 L ˜ g ρ i L ˜ g 2 L ˜ g ρ i L ˜ g 2 L ˜ g I m Z ¯ + i = 1 n α l 2 4 ξ i
According to the Lemma of linear matrix inequality in [41], ρ i L ˜ g 2 ρ i L ˜ g 2 L ˜ g = L ˜ g 2 > 0 , ρ i L ˜ g 2 L ˜ g > 0 , so the matrix ρ i L ˜ g 2 ρ i L ˜ g 2 L ˜ g ρ i L ˜ g 2 L ˜ g ρ i L ˜ g 2 L ˜ g > 0 . Therefore, the inequality (44) can be rewritten as follows
V 1 ˙ Z ¯ T ρ i L ˜ g 2 L ˜ g 2 L ˜ g 2 L ˜ g 2 0 n L ˜ g L ˜ g L ˜ g I m Z ¯ + i = 1 n α l 2 4 ξ i = Z ¯ T ρ i Θ Ξ I m Z ¯ + Δ
where
Θ = L ˜ g 2 L ˜ g 2 L ˜ g 2 L ˜ g 2 ,   Ξ = 0 n L ˜ g L ˜ g L ˜ g ,   Δ = i = 1 n α l 2 4 ξ i .
Set ρ i > 1 λ min Θ λ max Ξ + 2 λ max P , where λ min Θ , λ max Ξ , λ max P denote the smallest eigenvalue of the matrix Θ , the largest eigenvalue of the matrix Ξ , and the largest eigenvalue of the matrix P, respectively, and = min σ 1 Y 1 , , σ n Y n .
Then, we can obtain
V 1 ˙ Z ¯ T P I m Z ¯ + Δ = V 1 + Δ
In accordance with Theorem 1, the subsequent inequality can be presented in the following manner
V 1 ( t ) V 1 ( 0 ) e t + Δ 1 e t
Consequently, by choosing suitable parameters as indicated in [44], one can achieve the performance of adaptive leader-following formation.
  • Part C. Proof of collision avoidance
The collision avoidance performance is analyzed only for the i-th ( i Π ) USV and its collision avoidance neighbor the j-th ( j Π c ) USV, and the others can be analyzed by the same way in [42].
Define the energy function holds quadratic form as follows
S ( t ) = 1 2 v i T ( t ) v i ( t ) + 1 2 d i k T ( t ) d i k ( t )
Taking the time derivative of (48), we can get
S ˙ ( t ) = ω i v i T ( t ) ( p , v ) φ i j ( p , v ) + d i k T ( t ) v i ( t ) v k ( t ) v i T ( t ) F ^ i k i v i T ( t ) e i x 1 ( t ) + e i x 2 ( t )
where d i k ( t ) is the relative position variable.
Because of the designing of the repulsive potential function φ i j ( p , v ) , φ i d i k ( t ) will approach infinity by designing the parameter ω ¯ i > 0 if the i-th USV is closing to the j-th USV. Therefore, if the i-th USV and the j-th USV are close to each other, and design the appropriate parameter, the inequality can be met in the following manner:
v i T ( t ) ( p , v ) φ i j ( p , v ) > 1 2 d i k T ( t ) d i k ( t ) + 1 2 v i T ( t ) v i ( t ) 1 ω i d i k T ( t ) v i ( t ) v k ( t ) + 1 ω i v i T ( t ) F ^ i + k i ω i v i T ( t ) e i x 1 ( t ) + e i x 2 ( t )
Substituting the aforementioned inequality (50) into (49), the subsequent inequality can be derived
S ˙ ( t ) > ω i S ( t )
Based on the principles outlined in Theorem 2, the subsequent inequality can be derived
d i k T ( t ) d i k ( t ) > 2 e ω i t t 0 S ( t ) v i T ( t ) v i ( t )
Since the previously indicated term v i T ( t ) v i ( t ) is continuous and bounded at this dwell time. By designing the gain parameter ω i > 0 large enough, we get can the inequality 2 e ω i t t 0 S ( t ) v i T ( t ) v i ( t ) > ( 2 r 1 ) 2 and the result can be obtained as d i k ( t ) > 2 r 1 . As a result, the collision avoidance performance can be guaranteed by the proposed formation control with collision avoidance strategy.
  • Part D. Proof of the decentralized formation controller
If the feasibility condition c v c ˙ p = 0 is satisfied, it then allows the time derivative of Equations (28) and (29) to be expressed as follows
e ˙ z p = e z v
e ˙ z v = F a z + G a z u a z + Δ a z c ˙ v .
where e z v may be regarded as the virtual control input in the system (51). The stability of the system (51) can be assured by constructing the virtual control input ζ = k ζ e z p .
A positive Lyapunov function is considered as
V z p = 1 2 e z p T e z p
A new error is defined as
e ζ = e z v ζ
By computing the time derivative of Equation (53) and incorporating the control input from Equation (30) into this calculation, one can obtain
e ˙ ζ = σ e ζ + F a z + Δ a z + H ˜ a z u a z e ζ κ ^ 11 ω 1 T ω 1 2 ς 1 2 e ζ κ ^ 12 2 ς 2 2 .
V z is given as
V z = 1 2 e z p T e z p + 1 2 e ζ T e ζ + κ ˜ 11 2 2 l 11 + κ ˜ 12 2 2 l 12 + Tr H ˜ a z T H ˜ a z 2 l 13
with κ ˜ 11 = κ 11 κ ^ 11 ,   κ ˜ 12 = κ 12 κ ^ 12 and H ˜ a z = H a z H ^ a z .
V ˙ z k ζ e z p T e z p σ e ζ T e ζ + e ζ T F a z + e ζ T Δ a z + e ζ T H ˜ a z u a z e ζ T e ζ κ ^ 11 ω 1 T ω 1 2 ς 1 2 e ζ T e ζ κ ^ 12 2 ς 2 2 + κ ˜ 11 κ ˜ ˙ 11 l 11 + κ ˜ 12 κ ˜ ˙ 12 l 12 + Tr H ˜ a z T H ˜ ˙ a z l 13
Given that F a z is an undefined function, as per the universal approximation theorem cited in [45], for any arbitrarily small constant ϵ 1 , there exists a fuzzy logic system θ 1 * T ω 1 that can be represented as F a z = θ 1 * T ω 1 + ϵ 1 . In this context, θ 1 * T represents the optimal weight, ω 1 is the fuzzy basis vector, ϵ 1 denotes the fuzzy system’s approximation error, and it is established that ϵ 1 ϵ ¯ 1 .
Regarding Young’s inequality, it can be deduced that
e ζ T θ 1 * T ω 1 e ζ T e ζ κ 11 ω 1 T ω 1 2 ς 1 2 + ς 1 2 2
e ζ T ϵ 1 + Δ a z e ζ T e ζ κ 12 2 ς 2 2 + ς 2 2 2
where κ 11 = θ 1 * T θ 1 * and κ 12 = ϵ ¯ 1 + Δ ¯ a z T ϵ ¯ 1 + Δ ¯ a z .
Substituting (58) and (59) into (57), we can obtain
V ˙ i z k ζ e z p T e z p σ e ζ T e ζ + e ζ T e ζ κ 11 ω 1 T ω 1 2 ς 1 2 + e ζ T e ζ κ 12 2 ς 2 2 e ζ T e ζ κ ^ 11 ω 1 T ω 1 2 ς 1 2 e ζ T e ζ κ ^ 12 2 ς 2 2 κ ˜ 12 k 12 κ ^ 12 + e ζ T e ζ 2 ς 2 2 κ ˜ 11 k 11 κ ^ 11 + e ζ T e ζ ω 1 T ω 1 2 ς 1 2 + e ζ T H ˜ a z u a z Tr H ˜ a z T k 13 H ^ a z + e ζ u a z + ς 1 2 2 + ς 2 2 2 k ζ e z p T e z p σ e ζ T e ζ k 11 2 κ ˜ 11 2 k 12 2 κ ˜ 12 2 k 13 2 H ˜ a z F 2 + k 11 2 κ 11 2 + k 12 2 κ 12 2 + k 13 2 H a z F 2 + ς 1 2 2 + ς 2 2 2 ϑ z V z + v z
where θ z = min 2 k ζ , 2 σ , l 11 k 11 , l 12 k 12 , l 13 k 13 > 0 , v z = k 11 2 κ 11 2 + k 12 2 κ 22 2 + k 13 2 H a z F 2 + ς 1 2 2 + ς 2 2 2 .
Based on Equation (60), V ˙ z is derived as follows
V ˙ z ϑ z V z + v z
According to the boundedness theorem, the closed-loop system solution is uniformly eventually bounded. Using e z p and e ζ , it is possible to deduce that the UAV’s velocity error and altitude tracking error are uniformly eventually bounded.

4. Simulation Result

In this section, a simulation analysis for the heterogeneous systems with collision avoidance performance under parameter uncertainty and external disturbance is provided to illustrate the effect of the proposed approach. One UAV and four USVs make up the heterogeneous system. The simulation experiments described in this section were executed on the MATLAB R2020a simulation platform.
The communication topology is defined in Figure 3, where agent l represents the virtual leader, agent 0 represents the UAV, and agents 1–4 represent the four USVs. If B g = diag { 0 , 0 , 0 , 0 , 1 } is the adjacency weight matrix between agents and the virtual leader, then the adjacency matrix A g and the Laplacian matrix L g are as follows:
A g = 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 g = 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
The system parameters of a UAV and USVs are given in Table 1 and Table 2, respectively. APF design parameters r i n = 0.5 , r out = 2 , ω ¯ i = 0.3 , ω i = 1 , where i = 1 , , 4 . The external disturbances are given as Δ a x y = [ 0.2 cos ( 0.5 t ) , 0.8 cos ( t ) ] T , Δ a z = 0.3 cos ( 0.2 t ) , w b i x y = [ 1.1 cos ( 0.5 t ) , 0.2 sin ( 2 t ) ] T . The total duration of the simulation runs T z = 80   s and the sampling time T s = 0.6   s .
The initial position state vector of the UAV and USVs are defined as x 0 ( 0 ) = [ 6 , 2 , 10 ] T , x 1 ( 0 ) = [ 8 , 0 , 0 ] T ,   x 2 ( 0 ) = [ 4 , 0 , 0 ] T , x 3 ( 0 ) = [ 4 , 4 , 0 ] T ,   x 4 ( 0 ) = [ 8 , 4 , 0 ] T ,   δ 0 = [ 0 , 0 ] T ,   δ 1 = [ 0 , 3 ] T ,   δ 2 = [ 3 , 0 ] T ,   δ 3 = [ 3 , 0 ] T ,   δ 4 = [ 0 , ] T ,   = 3 2 .
The simulation results by using the proposed control method for heterogeneous systems with collision avoidance strategy are shown in Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11. The trajectories of the four USVs and the UAV in the 3D environment are depicted in Figure 4. The black line represents the UAV’s trajectory, while the other color lines represent the trajectories of the four USVs. From the given starting point, the USVs maintain a safe distance and establish a designated configuration between each other. Figure 5 shows the trajectories of heterogeneous systems without collision avoidance. At the sampling time T s = 0.6   s , Figure 5 shows that three USVs have been collided. Simulation trajectory plots demonstrate the importance of collision avoidance for heterogeneous systems.
Figure 6 and Figure 7 represent the position tracking errors and velocity errors in x-label and y-label, respectively. Under the influence of external disturbances and parameter uncertainties, the formation trajectory tracking errors are capable of converging to a minimal residual set. Figure 8 and Figure 9 depict the surge force and yaw torque. The outcomes of the simulations indicate that the USVs are able to continue tracking the virtual leader even after activating the improved APF for collision avoidance.
Figure 10 and Figure 11 further display the distances between USVs without and with considering the collision avoidance. Figure 10 shows that the USVs have collided, because of the distance d i j min < 2 r 1 = 1 . Figure 11 shows that the distances of USVs can always maintain d i j min > 2 r 1 = 1 at any time. The simulation diagrams further confirm that the proposed formation control with collision avoidance strategy has an excellent performance. The performance of ESO in estimating F i is depicted in Figure 12, which clearly demonstrates the effectiveness of the proposed solution for the accurate compensation of uncertainties and external disturbances within heterogeneous systems.

5. Conclusions

In this paper, the cooperative formation trajectory tracking problem for heterogeneous systems with collision avoidance performance under parameter uncertainty and external disturbance is investigated. Under the unified dynamic model of a UAV and USVs in the XY plane, a formation control protocol based on ESO is proposed. The collision avoidance control approach for USV formation is then constructed using APF theory. Furthermore, the APF approach is improved by rebuilding the repulsive potential field to bring the unmanned ship’s collision avoidance behavior more in compliance with COLREGs. According to the Lyapunov stability theory, the time-varying formation errors and tracking errors are uniformly ultimately limited. Simulation results verify the effectiveness of the proposed scheme.
This paper focuses solely on heterogeneous formation and collision avoidance between a single UAV and multiple USVs. Going forward, cooperative collision avoidance involving multiple UAVs and USVs within intricate environments deserves further study. In discussions of formation control issues for multi-agent or multi-vehicle systems, it is common to assume that the controlled entities operate within an ideal network communication environment. However, in practical scenarios, such as with USVs, the communication environment is often fraught with issues like noise, time delays, and packet loss. Consequently, the problem of cooperative collision avoidance control for swarm systems under adverse communication conditions also merits thorough investigation.

Author Contributions

Conceptualization, Y.H. and J.N.; methodology, Y.H. and J.N.; software, Y.H.; validation, Y.H.; writing—original draft preparation, Y.H.; writing—review and editing, Y.H., W.L., J.N. and Z.L.; supervision, W.L. and J.N.; project administration, W.L. and Z.L.; funding acquisition, W.L. 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 52271304,and 61976033, and in 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 conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
UAVUnmanned Aerial Vehicle
UGVUnmanned Ground Vehicle
USVUnmanned Surface Vehicle
AUVAutonomous Underwater Vehicle
ESOExtended State Observer
APFArtificial Potential Field
COLREGInternational Regulations for Preventing Collisions at Sea
RBFRadial Basis Function

References

  1. Yan, X.; Jiang, D.; Miao, R.; Li, Y. Formation control and obstacle avoidance algorithm of a multi-USV system based on virtual structure and artificial potential field. J. Mar. Sci. Eng. 2021, 9, 161. [Google Scholar] [CrossRef]
  2. Shan, Q.; Wang, X.; Li, T.; Chen, C.P. Finite-time control for USV path tracking under input saturation with random disturbances. Appl. Ocean Res. 2023, 138, 103628. [Google Scholar] [CrossRef]
  3. Guerrero-Ibañez, J.; Contreras-Castillo, J.; Zeadally, S. Deep learning support for intelligent transportation systems. Trans. Emerg. Telecommun. Technol. 2021, 32, e4169. [Google Scholar] [CrossRef]
  4. Li, X.; Tan, J.; Liu, A.; Vijayakumar, P.; Kumar, N.; Alazab, M. A Novel UAV-Enabled Data Collection Scheme for Intelligent Transportation System Through UAV Speed Control. IEEE Trans. Intell. Transp. Syst. 2021, 22, 2100–2110. [Google Scholar] [CrossRef]
  5. Menouar, H.; Guvenc, I.; Akkaya, K.; Uluagac, A.S.; Kadri, A.; Tuncer, A. UAV-Enabled Intelligent Transportation Systems for the Smart City: Applications and Challenges. IEEE Commun. Mag. 2017, 55, 22–28. [Google Scholar] [CrossRef]
  6. Sun, Y.; Zhang, D.; Wang, Y.; Zong, Z.; Wu, Z. Model Experimental Study on a T-Foil Control Method with Anti-Vertical Motion Optimization of the Mono Hull. J. Mar. Sci. Eng. 2023, 11, 1551. [Google Scholar] [CrossRef]
  7. Ke, C.; Chen, H. Cooperative path planning for air–sea heterogeneous unmanned vehicles using search-and-tracking mission. Ocean Eng. 2022, 262, 112020. [Google Scholar] [CrossRef]
  8. Ren, Y.; Zhang, L.; Ying, Y.; Li, S.; Tang, Y. Model-Parameter-Free Prescribed Time Trajectory Tracking Control for Under-Actuated Unmanned Surface Vehicles with Saturation Constraints and External Disturbances. J. Mar. Sci. Eng. 2023, 11, 1717. [Google Scholar] [CrossRef]
  9. Peng, Z.; Wang, J.; Wang, D. Distributed maneuvering of autonomous surface vehicles based on neurodynamic optimization and fuzzy approximation. IEEE Trans. Control Syst. Technol. 2017, 26, 1083–1090. [Google Scholar] [CrossRef]
  10. Zhou, Z.; Li, M.; Hao, Y. A Novel Region-Construction Method for Multi-USV Cooperative Target Allocation in Air–Ocean Integrated Environments. J. Mar. Sci. Eng. 2023, 11, 1369. [Google Scholar] [CrossRef]
  11. Fu, H.; Yao, W.; Cajo, R.; Zhao, S. Trajectory Tracking Predictive Control for Unmanned Surface Vehicles with Improved Nonlinear Disturbance Observer. J. Mar. Sci. Eng. 2023, 11, 1874. [Google Scholar] [CrossRef]
  12. Gu, N.; Wang, D.; Peng, Z.; Liu, L. Distributed containment maneuvering of uncertain under-actuated unmanned surface vehicles guided by multiple virtual leaders with a formation. Ocean Eng. 2019, 187, 105996. [Google Scholar] [CrossRef]
  13. Tan, G.; Zhuang, J.; Zou, J.; Wan, L. Multi-type task allocation for multiple heterogeneous unmanned surface vehicles (USVs) based on the self-organizing map. Appl. Ocean Res. 2022, 126, 103262. [Google Scholar] [CrossRef]
  14. Peng, Z.; Wang, D.; Chen, Z.; Hu, X.; Lan, W. Adaptive dynamic surface control for formations of autonomous surface vehicles with uncertain dynamics. IEEE Trans. Control Syst. Technol. 2012, 21, 513–520. [Google Scholar] [CrossRef]
  15. Miao, R.; Wang, L.; Pang, S. Coordination of distributed unmanned surface vehicles via model-based reinforcement learning methods. Appl. Ocean Res. 2022, 122, 103106. [Google Scholar] [CrossRef]
  16. 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. 2022, 10, 564–574. [Google Scholar] [CrossRef]
  17. Peng, Z.; Wang, J.; Wang, D.; Han, Q.L. An overview of recent advances in coordinated control of multiple autonomous surface vehicles. IEEE Trans. Ind. Inform. 2020, 17, 732–745. [Google Scholar] [CrossRef]
  18. 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]
  19. Liu, W.; Ye, H.; Yang, X. Model-Free Adaptive Sliding Mode Control Method for Unmanned Surface Vehicle Course Control. J. Mar. Sci. Eng. 2023, 11, 1904. [Google Scholar] [CrossRef]
  20. 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]
  21. 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]
  22. 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]
  23. 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]
  24. Wang, N.; Ahn, C.K. Coordinated Trajectory-Tracking Control of a Marine Aerial-Surface Heterogeneous System. IEEE/ASME Trans. Mechatronics 2021, 26, 3198–3210. [Google Scholar] [CrossRef]
  25. 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]
  26. 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]
  27. Li, L.; Wu, D.; Huang, Y.; Yuan, Z.M. A path planning strategy unified with a COLREGS collision avoidance function based on deep reinforcement learning and artificial potential field. Appl. Ocean Res. 2021, 113, 102759. [Google Scholar] [CrossRef]
  28. Sun, X.; Wang, G.; Fan, Y. Collision avoidance guidance and control scheme for vector propulsion unmanned surface vehicle with disturbance. Appl. Ocean Res. 2021, 115, 102799. [Google Scholar] [CrossRef]
  29. Ghommam, J.; Saad, M.; Mnif, F.; Zhu, Q.M. Guaranteed Performance Design for Formation Tracking and Collision Avoidance of Multiple USVs With Disturbances and Unmodeled Dynamics. IEEE Syst. J. 2021, 15, 4346–4357. [Google Scholar] [CrossRef]
  30. Dai, S.L.; He, S.; Lin, H.; Wang, C. Platoon Formation Control With Prescribed Performance Guarantees for USVs. IEEE Trans. Ind. Electron. 2018, 65, 4237–4246. [Google Scholar] [CrossRef]
  31. Peng, Z.; Wang, D.; Li, T.; Han, M. Output-Feedback Cooperative Formation Maneuvering of Autonomous Surface Vehicles With Connectivity Preservation and Collision Avoidance. IEEE Trans. Cybern. 2020, 50, 2527–2535. [Google Scholar] [CrossRef] [PubMed]
  32. 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]
  33. Xu, X.; Pan, W.; Huang, Y.; Zhang, W. Dynamic Collision Avoidance Algorithm for Unmanned Surface Vehicles via Layered Artificial Potential Field with Collision Cone. J. Navig. 2020, 73, 1306–1325. [Google Scholar] [CrossRef]
  34. Lyu, H.; Yin, Y. COLREGS-constrained real-time path planning for autonomous ships using modified artificial potential fields. J. Navig. 2019, 72, 588–608. [Google Scholar] [CrossRef]
  35. Song, J.; Hao, C.; Su, J. Path planning for unmanned surface vehicle based on predictive artificial potential field. Int. J. Adv. Robot. Syst. 2020, 17, 172988142091846. [Google Scholar] [CrossRef]
  36. Tan, G.; Zhuang, J.; Zou, J.; Wan, L.; Sun, Z. Artificial potential field-based swarm finding of the unmanned surface vehicles in the dynamic ocean environment. Int. J. Adv. Robot. Syst. 2020, 17, 172988142092530. [Google Scholar] [CrossRef]
  37. Sang, H.; You, Y.; Sun, X.; Zhou, Y.; Liu, F. The hybrid path planning algorithm based on improved A* and artificial potential field for unmanned surface vehicle formations. Ocean Eng. 2021, 223, 108709. [Google Scholar] [CrossRef]
  38. 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]
  39. Song, Y.; He, L.; Zhang, D.; Qian, J.; Fu, J. Neuroadaptive Fault-Tolerant Control of Quadrotor UAVs: A More Affordable Solution. IEEE Trans. Neural Netw. Learn. Syst. 2019, 30, 1975–1983. [Google Scholar] [CrossRef]
  40. 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]
  41. Wen, G.; Chen, C.L.P.; Liu, Y.J. Formation Control With Obstacle Avoidance for a Class of Stochastic Multiagent Systems. IEEE Trans. Ind. Electron. 2018, 65, 5847–5855. [Google Scholar] [CrossRef]
  42. Wen, G.; Chen, C.L.P.; Dou, H.; Yang, H.; Liu, C. Formation control with obstacle avoidance of second-order multi-agent systems under directed communication topology. Sci. China Inf. Sci. 2019, 62, 1–14. [Google Scholar] [CrossRef]
  43. Shi, Q.; Li, T.; Li, J.; Chen, C.P.; Xiao, Y.; Shan, Q. Adaptive leader-following formation control with collision avoidance for a class of second-order nonlinear multi-agent systems. Neurocomputing 2019, 350, 282–290. [Google Scholar] [CrossRef]
  44. Wen, G.; Chen, C.L.P.; Feng, J.; Zhou, N. Optimized Multi-Agent Formation Control Based on an Identifier–Actor–Critic Reinforcement Learning Algorithm. IEEE Trans. Fuzzy Syst. 2018, 26, 2719–2731. [Google Scholar] [CrossRef]
  45. Kim, Y.H.; Ahn, S.C.; Kwon, W.H. Computational complexity of general fuzzy logic control and its simplification for a loop controller. Fuzzy Sets Syst. 2000, 111, 215–224. [Google Scholar] [CrossRef]
Figure 1. The kinematic model of USV.
Figure 1. The kinematic model of USV.
Jmse 11 02332 g001
Figure 2. Repulsive potential field partition.
Figure 2. Repulsive potential field partition.
Jmse 11 02332 g002
Figure 3. Communication topology graph.
Figure 3. Communication topology graph.
Jmse 11 02332 g003
Figure 4. Formation evolution of heterogeneous systems with collision avoidance.
Figure 4. Formation evolution of heterogeneous systems with collision avoidance.
Jmse 11 02332 g004
Figure 5. Formation evolution of heterogeneous systems without collision avoidance.
Figure 5. Formation evolution of heterogeneous systems without collision avoidance.
Jmse 11 02332 g005
Figure 6. Position tracking errors of the UAV and USVs with collision avoidance.
Figure 6. Position tracking errors of the UAV and USVs with collision avoidance.
Jmse 11 02332 g006
Figure 7. Velocity errors of the UAV and USVs with collision avoidance.
Figure 7. Velocity errors of the UAV and USVs with collision avoidance.
Jmse 11 02332 g007
Figure 8. The surge force of heterogeneous system.
Figure 8. The surge force of heterogeneous system.
Jmse 11 02332 g008
Figure 9. The yaw torque of heterogeneous system.
Figure 9. The yaw torque of heterogeneous system.
Jmse 11 02332 g009
Figure 10. The distance between the USVs without collision avoidance.
Figure 10. The distance between the USVs without collision avoidance.
Jmse 11 02332 g010
Figure 11. The distance between the USVs with collision avoidance.
Figure 11. The distance between the USVs with collision avoidance.
Jmse 11 02332 g011
Figure 12. The estimated performance of ESO.
Figure 12. The estimated performance of ESO.
Jmse 11 02332 g012
Table 1. The model parameters of UAV.
Table 1. The model parameters of UAV.
ParameterValueUnit
m a 2kg
g a 9.8 m · s 2
J a x , J a y , J a z 1.5 N · s 2 · rad 1
d a x , d a y , d a z 0.012 N · s 2 · rad 1
Table 2. The model parameters of USV.
Table 2. The model parameters of USV.
ParameterValueUnit
m μ b i 25.8kg
m v b i 33.8kg
m r b i 2.76kg
d μ b i 0.725kg · s 1
d v b i 0.89kg · s 1
d r b i 1.9 kg · m 2 · s 1
d μ b i 1 1.33 kg · s 1
d v b i 1 36.47 kg · s 1
d r b i 1 0.75 kg · m 2 · s 1
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

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. https://doi.org/10.3390/jmse11122332

AMA Style

Huang Y, Li W, Ning J, Li Z. Formation Control for UAV-USVs Heterogeneous System with Collision Avoidance Performance. Journal of Marine Science and Engineering. 2023; 11(12):2332. https://doi.org/10.3390/jmse11122332

Chicago/Turabian Style

Huang, Yuyang, Wei Li, Jun Ning, and Zhihui Li. 2023. "Formation Control for UAV-USVs Heterogeneous System with Collision Avoidance Performance" Journal of Marine Science and Engineering 11, no. 12: 2332. https://doi.org/10.3390/jmse11122332

APA Style

Huang, Y., Li, W., Ning, J., & Li, Z. (2023). Formation Control for UAV-USVs Heterogeneous System with Collision Avoidance Performance. Journal of Marine Science and Engineering, 11(12), 2332. https://doi.org/10.3390/jmse11122332

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