Consensus-based adaptive guaranteed-performance formation control for second-order multi-agent systems

Abstract

Adaptive guaranteed-performance formation analysis and design problems for second-order multi-agent systems are studied, where the global information is not required, which means the main results of this article are fully distributed. First, an adaptive guaranteed-performance formation control protocol is presented for second-order multi-agent systems, where the control input is constructed using neighboring state errors and adaptively adjustable interaction weights. Then, an adaptive guaranteed-performance formation control is proposed based on Riccati inequalities. Furthermore, the guaranteed-performance cost is determined and the adjusting approach of the formation control gain is presented in terms of linear matrix inequalities. Finally, a numerical simulation is provided to demonstrate the effectiveness of the theoretical results.

Keywords

Guaranteed-performance control multi-agent systems consensus approach formation control adaptive control

Introduction

In the past few years, distributed cooperative control of multi-agent systems has aroused great attention from both scientific and engineering communities.^1–6 Multi-agent cooperative control problems in previous works^7–17 were solved by consensus-based approach, which requires that the specific state variables of all agents reach an agreement. As an important research topic of cooperative control of multi-agent systems, formation control has been widely applied in transportation, surveillance and rescue, and so on. As a point of fact, formation control problems have been extensively studied during the past decades, and some typical formation control approaches have been proposed, such as leader–follower, virtual structure, and behavior-based ones.^18–20 Recently, motivated by the development of consensus theory, the formation control problems were solved by the consensus-based approach as in previous works.^21–25

Multi-agent systems in previous works^26–28 achieved formation without considering the control performance, however, multi-agent systems in practical tasks are usually required to achieve certain performance, such as the shortest time and the least energy consumption. For example, Xu et al.²⁹ optimized the acceleration control parameters using the particle swarm optimization algorithm to reduce the energy consumption in the multiple robotic fishes’ formation process. The formation performance usually can be modeled as certain cost functions. A quadratic performance function is constructed to achieve guaranteed-performance time-varying formation in Wang et al.³⁰ A cost function was constructed by synchronization errors and control input to describe the formation performance of the multi-agent system in Cheng and Ugrinovskii.³¹ As a matter of fact, the global information of the multi-agent system is required in Wang et al.³⁰ and Cheng and Ugrinovskii,³¹ which means their main results are not fully distributed.

A distributed formation control algorithm is proposed for multi-vehicle nonlinear system model with simple kinetics characteristics coupling time-delay and jointly connected topologies in Xue et al.³² and Savkin et al.³³ proposed a distributed motion control algorithm for groups of robots. It should be pointed that the formation control results in Xue et al.³² and Savkin et al.³³ are fully distributed, while the formation performance was not considered. To the best of our knowledge, fully distributed guaranteed-performance formation control problems for second-order multi-agent systems with adaptive control scheme are still open.

Motivated by the facts stated above, this article investigated adaptive guaranteed-performance formation control for second-order multi-agent systems. First, an adaptive guaranteed-performance formation control protocol is proposed. Then, an adaptive guaranteed-performance formation control criterion is presented based on Riccati inequalities, and the guaranteed-performance cost is determined. Moreover, a formation control gain regulation approach is proposed.

Compared with the related works on formation control of second-order multi-agent systems, the new features of this article are threefold. First, a translation-adaptive strategy is proposed to deal with the fully distributed formation control problems for second-order multi-agent systems. However, the global information is required in Wang et al.³⁰ and Cheng and Ugrinovskii.³¹ Second, the multi-agent systems in this article can achieve guaranteed-performance formation control, while the formation performance was not considered in Xue et al.³² and Savkin et al.³³ Third, the formation control gain regulation strategy is proposed based on linear matrix inequality (LMI) techniques, while the formation control gains in previous works^26–28 cannot be adjusted.

The remainder of this article is organized as follows. In section “Problem description,” the problem description and the formation control protocol are proposed. The main results are given under section “Main results.” Section “Numerical simulation” gives numerical simulations. Concluding remarks are presented in the final section.

Notations: let $R^{d}$ stand for the d-dimension real column vector space and $R^{d \times d}$ represent the set of $d \times d$ dimensional real matrices. For simplicity, let 0 denote zero matrices of any size with zero vector and zero number as special cases. $I_{N}$ stands for the identity matrix of dimension $N$ . $1_{N}$ represents a column vector of dimension $N$ , whose entries are equal to 1. $P^{T}$ stands for the transpose of the symmetric matrix $P$ . $Q^{T} = Q > 0$ and $Q^{T} = Q < 0$ mean that $Q$ is positive semidefinite and negative semidefinite, respectively. The symbol ⊗ denotes the Kronecker product. The notation * stands for the symmetric terms of a symmetric matrix.

Problem description

Consider a second-order multi-agent system consisting of $N$ homogeneous agents, whose dynamics is described as follows

{\begin{matrix} {\overset{\cdot}{x}}_{i} (t) = v_{i} (t) \\ {\overset{\cdot}{v}}_{i} (t) = u_{i} (t) \end{matrix}

(1)

where $i \in I_{N}$ with $I_{N} = {1, 2, \dots, N}$ . $x_{i} (t) \in R^{d}$ , $v_{i} (t) \in R^{d}$ and $u_{i} (t) \in R^{d}$ are the position state, the velocity state, and the control input of agent $i$ , respectively. The interaction topology of multi-agent system (1) can be described by a connected undirected graph $G = G (V, E)$ , where $V = {V_{1}, V_{2}, \dots, V_{N}}$ denotes the set of nodes, $E \subseteq {(V_{i}, V_{j}), i \neq j, V_{i}, V_{j} \in V}$ denotes the set of edges. In the graph $G$ , each node represents an agent in the system. The edge between two nodes represents the interaction relationship between agents, and the length of edge represents the strength of the interaction between agents. $N_{i}$ stands for the neighboring set of agent $i$ . The weight matrix is represented by a symmetric adjacency matrix $W (t) = [l_{ij} w_{ij} (t)]_{N \times N}$ , where $l_{ii} = 0$ and $l_{ij} = 0$ if agent $j$ is not a neighbor of agent $i$ and $l_{ij} = 1$ otherwise, $w_{ij} (t)$ with $w_{ij} (t) > 0$ denotes the interaction weight from agent $j$ to agent $i$ . If $l_{ij} = 0$ , then $w_{ij} (t)$ denotes a virtual interaction weight from agent $j$ to agent $i$ . $D (t) = diag {d_{1} (t), d_{2} (t), \dots, d_{N} (t)}$ with $d_{i} (t) = \sum_{j = 1, j \neq i}^{N} l_{ij} w_{ij} (t)$ stands for the in-degree matrix of $G$ . $L_{w (t)} = D (t) - W (t)$ denotes the Laplacian matrix of $G$ . The properties of the Laplacian matrix for an undirected graph are shown in the following lemma.

Lemma 1

Let $G$ denote an undirected graph, then its Laplacian matrix $L \in R^{N \times N}$ has at least nonzero eigenvalue, and the corresponding eigenvector is $1_{N}$ . If $G$ is connected, $0$ is a simple eigenvalue of $L$ and all the other $N - 1$ eigenvalues are positive.³⁴

Definition 1

Multi-agent system (1) is said to achieve consensus if the following conditions are satisfied

\begin{matrix} lim_{t \to \infty} (x_{i} (t) - x_{j} (t)) = 0 \\ lim_{t \to \infty} (v_{i} (t) - v_{j} (t)) = 0 \end{matrix}

(2)

where $i, j \in I_{N}$ .

Definition 2

For a given vector $ι = [ι_{1}^{T}, ι_{2}^{T}, \dots, ι_{N}^{T}]^{T}$ with $ι_{i} \in R^{d}$ , multi-agent system (1) achieves formation $ι$ if the following conditions are satisfied

\begin{matrix} lim_{t \to \infty} (x_{i} (t) - x_{j} (t)) = ι_{i} - ι_{j} \\ lim_{t \to \infty} (v_{i} (t) - v_{j} (t)) = 0 \end{matrix}

(3)

where $i, j \in I_{N}$ , the vector $ι$ is called a formation vector.

Let $δ_{ij} (t) = [δ {ε_{ij}}^{T} (t), δ {v_{ij}}^{T} (t)]^{T}$ , where $δ ε_{ij} (t) = ε_{j} (t) - ε_{i} (t)$ with $ε_{i} (t) = x_{i} (t) - ι_{i}$ , and $δ v_{ij} (t) = δ v_{j} (t) - δ v_{i} (t)$ . Since ${\overset{\cdot}{ε}}_{i} (t) = {\overset{\cdot}{x}}_{i} (t) = v_{i} (t)$ , it can be obtained that

{\begin{matrix} {\overset{\cdot}{ε}}_{i} (t) = v_{i} (t) \\ {\overset{\cdot}{v}}_{i} (t) = u_{i} (t) \end{matrix}

(4)

and formation control criterion (3) can be rewritten as

\begin{matrix} lim_{t \to \infty} (ε_{i} (t) - ε_{j} (t)) = 0 \\ lim_{t \to \infty} (v_{i} (t) - v_{j} (t)) = 0 \end{matrix}

(5)

Remark 1

When the formation is achieved, it can be obtained by Definition 2 that $lim_{t \to \infty} (x_{i} (t) - x_{j} (t)) = ι_{i} - ι_{j}$ , which does not mean that the position states of the agents accurately converge to the corresponding formation vector under the formation control protocol. Moreover, if $ι = 0_{d \times N}$ , the formation control criterion (3) can be rewritten as the consensus criterion (2). Therefore, the consensus problem is a special case of formation problem. More generally, by equations (4) and (5), the formation control criterion (3) is written in a form similar to the consensus criterion (2); that is, a formation control problem can be described as a consensus problem.

Definition 3

Multi-agent system (1) is said to achieve adaptive guaranteed-performance formation by protocol (10) if there exist $K_{u}$ and $K_{w}$ such that

lim_{t \to \infty} δ_{ij} (t) = 0 (i, j = 1, 2, \dots, N)

(6)

and $J_{r} < J_{r}^{*}$ for any bounded initial states $x_{i} (0) (i = 1, 2, \dots, N)$ , where $J_{r}^{*}$ is said to be the guaranteed-performance cost.

The states of second-order multi-agent system (1) include position state and velocity state, so a control input $u_{i} (t)$ is proposed as follows

u_{i} (t) = u_{i ε} (t) + u_{iv} (t)

(7)

with

\begin{matrix} u_{i ε} (t) = K_{u 1} \sum_{j \in N_{i}} w_{ij} (ε_{j} (t) - ε_{i} (t)) \\ u_{iv} (t) = K_{u 2} \sum_{j \in N_{i}} w_{ij} (v_{j} (t) - v_{i} (t)) \end{matrix}

(8)

where $K_{u 1} \in R^{d \times d}$ and $K_{u 2} \in R^{d \times d}$ are the control gain matrices. Let $K_{u} = [K_{u 1}, K_{u 2}]$ , then the control input $u_{i} (t)$ can be rewritten as

u_{i} (t) = K_{u} \sum_{j \in N_{i}} w_{ij} (t) δ_{ij} (t)

(9)

Based on equation (9), a guaranteed-performance formation control protocol of second-order multi-agent system (1) is proposed as follows

{\begin{matrix} u_{i} (t) = K_{u} \sum_{j \in N_{i}} w_{ij} δ_{ij} (t) \\ {\overset{\cdot}{w}}_{ij} (t) = {δ_{ij}}^{T} (t) K_{w} δ_{ij} (t) \\ J_{r} = \frac{1}{N} \sum_{i = 1}^{N} \sum_{j = 1}^{N} \int_{0}^{+ \infty} {δ_{ij}}^{T} (t) Q δ_{ij} (t) dt \end{matrix}

(10)

where $i \in I_{N}$ , $K_{w}$ with $K_{w} = K_{w}^{T} \geq 0$ is a gain matrix and $Q = Q^{T} > 0$ . It should be pointed out that $w_{ij} (t)$ is constantly positive and adaptively adjusted, and it is consistently increasing unless the states of agent $i$ and agent $j$ achieve consensus. Let $γ_{ij}$ be the upper bound of $w_{ij} (t)$ .

Main results

Let $χ (t) = [χ_{1}^{T} (t), χ_{2}^{T} (t), \dots, χ_{N}^{T} (t)]^{T}$ with ${χ_{i}}^{T} (t) = [x_{i}^{T} (t), v_{i}^{T} (t)]^{T}$ , then multi-agent system (1) can be described in global form as

\overset{\cdot}{χ} (t) = (I_{N} \otimes A) χ (t) - (I_{N} \otimes B) u (t)

(11)

where

A = [\begin{matrix} 0 & I_{d} \\ 0 & 0 \end{matrix}], B = [\begin{matrix} 0 \\ I_{d} \end{matrix}]

(12)

Let $l = [l_{1}^{T}, l_{2}^{T}, \dots l_{N}^{T}]$ with $l_{i} = [ι_{i}^{T}, 0]^{T}$ represent the formation vector of system (11) and $ω (t) = [{ω_{1}}^{T} (t), {ω_{2}}^{T} (t), \dots, {ω_{N}}^{T} (t)]^{T}$ with $ω_{i}^{T} (t) = [χ_{i}^{T} (t) - ι_{i}^{T}]$ , then $ω (t) = χ (t) - l$ .

Due to formation vector $l$ being a constant vector, $\overset{\cdot}{l} = 0$ . Since $ω (t) = χ (t) - l$ , it can be obtained by equations (9) and (11) that

\overset{\cdot}{ω} (t) = (I_{N} \otimes A - L_{w (t)} \otimes B K_{u}) ω (t) + (I_{N} \otimes A) l

(13)

By the special structure of $I_{N}$ , $A$ , and $l$ , one has

(I_{N} \otimes A) l = 0

(14)

So, equation (13) can be converted into

\overset{\cdot}{ω} (t) = (I_{N} \otimes A - L_{w (t)} \otimes B K_{u}) ω (t)

(15)

where $L_{w (t)}$ denotes the Laplacian matrix of $G$ .

From Lemma 1, one has $L_{w (t)} 1_{N} = 0$ . Because $G$ is connected, there exists an orthonormal matrix $U = [1_{N} / \sqrt{N}, \bar{U}]$ satisfies that

U^{T} L_{w (0)} U = [\begin{matrix} 0 & 0^{T} \\ 0 & Λ_{λ} \end{matrix}]

(16)

where $Λ_{λ} = diag {λ_{2}, λ_{3}, \dots, λ_{N}}$ with $0 < λ_{2} \leq \dots \leq λ_{N}$ . Let $κ (t) = (U^{T} \otimes I_{2 d}) ω (t) = [κ_{1}^{T} (t), κ_{e}^{T} (t)]^{T}$ , where $κ_{1} (t) \in R^{2 d}$ , $κ_{e} (t) = [κ_{2}^{T} (t), κ_{3}^{T} (t), \dots, κ_{N}^{T} (t)]^{T} \in R^{2 d (N - 1)}$ , then equation (15) can be rewritten as

{\overset{\cdot}{κ}}_{1} (t) = A κ_{1} (t)

(17)

\begin{matrix} {\overset{\cdot}{κ}}_{e} (t) = (I_{N - 1} \otimes A - Λ_{λ} \otimes B K_{u}) κ_{e} (t) \\ = (I_{N - 1} \otimes A - {\bar{U}}^{T} L_{w (t)} \bar{U} \otimes B K_{u}) κ_{e} (t) \end{matrix}

(18)

Define

{\tilde{κ}}_{1} (t) = U e_{1} \otimes κ_{1} (t) = \frac{1}{\sqrt{N}} 1_{N} \otimes κ_{1} (t)

(19)

{\tilde{κ}}_{e} (t) = \sum_{i = 2}^{N} U e_{i} \otimes κ_{i} (t)

(20)

where $e_{i} (i \in I_{N})$ denotes N-dimensional column vectors with the i-th element 1 and 0 elsewhere.

It can be derived from equations (19) and (20) that

{\tilde{κ}}_{1} (t) = (U \otimes I_{2 d}) {[κ_{1}^{T} (t), 0^{T}]}^{T}

(21)

{\tilde{κ}}_{e} (t) = (U \otimes I_{2 d}) {[0^{T}, κ_{e}^{T} (t)]}^{T}

(22)

Since $U \otimes I_{2 d}$ is nonsingular, one can obtain that ${\tilde{κ}}_{1} (t)$ and ${\tilde{κ}}_{e} (t)$ are linearly independent. Due to

(U^{T} \otimes I_{2 d}) ({\tilde{κ}}_{1} (t) + {\tilde{κ}}_{e} (t)) = {[κ_{1}^{T} (t), κ_{e}^{T} (t)]}^{T}

(23)

and

(U^{T} \otimes I_{2 d}) ω (t) = {[κ_{1}^{T} (t), κ_{e}^{T} (t)]}^{T}

(24)

one can see that $ω (t) = {\tilde{κ}}_{1} (t) + {\tilde{κ}}_{e} (t)$ . Then, it can be obtained that the multi-agent system (1) achieves formation control if and only if $lim_{t \to \infty} κ_{e} (t) = 0$ .

Theorem 1

For any given translation factor $γ > 0$ , multi-agent system (1) is said to achieve adaptive guaranteed-performance formation by protocol (10) if there exists a matrix P with $P = P^{T} > 0$ such that $2 Q + PA + A^{T} P - γ PB B^{T} P \leq 0$ . In this case, $K_{u} = B^{T} P$ , $K_{w} = PB B^{T} P$ . The guaranteed-performance cost is

\begin{matrix} J_{r}^{*} = ω^{T} (0) ((I_{N} - \frac{1}{N} 1_{N} 1_{N}^{T}) \otimes P) ω (0) \\ + γ \int_{0}^{+ \infty} ω^{T} (t) ((I_{N} - \frac{1}{N} 1_{N} 1_{N}^{T}) \otimes PB B^{T} P) ω (t) dt \end{matrix}

(25)

Proof

Design a Lyapunov function candidate as follows

\begin{matrix} V (t) = {κ_{e}}^{T} (t) (I_{N - 1} \otimes P) κ_{e} (t) \\ + \sum_{i = 1}^{N} \sum_{j \in N_{i}} \frac{{(w_{ij} (t) - w_{ij} (0))}^{2}}{2} \\ + \frac{γ}{2 N} \sum_{i = 1}^{N} \sum_{j = 1, j \neq i}^{N} (γ_{ij} - w_{ij} (t)) \end{matrix}

(26)

Due to $P = P^{T} > 0$ and $γ_{ij} \geq w_{ij} (t)$ , it can be shown that $V (t) \geq 0$ . Since, the topologic graph $G$ is undirected, one can obtain that $L_{w (t)}$ is a symmetric matrix. By equation (18), one has

\begin{matrix} \overset{\cdot}{V} (t) = {κ_{e}}^{T} (t) (I_{N - 1} \otimes (PA + A^{T} P)) κ_{e} (t) \\ - {κ_{e}}^{T} (t) (2 {\bar{U}}^{T} L_{w (t)} \bar{U} \otimes PB B^{T} P) κ_{e} (t) \\ + \sum_{i = 1}^{N} \sum_{j \in N_{i}} (w_{ij} (t) - w_{ij} (0)) {\overset{\cdot}{w}}_{ij} (t) \\ - \frac{γ}{2 N} \sum_{i = 1}^{N} \sum_{j = 1, j \neq i}^{N} {\overset{\cdot}{w}}_{ij} (t) \end{matrix}

(27)

Since $U$ is an orthonormal matrix, it can be obtained that $\bar{U} {\bar{U}}^{T} = L_{N}$ , where $L_{N}$ denotes the Laplacian matrix of a topologic graph $G^{*}$ and the weight of all the edges of $G^{*}$ is $1 / N$ . From protocol (10) and equation (16), one can show that

\begin{matrix} \sum_{i = 1}^{N} \sum_{j \in N_{i}} (w_{ij} (t) - w_{ij} (0)) {\overset{\cdot}{w}}_{ij} (t) - \frac{γ}{2 N} \sum_{i = 1}^{N} \sum_{j = 1, j \neq i}^{N} {\overset{\cdot}{w}}_{ij} (t) \\ = {κ_{e}}^{T} ((2 {\bar{U}}^{T} L_{w (t)} \bar{U} - 2 {\bar{U}}^{T} L_{w (0)} \bar{U} - γ I_{N - 1}) \otimes K_{w}) κ_{e} \end{matrix}

(28)

Due to $K_{w} = PB B^{T} P$ , it can be derived that

\begin{matrix} \overset{\cdot}{V} (t) = {κ_{e}}^{T} (I_{N - 1} \otimes (PA + A^{T} P)) κ_{e} \\ - {κ_{e}}^{T} ((2 {\bar{U}}^{T} L_{w (0)} \bar{U} - γ I_{N - 1}) \otimes K_{w}) κ_{e} \end{matrix}

(29)

Then, it means that

\overset{\cdot}{V} (t) \leq \sum_{i = 2}^{N} {κ_{i}}^{T} (t) (PA + A^{T} P - (2 λ_{i} + γ) PB B^{T} P) κ_{i} (t)

(30)

If $2 Q + PA + A^{T} P - γ PB B^{T} P \leq 0$ , then one has $PA + A^{T} P - γ PB B^{T} P < 0$ . Hence, one can obtain that $\overset{\cdot}{V} (t) \leq - ε ∥ κ_{e} ∥^{2}$ , where $ε$ is a positive constant.

Thus, $lim_{t \to \infty} κ_{e} (t) = 0$ ; that is, multi-agent system (1) achieves adaptively guaranteed-performance formation control by protocol (10). In the following, the guaranteed-performance cost $J_{r}^{*}$ is determined.

Due to $ω^{T} (t) (L_{N} \otimes I) ω (t) = {κ_{e}}^{T} (t) κ_{e} (t)$ , one gets

ω^{T} (t) (L_{N} \otimes Q) ω (t) = \sum_{i = 2}^{N} {κ_{i}}^{T} (t) Q κ_{i} (t)

(31)

Let $\nabla = \sum_{i = 1}^{N} \sum_{j = 1}^{N} {(ω_{j} (t) - ω_{i} (t))}^{T} Q (ω_{j} (t) - ω_{i} (t))$ and $h > 0$ , then one can obtain that

J_{r}^{h} \overset{Δ}{=} \frac{1}{N} \int_{0}^{h} \nabla dt = \sum_{i = 2}^{N} \int_{0}^{h} 2 {κ_{i}}^{T} (t) Q κ_{i} (t) dt

(32)

Due to $lim_{t \to + \infty} (γ_{ik} - w_{ik} (t)) = 0$ , one can infer that

lim_{h \to + \infty} \sum_{i = 1}^{N} \sum_{j = 1, j \neq i}^{N} (γ_{ik} - w_{ik} (h)) = 0

(33)

Due to $λ_{i} > 0$ , by equation (30), it holds that

\overset{\cdot}{V} (t) \leq \sum_{i = 2}^{N} {κ_{i}}^{T} (t) (PA + A^{T} P - γ PB B^{T} P) κ_{i} (t)

(34)

Hence, it can be derived that

\begin{matrix} J_{r}^{h} \leq \sum_{i = 2}^{N} \int_{0}^{h} 2 {κ_{i}}^{T} (t) Q κ_{i} (t) dt \\ \leq - \int_{0}^{h} \overset{\cdot}{V} (t) dt \\ = V (0) - V (h) \end{matrix}

(35)

Due to

lim_{h \to + \infty} V (h) = lim_{h \to + \infty} {κ_{e}}^{T} (h) (I_{N - 1} \otimes P) κ_{e} (h) + lim_{h \to + \infty} \sum_{i = 1}^{N} \sum_{j \in N_{i}} \frac{{(w_{ij} (h) - w_{ij} (0))}^{2}}{2} > 0

(36)

Then, it can be shown that

lim_{h \to + \infty} J_{r}^{h} \leq lim_{h \to + \infty} V (0) = {κ_{e}}^{T} (0) (I_{N - 1} \otimes P) κ_{e} (0) + \frac{γ}{2 N} \sum_{i = 1}^{N} \sum_{j = 1, j \neq i}^{N} (γ_{ij} - w_{ij} (0))

(37)

Since $\bar{U} {\bar{U}}^{T} = L_{N}$ and $κ_{e} (t) = [0_{2 (N - 1) d \times 2 d}, I_{2 (N - 1) d}] (U^{T} \otimes I_{2 d}) ω (t)$ , it can be obtained that

{κ_{e}}^{T} (0) (I_{N - 1} \otimes P) κ_{e} (0) = ω^{T} (0) ((I_{N} - \frac{1}{N} 1_{N} 1_{N}^{T}) \otimes P) ω (0)

(38)

Due to $lim_{t \to + \infty} (γ_{ik} - w_{ik} (t)) = 0$ , one can conclude that

\sum_{i = 1}^{N} \sum_{j = 1, j \neq i}^{N} (γ_{ik} - w_{ik} (0)) = \sum_{i = 1}^{N} \sum_{j = 1, j \neq i}^{N} \int_{0}^{+ \infty} {\overset{\cdot}{w}}_{ik} (t) dt = 2 N \int_{0}^{+ \infty} ω^{T} (t) (L_{N} \otimes K_{w}) ω (t) dt

(39)

From Theorem 1, $K_{u} = B^{T} P$ and $K_{w} = PB B^{T} P$ . $P$ can be determined by LMI techniques if $γ$ and $Q$ are given.

Introduce a positive constant $μ$ such that $μ I \geq P$ , then one can obtain that $PB B^{T} P \leq μ^{2} B B^{T}$ if $λ_{max} (B B^{T}) \leq 1$ , where $λ_{max} (B B^{T})$ denotes maximum eigenvalues of $B B^{T}$ . By LMI techniques, an adaptive guaranteed-performance formation control criterion based on Schur complement Lemma in Boyd et al.³⁵ is proposed as follows.

Corollary 1

For any given $μ > 0$ , multi-agent system (1) achieves adaptive guaranteed-performance formation by protocol (10) if $λ_{max} (B B^{T}) \leq 1$ and there exist $γ > 0$ and ${\tilde{P}}^{T} = \tilde{P} \geq μ^{- 1} I$ such that

[\begin{matrix} A \tilde{P} + \tilde{P} A^{T} - γ B B^{T} & 2 \tilde{P} Q \\ - 2 Q \end{matrix}] < 0

(40)

In this case, $K_{u} = B^{T} {\tilde{P}}^{- 1}$ , $K_{w} = {\tilde{P}}^{- 1} B B^{T} {\tilde{P}}^{- 1}$ and

J_{r}^{*} = \sum_{i = 2}^{N} (μ {‖ κ_{i} (0) ‖}^{2} + γ μ^{2} \int_{0}^{+ \infty} {‖ B^{T} κ_{i} (t) ‖}^{2} dt)

(41)

Corollary 2

If multi-agent system (1) achieves adaptive guaranteed-performance formation by protocol (10), then

lim_{t \to \infty} (x_{i} (t) - ι_{i}) = lim_{t \to \infty} (e^{At} (\frac{1}{N} \sum_{i = 1}^{N} (x_{i} (0) - ι_{i})))

(42)

Remark 2

The translation-adaptive strategy is proposed in Theorem 1 to solve the guaranteed-performance constraints of the adaptive consensus strategy in Li et al.³⁶ We use the special properties of the Laplacian matrix of the entire topology to achieve eigenvalue translation, which means all eigenvalues equivalent.

Numerical simulation

The following is a numerical simulation example to demonstrate the effectiveness of the theoretical results.

Consider a second-order mobile robot swarms with five agents, which are labeled from 1 to 5. The dynamics of each agent can be described by multi-agent system (1) and

A = [\begin{matrix} 0 & I_{2} \\ 0 & 0 \end{matrix}], B = [\begin{matrix} 0 \\ I_{2} \end{matrix}]

Let $Q = 0.8 I_{4}$ , $γ = 5$ , then it can be obtained that

K_{u} = [\begin{matrix} 1.8861 & 0 & 2.0304 & 0 \\ 0 & 1.8861 & 0 & 2.0304 \end{matrix}]

K_{w} = [\begin{matrix} 3.5574 & 0 & 3.8295 & 0 \\ 0 & 3.5574 & 0 & 3.8295 \\ 3.8295 & 0 & 4.1225 & 0 \\ 0 & 3.8295 & 0 & 4.1225 \end{matrix}]

and $J_{r}^{*} = 21.7675$

Let $C_{i} = (x_{i} (t), y_{i} (t), v_{xi} (t), v_{yi} (t))$ with $i \in {1, 2, 3, 4, 5}$ denote the state vector of the ith agent. The initial state vector $(x_{0 i}, y_{0 i}, v_{x 0 i}, v_{y 0 i})$ and formation vector $ι = {[ι_{1}^{T}, ι_{2}^{T}, ι_{3}^{T}, ι_{4}^{T}, ι_{5}^{T}]}^{T}$ with $ι_{i} = (ι_{xi}, ι_{yi})$ are shown in Table 1. The interaction topology of the mobile robot swarms is shown in Figure 1. From Definition 2, the second-order mobile robot swarms is said to achieve formation $ι$ if and only if

lim_{t \to \infty} (((x_{i} (t), y_{i} (t)) - (x_{j} (t), y_{j} (t))) = (ι_{i} - ι_{j})

(43)

and

lim_{t \to \infty} (v_{xi} (t), v_{yi} (t)) = lim_{t \to \infty} (v_{xj} (t), v_{yj} (t))

(44)

with i, j \in {1, 2, 3, 4, 5}

Table 1.

Initial state vector and formation vector.

Serial number	$(x_{0 i}, y_{0 i}, v_{x 0 i}, v_{y 0 i})$	$(ι_{xi}, ι_{yi})$
1	(–2, 1, 0, 0)	(–4, 3)
2	(–1, 2, 0, 0)	(–2, 3)
3	(0, 3,0, 0)	(0, 3)
4	(1, 2,0, 0)	(2, 3)
5	(2, 3, 0, 0)	(4, 3)

Figure 1.

Interaction topology of the mobile robot swarms.

The position state trajectories and velocity state trajectories of the mobile robot swarms are shown in Figures 2 and 3, respectively. One can see that position state trajectory of each agent converges to the corresponding value and the velocity states achieve consensus. Figure 4 shows that the interaction weights among neighboring agents converge to a finite value. Figure 5 shows that the performance function converges to a value less than $J_{r}^{*}$ . Figure 6 depicts the movement tracks of the robot system, one can see that the expected formation of the system is formed by the formation control protocol designed; that is, the system achieves adaptive guaranteed-performance formation control.

Figure 2.

Position state trajectories of the mobile robot swarms.

Figure 3.

Velocity state trajectories of the mobile robot swarms.

Figure 4.

Interaction weights among neighboring agents.

Figure 5.

Performance function of the mobile robot swarms.

Figure 6.

Movement tracks of the mobile robot swarms.

Conclusion

Adaptive guaranteed-performance formation control problems for second-order multi-agent systems were studied. An adaptive guaranteed-performance formation control strategy is proposed in the current article. The interaction weights of the system can be adaptively adjusted by the designed formation control protocol and the formation control gain can be regulated by choosing some proper translation factors.

Based on this result, it is meaningful to further study the formation control of practical multi-agent systems. Another important research topic is to address the adaptive guaranteed-performance formation control of multi-agent systems with leader-follower structure. For future work, we will focus on the time-varying guaranteed-performance formation control problem with switching topology.

Footnotes

Handling Editor: James Baldwin

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work is supported by National Natural Science Foundation (NNSF) of China under Grant Nos 61374054, 61503012, 61703411, 61503009, 61333011, and 61421063, in part by the Innovation Foundation of Rocket Force University of Engineering under Grant 2015ZZDJJ03, in part by the Innovation Zone Project under Grant 18-163-11-ZT-004-005-01, and in part by the Scientific Research Project under Grant JJ20172B03058.

ORCID iD

Xueli Xie

References

H-B

Wen

G-H

Cheng

Y-Y

, et al. Distributed finite-time cooperative control of multiple high-order nonholonomic mobile robots. IEEE T Neur Net Lear 2016; 28: 2998–3006.

Liu

Z-Q

Wang

Y-L

Wang

T-B

. Incremental predictive control-based output consensus of networked unmanned surface vehicle formation systems. Inform Sci 2018; 457: 166–181.

Johnson

Dasgupta

. Robust Hurwitz stability and performance analysis of H-infinity controlled forward-velocity dynamics of UAVs in close formation flight using bounded phase conditions in a Kharitonov framework. J Inst Eng India: Series C 2014; 95: 223–231.

Wang

Zhang

J-Y

Zhang

D-B

, et al. Development of an experiment system for UAV formation flight. Robot 2017; 39: 160–166, 175.

Cheng

Y-Y

H-B

Y-G

, et al. Distributed finite-time attitude regulation for multiple rigid spacecraft via bounded control. Inform Sci 2016; 328: 144–157.

Zhang

Y-H

H-Y

Sun

, et al. Cooperative adaptive event-triggered control for multiagent systems with actuator failures. IEEE T Syst Man Cy A 2019; 49: 1759–1768.

J-X

Wang

Liu

, et al. Completely distributed guaranteed-performance consensualization for high-order multiagent systems with switching topologies. IEEE T Syst Man Cy A 2018; 49: 1338–1348.

J-X

Liu

, et al. Admissible output consensualization control for singular multi-agent systems with time delays. J Frankl Inst 2016; 353: 4074–4090.

Tan

Z-Y

Cai

Zhou

, et al. On performance of peer review for academic journals: analysis based on distributed parallel system. IEEE Access 2019; 7: 19024–19032.

10.

J-X

Wang

Yang

X-G

, et al. Limited-budget output consensus for descriptor multiagent systems with energy constraints. IEEE T Cybernetics 2019, http://arxiv.org/abs/1909.08345

11.

Cai

Q-X

, et al. On almost controllability of dynamical complex networks with noises. J Syst Sci Complex 2019; 32: 1125–1139.

12.

J-X

Fan

Z-L

Liu

, et al. Guaranteed-cost consensus for multiagent networks with Lipschitz nonlinear dynamics and switching topologies. Int J Robust Nonlin 2018; 28: 2841–2852.

13.

H-Y

Jia

Cai

, et al. Adaptive guaranteed-performance consensus control for multi-agent systems with an adjustable convergence speed. Discrete Dyn Nat Soc 2019; 2019: 5190301.

14.

Z-J

H-S

. A new perspective to graphical characterization of multiagent controllability. IEEE T Cybernetics 2016; 47: 1471–1483.

15.

Liu

X-Z

Z-J

Hou

. Stabilization of heterogeneous multiagent systems via harmonic control. Complexity 2018; 2018: 8265637.

16.

Zhang

Y-H

Sun

Liang

H-J

, et al. Event-triggered adaptive tracking control for multiagent systems with unknown disturbances. IEEE T Cybernetics. Epub ahead of print 27 September 2018. DOI: 10.1109/TCYB.2018.2869084

17.

Yang

X-J

W-Z

Wang

, et al. Fast spectral clustering learning with hierarchical bipartite graph for large-scale data. Pattern Recogn Lett. Epub ahead of print 22 June 2018. DOI: 10.1016/j.patrec.2018.06.024

18.

Consolini

Morbidi

Prattichizzo

, et al. Leader-follower formation control of nonholonomic mobile robots with input constraints. Automatica 2008; 44: 1343–1349.

19.

Balch

Arkin

. Behavior-based formation control for multirobot teams. IEEE T Robotic Autom 1998; 14: 926–939.

20.

Lewis

Tan

K-H

. High precision formation control of mobile robots using virtual structures. Auton Robots 1997; 4: 387–403.

21.

Liu

Lewis

, et al. Robust formation control for multiple quadrotors with nonlinearities and disturbances. IEEE T Cybernetics. Epub ahead of print 25 October 2018. DOI: 10.1109/TCYB.2018.2875559

22.

Thunberg

Goncalves

X-M

. Consensus and formation control on SE (3) for switching topologies. Automatica 2016; 66: 109–121.

23.

Peng

Z-X

Wen

G-G

Rahmani

, et al. Distributed consensus-based formation control for multiple nonholonomic mobile robots with a specified reference trajectory. Int J Syst Sci 2015; 46: 1447–1457.

24.

Cepeda-Gomez

Olgac

. Stability of formation control using a consensus protocol under directed communications with two time delays and delay scheduling. Int J Syst Sci 2016; 47: 433–449.

25.

Wang

Liu

G-B

J-X

, et al. Guaranteed cost formation control for multi-agent systems: consensus approach. In: Proceedings of the 2015 34th Chinese control conference (CCC), Hangzhou, China, 28–30 July 2015, pp.7309–7314. New York: IEEE.

26.

Cheng

Wang

H-L

Hou

Z-G

, et al. Reaching a consensus in networks of high-order integral agents under switching directed topologies. Int J Syst Sci 2016; 47: 1966–1981.

27.

W-F

Liu

Feng

. Output consensus of heterogeneous linear multi-agent systems by distributed event-triggered/self-triggered strategy. IEEE T Cybernetics 2016; 47: 1914–1924.

28.

Cheng

Y-Y

Jia

R-T

, et al. Robust finite-time consensus formation control for multiple nonholonomic wheeled mobile robots via output feedback. Int J Robust Nonlin 2018; 28: 2082–2096.

29.

, et al. Energy consumption optimization for the formation of multiple robotic fishes using particle swarm optimization. Energies 2018; 11: 082023.

30.

Wang

J-X

, et al. Robust time-varying formation design for multi-agent systems with disturbances: extended-state-observer method. Int J Robust Nonlin 2019, https://arxiv.org/abs/1909.08974

31.

Cheng

Ugrinovskii

. Guaranteed performance leader-follower control for multi-agent systems with linear IQC coupling. In: Proceedings of the 2013 American control conference, Washington, DC, 17–19 June 2013, pp.2625–2630. New York: IEEE.

32.

Xue

R-B

Song

J-M

Zhang

M-Q

, et al. Research on distributed multi-vehicle coordinated formation flight control with coupling time-delay and jointly-connected topologies. Acta Armamentarii 2015; 36: 492–502.

33.

Savkin

Wang

Baranzadeh

, et al. Distributed formation building algorithms for groups of wheeled mobile robots. Robot Auton Syst 2016; 75: 463–474.

34.

Godsil

Royle

. Algebraic graph theory. Cham: Springer, 2013.

35.

Boyd

Ghaoui

Feron

, et al. Linear matrix inequalities in system and control theory. Philadelphia, PA: Siam, 1994.

36.

Z-K

Wen

G-H

Duan

Z-S

, et al. Designing fully distributed consensus protocols for linear multi-agent systems with directed graphs. IEEE T Automat Contr 2014; 60: 1152–1157.