Next Article in Journal
Deep Reinforcement Learning Based Decision Making for Complex Jamming Waveforms
Previous Article in Journal
Co-Production of Dimethyl Carbonate, Dimethoxymethane and Dimethyl Ether from Methanol: Process Design and Exergy Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Collaborative Optimal Formation Control for Heterogeneous Multi-Agent Systems

College of Computer and Control Engineering, Qiqihar University, Qiqihar 161000, China
*
Author to whom correspondence should be addressed.
Entropy 2022, 24(10), 1440; https://doi.org/10.3390/e24101440
Submission received: 14 August 2022 / Revised: 25 September 2022 / Accepted: 28 September 2022 / Published: 10 October 2022

Abstract

:
In this paper, the distributed optimal control method is used to study the cooperative formation of heterogeneous multi-agents in the air–ground environment. The considered system consists of an unmanned aerial vehicle (UAV) and an unmanned ground vehicle (UGV). The optimal control theory is introduced into the formation control protocol, the distributed optimal formation control protocol is designed, and the stability is verified by graph theory. Furthermore, the cooperative optimal formation control protocol is designed, and the stability is analyzed using a block Kronecker product and matrix transformation theory. Through the comparison of simulation results, the introduction of optimal control theory shortens the formation time of the system and accelerates the convergence speed of the system.

1. Introduction

With the rapid development of science and technology, the world’s military powers attach great importance to the cooperation capability of an unmanned combat system. In recent decades, air–ground heterogeneous unmanned combat systems, which consist of unmanned aerial vehicles (UAVs) and unmanned ground vehicles (UGVs), have been favored by military powers due to their fast response speed, strong communication capability, strong payload capacity, and high target reconnaissance accuracy [1,2,3].
In the field of multi-agent systems, cooperative control has received extensive attention and research. Examples of such applications include the field of robot collaboration, UAV formation, cooperative transport, and combat reconnaissance [4]. Formation control is an application hotspot in the field of distributed cooperation. In general, formation control can be divided into two categories according to the presence or absence of leader agents: leader and leaderless [5,6,7,8]. Reference [9] used the leader–follower method to complete the trajectory tracking task of the UAV–UGV system but did not cooperate in the process of completing the formation. Furthermore, Reference [10] studied the cooperation problem of the UAV–UGV system by improving the artificial physics approach but did not study the formation problem. Reference [11] studied the time-varying formation control of cooperative heterogeneous multi-agent systems and combined formation with cooperation.
In the cooperative reconnaissance and cooperative strike of the unmanned combat system, based on forming an established formation, the unmanned combat system must also identify the complex and volatile battlefield environments and quickly cross barriers, such as fences and fortifications. Therefore, quickly making the whole formation reach the desired state is also an important concern of the formation problem. Reference [12] used a virtual-structure-based approach and multiple-impedance control to achieve the optimal formation of three mobile robots, and the mobile robots carried out the cooperative formation. However, this study was based on the study of homogeneous multi-agents. In Reference [13], the leader–follower strategy and the virtual leader strategy were integrated into an optimal control framework to study the optimal formation of multiple UAVs. However, this study did not investigate the cooperation of multiple UAVs. For the optimal formation of heterogeneous multi-agent, there are few related research results. Reference [14] is based on using reinforcement learning methods to achieve the optimal formation of heterogeneous multi-agent systems but does not study cooperative control. It should be noted that optimal formation control alone can only solve a relatively limited number of problems, and cooperative optimal formation control is still an open problem.
In addition, most of the existing research results are based on the same dynamic model, namely the homogeneous agent model [15,16,17,18]. Compared with homogeneous multi-agent systems, multi-agent systems composed of heterogeneous dynamic models are more flexible in practical applications. Therefore, it is of great significance to study heterogeneous multi-agent systems. A large number of valuable research results have been obtained for the heterogeneous cooperative problem [19,20,21].
In this paper, the cooperative optimal formation problem of heterogeneous multi-agent systems is studied on the unmanned aerial vehicle and unmanned ground vehicle model. There are threefold main innovations:
Firstly, a heterogeneous modeled UAV–UGV system is proposed, a cooperative architecture of heterogeneous multi-agent systems with equal number is designed, and a Laplacian matrix of communication topology is designed. In addition, a novel block Kronecker product is used to describe the UAV–UGV system. Based on this, distributed formation control is proposed.
The second contribution is to introduce the optimal control method into the formation control protocol, design the distributed formation optimal control protocol, and prove the stability using the method of graph theory.
The third contribution is to design a cooperative formation control protocol for the air–ground system based on the heterogeneous system model so that the UAV–UGV system can achieve the cooperative formation effect. Then, the optimal control is introduced into the cooperative formation control protocol, and the cooperative optimal formation control protocol is designed, which enables the UAV–UGV system formation to quickly achieve the expected effect.

2. Preliminaries

This section mainly introduces the preliminary knowledge of unmanned aerial vehicles and unmanned ground vehicles, including the use of graph theory to describe the internal relationship of the system and the state-space equation of the UAV system and the UGV system.

2.1. Graph Theory

A weighted undirected graph G = ( V , E , A ) consists of n vertices, where V = ( v 1 , v 2 , v n ) represents the set of all vertices in the undirected graph, and each vertex represents an agent. E = { e i j = ( v i , v j ) } v × v represents the edge set between vertices, and e i j = ( v i , v j ) represents the edge from v i vertex to v j vertex; an edge connection between two vertices indicates that there is an information interaction between these two vertices. The graph is undirected if it allows two-way communication; otherwise, it is directed.   A = [ a i j ] n × n represents the adjacency matrix indicating the relationship between agents, where a i j is the weight of the side e i j = ( v i , v j ) and where the diagonal elements of the matrix A are all 0. For i , j = 1 , 2 , 3 , , n ( i j ) , if the agents v i and v j can receive information from each other, then the elements in the adjacency matrix are a i j = a j i > 0 ; otherwise, the element in the adjacency matrix is 0 .
In an undirected graph, the degree represents the number of neighbors of a node, that is, the number of edges per node. D = d i g { d 1 , , d n } of undirected graph G is a diagonal matrix with d i = j = 1 n a i j .Then, the Laplacian matrix of G is defined as L = D A, which has at least one zero eigenvalue with 1 = [ 1 , 1 , 1 ] T as its corresponding right eigenvector. In addition, L has exactly one zero eigenvalue if and only if the directed graph G contains a directed spanning tree.

2.2. UGV Dynamics Model

Single UGV motion model:
{ p ˙ g i = v g i v ˙ g i = u g i
where   p g i = [ p g i x , p g i y , p g i z ] T stands for the position in ground space, v g i = [ v g i x , v g i , y , v g i z ] T is the velocity in the direction p g i , and u g i = [ u g i x , u g i y , u g i z ] T represents the input of agent i . If there are k UGVs, the above formula is converted to the states:
X ˙ G = A G X G + B G U G
where   X G = ( P G , V G ) T , P G = ( p 1 , p 2 , p 3 p k ) , p i = ( x i , y i , z i ) , i = 1 , 2 , , k ;
V G = ( v 1 , v 2 , v 3 v k ) , v i = ( v i x , v i y , v i z ) , i = 1 , 2 , , k
U G = ( u 1 , u 2 , u 3 u k ) , u i = ( u i x , u i y , u i z ) , i = 1 , 2 , , k
A G = [ 0 1 0 0 ] I k , B G = [ 0 1 ] I k .
Subscript G represents the state variable of the unmanned vehicle.
The expected formation state is h g = ( h g x , h g y , h g z ) T , the formation state is transformed into the position state, and the new error position state naturally appears, namely:   δ G = ( δ G x , δ G y , δ G Z ) T = ( P g i x h g x , P g i y h g y , P g i z h g z ) T .
Therefore, the problem of formation control becomes finding a protocol U G to drive the error vector δ G to zero, which means that
l i m t   δ g i δ g j = 0 , l i m t   = v g i = 0 .

2.3. UAV Dynamics Model

The motion model of a single UAV is:
{ x ¨ = g θ y ¨ = g ϕ z ¨ = f z / m g ϕ ¨ = M ϕ / I x θ ¨ = M θ / I y φ ¨ = M φ / I z
where   g is the acceleration of gravity; x , y , and z are the positions of the UAV in three coordinate systems; ϕ , θ , and φ are the roll angle, the pitch angle, and the yaw angle of the UAV, respectively; f z is the lift force in the direction of height; M ϕ , M θ , and M φ are the torques on the three axes of the body coordinate system; I x , I y , and I Z are the inertial matrices in the body coordinate system. For L UAVs, the above equations are converted to the state-space form as follows:
X ˙ A = A A X A + B A U A
where X A = ( P A , V A , Ω A , Ω ˙ A ) T , P A = ( p 1 , p 2 , p 3 p l ) , p i = ( x i , y i , z i ) , i = 1 , 2 , , l ;
V A = ( v 1 , v 2 , v 3 v l ) , v i = ( v i x , v i y , v i z ) , i = 1 , 2 , , l
Ω A = ( Ω 1 , Ω 2 , Ω 3 , Ω l ) , Ω i = ( g θ i , g ϕ i , 0 ) , i = 1 , 2 , , l
Ω ˙ A = ( Ω ˙ 1 , Ω ˙ 2 , Ω ˙ 3 , Ω ˙ l ) , Ω ˙ i = ( g θ ˙ i , g ϕ ˙ i , 0 ) , i = 1 , 2 , , l
A A = ( 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 ) I l , B A = ( 0 0 0 1 ) I l
U A = ( u 1 , u 2 , u 3 u l ) T , u i = ( u i x , u i y , u i z ) , i = 1 , 2 , , l .
Subscript A represents the state variable of the unmanned vehicle.
The expected formation state is h a = ( h a x , h a y , h a z ) T , the formation state is transformed into the position state, and the new error position state naturally appears, namely: δ A = ( δ A x , δ A y , δ A z ) T = ( P a i x h a x , P a i y h a y , P a i z h a z ) T .
Therefore, the problem of formation control becomes finding a protocol U A to drive the error vector δ A to zero, which means that
l i m t   δ a i δ a j = 0 , l i m t   = v a i = 0 , l i m t   Ω i = 0 , l i m t   Ω ˙ i = 0 .

2.4. Heterogeneous Multi-Agent System

To analyze heterogeneous multi-agent systems more conveniently, the UAV system and UGV system are written into the same state space and combined with the state-space model of the single agent above; the form of the heterogeneous multi-agent state-space model is defined as:
X ˙ = A X + B U
where X = ( X G T , X A T ) T , A = ( A G 0 0 A A ) , B = ( B G 0 0 B A ) ,   and   U = ( U G T , U A T ) T . The Laplace matrix is L = ( L A A L A G L G A L G G ) , where L A G , L G A represents information between heterogeneous agent systems. This paper takes the heterogeneous multi-agent system composed of three UGVs and three UAVs as the research object, and its Laplace matrix relationship is as follows:
L = ( 3 1 1 1 3 1 1 1 3 L A A 0 0 1 0 1 0 1 0 0 L G A 0 0 1 0 1 0 1 0 0 L A G 3 1 1 1 3 1 1 1 3 L G G )
The expected formation state is h = ( h a T , h g T ) T , the formation state is transformed into the position state, and the new error position state naturally appears, namely:
δ = ( δ A x , δ A y , δ A Z , δ G x , δ G y , δ G Z ) T = ( P a i x h a x , P a i y h a y , P a i z h a z , P g i x h g x , P g i y h g y , P g i z h g z ) T .
Therefore, the problem of formation control becomes finding a protocol U to drive the error vector δ to zero, which means that
l i m t δ i δ j = 0 , l i m t = v i = 0 , l i m t Ω i = 0 , l i m t Ω ˙ i = 0

3. Design of Control Protocol

To realize the formation of heterogeneous multi-agent systems of UAVs and UGVs, this section is based on the formation control protocol. Firstly, the optimal control law is applied to the single agent. Then, according to the combination of the optimal control law and the formation control protocol, a heterogeneous multi-agent system with distributed optimal formation control is realized. Finally, according to the motion equation of the heterogeneous multi-agent system, the cooperative formation control and cooperative optimal formation control of the heterogeneous multi-agent system are realized.
Lemma 1.
[22]. For an N * N Laplacian matrix L , N e L t , t > 0 is a random matrix with positive diagonal elements. If L has a unique zero eigenvalue, Rank ( N ) = N 1 , then its left eigenvector has v = [ v 1 , v 2 , v n ] T 0 and 1 N T v = 1 ,   L T v = 0 , where, t , e L t 1 N v T .

3.1. Formation Control

Formation control protocol for the UAVs:
u i a = α j N i   ( δ a j δ a i ) β v a i γ 1 Ω i γ 2 Ω ˙ i .
Formation control protocol for the UGVs:
u i g = α j N i   ( δ g j δ g i ) β v g i
where α , β , γ 1 , γ 2 represent the positive gain coefficients,
δ a i = ( δ a i x , δ a i y , δ a i z ) T ,   δ g i = ( δ g i x , δ g i y , δ g i z ) T ,   v g i = ( P ˙ g i x , P ˙ g i y , P ˙ g i z ) T ,   v a i = ( P ˙ a i x , P ˙ a i y , P ˙ a i z ) T
Ω i = ( g θ , g ϕ , 0 ) T , Ω ˙ i = ( g θ ˙ , g ϕ ˙ , 0 ) T
Protocols (9) and (10) shall be unified into the same type:
U = L d X ˜
Define the state-space form of the multi-agent system formation:
X ˜ ˙ = A X ˜ + B U
where
L d = ( α L s I 3 β I l I 3 γ 1 I l I 3 γ 2 I l I 3 0 0 0 0 0 0 α L s I 3 β I k I 3 ) ,
L S = ( 1 1 0 0 1 0 1 0 1 ) ,   X ˜ = ( δ a i T , v a i T , Ω T , Ω ˙ T , δ g i T , v g i T ) T , A = ( A A 0 0 A G ) ,   B = ( B A 0 0 B G ) , I is the identity matrix, and is the Kronecker product.
Theorem 1.
If Protocols (9) and (10) are satisfied α > 0 , β > 0 , γ 1 > 0 , γ 2 > 0 , β α , γ 1 > β , γ 2 > β , β γ 1 γ 2 > β 2 γ 2 α , Systems (2) and (5) can implement the formation defined in (3) and (6).
Proof of Theorem 1.
Substitute Formula (11) into Formula (12) to obtain: X ˜ ˙ = T d X ˜ where T d = [ 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 α L s I 3 β I l I 3 γ 1 I l I 3 γ 2 I l I 3 0 0 0 0 0 0 0 1 0 0 0 0 α L s I 3 β I k I 3 ]
According to the linear stability theorem, the parameters α , β , γ 1 , and γ 2 need to be selected so that T d has a zero eigenvalue and other eigenvalues have genuine negative parts. The parameters α and β need to meet the stability of UGV consistency, and the parameters γ 1 and γ 2 need to meet the stability of UAV consistency. After selecting parameters, T d can be converted to a Jordan standard type: T d = P J P 1 . Let v 1 T be the first row of P 1 and the left eigenvector have eigenvalue 0. Let w 1   be the first column of P and the right eigenvector have eigenvalue 0 . Therefore, v 1 T w 1 = 1 ; as time approaches infinity, the system’s state becomes:. According to Lemma 1, as time approaches infinity, Systems (2) and (5) asymptotically agree and the systems complete formation. □

3.2. Optimal Control

The solution of optimal control requires the states of all multi-agents. Before providing performance indicators, X ¯ and W ¯ are defined:
X ¯ = ( δ a i x , δ a i y , δ a i z , P ˙ a i x , P ˙ a i y , P ˙ a i z , g θ , g ϕ , 0 , g θ ˙ , g ϕ ˙ , 0 ) T , W ¯ = ( δ g i x , δ g i y , δ g i z , P ˙ g i x , P ˙ g i y , P ˙ g i z ) T Then, define the performance indicator function as:
J i = 0 [ X ¯ i T Q X ¯ i + u i a T R u i a ] d t
ω i = 0 [ w ¯ i T Q w ¯ i + u i g T γ u i g ] d t
As the UAV and the UGV are independent in different coordinate systems, the UAV weight must be set to Q = q I 12 , R = r * I 3 ,   where, q > 0 ,   r > 0 . The UGV weight must be set to T = λ * I 6 , γ = μ * I 3 , where λ > 0 , μ > 0 .
According to the optimal control theory, the optimal control law of a single agent UAV is: u a * = R 1 B A T P A X ¯ , where P A is the solution of Riccati’s in Equation (15)
A A T P A + P A A A P A B A R 1 B A T P A + Q = 0
The optimal control law for a single agent UGV is: u g * = γ 1 B G T P G w ¯ , where P G is the solution of Riccati’s in Equation (16)
A G T P G + P G A G P G B G γ 1 B G T P G + T = 0
Through the above calculation, the optimal control law u a * can be obtained. Let K = R 1 B A T P A , the dimension of matrix K is 3 12 , K expressed as K = [ k 1 , k 2 , k 3 , k 4 ] I 3 . Similarly, the optimal control law u g * can also be solved, let G = γ 1 B G T P G the dimension of matrix G is 3 6 , G expressed as G = [ g 1 , g 2 ] I 3 .

3.3. Distributed Optimal Formation Control

All UAVs have the same dynamics model, so all UAVs are homogeneous multi-agents. Similarly, all UGVs are homogeneous multi-agents. Therefore, optimal control laws can be extended to the formation control of UAV and UGV multi-agent systems.
u i a * = k 1 j N i   ( δ a j δ a i ) k 2 v a i k 3 Ω i k 4 Ω ˙ i
u i g * = g 1 j N i   ( δ g j δ g i ) g 2 v g i
Define the multi-agent system to be optimized:
X ˜ ˙ = A X ˜ + B U
where X ˜ = ( δ a i T , v a i T , Ω T , Ω ˙ T , δ g i T , v g i T ) T , A = ( A A 0 0 A G ) ,   B = ( B A 0 0 B G ) , k 1 , k 2 , k 3 , and k 4 are derived from the matrix K . g 1 and g 2 are derived from the matrix G .
Theorem 2.
If the unmanned ground vehicle system in (2) and the unmanned aerial vehicle system in (5) use Protocols (17) and (18), respectively, the formation can be completed, and Performance Functions (13) and (14) can be optimized.
Proof of Theorem 2.
Protocols (17) and (18) shall be unified into the same type:
U = L l X ˜
where
L l = ( k 1 L s I 3 k 2 I l I 3 k 3 I l I 3 k 4 I l I 3 0 0 0 0 0 0 g 1 L s I 3 g 2 I k I 3 ) ,
L S = ( 1 1 0 0 1 0 1 0 1 ) ,   X ˜ = ( δ a i T , v a i T , Ω T , Ω ˙ T , δ g i T , v g i T ) T , I is the identity matrix, and is the Kronecker product.
Let U = ( U 11 U 12 U 21 U 22 ) = ( k 1 L s k 2 I l k 3 I l k 4 I l 0 0 0 0 0 0 g 1 L s g 2 I k )
Substituting   U = ( U 11 U 12 U 21 U 22 ) X ˜ into Equation (19): X ˜ ˙ = T l X ˜
where
T l = [ 0 I 0 0 0 0 0 0 I 0 0 0 0 0 0 I 0 0 k 1 L s k 2 I k 3 I k 4 I 0 0 0 0 0 0 0 I 0 0 0 0 g 1 L s g 2 I ]
Elementary row and column transformation can be taken on T l :
T l = [ I 0 0 0 0 0 0 I 0 0 0 0 0 0 0 I 0 0 0 0 0 0 I 0 0 0 k 1 L s 0 0 0 0 0 0 0 0 g 1 L s ] = ( I 1 0 0 0 0 I 2 0 k 1 L g 1 L ) = E
If Rank ( L ) = N 1 , Rank ( [ I 1 I 2 ] ) = r , it can be seen that:   Rank ( T l ) = Rank ( E ) = r + N 1 .
where T l has only zero eigenvalues. Therefore, we must select the parameters k 1 , k 2 , k 3 , k 4 , g 1 , and   g 2 so that T l has zero eigenvalue, and all other eigenvalues have negative real parts. The parameters k 1 , k 2 , k 3   and   k 4 need to meet the stability of UAV consistency, and the parameters g 1 and g 2 need to meet the stability of UGV consistency. After determining the parameters, T l can be converted to a Jordan standard type: T l = P J P 1 . Let v 1 T be the first row of P 1 and the left eigenvector have eigenvalue 0. Let w 1 be the first column of P and the right eigenvector with eigenvalue 0 . Therefore, v 1 T w 1 = 1 ; when the time approaches infinity, the system’s state becomes: lim t   X ˜ = lim t   e T l t X ˜ ( 0 ) , e T l t X ˜ ( 0 ) ( w 1 v 1 T ) X ˜ ( 0 ) ( t ) . According to Lemma 1 it is then seen that the system can reach asymptotic consensus in cases where time tends toward infinity. □

3.4. Heterogeneous Cooperative Formation Control

For UAV:
u ˜ i a = α j N i   a i j ( δ j δ i ) β a i j ( v j v i ) γ 1 Ω i γ 2 Ω ˙ i
For UGV:
u ˜ i g = α j N i   a i j ( δ j δ i ) β a i j ( v j v i )
δ i = ( δ i a T , δ i g T ) T ,   v i = ( v i a T , v i g T ) T ,   Ω i = ( g θ , g ϕ , 0 ) T , Ω ˙ = ( g θ ˙ , g ϕ ˙ , 0 ) T
In combination with the Laplace matrix, Protocols (21) and (22) are rewritten as:
U ˜ = H X ˜
H = ( α L A A I 3 α L G A I 3 β L A A I 3 β L G A I 3 γ 1 I l I 3 0 γ 2 I l I 3 0 α L A G I 3 α L G G I 3 β L A G I 3 β L G G I 3 )
Define the state-space form of the heterogeneous multi-agent formation:
X ˜ ˙ = A X ˜ + B U ˜
where A = ( A A 0 0 A G ) ,   B = ( B A 0 0 B G ) , X ˜ = ( X ¯ T , W ¯ T ) T .
Theorem 3.
if Protocols (21) and (22) meet α > 0 , β > 0 , γ 1 > 0 , γ 2 > 0 , β α , γ 1 > β , γ 2 > β ,   β γ 1 γ 2 > β 2 γ 2 α , the heterogeneous system in (7) can be achieved, and (8) is defined in the formation, then formation control is realized.
Proof of Theorem 3.
Substitute Equation (23) into Equation (24) and obtain: X ˜ ˙ = T S X ˜
where T s = [ 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 α L A A I 3 β L A A I 3 γ 1 I l I 3 γ 2 I l I 3 α L A G I 3 β L A G I 3 0 0 0 0 0 1 α L G A I 3 β L G A I 3 0 0 α L G G I 3 β L G G I 3 ] , I is the identity matrix and is the Kronecker product.
Internal stability parameters γ 1 and γ 2 should guarantee the stability of a single UAV agent. Therefore, the UAV is written as:
( P ˙ A V ˙ A Ω ˙ A Ω ¨ A ) = Γ ( P A V A Ω A Ω ˙ A ) , where Γ = ( 0 I 0 0 0 0 I 0 0 0 0 I α I β I γ 1 I γ 2 I ) .
The characteristic polynomial is: det ( S I Γ ) = | S I 0 0 0 S I 0 0 0 S I α I β I γ 1 I S + γ 2 I |
where I is the identity matrix. For a single UAV:
| S I Γ | = s 4 + γ 2 s 3 + γ 1 s 2 + β s + α
According to the Routh–Hurwitz stability criterion:
α > 0 , β > 0 , γ 1 > 0 , γ 2 > 0 , γ 1 γ 2 > β , γ 1 γ 2 β > β 2 + γ 2 2
According to the linear stability theorem, the parameter α , β , γ 1 , and γ 2 should be selected so that there is a zero eigenvalue and other eigenvalues have negative genuine part parameters. The parameters α and β must meet the stability of UGV consistency, and the parameters γ 1 and γ 2 must meet the stability of UAV consistency. After selecting the parameters, T s can be converted to a Jordan standard type: T s = P J P 1 . Let v 1 T be the first row of P 1 and the left eigenvector have eigenvalue 0. Let w 1   be the first column of P and the right eigenvector have eigenvalue 0 . Therefore, v 1 T w 1 = 1 ; when the time approaches infinity, the system’s state becomes: lim t X ˜ = lim t e T s t X ˜ ( 0 ) , e T s t X ˜ ( 0 ) ( w 1 v 1 T ) X ˜ ( 0 ) ( t ) . According to Lemma 1, when the time approaches infinity, the system in (7) asymptotically agrees; that is, the system achieves cooperative formation.□

3.5. Heterogeneous Cooperative Optimal Formation Control

UAV systems and UGV systems are heterogeneous systems. Applying the optimal control law to the heterogeneous system can be expressed as:
u ˜ i a = k 1 j N i   a i j ( δ j δ i ) k 2 a i j ( v j v i ) k 3 Ω i k 4 Ω ˙ i
u ˜ i g = g 1 j N i   a i j ( δ j δ i ) g 2 a i j ( v j v i )
where k 1 , k 2 , k 3 ,and k 4 are derived from the matrix K . g 1 and g 2 are derived from the matrix G .
In combination with the Laplace matrix, Protocols ( 25 ) and ( 26 ) are rewritten as:
U ˜ = S X ˜
S = ( k 1 L A A I 3 k 2 L A A I 3 k 3 I l I 3 k 4 I l I 3 k 1 L A G I 3 k 2 L A G I 3 g 1 L G A I 3 g 2 L G A I 3 0 0 g 1 L G G I 3 g 2 L G G I 3 )
Define the state space of heterogeneous multi-agent formation:
X ˜ ˙ = A X ˜ + B U ˜
where A = ( A A 0 0 A G ) ,   B = ( B A 0 0 B G ) ,   X ˜ = ( X ¯ T , W ¯ T ) T
Theorem 4.
Heterogeneous systems can complete formation and optimize Performance Functions (13) and (14) if they use Protocols (25) and (26), respectively.
Proof of Theorem 4.
Each UAV and UGV is defined as a group of formation units, so their number is made the same, that is, L = K; the block Laplacian matrix of the system has the same number of rows and columns.
Then, the heterogeneous system in (7) becomes:
X ˙ = A ^ X + B ^ U
where A ^ = ( 0 I 0 0 0 0 0 0 0 0 0 0 0 0 0 I 0 0 0 0 0 0 I 0 0 0 0 0 0 I 0 0 0 0 0 0 ) , B ^ = ( 0 - - - - 0 I - - - - 0 0 - - - - 0 0 - - - - 0 0 - - - - 0 0 - - - - I ) .
According to State Equation (29):
( S I A ^ - - - - B ) = ( S I 0 0 0 0 0 0 0 S 0 0 0 0 I 0 0 0 S I 0 0 0 0 0 0 0 S I 0 0 0 0 0 0 0 S I 0 0 0 0 0 0 0 S 0 I )
It is easy to discover the eigenvalues of the matrix A ^ : λ 1 = λ 2 = = λ 6 = 0 .
Verify the rank of the matrix ( S I A ^ - - - - B ) with the above eigenvalues. Let S = λ 1 = λ 2 = = λ 6 = 0 :
Rank ( S I A ^ - - - - B )   = Rank ( 0 I 0 0 0 0 0 0 0 0 0 0 0 0 I 0 0 0 0 I 0 0 0 0 0 0 0 0 I 0 0 0 0 0 0 0 0 I 0 0 0 0 0 0 0 0 0 I )   = 6, according to the rank criterion of PBH, the system state space is controllable.
Let U ˜ = ( U 11 U 12 U 21 U 22 ) = ( k 1 L A A k 2 L A A k 3 I k 4 I k 1 L A G k 2 L A G g 1 L G A g 2 L G A 0 0 g 1 L G G g 2 L G G )
Substituting   U ˜ = ( U 11 U 12 U 21 U 22 ) X ˜ into Equation (28):   X ˜ ˙ = T * X ˜
where
T * = [ 0 I 0 0 0 0 0 0 I 0 0 0 0 0 0 I 0 0 k 1 L A A k 2 L A A k 3 I k 4 I g 1 L A G g 2 L A G 0 0 0 0 0 I k 1 L G A k 2 L G A 0 0 g 1 L G G g 2 L G G ]
The T * basic determinant change:
( λ + I I 0 0 0 0 I λ I 0 0 0 2 I 0 λ I I 0 0 C 1 C 2 k 3 I λ + k 4 I g 1 L A G g 2 L A G I 0 0 0 λ I C 3 k 2 L G A I 0 0 g 1 L G G g 2 L G G ) = λ I Λ ,
where
Λ = ( I I 0 0 0 0 I 0 I 0 0 0 2 I 0 I I 0 0 C 1 C 2 k 3 I k 4 I g 1 L A G g 2 L A G I 0 0 0 0 I C 3 I k 2 L G A 0 0 g 1 L G G g 2 L G G )
C 1 = I k 1 L A A g 2 L A G k 2 L A A k 4 I .
C 2 = k 2 L A A k 3 I I
C 3 = I k 1 L G A g 2 L G G
Now λ I T * λ I Λ , so matrix T * is similar to Λ. Therefore, there is a nonsingular transformation matrix Q, making Λ = Q T * Q 1 , X ˜ ˙ = Λ X ˜ ;   Λ is the matrix for which the sum of each row is zero. Therefore, there is at least one zero eigenvalue. Elementary row and column transformation can be taken on T * :
T * = ( I 0 0 0 0 0 0 I 0 0 0 0 0 0 0 0 I 0 0 0 0 0 0 I 0 0 k 1 L A A k 2 L A A g 1 L A G g 2 L A G 0 0 k 1 L G A k 2 L G A g 1 L G G g 2 L G G ) = ( I 1 0 0 0 0 I 2 0 k L g L ) = E
If Rank ( L ) = N 1 , Rank ( [ I 1 I 2 ] ) = r , it can be seen that:
Rank ( T * ) = Rank ( Λ ) = Rank ( E )   = r   + N 1
where the Λ and T * have only zero eigenvalues. Therefore, we must select the parameters k 1 , k 2 , k 3 , k 4 , g 1 , and   g 2 so that T * has zero eigenvalue and all other eigenvalues have negative real parts. The parameters k 1 , k 2 , k 3 ,   and   k 4 must meet the stability of UAV consistency, and the parameters g 1 and g 2 must meet the stability of UGV consistency. After determining the parameters, T * can be converted to a Jordan standard type: T * = P J P 1 . Let v 1 T be the first row of P 1 and the left eigenvector have eigenvalue 0. Let w 1 the first column of P and the right eigenvector have eigenvalue 0 . Therefore, v 1 T w 1 = 1 ; when the time approaches infinity, the system’s state becomes: lim t X ˜ = lim t e T * t X ˜ ( 0 ) , e T * t X ˜ ( 0 ) ( w 1 v 1 T ) X ˜ ( 0 ) ( t ) . According to Lemma 1, it is then seen that the system can reach asymptotic consensus in cases where time tends toward infinity.□

4. Simulations

The formation protocol, distributed optimal formation protocol, cooperative formation protocol, and cooperative optimal formation protocol individually designed in this paper are simulated and analyzed using Matlab2016a. The effectiveness of the designed control protocol is verified via simulation. To achieve a better cooperative formation task effect in the system, the speed values set by the UAV system and the UGV system are similar.
The initial state is as follows: the position of UAV1 is ( 30 , 50 , 50 ) m , the position of UAV2 is ( 90 , 30 , 50 ) m , and the position of UAV3 is ( 60 , 30 , 15 ) m . The speed of UAV1 is ( 1 , 1 , 2 ) m / s , the speed of UAV2 is ( 1 , 2 , 1 ) m / s , and the speed of UAV3 is ( 2 , 1 , 1 ) m / s . The attitude angle of UAV1 is ( 0 , 0 , 0 ) , the attitude angle of UAV2 is ( 0 , 0 , 0 ) , and the attitude angle of UAV3 is ( 0 , 0 , 0 )
The rate of the attitude angle of UAV1 is ( 0 , 0 , 0 ) / s , the rate of the attitude angle of UAV2 is ( 0 , 0 , 0 ) / s , and the rate of the attitude angle of UAV3 is ( 0 , 0 , 0 ) / s . The setting position of UAV1 is ( 10 , 10 , 30 ) m , the setting position of UAV2 is ( 0 , 5 , 30 ) m , and the setting position of UAV3 is ( 0 , 15 , 30 ) m . The UAV and UGV systems have a set speed of ( 1 , 0 , 0 ) m / s . The position of UGV1 is ( 90 , 50 , 0 ) m , the position of UGV2 is ( 65 , 10 , 0 ) m , and the position of UGV3 is ( 30 , 20 , 0 ) m . The speed of UGV1 is ( 1 , 1 , 0 ) m / s , the speed of UGV2 is ( 2 , 2 , 0 ) m / s , and the speed of UGV3 is ( 1 , 1 , 0 ) m / s . The setting position of UGV1 is ( 10 , 10 , 0 ) m , the setting position of UGV2 is ( 0 , 5 , 0 ) m , and the setting position of UGV3 is ( 0 , 15 , 0 ) m .
The parameters are α = 0.2 , β = 1.5 , γ 1 = 5 , γ 2 = 2 , k 1 = 2.3452 , k 2 = 6.4707 , k 3 = 7.7541 , k 4 = 4.5835 , g 1 = 0.4472 , and   g 2 = 1.0461 .
Simulations of the UAV system and the UGV system using formation control and distributed optimal formation control are shown in Figure 1 and Figure 2. By comparing the two figures in Figure 1 and Figure 2, it can be found that both protocols can be used to complete triangular formation at the same time, but when distributed optimal formation control is used, it can be observed that the error between the actual position and the set value is significantly reduced in a short period of time.
When using the formation control protocol, with parameters α = 0.2 , β = 1.5 , γ 1 = 5 , and   γ 2 = 2 . Figure 1 shows the change of position coordinates of the multi-agent system and shows the actual formation position of the UAV and UGV when t = 10   s . It can be seen that there is a significant error with the set value.
When using the distributed optimal formation control protocol, with parameters k 1 = 2.3452 , k 2 = 6.4707 , k 3 = 7.7541 ,   k 4 = 4.5835 , g 1 = 0.4472 , and   g 2 = 1.0461 . Figure 2 shows the position coordinate changes of the system and shows the actual formation status of the UAV and UGV when t = 10 s. It can be observed that the system can quickly complete the triangular formation and that the error between the actual position and the set value is small.
It can be seen from Figure 3 and Figure 4 that the state of each system variable changes with time when the formation control protocol is used. It can be seen from the figures that the state of each variable in the system is stable from 40 s to 50 s.
From Figure 5 and Figure 6, it can be seen that the system uses the distributed optimal formation control protocol to change the state of each variable with time. It is observed that the system reaches stability between 30 and 40 s.
The simulation of cooperative formation protocol and cooperative optimal formation protocol of heterogeneous multi-agents is shown in Figure 7 and Figure 8.
When using the cooperative formation control protocol, with parameters α = 0.2 , β = 1.5 , γ 1 = 5 , and   γ 2 = 2 . Figure 7 shows the change of position coordinates of the system and the expected formation state of the system at the 10th second. It can be seen that the actual formation position of the system has a significant error with the set value.
When using the cooperative optimal formation control protocol, with parameter k 1 = 2.3452 , k 2 = 6.4707 , k 3 = 7.7541 , k 4 = 4.5835 , g 1 = 0.4472 , and   g 2 = 1.0461 .
Figure 8 shows the change in position coordinates of the system and the expected formation state of the system at t=10 s. It can be observed that the system can quickly complete the triangular formation and that the error between the actual position and the set value is small. In addition, the cooperative optimal formation protocol speeds up the convergence rate of the system, which is of great help to the formation time of the system.
When using the formation control and the distributed optimal formation control, the formation states of the UAV system at different times are shown in Figure 9 and Figure 10.
Figure 9 shows the formation status of the UAV system at different times under the use of the formation control protocol. It can be observed that the system can complete the triangle formation but also that there is a relative reach error with the set position in the formation completion process and that there is still an error at the fiftieth second.
Figure 10 shows the formation status of the UAV system at different times under the use of the distributed optimal formation control protocol. It can be observed that the system can quickly complete the triangle formation and that the error between the actual position and the set value is very small at the thirtieth second.
When using the formation control and the distributed optimal formation control, the formation states of the UGV system at different times are shown in Figure 11 and Figure 12.
Figure 11 shows the formation status of the UGV system at different moments when the formation control protocol is used. It can be observed that the UGV system can complete the triangle formation, but there is still a certain error with the set value in the formation completion process. At the thirtieth second, the difference between the actual position and the set value gradually decreases.
Figure 12 shows the change in the formation shape of the UGV system at different moments when the distributed optimal formation control protocol is used; it can be observed that the UGV system can quickly complete the triangle formation and that the error between the actual position and the set value is small.
When using the heterogeneous cooperative formation control and the heterogeneous cooperative optimal formation control, the formation states of the UAV system at different times are shown in Figure 13 and Figure 14.
Figure 13 shows the formation status of the UAV system at different times under the use of the heterogeneous cooperative formation control protocol. It is observed that the system can complete the triangle formation. At the thirtieth second, the difference between the actual position and the set value gradually decreases.
Figure 14 shows the formation status of the UAV system at different times under the use of the heterogeneous cooperative optimal formation control. It can be observed that the system can quickly complete the triangle formation and that the error between the actual position and the set value is very small at the tenth second.
When using the heterogeneous cooperative formation control and the heterogeneous cooperative optimal formation control, the formation states of the UGV system at different times are shown in Figure 15 and Figure 16.
Figure 15 shows the formation status of the UGV system at different times under the use of the heterogeneous cooperative formation control protocol. It can be observed that the system can complete the triangle formation. At the thirtieth second, the difference between the actual position and the set value gradually decreases.
Figure 16 shows the formation status of the UGV system at different times under the use of the heterogeneous cooperative optimal formation control. It can be observed that the system can quickly complete the triangle formation and that the error between the actual position and the set value is very small at the tenth second.
In order to fully verify the theoretical results, the structure is complicated, and the communication Laplacian matrix is as follows:
L = ( 3 1 1 0 0 1 0 0 0 1 3 1 0 1 0 0 0 0 1 1 3 0 0 1 0 0 0 0 0 1 3 1 1 0 0 0 0 1 0 1 3 1 0 0 0 1 0 0 1 1 3 0 0 0 0 0 0 0 0 1 3 1 1 0 0 0 0 1 0 1 3 1 0 0 0 1 0 0 1 1 3 )
The added set of UAV states is as follows: the position of UAV4 is ( 40 , 50 , 10 ) m , the position of UAV5 is ( 60 , 30 , 20 ) m , and the position of UAV 6   is   ( 30 , 20 , 10 ) m . The speed of UAV4   is   ( 1 , 1 , 1 ) m / s , the speed of UAV5   is   ( 2 , 2 , 1 ) m / s , and the speed of UAV6 is ( 3 , 1 , 1 ) m / s .The attitude angle of UAV4 is ( 0 , 0 , 0 ) , the attitude angle of UAV5   is   ( 0 , 0 , 0 ) , and the attitude angle of UAV6 is ( 0 , 0 , 0 ) . The setting position of UAV4 is ( 10 , 10 , 10 ) m , the setting position of UAV5 is ( 0 , 5 , 10 ) m , and the setting position of UAV6 is ( 0 , 15 , 10 ) m . The UAV systems have a set speed of ( 1 , 0 , 0 ) m / s . The simulation results are shown in Figure 17, Figure 18, Figure 19 and Figure 20.
The experiments have verified the formation control and distributed optimization formation control, as shown in Figure 17 and Figure 18, respectively. By comparing Figure 17 and Figure 18, it can be found that both protocols can be used to complete triangular formation at the same time, but when distributed optimal formation control is used, it can be observed that the error between the actual position and the set value is significantly reduced in a short period of time.
When using the formation control protocol, with parameters α = 0.2 , β = 1.5 , γ 1 = 5 , γ 2 = 2 , Figure 17 shows the change of position coordinates of the complex system and the expected formation state of the complex system. However, it can be seen that there is a significant error with the set value.
When using the distributed optimal formation control protocol, with parameters k 1 = 2.3452 , k 2 = 6.4707 , k 3 = 7.7541 ,   k 4 = 4.5835 , g 1 = 0.4472 , and   g 2 = 1.0461 . Figure 18 shows the change of position coordinates of the complex system and the expected formation state of the complex system. It can be observed that the complex system can quickly complete the triangular formation and that the error between the actual position and the set value is small.
The experiments have verified the cooperative formation control and the cooperative optimal formation control, as shown in Figure 19 and Figure 20, respectively.
When using the cooperative formation control protocol, with parameters α = 0.2 , β = 1.5 , γ 1 = 5 , and   γ 2 = 2 , Figure 19 shows the change in position coordinates of the complex system and the expected formation state of the complex system. However, it can be seen that there is a significant error with the set value.
When using the cooperative optimal formation control protocol, with parameters k 1 = 2.3452 , k 2 = 6.4707 , k 3 = 7.7541 ,   k 4 = 4.5835 , g 1 = 0.4472 , and   g 2 = 1.0461 , Figure 20 shows the change in position coordinates of the complex system and the expected formation state of the complex system. It can be observed that the complex system can quickly complete the triangular formation and that the error between the actual position and the set value is small. In addition, the cooperative optimal formation protocol speeds up the convergence rate of the complex system, which is of great help to the formation time of the complex system.

5. Conclusions

In this paper, a heterogeneous multi-agent system has been established by analyzing the dynamics model of the unmanned ground vehicle and the unmanned aerial vehicle. Firstly, the formation control protocol is proposed based on the communication topology of a multi-agent system. Then, according to the internal state of a single agent, the optimal control law of a single agent system is designed using the optimal control theory, and the optimal control law is introduced into the system to achieve the distributed optimal formation. Finally, based on the cooperative architecture of the heterogeneous multi-agent system, the cooperative formation design of the heterogeneous multi-agent system is carried out, and the optimal control theory is introduced into the heterogeneous multi-agent system to realize the optimal cooperative formation of the heterogeneous system. The stability of the system is further analyzed by graph theory. The communication topology of the multi-agent system does not interfere with the protocol and the protocol can optimize the performance function while the system completes the formation task. The simulation results show that the optimal control can accelerate the convergence speed of the system and greatly help the system to quickly reach the desired formation state. In the next step, we plan to investigate the anomaly detection and recognition problems under heterogeneous multi-agent cooperation architecture, and we plan to apply the theoretical research results in practice to engineering applications.

Author Contributions

Conceptualization, Y.L. and M.L.; methodology, Y.L.; software, M.L. and J.L.; validation, Y.L., M.L. and J.L.; writing original draft preparation, Y.L.; resources, Y.G.; writing review and editing, M.L.; supervision, Y.G.; funding acquisition, Y.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (61872204), the Scientific Research Project of Heilongjiang Provincial Universities, China (Grant No.145109143).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author, Meichen Liu, upon reasonable request.

Acknowledgments

We appreciate all the authors for their contributions and the support of the foundation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cheung, C.; Grocholsky, B. UAV-UGV Collaboration with a PackBot UGV and Raven SUAV for Pursuit and Tracking of a Dynamic Target. In Unmanned Systems Technology X; SPIE: Bellingham, WA, USA, 2008; Volume 6962, pp. 381–390. [Google Scholar]
  2. Moseley, M.B.; Grocholsky, B.P.; Cheung, C.; Singh, S. Integrated long-range UAV/UGV collaborative target tracking. In Unmanned Systems Technology XI; SPIE: Bellingham, WA, USA, 2009; Volume 7332, pp. 30–40. [Google Scholar]
  3. Schulteis, T.M.; Price, J.G. Project stork UAV/UGV collaborative initiative. In Unmanned Ground Vehicle Technology VI; SPIE: Bellingham, WA, USA, 2004; Volume 5422, pp. 414–425. [Google Scholar]
  4. Cao, Y.; Yu, W.; Ren, W.; Chen, G. An overview of recent progress in the study of distributed multi-agent coordination. IEEE Trans. Ind. Inform. 2013, 9, 427–438. [Google Scholar] [CrossRef] [Green Version]
  5. Chen, F.; Dimarogonas, D.V. Leader–follower formation control with prescribed performance guarantees. IEEE Trans. Control. Netw. Syst. 2021, 8, 450–461. [Google Scholar] [CrossRef]
  6. Han, N.; Luo, X. Tracking and distributed formation control for leader-following heterogeneous multi-agent systems. In Proceedings of the 2016 35th Chinese Control Conference (CCC), Chengdu, China, 27–29 July 2016; pp. 7897–7901. [Google Scholar]
  7. Li, S.; Zhang, J.; Li, X.; Wang, F.; Luo, X.; Guan, X. Formation control of heterogeneous discrete-time nonlinear multi-agent systems with uncertainties. IEEE Trans. Ind. Electron. 2017, 64, 4730–4740. [Google Scholar] [CrossRef]
  8. Liang, S.; Wang, F.; Chen, Z.; Liu, Z. Formation Control for Discrete-Time Heterogeneous Multi-Agent Systems; John Wiley & Sons: Hoboken, NJ, USA, 2022. [Google Scholar]
  9. Brandão, A.S.; Sarcinelli-Filho, M.; Carelli, R. Leader-following control of a UAV-UGV formation. In Proceedings of the 2013 16th International Conference on Advanced Robotics (ICAR), Montevideo, Uruguay, 25–29 November 2013; pp. 1–6. [Google Scholar]
  10. Luo, Q.; Duan, H. An improved artificial physics approach to multiple UAVs/UGVs heterogeneous coordination. Sci. China Technol. Sci. 2013, 56, 2473–2479. [Google Scholar] [CrossRef]
  11. Rahimi, R.; Abdollahi, F.; Naqshi, K. Time-varying formation control of a collaborative heterogeneous multi agent system. Robot. Auton. Syst. 2014, 62, 1799–1805. [Google Scholar] [CrossRef]
  12. Abbaspour, A.; Alipour, K.; Jafari, H.Z.; Moosavian, S.A.A. Optimal formation and control of cooperative wheeled mobile robots. Comptes Rendus Mécanique 2015, 343, 307–321. [Google Scholar] [CrossRef]
  13. Ai, X.L.; Yu, J.Q.; Chen, Y.B.; Chen, F.Z.; Shen, Y.C. Optimal formation control with limited communication for multi-unmanned aerial vehicle in an obstacle-laden environment. Part G J. Aerosp. Eng. 2017, 231, 979–997. [Google Scholar] [CrossRef]
  14. Lin, W.; Zhao, W.; Liu, H. Robust optimal formation control of heterogeneous multi-agent system via reinforcement learning. IEEE Access 2020, 8, 218424–218432. [Google Scholar] [CrossRef]
  15. Huang, R.; Zhou, J.; Huang, H. Flexible formation control through time-varying weighting multi-agent flocking. In Proceedings of the 2019 Chinese Control And Decision Conference (CCDC), Nanchang, China, 3–5 June 2019; pp. 3424–3429. [Google Scholar]
  16. Yan, Z.-P.; Liu, Y.-B.; Yu, C.-B.; Zhou, J.-J. Leader-following coordination of multiple UUVs formation under two independent topologies and time-varying delays. J. Cent. South Univ. 2017, 24, 382–393. [Google Scholar] [CrossRef]
  17. Zhao, J.; Dai, F.; Song, Y. A Distributed Optimal Formation Control for Multi-Agent System of UAVs. In Proceedings of International Conference on Artificial Life and Robotics; ALife Robotics Corporation Ltd.: Oita, Japan, 2022. [Google Scholar]
  18. Zhou, J.; Zeng, D.; Lu, X. Multi-agent trajectory-tracking flexible formation via generalized flocking and leader-average sliding mode control. IEEE Access 2020, 8, 36089–36099. [Google Scholar] [CrossRef]
  19. Cui, Y.; Fei, M.; Du, D. Event-triggered cooperative compensation control for consensus of heterogeneous multi-agent systems. IET Control. Theory Appl. 2016, 10, 1573–1582. [Google Scholar] [CrossRef]
  20. Li, S.; Feng, G.; Wang, J.; Luo, X.; Guan, X. Adaptive control for cooperative linear output regulation of heterogeneous multi-agent systems with periodic switching topology. IET Control. Theory Appl. 2015, 9, 34–41. [Google Scholar] [CrossRef]
  21. Yuan, C.; He, H. Cooperative output regulation of heterogeneous multi-agent systems with a leader of bounded inputs. IET Control. Theory Appl. 2018, 12, 233–242. [Google Scholar] [CrossRef]
  22. Ren, W.; Beard, R.W. Distributed Consensus in Multi-Vehicle Cooperative Control; Springer: London, UK, 2008. [Google Scholar]
Figure 1. Formation control protocol for heterogeneous multi-agent systems.
Figure 1. Formation control protocol for heterogeneous multi-agent systems.
Entropy 24 01440 g001
Figure 2. Distributed optimal formation control protocol for heterogeneous multi-agent systems.
Figure 2. Distributed optimal formation control protocol for heterogeneous multi-agent systems.
Entropy 24 01440 g002
Figure 3. Position and velocity states change with time using the formation control protocol.
Figure 3. Position and velocity states change with time using the formation control protocol.
Entropy 24 01440 g003
Figure 4. Attitude angle changes with time using a formation control protocol.
Figure 4. Attitude angle changes with time using a formation control protocol.
Entropy 24 01440 g004
Figure 5. Changes in position and velocity states with time using the distributed optimal formation control protocol.
Figure 5. Changes in position and velocity states with time using the distributed optimal formation control protocol.
Entropy 24 01440 g005
Figure 6. Attitude angle changes with time using the distributed optimal formation control protocol.
Figure 6. Attitude angle changes with time using the distributed optimal formation control protocol.
Entropy 24 01440 g006
Figure 7. Distributed cooperative formation control for heterogeneous multi-agent systems.
Figure 7. Distributed cooperative formation control for heterogeneous multi-agent systems.
Entropy 24 01440 g007
Figure 8. Distributed cooperative optimal formation control for heterogeneous multi-agent systems.
Figure 8. Distributed cooperative optimal formation control for heterogeneous multi-agent systems.
Entropy 24 01440 g008
Figure 9. UAV system state formation based on formation control.
Figure 9. UAV system state formation based on formation control.
Entropy 24 01440 g009
Figure 10. UAV system state formation based on distributed optimal formation control.
Figure 10. UAV system state formation based on distributed optimal formation control.
Entropy 24 01440 g010
Figure 11. UGV system state formation based on formation control.
Figure 11. UGV system state formation based on formation control.
Entropy 24 01440 g011
Figure 12. UGV system state formation based on distributed optimal formation control.
Figure 12. UGV system state formation based on distributed optimal formation control.
Entropy 24 01440 g012
Figure 13. UAV system state formation based on heterogeneous cooperative formation control.
Figure 13. UAV system state formation based on heterogeneous cooperative formation control.
Entropy 24 01440 g013
Figure 14. UAV system state formation based on heterogeneous cooperative optimal formation control.
Figure 14. UAV system state formation based on heterogeneous cooperative optimal formation control.
Entropy 24 01440 g014
Figure 15. UGV system state formation based on heterogeneous cooperative formation control.
Figure 15. UGV system state formation based on heterogeneous cooperative formation control.
Entropy 24 01440 g015
Figure 16. UGV system state formation based on heterogeneous cooperative optimal formation control.
Figure 16. UGV system state formation based on heterogeneous cooperative optimal formation control.
Entropy 24 01440 g016
Figure 17. Formation control protocol for complex heterogeneous multi-agent systems.
Figure 17. Formation control protocol for complex heterogeneous multi-agent systems.
Entropy 24 01440 g017
Figure 18. Distributed optimal formation control protocol for complex heterogeneous multi-agent systems.
Figure 18. Distributed optimal formation control protocol for complex heterogeneous multi-agent systems.
Entropy 24 01440 g018
Figure 19. Distributed cooperative formation control for complex heterogeneous multi-agent systems.
Figure 19. Distributed cooperative formation control for complex heterogeneous multi-agent systems.
Entropy 24 01440 g019
Figure 20. Distributed cooperative optimal formation control for complex heterogeneous multi-agent systems.
Figure 20. Distributed cooperative optimal formation control for complex heterogeneous multi-agent systems.
Entropy 24 01440 g020
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Y.; Liu, M.; Lian, J.; Guo, Y. Collaborative Optimal Formation Control for Heterogeneous Multi-Agent Systems. Entropy 2022, 24, 1440. https://doi.org/10.3390/e24101440

AMA Style

Li Y, Liu M, Lian J, Guo Y. Collaborative Optimal Formation Control for Heterogeneous Multi-Agent Systems. Entropy. 2022; 24(10):1440. https://doi.org/10.3390/e24101440

Chicago/Turabian Style

Li, Yandong, Meichen Liu, Jiya Lian, and Yuan Guo. 2022. "Collaborative Optimal Formation Control for Heterogeneous Multi-Agent Systems" Entropy 24, no. 10: 1440. https://doi.org/10.3390/e24101440

APA Style

Li, Y., Liu, M., Lian, J., & Guo, Y. (2022). Collaborative Optimal Formation Control for Heterogeneous Multi-Agent Systems. Entropy, 24(10), 1440. https://doi.org/10.3390/e24101440

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