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

Asynchronous Distributed Consensus with Minimum Communication

Vishal Sawant, Debraj Chakraborty and Debasattam Pal Vishal Sawant, Debraj Chakraborty (Corresponding author) and Debasattam Pal are with the Department of Electrical Engineering, Indian Institute of Technology Bombay, Mumbai - 400076, India. Email: {vsawant, dc, debasattam}@ee.iitb.ac.in
Abstract

In this paper, the communication effort required in a multi-agent system (MAS) is minimized via an explicit optimization formulation. The paper considers a MAS of single-integrator agents with bounded inputs and a time-invariant communication graph. A new model of discrete asynchronous communication and a distributed consensus protocol based on it, are proposed. The goal of the proposed protocol is to minimize the aggregate number of communication instants of all agents, required to steer the state trajectories inside a pres-specified bounded neighbourhood within a pre-specified time. Due to information structure imposed by the underlying communication graph, an individual agent does not know the global parameters in the MAS, which are required for the above-mentioned minimization. To counter this uncertainty, the worst-case realizations of the global parameters are considered, which lead to min-max type optimizations. The control rules in the proposed protocol are obtained as the closed form solutions of these optimization problems. Hence, the proposed protocol does not increase the burden of run-time computation making it suitable for time-critical applications.

Index Terms:
Consensus, Distributed Optimization, Communication Cost, Asynchronous Communication

I Introduction

In the recent years, multi-agent systems (MASs) have received tremendous attention due to their extensive applications in unmanned aerial vehicles (UAVs) [3], sensor networks [14], power grids [18], industrial robotics [26] etc. A typical MAS consists of a group of agents which collaborate to achieve desired objectives. Often, the objective is to achieve consensus, i.e., to drive the states of all agents into an agreement. To achieve this objective, agents in MAS need to exchange information with each other over some communication network. It is well known that communication network is an expensive resource [2], [22]. Hence, with the intent of reducing communication effort, a few approaches have been proposed in [9], [27], [28] etc. However, except in our preliminary work [25], the minimization of communication effort via explicit optimization, has never been addressed.

The problem of achieving consensus of MAS with single-integrator agents was first extensively analyzed in [15]. After that, various consensus protocols for higher order agents were developed in [11], [23], [24] etc. These protocols are either based on discrete synchronous communication or on continuous communication. A global synchronization clock is necessary for the implementation of synchronous protocols, which is often a major practical constraint [16]. On the other hand, the energy consumption of agent’s transponder is proportional to the number and duration of transmissions [7]. Thus, continuous transmission limits the life of agent’s battery and hence, its flight time [4], [19]. Second, continuous transmissions over a shared bandwidth-limited channel by multiple agents may lead to congestion of communication channel [2], [22]. And finally, MASs such as a group of UAVs are often used for stealthy military applications [17], [34]. In such scenarios, it is of strategic advantage to keep radio transmissions to a minimum in order to avoid detection by the enemy. For all these applications, it is necessary to develop an asynchronous/intermittent communication based consensus protocol which minimizes communication effort, i.e., the number and/or duration of transmissions.

In order to reduce the required communication effort in a MAS, a few indirect approaches have been proposed in the literature. In self-triggered control [13] based consensus protocols [8], [32], the next communication instant is pre-computed based on the current state. Event-triggered control [13] based consensus protocols [9], [35] initiate communication only when a certain error state reaches a predetermined threshold. Intermittent communication based consensus protocols are investigated in [27], [28] and [29]. Consensus protocols based on asynchronous information exchange for a MAS with single-integrator agents are developed in [5], [6], [10] and [31]. Other work on asynchronous consensus include [1], [12], [30] and [33]. All of the above protocols result in the reduced communication effort as compared to the conventional continuous communication based protocols. However, there is no explicit minimization of communication effort and hence, the above protocols can result in sub-optimal communication performance.

To overcome this issue, in this paper, we develop a distributed protocol which minimizes the communication effort required to maintain the consensus of single-integrator agents, via explicit optimization. This protocol is based on discrete asynchronous communication. Our notion of consensus is less stringent than the conventional one [15], in that we only require the difference between neighbouring agents’ states (i.e., the local disagreement) to reduce below a pre-specified bound. In the proposed protocol, communication occurs at discrete time instants, namely update/communication instants, at which agents access the states of neighbouring agents and then based on that information, update their control. Thus, in the proposed protocol, the number of update instants is a good measure of communication effort. Hence, our basic objective is to minimize the aggregate number of update instants of all agents in the MAS. Now, intuitively, if the inter-update durations are increased, then the number of update instants decreases. However, increase in the inter-update durations increases the time required to achieve the consensus bound. Hence, the problem of minimization of the number of update instants is well-posed only when the consensus time is included in the formulation. Evidently, the time required for a MAS to achieve consensus depends on the initial configuration of the agents. Hence, to make the minimization of the number of update instants well-posed, we require that the local disagreements of all agents be steered below a pre-specified consensus bound within a pre-specified time, which is expressed as a function of the initial condition of the MAS.

Due to communication structure imposed by the network, the above-mentioned minimization is a decentralized optimal control [20] problem. Because of the said imposition, an individual agent does not have global information such as the number of agents in the MAS, the underlying communication graph, the complete initial condition of the MAS etc. To counter this lack of global information, we require that the consensus constraint be satisfied for any number of agents, any communication graph and any initial condition. Further, due to the imposed communication structure, an individual agent can not predict the control inputs of the neighbouring agents. In order to guard against the resulting uncertainty, the control inputs in the proposed protocol are obtained as a solution of certain max-min optimizations (Sections V-A and V-B). We obtain the closed form solutions of these optimizations (see (7)-(10)) and hence, extensive numerical computations are not required for their implementation. This makes the proposed protocol suitable for time-critical applications.

Our contributions in this paper are summarized as follows:

  1. 1.

    We develop a discrete asynchronous communication based distributed consensus protocol for a MAS with single-integrator agents (Section III).

  2. 2.

    The proposed protocol minimizes the aggregate number of update instants under the constraint of steering the local disagreements of all agents below a pre-specified limit within a pre-specified time (Theorems 15 and 16).

  3. 3.

    The control rules in the proposed protocol are solution of certain max-min optimizations (Sections V-A and V-B). We obtain the closed form expressions of the corresponding optimal controls (Lemmas 11 and 12).

The current paper is an extension of our work in [25] in three major ways. First, it was assumed in [25] that the initial local disagreements between agents are confined below a pre-specified bound. In the current paper, no such assumption on initial conditions has been made. Second, the protocol in [25] solves the minimization problem for a specific consensus time. On the other hand, the protocol developed in the current paper solves the minimization problem for a general pre-specified consensus time. Finally, in the current paper, the effect of pre-specified consensus time on optimal communication cost is analyzed. Such analysis was not presented in [25].

The remaining part of this paper is organized as follows. In Section II, the problem of minimizing the number of communication instants is formulated. In Section III, a distributed consensus protocol is proposed, which will be shown to be the solution for the special case of the formulated problem, in Sections IV and V. The protocol proposed in Section III is extended for the general case of the formulated problem, in Section VI. In Section VII, the simulation results are presented. The paper is concluded in Section VIII with future directions.

II Preliminaries and Problem formulation

II-A Graphs

A graph G=(V,E)G=(V,E) is a finite set of nodes VV connected by a set of edges E(V×V)E\subseteq(V\times V). An edge between nodes ii and jj is represented by an ordered pair (i,j)E(i,j)\in E. A graph GG is said to be simple if (i,i)E,iV(i,i)\not\in E,~\forall i\in V. A graph GG is said to be undirected if (i,j)E(i,j)\in E implies (j,i)E(j,i)\in E. In an undirected graph GG, if (i,j)E(i,j)\in E (and equivalently (j,i)E(j,i)\in E), then the nodes ii and jj are said to be neighbours of each other. A path between nodes ii and jj in an undirected graph GG is a sequence of edges (i,k1),(k1,k2),,(kr1,kr),(kr,j)E(i,k_{1}),(k_{1},k_{2}),\dots,(k_{r-1},k_{r}),(k_{r},j)\in E. An undirected graph GG is said to be connected if there exists a path between any two nodes in GG. Let nin_{i} denote the number of neighbours of node ii and |V||V| denote the cardinality of set VV. Then, the Laplacian matrix L|V|×|V|L\in{\mathbb{R}}^{|V|\times|V|} of a simple undirected graph G=(V,E)G=(V,E) is defined as

Li,j:={ni, if i=j1, if ij and (i,j)E0, if ij and (i,j)EL_{i,j}:=\begin{cases}n_{i},&\text{ if }~~~~i=j\\ -1,&\text{ if }~~~~i\neq j~~\text{ and }~~(i,j)\in E\\ 0,&\text{ if }~~~~i\neq j~~\text{ and }~~(i,j)\not\in E\end{cases}

II-B System description

Consider a multi-agent system (MAS) of nn single-integrator agents, labeled as a1a_{1}, a2,,ana_{2},\dots,a_{n}, with dynamics

x˙i(t)=ui(t),i=1,,n\dot{x}_{i}(t)=u_{i}(t),~~~~~i=1,\dots,n (1)

where xi(t)x_{i}(t)\in\mathbb{R} and ui(t)u_{i}(t)\in\mathbb{R} are the state and the control input of agent aia_{i}, respectively. Let X:=[x1,x2,,xn]X:=[x_{1},x_{2},\dots,x_{n}] and 𝐮:=[u1,u2,,un]{\bf{u}}:=[u_{1},u_{2},\dots,u_{n}] be the augmented state and control vectors of MAS (1), respectively. Define the set P:={1,2,,n}P:=\{1,2,\dots,n\}.

Let GG be a time-invariant simple undirected graph, whose nodes represent agents in MAS (1) whereas the edges represent the communication links between agents, over which they exchange information with their neighbours. Let SiS_{i} be the set of indices of neighbours of agent aia_{i}. Note that iSii\not\in S_{i} as GG is a simple graph. The cardinality of the set SiS_{i} is denoted by nin_{i}. Let LL denote the Laplacian matrix of GG.

We make the following assumptions about MAS (1):

  1. 1.

    The control input uiu_{i} of each agent belongs to the set

    𝒰:={u||u(t)|β,t0}\mathcal{U}:=\{u\in\mathcal{M}~|~|u(t)|\leq\beta,~\forall t\geq 0\}

    where \mathcal{M} denotes the set of measurable functions from [0,)[0,\infty) to \mathbb{R}.

  2. 2.

    The communication graph GG is connected.

  3. 3.

    The communication delay is zero.

II-C Consensus

In this paper, we will be using two notions of consensus, which we define next.

Definition 1.

MAS (1) is said to have achieved conventional consensus at instant t~\widetilde{t} if

t~:=inf{t^|xi(t)=xj(t),tt^,i,jP}<\widetilde{t}:=\inf\big{\{}\widehat{t}~\big{|}~x_{i}(t)=x_{j}(t),~~\forall t\geq\widehat{t},~~\forall i,j\in P\big{\}}<\infty

Define Z(t)=[z1(t),,zn(t)]:=LX(t)Z(t)=[z_{1}(t),\dots,z_{n}(t)]:=LX(t). Then, it follows from the definition of the Laplacian matrix that

zi(t)=jSi(xi(t)xj(t)),iPz_{i}(t)=\sum_{j\in S_{i}}\big{(}x_{i}(t)-x_{j}(t)\big{)},~~~~~\forall i\in P (2)

As ziz_{i} is the sum of differences of agent aia_{i}’s state with its neighbours, we call it the local disagreement of agent aia_{i}. It is well known [15] that for a MAS with connected, time-invariant communication graph, conventional consensus at instant t~\widetilde{t} is equivalent to zi(t)=0,tt~,iPz_{i}(t)=0,~\forall t\geq\widetilde{t},~\forall i\in P. However, in many practical applications, it is not necessary that each ziz_{i} becomes exactly zero. It is sufficient if each |zi||z_{i}| remains below a prespecified consensus bound. This motivates our next notion of consensus, namely α\alpha-consensus.

Definition 2.

Let α+\alpha\in{\mathbb{R}}^{+} be the prespecified consensus bound. MAS (1) is said to have achieved α\alpha-consensus at instant t~\widetilde{t} if

t~:=inf{t^||zi(t)|α,tt^,iP}<\widetilde{t}:=\inf\big{\{}\widehat{t}~\big{|}~|z_{i}(t)|\leq\alpha,~~\forall t\geq\widehat{t},~~\forall i\in P\big{\}}<\infty

II-D Communication model

In this paper, we consider a discrete communication model. Let tikt^{k}_{i} denote the kkth update instant (also referred to as the communication instant) of agent aia_{i}, at which it accesses the state information of its neighbours and then, based on that information, updates its control. Our communication model is asynchronous, i.e., the update instants tikt^{k}_{i}’s of two different agents need not coincide.

As the communication model is discrete, the number of update instants is a good measure of communication effort. Hence, we define the communication cost of agent aia_{i}, denoted by Ci(t)C_{i}(t), as the number of update instants of agent aia_{i} in the time interval [0,t][0,t]. Then, we define the aggregate communication cost of MAS (1), denoted by CMAS(t)C_{MAS}(t), as

CMAS(t):=iPCi(t)C_{MAS}(t):=\sum_{i\in P}C_{i}(t) (3)

II-E Problem formulation

Let X(0)X(0) be an initial condition of MAS (1), α\alpha be the prespecified consensus bound and TT be the prespecified consensus time. Our objective is to develop a protocol which minimizes the communication cost CMAS(T)C_{MAS}(T), under the constraint of achieving α\alpha-consensus of MAS (1) within time TT.

As discussed in Section I, due to information structure imposed by graph GG, an individual agent in MAS (1) has access only to its own information and that of its neighbours. Therefore, the proposed protocol needs to be distributed, i.e., based only on the local information. In addition, an individual agent does not have global information such as the number of agents nn in MAS (1), the structure of the communication graph GG, the complete initial condition X(0)X(0) etc. To address these uncertainties, the proposed protocol must be able to achieve α\alpha-consensus of MAS (1) within the prespecified time TT, for any nn, any connected GG and any X(0)nX(0)\in{\mathbb{R}}^{n}.

Since the input set 𝒰\mathcal{U} is magnitude bounded, for a fixed TT, it will not be possible to achieve α\alpha-consensus within time TT, for every X(0)nX(0)\in{\mathbb{R}}^{n}. Thus, the consensus time must be specified as a function of X(0)X(0). To highlight this dependence on X(0)X(0), we modify the notation of the prespecified consensus time from TT to T(X(0))T\big{(}X(0)\big{)}. Similarly, the communication costs CiC_{i} and CMASC_{MAS} depend on X(0)X(0). Thus, we modify their notations from Ci(t)C_{i}(t) and CMAS(t)C_{MAS}(t) to Ci(t,X(0))C_{i}\big{(}t,X(0)\big{)} and CMAS(t,X(0))C_{MAS}\big{(}t,X(0)\big{)}, respectively. Now, we formalize our objective as follows:

Problem 3.

Consider MAS (1) with initial condition X(0)nX(0)\in{\mathbb{R}}^{n} and connected communication graph GG. Let Ψ(n)\Psi(n) denote the set of connected graphs with nn nodes. Let T(X(0))T\big{(}X(0)\big{)} be the prespecified consensus time which is expressed as a function of X(0)X(0). Develop, if possible, a discrete asynchronous communication based protocol, i.e., admissible control 𝐮=[u1,,un]{\bf{u}}^{\ast}=[u_{1}^{\ast},\dots,u^{\ast}_{n}], adhering to graph GG, which is solution of the following optimization:

𝐮=\displaystyle{\bf{u}}^{\ast}=~ minui𝒰iPCMAS(T(X(0)),X(0))\displaystyle\min_{\begin{subarray}{c}u_{i}\in\mathcal{U}\\ \forall i\in P\end{subarray}}~~C_{MAS}\Big{(}T\big{(}X(0)\big{)},X(0)\Big{)}
s.t. |zi(t)|α,tT(X(0)),iP\displaystyle\text{ s.t. }~~~|z_{i}(t)|\leq\alpha,~~\forall t\geq T\big{(}X(0)\big{)},~~\forall i\in P (4)
n,GΨ(n),X(0)n\displaystyle\hskip 57.75905pt\forall n,~~\forall G\in\Psi(n),~~\forall X(0)\in{\mathbb{R}}^{n}

II-F Choosing T(X(0))T\big{(}X(0)\big{)}

In Problem 3, the time T(X(0))T\big{(}X(0)\big{)} can be specified as any function of X(0)X(0). However, in practical applications, it is desirable to set T(X(0))T\big{(}X(0)\big{)} to the minimum feasible value.

In [21], the time-optimal control rule is proposed which achieves conventional consensus of MAS (1) in minimum time. This control rule and the corresponding consensus time, denoted by uiu^{\star}_{i} and T(X(0))T^{\ast}\big{(}X(0)\big{)}, respectively, are presented below. Let X(0)=[x1(0),,xn(0)]X(0)=[x_{1}(0),\dots,x_{n}(0)] be an initial condition of MAS (1). Define xmin(0):=min{x1(0),,xn(0)}x^{min}(0):=\min\big{\{}x_{1}(0),\dots,x_{n}(0)\big{\}} and xmax(0):=max{x1(0),,xn(0)}x^{max}(0):=\max\big{\{}x_{1}(0),\dots,x_{n}(0)\big{\}}. Recall that ziz_{i} denote the local disagreement of agent aia_{i}. Let signsign denote the standard signum function. Then, the time-optimal consensus rule from [21] is

ui(t)=βsign(zi(t)),t0,iPu^{\star}_{i}(t)=-\beta\hskip 2.84544ptsign\big{(}z_{i}(t)\big{)},~~~~\forall t\geq 0,~~~~\forall i\in P (5)

with the corresponding consensus time

T(X(0))=xmax(0)xmin(0)2βT^{\ast}\big{(}X(0)\big{)}=\frac{x^{max}(0)-x^{min}(0)}{2\beta} (6)

Notice that the control rule (5) requires instantaneous access to the local disagreement ziz_{i}, and as a result, demands continuous communication. Then, it follows from the time optimality of T(X(0))T^{\ast}\big{(}X(0)\big{)} that a discrete communication based protocol cannot achieve α\alpha-consensus of MAS (1) within time T(X(0))T^{\ast}\big{(}X(0)\big{)}. This implies that Problem 3 is infeasible for T(X(0))T(X(0))T\big{(}X(0)\big{)}\leq T^{\ast}\big{(}X(0)\big{)} and we should assume T(X(0))>T(X(0))T\big{(}X(0)\big{)}>T^{\ast}\big{(}X(0)\big{)}. We further assume that T(X(0))2T(X(0))T\big{(}X(0)\big{)}\geq 2T^{\ast}\big{(}X(0)\big{)}. This assumption makes Problem 3 tractable and results in a particularly simple closed form solution for the control inputs. We solve Problem 3 for T(X(0))=2T(X(0))T\big{(}X(0)\big{)}=2T^{\ast}\big{(}X(0)\big{)} in Sections III-V and later extend it for T(X(0))>2T(X(0))T\big{(}X(0)\big{)}>2T^{\ast}\big{(}X(0)\big{)} in Section VI.

III Protocol \Romannum1\Romannum{1}: For T(X(0))=2T(X(0))T\big{(}X(0)\big{)}=2T^{\ast}\big{(}X(0)\big{)}

In this section, we present the protocol proposed for T(X(0))=2T(X(0))T\big{(}X(0)\big{)}=2T^{\ast}\big{(}X(0)\big{)}. We refer to this protocol as Protocol \Romannum1\Romannum{1}. Later, this protocol will be shown to be a solution of Problem 3 for T(X(0))=2T(X(0))T\big{(}X(0)\big{)}=2T^{\ast}\big{(}X(0)\big{)}.

Protocol \Romannum1\Romannum{1} has the following two elements:
11)  Computation of next update instant and control input
aa
)  At each update instant tik,k1t^{k}_{i},~k\geq 1, agent ai,iPa_{i},~i\in P, accesses the current states xj(tik)x_{j}(t^{k}_{i})’s of its neighbours aj,jSia_{j},~j\in S_{i}, and computes zi(tik)=jSi(xi(tik)xj(tik))z_{i}(t^{k}_{i})=\sum_{j\in S_{i}}\big{(}x_{i}(t^{k}_{i})-x_{j}(t^{k}_{i})\big{)}.
bb) After that, agent aia_{i} computes its next update instant tik+1t^{k+1}_{i} and control input uiu^{\ast}_{i} to be applied in the interval [tik,tik+1)[t^{k}_{i},t^{k+1}_{i}) as follows:
   \romannum1\romannum{1}) If |zi(tik)|α|z_{i}(t^{k}_{i})|\leq\alpha, then

tik+1\displaystyle t^{k+1}_{i} =tik+αβni\displaystyle=t^{k}_{i}+\frac{\alpha}{\beta n_{i}} (7)
ui(t)\displaystyle u^{\ast}_{i}(t) =zi(tik)αβ,t[tik,tik+1)\displaystyle=-\frac{z_{i}(t^{k}_{i})}{\alpha}\beta,\hskip 25.6073pt\forall t\in\big{[}t^{k}_{i},t^{k+1}_{i}\big{)} (8)

\romannum2\romannum{2}) If |zi(tik)|>α|z_{i}(t^{k}_{i})|>\alpha, then

tik+1\displaystyle t^{k+1}_{i} =tik+|zi(tik)|+α2βni\displaystyle=t^{k}_{i}+\frac{|z_{i}(t^{k}_{i})|+\alpha}{2\beta n_{i}} (9)
ui(t)\displaystyle u^{\ast}_{i}(t) =βsign(zi(tik)),t[tik,tik+1)\displaystyle=-\beta sign\big{(}z_{i}(t^{k}_{i})\big{)},~~~~\forall t\in\big{[}t^{k}_{i},t^{k+1}_{i}\big{)} (10)

cc) Then, agent aia_{i} broadcasts xi(tik)x_{i}(t^{k}_{i}) and ui(tik)u^{\ast}_{i}(t^{k}_{i}). This broadcast information is received by agent aia_{i}’s neighbours aj,jSia_{j},~j\in S_{i}, at the same instant tikt^{k}_{i}. The neighbours aj,jSia_{j},~j\in S_{i}, store this information with the reception time-stamp.

22) Accessing the states of neighbours at update instants
aa
) The time instant ti1=0t^{1}_{i}=0 is the first update instant of all agents ai,iPa_{i},~i\in P. At this instant, all agents broadcast their current states. Thus, each agent ai,iPa_{i},~i\in P, has a direct access to the current states xj(ti1)x_{j}(t^{1}_{i})’s of its neighbours aj,jSia_{j},~j\in S_{i}. For example, consider the communication graph GcG_{c} shown in Fig. 1 and the corresponding communication timeline shown in Fig. 2. At instant t=0t=0, agents a1a_{1}, a2a_{2} and a3a_{3} broadcast information to their neighbours.

112233
Figure 1: Communication graph GcG_{c}
0t1kt^{k}_{1}t3lt^{l}_{3}t2pt^{p}_{2}t1k+1t^{k+1}_{1}a1a_{1}a2a_{2}a3a_{3}
Figure 2: Timeline of communication over graph GcG_{c} in Fig. 1. Dark arrows indicate information transmissions.

bb) Let tik,k>1t^{k}_{i},~k>1, be any update instant of agent aia_{i}. Let tjl<tikt^{l}_{j}<t^{k}_{i} be the latest update instant of agent aj,jSia_{j},~j\in S_{i}, at which it had broadcast xj(tjl)x_{j}(t^{l}_{j}) and uj(tjl)u^{\ast}_{j}(t^{l}_{j}). At update instant tikt^{k}_{i}, agent aia_{i} accesses the stored information and retrieves xj(tjl)x_{j}(t^{l}_{j}) and uj(tjl)u^{\ast}_{j}(t^{l}_{j}) for all jSij\in S_{i}. For example, in the communication timeline shown in Fig. 2, at instants t1kt^{k}_{1} and t3lt^{l}_{3}, agent a2a_{2} receives information from its neighbours. Later, agent a2a_{2} uses this information at its update instant t2pt^{p}_{2}.
cc) It is known from (8) and (10) that every agent aj,jSia_{j},~j\in S_{i}, had applied control uj(t)=uj(tjl)u^{\ast}_{j}(t)=u^{\ast}_{j}(t^{l}_{j}) in the interval t[tjl,tik)t\in[t^{l}_{j},t^{k}_{i}). Using this fact, agent aia_{i} computes the current state xj(tik)x_{j}(t^{k}_{i}) of aj,jSia_{j},~\forall j\in S_{i}, as

xj(tik)=xj(tjl)+tjltikuj(tjl)𝑑tx_{j}\big{(}t^{k}_{i}\big{)}=x_{j}(t^{l}_{j})+\int_{t^{l}_{j}}^{t^{k}_{i}}\!u^{\ast}_{j}\big{(}t^{l}_{j}\big{)}~dt
Remark 4.

As per Assumption 3, the communication delay is zero. In addition, the time required for the computation of tik+1t^{k+1}_{i} and ui(tik)u^{\ast}_{i}(t^{k}_{i}) is negligible. This justifies the assumption that computation, transmission and reception of ui(tik)u^{\ast}_{i}(t^{k}_{i}) and xi(tik)x_{i}(t^{k}_{i}), happen at the same time instant tikt^{k}_{i}.

Remark 5.

It may appear from (8) and (10) that we have selected uiu^{\ast}_{i}’s which are constant over intervals [tik,tik+1)[t^{k}_{i},t^{k+1}_{i}), in order to simplify analysis. However, uiu^{\ast}_{i}’s in (8) and (10) are obtained as the solution of certain max-min optimizations (Sections V-A and V-B) and coincidently, they have this nice form.

IV α\alpha-consensus under Protocol \Romannum1\Romannum{1}

In this section, we show that Protocol \Romannum1\Romannum{1} achieves α\alpha-consensus of MAS (1) within time 2T(X(0))2T^{\ast}\big{(}X(0)\big{)}. However, before that, we present a necessary condition on feasible solutions of Problem 3, which will be utilized while proving the attainment of α\alpha-consensus within time 2T(X(0))2T^{\ast}\big{(}X(0)\big{)}.

IV-A Necessary condition on feasible solutions of Problem 3

A discrete communication based distributed protocol is said to be a feasible solution of Problem 3 if it satisfies constraint (4), i.e., achieves α\alpha-consensus of MAS (1) within time T(X(0))T\big{(}X(0)\big{)}, for any number of agents nn, any connected graph GG and any initial condition X(0)nX(0)\in{\mathbb{R}}^{n}. Recall that Ψ(n)\Psi(n) denotes the set of connected graphs with nn nodes. Then, the following lemma gives a necessary condition on the feasible solutions of Problem 3.

Lemma 6.

Consider MAS (1). Let T(X(0))T\big{(}X(0)\big{)} be the pre-specified consensus time. Let QQ be any discrete communication based distributed protocol. Under Protocol QQ, let t~i\widetilde{t}_{i} be the first time instant at which the local disagreement ziz_{i} of agent aia_{i} satisfies |zi(t~i)|α\big{|}z_{i}\big{(}~\!\widetilde{t}_{i}\big{)}\big{|}\leq\alpha. Then, Protocol QQ is a feasible solution of Problem 3 only if it ensures

|zi(t)|α,tt~i,n,GΨ(n),X(0)n\big{|}z_{i}(t)\big{|}\leq\alpha,\hskip 11.38092pt\forall t\geq\widetilde{t}_{i},\hskip 5.69046pt\forall n,\hskip 5.69046pt\forall G\in\Psi(n),\hskip 5.69046pt\forall X(0)\in{\mathbb{R}}^{n} (11)
Proof.

Recall that due to communication structure imposed by graph GG, an individual agent aia_{i} in MAS (1) does not know the complete initial condition X(0)X(0). Therefore, it does not know the exact value of T(X(0))T\big{(}X(0)\big{)} and how far t~i\widetilde{t}_{i} is from T(X(0))T\big{(}X(0)\big{)}. In fact, as in Example 9 presented below, there may exist a MAS of the form (1) in which t~i=T(X(0))\widetilde{t}_{i}=T\big{(}X(0)\big{)}. In such a case, violation of (11) results in the violation of constraint (4) in Problem 3. This contradicts the fact that Protocol QQ is a feasible solution of Problem 3. Hence, Protocol QQ must satisfy (11). This completes the proof. ∎

The following lemma shows that Protocol \Romannum1\Romannum{1} satisfies the necessary condition (11). The proof of this lemma relies on the derivation of control rules (7)-(10) in Protocol \Romannum1\Romannum{1}, which is deferred to Section V for better structure of the paper. Hence, we defer the proof of the said lemma to Section V-C.

Lemma 7.

Protocol \Romannum1\Romannum{1} satisfies the necessary condition (11).

IV-B α\alpha-consensus within time 2T(X(0))2T^{\ast}\big{(}X(0)\big{)}

The following theorem shows that Protocol \Romannum1\Romannum{1} achieves α\alpha-consensus of MAS (1) within time 2T(X(0))2T^{\ast}\big{(}X(0)\big{)}.

Theorem 8.

Consider MAS (1). Let α+\alpha\in{\mathbb{R}}^{+} be the pre-specified consensus bound. Let T(X(0))T^{\ast}\big{(}X(0)\big{)} be as defined in (6). Then, for every nn, every connected communication graph GG and every X(0)nX(0)\in{\mathbb{R}}^{n}, Protocol \Romannum1\Romannum{1} achieves α\alpha-consensus of MAS (1) in time less than or equal to 2T(X(0))2T^{\ast}\big{(}X(0)\big{)}.

Proof.

See the Appendix for the proof. ∎

According to Theorem 8, under Protocol \Romannum1\Romannum{1}, the duration 2T(X(0))2T^{\ast}\big{(}X(0)\big{)} is an upper bound on the time required to achieve α\alpha-consensus of MAS (1). Next, we show with the following example that there exists a MAS of the form (1) for which the α\alpha-consensus time under Protocol \Romannum1\Romannum{1} is arbitrarily close to 2T(X(0))2T^{\ast}\big{(}X(0)\big{)}.

Example 9.

Consider the connected graph GeG_{e} shown in Fig. 3, which has nn nodes. The subgraph GsG_{s} of GeG_{e} (inside the dashed square in Fig. 3) is a complete graph on n1n-1 nodes. The set of indices of neighbours of node 11 is S1={2,3,,r+1}S_{1}=\{2,3,\dots,r+1\}. Hence, the cardinality of S1S_{1} is rr.

Now, consider MAS (1) with communication graph GeG_{e}. Let the consensus bound and the control bound be α=3\alpha=3 and β=1\beta=1, respectively. Let the initial conditions of agents be x1(0)=0x_{1}(0)=0 and xi(0)=5,i=2,,nx_{i}(0)=5,~\forall i=2,\dots,n. Thus, xmin(0)=minixi(0)=0x^{min}(0)=\min_{i}x_{i}(0)=0 and xmax(0)=maxixi(0)=5x^{max}(0)=\max_{i}x_{i}(0)=5. Then, it follows from the definition of T(X(0))T^{\ast}\big{(}X(0)\big{)} in (6) that T(X(0))=2.5T^{\ast}\big{(}X(0)\big{)}=2.5 sec.

2233rrnnr{\footnotesize{r}}+{\footnotesize{+}}114411GsG_{s}
Figure 3: Communication graph GeG_{e} of the MAS in Example 9

The following theorem shows that under Protocol \Romannum1\Romannum{1}, the time required to achieve α\alpha-consensus of the MAS in Example 9 is arbitrarily close to 2T(X(0))2T^{\ast}\big{(}X(0)\big{)}.

Theorem 10.

Consider the MAS in Example 9. Let T(α,X(0))T\big{(}\alpha,X(0)\big{)} denote the time required by Protocol \Romannum1\Romannum{1} to achieve α\alpha-consensus of a MAS with initial condition X(0)X(0). Let ϵ>0\epsilon>0 be any real number. Then, there exist integers r=rϵr=r_{\epsilon} and n=nϵn=n_{\epsilon} such that the following holds.

(2T(X(0))ϵ)T(α,X(0))2T(X(0))\Big{(}2T^{\ast}\big{(}X(0)\big{)}-\epsilon\Big{)}\leq T\big{(}\alpha,X(0)\big{)}\leq 2T^{\ast}\big{(}X(0)\big{)} (12)
Proof.

See the Appendix for the proof. ∎

V Optimality of Protocol \Romannum1\Romannum{1}

The objective in Problem 3 is to minimize the number of update instants under the constraint of achieving α\alpha-consensus of MAS (1) within the pre-specified time. Intuitively, if the duration between the successive update instants is increased, then the number of update instants decreases. Motivated from this intuition, we obtain the solution of Problem 3 by maximizing the inter-update durations. We divide the solution process into two steps. First, we solve two maximization problems, one corresponding to |zi(tik)|α|z_{i}(t^{k}_{i})|\leq\alpha and the other corresponding to |zi(tik)|>α|z_{i}(t^{k}_{i})|>\alpha, in which the objective function is the inter-update duration. The control rules (7)-(8) and (9)-(10) in Protocol \Romannum1\Romannum{1} are solutions of these maximization problems, respectively. Later, we show how these two control rules together form the solution of Problem 3.

V-A Maximization of inter-update durations: For |zi(tik)|α\big{|}z_{i}\big{(}t^{k}_{i}\big{)}\big{|}\leq\alpha

Let zi(tik)z_{i}\big{(}t^{k}_{i}\big{)} be the local disagreement of agent aia_{i} at its update instant tikt^{k}_{i} such that |zi(tik)|α\big{|}z_{i}\big{(}t^{k}_{i}\big{)}\big{|}\leq\alpha. The goal of agent aia_{i} is to maximize the inter-update duration (tik+1tik)\big{(}t^{k+1}_{i}-t^{k}_{i}\big{)} by delaying its next update instant tik+1t^{k+1}_{i}. However, while doing so, agent aia_{i} must satisfy the necessary condition (11). For that purpose, agent aia_{i} needs to ensure that |zi(t)|α|z_{i}(t)|\leq\alpha for all t[tik,tik+1]t\in[t^{k}_{i},t^{k+1}_{i}]. Recall the definition of ziz_{i} in (2). Then, the evolution of ziz_{i} in the interval [tik,tik+1][t^{k}_{i},t^{k+1}_{i}] is given as

zi(t)=zi(tik)+nitiktui(τ)𝑑τjSi(tiktuj(τ)𝑑τ)z_{i}\big{(}t\big{)}=z_{i}\big{(}t^{k}_{i}\big{)}+n_{i}\int_{t^{k}_{i}}^{t}u_{i}(\tau)d\tau-\sum_{j\in S_{i}}\bigg{(}\int_{t^{k}_{i}}^{t}u_{j}(\tau)d\tau\bigg{)} (13)

This evolution depends on control inputs uiu_{i} and uj,jSiu_{j},~\forall j\in S_{i}, in the interval [tik,tik+1][t^{k}_{i},t^{k+1}_{i}]. Note that any instant t[tik,tik+1]t\in[t^{k}_{i},t^{k+1}_{i}] can be an update instant of a neighbouring agent aj,jSia_{j},~j\in S_{i}, at which aja_{j} updates its control. This updated value is not known to agent aia_{i} in advance, at instant tikt^{k}_{i}. Thus, while maximizing the inter-update duration (tik+1tik)\big{(}t^{k+1}_{i}-t^{k}_{i}\big{)}, agent aia_{i} needs to consider the worst-case realizations of the neighbouring inputs uj,jSiu_{j},~\forall j\in S_{i}, which result in the minimum value of the inter-update duration. This leads to the following max-min optimization:

ui=maxui𝒰,t~\displaystyle u^{\ast}_{i}=~\max\limits_{\begin{subarray}{c}u_{i}\in\mathcal{U},\\ \widetilde{t}\in\mathbb{R}\end{subarray}}~ minuj𝒰,jSit~tik\displaystyle\min\limits_{\begin{subarray}{c}u_{j}\in\mathcal{U},\\ \forall j\in S_{i}\end{subarray}}~~~\widetilde{t}-t^{k}_{i} (14)
tik+1=argmaxui𝒰,t~\displaystyle t^{k+1}_{i}=\arg~\max\limits_{\begin{subarray}{c}u_{i}\in\mathcal{U},\\ \widetilde{t}\in\mathbb{R}\end{subarray}}~ minuj𝒰,jSit~tik\displaystyle\min\limits_{\begin{subarray}{c}u_{j}\in\mathcal{U},\\ \forall j\in S_{i}\end{subarray}}~~~\widetilde{t}-t^{k}_{i} (15)
s.t.|zi(t)|α,t[tik,t~]\displaystyle~~~~\text{s.t.}~~~|z_{i}(t)|\leq\alpha,~~\forall t\in\big{[}t^{k}_{i},\widetilde{t}\hskip 2.84544pt\big{]} (16)

The following lemma shows that control law (7)-(8) is the solution of optimization (14)-(16).

Lemma 11.

Consider any agent aia_{i} in MAS (1). Let zi(tik)z_{i}\big{(}t^{k}_{i}\big{)} be the local disagreement of aia_{i} at its update instant tikt^{k}_{i} such that |zi(tik)|α\big{|}z_{i}\big{(}t^{k}_{i}\big{)}\big{|}\leq\alpha. Then, the control law (7)-(8) is the solution of optimization (14)-(16).

Proof.

See the Appendix for the proof. ∎

V-B Maximization of inter-update durations: For |zi(tik)|>α\big{|}z_{i}\big{(}t^{k}_{i}\big{)}\big{|}>\alpha

Let zi(tik)z_{i}\big{(}t^{k}_{i}\big{)} be the local disagreement of agent aia_{i} at its update instant tikt^{k}_{i} such that |zi(tik)|>α\big{|}z_{i}\big{(}t^{k}_{i}\big{)}\big{|}>\alpha. Then, in order to achieve α\alpha-consensus, it is necessary to steer |zi(tik)|\big{|}z_{i}\big{(}t^{k}_{i}\big{)}\big{|} below the consensus bound α\alpha. If that could be done in time-optimal manner, it is an added advantage. Recall that in Section II-F, we discussed the time-optimal consensus rule (5) which achieves conventional consensus of MAS (1) in minimum time. Motivated from this rule, we choose our control rule as

ui(t)=βsign(zi(tik)),t[tik,tik+1)u^{\ast}_{i}(t)=-\beta\hskip 2.84544ptsign\Big{(}z_{i}\big{(}t^{k}_{i}\big{)}\Big{)},~~~~~~~~\forall t\in\big{[}t^{k}_{i},t^{k+1}_{i}\big{)} (17)

As mentioned in Section V-A, the goal of agent aia_{i} is to maximize the inter-update duration (tik+1tik)\big{(}t^{k+1}_{i}-t^{k}_{i}\big{)} by delaying its next update instant tik+1t^{k+1}_{i}. However, while doing so, agent aia_{i} must satisfy the necessary condition (11).

Recall that |zi(tik)|>α\big{|}z_{i}\big{(}t^{k}_{i}\big{)}\big{|}>\alpha. Without loss of generality, assume that zi(tik)<αz_{i}\big{(}t^{k}_{i}\big{)}<-\alpha. Then, it follows from (17) that ui(t)=β,t[tik,tik+1)u^{\ast}_{i}(t)=\beta,~\forall t\in[t^{k}_{i},t^{k+1}_{i}). Let t^[tik,tik+1)\hat{t}\in[t^{k}_{i},t^{k+1}_{i}) be the first time instant at which zi(t^)=αz_{i}(\hat{t}~\!)=-\alpha. Then, in order to satisfy (11), agent aia_{i} needs to ensure that

zi(t)[α,α],t[t^,tik+1]z_{i}(t)\in[-\alpha,\alpha],~~~~~~~~\forall t\in[\hat{t},t^{k+1}_{i}]

which is equivalent to

zi(t)α,t[t^,tik+1]z_{i}(t)\leq\alpha,~~~~~~~~~\forall t\in[\hat{t},t^{k+1}_{i}]

As discussed in Section V-A, at instant tikt^{k}_{i}, agent aia_{i} does not know the future values of the neighbouring inputs in the interval [tik,tik+1)[t^{k}_{i},t^{k+1}_{i}). Thus, while maximizing the inter-update duration (tik+1tik)\big{(}t^{k+1}_{i}-t^{k}_{i}\big{)}, agent aia_{i} needs to consider the worst-case realizations of the neighbouring inputs uj,jSiu_{j},~\forall j\in S_{i}, which result in the minimum value of the inter-update duration. This leads to the following max-min optimization:

tik+1=argmaxt~\displaystyle t^{k+1}_{i}=\arg~\max\limits_{\widetilde{t}\in{\mathbb{R}}}~ minui=β,uj𝒰,jSit~tik\displaystyle\min\limits_{\begin{subarray}{c}u_{i}=\beta,\\ u_{j}\in\mathcal{U},\\ \forall j\in S_{i}\end{subarray}}~~~\widetilde{t}-t^{k}_{i} (18)
s.t.zi(t)α,t[tik,t~]\displaystyle~~~~\text{s.t.}~~~z_{i}(t)\leq\alpha,~~~\forall t\in\big{[}t^{k}_{i},\widetilde{t}\hskip 2.84544pt\big{]} (19)

The following lemma shows that control law (9)-(10) is the solution of optimization (18)-(19).

Lemma 12.

Consider any agent aia_{i} in MAS (1). Let zi(tik)z_{i}\big{(}t^{k}_{i}\big{)} be the local disagreement of aia_{i} at its update instant tikt^{k}_{i} such that zi(tik)<αz_{i}\big{(}t^{k}_{i}\big{)}<-\alpha. Then, the control law (9)-(10) is the solution of optimization (18)-(19).

Proof.

See the Appendix for the proof. ∎

Remark 13.

The optimization (18)-(19) and Lemma 12 correspond to the case zi(tik)<αz_{i}\big{(}t^{k}_{i}\big{)}<-\alpha. If zi(tik)>αz_{i}\big{(}t^{k}_{i}\big{)}>\alpha, the constraint (19) becomes zi(t)α,t[tik,t~]z_{i}(t)\geq-\alpha,~\forall t\in\big{[}t^{k}_{i},\widetilde{t}\hskip 2.27626pt\big{]}. Then, by following the arguments in the proof of Lemma 12, we can show that even for zi(tik)>αz_{i}\big{(}t^{k}_{i}\big{)}>\alpha, the control rule (9)-(10) is the solution of optimization (18)-(19).

V-C Feasibility of Protocol \Romannum1\Romannum{1}

In this section, we present the proof of Lemma 7 which claims that Protocol \Romannum1\Romannum{1} satisfies the necessary condition (11).

Proof.

of Lemma 7 : Consider MAS (1) with any nn, any connected communication graph GG and any initial condition X(0)nX(0)\in{\mathbb{R}}^{n}. Let aia_{i} be any agent in MAS (1) and t~i[tik,tik+1)\widetilde{t}_{i}\in[t^{k}_{i},t^{k+1}_{i}) be the first time instant under Protocol \Romannum1\Romannum{1} at which |zi(t~i)|α|z_{i}\big{(}\!\!~\widetilde{t}_{i}\big{)}|\leq\alpha. If |zi(tik)|α|z_{i}\big{(}t^{k}_{i}\big{)}|\leq\alpha, then it follows from (16) and Lemma 11 that |zi(t)|α,t[t~i,tik+1]|z_{i}(t)|\leq\alpha,~\forall t\in\big{[}~\!\widetilde{t}_{i},t^{k+1}_{i}\big{]} under Protocol \Romannum1\Romannum{1}. On the other hand, if |zi(tik)|>α|z_{i}\big{(}t^{k}_{i}\big{)}|>\alpha, then it follows from (19), Lemma 12 and Remark 13 that |zi(t)|α,t[t~i,tik+1]|z_{i}(t)|\leq\alpha,~\forall t\in\big{[}~\!\widetilde{t}_{i},t^{k+1}_{i}\big{]} under Protocol \Romannum1\Romannum{1}. Consequently, in both cases, Lemma 11 leads to |zi(t)|α,t>tik+1|z_{i}(t)|\leq\alpha,~\forall t>t^{k+1}_{i}. This proves that Protocol \Romannum1\Romannum{1} satisfies the necessary condition (11). ∎

By using Lemma 7, it is already shown in Theorem 8 that Protocol \Romannum1\Romannum{1} satisfies constraint (4) in Problem 3 for time 2T(X(0))2T^{\ast}\big{(}X(0)\big{)}, i.e., achieves α\alpha consensus of MAS (1) within time 2T(X(0))2T^{\ast}\big{(}X(0)\big{)}, for every nn, every connected GG and every X(0)nX(0)\in{\mathbb{R}}^{n}. In the next section, we prove the optimality of Protocol \Romannum1\Romannum{1}.

V-D Proof of optimality of Protocol \Romannum1\Romannum{1}

Consider MAS (1) with initial condition X(0)nX(0)\in{\mathbb{R}}^{n}. Let QQ be any protocol which is a feasible solution of Problem 3 for T(X(0))=2T(X(0))T\big{(}X(0)\big{)}=2T^{\ast}\big{(}X(0)\big{)}, i.e., a discrete communication based distributed protocol which achieves α\alpha-consensus of MAS (1) within time 2T(X(0))2T^{\ast}\big{(}X(0)\big{)}, for every nn, every connected GG and every X(0)nX(0)\in{\mathbb{R}}^{n}. Let CMASQ(2T(X(0)),X(0))C^{Q}_{MAS}\Big{(}2T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)} and CMAS(2T(X(0)),X(0))C^{\ast}_{MAS}\Big{(}2T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)} denote the value of the aggregate communication cost CMASC_{MAS} defined in (3), under Protocol QQ and Protocol \Romannum1\Romannum{1}, respectively.

Theorem 14.

Consider MAS (1) with initial condition X(0)nX(0)\in{\mathbb{R}}^{n}. Then, under Protocol \Romannum1\Romannum{1}, the following holds.

CMAS(2T(X(0)),X(0))CMASQ(2T(X(0)),X(0))C^{\ast}_{MAS}\Big{(}2T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)}\leq C^{Q}_{MAS}\Big{(}2T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)} (20)
Proof.

For the sake of contradiction, assume that

CMASQ(2T(X(0)),X(0))<CMAS(2T(X(0)),X(0))C^{Q}_{MAS}\Big{(}2T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)}<C^{\ast}_{MAS}\Big{(}2T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)} (21)

Recall that Ci(t,X(0))C_{i}\big{(}t,X(0)\big{)} denotes the number of update instants of agent aia_{i} in the interval [0,t)[0,t), corresponding to initial condition X(0)X(0). Let CiQ(2T(X(0)),X(0))C^{Q}_{i}\Big{(}2T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)} and Ci(2T(X(0)),X(0))C^{\ast}_{i}\Big{(}2T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)} denote the value of Ci(2T(X(0),X(0)))C_{i}\Big{(}2T^{\ast}\big{(}X(0),X(0)\big{)}\Big{)}, under Protocol QQ and Protocol \Romannum1\Romannum{1}, respectively. Then, (3) and (21) imply that there exists at least one agent, say aia_{i}, such that

CiQ(2T(X(0)),X(0))<Ci(2T(X(0)),X(0))C^{Q}_{i}\Big{(}2T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)}<C^{\ast}_{i}\Big{(}2T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)}

This implies that under Protocol QQ, at least one inter-update duration of agent aia_{i} is longer than that prescribed by control laws (7) and (9) in Protocol \Romannum1\Romannum{1}.

Recall from (15) and Lemma 11 that control law (7) gives the maximum inter-update duration under constraint (16). Similarly, it follows from (18) and Lemma 12 that control law (9) gives the maximum inter-update duration under constraint (19). Then, as one inter-update duration under Protocol QQ is longer than that prescribed by control laws (7) and (9), Protocol QQ must be violating either constraint (16) or constraint (19). Recall that violation of (16) or (19) by Protocol QQ results in the violation of necessary condition (11) on feasible protocols. This contradicts the fact that Protocol QQ is a feasible solution of Problem 3 and proves claim (20). ∎

VI Protocol \Romannum2\Romannum{2}: For T(X(0))>2T(X(0))T\big{(}X(0)\big{)}>2T^{\ast}\big{(}X(0)\big{)}

In Section V-D, we proved that Protocol \Romannum1\Romannum{1} is the solution of Problem 3 for T(X(0))=2T(X(0))T\big{(}X(0)\big{)}=2T^{\ast}\big{(}X(0)\big{)}. In this section, we extend Protocol \Romannum1\Romannum{1} for T(X(0))>2T(X(0))T\big{(}X(0)\big{)}>2T^{\ast}\big{(}X(0)\big{)}. We refer to the extended protocol as Protocol \Romannum2\Romannum{2}.

VI-A Protocol \Romannum2\Romannum{2}

Let γ>1\gamma>1 be a real number and T(X(0))=2γT(X(0))T\big{(}X(0)\big{)}=2\gamma T^{\ast}\big{(}X(0)\big{)} be the pre-specified α\alpha-consensus time. Define β~:=βγ\widetilde{\beta}:=\dfrac{\beta}{\gamma}. Then, Protocol \Romannum2\Romannum{2} is same as Protocol \Romannum1\Romannum{1}, except the control bound β~\widetilde{\beta} in place of β\beta.

VI-B Optimality of Protocol \Romannum2\Romannum{2}

In this section, we first show that Protocol \Romannum2\Romannum{2} achieves α\alpha-consensus of MAS (1) within the pre-specified time T(X(0))=2γT(X(0))T\big{(}X(0)\big{)}=2\gamma T^{\ast}\big{(}X(0)\big{)}.

Theorem 15.

Consider MAS (1). Let α+\alpha\in{\mathbb{R}}^{+} be the pre-specified consensus bound. Then, for every nn, every connected communication graph GG and every X(0)nX(0)\in{\mathbb{R}}^{n}, Protocol \Romannum2\Romannum{2} achieves α\alpha-consensus of MAS (1) in time less than or equal to 2γT(X(0))2\gamma T^{\ast}\big{(}X(0)\big{)}. Moreover, there exist n~\widetilde{n}, connected graph G~\widetilde{G} and initial condition X~(0)n~\widetilde{X}(0)\in{\mathbb{R}}^{\widetilde{n}} for which α\alpha-consensus time under Protocol \Romannum2\Romannum{2} is arbitrarily close to 2γT(X(0))2\gamma T^{\ast}\big{(}X(0)\big{)}.

Proof.

Recall that the control bounds in Protocols \Romannum1\Romannum{1} and \Romannum2\Romannum{2} are β\beta and β~=βγ\widetilde{\beta}=\dfrac{\beta}{\gamma}, respectively. Thus, the dynamics of MAS (1) under Protocol \Romannum2\Romannum{2} is γ\gamma times slower than that of under Protocol \Romannum1\Romannum{1}. Then, the claim follows from the arguments in the proofs of Theorems 8 and 10. ∎

Now, we prove the optimality of Protocol \Romannum2\Romannum{2} for T(X(0))=2γT(X(0))T\big{(}X(0)\big{)}=2\gamma T^{\ast}\big{(}X(0)\big{)}. Consider MAS (1) with an initial condition X(0)nX(0)\in{\mathbb{R}}^{n}. Let QQ be any protocol which is a feasible solution of Problem 3 for T(X(0))=2γT(X(0))T\big{(}X(0)\big{)}=2\gamma T^{\ast}\big{(}X(0)\big{)}. Let CMASQ(2γT(X(0)),X(0))C^{Q}_{MAS}\Big{(}2\gamma T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)} and CMAS(2γT(X(0)),X(0))C^{\ast}_{MAS}\Big{(}2\gamma T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)} denote the value of CMASC_{MAS} defined in (3), under Protocol QQ and Protocol \Romannum2\Romannum{2}, respectively.

Theorem 16.

Consider MAS (1) with initial condition X(0)nX(0)\in{\mathbb{R}}^{n}. Then, under Protocol \Romannum2\Romannum{2}, the following holds.

CMAS(2γT(X(0)),X(0))CMASQ(2γT(X(0)),X(0))C^{\ast}_{MAS}\Big{(}2\gamma T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)}\leq C^{Q}_{MAS}\Big{(}2\gamma T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)}
Proof.

The claim follows from the arguments in the proof of Theorem 14. ∎

Next, we analyze the effect of T(X(0))T\big{(}X(0)\big{)} on CMASC^{\ast}_{MAS}. Intuitively, with the increase in T(X(0))T\big{(}X(0)\big{)}, agents in MAS (1) can afford to delay their next update instants, which would result in lower CMASC^{\ast}_{MAS}. However, the following theorem shows that the above intuition is not true after T(X(0))=2T(X(0))T\big{(}X(0)\big{)}=2T^{\ast}\big{(}X(0)\big{)}.

Theorem 17.

Consider MAS (1) with initial condition X(0)nX(0)\in{\mathbb{R}}^{n}. Then, under Protocol \Romannum2\Romannum{2}, for every real γ>1\gamma>1, the following holds.

CMAS(2γT(X(0)),X(0))=CMAS(2T(X(0)),X(0))C^{\ast}_{MAS}\Big{(}2\gamma T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)}=C^{\ast}_{MAS}\Big{(}2T^{\ast}\big{(}X(0)\big{)},X(0)\Big{)} (22)
Proof.

Let tik,βt^{k,\!~\beta}_{i} and tik,β~t^{k,\!~\widetilde{\beta}}_{i} denote the kkth update instants of agent aia_{i} under Protocols \Romannum1\Romannum{1} and \Romannum2\Romannum{2}, respectively. Recall that the dynamics of MAS (1) under Protocol \Romannum2\Romannum{2} is γ\gamma times slower than that of under Protocol \Romannum1\Romannum{1}. Then, it is easy to show that

tik,β~=γtik,β,iP,k1t^{k,\!~\widetilde{\beta}}_{i}=\gamma~\!t^{k,\!~\beta}_{i},~~~~~~~~\forall i\in P,~~~~~~~~~\forall k\geq 1

Hence, the timeline of agent aia_{i} under Protocol \Romannum2\Romannum{2} is just the γ\gamma-scaled version of its timeline under Protocol \Romannum1\Romannum{1}. As a result, the number of update instants of agent aia_{i} in the interval [0,2γT(X(0))]\big{[}0,\!~2\gamma T^{\ast}\big{(}X(0)\big{)}\big{]} under Protocol \Romannum2\Romannum{2} is equal to that of in the interval [0,2T(X(0))]\big{[}0,\!~2T^{\ast}\big{(}X(0)\big{)}\big{]} under Protocol \Romannum1\Romannum{1}. Then, the claim (22) follows from the definition of CMASC_{MAS} in (3). ∎

VII Simulation results

In this section, we present the simulation results obtained under Protocol \Romannum2\Romannum{2} for various T(X(0))T\big{(}X(0)\big{)}’s. Consider MAS (1) with n=6n=6 agents and the communication graph shown in Fig. 4.

112233445566
Figure 4: Communication graph for simulation

Let the consensus bound and the control bound be α=0.6\alpha=0.6 and β=1\beta=1, respectively. Let the initial condition of the MAS be X(0)=[7,2,4,3,1,5]X(0)=[7,2,4,3,1,5]. Then, it follows from the definition of T(X(0))T^{\ast}\big{(}X(0)\big{)} in (6) that T(X(0))=3T^{\ast}\big{(}X(0)\big{)}=3 sec.

We simulated the above MAS under Protocol \Romannum2\Romannum{2} for T(X(0))=2T(X(0))T\big{(}X(0)\big{)}=2T^{\ast}\big{(}X(0)\big{)}, T(X(0))=10T(X(0))T\big{(}X(0)\big{)}=10T^{\ast}\big{(}X(0)\big{)} and T(X(0))=20T(X(0))T\big{(}X(0)\big{)}=20T^{\ast}\big{(}X(0)\big{)}. The corresponding disagreement trajectories ziz_{i}’s are plotted in Fig. 5(a), 5(b) and 5(c), respectively. The corresponding α\alpha-consensus times and optimal communication costs are given in Table \Romannum1\Romannum{1}.

Notice that in each of the above three cases, Protocol \Romannum2\Romannum{2} achieves α\alpha-consensus within the prespecified time T(X(0))T\big{(}X(0)\big{)}, as proved in Theorem 15. Further, the α\alpha-consensus times corresponding to T(X(0))=10T(X(0))T\big{(}X(0)\big{)}=10T^{\ast}\big{(}X(0)\big{)} and T(X(0))=20T(X(0))T\big{(}X(0)\big{)}=20T^{\ast}\big{(}X(0)\big{)} are 55 and 1010 times of that corresponding to T(X(0))=2T(X(0))T\big{(}X(0)\big{)}=2T^{\ast}\big{(}X(0)\big{)}, respectively. This observation is in accordance with the arguments in the proof of Theorem 15. In addition, observe that the trajectories in Figures 5(b) and 5(c) are just the time-stretched versions of the corresponding trajectories in Fig. 5(a). The optimal communication costs CiC^{\ast}_{i}’s and CMASC^{\ast}_{MAS} corresponding to T(X(0))=2T(X(0))T\big{(}X(0)\big{)}=2T^{\ast}\big{(}X(0)\big{)}, T(X(0))=10T(X(0))T\big{(}X(0)\big{)}=10T^{\ast}\big{(}X(0)\big{)} and T(X(0))=20T(X(0))T\big{(}X(0)\big{)}=20T^{\ast}\big{(}X(0)\big{)} are same, as proved in Theorem 17.

Refer to caption
(a) T(X(0))=2T(X(0))T\big{(}X(0)\big{)}=2T^{\ast}\big{(}X(0)\big{)}
Refer to caption
(b) T(X(0))=10T(X(0))T\big{(}X(0)\big{)}=10T^{\ast}\big{(}X(0)\big{)}
Refer to caption
(c) T(X(0))=20T(X(0))T\big{(}X(0)\big{)}=20T^{\ast}\big{(}X(0)\big{)}
Figure 5: Comparison of local disagreement trajectories under Protocol \Romannum2\Romannum{2} for various T(X(0))T\big{(}X(0)\big{)}’s
TABLE I: Consensus time and communication cost under Protocol \Romannum2\Romannum{2} for various T(X(0))T\big{(}X(0)\big{)}’s
T(X(0))T\big{(}X(0)\big{)} α\alpha-consensus time [sec] C1C^{\ast}_{1} C2C^{\ast}_{2} C3C^{\ast}_{3} C4C^{\ast}_{4} C5C^{\ast}_{5} C6C^{\ast}_{6} CMASC^{\ast}_{MAS}
2T(X(0))2T^{\ast}\big{(}X(0)\big{)} 2.262.26 88 99 3030 3030 99 99 9595
10T(X(0))10T^{\ast}\big{(}X(0)\big{)} 11.311.3 88 99 3030 3030 99 99 9595
20T(X(0))20T^{\ast}\big{(}X(0)\big{)} 22.622.6 88 99 3030 3030 99 99 9595

VIII Conclusion and future work

In this paper, a distributed consensus protocol (Protocol \Romannum2\Romannum{2}) is proposed for a MAS of single-integrator agents with bounded inputs and time-invariant communication graph. We showed (Theorems 15 and 16) that the proposed protocol minimizes the aggregate number of communication instants required to achieve α\alpha-consensus within the pre-specified time T(X(0))2T(X(0))T\big{(}X(0)\big{)}\geq 2T^{\ast}\big{(}X(0)\big{)}. The control rules in the proposed protocol are obtained by maximizing the inter-update durations. We computed (Lemmas 11 and 12) the closed form solutions of these optimizations. Finally, we presented the simulation results which verify the theoretical claims. The work of extending the proposed protocol to MASs with complex dynamics is in progress.

Appendix

A. Intermediate results for the proof of Theorem 8

In this section, we develop a few intermediate results which will be used later in the proof of Theorem 8. Let xi(tik)x_{i}\big{(}t^{k}_{i}\big{)} and zi(tik)z_{i}\big{(}t^{k}_{i}\big{)} be the state and the local disagreement of agent aia_{i} respectively, at its update instant tikt^{k}_{i}. Recall that SiS_{i} is the set of indices of neighbours of agent aia_{i}, with cardinality nin_{i}.

Lemma 18.

Consider agent aia_{i} in MAS (1) with |zi(tik)|α|z_{i}\big{(}t^{k}_{i}\big{)}|\leq\alpha. Then, under Protocol \Romannum1\Romannum{1}, the following holds.
11xi(tik+1)=jSixj(tik)nix_{i}\big{(}t^{k+1}_{i}\big{)}=\dfrac{\sum_{j\in S_{i}}x_{j}\big{(}t^{k}_{i}\big{)}}{n_{i}}
22|zi(t)|α,ttik|z_{i}(t)|\leq\alpha,\hskip 9.98685pt\forall t\geq t^{k}_{i}

Proof.

At instant t>tikt>t^{k}_{i}, the state xix_{i} of agent aia_{i} is xi(t)=xi(tik)+tiktui(τ)𝑑τx_{i}(t)=x_{i}(t^{k}_{i})+\int_{t^{k}_{i}}^{t}u_{i}(\tau)d\tau. It is given that |zi(tik)|α|z_{i}(t^{k}_{i})|\leq\alpha. Then, control rule (7)-(8) in Protocol \Romannum1\Romannum{1} lead to

xi(tik+1)=xi(tik)+tiktik+αβnizi(tik)αβ𝑑τ=xi(tik)zi(tik)nix_{i}\big{(}t^{k+1}_{i}\big{)}=x_{i}\big{(}t^{k}_{i}\big{)}+\!\!\int_{t^{k}_{i}}^{t^{k}_{i}+\dfrac{\alpha}{\beta n_{i}}}\frac{-z_{i}(t^{k}_{i})}{\alpha}\beta d\tau=x_{i}\big{(}t^{k}_{i}\big{)}-\dfrac{z_{i}\big{(}t^{k}_{i}\big{)}}{n_{i}}

Then, it follows from the definition of ziz_{i} in (2) that

xi(tik+1)=xi(tik)jSi(xi(tik)xj(tik))nix_{i}\big{(}t^{k+1}_{i}\big{)}=x_{i}\big{(}t^{k}_{i}\big{)}-\dfrac{\sum_{j\in S_{i}}\Big{(}x_{i}\big{(}t^{k}_{i}\big{)}-x_{j}\big{(}t^{k}_{i}\big{)}\Big{)}}{n_{i}} (23)

Recall that the cardinality of SiS_{i} is nin_{i}. Then, by re-arranging the terms in (23), we get xi(tik+1)=jSixj(tik)nix_{i}\big{(}t^{k+1}_{i}\big{)}=\dfrac{\sum_{j\in S_{i}}x_{j}\big{(}t^{k}_{i}\big{)}}{n_{i}}. This proves the claim in Lemma 18.11.

Next, we prove that under Protocol \Romannum1\Romannum{1}, ziz_{i} satisfies |zi(t)|α,t[tik,tik+1]|z_{i}(t)|\leq\alpha,~\forall t\in[t^{k}_{i},t^{k+1}_{i}]. Then, the claim in Lemma 18.22 follows by repeating arguments in the subsequent inter-update intervals.

Let the control input uiu^{\ast}_{i} be as defined in (8). We know from the dynamics of xix_{i} in (1) and the definition of ziz_{i} in (2) that

zi(tik+1)=zi(tik)+nitiktik+1ui(t)𝑑tjSi(tiktik+1uj(t)𝑑t)z_{i}\big{(}t^{k+1}_{i}\big{)}=z_{i}\big{(}t^{k}_{i}\big{)}+n_{i}\int_{t^{k}_{i}}^{t^{k+1}_{i}}\!\!\!\!\!\!\!\!u^{\ast}_{i}(t)~dt-\sum_{j\in S_{i}}\Bigg{(}\int_{t^{k}_{i}}^{t^{k+1}_{i}}\!\!\!\!\!\!\!\!u_{j}(t)~dt\Bigg{)} (24)

By using (7) and (8), it is easy to show that

zi(tik)+nitiktik+1ui(t)𝑑t=0z_{i}\big{(}t^{k}_{i}\big{)}+n_{i}\int_{t^{k}_{i}}^{t^{k+1}_{i}}\!\!\!\!\!\!\!\!u^{\ast}_{i}(t)~dt=0 (25)

Further, recall that |uj(t)|β,t0,jP|u_{j}(t)|\leq\beta,~\forall t\geq 0,~\forall j\in P. Then, it follows from (7) that

|tiktik+1uj(t)𝑑t|=|tiktik+αβniuj(t)𝑑t|αni\Bigg{|}\int_{t^{k}_{i}}^{t^{k+1}_{i}}\!\!\!\!\!\!\!\!u_{j}(t)~dt\Bigg{|}=\Bigg{|}\int_{t^{k}_{i}}^{t^{k}_{i}+\dfrac{\alpha}{\beta n_{i}}}u_{j}(t)~dt\Bigg{|}\leq\frac{\alpha}{n_{i}} (26)

Then, (24), (25) and (26) together lead to

|zi(tik+1)|α\big{|}z_{i}\big{(}t^{k+1}_{i}\big{)}\big{|}\leq\alpha (27)

Now, we claim that |zi(t)|α,t[tik,tik+1]|z_{i}(t)|\leq\alpha,~\forall t\in[t^{k}_{i},t^{k+1}_{i}]. Assume for the sake of contradiction that there exists t^(tik,tik+1)\hat{t}\in(t^{k}_{i},t^{k+1}_{i}) such that |zi(t^)|>α|z_{i}\big{(}\hat{t}\hskip 1.13791pt\big{)}|>\alpha. Take uj(t)=ui(t),t[t^,tik+1),jSiu_{j}(t)=u^{\ast}_{i}(t),~\forall t\in\big{[}\hat{t},t^{k+1}_{i}\big{)},~\forall j\in S_{i}. Then, it follows from (24) that |zi(tik+1)|=|zi(t^)|>α|z_{i}\big{(}t^{k+1}_{i}\big{)}|=|z_{i}\big{(}\hat{t}\hskip 1.13791pt\big{)}|>\alpha, which contradicts (27) and proves our claim that |zi(t)|α,t[tik,tik+1]|z_{i}(t)|\leq\alpha,~\forall t\in[t^{k}_{i},t^{k+1}_{i}]. By repeating the above arguments in the subsequent inter-update intervals of agent aia_{i}, we get Lemma 18.22. This completes the proof. ∎

Let tikt^{k}_{i} be an update instant of agent aia_{i} in MAS (1). Define

ximin(tik)\displaystyle x_{i}^{min}\big{(}t^{k}_{i}\big{)} :=minj(Sii)xj(tik)\displaystyle:=\min_{j\in(S_{i}\cup\hskip 0.85355pti)}~x_{j}\big{(}t^{k}_{i}\big{)} (28)
ximax(tik)\displaystyle x_{i}^{max}\big{(}t^{k}_{i}\big{)} :=maxj(Sii)xj(tik)\displaystyle:=\max_{j\in(S_{i}\cup\hskip 0.85355pti)}~x_{j}\big{(}t^{k}_{i}\big{)} (29)

The following lemma shows that the state xix_{i} remains bounded between ximin(tik)x_{i}^{min}\big{(}t^{k}_{i}\big{)} and ximax(tik)x_{i}^{max}\big{(}t^{k}_{i}\big{)}, in the interval [tik,tik+1][t^{k}_{i},t^{k+1}_{i}].

Lemma 19.

Consider any agent aia_{i} in MAS (1). Under Protocol \Romannum1\Romannum{1}, the state of agent aia_{i} satisfies the following.

ximin(tik)xi(t)ximax(tik),t[tik,tik+1]x_{i}^{min}\big{(}t^{k}_{i}\big{)}\leq x_{i}(t)\leq x_{i}^{max}\big{(}t^{k}_{i}\big{)},~~~~~\forall t\in\big{[}t^{k}_{i},t^{k+1}_{i}\big{]} (30)
Proof.

Based on the value of zi(tik)z_{i}\big{(}t^{k}_{i}\big{)}, there are two cases, |zi(tik)|α|z_{i}\big{(}t^{k}_{i}\big{)}|\leq\alpha and |zi(tik)|>α|z_{i}\big{(}t^{k}_{i}\big{)}|>\alpha. We will prove the claim for the first case. The proof for the second case follows analogously.

Without loss of generality, assume that αzi(tik)0-\alpha\leq z_{i}\big{(}t^{k}_{i}\big{)}\leq 0. Then, it follows from (8) that ui(t)=zi(tik)αβ0,t[tik,tik+1)u_{i}(t)=\dfrac{-z_{i}\big{(}t^{k}_{i}\big{)}}{\alpha}\beta\geq 0,~\forall t\in[t^{k}_{i},t^{k+1}_{i}). Hence, x˙i(t)=ui(t)0,t[tik,tik+1)\dot{x}_{i}(t)=u_{i}(t)\geq 0,~\forall t\in[t^{k}_{i},t^{k+1}_{i}), which implies that xi(tik)xi(t),t[tik,tik+1]x_{i}\big{(}t^{k}_{i}\big{)}\leq x_{i}(t),~\forall t\in[t^{k}_{i},t^{k+1}_{i}]. Then, the definition of ximin(tik)x_{i}^{min}\big{(}t^{k}_{i}\big{)} in (28) leads to

ximin(tik)xi(tik)xi(t),t[tik,tik+1]x_{i}^{min}\big{(}t^{k}_{i}\big{)}\leq x_{i}\big{(}t^{k}_{i}\big{)}\leq x_{i}(t),~~~~~~\forall t\in\big{[}t^{k}_{i},t^{k+1}_{i}\big{]} (31)

Next, let q(Sii)q\in(S_{i}\cup i) be such that xq(tik)=ximax(tik)x_{q}\big{(}t^{k}_{i}\big{)}=x_{i}^{max}\big{(}t^{k}_{i}\big{)}. Then, based on the value of qq, there are following two cases:
Case 1 : q=iq=i
In this case, we must have xi(tik)=xj(tik),jSix_{i}\big{(}t^{k}_{i}\big{)}=x_{j}\big{(}t^{k}_{i}\big{)},~\forall j\in S_{i}. Otherwise, we will have zi(tik)>0z_{i}\big{(}t^{k}_{i}\big{)}>0, which violates the assumption αzi(tik)0-\alpha\leq z_{i}\big{(}t^{k}_{i}\big{)}\leq 0. Thus, zi(tik)=0z_{i}\big{(}t^{k}_{i}\big{)}=0. Then, (8) leads to ui(t)=0,t[tik,tik+1)u_{i}(t)=0,~\forall t\in[t^{k}_{i},t^{k+1}_{i}) and hence, xi(tik+1)=xi(tik)=ximax(tik)x_{i}\big{(}t^{k+1}_{i}\big{)}=x_{i}\big{(}t^{k}_{i}\big{)}=x_{i}^{max}\big{(}t^{k}_{i}\big{)}.
Case 2 : qSiq\in S_{i}
Recall from Lemma 18.11 that xi(tik+1)x_{i}\big{(}t^{k+1}_{i}\big{)} is the average of xj(tik)x_{j}\big{(}t^{k}_{i}\big{)}’s for all jSij\in S_{i}. Thus, xi(tik+1)maxjSixj(tik)=ximax(tik)x_{i}\big{(}t^{k+1}_{i}\big{)}\leq\max_{j\in S_{i}}x_{j}\big{(}t^{k}_{i}\big{)}=x_{i}^{max}\big{(}t^{k}_{i}\big{)}.

In both of the above cases, we have xi(tik+1)ximax(tik)x_{i}\big{(}t^{k+1}_{i}\big{)}\leq x_{i}^{max}\big{(}t^{k}_{i}\big{)}. Recall that x˙i(t)0,t[tik,tk+1]\dot{x}_{i}(t)\geq 0,~\forall t\in[t^{k}_{i},t^{k+1}]. Thus,

xi(t)xi(tik+1)ximax(tik),t[tik,tik+1]x_{i}(t)\leq x_{i}\big{(}t^{k+1}_{i}\big{)}\leq x^{max}_{i}\big{(}t^{k}_{i}\big{)},~~~~~\forall t\in\big{[}t^{k}_{i},t^{k+1}_{i}\big{]} (32)

Then, the claim (30) follows from (31) and (32). ∎

Let xmin(0)x^{min}(0) and xmax(0)x^{max}(0) be as in the definition of T(X(0))T^{\ast}\big{(}X(0)\big{)} in (6). The following lemma shows that the states of all agents in MAS (1) remain bounded between xmin(0)x^{min}(0) and xmax(0)x^{max}(0), for all t0t\geq 0.

Lemma 20.

Consider MAS (1) with initial condition X(0)nX(0)\in{\mathbb{R}}^{n}. Under Protocol \Romannum1\Romannum{1}, the following holds.

xmin(0)xi(t)xmax(0),t0,iPx^{min}(0)\leq x_{i}(t)\leq x^{max}(0),~~~~~\forall t\geq 0,~~~~~\forall i\in P (33)
Proof.

Recall that according to Protocol \Romannum1\Romannum{1}, the time instant ti1=0t^{1}_{i}=0 is the first update instant of every agent ai,iPa_{i},~i\in P. Let ximin(tik)x_{i}^{min}\big{(}t^{k}_{i}\big{)} and ximax(tik)x_{i}^{max}\big{(}t^{k}_{i}\big{)} be as defined in (28) and (29), respectively. Then, it follows from Lemma 19 that

ximin(0)xi(t)ximax(0),t[0,ti2],iPx_{i}^{min}(0)\leq x_{i}(t)\leq x_{i}^{max}(0),~~~~\forall t\in[0,t^{2}_{i}],~~~\forall i\in P (34)

Further, it is clear from the definitions of ximinx_{i}^{min} and xmin(0)x^{min}(0) that xmin(0)ximin(0),iPx^{min}(0)\leq x_{i}^{min}(0),~\forall i\in P. Similarly, ximax(0)xmax(0),iPx_{i}^{max}(0)\leq x^{max}(0),~\forall i\in P. Then, (34) leads to

xmin(0)xi(t)xmax(0),t[0,ti2],iPx^{min}(0)\leq x_{i}(t)\leq x^{max}(0),~~~~\forall t\in[0,t^{2}_{i}],\hskip 14.22636pt\forall i\in P

Now, by applying the above arguments in the subsequent inter-update intervals [ti2,ti3)[t^{2}_{i},t^{3}_{i}), [ti3,ti4),[t^{3}_{i},t^{4}_{i}),\dots of all agents, we get (33). This completes the proof. ∎

B. Proof of Theorem 8

For the sake of contradiction, assume that the claim is not correct, i.e., for some nn, GG and X(0)X(0), Protocol \Romannum1\Romannum{1} does not achieve α\alpha-consensus of MAS (1) within time 2T(X(0))2T^{\ast}\big{(}X(0)\big{)}. This implies that there exists at least one agent, say aia_{i}, and a time instant t^i>2T(X(0))\hat{t}_{i}>2T^{\ast}\big{(}X(0)\big{)}, such that |zi(t^i)|>α|z_{i}(\hat{t}_{i})|>\alpha. Recall from Lemma 7 that Protocol \Romannum1\Romannum{1} satisfies the necessary condition (11). Then, it follows from (11) that under Protocol \Romannum1\Romannum{1}, the following holds.

zi(t)<α,t[0,t^i]orzi(t)>α,t[0,t^i]z_{i}(t)<-\alpha,~\forall t\in[0,\hat{t}_{i}]\hskip 17.07182ptor\hskip 17.07182ptz_{i}(t)>\alpha,~\forall t\in[0,\hat{t}_{i}]

Without loss of generality, assume that zi(t)<α,t[0,t^i]z_{i}(t)<-\alpha,~\forall t\in[0,\hat{t}_{i}]. Then, it follows from control rule (10) that ui(t)=β,t[0,t^i]u_{i}(t)=\beta,~\forall t\in[0,\hat{t}_{i}] and hence, xi(t^i)=xi(0)+βt^ix_{i}\big{(}\hat{t}_{i}\big{)}=x_{i}(0)+\beta\hat{t}_{i}. As t^i>2T(X(0))\hat{t}_{i}>2T^{\ast}\big{(}X(0)\big{)}, we get

xi(t^i)>(xi(0)+2βT(X(0)))x_{i}\big{(}\hat{t}_{i}\big{)}>\Big{(}x_{i}(0)+2\beta T^{\ast}\big{(}X(0)\big{)}\Big{)} (35)

Let xmin(0)x^{min}(0) and xmax(0)x^{max}(0) be as in the definition of T(X(0))T^{\ast}\big{(}X(0)\big{)} in (6). Then, (6) and (35) lead to

xi(t^i)>(xi(0)+xmax(0)xmin(0))x_{i}\big{(}\hat{t}_{i}\big{)}>\Big{(}x_{i}(0)+x^{max}(0)-x^{min}(0)\Big{)} (36)

It is clear from the definition of xmin(0)x^{min}(0) that xi(0)xmin(0)x_{i}(0)\geq x^{min}(0). Then, (36) leads to xi(t^i)>xmax(0)x_{i}(\hat{t}_{i})>x^{max}(0). This contradicts Lemma 20 and proves the claim in Theorem 8.

C. Intermediate result for the proof of Theorem 10

Lemma 21.

Consider the MAS in Example 9 with the communication graph GeG_{e} shown in Fig. 3. Let rr be the number of neighbours of agent a1a_{1}. Then, for every μ+\mu\in{\mathbb{R}}^{+}, there exists an integer n=nμ,rn=n_{\mu,r} such that under Protocol \Romannum1\Romannum{1}, the following holds.

5μxi(t)5,tμ,i=2,3,,nμ,r5-\mu\leq x_{i}(t)\leq 5,~~~~~~~\forall t\geq\mu,~~~~~~~~\forall i=2,3,\dots,n_{\mu,r} (37)
Proof.

We first define the integer nμ,rn_{\mu,r} for which the claim (37) holds. Let μ+\mu\in{\mathbb{R}}^{+} be the given number. Define μ~:=μ2\widetilde{\mu}:=\dfrac{\mu}{2}. Let ceil be the standard ceiling function. Define the integers nμ1:=1+ceil(4μ~)n_{{\mu}_{1}}:=1+ceil\bigg{(}\dfrac{4}{\widetilde{\mu}}\bigg{)}, nμ2:=ceil(5μμ~)n_{{\mu}_{2}}:=ceil\bigg{(}\dfrac{5}{\mu-\widetilde{\mu}}\bigg{)}, nr1:=1+(16r5r+3)n_{r_{1}}:=1+\bigg{(}\dfrac{16r}{5r+3}\bigg{)} and nr2:=3r7n_{r_{2}}:=3r-7. Then, define

nμ,r:=max{8,nμ1,nμ2,nr1,nr2}+1n_{{\mu},r}:=\max\big{\{}8,n_{{\mu}_{1}},n_{{\mu}_{2}},n_{r_{1}},n_{r_{2}}\big{\}}+1 (38)

Now, consider the MAS in Example 9 with n=nμ,rn=n_{{\mu},r}. We will show that the claim (37) holds for this MAS. For the sake of clarity, we divide the remaining proof into five parts as follows.
1)1) Computation of ti2t^{2}_{i} and xi(ti2)x_{i}\big{(}t^{2}_{i}\big{)} for i1i\geq 1

According to Protocol \Romannum1\Romannum{1}, the first update instant of every agent is ti1=0t^{1}_{i}=0. At this instant, the states of agents are x1(0)=0x_{1}(0)=0 and xi(0)=5,i=2,,nμ,rx_{i}(0)=5,~\forall i=2,\dots,n_{{\mu},r}. Then, the local disagreement of agent a1a_{1} at instant t=0t=0 is z1(0)=5rz_{1}(0)=-5r. Note that r1r\geq 1. Thus, z1(0)<3=αz_{1}(0)<-3=-\alpha. Then, it follows from control rule (9)-(10) that

t12=\displaystyle t^{2}_{1}=~\!\! |z1(0)|+α2βn1=5r+32r\displaystyle\dfrac{|z_{1}(0)|+\alpha}{2\beta n_{1}}=\dfrac{5r+3}{2r} (39)
u1(t)=\displaystyle u^{\ast}_{1}(t)=~\!\! β=1,t[0,t12)\displaystyle\beta=1,~~~~~~~\forall t\in\big{[}0,t^{2}_{1}\big{)} (40)

As a result, x1(t12)=x1(0)+t12=5r+32rx_{1}\big{(}t^{2}_{1}\big{)}=x_{1}(0)+t^{2}_{1}=\dfrac{5r+3}{2r}.

Now, consider any agent al,l=2,,r+1a_{l},~l=2,\dots,r+1. Note that nl=nμ,r1n_{l}=n_{{\mu},r}-1. Thus, the local disagreement of agent ala_{l} at instant t=0t=0 is zl(0)=(50)+(nμ,r2)(55)=5>αz_{l}(0)=(5-0)+(n_{{\mu},r}-2)(5-5)=5>\alpha. Then, it follows from control rule (9)-(10) that

tl2=\displaystyle t^{2}_{l}= |zl(0)|+α2βnl=4nμ,r1\displaystyle\dfrac{|z_{l}(0)|+\alpha}{2\beta n_{l}}=\dfrac{4}{n_{{\mu},r}-1} (41)
ul(t)=\displaystyle u^{\ast}_{l}(t)= β=1,t[0,tl2)\displaystyle-\beta=-1,~~~~~\forall t\in\big{[}0,t^{2}_{l}\big{)} (42)

As a result, xl(tl2)=xl(0)tl2=54nμ,r1x_{l}\big{(}t^{2}_{l}\big{)}=x_{l}(0)-t^{2}_{l}=5-\dfrac{4}{n_{{\mu},r}-1}. It follows from the definitions of nμ1n_{{\mu}_{1}} and nμ,rn_{{\mu},r} in (38) that 4nμ,r1<μ~\dfrac{4}{n_{{\mu},r}-1}<\widetilde{\mu}. Thus, xl(tl2)x_{l}\big{(}t^{2}_{l}\big{)} satisfies (5μ~)<xl(tl2)<5(5-\widetilde{\mu})<x_{l}\big{(}t^{2}_{l}\big{)}<5. Recall from (42) that x˙l(t)<0,t[0,tl2)\dot{x}_{l}(t)<0,~\forall t\in[0,t^{2}_{l}). Recall also that xl(0)=5x_{l}(0)=5. Hence,

(5μ~)<xl(t)<5,t[0,tl2],l=2,,r+1(5-\widetilde{\mu})<x_{l}(t)<5,~~\forall t\in\big{[}0,t^{2}_{l}\big{]},~~\forall l=2,\dots,r+1 (43)

Next, consider any agent aj,j=r+2,,nμ,ra_{j},~j=r+2,\dots,n_{{\mu},r}. Note that nj=nμ,r2n_{j}=n_{{\mu},r}-2. Thus, the local disagreement of agent aja_{j} at instant t=0t=0 is zj(0)=(55)(nμ,r2)=0z_{j}(0)=(5-5)(n_{{\mu},r}-2)=0. Then, it follows from control rule (7)-(8) that

tj2=\displaystyle t^{2}_{j}=~\!\! αβnj=3nμ,r2\displaystyle\dfrac{\alpha}{\beta n_{j}}=\dfrac{3}{n_{{\mu},r}-2} (44)
uj(t)=\displaystyle u^{\ast}_{j}(t)=~\!\! 0,t[0,tj2)\displaystyle 0,~~~~~~\forall t\in\big{[}0,t^{2}_{j}\big{)} (45)

As a result, the state xjx_{j} satisfies

xj(t)=xj(0)=5,t[0,tj2],j=r+2,,nμ,rx_{j}(t)=x_{j}(0)=5,~~~~\forall t\in\big{[}0,t^{2}_{j}\big{]},~~~~\forall j=r+2,\dots,n_{{\mu},r} (46)

22) Ordering of update instants ti2t^{2}_{i}’s

Let ti2t^{2}_{i}’s be as given in (41) and (44). Recall from (38) that nμ,r>8n_{{\mu},r}>8. Then, it is easy to show that 3nμ,r2<4nμ,r1\dfrac{3}{n_{{\mu},r}-2}<\dfrac{4}{n_{{\mu},r}-1}. Then, it follows from (41) and (44) that

tr+22=tr+32==tnμ,r2<t22=t32==tr+12t^{2}_{r+2}=t^{2}_{r+3}=\dots=t^{2}_{n_{{\mu},r}}<t^{2}_{2}=t^{2}_{3}=\dots=t^{2}_{r+1} (47)

Let t12t^{2}_{1} be as given in (39). It follows from the definitions of nr1n_{r_{1}} and nμ,rn_{{\mu},r} in (38) that 4nμ,r1<5r+32r\dfrac{4}{n_{{\mu},r}-1}<\dfrac{5r+3}{2r}. Then, (39) and (41) lead to

t22=t32==tr+12<t12t^{2}_{2}=t^{2}_{3}=\dots=t^{2}_{r+1}<t^{2}_{1} (48)

33) Computation of ti3t^{3}_{i} for i2i\geq 2

Consider any agent al,l=2,,r+1a_{l},~l=2,\dots,r+1. Recall from (47) and (48) that tr+22<tl2<t12t^{2}_{r+2}<t^{2}_{l}<t^{2}_{1}. Now, we compute zl(tr+22)z_{l}\big{(}t^{2}_{r+2}\big{)}. According to (40) and (44), the state of agent a1a_{1} at instant tr+22t^{2}_{r+2} is x1(tr+22)=3nμ,r2x_{1}\big{(}t^{2}_{r+2}\big{)}=\dfrac{3}{n_{{\mu},r}-2}. Similarly, according to (42) and (44), the state of agent al,l=2,,r+1a_{l},~l=2,\dots,r+1, at instant tr+22t^{2}_{r+2} is xl(tr+22)=53nμ,r2x_{l}\big{(}t^{2}_{r+2}\big{)}=5-\dfrac{3}{n_{{\mu},r}-2}. Recall from (46) that xj(tr+22)=5,j=r+2,,nμ,rx_{j}\big{(}t^{2}_{r+2}\big{)}=5,~\forall j=r+2,\dots,n_{{\mu},r}. Thus, for l=2,,r+1l=2,\dots,r+1, we get

zl(tr+22)=k=1,klnμ,r(xl(tr+22)xk(tr+22))=53(nμ,rr+1)nμ,r2z_{l}\big{(}t^{2}_{r+2}\big{)}=\!\!\!\sum_{k=1,k\neq l}^{n_{{\mu},r}}\!\!\!\Big{(}x_{l}\big{(}t^{2}_{r+2}\big{)}-x_{k}\big{(}t^{2}_{r+2}\big{)}\Big{)}\!=5-\dfrac{3(n_{{\mu},r}-r+1)}{n_{{\mu},r}-2} (49)

Then, it follows from definitions of nr2n_{r_{2}} and nμ,rn_{{\mu},r} in (38) that zl(tr+22)[2,3][α,α]z_{l}\big{(}t^{2}_{r+2}\big{)}\in\![2,3]\!\subset\![-\alpha,\alpha]. ​After that, Lemmas 6 and 7 give

|zl(t)|α,ttr+22,l=2,,r+1|z_{l}(t)|\leq\alpha,~~\forall t\geq t^{2}_{r+2},~~\forall l=2,\dots,r+1 (50)

Recall from (47) that tr+22<tl2,l=2,,r+1t^{2}_{r+2}<t^{2}_{l},\forall l=2,\dots,r+1. Hence,

|zl(tl2)|α,l=2,,r+1|z_{l}\big{(}t^{2}_{l}\big{)}|\leq\alpha,~~~\forall l=2,\dots,r+1 (51)

Then, (7) and (41) together lead to

tl3=(tl2+αβnl)=7nμ,r1,l=2,,r+1t^{3}_{l}=\bigg{(}t^{2}_{l}+\dfrac{\alpha}{\beta n_{l}}\bigg{)}=\dfrac{7}{n_{{\mu},r}-1},~\forall l=2,\dots,r+1 (52)

Next, consider any agent aj,j=r+2,,nμ,ra_{j},~j=r+2,\dots,n_{{\mu},r}. Recall that zj(0)=0z_{j}(0)=0. Then, it follows from Lemma 18.22 that

|zj(t)|α,t0,j=r+2,,nμ,r|z_{j}(t)|\leq\alpha,~~\forall t\geq 0,~~\forall j=r+2,\dots,n_{{\mu},r} (53)

Subsequently, (7) and (44) together lead to

tj3=tj2+αβnj=6nμ,r2,j=r+2,,nμ,rt^{3}_{j}=t^{2}_{j}+\dfrac{\alpha}{\beta n_{j}}=\dfrac{6}{n_{{\mu},r}-2},~~~~\forall j=r+2,\dots,n_{{\mu},r} (54)

Recall from (38) that nμ,r>8n_{{\mu},r}>8. Then, it is easy to show that 4nμ,r1<6nμ,r2<7nμ,r1\dfrac{4}{n_{{\mu},r}-1}<\dfrac{6}{n_{{\mu},r}-2}<\dfrac{7}{n_{{\mu},r}-1}. Then, (41), (52) and (54) imply that

tr+12<tr+23==tnμ,r3<t23=t33==tr+13t^{2}_{r+1}<t^{3}_{r+2}=\dots=t^{3}_{n_{{\mu},r}}<t^{3}_{2}=t^{3}_{3}=\dots=t^{3}_{r+1} (55)

Let t12t^{2}_{1} be as given in (39). It follows from the definitions of nr1n_{r_{1}} and nμ,rn_{{\mu},r} in (38) that 7nμ,r1<5r+32r\dfrac{7}{n_{{\mu},r}-1}<\dfrac{5r+3}{2r}. Then, (39) and (52) together lead to

t23=t33==tr+13<t12t^{3}_{2}=t^{3}_{3}=\dots=t^{3}_{r+1}<t^{2}_{1} (56)

44) Computation of xj(tj3)x_{j}\big{(}t^{3}_{j}\big{)} for jr+2j\geq r+2

Consider any agent aj,j=r+2,,nμ,ra_{j},~j=r+2,\dots,n_{{\mu},r}. Recall that |zj(t)|α,t0|z_{j}(t)|\leq\alpha,\forall t\geq 0 and nj=nμ,r2n_{j}=n_{{\mu},r}-2. Then, Lemma 18.11 leads to

xj(tj3)=kSjxk(tj2)nμ,r2=k=2,kjnμ,rxk(tj2)nμ,r2x_{j}\big{(}t^{3}_{j}\big{)}=\dfrac{\sum_{k\in S_{j}}~x_{k}\big{(}t^{2}_{j}\big{)}}{n_{{\mu},r}-2}=\dfrac{\sum_{k=2,k\neq j}^{n_{{\mu},r}}~x_{k}\big{(}t^{2}_{j}\big{)}}{n_{{\mu},r}-2}

It follows from (43), (46) and (47) that (5μ~)<xk(tj2)5,k=2,,nμ,r(5-\widetilde{\mu})<x_{k}\big{(}t^{2}_{j}\big{)}\leq 5,~\forall k=2,\dots,n_{{\mu},r}. Hence, the average of xk(tj2)x_{k}\big{(}t^{2}_{j}\big{)}’s, i.e., xj(tj3)x_{j}\big{(}t^{3}_{j}\big{)}, satisfies (5μ~)<xj(tj3)5(5-\widetilde{\mu})<x_{j}\big{(}t^{3}_{j}\big{)}\leq 5. Recall that x˙j(t)=uj(t)\dot{x}_{j}(t)=u_{j}(t) is constant in the interval [tj2,tj3][t^{2}_{j},t^{3}_{j}]. Recall also that μ~=μ2\widetilde{\mu}=\dfrac{\mu}{2}. Hence,

5μ<5μ~<xj(t)5,t[tj2,tj3],j=r+2,,nμ,r5-\mu<5-\widetilde{\mu}<x_{j}(t)\leq 5,~\forall t\in\big{[}t^{2}_{j},t^{3}_{j}\big{]},~\forall j=r+2,\dots,n_{{\mu},r} (57)

55) Computation of xl(tl3)x_{l}\big{(}t^{3}_{l}\big{)} for l=2,,r+1l=2,\dots,r+1

Consider any agent al,l=2,,r+1a_{l},~l=2,\dots,r+1. Recall from (51) that |zl(tl2)|α|z_{l}\big{(}t^{2}_{l}\big{)}|\leq\alpha. Recall also that nl=nμ,r1n_{l}=n_{{\mu},r}-1. Then, Lemma 18.11 leads to,

xl(tl3)=x1(tl2)nμ,r1+k=2,klnμ,rxk(tl2)nμ,r1x_{l}\big{(}t^{3}_{l}\big{)}=\dfrac{x_{1}\big{(}t^{2}_{l}\big{)}}{n_{{\mu},r}-1}+\dfrac{\sum_{k=2,k\neq l}^{n_{{\mu},r}}x_{k}\big{(}t^{2}_{l}\big{)}}{n_{{\mu},r}-1} (58)

Recall from (43) that (5μ~)<xk(tl2)<5,k=2,,r+1(5-\widetilde{\mu})<x_{k}\big{(}t^{2}_{l}\big{)}<5,~\forall k=2,\dots,r+1. Further, it follows from (47), (55) and (57) that (5μ~)<xk(tl2)5,k=r+2,,nμ,r(5-\widetilde{\mu})<x_{k}\big{(}t^{2}_{l}\big{)}\leq 5,~\forall k=r+2,\dots,n_{{\mu},r}. Then, the definitions of nμ2n_{{\mu}_{2}}, nμ,rn_{{\mu},r} and μ~\widetilde{\mu} lead to (5μ)<xl(tl3)5(5-\mu)<x_{l}\big{(}t^{3}_{l}\big{)}\leq 5 for l=2,,r+1l=2,\dots,r+1. Recall that x˙l(t)=ul(t)\dot{x}_{l}(t)=u_{l}(t) is constant in the interval [tl2,tl3][t^{2}_{l},t^{3}_{l}\big{]}. Hence,

5μ<xl(t)5,t[tl2,tl3],l=2,,r+15-\mu<x_{l}(t)\leq 5,~~\forall t\in\big{[}t^{2}_{l},t^{3}_{l}\big{]},~~\forall l=2,\dots,r+1 (59)

Now, by repeating the arguments which lead to (57) and (59), we can show that

5μ<xi(t)5,ttr+22=3nμr2,i=2,,nμ,r5-\mu<x_{i}(t)\leq 5,~~\forall t\geq t^{2}_{r+2}=\dfrac{3}{n_{{\mu}_{r}}-2},~~\forall i=2,\dots,n_{{\mu},r}

It follows from the definitions of nμ1n_{{\mu}_{1}} and nμ,rn_{{\mu},r} in (38) that μtr+22\mu\geq t^{2}_{r+2}. Hence, 5μ<xi(t)5,tμ5-\mu<x_{i}(t)\leq 5,~\forall t\geq\mu for i=2,,nμ,ri=2,\dots,n_{{\mu},r}. This proves the claim (37) and completes the proof. ∎

D. Proof of Theorem 10

Let ϵ>0\epsilon>0 be the given real number. Define μ^:=ϵ2\widehat{\mu}:=\dfrac{\epsilon}{2} and r^:=max{ceil(2αϵ),ceil(α5ϵ)}\widehat{r}:=\max\bigg{\{}ceil\bigg{(}\dfrac{2\alpha}{\epsilon}\bigg{)},ceil\bigg{(}\dfrac{\alpha}{5-\epsilon}\bigg{)}\bigg{\}}. Let nμ^,r^n_{\widehat{\mu},\widehat{r}} be as defined in (38) corresponding to μ^\widehat{\mu} and r^\widehat{r}. Recall from (50) that |zl(t)|α,ttr+22,l=2,,r+1|z_{l}(t)|\leq\alpha,~\forall t\geq t^{2}_{r+2},~\forall l=2,\dots,r+1, where the expression of tr+22t^{2}_{r+2} is given in (44). Then, it follows from (44), the definitions of nμ1n_{{\mu}_{1}} and nμ,rn_{{\mu},r} in (38) and the definition of nμ^,r^n_{\widehat{\mu},\widehat{r}} that μ^tr+22\widehat{\mu}\geq t^{2}_{r+2}. Thus,

|zi(t)|α,tμ^,i=2,,r+1|z_{i}(t)|\leq\alpha,~~\forall t\geq\widehat{\mu},~~\forall i=2,\dots,r+1 (60)

Further, recall from (53) that

|zi(t)|α,t0,i=r+2,,nμ^,r^|z_{i}(t)|\leq\alpha,~~\forall t\geq 0,~~\forall i=r+2,\dots,n_{\widehat{\mu},\widehat{r}} (61)

Then, (60) and (61) imply that the α\alpha-consensus time is equal to the maximum of μ^{\widehat{\mu}} and the time required to steer z1z_{1} inside [α,α][-\alpha,\alpha].

Let t^1\widehat{t}_{1} be the first time instant at which |z1(t^1)|=α|z_{1}\big{(}\hskip 0.85355pt\widehat{t}_{1}\big{)}|=\alpha. Recall that the initial conditions of the agents in Example 9 are x1(0)=0x_{1}(0)=0 and xi(0)=5,i=2,,nμ^,r^x_{i}(0)=5,~\forall i=2,\dots,n_{\widehat{\mu},\widehat{r}}. Thus, z1(0)=5r<α=3z_{1}(0)=-5r<-\alpha=-3. Then, clearly, z1(t^1)=αz_{1}\big{(}\hskip 0.85355pt\widehat{t}_{1}\big{)}=-\alpha. We know from Lemma 21 that xi(t)(5u^),tμ^,i=2,,nμ^,r^x_{i}(t)\geq(5-\widehat{u}),~\forall t\geq\widehat{\mu},~\forall i=2,\dots,n_{\widehat{\mu},\widehat{r}}. Then, z1(t^1)=αz_{1}\big{(}\hskip 0.85355pt\widehat{t}_{1}\big{)}=-\alpha implies that x1(t^1)(5μ^αr^)x_{1}\big{(}\hskip 0.85355pt\widehat{t}_{1}\big{)}\geq\bigg{(}5-\widehat{\mu}-\dfrac{\alpha}{\widehat{r}}\bigg{)}. Further, the control rule (10) and z1(t^1)=αz_{1}\big{(}\hskip 0.85355pt\widehat{t}_{1}\big{)}=-\alpha together imply that x˙1(t)=β=1,t[0,t^1)\dot{x}_{1}(t)=\beta=1,~\forall t\in\big{[}0,\widehat{t}_{1}\big{)}. Recall that x1(0)=0x_{1}(0)=0. Hence, t^1(5μ^αr^)\widehat{t}_{1}\geq\bigg{(}5-\widehat{\mu}-\dfrac{\alpha}{\widehat{r}}\bigg{)}. Then, it follows from the definitions of μ^\widehat{\mu} and r^\widehat{r} that t^1(5ϵ)\widehat{t}_{1}\geq(5-\epsilon) and t^1μ^\widehat{t}_{1}\geq\widehat{\mu}. Subsequently, (60) and (61) imply that the α\alpha-consensus time T(α,X(0))T\big{(}\alpha,X(0)\big{)} satisfies T(α,X(0))(5ϵ)=(2T(X(0))ϵ)T\big{(}\alpha,X(0)\big{)}\geq(5-\epsilon)=\Big{(}2T^{\ast}\big{(}X(0)\big{)}-\epsilon\Big{)}. This proves the first inequality in (12).

The second inequality in (12) follows from Theorem 8. This completes the proof.

E. Proof of Lemma 11

It is given that |zi(tik)|α|z_{i}(t^{k}_{i})|\leq\alpha. Then, it follows from Lemma 18.22 that control rule (7)-(8) satisfies constraint (16). Let tik+1t^{k+1}_{i} and uiu^{\ast}_{i} be as defined in (7) and (8), respectively. Then, the inter-update duration (tik+1tik)\big{(}t^{k+1}_{i}-t^{k}_{i}\big{)} is αβni\dfrac{\alpha}{\beta n_{i}}, for all k1k\geq 1. Next, we show that αβni\dfrac{\alpha}{\beta n_{i}} is the max-min value of optimization (14)-(16). This will prove the optimality of control rule (7)-(8) and complete the proof.

For the sake of contradiction, assume that the max-min value of optimization (14)-(16) is greater than αβni\dfrac{\alpha}{\beta n_{i}}. Let that value be Ti:=(tik+1tik)>αβniT_{i}:=\big{(}t^{k+1}_{i}-t^{k}_{i}\big{)}>\dfrac{\alpha}{\beta n_{i}}. Let u~i𝒰\widetilde{u}_{i}\in\mathcal{U} be the corresponding optimal control. Then, it follows from the evolution of ziz_{i} given in (13) that

zi(tik+1)=zi(tik)+nitiktik+Tiu~i(t)𝑑tjSi(tiktik+Tiuj(t)𝑑t)z_{i}\big{(}t^{k+1}_{i}\big{)}=z_{i}(t^{k}_{i})+n_{i}\int_{t^{k}_{i}}^{t^{k}_{i}+T_{i}}\!\!\!\!\!\!\widetilde{u}_{i}(t)~dt-\sum_{j\in S_{i}}\Bigg{(}\int_{t^{k}_{i}}^{t^{k}_{i}+T_{i}}\!\!\!\!\!\!\!u_{j}(t)~dt\Bigg{)} (62)

Define z~i(tik+1):=zi(tik)+nitiktik+Tiu~i(t)𝑑t\widetilde{z}_{i}\big{(}t^{k+1}_{i}\big{)}:=z_{i}\big{(}t^{k}_{i}\big{)}+n_{i}\int_{t^{k}_{i}}^{t^{k}_{i}+T_{i}}\widetilde{u}_{i}(t)~dt. Then, depending on the value of z~i(tik+1)\widetilde{z}_{i}\big{(}t^{k+1}_{i}\big{)}, there are following three cases.
Case 1 : z~i(tik+1)|>α\widetilde{z}_{i}\big{(}t^{k+1}_{i}\big{)}\big{|}>\alpha
Take uj(t)=0,t[tik,tik+Ti],jSiu_{j}(t)=0,~\forall t\in[t^{k}_{i},t^{k}_{i}+T_{i}],~\forall j\in S_{i}. Then, it follows from (62) that |zi(tik+1)|>α\big{|}z_{i}\big{(}t^{k+1}_{i}\big{)}\big{|}>\alpha.
Case 2 : 0z~i(tik+1)α0\leq\widetilde{z}_{i}\big{(}t^{k+1}_{i}\big{)}\leq\alpha
Take uj(t)=β,t[tik,tik+Ti],jSiu_{j}(t)=-\beta,~\forall t\in[t^{k}_{i},t^{k}_{i}+T_{i}],~\forall j\in S_{i}. Then, it follows from (62) that

zi(tik+1)=z~i(tik+1)+jSiβ(tik+Titik)=z~i(tik+1)+βniTiz_{i}\big{(}t^{k+1}_{i}\big{)}=\widetilde{z}_{i}\big{(}t^{k+1}_{i}\big{)}+\sum_{j\in S_{i}}\beta\big{(}t^{k}_{i}+T_{i}-t^{k}_{i}\big{)}=\widetilde{z}_{i}\big{(}t^{k+1}_{i}\big{)}+\beta n_{i}T_{i}

Recall that Ti>αβniT_{i}>\dfrac{\alpha}{\beta n_{i}}, which is equivalent to βniTi>α\beta n_{i}T_{i}>\alpha. Then, the fact z~i(tik+1)0\widetilde{z}_{i}\big{(}t^{k+1}_{i}\big{)}\geq 0 implies that |zi(tik+1)|>α\big{|}z_{i}\big{(}t^{k+1}_{i}\big{)}\big{|}>\alpha.
Case 3 : αz~i(tik+1)<0-\alpha\leq\widetilde{z}_{i}\big{(}t^{k+1}_{i}\big{)}<0
Take uj(t)=β,t[tik,tik+1],jSiu_{j}(t)=\beta,~\forall t\in[t^{k}_{i},t^{k+1}_{i}],~\forall j\in S_{i}. Then, by following the arguments in Case 22, we get |zi(tik+1)|>α\big{|}z_{i}\big{(}t^{k+1}_{i}\big{)}\big{|}>\alpha.

In each of the above cases, for control u~i\widetilde{u}_{i}, there exist uj𝒰,jSiu_{j}\in\mathcal{U},\forall j\in S_{i}, such that |zi(tik+1)|>α\big{|}z_{i}\big{(}t^{k+1}_{i}\big{)}\big{|}>\alpha. This violates constraint (16) and in result, contradicts the assumption that u~i\widetilde{u}_{i} is a solution of optimization (14)-(16). This contradiction proves our claim that αβni\dfrac{\alpha}{\beta n_{i}} is the max-min value of optimization (14)-(16) and completes the proof.

F. Proof of Lemma 12

We first show that control rule (9)-(10) satisfies constraint (19). Let tik+1t^{k+1}_{i} and uiu^{\ast}_{i} be as defined in (9) and (10), respectively. Define Ti:=tik+1tik=|zi(tik)|+α2βniT_{i}:=t^{k+1}_{i}-t^{k}_{i}=\dfrac{|z_{i}\big{(}t^{k}_{i}\big{)}|+\alpha}{2\beta n_{i}}. It is given that zi(tik)<αz_{i}\big{(}t^{k}_{i}\big{)}<-\alpha. Then, (10) leads to ui(t)=β,t[tik,tik+Ti]u^{\ast}_{i}(t)=\beta,~\forall t\in[t^{k}_{i},t^{k}_{i}+T_{i}]. Recall the evolution of ziz_{i} given in (13). By putting the expressions of tik+1t^{k+1}_{i} and uiu^{\ast}_{i} in (13), we get

zi(tik+1)=zi(tik)+|zi(tik)|+α2jSi(tiktik+Tiuj(t)𝑑t)z_{i}\big{(}t^{k+1}_{i}\big{)}=z_{i}\big{(}t^{k}_{i}\big{)}+\dfrac{\big{|}z_{i}\big{(}t^{k}_{i}\big{)}\big{|}+\alpha}{2}-\sum_{j\in S_{i}}\Bigg{(}\int_{t^{k}_{i}}^{t^{k}_{i}+T_{i}}\!\!u_{j}(t)~dt\Bigg{)} (63)

Recall that the cardinality of SiS_{i} is nin_{i} and |uj(t)|β,t0,jP|u_{j}(t)|\leq\beta,~\forall t\geq 0,~\forall j\in P. Then, it follows from the definition of TiT_{i} that

|jSi(tiktik+Tiuj(t)𝑑t)||zi(tik)|+α2\Bigg{|}\sum_{j\in S_{i}}\Bigg{(}\int_{t^{k}_{i}}^{t^{k}_{i}+T_{i}}\!\!u_{j}(t)~dt\Bigg{)}\Bigg{|}\leq\dfrac{\big{|}z_{i}\big{(}t^{k}_{i}\big{)}\big{|}+\alpha}{2} (64)

Now, (63) and (64) together imply that zi(tik+1)(zi(tik)+|zi(tik)|+α)z_{i}\big{(}t^{k+1}_{i}\big{)}\leq\Big{(}z_{i}\big{(}t^{k}_{i}\big{)}+|z_{i}\big{(}t^{k}_{i}\big{)}|+\alpha\Big{)}. As zi(tik)<0z_{i}\big{(}t^{k}_{i}\big{)}<0, i.e., zi(tik)+|zi(tik)|=0z_{i}\big{(}t^{k}_{i}\big{)}+|z_{i}\big{(}t^{k}_{i}\big{)}|=0, we have zi(tik+1)αz_{i}\big{(}t^{k+1}_{i}\big{)}\leq\alpha. Further, it follows from the definition of ziz_{i} in (2) that z˙i(t)=(niui(t)jSiuj(t))=(niβjSiuj(t)),t[tik,tik+1]\dot{z}_{i}(t)=\Big{(}n_{i}u^{\ast}_{i}(t)-\sum_{j\in S_{i}}u_{j}(t)\Big{)}=\Big{(}n_{i}\beta-\sum_{j\in S_{i}}u_{j}(t)\Big{)},~\forall t\in[t^{k}_{i},t^{k+1}_{i}]. Recall again that the cardinality of SiS_{i} is nin_{i} and |uj(t)|β,t0,jP|u_{j}(t)|\leq\beta,~\forall t\geq 0,~\forall j\in P. Hence, z˙i(t)0,t[tik,tik+1]\dot{z}_{i}(t)\geq 0,~\forall t\in[t^{k}_{i},t^{k+1}_{i}] which implies that zi(t)zi(tik+1),t[tik,tik+1]z_{i}(t)\leq z_{i}\big{(}t^{k+1}_{i}\big{)},~\forall t\in[t^{k}_{i},t^{k+1}_{i}]. Then, the above-proved fact zi(tik+1)αz_{i}\big{(}t^{k+1}_{i}\big{)}\leq\alpha implies that control law (9)-(10) satisfies constraint (19).

Next, we show that control rule (9)-(10) achieves the max-min value of optimization (18)-(19). Recall that under control rule (9)-(10), the inter-update duration is Ti=|zi(tik)|+α2βniT_{i}=\dfrac{|z_{i}\big{(}t^{k}_{i}\big{)}|+\alpha}{2\beta n_{i}}. We claim that TiT_{i} is the max-min value of optimization (18)-(19). For the sake of contradiction, assume that the max-min value of optimization (18)-(19) is (Ti+ϵ)(T_{i}+\epsilon) for some ϵ>0\epsilon>0. Then, (tik+1tik)=Ti+ϵ\big{(}t^{k+1}_{i}-t^{k}_{i}\big{)}=T_{i}+\epsilon. Let the inputs of the neighbouring agents be uj(t)=β,t[tik,tik+1),jSiu_{j}(t)=-\beta,~\forall t\in[t^{k}_{i},t^{k+1}_{i}),~\forall j\in S_{i}. Then, it follows from the evolution of ziz_{i} given in (13) that

zi(tik+1)=zi(tik)+nitiktik+Ti+ϵβ𝑑t+jSi(tiktik+Ti+ϵβ𝑑t)z_{i}\big{(}t^{k+1}_{i}\big{)}=z_{i}\big{(}t^{k}_{i}\big{)}+n_{i}\int_{t^{k}_{i}}^{t^{k}_{i}+T_{i}+\epsilon}\!\!\!\!\!\!\beta~dt+\sum_{j\in S_{i}}\Bigg{(}\int_{t^{k}_{i}}^{t^{k}_{i}+T_{i}+\epsilon}\beta~dt\Bigg{)}

Now, by putting the expression of TiT_{i} and rearranging the terms, we get zi(tik+1)=(zi(tik)+|zi(tik)|+α+2βniϵ)z_{i}\big{(}t^{k+1}_{i}\big{)}=\Big{(}z_{i}\big{(}t^{k}_{i}\big{)}+|z_{i}\big{(}t^{k}_{i}\big{)}|+\alpha+2\beta n_{i}\epsilon\Big{)}. As zi(tik)<0z_{i}\big{(}t^{k}_{i}\big{)}<0, we have zi(tik)+|zi(tik)|=0z_{i}\big{(}t^{k}_{i}\big{)}+|z_{i}\big{(}t^{k}_{i}\big{)}|=0. Thus, zi(tik+1)=(α+2βniϵ)>αz_{i}\big{(}t^{k+1}_{i}\big{)}=\big{(}\alpha+2\beta n_{i}\epsilon\big{)}>\alpha. This violates constraint (19) and contradicts the fact that (Ti+ϵ)(T_{i}+\epsilon) is the max-min value of optimization (18)-(19). Hence, the inter-update duration TiT_{i} under control rule (9)-(10) is the max-min value of optimization (18)-(19). This completes the proof.

References

  • [1] J. Almeida, C. Silvestre and A. Pascoal, Continuous-time consensus with discrete-time communications, Systems and Control Letters, vol. 61, no. 7, pp. 788-796, 2012.
  • [2] J. Baillieul and P. Antsaklis, Control and communication challenges in networked real-time systems, Proceedings of the IEEE, vol. 95, no. 1, pp. 9-28, 2007.
  • [3] E. Barmpounakis, E. Vlahogianni and J. Golias, Unmanned aerial aircraft systems for transportation engineering: current practice and future challenges, International Journal of Transportation Science and Technology, vol. 5, no. 3, pp. 111-122, 2016.
  • [4] E. Bertran and A. Cerda, On the tradeoff between electrical power consumption and flight performance in fixed-wing UAV autopilots, IEEE Transactions on Vehicular Technology, vol. 56, no. 11, pp. 8832-8840, 2016.
  • [5] M. Cao, A. Morse and B. Anderson, Coordination of an asynchronous multi-agent system via averaging, International Federation of Automatic Control (IFAC) World Congress, Prague, Czech Republic, 2005.
  • [6] M. Cao, A. Morse and B. Anderson, Agreeing asynchronously, IEEE Transactions on Automatic Control, vol. 53, no. 8, pp. 1826-1838, 2008.
  • [7] L. Corral, I. Fronza, N. Ioini and A. Ibershimi, A measurement tool to track drones battery consumption during flights, Proceedings of International Conference on Mobile Web and Information Systems, Vienna, Austria, 2016.
  • [8] D. Dimarogonas, E. Frazzoli, and K. Johansson, Distributed self-triggered control for multi-agent systems, Proceedings of the IEEE Conference on Decision and Control, Atlanta, USA, 2010.
  • [9] D. Dimarogonas, E. Frazzoli and K. Johansson, Distributed event-triggered control for multi-agent systems, IEEE Transactions on Automatic Control, vol. 57, no. 5, pp. 1291-1297, 2012.
  • [10] L. Fang and P. Antsaklis, Information consensus of asynchronous discrete-time multi-agent systems, Proceedings of the American Control Conference, Portland, USA, 2005.
  • [11] J. Fax and R. Murray, Information flow and cooperative control of vehicle formations, IEEE Transactions on Automatic Control, vol. 49, no. 9, pp. 1465-1476, 2004.
  • [12] S. Giannini, D. Paola, A. Petitti and A. Rizzo, On the convergence of the max-consensus protocol with asynchronous updates, Proceedings of the IEEE Conference on Decision and Control, Florence, Italy, 2013.
  • [13] W. Heemels, K. Johansson and P. Tabuada, An introduction to event-triggered and self-triggered control, Proceedings of the IEEE Conference on Decision and Control, Hawaii, USA, 2012.
  • [14] K. Hla, Y. Choi and J. Park, The multi-agent system solutions for wireless sensor network applications, Agent and Multi-agent Systems: Technologies and Applications, pp. 454-463, 2008.
  • [15] A. Jadbabaie and A. Morse, Coordination of groups of mobile autonomous agents using nearest neighbor rules, IEEE Transactions on Automatic Control, vol. 48, no. 6, pp. 988-1001, 2003.
  • [16] A. Joshi, N. Limbu, I. Ahuja, A. Mulla, H. Chung, and D. Chakraborty, Implementation of distributed consensus on an outdoor testbed, Proceedings of the European Control Conference, Aalborg, Denmark, 2016.
  • [17] H. Liu, S. Chen, L. Shen and J. Chen, Tactical trajectory planning for stealth unmanned aerial vehicle to win the radar game, Defence Science Journal, vol. 62, no. 6, pp. 375-381, 2012.
  • [18] G. Merabet, M. Essaaidi, H. Talei, M. Abid, N. Khalil, M. Madkour and D. Benhaddou, Applications of Multi-Agent Systems in Smart Grids: A Survey, International Conference on Multimedia Computing and Systems (ICMCS), Marrakesh, Morocco, 2014.
  • [19] J. Meyer, F. Du-plessis and W. Clarke, Design considerations for long endurance unmanned aerial vehicles, IntechOpen, 2009.
  • [20] K. Movric and F. Lewis, Cooperative optimal control for multi-agent systems on directed graph topologies, IEEE Transactions on Automatic Control, vol. 59, no. 3, pp. 769-774, 2014.
  • [21] A. Mulla, D. Patil and D. Chakraborty, Distributed computation of minimum time consensus for multi-agent systems, Proceedings of the IEEE Conference on Decision and Control, Los Angeles, USA, 2014.
  • [22] G. Nair, F. Fagnani, S. Zampieri and R. Evans, Feedback control under data rate constraints: an overview, Proceedings of the IEEE, vol. 95, no. 1, pp. 108-137, 2007.
  • [23] W. Ren, Consensus tracking under directed interaction topologies: algorithms and experiments, IEEE Transactions on Control Systems Technology, vol. 18, no. 1, pp. 230-237, 2010.
  • [24] W. Ren and R. Beard, Distributed Consensus in Multi-vehicle Cooperative Control, Springer, 2008.
  • [25] V. Sawant, D. Chakraborty and D. Pal, Asynchronous consensus of continuous-time multi-agent systems with minimum communication, Proceedings of the IEEE Conference on Decision and Control, Nice, France, 2019.
  • [26] J. Shamma, Cooperative Control of Distributed Multi-agent Systems, Wiley, 2007.
  • [27] G. Wen, Z. Duan, W. Ren and G. Chen, Distributed consensus of multi-agent systems with general linear node dynamics and intermittent communications, International Journal of Robust and Nonlinear Control, vol. 24, no. 16, pp. 2438-2457, 2013.
  • [28] G. Wen, Z. Duan, W. Yu and G. Chen, Consensus in multi-agent systems with communication constraints, International Journal of Robust and Nonlinear Control, vol. 22, no. 2, pp. 170-182, 2012.
  • [29] G. Wen, Z. Duan, W. Yu and G. Chen, Consensus of second-order multi-agent systems with delayed nonlinear dynamics and intermittent communications, International Journal of Control, vol. 86, no. 2, pp. 322-331, 2012.
  • [30] W. Xia amnd M. Cao, Sarymsakov matrices and asynchronous implementation of distributed coordination algorithms, IEEE Transactions on Automatic Control, vol. 59, no. 8, pp. 2228-2233, 2014.
  • [31] F. Xiao and L. Wang, Asynchronous consensus in continuous-time multi-agent systems with switching topology and time-varying delays, IEEE Transactions on Automatic Control, vol. 53, no. 8, 1804-1816, 2008.
  • [32] X. Yi, J. Wei and K. Johansson, Self-triggered control for multi-agent systems with quantized communication or sensing, Proceedings of the IEEE Conference on Decision and Control, Las Vegas, USA, 2016.
  • [33] J. Zhan and X. Li, Asynchronous consensus of multiple double-integrator agents with arbitrary sampling intervals and communication delays, IEEE Transactions on Circuits and Systems-\Romannum1\Romannum{1}, vol. 62, no. 9, pp. 2301-2311, 2015.
  • [34] Z. Zhao, Y. Niu, Z. Ma and X. Ji, A fast stealth trajectory planning algorithm for stealth UAV to fly in multi-radar network, Proceedings of the IEEE International Conference on Real-time Computing and Robotics, Angkor Wat, Cambodia, 2016.
  • [35] W. Zhu, Z. Jiang and G. Feng, Event-based consensus of multi-agent systems with general linear models, Automatica, vol. 50, no. 2, pp. 552-558, 2014.
Vishal Sawant received the Bachelor of Engineering degree from Mumbai University, India, in 2010 and the M.Tech. degree from the Department of Electrical Engineering, IIT Bombay, India, in 2012. He then worked as a software engineer for two years. In July 2014, he started Ph.D. degree at the Department of Electrical Engineering, IIT Bombay, India. From May 2020 to September 2021, he worked as a postdoctoral researcher at the Cognitive Robotics group, TU Delft, The Netherlands. Currently, he is a postdoctoral researcher at the Automation and Control Section, Aalborg University, Denmark. His research interests include multi-agent systems, cyber-physical security and distributed optimization.
Debraj Chakraborty received the Bachelor’s degree from Jadhavpur University, Kolkata, India, in 2001 and the M. Tech. degree from Indian Institute of Technology Kanpur, Kanpur, India, in 2003, both in electrical engineering, and the Ph. D. degree in electrical and computer engineering from University of Florida, Gainesville, FL, USA, in 2007. He joined the Indian Institute of Technology Bombay, Mumbai, India, in 2007, where he is currently a Professor in the Department of Electrical Engineering. His research interests include optimal control, linear systems, and multiagent systems.
Debasattam Pal received his Bachelor of Engineering (B.E.) degree from the Department of Electrical Engineering of Jadavpur University, Kolkata, in 2005. He received his M.Tech. and Ph.D. degrees from the Department of Electrical Engineering, IIT Bombay, in the years 2007 and 2012, respectively. He then worked as an Assistant Professor in IIT Guwahati from July, 2012 to May, 2014. He joined IIT Bombay in June, 2014, where he is currently an Associate Professor in the EE Department. His main area of research is systems and control theory. More specifically, his areas of interest are: multidimensional systems theory, algebraic analysis of systems, dissipative systems, optimal control and computational algebra.