1. Introduction
The collaborative control of multi-unmanned surface vehicle (USV) systems has been increasingly applied in military fields in recent years, such as coordinated escort, coordinated mine sweeping, and coordinated encirclement of enemy ships [
1]. Similarly, in civilian fields, it can expand the scope of marine operations, improve operational efficiency, and complete complex tasks that cannot be achieved by a single ship, including coordinated search and rescue, coordinated reconnaissance, long-distance supply, maritime patrol, and coordinated transportation, among others. The formation problem is an important issue in the control of multi-USV systems [
2,
3,
4], which requires USVs to maintain a certain relative spatial position while in motion. This plays an important role when USV clusters complete common tasks such as oil and gas exploration and maritime search and rescue [
5].
The field of formation control application is broad, which has led to its rapid development [
6,
7,
8]. Formation control can be classified into two types of problems [
9] based on whether the formation spatial position changes over time: (1) formation problem, which controls multiple agents to achieve a pre-specified formation that does not change over time, and (2) formation tracking problem, which requires agents to track a pre-specified formation trajectory that may change over time. In this case, the formation may also change over time, and changing the formation based on actual motion trajectory is a more practical approach. Several commonly used formation methods include the leader-following method [
10], behavior-based method [
11], virtual structure method [
12], etc. A unified formation method framework based on consensus for second-order multi-agent systems was proposed in reference [
13], which makes the aforementioned formation methods special cases of this framework. The leader-following method combined with the line-of-sight (LOS) strategy has been used to study the safety formation problem of multi-USV systems considering unmodeled dynamics, environmental disturbances, input saturation, and output constraints [
14]. In actual complex maritime missions, such as cooperative surveillance, it is not enough to have multiple USVs tracking a target. Usually, multiple leader USVs form a safe area through formation, and USV followers enter the area to execute hazardous tasks [
15]. The USV leaders maintain the desired formation, and the follower USVs converge to the convex hull formed by the leader USVs. This problem is known as the formation containment problem [
16]. The leaders’ formation can be controlled by setting the distance and angle between the leader and the target from the shore, while the follower’s formation is controlled by the leader, and its relative position in the convex hull cannot be directly determined. However, for safety and the efficient completion of various tasks, the position of the follower in the convex hull is essential. Due to the constraints of the containment control law, the formation of the followers cannot be arbitrarily set by the control protocol. To design a formation for the followers, a topology structure can be used. Reference [
17] proved that the topology network enables the followers of the multi-agent system to converge to the convex hull and maintain a convex combination of some leader’s formation within the convex hull; based on this theory, two kinds of follower regular polygon formations were designed. However, to our knowledge, the relevant problems of designing formations to make the system have convergence rate advantages have not been widely studied.
In the problem of formation containment of cooperative control of multiple USVs, saving communication cost has significance. Reference [
18] applies a dynamic event-triggered mechanism (DET) to the formation containment-tracking control of discrete multi-agent systems, saving communication cost and communication resource utilization. In addition to the event-triggered method to improve communication bandwidth utilization, communication resource utilization can also be improved by designing communication topology [
17], and simplifying the communication link between USVs can also play a role in improving communication resource utilization [
17]. Therefore, this article selects a low communication resource connection method with a simple structure, low cost, and high effectiveness, tree topology [
19,
20], as the communication connection topology between various levels of USVs.
Due to its advantages such as short communication line length, low cost, ease of node expansion, convenient search, and isolation of fault paths [
19,
21,
22], tree topology networks are often used in the topology structure of tracking and formation control for various multi-agent system clusters [
23,
24]. Moreover, the tree topology structure is also a directed acyclic graph (DAG).
A directed acyclic graph (DAG) is a kind of directed graph with topological order and no directed cycle. Algebraic connectivity is a crucial index in graph theory and is defined as the smallest non-zero eigenvalue of the Laplacian matrix of the communication topology of a multi-agent system [
25,
26]. In the containment problem of multiple USVs, if there are multiple follower USVs in the convex hull, their algebraic connectivity is often related to the degree of formation density of the follower USVs in the convex hull; the larger the algebraic connectivity, the denser the follower USVs. In more complex and harsh actual sea areas, if the formation of follower USVs in the convex hull is too dense, it is easy to cause collisions between agents, and it is not conducive to obstacle avoidance, so the formation of followers in the convex hull should be as dispersed and sparse as possible. Therefore, a relatively small algebraic connectivity is needed when forming a formation for the follower USVs. In addition, algebraic connectivity can also reflect the convergence rate of the collaborative multi-USV system [
27,
28]. When the formation of follower USVs guarantees relative dispersion and safety, it is best to have the fastest possible convergence rate.
Adding a reverse edge on a DAG may decrease its algebraic connectivity [
28,
29,
30]. However, there is currently no effective and intuitive standard for judging whether adding a reverse edge to a general DAG affects its algebraic connectivity [
28]. Nevertheless, some general rules that affect the convergence rate have been summarized for some special DAGs, such as 1D chains and 2D grids of geographic areas commonly used in mobile ad hoc networks (MANETs) in the communication field [
31,
32]. The interference of the reverse edge depends only on the reverse range, but it is independent of the specific location of the reverse edge in the graph and the size of the graph [
30]. For DAGs with nodes of in-degree 1, reference [
28] defines a reverse coefficient to determine whether a reverse edge is an interfering edge by comparing it with one through calculation. Reference [
29] defined a single-chain DAG local structure as a stem; for example, the DAG in
Figure 1, a local structure like
is a stem.
The reference presented an accurate condition for the interference reverse edge of this structure, took the dominant convergence rate as a function, and quantitatively studied its monotonic relationship with the reverse range nodes’ in-degree. In addition, due to the stem presenting a network structure with a single-chain topology with topological order, as a relatively simple and connected topology, it has advantages such as simplicity, reliability, and low latency in communication, making it suitable for small and rapidly deployable networks [
33,
34,
35,
36,
37], such as the USV communication network. Therefore, it can be used for communication topology connection between followers in multi-USV system cooperative tasks.
Inspired by the above studies, this paper designed a special type of DAG communication topology network for USV clusters based on the tree-like topology structure, which is advantageous in saving communication costs. The algebraic connectivity of the communication topology was also taken into consideration, and a topology reconstruction method based on reverse edge was proposed for optimizing the formation of followers in the formation containment tracking problem. Specifically, this paper has achieved the following:
- (1)
A simple form of a USV (unmanned surface vehicle) follower network is considered in this article, and a stem structure is used to describe the communication topology between followers. As the algebraic connectivity would be decreased by adding a reverse edge to the stem structure, this property is utilized by this article to add a communication link to USV followers in order to reduce the algebraic connectivity, thereby making the follower formations more dispersed and more conducive to USV distributed security operations.
- (2)
The changes in algebraic connectivity by adding different reverse links to the stem-structured communication topology of the USV follower cluster are analyzed in this article. A method for optimal addition of a reverse link is proposed based on the above rules, considering the appropriate algebraic connectivity of the communication topology, in order to optimize the problem of slower convergence rate caused by adding a reverse link to the stem structure. The formation of the USV followers is improved by this method while minimizing the loss of convergence rate to the convex hull of the leaders.
- (3)
Compared to the more complex topology reconstruction methods that require extensive changes to the communication links, the advantage of the method described in this article is that it only involves minor modifications to the original DAG graph, i.e., adding a reverse edge to the stem structure, which can significantly change the formation of USV followers while ensuring a good convergence rate.
The structure of the rest of this paper is as follows:
Section 2 presents some preliminary knowledge to support the conclusion.
Section 3 proposes a method of adding reverse edges for a special class of communication topologies to make the formation of multiple USVs more dispersed, while minimizing the convergence rate of system loss.
Section 4 conducts some simulation experiments. Finally,
Section 5 presents the conclusion of this paper.
2. Prelimination and Problem Formulation
2.1. Graph Theory Stem and Reverse Edge
Consider a network with a number of nodes , represented by a directed graph . Here, is the set of nodes, is the set of edges, is the adjacency matrix, and its element . In particular, the existence of an edge from node j to node i is represented as , and its absence is represented as . Additionally, this article assumes that there are no circular connections where a node is connected to itself, i.e., . The in-degree of node is defined as . Let represent the in-degree of graph G. represent the in-degree of the node , in-degree means the number of the agents that connect agent . All other off-diagonal elements are 0. The Laplacian matrix of graph G is defined as .
The state represented as a single agent object is described by the following expression for the state of all agents in the network, since there are
N agents in the entire network. The
in this paper means the USV
i position coordinates:
where
. Writing the equation in vector form yields:
In order for individual agents to achieve what is known as consensus, they need to be part of a connected network. The simplest requirement for a connected network is that it has a spanning tree. That is, if a directed graph has a root node with directed paths to all other nodes, the directed graph is said to contain a spanning tree.
In addition, the matrix L has an eigenvector with a non-negative eigenvalue, and the smallest eigenvalue is 0. The non-zero smallest eigenvalue is called the algebraic connectivity, which has rich meanings in specific application contexts.
Definition 1 ([
30])
. For a network with a directed graph that contains a spanning tree, the minimum real part of the non-zero smallest eigenvalue of its Laplacian matrix, i.e., , is called the dominant convergence rate, if all eigenvalues have real parts. Below, we introduce a special network structure called DAG (directed acyclic graph).
Definition 2 ([
31])
. DAG is a directed graph without directed cycles. It consists of a set of nodes and a set of edges, where each edge connects one node to another, ensuring that there is no node i that can start from a series of edges and eventually loop back to i. Definition 3 ([
29])
. Every DAG always has a topological ordering, which is a sequence of nodes where the starting node of each edge comes before its ending node. A topological ordering is only possible for a graph if it is a DAG. However, in general, this ordering is not unique. Typical examples of DAGs include singly linked chains, grids, trees, mushrooms, etc. Among them, there is a local DAG structure called stem, which has a singly linked network structure and is located in a local region of a large network. More importantly, it can be used for local quantitative analysis and can reduce the workload of analyzing large networks.
Remark 1. The following graph in Figure 2 is a directed acyclic graph (DAG) with topological ordering. The ordering refers to the sequence of the nodes, ensuring that for each edge, the start node always precedes the end node. One topological ordering of a graph is . Definition 4 ([
30])
. For a DAG with a topological ordering of N nodes, suppose there is a path from node ℏ to node ℓ and . If there is no edge from i to j and , , the number of nodes on the stem, i.e., , is called the length of the stem. Nodes ℏ and ℓ are, respectively, called the start node and end node, and the ordered set of nodes from the start node to the end node according to the topological ordering is called the reverse range. Remark 2. The following graph in Figure 1 is a directed acyclic graph. For example, is a stem structure, where start node and the end node . Choose any : , and . The possible values for are as follows: and . By definition, there is no direct edge from i to j in the graph. Therefore, is a stem structure. The greatest advantage of the stem structure lies in its ability to perform both quantitative and qualitative analyses, as stated in the following lemma:
Lemma 1 ([
30])
. Consider a graph with N nodes and a topological ordering, where the Laplacian matrix is a lower triangular matrix. Assume that the graph has a stem with a start node h and another node l, with . The in-degree requirements for all nodes in the stem structure are as follows:The in-degree of node h is .
The in-degree of node l is .
The in-degrees of nodes from to form a set , where . In this sense, there are nodes with an incoming degree of , for , satisfying . The following relationships hold:
The impact of adding a reverse edge from node l to node h is as follows:
- (1)
Only eigenvalues in the Laplacian matrix are changed, and the remaining eigenvalues remain unchanged.
- (2)
The changed eigenvalues are determined by the following polynomial equations: The equations always have a real root . In particular, can be explicitly expressed as .
- (3)
is smaller than the remaining r roots in Equation (4). - (4)
For , when all other independent variables are fixed, the function is strictly monotonically increasing with respect to the independent variable .
- (5)
For , when other independent variables are fixed, if , the function is strictly monotonically increasing with respect to the independent variable . If , then the function is strictly monotonically decreasing with respect to the independent variable .
- (6)
If , then is the dominant convergence rate.
2.2. The Formation Containment Tracking Control and USV Model
Below is an introduction to the basic description of formation containment tracking tasks. This task essentially divides the multi-agent system into three roles: the tracking leader is responsible for determining the macroscopic motion trajectory of the entire system, the formation leader is responsible for tracking the tracking leader according to the formation set by the motion control protocol, and the followers converge into the convex hull formed by the formation leader according to the motion control protocol.
The main content of this section is to provide motion control protocols for both the formation leader and followers based on the kinematic descriptions of the three roles mentioned above, i.e., Equations (
11) and (
12). The parameters of the control protocols are determined according to Algorithm 1 in reference [
16]. This paper only uses this method to design control protocols.
In order to complete the formation containment tracking tasks, the communication connection methods for USVs with different roles have the following requirements.
Assumption 1 ([
16])
. The communication topology between the tracking leader and the formation leader satisfies the following conditions: it includes a spanning tree, and the root of the tree is the tracking leader. Assumption 2 ([
16])
. The topology relationship between the formation leader and the followers satisfies the following condition: for any follower, there is at least one formation leader that has a directed path to it. Considering the entire multi-USV system, according to the order of one tracking leader, N formation leaders, and M followers, the Laplacian matrix can be divided into , where , , , .
Based on Assumptions 1 and 2, there exists a diagonal matrix such that and . There also exists a diagonal matrix such that and . Feasible matrices and can be calculated as and . When applied to a multi-USV system, due to the fact that the research object of this paper is small unmanned boats, and the spatial distance between each USV in the formation is much greater than the length of the USV itself, the motion angle deviation of each USV can be ignored. Therefore, an approximate first-order linear point-mass motion model can be used to describe each USV.
The dynamic descriptions of the formation leaders and followers are as follows:
where the position coordinates of the single USV is represented by
and
as the control input. The dynamic description of the tracking leader USV is given by:
where
represents the position coordinates of the tracking leader USV, and
represents the control input. The time-varying non-zero control input
is unknown to all the formation leaders and followers.
Assumption 3. The position control input of the tracking leader is bounded, and there exists a positive real number γ such that .
The expected time-varying formation for the formation leader is given by
, where
is piecewise continuously differentiable, which requires that
and
are both bounded.
The followers should converge to the convex hull formed by the formation leader, and there exists a set of non-negative numbers
such that
. The specific convergence position satisfies the following equation.
The time-varying formation tracking position error between a formation leader and its neighbors is defined as:
where
if the neighbor of the formation leader includes the tracking leader; otherwise,
.
For the follower, the relative error between itself and its neighbors is defined as:
Consider the following formation containment tracking control protocol:
where
represents the time-varying formation tracking compensation input determined by
.
and
are constant gain matrices,
and
are positive constants, and
is a nonlinear function.
2.3. Topology Design
Due to the fact that neither Assumption 1 nor Assumption 2 require any specific topology between the tracking leader and the followers, in general, the tracking leader will not be connected to the followers. As for the topology between the tracking leader and the formation leader, only a spanning tree with the tracking leader as the root is required, and the formation leader’s formation, i.e., the relative position between the formation leader and the tracking leader, is determined by designing the relative position and other parameters in the control protocol. Therefore, in the third part of this paper, the topology between the tracking leader and the formation leader is ignored, and only the topology between the formation leader and the followers is designed.
Assuming that the graph to be considered is composed of leaders with indices 1, 2, …, M and followers with indices M + 1, M + 2, …, N, and their connection satisfies Assumption 2, the Laplacian matrix of the graph becomes , where and .
Then, the followers will converge to the convex hull formed by the formation leaders, and the specific convergence position is given by the following two lemmas.
Lemma 2 ([
9])
. Under Assumption 2, all eigenvalues of have positive real parts, each element of is non-negative, and the sum of each row of is equal to 1. Lemma 3 ([
9])
. Assuming every single USV is stable, which means that matrix A,B in Equation (12) satisfies that the controllability matrix is full rank, and that the graph connectivity satisfies Assumption 2, the control protocol can solve the containment control problem for the followers. Then, the specific convergent coordinates of the followers can be obtained from Equation (13).where represents the coordinates of the followers and represents the coordinates of the formation leader. The main goal of this article is to briefly describe the topology structure between the formation leader and followers in the design of the formation. Goal: To propose a topology structure based on a tree-like topology, where the followers are interconnected with a special single-chain DAG structure. The topology includes M formation leaders and followers. Using the relevant theory of algebraic connectivity affected by adding a reverse edge to the stem structure in the DAG, we design a graph that describes the communication connection between followers and formation leaders. The graph should satisfy the following requirements:
- (1)
Followers’ formation in the convex hull should be relatively dispersed, and the algebraic connectivity should lower than that of the original tree-like topology.
- (2)
Adding a reverse edge can decrease algebraic connectivity while maximizing the convergence rate of the system.
3. Main Results
First, according to the requirements for the prototype topology design in the target, we present the following topology as the prototype diagram for this paper. The communication between USV tracking leader, formation leader, and followers is carried out through a tree-shaped network, with a directed spanning tree rooted at the tracking leader and formation leader connected to several followers in sequence. The communication topology between USV followers is connected by a single-chain DAG local structure stem, the advantages of which have been briefly explained earlier in this paper and will not be repeated here. The specific topology links are described in detail in Assumptions 4 and 5.
Assumption 4. Graph has a total of n formation leaders numbered as ; at the same time, there are n followers numbered as . The followers are connected in sequence according to the topological order; a formation leader connects d followers starting from the follower whose number is larger than its own number by n in sequence according to the topological order, which construct a "minimal unit" of a formation leader, where . If a formation leader does not have enough followers to connect in the ordering, then it connects the remaining followers from the follower in the ordering. See Figure 3. Assumption 5. Graph has a total of n formation leaders numbered as ; at the same time, there are n followers numbered as . The followers are connected in sequence according to the topological order; a formation leader connects d followers as a "minimal unit", then the n "minimal units" construct the general typology structure of G″. See Figure 4. Remark 3. The structure of G″ is constructed by n "minimal units"; the ellipsis in Figure 3 and Figure 4 means there are edges connecting other followers, and d can be any constant less than n and larger than or equal to 2. For example, when , the topology structure is as in Figure 5. The entire graph G″ is a DAG, and the numbering of the nodes connected to each other follows the topological order. Therefore, the Laplacian matrix of graph G″ is a lower triangular matrix.
Due to the fact that the addition of a reverse edge is limited to followers only (the formation leader cannot receive signals from followers), according to Lemma 1, in order to study the effect of adding reverse edges on the algebraic connectivity of stem structures, it is necessary to pay attention to the in-degree of the nodes in the stem structure, i.e. nodes numbered from
to
must satisfy the requirements of Lemma 1 for in-degree. By observation, the in-degree of follower node
is
d, while the in-degree of all other follower nodes is
. If we choose follower node
as the head node and another follower node as the tail node, the in-degree between the follower nodes does not satisfy the condition of Equation (
3), i.e., the in-degree of the tail node must not be greater than the in-degree of the head node and all intermediate nodes. Obviously, this case cannot be directly analyzed by Lemma 1.
In order to use Lemma 1 smoothly, a virtual node with an in-degree of 1, denoted as , was added to the original DAG G″. This node can be directly connected to any follower node except for follower node (since node is the first among all follower nodes in the topological order, it cannot serve as the tail node). When a follower node is needed as the tail node, a virtual node with an in-degree of 1 can be added after it, and it can be used as the tail node for the reverse edge analysis. Its specific characteristics are given by Lemma 4.
Lemma 4. When , , the algebraic connectivity for the DAG graph G″ under Assumption 4. If a virtual follower node with is added, connected to any non-follower node other than node , the resulting new graph G′ remains a DAG and its Laplacian matrix is a lower triangular matrix. The topology between followers is a stem, as shown in Figure 6. Then, the algebraic connectivity for the new graph G′. Proof. Due to the Laplacian matrix
of the original DAG graph
, which is given by
, where
as
is a lower triangular matrix, it follows that the diagonal elements are the eigenvalues of
. Its algebraic connectivity is
. After adding a virtual follower to obtain graph
, the Laplacian matrix of the new graph
is given by
, where
; similarly,
is also a lower triangular matrix, and its diagonal elements are the eigenvalues of
. Its algebraic connectivity is
. □
Remark 4. As stated in item (6) of Lemma 1, if the in-degree of the tail node of a reverse edge is always 1, then the algebraic connectivity will always lie in the range of . This makes it easier to compare the effect of adding different reverse edges on the algebraic connectivity in a unified range when considering a stem structure.
Remark 5. The virtual follower node only receives information from the follower node it is connected to in the graph. Therefore, during the movement process, the trajectory of the virtual follower node coincides completely with that of the follower node it represents. In practical applications, the virtual follower node and the follower node it represents stand for the same unmanned transportation device.
Remark 6. The virtual node does not actually change the algebraic connectivity or the convergence point of the original graph. The purpose of adding this node is only to facilitate the analysis of changes in the algebraic connectivity of the graph using Lemma 1.
Since the addition of the reverse edge is limited to the followers only, according to clause (1) of Lemma 1, only the eigenvalues corresponding to the nodes within the scope of the reverse edge will change. Therefore, the eigenvalues corresponding to the leader part will not change. In the subsequent analysis of the algebraic connectivity, the formation leader can be ignored, as shown in
Figure 7.
In a stem structure with a total of t nodes, there are a total of possible reverse edges that can be added. When the tail node is located at the node with the highest index, the number of possible reverse edges that can be added is .
For convenience, this article first analyzes the case where the tail node of a reverse edge is located at the highest indexed follower node, and the initial node is one of the remaining follower nodes, the reverse edge that can reduce the algebraic connectivity and minimize the loss of algebraic connectivity is identified. The cases where other followers are the tail node of the reverse edge will be discussed later, and their sizes will be analyzed.
Now, we begin to analyze the case where the tail node of the reverse edge is the follower node with the largest topological order: connect the virtual node
to the follower node with the largest topological order
. According to Lemma 4, we take
as the tail node of the reverse edge, as shown in
Figure 8. Since the connection between followers in graph
is a stem structure, adding a reverse edge with the tail node at
and the head node selected from the set of candidate nodes
must correspond to a reverse range that is a stem structure. The impact of adding a reverse edge on the rate in a DAG with a reverse range that is a stem structure is given by the following proposition.
Proposition 1. Consider a graph G with nodes that has a topological order and whose Laplacian matrix is a lower triangular matrix. The graph contains a stem consisting of nodes, each numbered in topological order as . The node is a virtual node, and the in-degrees of these nodes have the following characteristics:
The in-degree of node is a constant d, where .
The in-degree of the virtual node is a constant 1.
The in-degree of nodes is a constant .
. Now we fix node as the tail node of the reverse range and select one of the nodes as the head node of the reverse range. Depending on the choice of the head node, we have the following conclusions about the effect of adding a reverse edge on the dominant convergence rate :
- (1)
When selecting as the starting node, the larger the serial number of the selected starting node, the smaller the dominant convergence rate . Among them, selecting node as the starting node results in the maximum algebraic connectivity.
- (2)
When the constant d satisfies and , choosing the node as the head node for the reverse edge will result in a strictly larger algebraic connectivity than choosing as the head node.
Proof. - (1)
It can be inferred from Lemma 1 that the dominant convergence rate
is a function of the first node’s in-degree
, the in-degrees of the middle nodes
, and the number of nodes with each in-degree
. Combining this corollary with the topology shown in
Figure 4, when the tail node is fixed and the head node is chosen from nodes
, the function
becomes a function of the independent variables
,
, and
, which are as follows:
. There is only one type of middle node with in-degree
, and the number of middle nodes with in-degree
varies depending on the choice of the first node. As the index of the first node selected increases and the value of
decreases, while
and
remain unchanged. Since
and
,
strictly increases monotonically with respect to
. Therefore, the dominant convergence rate decreases as
decreases.
- (2)
Before adding the reverse edge, as shown in
Figure 8, the in-degrees of node
and the other nodes in the candidate set
are all different, so we cannot directly compare their convergence rates. In order to compare the convergence rate of selecting node
with selecting other head nodes, according to the second statement of Lemma 1, the changed eigenvalues in the reverse scope can be given by the polynomial Equation
4. Applying Equation (
4) to the case where node
is selected as the head node in
Figure 4:
Applying Equation (
4) to the case where node
is selected as the head node in
Figure 4:
Equations (
14) and (
15) present the relationship between the dominant convergence rate
and the positive independent variables
d and
n in the form of implicit functions. To distinguish them, we denote the function
in the implicit function (
14) as
, and the function
in the implicit function (
15) as
, as shown in the following Equations (
16) and (
17):
Dividing both sides of Equations (
16) and (
17) by each other, we obtain:
Construct a function
with
and
as variables:
Transforming Equation (
18), we obtain:
Firstly, we study the monotonicity of
for
.
Since d and n are positive integers, it is obvious that and is decreasing for .
The second step is to study the relationship between
and
in Equation (
20). We can take their difference means subtracting one from the other:
We obtain a quadratic function with parameter d. We can find the axis of symmetry and two roots of this function and . Now we compare the values of and by discussing the possible values of d.
When
, we can obtain
. When
is in the domain of
, it can be observed that when
,
. In this case,
. By comparing this with Equation (
20), we can see that
.
When
, we can obtain the axis of symmetry of the graph of the function
as
and the point of intersection of the function with the
axis on the left side as
. Now, when
is in the domain of
, the horizontal coordinate of this intersection point is
. Obviously, if
, then
holds for all
. By comparing Equation (
20), we can see that
.
It can be seen that when the positive integer and , the convergence rate of selecting node as the head node is faster than selecting node .
□
Considering both the conclusions of Proposition 1 and a comprehensive analysis, it can be concluded that when the positive integer and , selecting the node with the latest topological order as the tail node in the stem structure can reduce the algebraic connectivity of all reverse edges and improve the formation structure. At the same time, selecting the node with the earliest topological order as the starting node for the reverse edge can enable the system to achieve the fastest convergence rate.
Next, we compare the cases of adding reverse edges when other follower nodes in the stem structure, i.e., including nodes to , are selected as the tail node.
Here, we compare the number of nodes crossed by the reverse edges and compare it with the algebraic connectivity obtained from the reverse edges with node as the tail node and nodes to as the starting nodes. The impact of adding a reverse edge on algebraic connectivity is largely determined by the range of nodes crossed by the reverse edge. Since the stem structure is a single-chain structure, the range of nodes crossed can only be a single chain. The specific variables include the in-degree of the starting node, the in-degree of the intermediate nodes, and the number of intermediate nodes with each in-degree. In the model studied in this paper, there is only one type of in-degree for intermediate nodes crossed by the reverse edge, which is . When the number of crossed nodes is the same, the number of intermediate nodes with this unique in-degree is also the same. Therefore, the choice of the tail node for the reverse edges does not affect the size of algebraic connectivity, and the difference in algebraic connectivity is only determined by the difference in the in-degree of the starting node,.
In the model analyzed in this paper, the difference in the in-degree mainly depends on whether the node with the earliest topological order is selected as the starting node. When node is the starting node, the in-degree of the starting node is , while for other nodes, the in-degree of the starting node is . Next, we compare the impact of adding reverse edges to algebraic connectivity when the number of crossed nodes is the same for different starting nodes.
When the number of crossed nodes is , there are only two possible algebraic connectivity values: when the starting node is the follower node with the smallest index, the in-degree of the starting node is , and the algebraic connectivity is . For all other nodes as starting nodes, the in-degree of the starting node is , and the algebraic connectivity obtained is . To compare the size relationship between these two algebraic connectivity values, the following proposition is given:
Corollary 1. The conditions are consistent with Proposition 1. When the number of nodes crossed by the reverse edge is the same t, the tail node has an in-degree of 1, and the two algebraic connectivity values with a difference of 1 in the in-degree of the head node, i.e., and , is strictly greater than .
Proof. We use point (4) of Lemma 1 and conclude that: for , when the other independent variables are fixed, the function is strictly monotonically increasing with respect to the independent variable . Therefore, when the other independent variables are fixed, is strictly monotonically increasing with respect to the in-degree of the first node, so is strictly greater than . □
Combining all the above conclusions, adding a reverse edge on the stem structure of the followers can improve their formation within the convex hull. From the perspective of achieving the maximum convergence rate, selecting a reverse edge in which the first node is the follower node with the smallest index and the last node is the follower node with the largest index can achieve the maximum convergence rate while ensuring improvement in the formation.
Remark 7. In graph , selecting the reverse edge with the first node as has a faster convergence rate than selecting other follower nodes as the first node. In addition, by increasing d, the algebraic connectivity of graph G can be made closer to the algebraic connectivity of graph , allowing the system to achieve almost the same convergence rate as before adding the reverse edge. However, as the constant d increases, the formation of followers becomes denser. Therefore, d should not be too large.
Corollary 2. For the follower USV cluster with the communication topology shown in Figure 8, adding such a reverse communication link can improve the formation of the follower fleet while maximizing the dominant convergence rate closest to before adding the reverse edge: the head node of the reverse link is the follower node with the smallest topological order, and the tail node is the follower node with the largest topological order. Proof. This conclusion can be easily proved by the Proofs of Proposition 1 and Corollary 1. □
4. Simulation
First, we obtain the parameters of the USV control protocol based on Algorithm 1 in [
16]. While completing the formation containment tracking task, we pay attention to observing the formation of the USV followers. The specific results are as follows:
When
and
, we have the following designed topology, where USV 0 is the tracking leader; USV 1, 2, 3, 4, and 5 are the formation leaders; USV 6, 7, 8, 9, and 10 are the followers; and USV 11 is a fictitious follower. As shown in
Figure 9, the connection between the tracking leader and the formation leaders is represented by dashed lines because the analysis of algebraic connectivity does not consider their connection relationship.
Next, based on the given topology, we design the parameters. In the first step, we provide the formation vector
for the formation leaders, and plan to maintain a regular pentagon formation centered around the tracking leader, i.e.,
. It can be clearly solved that
. In the second step, we provide the motion parameters of each USV particle model:
,
; it is not difficult to solve for the positive definite matrix
, and it is obvious that the system
is stable. In the third step, it follows from Lemma 3 in Reference [
16] that the positive definite matrix
,
,
,
. It can be calculated that
,
,
, and
, taking parameters
and
. From the fourth step, take the parameters
,
, and
. Then, the multi-USV system can accomplish formation containment tracking tasks, and each agent is relatively stationary to each other, maintaining a certain spatial relative position, as shown in
Figure 10.
This paper proposes a method of adding a reverse edge to optimize the convergence position of followers in the formation. To demonstrate the effectiveness of this method, we use the preset formation with an added reverse edge shown on the left side of
Figure 11 as a comparative experiment. The preset formation is presented as a pentagon similar to the leader’s formation, and it should be noted that, according to Remark 6, the convergence positions of the virtual node 11 and follower node 10 always remain consistent. Based on this, we remove the virtual node 11, connect follower 10 directly to follower 6, and add a reverse edge to adjust the convergence positions of other followers. By doing so and comparing it with the formation of followers on the right side of
Figure 11 without the addition of reverse edges, we can obtain the following advantages:
- (1)
The distance between each USV is relatively far, which can reduce the possibility of collision between them.
- (2)
After adding reverse edges, the convergence positions of the followers are more evenly distributed and maintain a fixed angular difference with the formation leader, which has better obstacle avoidance ability in complex sea areas. In contrast, in the case where reverse edges are not added (on the right side of
Figure 11), there are two followers, 6 and 10, who are very close to the formation leader and even appear on the edge of the convex hull. This obviously does not benefit the safety and stability of the formation.
As can be seen, the communication topology designed in this article exhibits rotational symmetry in the relationship between the convergence positions of followers with the addition of reverse edges and the position of the formation leader. Thus, the formation of the followers is similar to that of the formation leader. If the formation of the leader is in the shape of a regular polygon, the followers can converge into a similar polygon shape inside the convex hull uniformly.
Next, we ignore the tracking leader and compare the effects of adding different reverse edges to the stem structure of the original follower network on the algebraic connectivity that characterizes the convergence rate.
This article proposes a method of adding a reverse edge to optimize the formation of the team. As shown in
Figure 12, in the follower network, we choose any two nodes on the stem and take the node with a smaller index as the head node while artificially setting a node with an in-degree of 1 as the tail node of the reverse edge after the node with a larger index. This approach results in an algebraic connectivity of
. When we add reverse edges, the algebraic connectivity decreases.
Figure 13 shows the effects of different reverse edges on the algebraic connectivity.
Considering the follower network part of the communication topology, as shown in
Figure 12, the network structure appears as a stem. We choose any two points on the stem, and according to the topological order, the node with a smaller index is taken as the head node. For ease of analysis, a virtual node with an in-degree of 1 is artificially added as the tail node of the reverse edge after the node with a larger index, in accordance with the requirements of Proposition 1. This makes the algebraic connectivity
. When the reverse edge is added, the algebraic connectivity decreases. The change in the algebraic connectivity with the addition of different reverse edges is shown in
Figure 13.
In
Figure 13, the vertical axis represents the algebraic connectivity of the communication topology, and the horizontal axis represents the first node number of the reverse edge where
, and
, respectively, represent a node with an in-degree of 1 that is virtually added as the tail node of the reverse edge after nodes 10, 9, 8, and 7.
By observing
Figure 13 and only looking at the images where the first node of the reverse edge is 7, 8, or 9, the algebraic connectivity after adding the reverse edge can be expressed as the function
, where
d is a positive constant representing the unique in-degree of the middle node in the reverse edge, and
r represents the number of nodes in the middle node with an in-degree of
d. We can verify some conclusions:
- (1)
It can be seen that when , the larger the value of r, the smaller the algebraic connectivity .
- (2)
The impact of the reverse edge on the algebraic connectivity is only related to the range covered by the reverse edge and the in-degree of the nodes within that range, and is independent of the location of the reverse edge.
Remark 8. The fact that the location of the reverse edge does not affect the algebraic connectivity can be illustrated by the following example: for instance, when the virtual node after 8 is connected to 7, the virtual node after 9 is connected to 8, and the virtual node after 10 is connected to 9, even though the first and last nodes are not the same; their algebraic connectivity is the same because the length covered by the reverse edge and the in-degree of the nodes it covers are the same.
The following experiments compare the convergence rates when the tail node of the reverse edge is fixed and the first node is either the node with the earliest topological order or the follower node with the second earliest topological order. Comparing the algebraic connectivity function values for the two cases is not straightforward because
for both cases has a constant in-degree
d for the middle node, but different independent variables, namely the in-degree of the first node
and the number of middle nodes with in-degree
d:
r(which is the number of middle nodes in this structure), making direct comparison difficult.To investigate this problem more generally, we consider the difference between the algebraic connectivity of the follower node with the smallest sequence number and the follower node with the second smallest sequence number, when the unique in-degree
d of the middle node in the reverse edge is in the range
and the number of middle nodes
r is in the range
. As shown in
Figure 14, when
d is close to 1, the difference between the two is negative. However, when
, the difference is greater than zero, and as
d or
r increases, the difference between the two algebraic connectivity values approaches zero, but always remains strictly greater than zero.
Therefore, it can be verified that when , adding a fixed tail node reverse edge and selecting the reverse edge of the minimum topological order node as the first node will lead to a strictly higher dominant convergence rate than selecting the reverse edge of the second smallest topological order node as the first node.
The result in Proposition 1 can therefore be easily verified.