Next Article in Journal
A Multicriteria Simheuristic Approach for Solving a Stochastic Permutation Flow Shop Scheduling Problem
Next Article in Special Issue
Modeling and Control of IPMC-Based Artificial Eukaryotic Flagellum Swimming Robot: Distributed Actuation
Previous Article in Journal
PM2.5 Concentration Prediction Based on CNN-BiLSTM and Attention Mechanism
Previous Article in Special Issue
Optimal Coronavirus Optimization Algorithm Based PID Controller for High Performance Brushless DC Motor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Containment Control of First-Order Multi-Agent Systems under PI Coordination Protocol †

1
Key Laboratory of Advanced Process Control for Light Industry (Ministry of Education), Institute of Automation, Jiangnan University, Wuxi 214122, China
2
School of Automation, Nanjing University of Science and Technology, Nanjing 210094, China
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in the 2020 Chinese Automation Congress (CAC).
Algorithms 2021, 14(7), 209; https://doi.org/10.3390/a14070209
Submission received: 15 June 2021 / Revised: 8 July 2021 / Accepted: 12 July 2021 / Published: 14 July 2021
(This article belongs to the Special Issue Algorithms for PID Controller 2021)

Abstract

:
This paper investigates the containment control problem of discrete-time first-order multi-agent system composed of multiple leaders and followers, and we propose a proportional-integral (PI) coordination control protocol. Assume that each follower has a directed path to one leader, and we consider several cases according to different topologies composed of the followers. Under the general directed topology that has a spanning tree, the frequency-domain analysis method is used to obtain the sufficient convergence condition for the followers achieving the containment-rendezvous that all the followers reach an agreement value in the convex hull formed by the leaders. Specially, a less conservative sufficient condition is obtained for the followers under symmetric and connected topology. Furthermore, it is proved that our proposed protocol drives the followers with unconnected topology to converge to the convex hull of the leaders. Numerical examples show the correctness of the theoretical results.

1. Introduction

Distributed coordination control of multi-agent systems has been an active research area for its widespread potential engineering applications including cooperative surveillance, sensor networks, spacecraft formation flying, etc. The most fundamental problem studied in coordination control is consensus problem [1,2,3,4,5,6,7,8,9,10,11] which means that each agent reaches an agreement based on the information from its relative agents. According to the different number of leaders in the multi-agent system, the consensus problem is generally divided into leaderless case [12,13], single-leader–follower case [14,15,16], and multiple-leader–follower case.
The containment control problem [17,18,19,20,21,22,23,24] is considered as a special multiple-leader–follower consensus problem for multi-agent system with multiple leaders and followers, and it requires all followers to converge into the convex hull spanned by leaders. Many valuable results on the containment control problem of first-order multi-agent systems, which are investigated in this paper, have been obtained in recent decades. Basic results for realizing containment of continuous-time first-order multi-agent systems with stationary leaders have been given by Liu in [25], and the convergence conditions under fixed and switching topologies are dependent on the topology structure. Wang and their colleagues [26] investigated the containment problem of first-order multi-agent system with communication noises, and designed a time-varying gain to reduce noises. A PD-type control protocol was introduced and a parameter condition was given to guarantee the containment achievement under input delays in Rong’s work [27]. Mu and partners [28] provided necessary and sufficient criteria for containment convergence if the communication data rates are limited. Miao and colleagues [29] applied the event-triggered scheme into containment control algorithms, and the convergence problem with single time delay and multiple time delays were discussed, respectively. Considering the discrete-time first-order multi-agent systems. The authors in [30] adopted the asynchronous containment control protocol and got sufficient conditions for reaching containment. Additionally, containment problem of multi-agent systems with second-order dynamics and general linear dynamics has also attracted extensive attention [31,32,33,34].
According to aforementioned works, we note that general containment control protocols make each follower finally converge to different point in the convex hull spanned by leaders. General containment problem only focuses on the containment and ignores the consensus, which is called rendezvous in this paper, of the followers. If we take both containment and consensus into consideration, a different control protocol is required, and we define this problem as a containment-rendezvous problem in this paper. Actually, the average-tracking problem of reference signals is a typical containment-rendezvous problem, if we regard the reference signals as leaders. Shan and Liu [35] considered the average-tracking problem of first-order multi-agent systems with unmatched reference signals, and used the frequency-domain analysis method, which will be adopted in this paper, to get the convergence condition. On the basis of average-consensus algorithm, Chung [36] designed a containment-rendezvous algorithm that made all followers track one point in the convex hull of first-order agents, but the results depended on the balance property of communication topology.
Inspired by above works, this paper focuses on the containment-rendezvous problem of discrete-time first-order multi-agent systems under general communication topology. On the basis of the average-consensus algorithm, a proportional-integral (PI) coordination control protocol is proposed for reaching containment-rendezvous. We first analyze the protocol for the followers under general connected topology, and obtain a sufficient convergence condition. Besides, we consider the followers with symmetric and connected topology and get a less conservative convergence condition which has been simply discussed in our initial work [37]. Furthermore, we investigate the followers under unconnected topology that is regarded as a union of several connected parts, and obtain the convergence condition according to that of general connected topology.
This paper is organized as follows. We give some basic concepts about graph theory, the multi-agent systems, and the coordination control protocol in Section 2. The convergence of the protocol under general connected topology, symmetric and connected topology and unconnected topology of the followers is proved in Section 3. In Section 4 and Section 5, the numerical simulations and the conclusion are presented, respectively.

2. Problem Formulation

2.1. Graph Theory

Consider a multi-agent system with n agents is denoted by a graph G ( V , E ) , where V = { 1 , 2 , , n } and E V × V stand for the vertex set and edge set, respectively. An edge ( i , j ) E represents that agent j is able to access information of agent i and means that vertex i is a neighbor of vertex j. If agent j has no neighbor, it is called a leader, otherwise, it is a follower. The index set is denoted as N j = { i V : ( i , j ) E , i j } . A directed path from i to j is a sequence of edges in a graph of the form ( i , h 0 ) , ( h 1 , h 2 ) , , ( h k , j ) , where h k V . The adjacency matrix is a nonnegative matrix A = [ a i j ] R n × n defined as a j i > 0 if ( i , j ) E , and a j i = 0 otherwise. Furthermore, self edges are not allowed in this paper, i.e., a i i = 0 . The Laplacian matrix is defined as L = [ l i j ] R n × n , where l i i = j = 1 , j i n a i j and l i j = a i j , i j .
A directed graph is called a directed tree if each node in graph has exactly one parent except for one node which is called the root, and the root has directed paths to each other node. A directed spanning tree of a directed graph is a direct tree that contains all nodes of the directed graph. A directed graph has a spanning tree if there exists a directed spanning tree as a subset of the directed graph.

2.2. Agents’ Dynamics and Coordination Protocol

Investigate a discrete-time multi-agent system consisting of m leaders and n m followers labelled by 1 , , m and m + 1 , , n , respectively. The dynamic model of agent i is given by
x i ( k + 1 ) = x i ( k ) , i = 1 , , m , x i ( k + 1 ) = x i ( k ) + u i ( k ) , i = m + 1 , , n ,
where x i ( k ) R p and u i ( k ) R p are the state and control input of agent i, respectively.
According to the PI control strategy, we use the following PI coordination control algorithm for the first-order agents,
u i ( k ) = γ 1 j = 1 n a i j ( x j ( k ) x i ( k ) ) + γ 2 r i ( k ) , r i ( k + 1 ) = r i ( k ) + j = m + 1 n a i j ( x j ( k ) x i ( k ) ) , i = m + 1 , , n ,
where a i j is the ( i , j ) entry of the adjacency matrix A, r i ( k ) represents the integral term, γ 1 and γ 2 are positive gain parameters to be decided for the proportional term and integral term, respectively.
The states of leaders (1) remain static since the inputs of leaders are always zero, so we only investigate the followers’ dynamics here. With the protocol (2), the dynamics of follower i are rewritten as
x i ( k + 1 ) = x i ( k ) + γ 1 j = 1 n a i j ( x j ( k ) x i ( k ) ) + γ 2 r i ( k ) , r i ( k + 1 ) = r i ( k ) + j = m + 1 n a i j ( x j ( k ) x i ( k ) ) , i = m + 1 , , n .
In order to analyze the convergence performance of system (3), we introduce two topologies, one of which is named as leader–follower topology composed of the leaders and followers, and the other one is named as follower topology composed of the followers. The Laplacian matrix of leader–follower topology is L given by
L = 0 m × m 0 m × ( n m ) L 1 L 2 ,
where L 1 represents the topology between leaders and followers, and L 2 represents the topology among followers. Meanwhile, the Laplacian matrix corresponding to the follower topology is L F formulated as
L F = L 2 D ,
where the diagonal matrix D is defined as
D = d i a g { j = 1 m a m + 1 , j , , j = 1 m a n , j } .
Generally, we need the basic assumption on the leader–follower topology of system (3) as follows.
Assumption 1.
For each of the followers in the leader–follower topology, there is at least one leader that has a directed path to the follower.
On the basis of Assumption 1, we make further assumptions on the follower topology as follows.
Assumption 2.
The follower topology has a directed spanning tree.
Assumption 3.
The edges of the follower topology are bidirectional, i.e., Laplacian matrix L F is symmetric, and it has a directed spanning tree.
Then we have the following lemma.
Lemma 1
([38]). The Laplacian matrix L F of the follower topology has a simple eigenvalue 0 with the corresponding eigenvector [ 1 , 1 , , 1 ] T and all the other eigenvalues have positive real parts if and only if the topology satisfies Assumption 2.

3. Main Results

System (3) is reformulated in a vector form as
X F ( k + 1 ) = ( ( I n m γ 1 L 2 ) I p ) X F ( k ) γ 1 ( L 1 I p ) X L ( k ) + γ 2 R ( k ) , R ( k + 1 ) = R ( k ) ( L F I p ) X F ( k ) ,
where X L ( k ) = [ x 1 T ( k ) , , x m T ( k ) ] T , X F ( k ) = [ x m + 1 T ( k ) , , x n T ( k ) ] T , R ( k ) = [ r m + 1 T ( k ) , , r n T ( k ) ] T .
Define R ^ ( k ) = R ( k ) γ 1 γ 2 ( L 1 I p ) X L ( k ) , and we get
X F ( k + 1 ) = ( I n m γ 1 L 2 I p ) X F ( k ) + γ 2 R ^ ( k ) , R ^ ( k + 1 ) = R ^ ( k ) ( L F I p ) X F ( k ) .
Let Y ( k ) = [ X F T ( k ) , R ^ T ( k ) ] T , and the system (5) is expressed in a compact form as
Y ( k + 1 ) = ( I n m γ 1 L 1 γ 2 I n m L F I n m I p ) Y ( k ) .
To continue the convergence analysis of system (6), some useful lemmas are listed firstly.
Lemma 2
([39]). Let P ( z ) be a polynomial of order two with complex coefficients in the form of P ( z ) = z 2 + ( p 1 + j q 1 ) z + p 2 + j q 2 , where j is the imaginary unit. The polynomial P ( z ) has all its zeros in the open left half of the z-complex plane if and only if p 1 > 0 and p 1 2 p 2 + p 1 q 1 q 2 q 2 2 > 0 .
Lemma 3
([35]). If Q 4 is invertible, det Q 1 Q 2 Q 3 Q 4 = det ( Q 4 ) det ( Q 1 Q 2 Q 4 1 Q 3 ) .
Lemma 4
([40]). Let Q C n × n , Q = Q * 0 and T = d i a g { t i , t i C } , then
λ ( Q T ) ρ ( Q ) C o ( 0 { t i } ) ,
where λ ( · ) denotes matrix eigenvalue, ρ ( · ) denotes matrix spectral radius and C o ( · ) denotes the convex hull.
Next, we will obtain the convergence conditions of the system (6) according to different follower topologies.

3.1. General Directed Follower Topology

Theorem 1.
Consider the multi-agent system (3) with leader–follower and follower topologies satisfying Assumptions 1 and 2, respectively. With condition i = m + 1 n q i r i ( 0 ) = 0 , all the followers reach containment-rendezvous asymptotically that the followers converge to an agreement value in the convex hull spanned by the leaders, if γ 1 and γ 2 satisfy
γ 1 γ 2 > 0 ,
( γ 1 γ 2 ) 2 ( γ 2 2 γ 1 + 4 R e ( κ i ) | κ i | 2 ) 4 γ 2 I m ( κ i ) 2 | κ i | 4 > 0 ,
for i = m + 1 , , n , and
| λ ( γ 2 D Θ 1 ( e j ω ) ) | < 1
hold with ω ( 0 , π ] , where Θ ( e j ω ) = ( e j ω 1 ) 2 I + γ 1 ( e j ω 1 ) ( L F + D ) + γ 2 ( L F + D ) and κ i , i = m + 1 , , n represent the eigenvalues of the matrix L F + D .
Proof. 
According to the properties of the Kronecker product, we set the agents’ state dimension as p = 1 in the following proof, and system (6) is written as
Y ^ ( k + 1 ) = I n m γ 1 ( L F + D ) γ 2 I n m L F I n m Y ^ ( k ) .
Meanwhile, we divide the proof into two steps including convergence analysis and the analysis of final rendezvous state.
Step 1: To analyze the convergence performance of system (8), we investigate the characteristic equation of (8) as follows,
det ( z I I γ 1 ( L F + D ) γ 2 I L F I ) = 0 ,
and it can be reformulated from Lemma 3 as
det ( ( z 1 ) 2 I + γ 1 ( z 1 ) ( L F + D ) + γ 2 L F ) = 0 .
Evidently, it is obtained from Lemma 1 that the Equation (10) has a root at z = 1 .
Before analyzing Equation (10), we first pay attention to the following equation
det ( ( z 1 ) 2 I + γ 1 ( z 1 ) ( L F + D ) + γ 2 ( L F + D ) ) = 0 , .
Obviously, (11) is equivalent to
z 2 + ( γ 1 κ i 2 ) z + ( γ 2 γ 1 ) κ i + 1 = 0 , i = m + 1 , , n .
Instead of studying Equation (12) directly, we apply the bilinear transformation s = z + 1 z 1 to it and get
γ 2 κ i s 2 + 2 ( γ 1 γ 2 ) κ i s + ( γ 2 2 γ 1 ) κ i + 4 = 0 .
Thus, Equation (12) has all roots within the unit circle, if and only if Equation (13) has all roots in the open left half complex plane. Then, we reformulate Equation (13) as
s 2 + ( p 1 + j q 1 ) s + p 2 i + j q 2 i = 0 , i = m + 1 , , n ,
where p 1 = 2 ( γ 1 γ 2 ) γ 2 , q 1 = 0 , p 2 i = γ 2 2 γ 1 γ 2 + 4 R e ( κ i ) γ 2 | κ i | 2 , q 2 i = 4 I m ( κ i ) γ 2 | κ i | 2 .
Based on Lemma 2, the polynomial has all zeros in the open left half complex plane if and only if the gains γ 1 and γ 2 satisfy
p 1 > 0 , p 1 2 p 2 i + p 1 q 1 q 2 i q 2 i 2 > 0 , i = m + 1 , , n .
which is equivalent to
γ 1 γ 2 > 0 , ( γ 1 γ 2 ) 2 ( γ 2 2 γ 1 + 4 R e ( κ i ) | κ i | 2 ) 4 γ 2 I m 2 ( κ i ) | κ i | 4 > 0 , i = m + 1 , n .
Thus, the roots of Equation (11) all lie inside the unit circle if and only if there exist gain parameters γ 1 and γ 2 satisfying condition (16).
Back to Equation (10), we reformulate it as
det ( I γ 2 D Θ 1 ( z ) ) = 0 ,
where Θ ( z ) = ( z 1 ) 2 I + γ 1 ( z 1 ) ( L F + D ) + γ 2 ( L F + D ) . Under condition (16), it is obvious that if condition | λ ( γ 2 D Θ 1 ( e j ω ) ) | < 1 holds with ω ( 0 , π ] , the roots of Equation (17) lie inside the unit circle except for one root at z = 1 . Thus, the characteristic Equation (9) has all roots within or on the unit circle and we have finally proved the asymptotic convergence of the system under Assumption 1, i.e.
lim k Y ^ ( k ) = lim k [ X F T ( k ) , R ^ T ( k ) ] T = [ X d T , R ^ d T ] T ,
where X d R n m and R ^ d R n m are constant vectors.
Step 2: We will prove that all followers reach an agreement value in convex hull spanned by leaders.
According to Equation (4), we obtain with p = 1
lim k L F X F ( k ) = 0 .
Under Assumption 1,we get from Lemma 1 and Equation (19)
lim k x i ( k ) = x d , i = m + 1 , , n ,
where x d R is a constant. From Equation (20), it is clear that followers converge to the same value. Taking the z transformation of Equation (4) with p = 1 , we get
z X F ( z ) z X F ( 0 ) = ( I γ 1 L 2 ) X F ( z ) + γ 2 R ( z ) z z 1 γ 1 L 1 X L . z R ( z ) z R ( 0 ) = R ( z ) L F X F ( z )
Re-express Equation (21) as
X F ( z ) = [ ( z 1 ) I + γ 1 L 2 + γ 2 L F z 1 ] 1 [ z X F ( 0 ) + z z 1 γ 2 R ( 0 ) z z 1 γ 1 L 1 X L ] = [ ( z 1 ) I + γ 1 D + ( γ 1 + γ 2 z 1 ) L F ] 1 [ z X F ( 0 ) + z z 1 γ 2 R ( 0 ) z z 1 γ 1 L 1 X L ] .
Since the convergence of the system has been proved, using the final value theorem yields
lim z 1 ( z 1 ) X F ( z ) = [ ( z 1 ) I + γ 1 D + ( γ 1 + γ 2 z 1 ) L F ] 1 [ z ( z 1 ) X F ( 0 ) + z γ 2 R ( 0 ) z γ 1 L 1 X L ] = x d [ 1 , 1 , , 1 ] T .
Letting W ( z ) = ( z 1 ) X F ( z ) , we have lim z 1 W ( z ) = x d [ 1 , 1 , , 1 ] T and
[ ( z 1 ) I + γ 1 D + ( γ 1 + γ 2 z 1 ) L F ] W ( z ) = [ z ( z 1 ) X F ( 0 ) + γ 2 z R ( 0 ) z γ 1 L 1 X L ] .
Multiplying Q = [ q m + 1 , q m + 2 , , q n ] on both sides of (24), we get
Q [ ( z 1 ) I + γ 1 D ] W ( z ) = Q [ z ( z 1 ) X F ( 0 ) + γ 2 z R ( 0 ) z γ 1 L 1 X L ] ,
where Q is the left eigenvector of L F corresponding to eigenvalue 0. Then, we take the limit of Equation (25) as z 1 and have
Q γ 1 D x d [ 1 , 1 , , 1 ] T = Q [ γ 2 R ( 0 ) γ 1 L 1 X L ] .
Equation (26) is rewritten as
γ 1 x d i = m + 1 n j = 1 m q i a i j = γ 1 i = m + 1 n j = 1 m q i a i j x j + γ 2 i = m + 1 n q i r i ( 0 ) ,
and we finally get
x d = i = m + 1 n j = 1 m q i a i j x j i = m + 1 n j = 1 m q i a i j ,
with i = m + 1 n q i r i ( 0 ) = 0 . Apparently, x d is in the convex hull formed by leaders.
When p 1 , the proof is the same while the final state x d is a constant vector instead of a constant. Hence, we can draw the conclusion that all followers will eventually reach an agreement value in the convex hull spanned by leaders under our proposed protocol. Theorem 1 is proved. □

3.2. Symmetric Follower Topology

By means of some existing results [40,41], we present the convergence condition, which is less conservative than (7), of system (3) under symmetric follower topology.
Theorem 2.
The leader–follower and follower topologies of multi-agent system (3) satisfy Assumptions 1 and 3, respectively. With i = m + 1 n r i ( 0 ) = 0 , all the followers converge to an agreement value i = m + 1 n j = 1 m a i j x j i = m + 1 n j = 1 m a i j that lies in convex hull spanned by the leaders, if γ 1 and γ 2 satisfy
γ 2 γ 1 2 θ i ,
and
ρ ( L F ) ( γ 2 2 γ 1 ) 2 γ 1 θ i + 4 > 0 ,
for i = m + 1 , , n , where θ i = j = 1 m a i j
Proof. 
We have already known that the characteristic Equation (9) has a root as z = 1 , so we only consider the situation that z 1 . From Lemma 3, Equation (9) equals to
det ( ( z 1 ) ( ( z 1 ) I + γ 1 D ) + ( γ 1 ( z 1 ) + γ 2 ) L F ) = 0 .
When z 1 , (30) is reformulated as
det ( I + G ( z ) L F ) = 0 ,
where
G ( z ) = γ 1 ( z 1 ) + γ 2 ( z 1 ) ( z 1 + γ 1 θ m + 1 ) γ 1 ( z 1 ) + γ 2 ( z 1 ) ( z 1 + γ 1 θ n ) ,
and θ i = j = 1 m a i j , i = m + 1 , , n is the ( i , i ) entry of the diagonal matrix D.
Let F ( z ) = det ( I + G ( z ) L F ) . According to the generalized Nyquist stability criterion [41], the zeros of F ( z ) are in the unit circle when G ( z ) does not have poles out of the unit circle, if λ ( G ( e j ω ) L F ) does not enclose the point ( 1 , j 0 ) for ω [ π , π ] . Because of the symmetry of the Laplacian matrix L F , it follows from Lemma 4 that
λ ( G ( e j ω ) L F ) ρ ( L F ) C o ( 0 g i ( e j ω ) ) , i = m + 1 , , n ,
where
g i ( e j ω ) = γ 1 ( e j ω 1 ) + γ 2 ( e j ω 1 ) ( e j ω 1 + γ 1 θ i ) .
For calculating conveniently, we assume that γ 2 γ 1 . In order to ensure that all the poles of G ( e j ω ) are within or on the unit circle, we have
0 γ 1 θ i 2 .
Equation (33) is rewritten as
g i ( e j ω ) = γ 1 cos ω + γ 2 γ 1 + j γ 1 sin ω a + j b ,
where a = 2 cos 2 ω + ( γ 1 θ i 2 ) cos ω γ 1 θ i and b = sin ω ( 2 cos ω + γ 1 θ i 2 ) .
To analyze the intersections of g i ( e j ω ) on the real axis, we get
I m ( g i ( e j ω ) ) = ( 2 ( γ 1 γ 2 ) cos ω + ( γ 1 γ 2 ) ( γ 1 θ i 2 ) γ 1 2 θ i ) sin ω = 0 .
For ω ( 0 , π ] , Equation (36) has only one solution ω = π and
g i ( e j π ) = γ 2 2 γ 1 4 2 γ 1 θ i .
It is evident that ρ ( L F ) C o ( 0 g i ( e j ω ) ) does not enclose the point ( 1 , j 0 ) , if all the points ( ρ ( L F ) g i ( e j π ) , j 0 ) are on the right side of the point ( 1 , j 0 ) for i = m + 1 , , n , i.e.,
ρ ( L F ) γ 2 2 γ 1 4 2 γ 1 θ i > 1 .
Hence, λ ( G ( e j ω ) L F ) does not enclose the point ( 1 , j 0 ) and G ( z ) has no poles out of the unit circle, if we choose the gain parameters γ 1 and γ 2 satisfying γ 2 γ 1 2 θ i , and ρ ( L F ) ( γ 2 2 γ 1 ) 2 γ 1 θ i + 4 > 0 , for i = m + 1 , , n . Thus, F ( z ) has all zeros within the unit circle, which means that the roots of characteristic Equation (9) are within or on the unit circle. Hence, the convergence of the system is proved.
The analysis of final rendezvous state is omitted here for it is almost same as the proof of Theorem 1. The only difference is that the left eigenvector of the symmetric Laplacian Matrix L F corresponding to eigenvalue 0 becomes [ 1 , , 1 ] , and the final value is
x d = i = m + 1 n j = 1 m a i j x j i = m + 1 n j = 1 m a i j .
Theorem 2 is proved. □

3.3. Unconnected Follower Topology

Considering the follower topology is unconnected, i.e., it has no spanning tree, the containment-rendezvous cannot be achieved by our proposed PI coordination control protocol. In this case, we can divide the follower topology into several connected parts, and obtain the convergence conditions based on the results in Theorem 1.
In order to analyze the convergence behavior, we divide the unconnected follower topology into N connected parts, each of which has a spanning tree or has only one agent, and we labelled the parts as 1 , , N . For each part, the state and integral vectors of followers are defined as X F l ( k ) and R l ( k ) , l = 1 , 2 , , N . Hence, the dynamic models of followers are formulated as
X F l ( k + 1 ) = ( I n m γ 1 L 2 l I p ) X F l ( k ) γ 1 ( L 1 l I p ) X L ( k ) + γ 2 R l ( k ) , R l ( k + 1 ) = R l ( k ) ( L F l I p ) X F l ( k ) , l = 1 , , N .
where L 1 l , L 2 l and L F l are same as the definitions of above L 1 , L 2 and L F , respectively.
For each part l, the characteristic equation is given by
det ( ( z 1 ) 2 I + γ 1 ( z 1 ) ( L F l + D l ) + γ 2 L F l ) = 0 ,
where D l is same as the definition of above D.
Similar to Theorem 1, we take into account the following equation
det ( ( z 1 ) 2 I + γ 1 ( z 1 ) ( L F l + D l ) + γ 2 ( L F l + D l ) ) = 0 .
Then, in the light of Theorem 1, we come to the following results.
Theorem 3.
The leader–follower topology of multi-agent system (3) satisfies Assumption 1 and the follower topology is unconnected. All the n m followers converge to the convex hull spanned by the m leaders, if for all N parts, the roots of Equation (41) lie in the unit circle, and
| λ ( γ 2 D l Θ l 1 ( e j ω ) ) | < 1
hold with ω ( 0 , π ] , where Θ ( e j ω ) = ( e j ω 1 ) 2 I + γ 1 ( e j ω 1 ) ( L F l + D l ) + γ 2 ( L F l + D l ) .
Proof. 
Divide the unconnected follower topology into N connected parts (40), and the state and the integral term of followers are expressed as X F ( k ) = [ X F 1 ( k ) , , X F N ( k ) ] T and R ( k ) = [ R 1 ( k ) , , R N ( k ) ] T . It is evident the the dynamics (40) of each part have completely same form as (4).
According to the proof of Theorem 1, the followers in one part reach the containment-rendezvous under condition in Theorem 3. In each part. the followers reach an agreement value in convex hull spanned by the leaders. Evidently, the convex hull spanned by the leaders in each part must be contained in the convex hull composed of all the leaders in the system. Hence, all followers converge to the convex hull composed of all the leaders. □

4. Simulations

In this section, we consider a first-order multi-agent system including 6 leaders and 5 followers labelled by 1 , , 6 and 7 , , 11 , respectively. The initial states of agent i are set to x 1 ( 0 ) = [ 4 , 4 ] T , x 2 ( 0 ) = [ 5 , 2 ] T , x 3 ( 0 ) = [ 6 , 7 ] T , x 4 ( 0 ) = [ 7 , 1 ] T , x 5 ( 0 ) = [ 9 , 3 ] T , x 6 ( 0 ) = [ 10 , 8 ] T , x 7 ( 0 ) = [ 5 , 1 ] T , x 8 ( 0 ) = [ 9 , 2 ] T , x 9 ( 0 ) = [ 2 , 6 ] T , x 10 ( 0 ) = [ 8 , 10 ] T and x 11 ( 0 ) = [ 10 , 5 ] T . All the leader–follower topologies considered in this section satisfy Assumption 1. The weight of each edge between a leader and a follower is set to 1 and the weight of each edge among followers is shown in the pictures. The simulation results with different follower topology of the agents are exhibited as follows.
General Topology. The general follower topology satisfying Assumptions 1 and 2 is shown in Figure 1.
Select the gain parameters as γ 1 = 0.2 and γ 2 = 0.05 satisfying the conditions in Theorem 1, and all followers reach an agreement value in the convex hull spanned by leaders (see Figure 2 and Figure 3).
Symmetric Topology. The symmetric follower topology satisfying Assumptions 1 and 3 is shown in Figure 4.
Since the condition is easier to calculate, we are able to give a more specific range of the parameters. Under the given topology, we have γ 1 0.67 and we choose γ 1 as 0.2 here. With the condition γ 1 = 0.2 , we then have γ 2 > 0.1 . Finally, we choose γ 1 = 0.2 and γ 2 = 0.11 to guarantee that the conditions in Theorem 2 hold. Then, all followers reach the containment-rendezvous asymptotically (see Figure 5 and Figure 6).
Unconnected Topology. The unconnected follower topology shown in Figure 7 satisfies Assumption 1. It is evident that the topology can be divided into two connected follower topologies as shown in Figure 8.
The gain parameters are set as γ 1 = 0.2 and γ 2 = 0.1 satisfying the requirement in Theorem 3. It is seen from Figure 9 and Figure 10, all followers are divided into two group and followers in each group reach own agreement value in the convex hull spanned by leaders.

5. Conclusions

In this paper, containment-rendezvous problem of discrete-time first-order multi-agent systems is analyzed. The proposed control protocol includes a proportional term and an integral term. The proportional term ensures the realization of the containment, and the integral term guarantees the rendezvous. According to the frequency-domain analysis and numerical example, the effectiveness of our proposed protocol under the general connected follower topology is proved. For the symmetric and connected follower topology, a simpler convergence condition is presented. The containment control problem under unconnected follower topology is further discussed. Notably, the unconnected follower topology can be divided into several connected ones, so the followers still converge to the containment formed by all the leaders. Since the work in our paper is only a theoretic research and do not consider the trajectory of the agents, compared with some practical works [42,43,44], we will continue to study this question in a more practical way in our future work.

Author Contributions

Theoretical analysis, M.H.; simulation, M.H.; writing—original draft preparation, M.H., writing—review and edit, C.L.; supervision, C.L., funding acquisition, C.L., L.S. All authors have read and agreed to the published version of the manuscript.

Funding

Project supported by National Natural Science Foundation of China under Grants 61973139 and 61473138, the Fundamental Research Funds for the Central Universities under Grant JUSRP22014, the Natural Science Foundation of Jiangsu Province under BK20191286 and the Fundamental Research Funds for the Central Universities under 30920021139.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Olfati-Saber, R.; Murray, R.M. Consensus problems in networks of agents with switching topology and time-delays. IEEE Trans. Automat. Control 2004, 49, 1520–1533. [Google Scholar] [CrossRef] [Green Version]
  2. Olfati-Saber, R. Distributed Kalman filter with embedded consensus filters. In Proceedings of the 44th IEEE Conference on Decision and Control, Seville, Spain, 15 December 2005. [Google Scholar]
  3. Freeman, R.A.; Yang, P.; Lynch, K.M. Stability and convergence properties of dynamic average consensus estimators. In Proceedings of the 45th IEEE Conference on Decision and Control, San Diego, CA, USA, 13–15 December 2006; pp. 338–343. [Google Scholar]
  4. Olfati-Saber, R.; Fax, J.A.; Murray, R.M. Consensus and cooperation in networked multi-agent systems. Proc. IEEE 2007, 95, 215–233. [Google Scholar] [CrossRef] [Green Version]
  5. Ren, W.; Beard, R.W. Distributed Consensus in Multi-Vehicle Cooperative Control; Springer: London, UK, 2008; Volume 27, pp. 71–82. [Google Scholar]
  6. Liu, D.J.; Liu, C.L. Consensus problem of discrete-time second-order multi-agent network with communication delays. In Proceedings of the 2009 Third International Symposium on Intelligent Information Technology Application, Nanchang, China, 21–22 November 2009; pp. 340–344. [Google Scholar]
  7. Li, S.; Guo, Y. Distributed consensus filter on directed switching graphs. Int. J. Robust Nonlinear Control 2015, 25, 2019–2040. [Google Scholar] [CrossRef]
  8. Rezaei, M.H.; Kabiri, M.; Menhaj, M.B. Adaptive consensus for high-order unknown nonlinear multi-agent systems with unknown control directions and switching topologies. Inf. Sci. 2018, 459, 224–237. [Google Scholar] [CrossRef]
  9. Rezaei, M.H.; Menhaj, M.B. Adaptive output stationary average consensus for heterogeneous unknown linear multi-agent systems. IET Control Theory Appl. 2018, 12, 847–856. [Google Scholar] [CrossRef]
  10. Liu, C.L.; Gu, X.Y.; Shan, L. Average-consensus tracking of ramp inputs via second-order multi-agent systems. In Proceedings of the 2019 Chinese Control And Decision Conference (CCDC), Nanchang, China, 3–5 June 2019; pp. 3442–3447. [Google Scholar]
  11. Zheng, M.; Liu, C.L.; Liu, F. Average-consensus tracking of sensor network via distributed coordination control of heterogeneous multi-agent systems. IEEE Control Syst. Lett. 2019, 3, 132–137. [Google Scholar] [CrossRef]
  12. Ren, W.; Beard, R.W.; Atkins, E.M. A survey of consensus problems in multi-agent coordination. In Proceedings of the 2005 American Control Conference, Portland, OR, USA, 8–10 June 2005; pp. 1859–1864. [Google Scholar]
  13. Rezaei, M.H.; Menhaj, M.B. Stationary average consensus protocol for a class of heterogeneous high-order multi-agent systems with application for aircraft. Int. J. Syst. Sci. 2018, 49, 284–298. [Google Scholar] [CrossRef]
  14. Ren, W.; Beard, R.W. Consensus seeking in multiagent systems under dynamically changing interaction topologies. IEEE Trans. Autom. Control 2005, 50, 1859–1864. [Google Scholar]
  15. Jiang, W.; Wen, G.; Meng, Y.; Rahmani, A. Distributed adaptive time-varying formation tracking for linear multi-agent systems: A dynamic output approach. In Proceedings of the 2017 36th Chinese Control Conference (CCC), Dalian, China, 26–28 July 2017; pp. 8571–8576. [Google Scholar]
  16. Jiang, W.; Peng, Z.; Rahmani, A.; Hu, W.; Wen, G. Distributed consensus of linear MASs with an unknown leader via a predictive extended state observer considering input delay and disturbances. Neurocomputing 2018, 315, 465–475. [Google Scholar] [CrossRef] [Green Version]
  17. Cao, Y.C.; Ren, W.; Egerstedt, M. Distributed containment control with multiple stationary or dynamic leaders in fixed and switching directed networks. Automatica 2012, 48, 1586–1597. [Google Scholar] [CrossRef]
  18. Li, B.; Cheng, Z.Q.; Liu, Z.X.; Zhang, Q. Containment control of discrete-time multi-agent systems with multiple stationary leaders and time-delays. In Proceedings of the 2015 34th Chinese Control Conference(CCC), Hangzhou, China, 28–30 July 2015; pp. 7062–7066. [Google Scholar]
  19. Hua, Y.Z.; Dong, X.W.; Han, L.; Li, Q.D.; Zhang, R. Formation-containment tracking for general linear multi-agent systems with a tracking-leader of unknown control input. Syst. Control Lett. 2018, 122, 67–76. [Google Scholar] [CrossRef]
  20. Jiang, W. Fully Distributed Time-Varying Formation and Containment Control for Multi-Agent/Multi-Robot Systems. Ph.D. Thesis, Ecole Centrale de Lille, Villeneuve-d’Ascq, France, 2018. [Google Scholar]
  21. Wang, F.Y.; Liu, Z.X.; Chen, Z.Q. Containment control for second-order nonlinear multi-agent systems with aperiodically intermittent position measurements. J. Frankl. Inst. 2019, 356, 8706–8725. [Google Scholar] [CrossRef]
  22. Shi, L.; Xiao, Y.; Shao, J.L.; Zheng, W.X. Containment control of asynchronous discrete-time general linear multiagent systems with arbitrary network topology. IEEE Trans. Cybern. 2020, 50, 2546–2556. [Google Scholar] [CrossRef]
  23. Ding, Y.; Ren, W. Sampled-data containment control for double-integrator agents with dynamic leaders with nonzero inputs. Syst. Control Lett. 2020, 139, 104673. [Google Scholar] [CrossRef]
  24. Wang, F.Y.; Ni, Y.H.; Liu, Z.X. Containment control for general second-order multiagent systems with switched dynamics. IEEE Trans. Autom. Control 2020, 50, 550–560. [Google Scholar] [CrossRef]
  25. Liu, H.Y.; Xie, G.M.; Wang, L. Necessary and sufficient conditions for containment control of networked multi-agent systems. Automatica 2012, 48, 1415–1422. [Google Scholar] [CrossRef]
  26. Wang, Y.P.; Cheng, L.; Hou, Z.G.; Min, T.; Wang, M. Containment control of multi-agent systems in a noisy communication environment. Automatica 2014, 50, 1922–1928. [Google Scholar] [CrossRef]
  27. Rong, L.N.; Shen, H. Distributed PD-type protocol based containment control of multi-agent systems with input delays. J. Frankl. Inst. 2015, 352, 3600–3611. [Google Scholar] [CrossRef]
  28. Mu, X.W.; Liu, K. Containment control of single-integrator network with limited communication data rate. IEEE Trans. Autom. Control 2016, 61, 2232–2238. [Google Scholar] [CrossRef]
  29. Miao, G.Y.; Cao, J.D.; Alsaedi, A.; Alsaadif, F.E. Event-triggered containment control for multi-agent systems with constant time delays. J. Frankl. Inst. 2017, 354, 6956–6977. [Google Scholar] [CrossRef]
  30. Shao, J.L.; Shi, L.; Gong, L.S. Analysis of asynchronous containment control problem for discrete-time multi-agent systems. In Proceedings of the 2017 IEEE 2nd Advanced Information Technology, Electronic and Automation Control Conference (IAEAC), Chongqing, China, 25–26 March 2017; pp. 247–251. [Google Scholar]
  31. Han, T.; Li, J.; Guan, Z.H.; Cai, C.X.; Zhang, D.X.; He, D.X. Containment control of multi-agent systems via a disturbance observer-based approach. J. Frankl. Inst. 2019, 35, 2919–2933. [Google Scholar] [CrossRef]
  32. Liu, H.Y.; Chong, L.; Tan, M.; Hou, Z.G. Containment control of continuous-time linear multi-agent systems with aperiodic sampling. Automatica 2015, 57, 78–84. [Google Scholar] [CrossRef]
  33. Liu, T.F.; Qi, L.; Jiang, Z.P. Distributed containment control of multi-agent systems with velocity and acceleration saturations. Automatic 2020, 117, 108992. [Google Scholar] [CrossRef]
  34. Wang, D.; Wang, D.; Wang, W. Necessary and sufficient conditions for containment control of multi-agent systems with time delay. Automatica 2019, 103, 418–423. [Google Scholar] [CrossRef]
  35. Shan, L.; Liu, C.L. Average-consensus tracking of multi-agent systems with additional interconnecting agents. J. Frankl. Inst. 2018, 355, 8957–8970. [Google Scholar] [CrossRef]
  36. Chung, Y.F.; Kia, S.S. Distributed dynamic containment control over a strongly connected and weight-balanced digraph. IFAC-PapersOnLine 2019, 52, 25–30. [Google Scholar] [CrossRef]
  37. Huang, M.Y.; Liu, C.L. Containment-rendezvous control for first-order multi-agent systems. In Proceedings of the 2020 Chinese Automation Congress (CAC), Shanghai, China, 6–8 November 2020; pp. 4885–4890. [Google Scholar]
  38. Lin, Z.Y.; Francis, B.; Maggiore, M. Necessary and sufficient graphical conditions for formation control of unicycles. IEEE Trans. Autom. Control 2005, 50, 121–127. [Google Scholar]
  39. Parks, P.C.; Hahn, V. Stability Theory; Prentice Hall: Englewood Cliffs, NJ, USA, 1993. [Google Scholar]
  40. Lestas, I.; Vinnicombe, G. Scalable robustness for consensus protocols with heterogeneous dynamics. IFAC Proc. Vol. 2005, 38, 185–190. [Google Scholar] [CrossRef] [Green Version]
  41. Desoer, C.A.; Wang, Y.T. On the generalized Nyquist stability criterion. IEEE Trans. Autom. Control 1980, 50, 187–196. [Google Scholar] [CrossRef]
  42. Olfati-Saber, R. Flocking for multi-agent dynamic systems: Algorithms and theory. IEEE Trans. Autom. Control 2006, 51, 401–420. [Google Scholar] [CrossRef] [Green Version]
  43. Nguyen, T.; La, H.M.; Le, T.D.; Jafari, M. Formation control and obstacle avoidance of multiple rectangular agents with limited communication ranges. IEEE Trans. Control Netw. Syst. 2016, 4, 680–691. [Google Scholar] [CrossRef]
  44. Tran, V.P.; Garratt, M.; Petersen, I.R. Switching time-invariant formation control of a collaborative multi-agent system using negative imaginary systems theory. Control Eng. Pract. 2020, 95, 104245. [Google Scholar] [CrossRef]
Figure 1. General topology of the agents.
Figure 1. General topology of the agents.
Algorithms 14 00209 g001
Figure 2. Trajectories of agents with general topology.
Figure 2. Trajectories of agents with general topology.
Algorithms 14 00209 g002
Figure 3. States of followers with general topology.
Figure 3. States of followers with general topology.
Algorithms 14 00209 g003
Figure 4. Symmetric topology of agents.
Figure 4. Symmetric topology of agents.
Algorithms 14 00209 g004
Figure 5. Trajectories of agents with symmetric topology.
Figure 5. Trajectories of agents with symmetric topology.
Algorithms 14 00209 g005
Figure 6. States of followers with symmetric topology.
Figure 6. States of followers with symmetric topology.
Algorithms 14 00209 g006
Figure 7. Unconnected topology of the agents.
Figure 7. Unconnected topology of the agents.
Algorithms 14 00209 g007
Figure 8. Two connected parts in the unconnected topology.
Figure 8. Two connected parts in the unconnected topology.
Algorithms 14 00209 g008
Figure 9. Trajectories of agents with unconnected topology.
Figure 9. Trajectories of agents with unconnected topology.
Algorithms 14 00209 g009
Figure 10. States of followers with unconnected topology.
Figure 10. States of followers with unconnected topology.
Algorithms 14 00209 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, M.; Liu, C.; Shan, L. Containment Control of First-Order Multi-Agent Systems under PI Coordination Protocol. Algorithms 2021, 14, 209. https://doi.org/10.3390/a14070209

AMA Style

Huang M, Liu C, Shan L. Containment Control of First-Order Multi-Agent Systems under PI Coordination Protocol. Algorithms. 2021; 14(7):209. https://doi.org/10.3390/a14070209

Chicago/Turabian Style

Huang, Mingyang, Chenglin Liu, and Liang Shan. 2021. "Containment Control of First-Order Multi-Agent Systems under PI Coordination Protocol" Algorithms 14, no. 7: 209. https://doi.org/10.3390/a14070209

APA Style

Huang, M., Liu, C., & Shan, L. (2021). Containment Control of First-Order Multi-Agent Systems under PI Coordination Protocol. Algorithms, 14(7), 209. https://doi.org/10.3390/a14070209

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