Next Article in Journal
A Numerical Study of Geomorphic and Oceanographic Controls on Wave-Driven Runup on Fringing Reefs with Shore-Normal Channels
Previous Article in Journal
Improved RRT Algorithm for AUV Target Search in Unknown 3D Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Swarm Control for Connectivity-Preserving and Collision-Avoiding Unmanned Surface Vehicles Subject to Multiple Constraints

College of Intelligent Systems Science and Engineering, Harbin Engineering University, Harbin 150001, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2022, 10(6), 827; https://doi.org/10.3390/jmse10060827
Submission received: 4 May 2022 / Revised: 15 June 2022 / Accepted: 15 June 2022 / Published: 17 June 2022
(This article belongs to the Section Ocean Engineering)

Abstract

:
This paper investigates swarm control for unmanned surface vessels subject to multiple constraints. These constraints can be summarized as model parameter uncertainty, the unavailability of velocity measurements, time-varying environmental disturbances, input saturation and output constraints. Firstly, to recover unmeasured velocity information, to identify unknown vehicle dynamics and to estimate time-varying environmental disturbances, a neural adaptive state observer is designed for each vessel. Secondly, to avoid complex calculations, a second-order linear tracking differentiator is employed to generate a smooth reference signal and to extract the time derivative of the kinematic control law. Thirdly, to solve the input saturation, an auxiliary dynamic system is introduced. Fourthly, the barrier Lyapunov function is used to achieve connectivity preservation, collision avoidance and swarm control. Meanwhile, by using the estimated velocities of vessels, an output feedback controller is designed. The stability of the closed-loop system is proved. The simulation results show the effectiveness of the proposed swarm control strategy.

1. Introduction

In recent years, the swarm control of multiple unmanned surface vehicles (USVs) has attracted extensive attention. The purpose of swarm control is to drive multiple agents to accomplish tasks uniformly, which is very different from the traditional control based on a single object [1]. Swarm control can be applied to a variety of tasks, such as arranging a group of USVs to form a safe area to protect other vessels [2,3], tracking floating pollutants [4], monitoring environment [5,6], etc. However, USVs will be affected by multiple constraints in operation. These constraints can be summarized as model parameter uncertainty, unavailability of velocity measurements, time-varying environmental disturbances, input saturation, and output constraints. These problems will reduce the performance of the system, and even lead to instability. Therefore, it is necessary to analyze these constraints and to design a swarm control strategy.
In applications, in addition to position information, velocity information is also very important. Generally, position information can be easily obtained by GNSS equipment, while velocity information may not be accurately computed by the GNSS receiver. In the high-speed operation of USVs, it is more difficult to obtain velocity information that is sufficiently timely for use in a control algorithm. Meanwhile, velocity obtained from position measurements by numerical differentiation techniques are often unfeasible due to inevitable measurement noises [7]. Therefore, it is significant to study output feedback schemes independent of the velocity measurement [8].
Some control strategies without velocity information have been proposed in the literature [9,10,11,12,13] such as the back-stepping observer approach [9], high-gain observer [10,11,12], passive nonlinear observer [13], etc. These control schemes have a common premise, that is, it is assumed that the vessel model parameter is known. However, due to imprecise measurements and external disturbances, accurate model parameters often cannot be obtained. In general, the hydrodynamic parameters of USVs are also time varying. It is emphasized that external disturbances and uncertainties can make the entire system unstable [14]. This increases the difficulty of measurement. In order to deal with uncertain dynamics, neural network (NN) approaches have been considered due to their inherent advantages including excellent approximation and learning performance [1,7,15,16,17,18,19,20,21,22].
Due to the power and response speed constraints of the vehicle’s actuator, and the output constraints of the propeller, the vessel’s control forces and moment cannot be infinite. Therefore, the amplitude of the control signal is usually limited to a certain range. Moreover, due to the requirements of different operational scenarios (e.g., acceleration from stationary to maximum speed, high sea state operation, etc.), the ship will have input saturation constraints. Otherwise, control input overshoot may occur, which will affect the performance of the control system and even lead to system instability. In [23], an adaptive steering control for uncertain ship dynamics with input constraints was designed. In [24], a dynamic surface control scheme was proposed for a class of uncertain strict-feedback nonlinear systems subject to input saturation. Furthermore, auxiliary dynamic systems have been used extensively for nonlinear systems design due to their efficiency and design flexibility [16,25,26,27].
In addition to input constraints, another challenge to swarm control for USVs is output (position) constraints. Prevailing collision avoidance approaches can be divided into two methodologies, potential functions [28,29,30,31,32] and prescribed performance functions [27,33,34,35,36]. Potential function approaches introduce a potential energy function that increases with the decrease in the distance between USVs. Then, a control low based on the gradient potential energy function is designed to minimize the potential energy of the system [28,29,30,31,32]. However, this method requires the design of another potential energy function in the Lyapunov function. Thus, there may be a conflict when choosing parameters that must simultaneously ensure the stability criteria are met and provide the required collision avoidance behavior. For the second method, the barrier Lyapunov function (BLF) is an effective way to prevent the violation of constraints for many practical systems. Many successful applications of this approach exist in the literature [27,33,34,35,36].
Motivated by the aforementioned observations, this work investigates the swarm control of USVs considering model uncertainties, the unavailability of velocity measurements, time-varying environmental disturbances, input saturation and output constraints. In summary, the main contributions are as follows.
  • A neural adaptive state observer is designed to recover velocity information and to estimate composite disturbances including model uncertainty and time-varying environmental disturbances.
  • An auxiliary dynamic system is introduced to deal with input saturation. A modified BLF is provided to achieve connectivity preservation, collision avoidance and swarm control.
  • In combination with the observer, an output feedback controller is proposed for the follower USVs based on a second-order linear tracking differentiator, an adaptive law, a modified BLF and graph theory. Meanwhile, the stability of the closed-loop system is proved via Lyapunov theory.
This paper is organized as follows. A table of notations and some variables used in the paper is presented (Table 1). Section 2 describes some preliminaries and mathematical modeling. Section 3 provides the neural adaptive state observer design. Section 4 presents the output feedback controller design and analyzes the stability by the Lyapunov method. Section 5 compares the results of simulations to verify the effectiveness of the proposed scheme. Section 6 concludes this paper.

2. Preliminaries and Mathematical Modeling

2.1. Algebraic Graph Theory

Graph theory is used to describe the communication topology of n follower USVs and a virtual leader vehicle (denoted by 0). A directed graph G = ( V , E ) consists of a vertex set V = { 0 , 1 , 2 , , n } and the set of edges E { ( i , j ) V × V } . If ( i , j ) E , node j is an adjacent node of node i. N i = { j V , ( i , j ) E } represents the set of all adjacent nodes of node i, as can be seen in [37].
Consider a directed graph G composed of n nodes; the adjacency matrix A = [ a i j ] n × n represents the link relationship between nodes, where a i j = 1 , if ( i , j ) E ; a i j = 0 , otherwise. If a i j = a j i , the graph is undirected; otherwise, it is directed. The Laplacian matrix is defined as L = D A , where D = diag { d 1 , d 2 , , d n } with d i = j = 1 n a i j .
In particular, a diagonal matrix A 0 = diag { a i 0 } is defined as a leader adjacency matrix, where a i 0 = 1 , if and only if the ith USV receives information from the virtual leader vehicle; a i 0 = 0 , otherwise. Finally, the information exchange matrix is defined as B = L + A 0 .
Assumption A1.
The graph G is directed, and there is at least one spanning tree from root node to the leader node, i.e., B is a positive definite matrix.

2.2. Barrier Lyapunov Function

Consider a continuous system
x ˙ = f ( x ) , x D ,
where D is an open region containing the origin. If a continuously differentiable, positive definite function V ( x ) satisfies lim x D V = + , and V ( x ( t ) ) b , t 0 , where x ( 0 ) D , b > 0 is a constant. Then, V ( x ) is a barrier Lyapunov function, see [33,38].
A barrier Lyapunov function candidate is as follows
V = 1 2 ln k 2 k 2 z 2 ,
where k > 0 is a constant, k < z < k .
The following lemma formalizes a result on the use of a BLF candidate for constraint satisfaction.
Lemma 1.
For any positive constant k, and any z R satisfying k < z < k , we have
ln k 2 k 2 z 2 < z 2 k 2 z 2 .
Proof. 
Let p 1 = z 2 k 2 z 2 ln k 2 k 2 z 2 , we have
p 1 = z 2 k 2 z 2 ln ( 1 z 2 k 2 z 2 ) .
As k < z < k , one has z 2 < k 2 . Then, we have the inequalities as follows:
0 z 2 k 2 z 2 < k 2 k 2 z 2 = 1 z 2 k 2 z 2 ,
0 z 2 k 2 z 2 < 1 2 ,
Let p 2 = z 2 k 2 z 2 , p 3 = ln ( 1 p 2 ) + p 2 . The derivative of p 3 is
p 3 ˙ = 2 p 2 1 p 2 > 0 .
This shows that p 3 is continuously increasing and the minimum of p 3 is p 3 m i n = 0 . Thus, we have ln k 2 k 2 z 2 < z 2 k 2 z 2 . □

2.3. Neural Network

For any real continuous function f ( ς ) : R n R k on a sufficiently large compact set Ω , there exists a radial basis function neural network (NN) such that
f ( ς ) = W T h ( ς ) + ε ,
where input vector ς Ω R n . W = [ ω 1 , ω 2 , , ω l ] T R k × l represents the output weight vector, satisfying W W * with W * > 0 being a positive constant. The NN node number l > 1 . h ( ς ) = [ h 1 ( ς ) , , h l ( ς ) ] T . Since the neural function is bounded, there exists a positive constant h max > 0 such that h ( ς ) h max . ε R k is the approximation error, satisfying ε ε * with ε * > 0 being a positive constant, see [22].
Define W ^ as an estimate of W, and let the estimation error be denoted by W ˜ = W ^ W . The optimal weight vector W is rewritten as
W = arg min W ^ R m × n sup ς Ω f ( ς ) W ^ T h ( ς ) ,
In this paper, we choose h i ( ς ) as a Gaussian function, as follows
h i ( ς ) = exp ( ς c i ) T ( ς c i ) d i 2 ,
where c i = [ c i 1 , c i 2 , , c i n ] T is the center of the receptive field, and d i is the width of the Gaussian function.

2.4. USVs Modeling

Consider a group of USVs consisting of a virtual leader vehicle (subscript 0) and n follower USVs (subscripts 1 , 2 , , n ). The 3-degree-of-freedom (DOFs) kinematics and dynamics equations of the ith USV can be expressed in vector form as [39]
η ˙ i = R i ψ i υ i ,
M i υ ˙ i + D i υ i = τ i + d i ,
where R i ( ψ i ) is the rotation matrix, given as
R i ( ψ i ) = cos ( ψ i ) sin ( ψ i ) 0 sin ( ψ i ) cos ( ψ i ) 0 0 0 1 ,
with properties: | R i ( ψ i ) | = 1 and R i T ( ψ i ) R i ( ψ i ) = I 3 . η i : = [ x i , y i , ψ i ] T is the position and yaw angle in the earth-fixed frame X E O E Y E (see Figure 1). υ i : = [ u i , v i , r i ] T is the velocity vector in the body-fixed frame X B O B Y B . The system inertia matrix M i R 3 × 3 is positive definite and constant. The damping matrix D i R 3 × 3 is also defined as symmetric and positive. The control input τ i : = [ τ i 1 , τ i 2 , τ i 3 ] T is the control input, which is produced by the a propeller and a rudder, etc. d i : = [ d i 1 , d i 2 , d i 3 ] T is a time-varying environmental disturbance.
In this paper, the input saturation can be described as
τ i = τ i , max , if τ i c > τ i , max τ i c , if τ i , min τ i c τ i , max , τ i , min , if τ i c < τ i , min
where τ i , max R 3 and τ i , min R 3 are the maximum and minimum control forces and moment of the ith vehicle, respectively. τ i c = [ τ i c 1 , τ i c 2 , τ i c 3 ] T is calculated by the controller.
The main goal of this paper is to design an output feedback controller for each USV to track the reference signal ( η d ) subject to model uncertainties and input constraints only using position measurements. Specifically, it is to achieve the following objectives.
(1) Formation objective: the ith USV is driven to a formation pattern with relative position and heading, and we have:
lim t + η i η d μ i σ i ,
where η i : = [ x i , y i , ψ i ] T . η d : = [ x d , y d , ψ d ] T is the desired reference point. μ i : = [ x i μ , y i μ , ψ i μ ] T represents the expected offset of the vehicle’s position and heading relative to the desired reference point, which is the position of the virtual leader. To maintain the USVs in a fixed formation, μ i is a constant vector. σ i > 0 is a constant.
(2) Collision avoidance: the position for ith USV must remain within the set
Ω η i = η i μ i < k i c ,
where k i c = k i b λ min ( H ) + η d + μ i , see Section 4.3 for details.
Assumption A2.
The reference signal η d is smooth and differentiable everywhere. Its first derivative η ˙ d and second derivative η ¨ d exist and are bounded.
Assumption A3.
The coincident disturbances, which mainly include model uncertainties and time-varying environmental disturbances, are bounded.
Assumption A4.
The position of each USV is available.
Assumption A5.
The initial positions of USVs are meeting the maximum collision avoidance distance, i.e., z i 1 ( 0 ) Ω z i 1 , i = 1 , , n . USVs are always within the communication range.

2.5. Environmental Disturbances Modeling

Unmodeled external forces and moments due to wind, ocean currents, and second waves are lumped together into an earth-fixed slowly varying bias term d i R 3 [13,40]. A widely used bias model for USVs is the first-order Markov process. In this paper, the environmental disturbances are modeled as
d i = R i T b ,
where b represents a first-order Markov process, given as
b ˙ = T 1 b + E b ϑ b ,
where T R 3 × 3 is a diagonal matrix of positive bias time constants, E b R 3 × 3 is a diagonal matrix scaling the amplitude of ϑ b , and ϑ b R 3 is a vector of zero-mean Gaussian white noise.
Remark 1.
The USV in this paper has a low speed and the part above the waterline (i.e., superstructure) is small, and only calm sea conditions are considered in this paper. Therefore, the impact of sea wind on the USV is ignored. Because the draft of the USV is very shallow, the ocean current has little impact on it. Therefore, the impact of the ocean current on the USV is not considered. The wave disturbance consists of a low-frequency part and high-frequency part. The oscillation motion caused by the high-frequency part shall not enter the feedback control circuit. In this paper, only the low-frequency part is considered during the control process.

3. Neural Adaptive State Observer Design

In engineering applications, the parameters M i and D i cannot always be measured accurately. Under such circumstances, according to Equations (11) and (12), an approximate model for the dynamics of the ith USV can be written as
η ˙ i = R i υ i ,
M ¯ i υ ˙ i = τ i f i ( υ i ) ,
where M ¯ i T = M ¯ i is the nominal inertial matrix, which is positive definite. A new variable ς i is used as the input vector of the NN. ς i = [ η i T , η i T ( t t d ) , η i T ( t 2 t d ) , τ i T ] T , and t d is a positive constant. In Equation (8), the f i ( υ i ) = [ f i 1 , f i 2 , f i 3 ] T is written as
f i ( υ i ) = W i T h i j ( ς i ) + ε i ,
where W i = [ W i 1 , W i 2 , W i 3 ] T R m × 3 . j = 1 , 2 , , m represents the jth neuron. ε i R 3 is the approximation error, satisfying ε i ε i , max with ε i , max R being a positive constant. h i j ( ς i ) R m × 1 , which satisfies h i j ( ς i ) h i j , max with h i j , max R being a positive constant.
Let η ^ i = [ x ^ i , y ^ i , ψ ^ i ] T represent an estimation of η i , and υ ^ i = [ u ^ i , v ^ i , r ^ i ] T represent an estimate of υ i . Define the position estimation error as η ˜ i = η ^ i η i . We design a neural adaptive state observer (NASO) as
η ^ ˙ i = R i υ ^ i K o i 1 η ˜ i ,
M ¯ i υ ^ ˙ i = τ i W ^ i T h i j ( ς i ) K o i 2 R i T η ˜ i ,
where K o i 1 R 3 × 3 and K o i 2 R 3 × 3 are positive definite diagonal gain matrices. W ^ i = [ W ^ i 1 , W ^ i 2 , W ^ i 3 ] T is an estimate of W i . We design a weight adaptive update law for W ^ i as
W ^ ˙ i k = γ i k h i j ( ς i ) p i k k i w W ^ i k , k = 1 , 2 , 3
where p i k = η ˜ i T R i = [ p i 1 , p i 2 , p i 3 ] , and γ i k > 0 and k i w > 0 are constants.
Define the velocity estimation error as υ ˜ i = υ ^ i υ i . Based on the above analysis, the dynamic error estimation equations of the NASO can be written as
η ˜ ˙ i = R i υ ˜ i K o i 1 η ˜ i ,
M ¯ i υ ˜ ˙ i = W ˜ i T h i j ( ς i ) K o i 2 R i T η ˜ i + ε i ,
To facilitate the stability analysis of the NASO, a new variable X i = [ η ˜ i T , υ ˜ i T ] T is defined. Equations (25) and (26) then become
X ˙ i = A i X i + B i W ˜ i T h i j ( ς i ) + ε i ,
η ˜ i = C i X i ,
where A i , B i and C i are defined as:
A i = K o i 1 R i K o i 2 M ¯ i 1 R i T 0 3 , B i = 0 3 M ¯ i 1 , C i = I 3 0 3 .
For convenience, a variable χ i = T i X i is introduced with T i = diag { R i T , I 3 } . Then, (27) can be written as
χ ˙ i = A i 0 + r i S T χ i + B i W ˜ i T h i j ( ς i ) + ε i ,
where r i is the yaw rate, S T = diag { S T , 0 3 } , and S and A i 0 are defined as:
S = 0 1 0 1 0 0 0 0 0 , A i 0 = K o i 1 I 3 K o i 2 M ¯ i 1 0 3 .
Similarly, Equation (24) is written as
W ^ ˙ i k = γ i k h i j ( ς i ) s i k k i w W ^ i k , k = 1 , 2 , 3
where s i k = χ i T C i T = [ s i 1 , s i 2 , s i 3 ] .
Using this notation, we can show that the following theorem holds:
Theorem 1.
The NASO estimation error is bounded; if the NASO is defined using Equations (27) and (28), the weight adaptive update law is defined using Equation (32), the parameters satisfy k i w γ i k h i j , max 2 > 0 and there exist positive definite symmetric matrices Q i , P i R 6 × 6 such that the linear matrix inequalities (LMIs) are satisfied:
A i 0 T P i + P i A i 0 + P i B i B i T P i + Q i + F i F i T + r i , max ( S T T P i + P i S T ) 0 ,
A i 0 T P i + P i A i 0 + P i B i B i T P i + Q i + F i F i T r i , max ( S T T P i + P i S T ) 0 ,
where F i = C i T P i B i ; r i , max is the upper bound of r i , satisfying r i r i , max with r i , max R being a positive constant.
Proof. 
We choose the Lyapunov function V i o as
V i o = 1 2 χ i T P i χ i + 1 2 k = 1 3 1 γ i k W ˜ i k T W ˜ i k ,
The time derivative of V i o is
V ˙ i o = 1 2 χ i T P i A i 0 + A i 0 T P i + r i P i S T + r i S T T P i χ i + χ i T P i B i W ˜ i T h i j ( ς i ) + ε i + k = 1 3 1 γ i k W ˜ i k T W ˜ ˙ i k
Let λ i = 1 2 χ i T P i A i 0 + A i 0 T P i + r i P i S T + r i S T T P i χ i + χ i T P i B i W ˜ i T h i j ( ς i ) + ε i , then
V ˙ i o = λ i + k = 1 3 1 γ i k W ˜ i k T W ˜ ˙ i k = λ i + k = 1 3 1 γ i k W ˜ i k T ( W ^ ˙ i k W ˙ i k ) = λ i + k = 1 3 W ˜ i k T h i j ( ς i ) s i k k = 1 3 k i w γ i k W ˜ i k T W ^ i k = λ i + W ˜ i T h i j ( ς i ) χ i T C i T k = 1 3 k i w γ i k W ˜ i k T W ^ i k
Substituting Equations (33) and (34) into (37), then
V ˙ i o 1 2 χ i T P i A i 0 + A i 0 T P i + r i P i S T + r i S T T P i χ i + χ i T F i W ˜ i T h i j ( ς i ) + χ i T P i B i ε i k = 1 3 k i w γ i k W ˜ i k T W ^ i k
Using Young’s inequality [41], we have
χ i T F i W ˜ i T h i j ( ς i ) 1 2 χ i T F i F i T χ i + h i j , max 2 2 W ˜ i T W ˜ i χ i T P i B i ε i 1 2 χ i T P i B i B i T P i χ i + 1 2 ε i , max 2 W ˜ i k T W ^ i k = W ˜ i k T ( W ˜ i k + W i k ) W ˜ i k 2 + 1 2 W ˜ i k 2 + 1 2 W i k 2 1 2 W ˜ i k 2 + 1 2 W i k 2
Substituting Equation (39) into (38), then
V ˙ i o 1 2 χ i T P i A i 0 + A i 0 T P i + r i P i S T + r i S T T P i + F i F i T + P i B i B i T P i χ i + h i j , max 2 2 W ˜ i T W ˜ i + 1 2 ε i , max 2 + k = 1 3 k i w 2 γ i k W ˜ i k 2 + W i k 2
Substituting Equations (33) and (34) into (40), then
V ˙ i o 1 2 λ min ( Q i ) χ i 2 k = 1 3 1 2 γ i k k i w γ i k h i j , max 2 W ˜ i k 2 + k = 1 3 k i w 2 γ i k W i k 2 + 1 2 ε i , max 2 a i 1 V i o + a i 2
As a result of the above control low, V ˙ i o a i 1 V i o + a i 2 and this results in the ultimately uniformly bounded regulation of the state.
Then, Equation (41) can be written as
V ˙ i o ( t ) a i 1 V i o ( t ) + a i 2 ,
where a i 1 = min λ min ( Q i ) λ max ( P i ) , k i w γ i k h i j , max 2 > 0 , and a i 2 = k = 1 3 k i w 2 γ i k W i k 2 + 1 2 ε i , max 2 , 0 < a i 2 a i 2 , max with a i 2 , max R being a positive constant.
The coefficients in the observer design process are determined by Equation (42). Since the observer design process needs to meet the Lyapunov theory, the selection of parameters needs to meet Equation (42).
Then, Equation (42) becomes
V i o ( t ) V i o ( 0 ) a i 2 a i 1 e a i 1 t + a i 2 a i 1 ,
From Equation (43), when t , V i o a i 2 a i 1 . Thus, the signals in closed-loop are bounded. Then, the states χ i are bounded. Noticing T i T = T i 1 , T i T 1 and using X i = T i T χ i , the estimation error signal X i is bounded.
This completes the proof. □

4. Output Feedback Controller Design

Firstly, a neural network adaptive state observer is designed, which can reconstruct the velocities and estimate the coincidence disturbances. Secondly, an ADS is designed to deal with the problem of input saturation. Thirdly, in order to solve the problem that the first-order filter is sensitive to noise in the traditional dynamic surface control technology, a second-order linear tracking differentiator (SOLTD) is introduced. Then, an output feedback controller is designed, as shown in Figure 2.

4.1. Auxiliary Dynamic System

In this subsection, an auxiliary dynamic system (ADS) is introduced to solve the input saturation problem, as can be seen in Figure 2. For ADS, the system input is the deviation ( Δ τ i ), the output is the velocity tracking error compensation ( β i 1 ) for ith USV. The ADS is designed as
β ˙ i 1 = L i 1 β i 1 + M ¯ i 1 Δ τ i
where L i 1 R 3 × 3 is a positive definite diagonal matrix. Δ τ i = τ i τ i c , Δ τ i Δ τ i , max , Δ τ i , max > 0 is a constant.

4.2. Output Feedback Controller Design

In this subsection, an output feedback controller for multiple USVs is designed using dynamic surface control technology. The design process is divided into the following steps.
Step 1: According to the communication topology between USVs, the first tracking error of the ith USV in the earth-fixed frame is defined as
z i 1 = j N i a i j η i η j μ i j + a i 0 ( η i η d μ i ) ,
where N i , a i j and a i 0 are defined in Section 2.1. η i , η d and μ i are explained in Equation (15), η j has similar definition, μ i j = μ i μ j .
For the time derivative of z i 1 , we obtain
z ˙ i 1 = a i d R i υ i j N i a i j R j υ j a i 0 η ˙ d ,
where a i d = d i + a i 0 , d i and a i 0 are defined in Section 2.1.
The kinematic control low ( α i ) of the ith USV is designed as
α i = R i T a i d K i 1 z i 1 κ i a z i 3 + j N i a i j R j υ ^ j + a i 0 η ˙ d
where K i 1 R 3 × 3 is a positive definite diagonal matrix. κ i a = 2 a i d + d i 2 > 0 is a constant. z i 3 = [ z i 3 , 1 , z i 3 , 2 , z i 3 , 3 ] T with z i 3 , l = z i 1 , l k i b , l 2 z i 1 , l 2 , l = 1 , 2 , 3 , z i 1 = [ z i 1 , 1 , z i 1 , 2 , z i 1 , 3 ] T . Define a compact set Ω z i 1 = z i 1 k i b < z i 1 < k i b , where k i b = [ k i b , 1 , k i b , 2 , k i b , 3 ] T . Similarly, R j is the rotation matrix of jth USV. υ ^ j is the velocity vector estimation value.
To avoid the calculation of the time derivative of α i , and considering the noise sensitivity of the first-order low-pass filter, a SOLTD is introduced as follows
υ ˙ i r = υ i r d υ ˙ i r d = ι i 2 ( υ i r α i ) 2 ι i υ i r d
where ι i > 0 is a time constant, υ i r R 3 is the output vector of SOLTD.
Step 2: The second tracking error of the ith USV in earth-fixed frame is defined as
z i 2 = υ ^ i υ i r β i 1 .
Using the Equations (23) and (44), the time derivative of z i 2 , we have
M ¯ i z ˙ i 2 = K o i 2 R i T η ˜ i W ^ i h i j ( ς i ) + τ i c M ¯ i ( υ i r d L i 1 β i 1 )
The dynamic equation of the ith USV is designed as
τ i c = K i 2 z i 2 + W ^ i h i j ( ς i ) + M ¯ i ( υ i r d L i 1 β i 1 ) .
where K i 2 R 3 × 3 is a positive definite diagonal matrix.
Substituting Equations (47) and (51) into (46) and (50), we obtain
z ˙ i 1 = K i 1 z i 1 k i a z i 3 + a i d R i ( υ ˜ i + β i 1 + z ^ i 2 + q i ) + d i R j υ ˜ j ,
M ¯ i z ˙ i 2 = K i 2 z i 2 K o i 2 R i T η ˜ i ,
where q i = υ i r α i , q i q i , max , q i , max > 0 is a constant.

4.3. Stability Analysis

Theorem 2.
Consider a closed-loop system subject to model parameter uncertainties, velocity measurements not being available, time-varying environmental disturbances, input saturation and output constraints. Let the USV dynamics be given by Equations (19) and (20), NASO by Equations (22) and (23), weight adaptive update law by Equation (24), ADS by (44), and kinematic and dynamic equations by Equations (47) and (51). Suppose that the closed-loop system satisfies Assumptions 1–5. Then, the following statements hold.
(i)
All signals in the closed-loop system are uniformly ultimately bounded.
(ii)
All USVs track the reference signal with a bounded tracking error.
(iii)
The output position of each USV satisfies output constraints.
Proof. 
Consider the above closed-loop system, and the Lyapunov function V i is chosen as
V i = 1 2 i = 1 n 2 V i o + l = 1 3 ln k i b , l 2 k i b , l 2 z i 1 , l 2 + z i 2 T M ¯ i z i 2 + β i 1 T β i 1
For the time derivative of V i , we obtain
V ˙ i = i = 1 n { V ˙ i o + z i 3 T z ˙ i 1 + z i 2 T M ¯ i z ˙ i 2 + β i 1 T β ˙ i 1 } = i = 1 n { V ˙ i o z i 3 T K i 1 z i 1 k i a z i 3 T z i 3 + z i 3 T [ a i d R i ( υ ˜ i + β i 1 + z ^ i 2 + q i ) + d i R j υ ˜ j ] z i 2 T K i 2 z i 2 z i 2 T K o i 2 R i T η ˜ i β i 1 T L i 1 β i 1 + β i 1 T M ¯ i 1 Δ τ i }
Using Young’s inequality [41], we have
z i 3 T a i d R i υ ˜ i a i d 2 z i 3 T z i 3 + a i d 2 υ ˜ i 2 z i 3 T a i d R i β i 1 a i d 2 z i 3 T z i 3 + a i d 2 β i 1 2 z i 3 T a i d R i z i 2 a i d 2 z i 3 T z i 3 + a i d 2 z i 2 2 z i 3 T a i d R i q i a i d 2 z i 3 T z i 3 + a i d 2 q i , max 2 z i 3 T d i R j υ ˜ j d i 2 z i 3 T z i 3 + d i 2 υ ˜ j 2 z i 2 T K o i 2 R i T η ˜ i λ max ( K o i 2 ) 2 z i 2 2 + η ˜ i 2 β i 1 T M ¯ i 1 Δ τ i λ max ( M ¯ i 1 ) 2 β i 1 2 + Δ τ i , max 2
Substituting Equation (56) into (55), we obtain:
V ˙ i i = 1 n { V ˙ i o λ min ( K i 1 ) l = 1 3 ln k i b , l 2 k i b , l 2 z i 1 , l 2 λ min ( K i 2 ) a i d + λ max ( K o i 2 ) 2 z i 2 2 λ min ( L i 1 ) a i d + λ max ( M ¯ i 1 ) 2 β i 1 2 + λ max ( K o i 2 ) 2 η ˜ i 2 + a i d + d i 2 υ ˜ i 2 + a i d 2 q i , max 2 + λ max ( M ¯ i 1 ) 2 Δ τ i , max 2 } i = 1 n { V ˙ i o b i 1 l = 1 3 ln k i b , l 2 k i b , l 2 z i 1 , l 2 b i 2 z i 2 2 b i 3 β i 1 2 + b i 4 η ˜ i 2 + b i 5 υ ˜ i 2 + b i 6 } i = 1 n { V ˙ i o b i 1 l = 1 3 ln k i b , l 2 k i b , l 2 z i 1 , l 2 b i 2 z i 2 2 b i 3 β i 1 2 + λ max ( A i 1 ) χ i 2 + b i 6 } .
where we select appropriate parameters K i 1 , K i 2 , K o i 2 and L i 1 to meet:
b i 1 = λ min ( K i 1 ) > 0 ,
b i 2 = λ min ( K i 2 ) a i d + λ max ( K o i 2 ) 2 > 0 ,
b i 3 = λ min ( L i 1 ) a i d + λ max ( M ¯ i 1 ) 2 > 0 ,
b i 4 = λ max ( K o i 2 ) 2 > 0 ,
b i 5 = a i d + d i 2 > 0 ,
and b i 6 = a i d 2 q i , max 2 + λ max ( M ¯ i 1 ) 2 Δ τ i , max 2 , A i 1 = diag { b i 4 , b i 5 } .
Substituting Equation (40) into (57), we obtain
V ˙ i i = 1 n { a i 1 V i o b i 1 l = 1 3 ln k i b , l 2 k i b , l 2 z i 1 , l 2 b i 2 z i 2 2 b i 3 β i 1 2 + λ max ( A i 1 ) χ i 2 + a i 2 + b i 6 } i = 1 n { a i 1 λ max ( A i 1 ) λ min ( P i ) V i o b i 1 l = 1 3 ln k i b , l 2 k i b , l 2 z i 1 , l 2 b i 2 z i 2 2 b i 3 β i 1 2 + a i 2 + b i 6 } c i 1 V i + c i 2
As a result of the above control low, V ˙ i c i 1 V i + c i 2 and this results in the ultimately uniformly bounded regulation of the state.
Then, Equation (63) can be written as
V ˙ i o ( t ) c i 1 V i o ( t ) + c i 2 ,
where c i 1 = min a i 1 λ max ( A i 1 ) λ min ( P i ) , 2 b i 1 , 2 b i 2 λ max ( M ¯ i ) , 2 b i 3 > 0 , and c i 2 = i = 1 n a i 2 + b i 6 ,
0 < c i 2 c i 2 , max with c i 2 , max R being a positive constant.
The coefficients in the controller design process are determined by Equation (64). Since the controller design process needs to meet the Lyapunov theory, the selection of parameters needs to meet Equation (64).
From the Equation (64), when t , V i o c i 2 c i 1 . Then, Equation (64) becomes
V i o ( t ) V i o ( 0 ) c i 2 c i 1 e c i 1 t + c i 2 c i 1 ,
According to the definition of the Lyapunov function V i and Theorem 1, it can be concluded that χ i , W ˜ i , z i 2 , β i 1 and ln k i b , l 2 k i b , l 2 z i 1 , l 2 are bounded, where ln k i b , l 2 k i b , l 2 z i 1 , l 2 implies that z i 1 always remains within the set Ω z i 1 . Therefore, all signals in the closed-loop system are uniformly ultimately bounded.
Therefore, conclusion (i) is valid.
We now prove (ii) of Theorem 2. The tracking error of the ith USV in the earth-fixed frame is defined as δ i = η i η d μ i . According to the definition of z i 1 in Equation (45), we have
z 1 = H I 3 δ ,
where z 1 = [ z 11 T , z 21 T , , z n 1 T ] T and δ = [ δ 1 T , δ 2 T , , δ n T ] T . H is defined in Section 2.1.
According to Assumption 1, all eigenvalues of matrix H have positive real parts. Thus, we obtain
δ z 1 λ min ( H )
According to Equation (63), the variable z i 1 satisfies z i 1 c i 2 , max / b i 1 . Therefore, δ is bounded. Thus, all USVs track the reference signal with a bounded tracking error. These prove that ii) holds.
From the Lemma 1 and BLF candidate ln k i b , l 2 k i b , l 2 z i 1 , l 2 , we obtain z i 1 , l < k i b , l , i.e., k i b , l < z i 1 , l < k i b , l .
Using Equation (66) and δ i = η i η d μ i , we have
η i k i b λ min ( H ) + η d + μ i
Let k i c = k i b λ min ( H ) + η d + μ i . That is, the output of each USV is always kept in set Ω η i = η i k i c < η i < k i c . Therefore, conclusion iii) is valid.
This completes the proof. □

5. Simulation Results

Simulations are carried out with Matlab 2018a. The simulations are run on a PC with a dual-core 2.30 GHz Intel(R) Core(TM) i5-8300H CPU and 8 GB of RAM.
In this section, we simulate a USV swarm consisting of one virtual leader vehicle (indexed by 0) and four follower USVs (indexed by 1 , 2 , 3 , 4 , 5 , 6 ) to demonstrate the effectiveness of the proposed method. The directed communication graph is shown in Figure 3.
In simulations, the model of surface ship Cybership II is used [42]. The time-varying environmental disturbances are modeled as a first-order Markov process [40]. The control forces and moment are limited as τ i 1 , max = τ i 1 , min = 2 N , τ i 2 , max = τ i 2 , min = 2 N and τ i 3 , max = τ i 3 , min = 1.5 Nm . The constraint k i b is set as k i b = [ 0.6 m , 0.6 m , 0.2 r a d ] T . The desired reference trajectory η d is generated as Equation (69). Some parameters are set as shown in Table 2.
η d = [ t 10 , 0 , 0 ] , 0 t < 50 ; [ 5 sin ( 1 50 ( t 50 ) ) + 5 , 5 cos ( 1 50 ( t 50 ) ) + 5 , 1 50 ( t 50 ) ] , 50 t < 50 + 50 π ; [ 5 sin ( 1 50 ( t 50 50 π ) ) + 5 , 5 cos ( 1 50 ( t 50 50 π ) ) + 15 , 1 50 ( t 50 50 π ) ] , 50 + 50 π t < 50 + 100 π ; [ 1 10 ( t 50 100 π ) , 20 , 0 ] , 50 + 100 π t 400 ;

5.1. Performance of Proposed Control Strategy

In this subsection, the simulation results are given to verify the performance of the proposed control strategy. The parameters of environmental disturbance are selected as T = d i a g { 10 3 , 10 3 , 10 3 } , b ( 0 ) = [ 0 , 0 , 0 ] T , E b = d i a g { 10 5 , 10 5 , 10 5 } , ϑ b = [ 0.1 , 0.1 , 0.1 ] T . The parameters of the observer are selected as K o i 1 = d i a g { 30 , 30 , 30 } , K o i 2 = d i a g { 30 , 30 , 30 } , γ i = 1000 and k i = 0.2 . The parameters of the controller are chosen as ι i = 0.001 , K i 1 = d i a g { 0.6 , 0.6 , 0.6 } and K i 2 = d i a g { 1.1 , 1.1 , 1.1 } .
The simulation results are shown in Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8. Figure 4 shows the trajectories of six USVs under time-varying environmental disturbance and input saturation. It can be shown that, after adjustment, the position and heading of six USVs can reach the present point, maintain the desired relative position between each other, and realize the swarm control of multiple USVs in a fixed formation, even if there is a position and heading deviation at the initial time. Figure 5 describes the tracking errors of four USVs, including the position and heading. It can be seen from Figure 5 that the tracking error is convergent. Figure 6 shows the control input of six USVs. According to the analysis of Figure 6, the control input is limited within the range of input constraints. Therefore, the input saturation can be realized using the ADS designed in this paper.
Taking USV 1 as an example, the estimation effects of the observer are shown in Figure 7 and Figure 8. Figure 7 depicts the comparison effect between the real value and the estimated value of the velocity. The solid lines and dashed line represent the real value and estimated value, respectively. It can be seen from Figure 7 that the designed NASO is effective and can realize the estimation of speed. Figure 8 describes the real values and estimated values of the total unknown function acting on the first USV. It can be seen from Figure 8 that the NASO is stable despite boundary estimate errors.

5.2. Comparison Group

To further illustrate the effectiveness of the proposed control strategy, a neural adaptive dynamic surface control (NADSC) approach without ADS is considered, whose motion is described by
z i 1 = j N i a i j η i η j μ i j + a i 0 ( η i η d μ i ) α i = R i T a i d K i 1 z i 1 κ i a z i 3 + j N i a i j R j υ ^ j + a i 0 η ˙ d υ ˙ i r = υ i r d = ι i 2 ( υ i r α i ) z i 2 = υ ^ i υ i r τ i c = K i 2 z i 2 + W ^ i h i j ( ς i ) + M ¯ i υ i r d
where the parameters are set to be the same as the proposed controller.
The simulation results are shown in Figure 9, Figure 10 and Figure 11. Figure 9 describes the tracking error of six USVs without ADS. Compared with Figure 5, it can be concluded that the tracking error fluctuates greatly at 0–10 s. Figure 10 shows the control input of six USVs without ADS. Compared with Figure 6, it can be concluded that the control input fluctuates greatly at 0–10 s. Figure 11 depicts the control inputs of six USVs without input constraints. It can be seen from Figure 11 that the maximum control force and moment is τ i = [ τ i 1 , τ i 2 , τ i 3 ] T = [ 200 N , 200 N , 150 Nm ] T . Compared with Figure 6, the maximum control force and torque of the USV were exceeded.
To summarize, the proposed approach combining the NASO, ADS and LBF successfully handled the output feedback swarm control problem with satisfactory results.

6. Conclusions

This paper investigated swarm control for USVs in the presence of model uncertainty, the unavailability of velocity measurements, unknown environmental disturbances, input saturation and output constraints. NASO is designed to estimate the unknown model uncertainty, unmeasured velocity and unknown environmental disturbance. ADS is introduced to mitigate the input saturation problem. An output feedback controller is designed, which is composed of NASO, ADS, SOLTD and BLF. The stability of the system is proved via the Lyapunov method. Finally, the effectiveness of the proposed control strategy is verified in simulation. However, in the process of designing the output feedback controller, we do not take the problem of obstacle avoidance into account. Then, this will be the focus of our research.

Author Contributions

Conceptualization, G.X. and X.S.; methodology, X.S.; software, X.S. and X.X.; formal analysis, X.S. and X.X.; investigation, X.S.; supervision, G.X. and X.X.; funding acquisition, G.X. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the 7th Generation Ultra Deep Water Drilling Unit Innovation Project.

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:
USVunmanned surface vehicle
BLFbarrier Lyapunov function
NNneural network
NASOneural adaptive state observer
LMIslinear matrix inequalities
SOLTDsecond-order linear tracking differentiator
ADSauxiliary dynamic system
NADSCneural adaptive dynamic surface control

References

  1. Peng, Z.; Wang, J.; Wang, D. Distributed containment maneuvering of multiple marine vessels via neurodynamics-based output feedback. IEEE Trans. Ind. Electron. 2017, 64, 3831–3839. [Google Scholar] [CrossRef]
  2. Mahacek, P.; Kitts, C.A.; Mas, I. Dynamic Guarding of Marine Assets Through Cluster Control of Automated Surface Vessel Fleets. IEEE/ASME Trans. Mechatron. 2012, 17, 65–75. [Google Scholar] [CrossRef]
  3. Raboin, E.; Švec, P.; Nau, D.S.; Gupta, S.K. Model-predictive asset guarding by team of autonomous surface vehicles in environment with civilian boats. Auton. Robot. 2014, 38, 261–282. [Google Scholar] [CrossRef]
  4. Kim, Y.H.; Lee, S.W.; Yang, H.S.; Shell, D.A. Toward autonomous robotic containment booms: Visual servoing for robust inter-vehicle docking of surface vehicles. Intell. Serv. Robot. 2012, 5, 1–18. [Google Scholar] [CrossRef]
  5. Zoss, B.M.; Mateo, D.; Kuan, Y.K.; Tokić, G.; Ghamanbaz, M.; Goh, L.; Vallegra, F.; Bouffanais, R.; Yue, D.K.P. Distributed system of autonomous buoys for scalable deployment and monitoring of large waterbodies. Auton. Robot. 2018, 42, 1669–1689. [Google Scholar] [CrossRef]
  6. Chamanbaz, M.; Notarstefano, G.; Bouffanais, R. A randomized distributed ellipsoid algorithm for uncertain feasibility problems. In Proceedings of the 2017 IEEE 56th Annual Conference on Decision and Control, Melbourne, Australia, 12–15 December 2017; pp. 1305–1310. [Google Scholar]
  7. Zhang, L.J.; Jia, H.-M.; Qi, X. NNFFC-adaptive output feedback trajectory tracking control for a surface ship at high speed. Ocean. Eng. 2011, 38, 1430–1438. [Google Scholar] [CrossRef]
  8. Li, J.; Ren, W.; Xu, S. Distributed containment control with multiple dynamic leaders for double-integrator dynamics using only position measurements. IEEE Trans. Autom. Control 2012, 57, 1553–1559. [Google Scholar] [CrossRef]
  9. Ihle, I.A.F.; Skjetne, R.; Fossen, T.I. Output feedback control for maneuvering systems using observer backstepping. In Proceedings of the 13th Mediterranean Conference on Control and Automation, Limassol, Cyprus, 27–28 June 2005; pp. 1512–1517. [Google Scholar]
  10. Dai, S.L.; Wang, M.; Wang, C.; Li, L. Learning from adaptive neural network output feedback control of uncertain ocean surface ship dynamics. Int. J. Adapt. Control. Signal Process. 2014, 28, 341–365. [Google Scholar] [CrossRef]
  11. Wei, H.; Zhao, Y.; Sun, C. Adaptive Neural Network Control of a Marine Vessel With Constraints Using the Asymmetric Barrier Lyapunov Function. IEEE Trans. Cybern. 2017, 47, 1641–1651. [Google Scholar]
  12. Du, J.; Hu, X.; Liu, H.; Chen, C.L. Adaptive Robust Output Feedback Control for a Marine Dynamic Positioning System Based on a High-Gain Observer. IEEE Trans. Neural Netw. Learn. Syst. 2015, 26, 2775–2786. [Google Scholar] [CrossRef]
  13. Xia, G.; Shao, X.; Zhao, A. Robust Nonlinear Observer Furthermore, Observer-Backstepping Control Design For Surface Ships. Asian J. Control. 2015, 17, 1377–1393. [Google Scholar] [CrossRef]
  14. Zhang, J.; Yu, S.; Yan, Y. Fixed-time output feedback trajectory tracking control of marine surface vessels subject to unknown external disturbances and uncertainties. ISA Trans. 2019, 93, 145–155. [Google Scholar] [CrossRef]
  15. Li, X.; Shi, P.; Wang, Y. Distributed cooperative adaptive tracking control for heterogeneous systems with hybrid nonlinear dynamics. Nonlinear Dyn. 2018, 95, 2131–2141. [Google Scholar] [CrossRef]
  16. Xia, G.; Sun, C.; Zhao, B.; Xia, X.; Sun, X. Neuroadaptive Distributed Output Feedback Tracking Control for Multiple Marine Surface Vessels With Input and Output Constraints. IEEE Access 2019, 7, 123076–123085. [Google Scholar] [CrossRef]
  17. Peng, Z.; Wang, D.; Liu, H.H.T.; Sun, G.; Wang, H. Distributed robust state and output feedback controller designs for rendezvous of networked autonomous surface vehicles using neural networks. Neurocomputing 2013, 115, 130–141. [Google Scholar] [CrossRef]
  18. Lin, X.; Wang, Y.; Liu, Y. Neural-network-based robust terminal sliding-mode control of quadrotor. Asian J. Control 2022, 24, 427–438. [Google Scholar] [CrossRef]
  19. Yan, Z.; Wang, M.; Xu, J. Global adaptive neural network control of underactuated autonomous underwater vehilces with parametric modeling uncertainty. Asian J. Control 2019, 21, 1342–1354. [Google Scholar] [CrossRef]
  20. Seshagiri, S.; Khalil, H.K. Output feedback control of nonlinear systems using RBF neural networks. IEEE Trans. Neural Netw. 2000, 11, 69–79. [Google Scholar] [CrossRef] [Green Version]
  21. Lu, Y.; Zhang, G.; Sun, Z.; Zhang, W. Adaptive cooperative formation control of autonomous surface vessels with uncertain dynamics and external disturbances. Ocean Eng. 2018, 167, 36–44. [Google Scholar] [CrossRef]
  22. Tee, K.P.; Ge, S.S. Control of fully actuated ocean surface vessels using a class of feedforward approximators. IEEE Trans. Control. Syst. Technol. 2006, 14, 750–756. [Google Scholar] [CrossRef]
  23. Kahveci, N.E.; Ioannou, P.A. Adaptive steering control for uncertain ship dynamics and stability analysis. Automatica 2013, 49, 685–697. [Google Scholar] [CrossRef]
  24. Chen, M.; Tao, G.; Jiang, B. Dynamic Surface Control Using Neural Networks for a Class of Uncertain Nonlinear Systems With Input Saturation. IEEE Trans. Neural Netw. Learn. Syst. 2015, 29, 2086–2097. [Google Scholar] [CrossRef] [PubMed]
  25. Huang, B.; Song, S.; Zhu, C.; Li, J.; Zhou, B. Finite-time distributed formation control for multiple unmanned surface vehicles with input saturation. Ocean Eng. 2021, 233, 109158. [Google Scholar] [CrossRef]
  26. Zhou, W.; Wang, Y.; Ahn, C.K.; Cheng, J.; Chen, C. Adaptive Fuzzy Backstepping-Based Formation Control of Unmanned Surface Vehicles With Unknown Model Nonlinearity and Actuator Saturation. IEEE Trans. Veh. Technol. 2020, 69, 14749–14764. [Google Scholar] [CrossRef]
  27. Xia, G.; Xia, X.; Zhao, B.; Sun, C.; Sun, X. Distributed Tracking Control for Connectivity-Preserving and Collision-Avoiding Formation Tracking of Underactuated Surface Vessels with Input Saturation. Appl. Sci. 2020, 10, 3372. [Google Scholar] [CrossRef]
  28. Do, K.D. Formation Tracking Control of Unicycle-Type Mobile Robots With Limited Sensing Ranges. IEEE Trans. Control. Syst. Technol. 2008, 16, 527–538. [Google Scholar] [CrossRef] [Green Version]
  29. Do, K.D. Practical formation control of multiple underactuated ships with limited sensing ranges. Robot. Auton. Syst. 2011, 59, 457–471. [Google Scholar] [CrossRef]
  30. Hu, H.; Yoon, S.Y.; Lin, Z. Coordinated Control of Wheeled Vehicles in the Presence of a Large Communication Delay Through a Potential Functional Approach. IEEE Trans. Intell. Transp. Syst. 2014, 15, 2261–2272. [Google Scholar] [CrossRef]
  31. Yoo, S.J.; Kim, T.-H. Distributed formation tracking of networked mobile robots under unknown slippage effects. Automatica 2015, 54, 100–106. [Google Scholar] [CrossRef]
  32. Huang, Y.; Ding, H.; Zhang, Y.; Wang, H.; Cao, D.; Xu, N.; Hu, C. A Motion Planning and Tracking Framework for Autonomous Vehicles Based on Artificial Potential Field Elaborated Resistance Network Approach. IEEE Trans. Ind. Electron. 2020, 67, 1376–1386. [Google Scholar] [CrossRef]
  33. Tee, K.P.; Ge, S.S. Control of nonlinear systems with partial state constraints using a barrier Lyapunov function. Int. J. Control 2011, 84, 2008–2023. [Google Scholar] [CrossRef]
  34. Ghommam, J.; Saad, M. Adaptive Leader–Follower Formation Control of Underactuated Surface Vessels Under Asymmetric Range and Bearing Constraints. IEEE Trans. Veh. Technol. 2018, 67, 852–865. [Google Scholar] [CrossRef]
  35. Park, B.S.; Yoo, S.J. An Error Transformation Approach for Connectivity-Preserving and Collision-Avoiding Formation Tracking of Networked Uncertain Underactuated Surface Vessels. IEEE Trans. Cybern. 2019, 49, 2955–2966. [Google Scholar] [CrossRef]
  36. Yu, J.; Zhao, L.; Yu, H.; Lin, C. Barrier Lyapunov functions-based command filtered output feedback control for full-state constrained nonlinear systems. Automatica 2019, 105, 71–79. [Google Scholar] [CrossRef]
  37. Zhang, X.; Liu, L.; Feng, G. Leader–follower consensus of time-varying nonlinear multi-agent systems. Automatica 2015, 52, 8–14. [Google Scholar] [CrossRef]
  38. Liu, Y.; Lu, S.; Tong, S.; Chen, X.; Chen, C.L.P.; Li, D. Adaptive control-based Barrier Lyapunov Functions for a class of stochastic nonlinear systems with full state constraints. Automatica 2018, 87, 83–93. [Google Scholar] [CrossRef]
  39. Fossen, T.I.; Strand, J.P. Passive nonlinear observer design for ships using Lyapunov methods: Full-Scale experiments with a supply vessel. Automatica 1999, 35, 3–16. [Google Scholar] [CrossRef]
  40. Fossen, T.I. How to incorporate wind, waves and ocean currents in the marine craft equations of motion. In Proceedings of the 9th IFAC Conference on Manoeuvring and Control of Marine Craft, Arenzano, Italy, 19–21 September 2012; pp. 126–131. [Google Scholar]
  41. Park, B.S.; Yoo, S.J.; Park, J.B.; Choi, Y.H. Adaptive neural sliding mode control of nonholonomic wheeled mobile robots with model uncertainty. IEEE Trans. Control. Syst. Technol. 2009, 17, 207–214. [Google Scholar] [CrossRef]
  42. Skjetne, R.; Kokotović, P.V. Adaptive maneuvering, with experiments, for a model ship in a marine control laboratory. Automatica 2005, 41, 289–298. [Google Scholar] [CrossRef]
Figure 1. Earth-fixed frame and body-fixed frame.
Figure 1. Earth-fixed frame and body-fixed frame.
Jmse 10 00827 g001
Figure 2. The structure diagram of swarm control for USVs.
Figure 2. The structure diagram of swarm control for USVs.
Jmse 10 00827 g002
Figure 3. Directed communication topology.
Figure 3. Directed communication topology.
Jmse 10 00827 g003
Figure 4. Trajectories of six USVs.
Figure 4. Trajectories of six USVs.
Jmse 10 00827 g004
Figure 5. Tracking errors of six USVs.
Figure 5. Tracking errors of six USVs.
Jmse 10 00827 g005
Figure 6. Control input of six USVs.
Figure 6. Control input of six USVs.
Jmse 10 00827 g006
Figure 7. Velocity estimation performance of NASO with solid lines for the real states and dashed lines for their estimation.
Figure 7. Velocity estimation performance of NASO with solid lines for the real states and dashed lines for their estimation.
Jmse 10 00827 g007
Figure 8. Total unknown function estimation performance of NASO with solid lines for the real states and dashed lines for their estimation.
Figure 8. Total unknown function estimation performance of NASO with solid lines for the real states and dashed lines for their estimation.
Jmse 10 00827 g008
Figure 9. Tracking errors of six USVs using NADSC without ADS.
Figure 9. Tracking errors of six USVs using NADSC without ADS.
Jmse 10 00827 g009
Figure 10. Control input of six USVs using NADSC without ADS.
Figure 10. Control input of six USVs using NADSC without ADS.
Jmse 10 00827 g010
Figure 11. Control input of six USVs using NADSC without ADS and input constraint.
Figure 11. Control input of six USVs using NADSC without ADS and input constraint.
Jmse 10 00827 g011
Table 1. Notations and variables used in this paper.
Table 1. Notations and variables used in this paper.
VariableDefinition
A B Set of elements belonging to A but not belonging to B
· Absolute value of a scalar
· Euclidean norm
R m × n m × n dimensional Euclidean space
( · ) T Transpose of a matrix
( · ) 1 Inverse of a matrix
Kronecker product of matrix
A R n × n Adjacency matrix defined as A = [ a i j ] n × n with a i j = a j i
d i Defined as d i = j = 1 n a i j
D R n × n Degree matrix defined as D = diag { d 1 , d 2 , , d n }
diag { d i } A block-diagonal matrix with d i being the ith diagonal element
L R n × n Laplacian matrix defined as L = D A
B R n × n Information exchange matrix defined as B = L A 0
λ min ( · ) Minimum of eigenvalues of a matrix
λ max ( · ) Maximum of eigenvalues of a matrix
I n n × n dimensional identity matrix
0 n n × n dimensional zero matrix
Table 2. The parameters of six USVs.
Table 2. The parameters of six USVs.
μ i Parameters/ ( m , m , rad ) η i Parameters/ ( m , m , rad )
μ 1 [ 1.2 , 0 , 0 ] T η 1 [ 1.4 , 0.2 , π / 5 ] T
μ 2 [ 0.6 , 1.2 , 0 ] T η 2 [ 0.7 , 1.0 , π / 5 ] T
μ 3 [ 0.6 , 1.2 , 0 ] T η 3 [ 0.7 , 1.0 , π / 5 ] T
μ 4 [ 1.2 , 0 , 0 ] T η 4 [ 1.3 , 0.2 , π / 5 ] T
μ 5 [ 0.6 , 1.2 , 0 ] T η 5 [ 0.5 , 1.5 , π / 5 ] T
μ 6 [ 0.6 , 1.2 , 0 ] T η 6 [ 0.7 , 1.5 , π / 5 ] T
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xia, G.; Sun, X.; Xia, X. Swarm Control for Connectivity-Preserving and Collision-Avoiding Unmanned Surface Vehicles Subject to Multiple Constraints. J. Mar. Sci. Eng. 2022, 10, 827. https://doi.org/10.3390/jmse10060827

AMA Style

Xia G, Sun X, Xia X. Swarm Control for Connectivity-Preserving and Collision-Avoiding Unmanned Surface Vehicles Subject to Multiple Constraints. Journal of Marine Science and Engineering. 2022; 10(6):827. https://doi.org/10.3390/jmse10060827

Chicago/Turabian Style

Xia, Guoqing, Xianxin Sun, and Xiaoming Xia. 2022. "Swarm Control for Connectivity-Preserving and Collision-Avoiding Unmanned Surface Vehicles Subject to Multiple Constraints" Journal of Marine Science and Engineering 10, no. 6: 827. https://doi.org/10.3390/jmse10060827

APA Style

Xia, G., Sun, X., & Xia, X. (2022). Swarm Control for Connectivity-Preserving and Collision-Avoiding Unmanned Surface Vehicles Subject to Multiple Constraints. Journal of Marine Science and Engineering, 10(6), 827. https://doi.org/10.3390/jmse10060827

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