This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

On Topological Conditions for Enabling Transient Control in Leader-follower Networks

Fei Chen fchen@kth.se    Dimos V. Dimarogonas dimos@kth.se Division of Decision and Control Systems, KTH Royal Institute of Technology, SE-100 44, Stockholm, Sweden
Abstract

We derive necessary and sufficient conditions for leader-follower multi-agent systems such that we can further apply prescribed performance control to achieve the desired formation while satisfying certain transient constraints. A leader-follower framework is considered in the sense that a group of agents with external inputs are selected as leaders in order to drive the group of followers in a way that the entire system can achieve target formation within certain prescribed performance transient bounds. We first derive necessary conditions on the leader-follower graph topology under which the target formation together with the prescribed performance guarantees can be fulfilled. Afterwards, the derived necessary conditions are extended to necessary and sufficient conditions for leader-follower formation control under transient constraints. Finally, the proposed results are illustrated with simulation examples.

keywords:
Leader-follower network, Multi-agent systems, Prescribed performance control, Necessary and sufficient condition
thanks: This work was supported by the ERC Consolidator Grant LEAFHOUND, the Swedish Research Council (VR) and the Knut och Alice Wallenberg Foundation (KAW).

,

1 Introduction

The problem of decentralized formation control (Fax & Murray, 2004) of multi-agent systems has been popular due to its wide applications in robotics and manufacturing. Some approaches further consider transient constraints for more complex tasks or certain prescribed system behavior, e.g., connectivity maintenance (Zavlanos et al., 2011) and collision avoidance (Chang et al., 2003). However, most existing approaches apply the designed local control strategies to all agents, which might sometimes be costly or redundant, motivating the need for leader-follower approaches. We thus propose here a systematic way on how to apply transient control in a leader-follower framework, in which the set of agents with advanced capabilities are selected as leaders to guide the whole agent group to the global equilibrium while fulfilling the transient constraints.

Control of leader-follower multi-agent systems under transient constraints has been pursued, e.g., Katsoukis & Rovithakis (2016); Dai et al. (2019); He et al. (2018); Wei et al. (2018). This standard leader-follower setting typically assumes a single leader which is treated as a reference for the remaining followers and all the follower controllers are designed for the purpose of tracking the leader. It has been extensively studied due to the close relation with the vast literature on multi-agent systems, specifically within a standard leader-follower setting. We instead consider a different leader-follower framework which considers arbitrary number of leaders with external inputs in addition to a first-order formation protocol. The remaining agents are followers only obeying the first-order formation protocol and have no knowledge or control design freedom for additional objectives, such as the transient behavior pursued here. Therefore, the followers only follow standard cooperative control protocols regardless of the objective of the whole team, while only the leader controllers are designed accordingly to guide the followers such that the target objective is achieved. Recent research in such a leader-follower setting mainly focuses on controllability of leader-follower networks. To name a few, Tanner (2004) derived conditions on the network topology, which ensures that the network can be controlled by particular members acting as leaders, while an extension of controllability conditions to leader-follower networks with double integrator dynamics was addressed by Goldin & Raisch (2010). Egerstedt et al. (2012); Rahmani et al. (2009) investigated necessary conditions for the controllability of the corresponding leader-follower networks using equitable partitions of graphs.  Sun et al. (2017) studied the controllability problem of leader-follower multi-agent systems defined by undirected signed graphs through almost equitable partitions and provided a necessary condition for the controllability of the network. The classes of essentially controllable, completely uncontrollable, and conditionally controllable graphs were discussed by Aguilar & Gharesifard (2014).  Yazıcıoğlu et al. (2016); Zhang et al. (2011, 2014) derived lower or upper bounds for controllable subspaces of leader-follower networks. Another category in the leader-follower setup considers leader selection problems to guarantee controllability (Yazicioğlu & Egerstedt, 2013) or to maximize a system performance metric (Fitch & Leonard, 2016; Patterson et al., 2016). Furthermore,  Pirani et al. (2017) studied the robustness of leader–follower consensus dynamics to disturbances and time delays.

As a comparison, we consider transient constraints rather than controllability in the leader-follower setting. Note that controllability is neither sufficient nor necessary to guarantee the fulfillment of the transient constraints, thus the topological conditions to guarantee transient constraints are expected to be different from those for controllability. In addition, the followers are solely indirectly guided through their dynamic couplings with the controlled leaders without any further control and knowledge of the transient control objective in hand, which makes the problem far more challenging than the existing work on transient constraints for (standard leader-follower) multi-agent systems, e.g., Katsoukis & Rovithakis (2016). Specifically, a main contribution of the paper in hand is that we consider control of leader-follower multi-agent systems under transient constraints, where the followers do not have any knowledge of the prescribed performance bounds, while only the leader controllers are designed to guide the followers such that the target formation is achieved within the transient constraints. Such leader-follower framework finds application in various domains, including multi-vehicle platooning. This involves designing leader vehicle controllers to ensure collision avoidance and maintain connectivity for the entire platoon. In addition, it encompasses multi-robot coordination, involving formation control or even more complex task planning for robot teams under spatiotemporal constraints, where only the controllers of the leader robots are developed. Furthermore, the framework enables collaborative manipulation, involving grasping and transporting objects in a collaborative leader-follower manner while satisfying the transient constraints. All these tasks are assigned to the available leader agents, demonstrating both a distributed and scalable nature. Additionally we investigate the topological conditions on the leader-follower networks such that the target objective can be achieved within the prescribed transient constraints, which offers advantages in various topics, including addressing leader selection problems amidst transient constraints, studying the network’s robustness in terms of agent failures, and exploring network reconfigurations. The technical challenges that arise are due to the consideration of transient constraints which are only known to the leaders, and also the combination of non-specific topologies, leader amount and leader positions. Preliminary results of leader-follower formation control with prescribed performance guarantees have been presented by Chen & Dimarogonas (2020b, a), where  Chen & Dimarogonas (2020a) have also briefly discussed necessary conditions on the leader-follower graph topology. In this paper, we present more details on the topological conditions while the necessary conditions  (Chen & Dimarogonas, 2020a) are further extended to necessary and sufficient conditions on the leader-follower graph topology such that target formation within certain prescribed performance bounds can be achieved. This introduces for the first time a new methodology for leader selection to guarantee target formation while satisfying the transient constraints within the leader-follower framework.

The contributions and novelty of the paper can be thus summarized as: i) we consider control of leader-follower multi-agent systems under transient constraints, where the followers do not have any knowledge of the prescribed performance bounds; ii) we derive necessary and sufficient conditions on the leader-follower graph topology such that the target formation together with the prescribed transient behavior can be fulfilled for the considered leader-follower multi-agent system.

2 Preliminaries and Problem Statement

2.1 Graph Theory

An undirected graph (Mesbahi & Egerstedt, 2010) is defined as 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}) with the vertices set 𝒱={1,2,,n}\mathcal{V}=\{1,2,\dots,n\} and the edges set ={(i,j)𝒱×𝒱j𝒩i}\mathcal{E}=\{(i,j)\in\mathcal{V}\times\mathcal{V}\mid j\in\mathcal{N}_{i}\} indexed by e1,,eme_{1},\dots,e_{m}. m=||m=|\mathcal{E}| is the number of edges and 𝒩i\mathcal{N}_{i} denotes the neighbourhood of agent ii. For an edge ek=(i,j)e_{k}=(i,j), 𝐯(ek)={i,j}\mathbf{v}(e_{k})=\{i,j\} is the set composed of the vertices of eke_{k}. The adjacency matrix 𝔸\mathbb{A} of 𝒢\mathcal{G} is the n×nn\times n symmetric matrix whose elements aija_{ij} are given by aij=1a_{ij}=1, if (i,j)(i,j)\in\mathcal{E}, and aij=0a_{ij}=0, otherwise. The degree of vertex ii is defined as d(i)=j𝒩iaijd(i)={\sum}_{j\in\mathcal{N}_{i}}a_{ij}. Then the degree matrix is Δ=diag(d(1),d(2),,d(n))\Delta=\mbox{diag}(d(1),d(2),\dots,d(n)). The graph Laplacian of 𝒢\mathcal{G} is L=Δ𝔸L=\Delta-\mathbb{A}. A path is a sequence of edges connecting two distinct vertices. A graph is connected if there exists a path between any pair of vertices. By assigning an orientation to each edge of 𝒢\mathcal{G} the incidence matrix D=D(𝒢)=[dij]n×mD=D(\mathcal{G})=[d_{ij}]\in\mathbb{R}^{n\times m} is defined. The rows of DD are indexed by the vertices and the columns are indexed by the edges with dij=1d_{ij}=1 if the vertex ii is the head of the edge (i,j)(i,j), dij=1d_{ij}=-1 if the vertex ii is the tail of the edge (i,j)(i,j) and dij=0d_{ij}=0 otherwise. The graph Laplacian of 𝒢\mathcal{G} is described as L=DDTL=DD^{T}. Le=DTDL_{e}=D^{T}D is the so-called edge Laplacian (Zelazo & Mesbahi, 2011) and we let cijc_{ij} denote the elements of LeL_{e}.

2.2 System Description and Proposed Control

In this work, we consider a multi-agent system under the communication graph 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}) with vertices 𝒱={1,2,,n}\mathcal{V}=\{1,2,\dots,n\}. Suppose that the first nfn_{f} agents are followers while the last nln_{l} agents are leaders with 𝒱={1,2,,nf}\mathcal{V}_{\mathscr{F}}=\{1,2,\dots,n_{f}\}, 𝒱={nf+1,nf+2,,nf+nl}\mathcal{V}_{\mathscr{L}}=\{n_{f}+1,n_{f}+2,\dots,n_{f}+n_{l}\} and n=nf+nl,𝒱𝒱=,𝒱𝒱=𝒱n=n_{f}+n_{l},\mathcal{V}_{\mathscr{F}}\cap\mathcal{V}_{\mathscr{L}}=\emptyset,\mathcal{V}_{\mathscr{F}}\cup\mathcal{V}_{\mathscr{L}}=\mathcal{V}. For an agent ii with d(i)=1d(i)=1, we call it an end agent or an end node. Moreover, if it further holds that i𝒱i\in\mathcal{V}_{\mathscr{L}} or i𝒱i\in\mathcal{V}_{\mathscr{F}}, we call it an end leader and an end follower, respectively. Here and later on, the subscript ``",``"``\mathscr{F}",``\mathscr{L}" stands for follower and leader vertices set, respectively. We also denote the follower-follower edge set as ={ek1km,k,ek=(i,j),i,j𝒱}\mathcal{E}_{\mathscr{F}\mathscr{F}}=\{e_{k}\mid 1\leq k\leq m,k\in\mathbb{Z},e_{k}=(i,j),i,j\in\mathcal{V}_{\mathscr{F}}\}, the leader-leader edge set as ={ek1km,k,ek=(i,j),i,j𝒱}\mathcal{E}_{\mathscr{L}\mathscr{L}}=\{e_{k}\mid 1\leq k\leq m,k\in\mathbb{Z},e_{k}=(i,j),i,j\in\mathcal{V}_{\mathscr{L}}\} and the leader-follower edge set as =()\mathcal{E}_{\mathscr{L}\mathscr{F}}=\mathcal{E}\setminus(\mathcal{E}_{\mathscr{F}\mathscr{F}}\cup\mathcal{E}_{\mathscr{L}\mathscr{L}}). In a similar manner, we can define 𝒱sup𝒱,𝒱sup𝒱,sup,sup,sup\mathcal{V}_{\mathscr{F}}^{\text{sup}}\subseteq\mathcal{V}_{\mathscr{F}},\mathcal{V}_{\mathscr{L}}^{\text{sup}}\subseteq\mathcal{V}_{\mathscr{L}},\mathcal{E}_{\mathscr{F}\mathscr{F}}^{\text{sup}}\subseteq\mathcal{E}_{\mathscr{F}\mathscr{F}},\mathcal{E}_{\mathscr{L}\mathscr{L}}^{\text{sup}}\subseteq\mathcal{E}_{\mathscr{L}\mathscr{L}},\mathcal{E}_{\mathscr{L}\mathscr{F}}^{\text{sup}}\subseteq\mathcal{E}_{\mathscr{L}\mathscr{F}} with respect to the subgraphs 𝒢sup=(𝒱sup,sup)\mathcal{G}^{\text{sup}}=(\mathcal{V}^{\text{sup}},\mathcal{E}^{\text{sup}}) of 𝒢\mathcal{G} over different superscripts, e.g., we use sup{f,p,,,}{\text{sup}}\in\{f,p,\star,\prime,\cdots\} later on.

Let pip_{i}\in\mathbb{R} be the position of agent ii, where we only consider the one dimensional case, without loss of generality. Specifically, the results can be extended to higher dimensions with Kronecker product. The target relative position-based formation is described as follows:

:={ppipj=pijdes,(i,j)},\mathcal{F}:=\{p\mid p_{i}-p_{j}=p^{des}_{ij},(i,j)\in\mathcal{E}\}, (1)

where pijdes:=pidespjdes,(i,j)p^{des}_{ij}:=p^{des}_{i}-p^{des}_{j},(i,j)\in\mathcal{E} is the desired relative position between agent ii and agent jj, which is constant and denoted as the difference between the absolute desired positions pides,pjdesp^{des}_{i},p^{des}_{j}\in\mathbb{R}. The state evolution of each agent i𝒱i\in\mathcal{V} is governed by the following dynamics:

p˙i=j𝒩i(pipjpijdes)+biui,\dot{p}_{i}=-\sum\limits_{j\in\mathcal{N}_{i}}(p_{i}-p_{j}-p^{des}_{ij})+b_{i}u_{i}, (2)

with bi=1b_{i}=1 if i𝒱i\in\mathcal{V}_{\mathscr{L}}, and bi=0b_{i}=0 if i𝒱i\in\mathcal{V}_{\mathscr{F}}. This means that followers are governed by the first-order formation protocol, while leaders are governed by the first-order formation protocol with an assigned external input. By stacking (2), the dynamics of the leader-follower multi-agent system are rewritten as:

Σ:p˙=L(ppdes)+Bu,\Sigma:\dot{p}=-L(p-p^{des})+Bu, (3)

where p=[p1,,pn]T,pdes=[p1des,,pndes]Tnp=[p_{1},\dots,p_{n}]^{T},p^{des}=[p^{des}_{1},\dots,p^{des}_{n}]^{T}\in\mathbb{R}^{n}, u=[unf+1,,unf+nl]Tnlu=[u_{n_{f}+1},\dots,u_{{n_{f}+n_{l}}}]^{T}\in\mathbb{R}^{n_{l}} and B=[0nf×nlInl].B=\left[\begin{smallmatrix}0_{n_{f}\times n_{l}}\\ I_{n_{l}}\end{smallmatrix}\right]. Denote p¯=[p¯1,,p¯m]T\bar{p}=[\bar{p}_{1},\dots,\bar{p}_{m}]^{T}, p¯des=[p¯1des,,p¯mdes]Tm\bar{p}^{des}=[\bar{p}^{des}_{1},\dots,\bar{p}^{des}_{m}]^{T}\in\mathbb{R}^{m} as the respective stack vector of relative positions and target relative positions between the pair of communication agents for the edge (i,j)=ek(i,j)=e_{k}\in\mathcal{E}, where p¯kpij=pipj,p¯kdespijdes=pidespjdes,k=1,2,,m\bar{p}_{k}\triangleq p_{ij}=p_{i}-p_{j},\bar{p}^{des}_{k}\triangleq p^{des}_{ij}=p^{des}_{i}-p^{des}_{j},k=1,2,\dots,m. Accordingly, x¯=p¯p¯des=[x¯1,,x¯m]T\bar{x}=\bar{p}-\bar{p}^{des}=[\bar{x}_{1},\dots,\bar{x}_{m}]^{T} is denoted as the shifted relative position vector. The incidence matrix DD can be decomposed by the rows into the first nfn_{f} and the remaining last nln_{l} rows, i.e., D=[DTDT]TD=\begin{bmatrix}D_{\mathscr{F}}^{T}&D_{\mathscr{L}}^{T}\end{bmatrix}^{T}. Multiplying with DTD^{T} on both sides of (3), the dynamics (3) can be reorganized into the edge space as

Σe:x¯˙=Lex¯+DTu,\Sigma_{e}:\dot{\bar{x}}=-L_{e}\bar{x}+D_{\mathscr{L}}^{T}u, (4)

The aim of prescribed performance control (PPC) (Bechlioulis & Rovithakis, 2008) is to prescribe the evolution of the relative position p¯i(t)\bar{p}_{i}(t) within some predefined region described as

p¯idesρx¯i(t)<p¯i(t)<p¯ides+ρx¯i(t),\bar{p}^{des}_{i}-\rho_{\bar{x}_{i}}(t)<\bar{p}_{i}(t)<\bar{p}^{des}_{i}+\rho_{\bar{x}_{i}}(t), (5)

or equivalently, to prescribe x¯i(t)\bar{x}_{i}(t) within ρx¯i(t)<x¯i(t)<ρx¯i(t).-\rho_{\bar{x}_{i}}(t)<\bar{x}_{i}(t)<\rho_{\bar{x}_{i}}(t). Here ρx¯i(t):++{0},i=1,2,,m\rho_{\bar{x}_{i}}(t):\mathbb{R}_{+}\rightarrow\mathbb{R}_{+}\setminus\{0\},i=1,2,\dots,m are positive, smooth and strictly decreasing performance functions and one example choice is ρx¯i(t)=(ρx¯i0ρx¯i)elx¯it+ρx¯i\rho_{\bar{x}_{i}}(t)=(\rho_{\bar{x}_{i0}}-\rho_{\bar{x}_{i\infty}})e^{-l_{\bar{x}_{i}}t}+\rho_{\bar{x}_{i\infty}} with ρx¯i0,ρx¯i\rho_{\bar{x}_{i0}},\rho_{\bar{x}_{i\infty}} and lx¯il_{\bar{x}_{i}} being positive parameters. We further define the modulated error as x¯^i(t)=x¯i(t)ρx¯i(t)\hat{\bar{x}}_{i}(t)=\frac{\bar{x}_{i}(t)}{\rho_{\bar{x}_{i}}(t)} and the corresponding prescribed performance region as 𝒟x¯i{x¯^i:x¯^i(1,1)}.\mathcal{D}_{\bar{x}_{i}}\triangleq\{\hat{\bar{x}}_{i}:\hat{\bar{x}}_{i}\in(-1,1)\}. Then the modulated error is transformed through a function Tx¯iT_{\bar{x}_{i}} that defines the smooth and strictly increasing mapping Tx¯i:𝒟x¯iT_{\bar{x}_{i}}:\mathcal{D}_{\bar{x}_{i}}\rightarrow\mathbb{R}, Tx¯i(0)=0T_{\bar{x}_{i}}(0)=0. Here, one example choice is Tx¯i(x¯^i)=ln(1+x¯^i1x¯^i).T_{\bar{x}_{i}}(\hat{\bar{x}}_{i})=\ln\left(\frac{1+\hat{\bar{x}}_{i}}{1-\hat{\bar{x}}_{i}}\right). The transformed error is then defined as εx¯i(x¯^i)=Tx¯i(x¯^i).\varepsilon_{\bar{x}_{i}}(\hat{\bar{x}}_{i})=T_{\bar{x}_{i}}(\hat{\bar{x}}_{i}). Differentiating εx¯i(x¯^i)\varepsilon_{\bar{x}_{i}}(\hat{\bar{x}}_{i}) with respect to time, we derive ε˙x¯i(x¯^i)=𝒥Tx¯i(x¯^i,t)(x¯˙i+αx¯i(t)x¯i),\dot{\varepsilon}_{\bar{x}_{i}}(\hat{\bar{x}}_{i})=\mathcal{J}_{T_{\bar{x}_{i}}}(\hat{\bar{x}}_{i},t)(\dot{\bar{x}}_{i}+\alpha_{\bar{x}_{i}}(t)\bar{x}_{i}), where 𝒥Tx¯i(x¯^i,t)Tx¯i(x¯^i)x¯^i1ρx¯i(t)\mathcal{J}_{T_{\bar{x}_{i}}}(\hat{\bar{x}}_{i},t)\triangleq\frac{\partial T_{\bar{x}_{i}}(\hat{\bar{x}}_{i})}{\partial\hat{\bar{x}}_{i}}\frac{1}{\rho_{\bar{x}_{i}}(t)} and αx¯i(t)ρ˙x¯i(t)ρx¯i(t)\alpha_{\bar{x}_{i}}(t)\triangleq-\frac{\dot{\rho}_{\bar{x}_{i}}(t)}{\rho_{\bar{x}_{i}}(t)} are the normalized Jacobian of the transformed function Tx¯iT_{\bar{x}_{i}} and the normalized derivative of the performance function, respectively. The basic idea of PPC is to verify that the transformed error εx¯i(x¯^i)\varepsilon_{\bar{x}_{i}}(\hat{\bar{x}}_{i}) is bounded, which in turn results in the satisfaction of (5(Bechlioulis & Rovithakis, 2008). For the edge dynamics (4), the following PPC control strategy is proposed by Chen & Dimarogonas (2020b, a):

uj=iΦjgx¯i𝒥Tx¯i(x¯^i,t)εx¯i(x¯^i),j𝒱,u_{j}=-\sum\limits_{i\in\Phi_{j}}g_{\bar{x}_{i}}\mathcal{J}_{T_{\bar{x}_{i}}}(\hat{\bar{x}}_{i},t)\varepsilon_{\bar{x}_{i}}(\hat{\bar{x}}_{i}),\hskip 14.22636ptj\in\mathcal{V}_{\mathscr{L}}, (6)

where Φj={i|(j,k)=i,k𝒩j}\Phi_{j}=\{i|(j,k)=i,k\in\mathcal{N}_{j}\}, i.e., the set of all the edges that include agent j𝒱j\in\mathcal{V}_{\mathscr{L}} as a node, and gx¯ig_{\bar{x}_{i}} is a positive scalar gain to be appropriately tuned. It is concluded by Chen & Dimarogonas (2020b, a) that the leader-follower multi-agent system (3) can achieve the target formation \mathcal{F} as in (1) while satisfying (5) by applying (6) under some assumptions on the graph topology, which will be further discussed in this paper. We refer the readers to Chen & Dimarogonas (2020b, a) for the detailed convergence analysis.

2.3 Problem Statement

In this work, we focus on investigating the leader-follower graph topology and the following problem is formulated.

Problem 1.

Derive the necessary and sufficient conditions on the leader-follower communication graph 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}) such that the leader-follower multi-agent system (3) achieves the target formation \mathcal{F} as in (1) while satisfying (5) by applying the proposed control strategy (6).

3 Necessary Conditions on Graph Topology

In this section, we derive necessary conditions on the graph topology for both tree graphs and general graphs with cycles such that under these conditions we can design the leaders to achieve the target formation with prescribed performance guarantees.

We first discuss the tree graphs and then the results for general graphs with cycles are built based on the results of tree graphs. We first define a leaderless graph 𝒢f=(𝒱f,f)\mathcal{G}^{f}=(\mathcal{V}^{f},\mathcal{E}^{f}) with only followers, i.e., 𝒱f=𝒱f\mathcal{V}_{\mathscr{F}}^{f}=\mathcal{V}^{f}, 𝒱f=\mathcal{V}_{\mathscr{L}}^{f}=\emptyset and f=f\mathcal{E}_{\mathscr{F}\mathscr{F}}^{f}=\mathcal{E}^{f}. The insight here is that we would like to analyze how the leader-follower multi-agent system described by the graph 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}) behaves when it contains 𝒢f\mathcal{G}^{f} as an induced subgraph. The definition of subgraph and induced subgraph is given as follows:

Definition 2.

(Subgraph and induced subgraph.) A graph 𝒢=(𝒱,)\mathcal{G}^{\prime}=(\mathcal{V}^{\prime},\mathcal{E}^{\prime}) is a subgraph of the graph 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}) if 𝒱𝒱\mathcal{V}^{\prime}\subseteq\mathcal{V} and \mathcal{E}^{\prime}\subseteq\mathcal{E}; for any i𝒱i\in\mathcal{V}_{\mathscr{F}}^{\prime} we have i𝒱i\in\mathcal{V}_{\mathscr{F}} and for any i𝒱i\in\mathcal{V}_{\mathscr{L}}^{\prime} we have i𝒱i\in\mathcal{V}_{\mathscr{L}}. A subgraph 𝒢\mathcal{G}^{\prime} of 𝒢\mathcal{G} is an induced subgraph of 𝒢\mathcal{G}, denoted as 𝒢𝒢\mathcal{G}^{\prime}\subseteq\mathcal{G} if it further holds that for any edge (i,j),i,j𝒱(i,j)\in\mathcal{E},i,j\in\mathcal{V}^{\prime}, we have (i,j)(i,j)\in\mathcal{E}^{\prime}.

Definition 3.

(Path subgraph.) A path pp of the graph 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}) is called a path subgraph p=(𝒱p,p)p=(\mathcal{V}^{p},\mathcal{E}^{p}) of 𝒢\mathcal{G} such that 𝒱p𝒱\mathcal{V}^{p}\subseteq\mathcal{V} and p\mathcal{E}^{p}\subseteq\mathcal{E}, for any i𝒱pi\in\mathcal{V}_{\mathscr{F}}^{p} we have i𝒱i\in\mathcal{V}_{\mathscr{F}} and for any i𝒱pi\in\mathcal{V}_{\mathscr{L}}^{p} we have i𝒱i\in\mathcal{V}_{\mathscr{L}};

From now on, we denote x¯\bar{x} as the edge state of 𝒢f\mathcal{G}^{f} and LeL_{e} as the edge Laplacian of 𝒢f\mathcal{G}^{f}. We know that the edge dynamics of 𝒢f\mathcal{G}^{f} are simply described as x¯˙=Lex¯\dot{\bar{x}}=-L_{e}\bar{x} since the leader set 𝒱f\mathcal{V}^{f}_{\mathscr{L}} is empty. Denote each column (corresponding to an edge) of the incidence matrix of 𝒢f\mathcal{G}^{f} by the vector eie_{i}. Then (Le)ij=eiTej=cij=2(L_{e})_{ij}=e_{i}^{T}e_{j}=c_{ij}=2 if i=ji=j; cij=0c_{ij}=0 if ei,eje_{i},e_{j} share no nodes; cij=1c_{ij}=1 if ei,eje_{i},e_{j} share a single node and have the same direction with respect to the sharing node; cij=1c_{ij}=-1 if ei,eje_{i},e_{j} share a single node but have different direction with respect to the sharing node (Zelazo & Mesbahi, 2011). Based on these simple rules, we can derive the dynamics in the edge space. Example 4 elucidates such a derivation. We define the neighbors of edge eie_{i} as 𝒩(ei):={ej|eiTej|=1}.\mathcal{N}(e_{i}):=\{e_{j}\mid|e_{i}^{T}e_{j}|=1\}.

Example 4.

We consider a leaderless tree graph 𝒢f\mathcal{G}^{f} as in Fig. 1 where we assign arbitrary directions for all the edges in order to derive the corresponding incidence matrix DD and the edge Laplacian LeL_{e} (note that we still focus on undirected graphs). Then the edge dynamics of the system described by 𝒢f\mathcal{G}^{f} are x¯˙=Lex¯=DTDx¯\dot{\bar{x}}=-L_{e}\bar{x}=-D^{T}D\bar{x} with x¯=[x¯1x¯2x¯3x¯4]T\bar{x}=\begin{bmatrix}\bar{x}_{1}&\bar{x}_{2}&\bar{x}_{3}&\bar{x}_{4}\end{bmatrix}^{T} and D,LeD,L_{e} matrices as follows:

D=[11111000010000100001],Le=[2111121111211112].D=\left[\begin{smallmatrix}-1&1&1&1\\ 1&0&0&0\\ 0&-1&0&0\\ 0&0&-1&0\\ 0&0&0&-1\end{smallmatrix}\right],L_{e}=\left[\begin{smallmatrix}2&-1&-1&-1\\ -1&2&1&1\\ -1&1&2&1\\ -1&1&1&2\end{smallmatrix}\right].

So the dynamics of the edge e1e_{1} are x¯˙1=2x¯1+x¯2+x¯3+x¯4\dot{\bar{x}}_{1}=-2\bar{x}_{1}+\bar{x}_{2}+\bar{x}_{3}+\bar{x}_{4} according to first row of LeL_{e}. Based on above discussions, we can also derive x¯˙1\dot{\bar{x}}_{1} by directly checking the graph only. In Fig. 1, edge e1e_{1} has three neighboring edges e2,e3,e4e_{2},e_{3},e_{4} since they share the node indexed by 11. Then we denote x¯˙1=ax¯1+bx¯2+cx¯3+dx¯4\dot{\bar{x}}_{1}=a\bar{x}_{1}+b\bar{x}_{2}+c\bar{x}_{3}+d\bar{x}_{4} with [abcd]\begin{bmatrix}-a&-b&-c&-d\end{bmatrix} being the first row of LeL_{e}. We know that a=2a=-2 since the diagonal entry of LeL_{e} is always 22. b=c=d=1b=c=d=1 since edge e1e_{1} shares node 11 with edges e2,e3,e4e_{2},e_{3},e_{4} butt have different directions with respect to node 11 (e1e_{1} points outward node 11, e2,e3,e4e_{2},e_{3},e_{4} all point inward node 11, and we assign 1-1 to the corresponding entries of LeL_{e}).

Refer to caption
Figure 1: Example of dynamics for a tree graph via edge Laplacian.

Then the following lemma is proposed for tree graphs and acts as a basis for the later extended results for general graphs with cycles.

Lemma 5.

Consider the leader-follower multi-agent system Σ\Sigma described by the tree graph 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}). A necessary condition on 𝒢\mathcal{G} under which we can design the leaders using (6) to achieve the target formation \mathcal{F} as in (1) while satisfying (5) is that every leaderless graph 𝒢f=(𝒱f,f)\mathcal{G}^{f}=(\mathcal{V}^{f},\mathcal{E}^{f}), such that 𝒢f𝒢\mathcal{G}^{f}\subseteq\mathcal{G}, should satisfy

|𝒩(ei)|2,eif.|\mathcal{N}(e_{i})|\leq 2,\hskip 8.53581pt\forall e_{i}\in\mathcal{E}^{f}. (7)
{pf}

The proof uses contradiction based on the entries of the edge Laplacian LeL_{e} of 𝒢f\mathcal{G}^{f}. Suppose 𝒢f\mathcal{G}^{f} is a leaderless induced subgraph of 𝒢\mathcal{G} and there exists eife_{i}\in\mathcal{E}^{f} satisfying |𝒩(ei)|3|\mathcal{N}(e_{i})|\geq 3. Without loss of generality, let us assume that e1fe_{1}\in\mathcal{E}^{f} satisfies 𝒩(e1)={e2,e3,e4}\mathcal{N}(e_{1})=\{e_{2},e_{3},e_{4}\}, thus |𝒩(ei)|=3|\mathcal{N}(e_{i})|=3. Suppose that e2,e3,e4e_{2},e_{3},e_{4} all share a single node with e1e_{1} but with different directions as shown in Fig. 1. This can be assumed without loss of generality since we can assign arbitrary directions to the edges. Then the state evolution of e1e_{1} is derived as x¯˙1=2x¯1+x¯2+x¯3+x¯4\dot{\bar{x}}_{1}=-2\bar{x}_{1}+\bar{x}_{2}+\bar{x}_{3}+\bar{x}_{4} according to Example 1. We can see that when all x¯i,i=1,2,3,4\bar{x}_{i},i=1,2,3,4 are initialised arbitrarily close to the prescribed performance bound ρx¯i0=ρ0,i=1,2,3,4\rho_{\bar{x}_{i0}}=\rho_{0},i=1,2,3,4, then x¯˙1=ρ0>0\dot{\bar{x}}_{1}=\rho_{0}>0, thus x¯1\bar{x}_{1} will continue evolving to violate the performance bound. This leads to a contradiction since no matter how we design the leaders in 𝒢\mathcal{G}, x¯1\bar{x}_{1} will always increase to violate the bound. Hence, we can conclude that 𝒢\mathcal{G} should not contain a leaderless induced subgraph 𝒢f=(𝒱f,f)\mathcal{G}^{f}=(\mathcal{V}^{f},\mathcal{E}^{f}) with 𝒱Ff=𝒱f\mathcal{V}_{F}^{f}=\mathcal{V}^{f} such that there exists eife_{i}\in\mathcal{E}^{f} satisfying |𝒩(ei)|3|\mathcal{N}(e_{i})|\geq 3. Or in other words, every leaderless graph 𝒢f=(𝒱f,f)\mathcal{G}^{f}=(\mathcal{V}^{f},\mathcal{E}^{f}), such that 𝒢f𝒢\mathcal{G}^{f}\subseteq\mathcal{G}, should satisfy |𝒩(ei)|2,eif.|\mathcal{N}(e_{i})|\leq 2,\forall e_{i}\in\mathcal{E}^{f}.∎ Next, based on Lemma 1, we derive a necessary condition for graphs with cycles. We denote here (𝒞)\mathcal{E}(\mathcal{C}) as the edge set of the cycle 𝒞\mathcal{C} with cardinality |(𝒞)||\mathcal{E}(\mathcal{C})|. We then perform the following graph decomposition which we call complete decomposition.

Definition 6.

(Complete decomposition.) A graph 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}) is decomposed with respect to the edge eie_{i}\in\mathcal{E} as 𝒢:=𝒞ei𝒳i𝒞ei𝒫i\mathcal{G}:=\cup_{\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}}\mathcal{C}_{e_{i}}\cup\mathcal{P}_{i}, where 𝒳i:={𝒞eiei𝒞ei}\mathcal{X}_{i}:=\{\mathcal{C}_{e_{i}}\mid e_{i}\in\mathcal{C}_{e_{i}}\} is the cycle set composed of all the cycles 𝒞ei\mathcal{C}_{e_{i}} in 𝒢\mathcal{G} that contain eie_{i} as an edge, and satisfy:

  • for every pair 𝒞eia,𝒞eib𝒳i\mathcal{C}_{e_{i}}^{a},\mathcal{C}_{e_{i}}^{b}\in\mathcal{X}_{i}, (𝒩(ei)𝒞eia)(𝒩(ei)𝒞eib)=(\mathcal{N}(e_{i})\cap\mathcal{C}_{e_{i}}^{a})\cap(\mathcal{N}(e_{i})\cap\mathcal{C}_{e_{i}}^{b})=\emptyset.

  • for every 𝒞ei𝒳i\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}, there does not exist a cycle 𝒞\mathcal{C} of 𝒢\mathcal{G} such that ei𝒞,(𝒩(ei)𝒞ei)(𝒩(ei)𝒞)e_{i}\in\mathcal{C},(\mathcal{N}(e_{i})\cap\mathcal{C}_{e_{i}})\cap(\mathcal{N}(e_{i})\cap\mathcal{C})\neq\emptyset, and |(𝒞)|<|(𝒞ei)||\mathcal{E}(\mathcal{C})|<|\mathcal{E}(\mathcal{C}_{e_{i}})|,

and where 𝒫i:={ekek𝒞ei,𝒞ei𝒳i}\mathcal{P}_{i}:=\{e_{k}\mid e_{k}\notin\mathcal{C}_{e_{i}},\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}\} is the set of the edges that do not belong to any cycle in 𝒳i\mathcal{X}_{i}. Then, we call this decomposition a complete decomposition of 𝒢\mathcal{G} with respect to the edge eie_{i}\in\mathcal{E}.

Next, we elucidate the complete decomposition with the following example.

Example 7.

Consider two leaderless graphs as shown in Fig.2. In the left figure, the complete decomposition with respect to the edge e1e_{1} is 𝒢:=𝒞e1𝒳1𝒞e1𝒫1\mathcal{G}:=\cup_{\mathcal{C}_{e_{1}}\in\mathcal{X}_{1}}\mathcal{C}_{e_{1}}\cup\mathcal{P}_{1}, where the cycle set 𝒳1\mathcal{X}_{1} includes only the cycle on the top, i.e., the cycle composed by the edges e1,e2,e3e_{1},e_{2},e_{3}, and 𝒫1={e4,e5,e6}\mathcal{P}_{1}=\{e_{4},e_{5},e_{6}\}. The complete decomposition with respect to the edge e3e_{3} is 𝒢:=𝒞e3𝒳3𝒞e3𝒫3\mathcal{G}:=\cup_{\mathcal{C}_{e_{3}}\in\mathcal{X}_{3}}\mathcal{C}_{e_{3}}\cup\mathcal{P}_{3}, where the cycle set 𝒳3\mathcal{X}_{3} includes the cycle composed by the edges e1,e2,e3e_{1},e_{2},e_{3} and the cycle composed by the edges e3,e4,e5,e6e_{3},e_{4},e_{5},e_{6}, and thus 𝒫3=\mathcal{P}_{3}=\emptyset. The second condition of the cycle set in Definition 3 implies that we do not consider the large cycle 𝒞\mathcal{C}^{\prime} composed of the edges e1,e2,e6,e5,e4e_{1},e_{2},e_{6},e_{5},e_{4} as in the cycle set 𝒳1\mathcal{X}_{1} or 𝒳3\mathcal{X}_{3}. For instance, consider e1e_{1} again, and let us assume 𝒞𝒳1\mathcal{C}^{\prime}\in\mathcal{X}_{1}. Then there exists a cycle 𝒞\mathcal{C} which is composed by e1,e2,e3e_{1},e_{2},e_{3} such that e1𝒞,(𝒩(e1)𝒞)(𝒩(e1)𝒞)={e2}e_{1}\in\mathcal{C},(\mathcal{N}(e_{1})\cap\mathcal{C}^{\prime})\cap(\mathcal{N}(e_{1})\cap\mathcal{C})=\{e_{2}\}\neq\emptyset since 𝒩(e1)={e2,e3,e4}\mathcal{N}(e_{1})=\{e_{2},e_{3},e_{4}\}. In addition |(𝒞)|=3<|(𝒞)|=5|\mathcal{E}(\mathcal{C})|=3<|\mathcal{E}(\mathcal{C}^{\prime})|=5, which contradicts the second condition of Definition 2. Therefore, we can conclude that 𝒞𝒳1\mathcal{C}^{\prime}\notin\mathcal{X}_{1}. In the right figure, the complete decomposition with respect to the edge e4e_{4} is 𝒢:=𝒞e4𝒳4𝒞e4𝒫4\mathcal{G}:=\cup_{\mathcal{C}_{e_{4}}\in\mathcal{X}_{4}}\mathcal{C}_{e_{4}}\cup\mathcal{P}_{4}, where the cycle set 𝒳4\mathcal{X}_{4} includes the cycle composed by the edges e1,e4,e5e_{1},e_{4},e_{5} and the cycle composed by the edges e3,e4,e6e_{3},e_{4},e_{6}; 𝒫4={e2,e7}\mathcal{P}_{4}=\{e_{2},e_{7}\} and thus 𝒩(e4)𝒫4={e7}\mathcal{N}(e_{4})\cap\mathcal{P}_{4}=\{e_{7}\}.

Refer to caption
Figure 2: Complete decomposition examples

The complete decomposition decomposes a large scale graph into certain cycles together with the remaining edges that do not belong to any cycle. We can then derive the necessary condition for 𝒢\mathcal{G} based on the decomposed cycle set 𝒳i\mathcal{X}_{i} with respect to the edge eie_{i}\in\mathcal{E} and the remaining edges in 𝒫i\mathcal{P}_{i}. This is done based on analysing the convergence of each edge, which is affected by its neighboring edges. Thus the insights here of the complete decomposition are: for the first requirement of Definition 6, we avoid considering the neighboring edges repeatedly; for the second requirement of Definition 6, we consider the smallest cycle that contains the edge in hand. The following theorem proposes a necessary condition on general graphs with cycles.

Theorem 8.

Consider the leader-follower multi-agent system Σ\Sigma described by the graph 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}). A necessary condition on 𝒢\mathcal{G} under which we can design the leaders using (6) to achieve the target formation \mathcal{F} as in (1) while satisfying (5) is that every leaderless graph 𝒢f=(𝒱f,f)\mathcal{G}^{f}=(\mathcal{V}^{f},\mathcal{E}^{f}), such that 𝒢f𝒢\mathcal{G}^{f}\subseteq\mathcal{G}, should satisfy

𝒞ei𝒳i{min(|(𝒞ei)|4,2)}+|Ei|2,eif.\sum\limits_{\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}}\left\{\min(|\mathcal{E}(\mathcal{C}_{e_{i}})|-4,2)\right\}+|E_{i}|\leq 2,\hskip 8.53581pt\forall e_{i}\in\mathcal{E}^{f}. (8)

where Ei={ekek𝒩(ei)𝒫i}E_{i}=\{e_{k}\mid e_{k}\in\mathcal{N}(e_{i})\cap\mathcal{P}_{i}\}, and 𝒳i\mathcal{X}_{i} is the cycle set obtained via the complete decomposition of 𝒢f\mathcal{G}^{f} with respect to eie_{i} as in Definition 6.

{pf}

The proof is based on the discussion of the decomposed cycles 𝒞ei𝒳i\mathcal{C}_{e_{i}}\in\mathcal{X}_{i} and the remaining edges in 𝒫i\mathcal{P}_{i}. We can resort to Lemma 5 to deal with 𝒫i\mathcal{P}_{i}. Here, we first discuss the result for one cycle, e.g., 𝒞e1a\mathcal{C}_{e_{1}}^{a} with respect to the edge e1e_{1} which has mm edges (𝒞e1a)={e1,e2,,em}\mathcal{E}(\mathcal{C}_{e_{1}}^{a})=\{e_{1},e_{2},\dots,e_{m}\} and m3m\geq 3 (need at least 3 edges to form a cycle). We will check how the number of edges of the cycle affects convergence. When m=3m=3, the state evolution of an arbitrary edge e1e_{1} is x¯˙1=2x¯1+x¯2+x¯3\dot{\bar{x}}_{1}=-2\bar{x}_{1}+\bar{x}_{2}+\bar{x}_{3}, and since ei,i=1,2,3e_{i},i=1,2,3 form a cycle, we have that x¯1+x¯2+x¯3=0\bar{x}_{1}+\bar{x}_{2}+\bar{x}_{3}=0. Hence x¯2+x¯3=x¯1\bar{x}_{2}+\bar{x}_{3}=-\bar{x}_{1} and x¯˙1=2x¯1x¯1=3x¯1\dot{\bar{x}}_{1}=-2\bar{x}_{1}-\bar{x}_{1}=-3\bar{x}_{1}. This means that the cycle that forms a triangle will show a higher convergence rate of 3-3 for the edge dynamics. When m=4m=4, the state evolution of an arbitrary edge e1e_{1} is x¯˙1=2x¯1+x¯2+x¯4\dot{\bar{x}}_{1}=-2\bar{x}_{1}+\bar{x}_{2}+\bar{x}_{4}, and since ei,i=1,2,3,4e_{i},i=1,2,3,4 form a cycle, we have that i=14x¯i=0\sum_{i=1}^{4}\bar{x}_{i}=0. Hence x¯2+x¯4=x¯1x¯3\bar{x}_{2}+\bar{x}_{4}=-\bar{x}_{1}-\bar{x}_{3} and x¯˙1=2x¯1x¯1x¯3\dot{\bar{x}}_{1}=-2\bar{x}_{1}-\bar{x}_{1}-\bar{x}_{3}. Consider the worst case, when x¯1\bar{x}_{1} is arbitrarily close to the prescribed performance bound and x¯3=x¯1\bar{x}_{3}=-\bar{x}_{1}, we then have x¯˙1=2x¯1x¯1+x¯1=2x¯1\dot{\bar{x}}_{1}=-2\bar{x}_{1}-\bar{x}_{1}+\bar{x}_{1}=-2\bar{x}_{1}. This means that the cycle still shows decay rate of 2-2 for the edge dynamics. When m=5m=5, the state evolution of an arbitrary edge e1e_{1} is x¯˙1=2x¯1+x¯2+x¯5\dot{\bar{x}}_{1}=-2\bar{x}_{1}+\bar{x}_{2}+\bar{x}_{5}, and since ei,i=1,2,3,4,5e_{i},i=1,2,3,4,5 form a cycle, we have that i=15x¯i=0\sum_{i=1}^{5}\bar{x}_{i}=0. Hence x¯2+x¯5=x¯1x¯3x¯4\bar{x}_{2}+\bar{x}_{5}=-\bar{x}_{1}-\bar{x}_{3}-\bar{x}_{4} and x¯˙1=2x¯1x¯1x¯3x¯4\dot{\bar{x}}_{1}=-2\bar{x}_{1}-\bar{x}_{1}-\bar{x}_{3}-\bar{x}_{4}. Consider the worst case, when x¯1\bar{x}_{1} is arbitrarily close to the performance bound and x¯3=x¯4=x¯1\bar{x}_{3}=\bar{x}_{4}=-\bar{x}_{1}, we then have x¯˙1=2x¯1x¯1+2x¯1=x¯1\dot{\bar{x}}_{1}=-2\bar{x}_{1}-\bar{x}_{1}+2\bar{x}_{1}=-\bar{x}_{1}. This means that the cycle with 55 edges still shows decay rate of 1-1 for the edge dynamics. This means that we still have the freedom to add one more edge that shares a node with e1e_{1}. For m6m\geq 6, the state evolution of an arbitrary edge e1e_{1} is x¯˙1=2x¯1+x¯2+x¯m\dot{\bar{x}}_{1}=-2\bar{x}_{1}+\bar{x}_{2}+\bar{x}_{m}, and since ei,i=1,2,,me_{i},i=1,2,\dots,m forms a cycle, we have that i=1mx¯i=0\sum_{i=1}^{m}\bar{x}_{i}=0. Then x¯2+x¯m=x¯1i=3m1x¯i\bar{x}_{2}+\bar{x}_{m}=-\bar{x}_{1}-\sum_{i=3}^{m-1}\bar{x}_{i} and x¯˙1=2x¯1x¯1i=3m1x¯i\dot{\bar{x}}_{1}=-2\bar{x}_{1}-\bar{x}_{1}-\sum_{i=3}^{m-1}\bar{x}_{i}. i=3m1x¯i-\sum_{i=3}^{m-1}\bar{x}_{i} cannot be greater than 3x¯13\bar{x}_{1} when x¯1\bar{x}_{1} is arbitrarily close to the performance bound since ei,i=1,2,,me_{i},i=1,2,\dots,m need to form a cycle satisfying i=1mx¯i=0\sum_{i=1}^{m}\bar{x}_{i}=0. This means that for the cycle with m6m\geq 6, in the worst case we have x¯˙1=0\dot{\bar{x}}_{1}=0, thus in the worst case x¯1\bar{x}_{1} will never evolve again to violate the performance bound. Till here, we can summarize that for a single cycle 𝒞e1a\mathcal{C}_{e_{1}}^{a}, in the worst case the decay rate of the edge dynamics is 2+min(|(𝒞e1a)|4,2)-2+\min(|\mathcal{E}(\mathcal{C}_{e_{1}}^{a})|-4,2). Then, based on the result for a single cycle, we build the result in the case that eie_{i} belongs to more than one cycles, i.e., eie_{i} is an edge of 𝒞ei𝒳i\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}. Then each cycle of 𝒞ei𝒳i\mathcal{C}_{e_{i}}\in\mathcal{X}_{i} will contribute a decay rate of min(|(𝒞ei)|4,2)\min(|\mathcal{E}(\mathcal{C}_{e_{i}})|-4,2) for eie_{i}. Since the cycles 𝒞ei𝒳i\mathcal{C}_{e_{i}}\in\mathcal{X}_{i} will not affect each other in contributing to the decay rate of eie_{i} since they are completely decomposed with respect to eie_{i}, then the total decay rate of eie_{i} that belongs to more than one cycle is 2+𝒞ei𝒳i{min(|(𝒞ei)|4,2)}-2+\sum\limits_{\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}}\left\{\min(|\mathcal{E}(\mathcal{C}_{e_{i}})|-4,2)\right\}, where 2-2 corresponds to the diagonal entry of the edge Laplacian Le-L_{e}. Finally, the remaining edges that affect the convergence of eie_{i} are the edges that share a node with eie_{i} but are not an edge of any 𝒞ei𝒳i\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}, i.e., Ei={ekek𝒩(ei)𝒫i}E_{i}=\{e_{k}\mid e_{k}\in\mathcal{N}(e_{i})\cap\mathcal{P}_{i}\}. In the worst case, each of these edges will contribute a decay rate of 11 to edge eie_{i}. Hence, the total decay rate of eie_{i} in the worst case is 2+𝒞ei𝒳i{min(|(𝒞ei)|4,2)}+|Ei|-2+\sum\limits_{\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}}\left\{\min(|\mathcal{E}(\mathcal{C}_{e_{i}})|-4,2)\right\}+|E_{i}| and should satisfy that 2+𝒞ei𝒳i{min(|(𝒞ei)|4,2)}+|Ei|0-2+\sum\limits_{\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}}\left\{\min(|\mathcal{E}(\mathcal{C}_{e_{i}})|-4,2)\right\}+|E_{i}|\leq 0, i.e., exactly the inequality in (8). This means that in the worst case when the term 𝒞ei𝒳i{min(|(𝒞ei)|4,2)}+|Ei|\sum\limits_{\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}}\left\{\min(|\mathcal{E}(\mathcal{C}_{e_{i}})|-4,2)\right\}+|E_{i}| starts to be greater than 2, eie_{i} will continue evolving to violate the prescribed performance bounds. Therefore, suppose that there exists eife_{i}\in\mathcal{E}^{f} satisfying 𝒞ei𝒳i{min(|(𝒞ei)|4,2)}+|Ei|3\sum\limits_{\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}}\left\{\min(|\mathcal{E}(\mathcal{C}_{e_{i}})|-4,2)\right\}+|E_{i}|\geq 3, then no matter how we design the leaders in 𝒢\mathcal{G}, in the worst case eie_{i} will always evolve to violate the prescribed performance bound, which leads to a contradiction. Hence, we can conclude that 𝒢\mathcal{G} should not contain a leaderless induced subgraph 𝒢f=(𝒱f,f)\mathcal{G}^{f}=(\mathcal{V}^{f},\mathcal{E}^{f}) such that there exists eife_{i}\in\mathcal{E}^{f} satisfying 𝒞ei𝒳i{min(|(𝒞ei)|4,2)}+|Ei|3\sum\limits_{\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}}\left\{\min(|\mathcal{E}(\mathcal{C}_{e_{i}})|-4,2)\right\}+|E_{i}|\geq 3. Or in other words, every leaderless graph 𝒢f=(𝒱f,f)\mathcal{G}^{f}=(\mathcal{V}^{f},\mathcal{E}^{f}) such that 𝒢f𝒢\mathcal{G}^{f}\subseteq\mathcal{G}, should satisfy (8).∎

Remark 9.

Lemma 5 and Theorem 8 indicate that when considering the leader-follower multi-agent system Σ\Sigma to achieve the target formation within the prescribed performance bounds, we should assign leaders such that (8) holds. They also propose a criterion in choosing leaders to achieve the target formation with prescribed performance guarantees, which can be further applied to solve leader selection problems. Note that Lemma 5 is a specific case of Theorem 8. That is, when the graph is a tree, then the term 𝒞ei𝒳i{min(|(𝒞ei)|4,2)}\sum\limits_{\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}}\left\{\min(|\mathcal{E}(\mathcal{C}_{e_{i}})|-4,2)\right\} vanishes and |Ei||E_{i}| is the number of edges that share a node with eie_{i}, i.e., |Ei|=|𝒩(ei)||E_{i}|=|\mathcal{N}(e_{i})|. This also leads to the condition |𝒩(ei)|2|\mathcal{N}(e_{i})|\leq 2 for tree graphs as in (7). Condition (7) of Lemma 5 indicates that for a leaderless tree graph, each edge can have at most 22 neighboring edges, which also means that each agent can have at most 33 neighboring agents. Furthermore, according to Theorem 8, for a general leaderless graph, agents that are belonging to some cycles can have more neighboring agents when including more cycles as long as condition (8) holds.

4 Necessary and Sufficient Conditions on Graph Topology

In the previous section, we derived necessary conditions on the graph topology that need to hold for both tree graphs and general graphs with cycles. We next derive sufficient conditions on the graph topology for both tree graphs and general graphs with cycles in this section. First of all, note that the conditions in Lemma 5 and Theorem 8 are necessary but not sufficient conditions since we only check the induced subgraphs (as in Definition 1) that are composed of followers. However, we have not taken the choice of leader vertices into account, which show the couplings between different induced subgraphs with only followers. In fact, for a general graph satisfying the necessary conditions in Lemma 5 or Theorem 8, we still may not find a suitable PPC law for the leaders to steer the leader-follower multi-agent system to achieve the target formation within the prescribed performance bounds because the conditions in Lemma 5 or Theorem 8 are not sufficient when taking the leaders into account. In the sequel, we take the choice of leader vertices into account and focus on finding a sufficient condition on the graph topology, which will show to head to a necessary and sufficient condition. Consider the leader-follower multi-agent system under the communication graph 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}), we first define the following follower-leader-follower (FLF) path and maximum follower-end subgraph as

Definition 10.

(FLF path.) A path subgraph p=(𝒱p,p)Pp=(\mathcal{V}^{p},\mathcal{E}^{p})\in P of 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}) as in Definition 3 with vi,vj𝒱pv_{i},v_{j}\in\mathcal{V}^{p} as the two end nodes, is a FLF path of 𝒢\mathcal{G} if 𝒱p={vi,vj}\mathcal{V}_{\mathscr{F}}^{p}=\{v_{i},v_{j}\}. PP is the collection of all the FLF paths of 𝒢\mathcal{G}. The neighborhood of the FLF path pPp\in P denoted as 𝔑(p)\mathfrak{N}(p) is the set of edges that share the node viv_{i} or vjv_{j} with pp, i.e., 𝔑(p)={ekvi𝐯(ek) or vj𝐯(ek),ekp}\mathfrak{N}(p)=\{e_{k}\mid v_{i}\in\mathbf{v}(e_{k})\text{ or }v_{j}\in\mathbf{v}(e_{k}),e_{k}\in\mathcal{E}\setminus\mathcal{E}^{p}\}.

Definition 11.

(Maximum follower-end subgraph.) A graph 𝒢=(𝒱,)\mathcal{G}^{\star}=(\mathcal{V}^{\star},\mathcal{E}^{\star}) is a maximum follower-end subgraph of the graph 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}), denoted as 𝒢𝒢\mathcal{G}^{\star}\preceq\mathcal{G} if the following conditions hold:

  • 𝒢𝒢\mathcal{G}^{\star}\subseteq\mathcal{G} (as in Definition 1);

  • every i𝒱i\in\mathcal{V}_{\mathscr{L}}^{\star} belongs to a FLF path of 𝒢\mathcal{G}^{\star};

  • there is no subgraph 𝒢=(𝒱,)\mathcal{G}^{\prime}=(\mathcal{V}^{\prime},\mathcal{E}^{\prime}) of 𝒢\mathcal{G} such that every i𝒱i\in\mathcal{V}_{\mathscr{L}}^{\prime} belongs to a FLF path satisfying |𝒱|>|𝒱||\mathcal{V}^{\prime}|>|\mathcal{V}^{\star}|.

Definition 11 is proposed in order to derive an induced subgraph 𝒢\mathcal{G}^{\star} of 𝒢\mathcal{G} (first requirement) such that the end leaders are ignored (second requirement) and includes as many agents as possible (third requirement). The insight here is that when a leader is placed in the graph as an end node, then this leader does not show any couplings between agents since it only connects to one agent and can freely move. For the maximum follower-end subgragh 𝒢\mathcal{G}^{\star}, we can analyse the convergence of eie_{i}\in\mathcal{E}_{\mathscr{F}\mathscr{F}}^{\star} by resorting to the ideas and arguments from Lemma 5 and Theorem 8. Note that we need to further consider the convergence of the dynamics of the leader-follower and leader-leader edge eie_{i}\in\mathcal{E}_{\mathscr{L}\mathscr{F}}^{\star}\cup\mathcal{E}_{\mathscr{L}\mathscr{L}}^{\star}, and the collection of all the FLF paths PP will traverse all the edges eie_{i}\in\mathcal{E}_{\mathscr{L}\mathscr{F}}^{\star}\cup\mathcal{E}_{\mathscr{L}\mathscr{L}}^{\star}, thus we next discuss the convergence results for FLF paths. When considering the maximum follower-end subgragh 𝒢\mathcal{G}^{\star}, we can also completely decompose 𝒢\mathcal{G}^{\star} with respect to a specific FLF path piPp_{i}\in P. The complete decomposition of a graph 𝒢\mathcal{G} with respect to a specific FLF path pip_{i} is defined in a similar manner to Definition 6, and is given as follows:

Definition 12.

(Complete decomposition w.r.t. FLF path.) A graph 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}) is decomposed with respect to the FLF path piPp_{i}\in P as 𝒢:=𝒞pi𝒴i𝒞pi𝒬i\mathcal{G}:=\cup_{\mathcal{C}_{p_{i}}\in\mathcal{Y}_{i}}\mathcal{C}_{p_{i}}\cup\mathcal{Q}_{i}, where 𝒴i:={𝒞pipi𝒞pi}\mathcal{Y}_{i}:=\{\mathcal{C}_{p_{i}}\mid p_{i}\subseteq\mathcal{C}_{p_{i}}\} is the cycle set composed of all the cycles 𝒞pi\mathcal{C}_{p_{i}} in 𝒢\mathcal{G} that contain pip_{i} as per Definition 3 and satisfy:

  • for every pair of 𝒞pia,𝒞pib𝒴i\mathcal{C}_{p_{i}}^{a},\mathcal{C}_{p_{i}}^{b}\in\mathcal{Y}_{i}, (𝔑(pi)𝒞pia)(𝔑(pi)𝒞pib)=(\mathfrak{N}(p_{i})\cap\mathcal{C}_{p_{i}}^{a})\cap(\mathfrak{N}(p_{i})\cap\mathcal{C}_{p_{i}}^{b})=\emptyset.

  • for every 𝒞pi𝒴i\mathcal{C}_{p_{i}}\in\mathcal{Y}_{i}, there does not exist a cycle 𝒞\mathcal{C} of 𝒢\mathcal{G} such that pi𝒞,(𝔑(pi)𝒞pi)(𝔑(pi)𝒞)p_{i}\subseteq\mathcal{C},(\mathfrak{N}(p_{i})\cap\mathcal{C}_{p_{i}})\cap(\mathfrak{N}(p_{i})\cap\mathcal{C})\neq\emptyset, and |(𝒞)|<|(𝒞pi)||\mathcal{E}(\mathcal{C})|<|\mathcal{E}(\mathcal{C}_{p_{i}})|,

and where 𝒬i:={ekek𝒞pi,𝒞pi𝒴i}\mathcal{Q}_{i}:=\{e_{k}\mid e_{k}\notin\mathcal{C}_{p_{i}},\mathcal{C}_{p_{i}}\in\mathcal{Y}_{i}\} is the set of the edges that do not belong to any cycle in 𝒴i\mathcal{Y}_{i}. Then, we call this decomposition a complete decomposition of 𝒢\mathcal{G} with respect to the FLF path piPp_{i}\in P.

Now, we derive the following theorem regarding the necessary and sufficient conditions on the graph topology of the leader-follower multi-agent system in order to achieve the target formation \mathcal{F} as in (1) under the control (6) while satisfying the prescribed performance bounds (5).

Theorem 13.

Consider the leader-follower multi-agent system Σ\Sigma described by the graph 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}), and let 𝒢𝒢\mathcal{G}^{\star}\preceq\mathcal{G}. A necessary and sufficient condition on 𝒢\mathcal{G} under which we can design the leaders using (6) to achieve the target formation \mathcal{F} as in (1) while satisfying (5) is that:

  • for any eie_{i}\in\mathcal{E}^{\star}_{\mathscr{F}\mathscr{F}},

    𝒞ei𝒳i{min(|(𝒞ei)|4,2)}+|Ei|2;\sum\limits_{\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}}\left\{\min(|\mathcal{E}(\mathcal{C}_{e_{i}})|-4,2)\right\}+|E_{i}|\leq 2; (9)
  • for any piPp_{i}\in P^{\star}, there either exists 𝒞pi𝒴i\mathcal{C}_{p_{i}}\in\mathcal{Y}_{i} such that |(𝒞pi)|<2|(pi)||\mathcal{E}(\mathcal{C}_{p_{i}})|<2|\mathcal{E}(p_{i})|, or

    𝒞pi𝒴i{min(|(𝒞pi)|2|(pi)|2,2)}+|Fi|2,\sum\limits_{\mathcal{C}_{p_{i}}\in\mathcal{Y}_{i}}\left\{\min(|\mathcal{E}(\mathcal{C}_{p_{i}})|-2|\mathcal{E}(p_{i})|-2,2)\right\}+|F_{i}|\leq 2, (10)

where Ei={ekek𝒩(ei)𝒫i}E_{i}=\{e_{k}\mid e_{k}\in\mathcal{N}(e_{i})\cap\mathcal{P}_{i}\}, 𝒳i\mathcal{X}_{i} is the cycle set obtained via the complete decomposition of 𝒢\mathcal{G}^{\star} with respect to eie_{i}\in\mathcal{E}^{\star}_{\mathscr{F}\mathscr{F}} as in Definition 6, PP^{\star} is the collection of all FLF paths of 𝒢\mathcal{G}^{\star}, Fi={ekek𝔑(pi)𝒬i}F_{i}=\{e_{k}\mid e_{k}\in\mathfrak{N}(p_{i})\cap\mathcal{Q}_{i}\}, and 𝒴i\mathcal{Y}_{i} is the cycle set obtained via the complete decomposition of 𝒢\mathcal{G}^{\star} with respect to piPp_{i}\in P^{\star} as in Definition 12.

{pf}

We first prove the sufficient part. That is, assuming that the conditions (9) or (10) hold, we show that we can design the leaders using (6) to achieve the target formation \mathcal{F} as in (1) while satisfying (5). Firstly, if 𝒞ei𝒳i{min(|(𝒞ei)|4,2)}+|Ei|2\sum\limits_{\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}}\left\{\min(|\mathcal{E}(\mathcal{C}_{e_{i}})|-4,2)\right\}+|E_{i}|\leq 2 holds for every edge eie_{i}\in\mathcal{E}^{\star}_{\mathscr{F}\mathscr{F}}, then in the worst case, when x¯i\bar{x}_{i} is arbitrarily close to the prescribed performance bound, it will not evolve to violate the performance bound. Therefore, there always exists a positive, smooth and strictly decreasing performance function ρx¯i(t)\rho_{\bar{x}_{i}}(t) that can be predefined such that the dynamics of the edge eie_{i}\in\mathcal{E}^{\star}_{\mathscr{F}\mathscr{F}} will evolve while satisfying (5). Next, suppose that (10) holds, we first denote an arbitrary FLF path piPp_{i}\in P^{\star} by the consecutive sequence of nodes as v1,,vnpv_{1},\cdots,v_{n_{p}} or the consecutive sequence of edges as e1,,enp1e_{1},\cdots,e_{n_{p}-1} (a path with npn_{p} nodes and np1n_{p}-1 edges) with v1,vnpv_{1},v_{n_{p}} as the two followers and e1,enp1e_{1},e_{n_{p}-1} as the two leader-follower edges of path piPp_{i}\in P^{\star}. If there exists a cycle 𝒞pi1𝒴i\mathcal{C}_{p_{i}}^{1}\in\mathcal{Y}_{i} that contains the FLF path pip_{i} as a subgraph and satisfies |(𝒞pi)|<2|(pi)||\mathcal{E}(\mathcal{C}_{p_{i}})|<2|\mathcal{E}(p_{i})|, then we know that the nodes v1,,vnpv_{1},\cdots,v_{n_{p}} and the edges e1,,enp1e_{1},\cdots,e_{n_{p}-1} of pip_{i} are also the respective nodes and edges of 𝒞pi1\mathcal{C}_{p_{i}}^{1}. Without loss of generality, suppose that the nodes of the cycle 𝒞pi1\mathcal{C}_{p_{i}}^{1} are v1,,vnp,vnp+1,,vncv_{1},\cdots,v_{n_{p}},v_{n_{p}+1},\cdots,v_{n_{c}} and the edges of the cycle 𝒞pi1\mathcal{C}_{p_{i}}^{1} are e1,,enp1,enp,,ence_{1},\cdots,e_{n_{p}-1},e_{n_{p}},\cdots,e_{n_{c}} (a cycle with ncn_{c} nodes and ncn_{c} edges), then we denote the remaining edges of 𝒞pi1\mathcal{C}_{p_{i}}^{1} besides e1,,enp1e_{1},\cdots,e_{n_{p}-1} as enp,,ence_{n_{p}},\cdots,e_{n_{c}}. Since e1,,enp1e_{1},\cdots,e_{n_{p}-1} and enp,,ence_{n_{p}},\cdots,e_{n_{c}} form the cycle 𝒞pi1\mathcal{C}_{p_{i}}^{1}, we have that i=1np1x¯i+j=npncx¯j=0\sum_{i=1}^{n_{p}-1}\bar{x}_{i}+\sum_{j=n_{p}}^{n_{c}}\bar{x}_{j}=0 or i=1np1x¯i=j=npncx¯j\sum_{i=1}^{n_{p}-1}\bar{x}_{i}=-\sum_{j={n_{p}}}^{n_{c}}\bar{x}_{j}. Since we have |(𝒞pi)|<2|(pi)||\mathcal{E}(\mathcal{C}_{p_{i}})|<2|\mathcal{E}(p_{i})|, i.e., nc<2(np1)n_{c}<2(n_{p}-1), then the number of edges of e1enp1e_{1}\cdots e_{n_{p}-1} (np1n_{p}-1 edges) are more than that of the remaining edges enpence_{n_{p}}\cdots e_{n_{c}} (ncnp+1n_{c}-n_{p}+1 edges). Therefore, i=1np1x¯i\sum_{i=1}^{n_{p}-1}\bar{x}_{i} will never evolve to violate the performance bound (the performance bound for i=1np1x¯i\sum_{i=1}^{n_{p}-1}\bar{x}_{i} is the summation of the performance bounds of each x¯i\bar{x}_{i}) since it is constrained by enpence_{n_{p}}\cdots e_{n_{c}} to form the cycle 𝒞pi1\mathcal{C}_{p_{i}}^{1}. Otherwise, if all 𝒞pi𝒴i\mathcal{C}_{p_{i}}\in\mathcal{Y}_{i} satisfy |(𝒞pi)|2|(pi)||\mathcal{E}(\mathcal{C}_{p_{i}})|\geq 2|\mathcal{E}(p_{i})|, we can further calculate the dynamics of pip_{i} as i=1np1x¯˙i=(x¯1+x¯np1)+{jej𝔑(pi)}x¯j\sum_{i=1}^{n_{p}-1}\dot{\bar{x}}_{i}=-(\bar{x}_{1}+\bar{x}_{n_{p}-1})+\sum\nolimits_{\{j\mid e_{j}\in\mathfrak{N}(p_{i})\}}\bar{x}_{j}. Similar to the arguments in Theorem 8 when discussing the cycles, the first term of (x¯1+x¯np1)-(\bar{x}_{1}+\bar{x}_{n_{p}-1}) contributes the decay rate of 2-2 in the dynamics of pip_{i}. The second term corresponds to the neighboring edges of pip_{i}, which include the edges belonging to some cycles that contain pip_{i} and the remaining edges defined by FiF_{i}. We first consider the case when there exists only one cycle 𝒞pi1\mathcal{C}_{p_{i}}^{1} that contains pip_{i} as a subgraph, e1enp1e_{1}\cdots e_{n_{p}-1} are the edges of pip_{i}, 𝔑(pi)𝒞pi1={enp,enc}\mathfrak{N}(p_{i})\cap\mathcal{C}_{p_{i}}^{1}=\{e_{n_{p}},e_{n_{c}}\} are the two neighboring edges of pip_{i} in the cycle 𝒞pi1\mathcal{C}_{p_{i}}^{1}, and the remaining edges are enp+1,,enc1e_{n_{p}+1},\cdots,e_{n_{c}-1}. Then the dynamics of pip_{i} are rewritten as i=1np1x¯˙i=(x¯1+x¯np1)+x¯np+x¯nc+{jekFi}x¯k\sum_{i=1}^{n_{p}-1}\dot{\bar{x}}_{i}=-(\bar{x}_{1}+\bar{x}_{n_{p}-1})+\bar{x}_{n_{p}}+\bar{x}_{n_{c}}+\sum\nolimits_{\{j\mid e_{k}\in F_{i}\}}\bar{x}_{k}. We know that i=1np1x¯i+x¯np+x¯nc+j=np+1nc1x¯j=0\sum_{i=1}^{n_{p}-1}\bar{x}_{i}+\bar{x}_{n_{p}}+\bar{x}_{n_{c}}+\sum_{j=n_{p}+1}^{n_{c}-1}\bar{x}_{j}=0 as they form the cycle 𝒞pi1\mathcal{C}_{p_{i}}^{1}, thus x¯np+x¯nc=i=1np1x¯ij=np+1nc1x¯j\bar{x}_{n_{p}}+\bar{x}_{n_{c}}=-\sum_{i=1}^{n_{p}-1}\bar{x}_{i}-\sum_{j=n_{p}+1}^{n_{c}-1}\bar{x}_{j} and the dynamics of pip_{i} are i=1np1x¯˙i=(x¯1+x¯np1)i=1np1x¯ij=np+1nc1x¯j+{jekFi}x¯k\sum_{i=1}^{n_{p}-1}\dot{\bar{x}}_{i}=-(\bar{x}_{1}+\bar{x}_{n_{p}-1})-\sum_{i=1}^{n_{p}-1}\bar{x}_{i}-\sum_{j=n_{p}+1}^{n_{c}-1}\bar{x}_{j}+\sum\nolimits_{\{j\mid e_{k}\in F_{i}\}}\bar{x}_{k}. We have that the number of edges of e1enp1e_{1}\cdots e_{n_{p}-1} is n1=|(pi)|=np1n_{1}=|\mathcal{E}(p_{i})|=n_{p}-1 and the number of edges of enp+1enc1e_{n_{p}+1}\cdots e_{n_{c}-1} is n2=|(𝒞pi1)||(pi)|2=ncnp1n_{2}=|\mathcal{E}(\mathcal{C}_{p_{i}}^{1})|-|\mathcal{E}(p_{i})|-2=n_{c}-n_{p}-1. When n2n12n_{2}-n_{1}\leq 2, the term i=1np1x¯ij=np+1nc1x¯j-\sum_{i=1}^{n_{p}-1}\bar{x}_{i}-\sum_{j=n_{p}+1}^{n_{c}-1}\bar{x}_{j} contributes a decay rate of n2n1n_{2}-n_{1} in the dynamics of pip_{i}, while n2n1>2n_{2}-n_{1}>2, the term i=1np1x¯ij=np+1nc1x¯j-\sum_{i=1}^{n_{p}-1}\bar{x}_{i}-\sum_{j=n_{p}+1}^{n_{c}-1}\bar{x}_{j} contributes a decay rate of 22 in the dynamics of pip_{i} due to the constraints to form the cycle 𝒞pi1\mathcal{C}_{p_{i}}^{1}. In addition the term {jekFi}x¯k\sum\nolimits_{\{j\mid e_{k}\in F_{i}\}}\bar{x}_{k} contributes a decay rate of |Fi||F_{i}| in the dynamics of pip_{i}. Hence, in the worst case the decay rate of the dynamics of pip_{i} is i=1np1x¯˙i=2+min(n2n1,2)+|Fi|=2+min(|(𝒞pi1)|2|(pi)|2,2)+|Fi|0\sum_{i=1}^{n_{p}-1}\dot{\bar{x}}_{i}=-2+\min(n_{2}-n_{1},2)+|F_{i}|=-2+\min(|\mathcal{E}(\mathcal{C}_{p_{i}}^{1})|-2|\mathcal{E}(p_{i})|-2,2)+|F_{i}|\leq 0 according to condition (10). Next, we extend the analysis to the cases when there are more than one cycle that contain pip_{i} as a subgraph. Similarly then the decay rate of the dynamics of pip_{i} is i=1np1x¯˙i=2+𝒞pi𝒴i{min(|(𝒞pi)|2|(pi)|2,2)}+|Fi|0\sum_{i=1}^{n_{p}-1}\dot{\bar{x}}_{i}=-2+\sum\limits_{\mathcal{C}_{p_{i}}\in\mathcal{Y}_{i}}\left\{\min(|\mathcal{E}(\mathcal{C}_{p_{i}})|-2|\mathcal{E}(p_{i})|-2,2)\right\}+|F_{i}|\leq 0 according to condition (10). This means that in the worst case, when condition (10) holds, we can always ensure that when each edge of piPp_{i}\in P^{\star} is arbitrarily close to the prescribed performance bound, it will not evolve to violate the performance bound. Then, for the path pip_{i} which only has two followers as the end nodes and the remaining nodes are all leaders, we can apply the control law (6) for the leaders in pip_{i} to ensure that every edge in pip_{i} satisfies the prescribed performance bound (5). This is generally due to the fact that each edge of a path graph has at most 2 neighboring edges. We also refer the readers to (Chen & Dimarogonas, 2020b) for the convergence analysis of the specific path graphs. Finally, since PP^{\star} is the collection of all FLF paths, it includes all the edges eie_{i}\in\mathcal{E}^{\star}_{\mathscr{L}\mathscr{L}} and all the edges eie_{i}\in\mathcal{E}^{\star}_{\mathscr{L}\mathscr{F}}. Note that since condition (9) is for all the edges eie_{i}\in\mathcal{E}^{\star}_{\mathscr{F}\mathscr{F}}, while condition (10) considers all the edges eie_{i}\in\mathcal{E}^{\star}_{\mathscr{L}\mathscr{L}}\cup\mathcal{E}^{\star}_{\mathscr{L}\mathscr{F}}, these two conditions are independent of each other. Till now, we have considered all the edges ei=e_{i}\in\mathcal{E}^{\star}=\mathcal{E}^{\star}_{\mathscr{F}\mathscr{F}}\cup\mathcal{E}^{\star}_{\mathscr{L}\mathscr{F}}\cup\mathcal{E}^{\star}_{\mathscr{L}\mathscr{L}} and the dynamics of all the edges satisfy the prescribed performance bound (5). Therefore, these two conditions are sufficient such that the leader-follower multi-agent system can achieve the target formation \mathcal{F} as in (1) while satisfying the prescribed performance bound (5).

For the necessity part, we use contradiction. That is, suppose that either condition (9) or condition (10) does not hold. Firstly, suppose that there exists an edge eiFFe_{i}\in\mathcal{E}^{\star}_{FF} such that 𝒞ei𝒳i{min(|(𝒞ei)|4,2)}+|Ei|3\sum\limits_{\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}}\left\{\min(|\mathcal{E}(\mathcal{C}_{e_{i}})|-4,2)\right\}+|E_{i}|\geq 3. Then, if eie_{i} and its neighboring edges are arbitrarily close to the prescribed performance bound, eie_{i} will continue to evolve to violate the performance bound. This results in a contradiction. Similarly, suppose that there exists a FLF path piPp_{i}\in P^{\star} such that condition (10) does not hold. Then, if the edges of pip_{i} and the neighboring edges of pip_{i} are arbitrarily close to the prescribed performance bound, pip_{i} will continue to evolve to violate the performance bound, which also leads to a contradiction. Note that the performance bound for pip_{i} is the summation of the performance bounds of each edge of pip_{i}. Finally, we can conclude that in order to achieve the target formation \mathcal{F} as in (1) while satisfying (5) by applying the control law (6), conditions (9) and (10) should hold. Therefore, conditions (9) and (10) are also necessary.∎

Remark 14.

Theorem 13 presents a necessary and sufficient condition on the leader-follower graph topology such that we can further design the leaders in order to steer the entire system to achieve the target formation within the transient constraints bounds. This also proposes a methodology for leader selection problems. Conditions (9) and (10) indicate the trade-offs among the leader-follower topology, number of cycles, number and positions of the leaders. Note that if we generalize the follower-follower edge also as a FLF path, i.e., a path that starts with a follower, ends with a follower, but crosses no leaders, then, the first condition in Theorem 13 can be regarded as a special case of the second condition, since |(𝒞ei)|<2|(ei)|=2|\mathcal{E}(\mathcal{C}_{e_{i}})|<2|\mathcal{E}(e_{i})|=2 will never hold for eiFFe_{i}\in\mathcal{E}^{\star}_{FF}, and min(|(𝒞ei)|2|(ei)|2,2)\min(|\mathcal{E}(\mathcal{C}_{e_{i}})|-2|\mathcal{E}(e_{i})|-2,2) is exactly min(|(𝒞ei)|4,2)\min(|\mathcal{E}(\mathcal{C}_{e_{i}})|-4,2). Similar to the discussions on Lemma 5 and Theorem 8 in the previous section, Theorem 13 holds for both tree graphs and general graphs with cycles.

5 Simulations and Examples

In this section, simulation examples are presented in order to verify the results. We first consider a multi-vehicle platooning example, and later on a multi-robot coordination example with different choice of leader robots. The simulations’ communication graphs are shown in Fig. 3, where the leaders and followers are represented by grey and white nodes, respectively.

Refer to caption
Figure 3: Leader-follower communication graphs.

Multi-vehicle platooning: we consider a vehicle platoon with 99 vehicles as shown in Fig. 3.CC with leader vehicles 𝒱={v3,v4,v7,v9}.\mathcal{V}_{\mathscr{L}}=\{v_{3},v_{4},v_{7},v_{9}\}. The vehicles are initialized with p=[0,20,60,105,125,145,185,205,250]Tp=[0,20,60,105,125,145,185,205,250]^{T}. All edges are subject to the transient constraints as in (5) with ρx¯i(t)=20et+0.1\rho_{\bar{x}_{i}}(t)=20e^{-t}+0.1 and p¯ides=30\bar{p}^{des}_{i}=30. The transient constraints under consideration encapsulate two primary tasks: ensuring collision avoidance with a minimum distance of 10 and maintaining connectivity within a maximum distance of 50. Furthermore, these constraints prescribe an additional requirement on the transient behavior, particularly concerning the convergence rate and overshoot towards the target platoon formation. Therefore, the vehicle platoon aims to achieve the formation (keeping a safe distance of 3030 between neighboring vehicles) while satisfying the time-varying constraints by only designing the controllers of the leader vehicles. We can check that the necessary and sufficient conditions proposed in Theorem 13 are trivially satisfied for multi-vehicle platooning since each vehicle has at most 22 neighboring vehicles. The simulation results when applying the PPC strategy (6) are depicted in Fig. 4. The top figure plots the platoon evolutions (dashed curve) within the end leader (indexed by 99) frame. The initial platoon is shown in black while the final platoon is shown in blue which achieves a distance of 3030 between neighboring vehicles. We can verify that all the trajectories of the relative positions of the neighboring vehicles evolve within the transient performance bounds as shown in the bottom figure. Therefore, we can conclude that the target multi-vehicle platoon is achieved while satisfying the transient constraints.

Refer to caption
Figure 4: Platoon displacements (top) and trajectories of relative positions (bottom).

Multi-robot coordination: we then consider a multi-robot coordination example in 22 dimensions under more complex communication graphs, and we show that different choices of leader robots influence the fulfillment of the transient constraints. Suppose all robots are initialized at the origin and the target relative position-based formation is p53des=[9.8,0]T,p52des=[0,9.8]T,p62des=[9.8,9.8]T,p76des=[9.8,9.8]T,p54des=[0,9.8]T,p84des=[9.8,9.8]T,p98des=[9.8,9.8]T,p51des=[9.8,0]T,p510des:=p5desp10des=[9.8,9.8]T,p511des:=p5desp11des=[9.8,9.8]Tp^{des}_{53}=[-9.8,0]^{T},p^{des}_{52}=[0,-9.8]^{T},p^{des}_{62}=[9.8,9.8]^{T},p^{des}_{76}=[9.8,-9.8]^{T},p^{des}_{54}=[0,9.8]^{T},p^{des}_{84}=[9.8,-9.8]^{T},p^{des}_{98}=[9.8,9.8]^{T},p^{des}_{51}=[9.8,0]^{T},p^{des}_{510}:=p^{des}_{5}-p^{des}_{10}=[9.8,-9.8]^{T},p^{des}_{511}:=p^{des}_{5}-p^{des}_{11}=[9.8,9.8]^{T}. We choose, without loss of generality, the same ρx¯i\rho_{\bar{x}_{i}} for all edges as ρx¯i(t)=15et+0.1\rho_{\bar{x}_{i}}(t)=15e^{-t}+0.1 and the same transformed function Tx¯iT_{\bar{x}_{i}} as Tx¯i(x¯^i)=ln(1+x¯^i1x¯^i)T_{\bar{x}_{i}}(\hat{\bar{x}}_{i})=\ln\left(\frac{1+\hat{\bar{x}}_{i}}{1-\hat{\bar{x}}_{i}}\right). Note that the presented example generally covers a multi-robot dispersion task and the transient constraints encompass both the maintenance of connectivity and the transient performance requirements. Next, we focus on verifying the results derived in Theorem 13, and we consider different positions of the leader robots, which are depicted in Fig. 3.AA and Fig. 3.BB respectively. In Fig. 3.AA, the graph is described by 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}) with 𝒱={v1,v4,v5,v6,v10,v11}\mathcal{V}_{\mathscr{F}}=\{v_{1},v_{4},v_{5},v_{6},v_{10},v_{11}\} and 𝒱={v2,v3,v7,v8,v9}\mathcal{V}_{\mathscr{L}}=\{v_{2},v_{3},v_{7},v_{8},v_{9}\}. The edge set is ={e1,,e11}\mathcal{E}=\{e_{1},\dots,e_{11}\}. We first construct the maximum follower-end subgraph 𝒢=(𝒱,)\mathcal{G}^{\star}=(\mathcal{V}^{\star},\mathcal{E}^{\star}) of 𝒢\mathcal{G} according to Definition 11, and obtain that 𝒱=𝒱{v8,v9}\mathcal{V}^{\star}=\mathcal{V}\setminus\{v_{8},v_{9}\}, ={e6,e7}\mathcal{E}^{\star}=\mathcal{E}\setminus\{e_{6},e_{7}\}. With respect to 𝒢\mathcal{G}^{\star}, we have ={e5,e8,e9,e10}\mathcal{E}^{\star}_{\mathscr{F}\mathscr{F}}=\{e_{5},e_{8},e_{9},e_{10}\}. For these follower-follower edges, we obtain 𝒞ei𝒳i{min(|(𝒞ei)|4,2)}+|Ei|=|Ei|=5>2,i{5,8,9,10}\sum\limits_{\mathcal{C}_{e_{i}}\in\mathcal{X}_{i}}\left\{\min(|\mathcal{E}(\mathcal{C}_{e_{i}})|-4,2)\right\}+|E_{i}|=|E_{i}|=5>2,i\in\{5,8,9,10\}, which means that condition (9) as in Theorem 13 does not hold. Furthermore, there exist 22 FLF paths in 𝒢\mathcal{G}^{\star}, i.e., p1=v5v2v6p_{1}=v_{5}v_{2}v_{6} and p2=v5v3v7v6p_{2}=v_{5}v_{3}v_{7}v_{6}. For both FLF paths, we then operate complete decomposition according to Definition 12. We can verify that condition (10) holds for p2p_{2} since |(𝒞p2)|=5<2|(p2)|=6|\mathcal{E}(\mathcal{C}_{p_{2}})|=5<2|\mathcal{E}(p_{2})|=6. However, for the FLF path p1p_{1}, we have that 𝒞p1𝒴1{min(|(𝒞p1)|2|(p1)|2,2)}+|F1|=min(542,2)+4=3>2\sum\limits_{\mathcal{C}_{p_{1}}\in\mathcal{Y}_{1}}\left\{\min(|\mathcal{E}(\mathcal{C}_{p_{1}})|-2|\mathcal{E}(p_{1})|-2,2)\right\}+|F_{1}|=\min(5-4-2,2)+4=3>2 since F1={e5,e8,e9,e10}F_{1}=\{e_{5},e_{8},e_{9},e_{10}\}. Then, the simulation results when applying the PPC strategy (6) for the leader-follower topology as in Fig. 3.AA are shown in Fig. 5(a). We can see that the performance bounds are violated since the necessary and sufficient conditions (9), (10) in Theorem 13 do not hold. As a comparison, we change the positions of the leader robots as shown in Fig. 3.BB, where the new graph 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}) has the respective follower and leader sets as 𝒱={v1,v2,v3,v4,v10,v11}\mathcal{V}_{\mathscr{F}}=\{v_{1},v_{2},v_{3},v_{4},v_{10},v_{11}\} and 𝒱={v5,v6,v7,v8,v9}\mathcal{V}_{\mathscr{L}}=\{v_{5},v_{6},v_{7},v_{8},v_{9}\}. The edge set is ={e1,,e11}\mathcal{E}=\{e_{1},\dots,e_{11}\}. Similarly, we first construct the maximum follower-end subgraph 𝒢=(𝒱,)\mathcal{G}^{\star}=(\mathcal{V}^{\star},\mathcal{E}^{\star}) of 𝒢\mathcal{G} according to Definition 11, and derive that 𝒱=𝒱{v8,v9}\mathcal{V}^{\star}=\mathcal{V}\setminus\{v_{8},v_{9}\}, ={e6,e7}\mathcal{E}^{\star}=\mathcal{E}\setminus\{e_{6},e_{7}\}. With respect to 𝒢\mathcal{G}^{\star}, we now have =\mathcal{E}^{\star}_{\mathscr{F}\mathscr{F}}=\emptyset since there is no follower-follower edge. Then, we check whether condition (10) holds for all FLF paths. For the FLF paths viv5vjv_{i}v_{5}v_{j} in 𝒢\mathcal{G}^{\star} such that i,j{1,4,10,11}i,j\in\{1,4,10,11\} and iji\neq j, note that they do not have neighboring edges in 𝒢\mathcal{G}^{\star}. For example, the FLF path p1=v1v5v11p_{1}=v_{1}v_{5}v_{11} does not have any neighboring edges in 𝒢\mathcal{G}^{\star} according to Definition 10. We can conclude that condition (10) holds for p1p_{1} since 𝒞p1𝒴1{min(|(𝒞p1)|2|(p1)|2,2)}+|F1|=|F1|=0<2\sum\limits_{\mathcal{C}_{p_{1}}\in\mathcal{Y}_{1}}\left\{\min(|\mathcal{E}(\mathcal{C}_{p_{1}})|-2|\mathcal{E}(p_{1})|-2,2)\right\}+|F_{1}|=|F_{1}|=0<2. Next, for the FLF paths viv5vjv_{i}v_{5}v_{j} in 𝒢\mathcal{G}^{\star} such that i{1,4,10,11}i\in\{1,4,10,11\} and j{2,3}j\in\{2,3\}, note that they have a single neighboring edge e3e_{3} if j=2j=2 or e11e_{11} if j=3j=3 in 𝒢\mathcal{G}^{\star}. For example, the FLF path p2=v1v5v2p_{2}=v_{1}v_{5}v_{2} has only one neighboring edge in 𝒢\mathcal{G}^{\star}, i.e., e3e_{3}. We can conclude that condition (10) holds for p2p_{2} since 𝒞p2𝒴2{min(|(𝒞p2)|2|(p2)|2,2)}+|F2|=|F2|=1<2\sum\limits_{\mathcal{C}_{p_{2}}\in\mathcal{Y}_{2}}\left\{\min(|\mathcal{E}(\mathcal{C}_{p_{2}})|-2|\mathcal{E}(p_{2})|-2,2)\right\}+|F_{2}|=|F_{2}|=1<2. The remaining FLF paths are p3=v2v6v7v3p_{3}=v_{2}v_{6}v_{7}v_{3} and p4=v2v5v3p_{4}=v_{2}v_{5}v_{3}. They both belong to the cycle composed by the edges e1,e2,e3,e4,e11e_{1},e_{2},e_{3},e_{4},e_{11}. We can verify that condition (10) holds for p3p_{3} since |(𝒞p3)|=5<2|(p3)|=6|\mathcal{E}(\mathcal{C}_{p_{3}})|=5<2|\mathcal{E}(p_{3})|=6. Moreover, for the FLF path p4p_{4}, we have that 𝒞p4𝒴4{min(|(𝒞p4)|2|(p4)|2,2)}+|F4|=min(542,2)=1<2\sum\limits_{\mathcal{C}_{p_{4}}\in\mathcal{Y}_{4}}\left\{\min(|\mathcal{E}(\mathcal{C}_{p_{4}})|-2|\mathcal{E}(p_{4})|-2,2)\right\}+|F_{4}|=\min(5-4-2,2)=-1<2 since F4=F_{4}=\emptyset. Then, the simulation results when applying the PPC strategy (6) for the leader-follower topology as in Fig. 3.BB are shown in Fig. 5(b). We can observe that all the trajectories evolve within the prescribed performance bounds since the necessary and sufficient conditions (9), (10) in Theorem 13 hold. We can conclude that for the leader-follower network as shown in Fig. 3.BB, since the necessary and sufficient conditions (9), (10) in Theorem 13 hold, the target formation is achieved within the prescribed performance bounds (as shown in Fig. 5(b)) when applying the PPC strategy (6).

Refer to caption
(a) Trajectories of the shifted relative positions under graph as in Fig. 3.AA.
Refer to caption
(b) Trajectories of the shifted relative positions under graph as in Fig. 3.BB.
Figure 5: Trajectories of the shifted relative positions.

6 Conclusions

In this paper, we have investigated conditions on the graph topology such that relative position-based formation within prescribed transient bounds can be achieved for leader-follower multi-agent systems. For both tree graphs and general graphs with cycles, necessary conditions on the graph topology are firstly proposed under which we can potentially design the leaders to achieve the target formation with prescribed performance guarantees, which are further extended to necessary and sufficient conditions on the leader-follower graph topology such that the target formation can be achieved within the prescribed performance bounds. Future research includes considering other transient approaches and also investigating leader selection problems.

References

  • (1)
  • Aguilar & Gharesifard (2014) Aguilar, C. O. & Gharesifard, B. (2014), ‘Graph controllability classes for the laplacian leader-follower dynamics’, IEEE transactions on automatic control 60(6), 1611–1623.
  • Bechlioulis & Rovithakis (2008) Bechlioulis, C. P. & Rovithakis, G. A. (2008), ‘Robust adaptive control of feedback linearizable mimo nonlinear systems with prescribed performance’, IEEE Transactions on Automatic Control 53(9), 2090–2099.
  • Chang et al. (2003) Chang, D. E., Shadden, S. C., Marsden, J. E. & Olfati-Saber, R. (2003), Collision avoidance for multiple agent systems, in ‘42nd IEEE International Conference on Decision and Control (IEEE Cat. No. 03CH37475)’, Vol. 1, IEEE, pp. 539–543.
  • Chen & Dimarogonas (2020a) Chen, F. & Dimarogonas, D. V. (2020a), Further results on leader-follower multi-agent formation control with prescribed performance guarantees, in ‘59th IEEE Conference on Decision and Control’, IEEE.
  • Chen & Dimarogonas (2020b) Chen, F. & Dimarogonas, D. V. (2020b), ‘Leader-follower formation control with prescribed performance guarantees’, IEEE Transactions on Control of Network Systems .
  • Dai et al. (2019) Dai, S.-L., He, S., Chen, X. & Jin, X. (2019), ‘Adaptive leader–follower formation control of nonholonomic mobile robots with prescribed transient and steady-state performance’, IEEE Transactions on Industrial Informatics 16(6), 3662–3671.
  • Egerstedt et al. (2012) Egerstedt, M., Martini, S., Cao, M., Camlibel, K. & Bicchi, A. (2012), ‘Interacting with networks: How does structure relate to controllability in single-leader, consensus networks?’, IEEE Control Systems 32(4), 66–73.
  • Fax & Murray (2004) Fax, J. A. & Murray, R. M. (2004), ‘Information flow and cooperative control of vehicle formations’, IEEE transactions on automatic control 49(9), 1465–1476.
  • Fitch & Leonard (2016) Fitch, K. & Leonard, N. E. (2016), Optimal leader selection for controllability and robustness in multi-agent networks, in ‘2016 European Control Conference (ECC)’, IEEE, pp. 1550–1555.
  • Goldin & Raisch (2010) Goldin, D. & Raisch, J. (2010), Controllability of second order leader-follower systems, in ‘2nd IFAC Workshop on Distributed Estimation and Control in Networked Systems 2010-NecSys 10’, pp. 233–238.
  • He et al. (2018) He, S., Wang, M., Dai, S.-L. & Luo, F. (2018), ‘Leader–follower formation control of usvs with prescribed performance and collision avoidance’, IEEE transactions on industrial informatics 15(1), 572–581.
  • Katsoukis & Rovithakis (2016) Katsoukis, I. & Rovithakis, G. A. (2016), Output feedback leader-follower with prescribed performance guarantees for a class of unknown nonlinear multi-agent systems, in ‘2016 24th Mediterranean Conference on Control and Automation (MED)’, IEEE, pp. 1077–1082.
  • Mesbahi & Egerstedt (2010) Mesbahi, M. & Egerstedt, M. (2010), Graph theoretic methods in multiagent networks, Vol. 33, Princeton University Press.
  • Patterson et al. (2016) Patterson, S., McGlohon, N. & Dyagilev, K. (2016), ‘Optimal k-leader selection for coherence and convergence rate in one-dimensional networks’, IEEE Transactions on Control of Network Systems 4(3), 523–532.
  • Pirani et al. (2017) Pirani, M., Shahrivar, E. M., Fidan, B. & Sundaram, S. (2017), ‘Robustness of leader–follower networked dynamical systems’, IEEE Transactions on Control of Network Systems 5(4), 1752–1763.
  • Rahmani et al. (2009) Rahmani, A., Ji, M., Mesbahi, M. & Egerstedt, M. (2009), ‘Controllability of multi-agent systems from a graph-theoretic perspective’, SIAM Journal on Control and Optimization 48(1), 162–186.
  • Sun et al. (2017) Sun, C., Hu, G. & Xie, L. (2017), ‘Controllability of multiagent networks with antagonistic interactions’, IEEE Transactions on Automatic Control 62(10), 5457–5462.
  • Tanner (2004) Tanner, H. G. (2004), On the controllability of nearest neighbor interconnections, in ‘Decision and Control, 2004. CDC. 43rd IEEE Conference on’, Vol. 3, IEEE, pp. 2467–2472.
  • Wei et al. (2018) Wei, C., Luo, J., Yin, Z. & Yuan, J. (2018), ‘Leader-following consensus of second-order multi-agent systems with arbitrarily appointed-time prescribed performance’, IET Control Theory & Applications 12(16), 2276–2286.
  • Yazıcıoğlu et al. (2016) Yazıcıoğlu, A. Y., Abbas, W. & Egerstedt, M. (2016), ‘Graph distances and controllability of networks’, IEEE Transactions on Automatic Control 61(12), 4125–4130.
  • Yazicioğlu & Egerstedt (2013) Yazicioğlu, A. Y. & Egerstedt, M. (2013), Leader selection and network assembly for controllability of leader-follower networks, in ‘American Control Conference (ACC), 2013’, IEEE, pp. 3802–3807.
  • Zavlanos et al. (2011) Zavlanos, M. M., Egerstedt, M. B. & Pappas, G. J. (2011), ‘Graph-theoretic connectivity control of mobile robot networks’, Proceedings of the IEEE 99(9), 1525–1540.
  • Zelazo & Mesbahi (2011) Zelazo, D. & Mesbahi, M. (2011), ‘Edge agreement: Graph-theoretic performance bounds and passivity analysis’, IEEE Transactions on Automatic Control 56(3), 544–555.
  • Zhang et al. (2011) Zhang, S., Camlibel, M. K. & Cao, M. (2011), Controllability of diffusively-coupled multi-agent systems with general and distance regular coupling topologies, in ‘Decision and Control and European Control Conference (CDC-ECC), 2011 50th IEEE Conference on’, IEEE, pp. 759–764.
  • Zhang et al. (2014) Zhang, S., Cao, M. & Camlibel, M. K. (2014), ‘Upper and lower bounds for controllable subspaces of networks of diffusively coupled agents’, IEEE Transactions on Automatic control 59(3), 745–750.