## A Novel Edge Laplacian-based Approach for Adaptive Formation Control of Uncertain Multi-agent Systems with Unified Relative Error Performance

Kun Li, Kai Zhao, Member, IEEE , Yongduan Song, Fellow, IEEE , and Lihua Xie, Fellow, IEEE

Abstract -For most existing prescribed performance formation control methods, performance requirements are not directly imposed on the relative states between agents but on the consensus error, which lacks a clear physical interpretation of their solution. In this paper, we propose a novel adaptive prescribed performance formation control strategy, capable of guaranteeing prescribed performance on the relative errors, for uncertain high-order multi-agent systems under a class of directed graphs. Due to the consideration of performance constraints for relative errors, a coupled nonlinear interaction term that contains global graphic information among agents is involved in the error dynamics, leading to a fully distributed control design more difficult and challenging. Here by proposing a series of nonlinear mappings and utilizing the edge Laplacian along with Lyapunov stability theory, the presented formation control scheme exhibits the following appealing features when compared to existing results: 1) different performance requirements can be guaranteed in a unified way by solely tuning the design parameters a priori, without the need for control redesign and stability reanalysis under the proposed fixed control protocol, making the design more user-friendly and the implementation less demanding; 2) the complex and burdensome verification process for the initial constraint, often encountered in existing prescribed performance controls, is completely obviated if the performance requirements are global; and 3) nonlinear interaction is completely decoupled and the asymptotic stability of the formation manifold is ensured via using the adaptive parameter estimate technique. Finally, simulations of various performance behaviors are performed to show the efficiency of the theoretical results.

response. In addition, for a platoon of vehicles to cooperatively transport a common load, the distances between any two vehicles in the team must be neither too small nor too large, as this can lead to collisions between the vehicles or a breakdown of the link.

Index Terms -Adaptive control, edge Laplacian, formation control, prescribed performance.

## I. INTRODUCTION

Cooperative control of multi-agent systems (MASs) has attracted considerable attention in recent decades due to its broad application potential. Representative problems include consensus [1], formation control [2], and converge control [3], to name a few. Wherein, formation control, the goal of which is to coordinate each agent using local information so that the whole team forms a prescribed spatial geometric pattern [4], has been widely studied as it is applicable to various practical tasks, such as transporting loads, exploring resources, and environmental monitoring [5].

In addition to the classical formation control objectives, MASs may also need to meet some additional steady and transient behaviors due to performance requirements or safety constraints in physical, e.g., completion time, maximum allowable overshoot, tolerable range of tracking accuracy, maintaining connectivity, and collision avoidance [6]. If we ignore such performance requirements, it may lead to the failure of cooperative formation tasks. Let us provide two specific examples to illustrate the necessity of performance and safety constraints in formation control. Consider the multi-robot hand used in neurosurgery, a large overshoot or slow response can present a risk to life. Therefore, the controller used in the robotic hand should not only have a desired steady-state response but also a desired transient

- K. Li, K. Zhao, and Y. D. Song are with School of Automation, Chongqing University, Chongqing 400044 China (e-mail: likun/zhaokai/ydsong@cqu.edu.cn).
- L. H. Xie is with School of Electrical and Electronic Engineering, Nanyang Technological University, 639798 Singapore (e-mail: elhxie@ntu.edu.sg).

Current efforts on such issues are primarily based on the artificial potential function [7], barrier Lyapunov function [8], funnel control (FC) [9]-[11], and prescribed performance bound (PPB) control [12][15]. Notably, only the FC and PPB methodologies can be used to guarantee the prescribed performance behavior of the MASs in terms of overshoot, convergence rate and steady-state error. Yet, it is essential to emphasize that, whereas such performance specifications are obtained, these FC and PPB-based works only take into account the performance of the sum of relative errors between agent i and all its neighbors (consensus errors) and not the relative errors between agent i and each neighbor, which does not reflect the inter-agent physical restrictions. Taking the cooperative formation control as an example, it is indeed more practical to specify the performance constraints on the relative position of agents, so that the interest of collision avoidance and connectivity maintenance can be tackled. On the other hand, from the user's perspective, it is more intuitive to define constraint functions based on the relative errors between agents, as the constraint requirements are normally shaped by the environmental boundaries.

In recent years, some typical results on the performance characteristics of relative errors between agents have been developed. In [16] and [17], consensus protocols with certain prescribed performance guarantees were proposed for single- and double-integrator dynamics, respectively. In [18], a decentralized formation control protocol with connectivity maintenance was presented for first-order nonlinear MASs. In [19], the relative position-based formation control for leader-follower MASs with prescribed performances was studied. In [20], a distributed dynamic average consensus algorithm was proposed for MASs under undirected communication topologies, so that the consensus with predefined and arbitrary convergence rate is guaranteed. Recently, a novel connectivity preservation and collisionfree formation control algorithm for double-integrator dynamics was proposed in [21].

By reviewing the above performance-constrained works [16][21], we can conclude that 1) the results only consider the singleor double-integrator dynamics under undirected graphs, which is impractical in the real applications since many engineering systems suffer from the high-order dynamics and modeling uncertainties; and 2) there exists an initial constraint on the relative errors in the algorithm implementation, which implies that the aforementioned methods only guarantee the semi-global results. To achieve the performance requirements, it is necessary to look for suitable initial values of system states and then to verify the initial constraint. However, due to the complexity and large-scale network structure of MASs, the verification process will become more and more demanding and grow explosively as the number of agents increases (the detailed discussion can be found in Section V-A). Moreover, when the considered system is restarted or the target formation is

changed, we have to check whether the original design parameters are applicable to the new initial data. If such a condition is not met, the corresponding control is no longer available, leading to the inflexibility of the control method. These limitations pose numerous obstacles to their application, which in turn motivates this work.

Motivated by the above observations, we will investigate the relative-error-performance formation problem for uncertain highorder MASs under directed graphs. As mentioned earlier, the difficulties in solving these problems are mainly in two aspects, 1) how to design a fully distributed formation protocol with prescribed performance guarantees for uncertain high-order MASs under directed graphs, and 2) how to ensure different performance behaviors between an agent and its neighbors within a fixed control framework according to its own environmental constraints. To address the above challenges, we present a unified prescribed performance solution for the formation control problem of uncertain high-order MASs. Compared with the existing results, the main contributions of this paper can be listed as follows:

- 1) Plug and Play Function: The proposed framework features that different performance requirements for the relative errors can be guaranteed by tuning the design parameters a priori and are completely decoupled from the system model, topology graph, and control gain selection. It provides a new formulation that achieves different prescribed performance for relative position errors of uncertain high-order nonlinear MASs in a unified way, without changing the control structure. Especially, the complex and explosive verification process for the initial constraints in existing prescribed performance controls is completely obviated if the performance requirement is in the global form;
- 2) User-Oriented Feature: Compared with the consensus errorbased PPB schemes [12]-[15], the proposed solution enables user to establish performance constraint boundaries directly based on environmental constraints, making it more favorable in practical applications. This physics-based performance control design approach provides a physical interpretation of the proposed approach; and
- 3) Broad Applicability: We propose an adaptive formation control scheme for uncertain high-order nonlinear MASs under a class of directed and undirected graphs. In this respect, existing methods in [16]-[21] that consider simple dynamics under undirected graphs are contained as special cases of our work.

Notations: In this paper, || · || is the Euclidean norm. The sets of real, positive real, and n -dimensional Euclidean space are R R , + , and R n , respectively. Let 1 n be unit column vector with n -dimensional. X ⊤ denotes the transpose of a matrix X . For X i ∈ R m n × , i = 1 , . . . , q , we denote diag [ X i ] = blockdiag { X ,. . . , X 1 q } ∈ R qm × qn . For a set of numbers V , | V | is its cardinality. Null ( A ) stands for the null space of matrix A λ . min ( A ) and λ max ( A ) denote the minimum and maximum eigenvalue of matrix A respectively.

## II. PRELIMINARIES

## A. Graph &amp; Edge Laplacian

We consider a team of N agents whose information exchange is described by a directed graph G ≜ ( V E , ) , where V = { 1 , ..., N } is the set of nodes and E ⊆ V × V with cardinality m is the set of edges. An edge ( i, j ) ∈ E indicates that agent j can obtain information from agent i , but not necessarily vice versa. A directed path in a directed graph is a sequence of edges, for example, ( i 1 , i 2 ) , ( i 2 , i 3 ) , . . . , ( i N -1 , i N ) , and a directed cycle is a directed path where the initial and final nodes are the same. If ( i, j ) ∈ E implies ( j, i ) ∈ E , the graph is said to be undirected . The adjacency matrix A = [ a ij ] ∈ R N × N with a ij = 1 , if

( i, j ) ∈ E , and a ij = 0 , otherwise. We denote the set of neighbors for agent i as N i = { j ∈ V| ( i, j ) ∈ E} . The degree matrix is △ = diag [ △ i ] ∈ R N × N with △ i = ∑ j ∈N i a ij being the sum of the i th row in A . Then, the graph Laplacian matrix of G is defined as L = △-A .

To ensure the specified performance for relative errors in the formation control, we choose an alternative graph-theoretic tool, i.e., the edge Laplacian [22]. Now we recall this concept as follows.

Label all edges of the digraph G with e , e 1 2 , . . . , e m . Suppose edge ( i, j ) in G corresponds to the k th directed edge, where k ∈ I = { 1 , . . . , m } . The incidence matrix E ( G ) ∈ R N × m for the digraph is a { 0 , ± } 1 -matrix with rows and columns indexed by nodes and edges of G , respectively, such that

$$\text{$n$} \text{$s$.} \\ \text{is}$$

From the definition of the incidence matrix, it follows that the null space of its transpose, Null ( E ( G ) ⊤ ) , contains span { 1 N } . Moreover, for the purpose of analysis, it is necessary to partition the incidence matrix as E ( G ) = E ⊙ ( G ) + E ⊗ ( G ) , where E ⊙ ( G ) ∈ R N × m corresponds to the so-called in-incidence matrix, whose elements are defined as

$$\iota _ { 1 }, \quad [ E _ { \odot } ] _ { i k } \coloneqq \begin{cases} - 1, & \text{if $i$ is the terminal node of edge $e_{k}$,} \\ 0, & \text{otherwise.} \end{cases}$$

and E ⊗ ( G ) ∈ R N × m corresponds to the so-called out-incidence matrix, whose elements are defined as

$$[ E _ { \otimes } ] _ { i k } \coloneqq \begin{cases} + 1, & \text{if $i$ is the initial node of edge $e_{k}$,} \\ 0, & \text{otherwise.} \end{cases}$$

According to the definitions of incidence matrix and in-incidence matrix, the edge Laplacian matrix L e ∈ R m m × and graph Laplacian L ∈ R N × N of a digraph G can be defined in terms of the incidence and in-incidence matrices as L e = E E ⊤ ⊙ and L = E E ⊙ ⊤ . For an undirected graph, L e = E E ⊤ and L = EE ⊤ , respectively. Moreover, using an appropriate permutation of the edge ordering, the incidence matrix is expressed as E = [ E t E c ] , where E t ∈ R N × ( N -1) denotes the full-column-rank incidence matrix corresponding to an arbitrary spanning tree G t ⊂ G and E c ∈ R N × ( m N -+1) represents the incidence matrix corresponding to the remaining edges not contained in G t . Defining R := [ I N -1 T ] , where T := ( E E ⊤ t t ) -1 E E ⊤ t c , one obtains representation of the incidence matrix of the graph, given by E = E R t . Furthermore, inspired by [23], the following facts will play a crucial role in the controller design and stability analysis presented later in this paper.

Lemma 1: [23] For a digraph that is a spanning tree, it holds that L s e = 1 2 ( E E ⊤ ⊙ + E E ⊤ ⊙ ) is positive definite; For a digraph that is a directed cycle, it holds that E E ⊤ t t is positive definite.

## B. Problem Statement

In this paper, we consider a group of N nonlinear agents that can be described as:

$$\begin{cases} \dot { x } _ { i, q } = x _ { i, q + 1 }, \ q = 1, \dots, n - 1, \\ \dot { x } _ { i, n } = u _ { i } + \varphi _ { i } ^ { \top } ( \overline { x } _ { i } ) \theta _ { i }, \ i = 1, \dots, N, \end{cases}$$

where x i,q ∈ R and u i ∈ R are the system state and control input of agent i , respectively, where we only consider the one dimension case in the theoretical analysis without loss of generality. The discussions in this work can be generalized to higher dimensions cases with Kronecker product. x i := [ x ⊤ i, 1 , x ⊤ i, 2 , · · · , x ⊤ i,n ] ⊤ ∈ R n ,

φ i : R n → R ν i × n is a known smooth nonlinear function, and θ i ∈ R ν i represents an unknown constant vector, where ν i ∈ R + .

Let e k := x i, 1 -x j, 1 ∈ R be the relative position between a pair of neighboring agents i and j for ( i, j ) ∈ E , and let e d k ∈ R denote the desired relative position for e k . The target relative positionbased formation is assumed to be admissible and is described as F := { | e e k = e d k , k ∈ I} , where e = E x ⊤ 1 ∈ R m denotes the stack vector of the relative states, and x 1 = [ x 1 1 , , x 2 1 , , . . . , x N, 1 ] ⊤ ∈ R N denotes the stack vector of the first-order components. Then we define the formation errors as ˜ e k = e k -e d k ∈ R for k ∈ I . It can be checked that the target formation is achieved if e k = 0 . The control objective in this paper is formally stated as follows.

Control Objective: By using the local available sensing measurements from the agent itself and its neighbor agents, we aim to design a distributed adaptive formation controller u i for (1) such that:

- O 1 : All signals in the closed-loop system are bounded, and the desired formation, described by the relative offsets { e d k } k ∈I , is achieved and maintained; and
- O 2 : Multiple prescribed performance requirements for the relative position errors { ˜ e k } k ∈I are guaranteed in a unified control framework without control redesign.

To this end, the following assumption is imposed.

Assumption 1: The sensing graph should either be a connected undirected graph or a directed graph that is a directed spanning tree or a directed cycle.

## III. SYSTEM AND PERFORMANCE TRANSFORMATION

In this section, we formulate the formation performance control problem by constructing a unified performance function. Motivated in part by a previous work [24], we propose a more comprehensive unified performance function to enhance the applicability of the algorithm.

## A. Unified Performance Function

Definition 1: Let P : ( -ι, ι ) → -∞ ∞ ( , ) be a unified performance function for ι ∈ R + , which satisfies the following conditions.

- 1) lim y → ι P ( y ) = ∞ , lim y →-ι P ( y ) = -∞ , and P (0) = 0 ;
- 2) The derivative of P is lower bounded by a positive constant q and tends to infinity as y →± ι ; and
- 3) P is continuously differentiable and P -( y ) = -P ( y ) .

Obviously, immediate examples of such P ( y ) include (but are not limited to) the following functions: P ( y ) = y √ 1 -y 2 with ι = 1 and P ( y ) = tan( y ) with ι = π 2 . Utilizing the unified performance function definition, the second objective O 2 can be mathematically stated as:

$$\mathcal { P } _ { k } ( - \underline { \delta } _ { k } \beta _ { k } ( t ) ) < \tilde { e } _ { k } ( t ) < \mathcal { P } _ { k } ( \overline { \delta } _ { k } \beta _ { k } ( t ) ) \text{ for } k \in I, \quad ( 2 )$$

where 0 &lt; δ ,δ k k ≤ 1 are user-chosen parameters. Additionally, β k ( ) t is a continuously decaying performance positive function with bounded derivative and a strictly positive limit as t →∞ , which is formulated as β k := ( β k 0 -β kf ) exp( -λ t k ) + β kf with 0 &lt; β kf &lt; β k 0 ≤ ι and λ k ∈ R + .

It can be readily seen that by choosing different design parameters δ k , δ k , and β k 0 , such a performance objective as defined in (2) covers different cases of interest, including lower/upper bounded onesided performance constraints, asymmetric performance constraints, and global performance constraints, for details see below.

- 1) If δ k β k 0 = δ k β k 0 = ι , then (2) can represent global performance constraints in the form of -∞ &lt; e ˜ (0) k &lt; ∞ ;
- 2) If δ k β k 0 = ι and 0 &lt; δ k &lt; 1 (resp. δ k β k 0 = ι and 0 &lt; δ k &lt; 1 ), then (2) can represent lower bounded one-sided

performance constraints in the form of -ϵ k &lt; e ˜ (0) k &lt; ∞ (resp. upper bounded one-sided performance constraints in the form of -∞ &lt; e ˜ (0) k &lt; ϵ k ), where ϵ k = P k ( -δ k β k 0 ) and ϵ k = P ( δ k β k 0 ) are some positive constants; and

- 3) If 0 &lt; δ β k k 0 &lt; ι and 0 &lt; δ β k k 0 &lt; ι , then (2) can also represent asymmetric performance constraints in the form of -ϵ k &lt; e ˜ (0) k &lt; ϵ k .

Overall, the designer has the flexibility to choose various design parameters (i.e., δ k , δ k , and β k 0 ) based on practical requirements. It is interesting to underscore that in the case of asymmetric performance constraints, the existing relative error PPB-based results [16]-[21] are relevant. Their algorithms are special cases of the algorithm developed in this paper.

## B. Transformed Error Dynamic Model

To deal with the constraints in (2), we will adopt the following nonlinear mapping:

$$s _ { k } \coloneqq \frac { \zeta _ { k } ( t ) } { ( \underline { \delta } _ { k } + \zeta _ { k } ( t ) ) ( \overline { \delta } _ { k } - \zeta _ { k } ( t ) ) },$$

$$\zeta _ { k } \coloneqq \frac { \eta _ { k } ( t ) } { \beta _ { k } ( t ) }, \, \eta _ { k } \coloneqq \mathcal { P } _ { k } ^ { - 1 } ( \tilde { e } _ { k } ( t ) )$$

with ζ k ( ) t denoting the modulated error and η k ( ) t being the normalized error. Furthermore, the following lemma can be easily derived from the aforementioned transformation.

Lemma 2: If ζ k (0) satisfies -δ k &lt; ζ k (0) &lt; δ k and s k ( ) t is bounded for ∀ t ≥ 0 , then there exist some positive constraints δ 1 k and δ 1 k so that -δ k &lt; -δ 1 k ≤ ζ k ( ) t ≤ δ 1 k &lt; δ k . Furthermore, (2) can be guaranteed.

Proof. The proof of this lemma is similar to [24], therefore it is omitted herein. ■

It is thus obvious that if an adaptive control law is designed such that s k ( ) t is bounded over t ∈ [0 , ∞ ) , the target O 2 for the tracking performance is achieved according to the Lemma 2. To facilitate the descriptions in control design and stability analysis, the unified performance function P ( y ) in this work is chosen as P ( y ) = y/ √ 1 -y 2 . Upon using the expression of s k as shown in (3), with the definitions of η k and β k in (4), it is seen that

$$\dot { s } _ { k } = \frac { \partial s _ { k } } { \partial \zeta _ { k } } \frac { \mathrm d \zeta _ { k } } { \mathrm d t } = J _ { k } \left ( \xi _ { k } \dot { \tilde { e } } _ { k } - \dot { \beta } _ { k } \zeta _ { k } \right ), \quad \quad ( 5 )$$

where J k = µ k β k &gt; 0 , µ k = δ k δ k + ζ 2 k ( δ k + ζ k ) 2 ( δ k -ζ k ) 2 &gt; 0 , ξ k = 1 ( 1+˜ e 2 ) √ ˜ e 2 +1 &gt; 0 are available for control design.

k k Based on the definition of the formation error ˜ e k , the transformed error dynamics for system (5) are rewritten in the following compact form

$$\dot { S } = J \left ( \xi \dot { \tilde { e } } + \dot { \beta } \zeta \right ) = W E ^ { \top } x _ { 2 } + D, \text{ \quad \ \ } \text{(6)}$$

where S = [ s , . . . , s 1 m ] ⊤ ∈ R m , J = diag [ J k ] ∈ R m m × , ξ = diag [ ξ k ] ∈ R m m × , ˙ β = diag [ -˙ β k ] , ζ = [ ζ 1 , . . . , ζ m ] ⊤ , W = Jξ , and D = Jβζ ˙ . Furthermore, given that each argument of W is positive, it follows that the diagonal matrix W is positive definite.

By replacing the equation of ˙ x i, 1 in (1) with ˙ s k in (5), the performance constraint objective is transformed to the stabilization problem of the following strict-feedback-like system:

$$\begin{cases} \lambda &. \\ \dot { S } = W E ^ { \top } x _ { 2 } + D, \\ \dot { x } _ { q } = x _ { q + 1 }, \ q = 2, \dots, n - 1, \\ \dot { x } _ { n } = u + \varphi ^ { \top } \theta, \end{cases}$$

where x q = [ x 1 ,q , x 2 ,q , . . . , x N,q ] ⊤ , u = [ u , u 1 2 , . . . , u N ] ⊤ , φ = diag [ φ ⊤ 1 , φ ⊤ 2 , . . . , φ ⊤ N ] , and θ = [ θ , θ 1 2 , . . . , θ N ] ⊤ .

## IV. CONTROLLER DESIGN AND STABILITY ANALYSIS

## A. Controller Design

In this subsection, based on the adaptive backstepping technique [25] we propose a novel adaptive formation design approach by utilizing the edge Laplacian. Firstly, we introduce the following change of coordinates:

$$z _ { 1 } = S, & & ( 8 ) & \underset { \text{defini} } { \dots \dots }$$

$$z _ { q } = x _ { q } - \alpha _ { q - 1 }, q = 2, \dots, n, \quad \quad ( 9 ) \quad T h \epsilon$$

where z 1 = [ s , . . . , s 1 m ] ⊤ and z q = [ z 1 ,q , . . . , z N,q ] ⊤ are intermediate variables, and α q -1 = [ α 1 ,q -1 , . . . , α N,q -1 ] ⊤ is the virtual controller, respectively. The design procedure contains n steps.

Step 1: From (1), (6) and (9), the derivative of z 1 is

$$\dot { z } _ { 1 } = W E ^ { \top } ( z _ { 2 } + \alpha _ { 1 } ) + D, \text{ \quad \ \ } ( 1 0 ) \text{ \quad when }$$

then the derivative of V 1 = 1 2 z ⊤ 1 z 1 along (10) yields,

$$\dot { V } _ { 1 } = z _ { 1 } ^ { \top } \left [ W E ^ { \top } ( z _ { 2 } + \alpha _ { 1 } ) + D \right ]. \quad \quad ( 1 1 ) \stackrel { \ell } { c _ { 2 } } > 1$$

Remark 1: We pause here to emphasize that, due to the consideration of relative error performance behaviors, the coupled nonlinear interaction term WE ⊤ and the additional coupling term Jβζ ˙ have been incorporated into the derivative of V 1 as presented in (11). This introduces an obstacle in designing the virtual controller α 1 based on local information to ensure the controllability of the subsystem dynamics (10). The main challenges arise from the following three aspects:

- 1) The transformed error dynamics, as described in (5), are significantly different from the conventional PPC control methodologies [12]-[15]. This distinction stems from the inclusion of a nonlinear term ξ k , which is coupled with the formation errors, leading to difficulties in ensuring the controllability of subsystems (10);
- 2) The term WE ⊤ represents a non-square matrix that contains global topological information, making the direct use of its inverse impractical for controller design; and
- 3) The third challenge arises from the fact that in the setup for directed graphs, the Laplacian matrix L is non-symmetric and positive semi-definite. Thus, it is challenging to ensure the controllability of the subsystem dynamics (10) using Laplacian matrix-based design methods.

To handle the above challenges, we perform the controller design of α 1 using the edge Laplacian technique. The virtual controller α 1 is designed as

$$\alpha _ { 1 } = - c _ { 1 } E _ { \odot } W z _ { 1 }, \ c _ { 1 } > 0. \quad \quad ( 1 2 ) \quad \dot { V } _ { 1 }$$

By substituting (12) into (11), the derivative of V 1 is given by

$$\dot { V } _ { 1 } = - c _ { 1 } z _ { 1 } ^$$

Remark 2: We highlight that each component of α 1 depends only on local information since E ⊙ represents the incoming edges on each node, i.e., the information available to each agent as defined by the directed graph. Moreover, equation (13) is applicable to both directed spanning-tree, directed-cycle and undirected graph topologies. In the following, we proceed by separately examining each of these topologies.

Case 1 (Directed spanning tree): In this case, we have G t = G . Therefore, E t = E R , = I N -1 , and E ⊙ t = E ⊙ , where E ⊙ t is the in-incidence matrix of the spanning tree. It is not difficult to obtain z 1 t = z 1 since R = I N -1 . Then we have

$$\dot { V } _ { 1 } & = - c _ { 1 } z _ { 1 t$$

where the third inequality follows from the fact that L s e is positive definite according to Lemma 1.

Then after using D = Jβζ ˙ in (6) and applying Young's inequality to last two term in the right-hand side of (14), we obtain

z

⊤

⊤

ϑ

⊤

2

1

2

1

t

WE z

t

2

≤

2

λ

max

(

E E

t

t

)

∥

Wz

1

t

∥

+ 2

ϑ

∥

z

2

∥

,

(15)

$$z _ { 1 t } ^ { \top } J \dot { \beta } \zeta \$$

where ϑ and ϵ are positive constants.

Now, given c 1 &gt; 0 , let ϑ &gt; 0 be such that c ′ 1 := [ c 1 λ min ( L s e ) -1 2 ϑλ max ( E E ⊤ t t ) ] λ min ( W 2 ) and c ′ 2 = c 2 -1 2 ϑ with c 2 &gt; 0 are positive, we bound the (14) as follows:

$$\dot { V } _ { 1 } & \leq - c _ { 1 } ^ { \$$

where c ′′ 1 := c ′ 1 -ϵ 2 λ max ( J 2 ) is positive for a sufficiently small ϵ &gt; 0 .

Case 2 (Directed cycle): In this case, we have E E ⊤ ⊙ + E E ⊤ ⊙ = E E ⊤ . By using this fact, we obtain

$$\begin{array} {$$

where the third inequality follows from the fact that E E ⊤ t t is positive definite according to Lemma 1.

Now, given c 1 &gt; 0 , let ϑ &gt; 0 be such that c ′ 1 := 1 2 [ c 1 λ min ( E E ⊤ t t ) -ϑλ max ( E E ⊤ t t ) ] λ min ( RWR RWR ⊤ ⊤ ) and c ′ 2 = c 2 -1 2 ϑ with c 2 &gt; 0 are positive. Note that RWR RWR ⊤ ⊤ and RJ R 2 ⊤ are positive definite due to R is full row-rank and W and J are positive definite. Then, similar to the case 1, by applying Young's inequality to last two term in the right-hand side of (18), we obtain

$$( 1 2 ) \quad & \dot { V } _ { 1 } \leq - c _ { 1 } ^ { \prime } \| z _ { 1 t } \| ^ { 2 } + \frac { 1 } { 2 \vartheta } \| z _ { 2 } \| ^ { 2 } + \frac { \epsilon } { 2 } z _ { 1 t } ^ { \top } R J J R ^ { \top } z _ { 1 t } + \frac { 1 } { 2 \epsilon } \zeta ^ { \top } \dot { \beta } \dot { \beta } \zeta \\ \text{by} \quad & \leq - c _ { 1 } ^ { \prime \prime } \| z _ { 1 t } \| ^ { 2 } + \frac { 1 } { 2 \vartheta } \| z _ { 2 } \| ^ { 2 } + \frac { 1 } { 2 \epsilon } \zeta ^ { \top } \dot { \beta } \dot { \beta } \zeta, \quad & ( 1 9 ) \\ \text{$1$} \infty$$

where c ′′ 1 := c ′ 1 -ϵ 2 λ max ( RJ R 2 ⊤ ) is positive for a sufficiently small ϵ &gt; 0 .

Case 3 (Undirected graph): In this setup, each agent can not only receive information from its neighbors but also broadcast its own information. Therefore, the virtual controller can be designed as follows:

$$\alpha _ { 1 } = - c _ { 1 } E W z _ { 1 },$$

where E represents the incidence matrix of the oriented undirected graph.

By substituting (20) into (11), the derivative of V 1 is given by

$$\nu _ { 1 } \, \text{souquoting} \, \infty \, \text{in} \, \infty \, \text{in} \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \incolon \\ \dot { V } _ { 1 } & = - c _ { 1 } z _ { 1 } ^ { \top } W E ^ { \top } E W z _ { 1 } + z _ { 1 } ^ { \top } W E ^ { \top } z _ { 2 } + z _ { 1 } ^ { \top } D \\ & = - c _ { 1 } z _ { 1 t } ^ { \top } R W R ^ { \top } E _ { t } ^ { \top } E _ { t } R W R ^ { \top } z _ { 1 t } + z _ { 1 } ^ { \top } W E ^ { \top } z _ { 2 } + z _ { 1 } ^ { \top } D \\ & \leq - c _ { 1 } \lambda _ { \min } ( E _ { t } ^ { \top } E _ { t } ) \| R W R ^ { \top } z _ { 1 t } \| ^ { 2 } + z _ { 1 t } ^ { \top } R W R ^ { \top } E _ { t } ^ { \top } z _ { 2 } \\ & \quad + z _ { 1 } ^ { \top } J \dot { \beta } \zeta, \quad ( 2 1 ) \quad & \quad \text{by apply} \,. \\ \dots \dots \, \quad \pi ^ { \top } \, \pi \, \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \} \,.$$

where E E ⊤ t t represents the edge Laplacian of a spanning tree G t ⊂ G and is symmetric positive definite.

It is seen from (21) that ˙ V 1 has a similar expression. Therefore, by employing the same analysis as in Case 2, we can easily obtain:

$$\dot { V } _ { 1 } \leq - c _ { 1 } ^ { \prime \prime } \| z _ { 1 t } \| ^ { 2 } + \frac { 1 } { 2 \vartheta } \| z _ { 2 } \| _ { 2 } ^ { 2 } + \frac { 1 } { 2 \epsilon } \zeta ^ { \top } \dot { \beta } \dot { \beta } \zeta, \quad ( 2 2 ) \quad \text{Herc}$$

where c ′′ 1 := c ′ 1 -ϵ 2 λ max ( RJ R 2 ⊤ ) &gt; 0 and c ′ 1 := [ c 1 λ min ( E E ⊤ t t ) -1 2 ϑλ max ( E E ⊤ t t ) ] λ min ( RWR RWR ⊤ ⊤ ) &gt; 0 .

In summary, we see from (17), (19), and (22) that the expression of ˙ V 1 is same for both the considered directed graph and undirected graph. With this in mind, we proceed with the controller design.

Step 2: Using (9), the derivative of z 2 is computed as

$$\dot { z } _ { 2 } = \dot { x } _ { 2 } - \dot { \alpha } _ { 1 } = z _ { 3 } + \alpha _ { 2 } - \dot { \alpha } _ { 1 }. \quad \quad ( 2 3 ) \ \text{soluti} \.$$

From (11) and (23), the derivative of V 2 = V 1 + 1 2 z ⊤ 2 z 2 is

$$\dot { V } _ { 2 } & \leq - \, c _ { 1 } ^ { \prime \prime } \| z _ { 1 t } \| ^ { 2 } + \frac { 1 } { 2 \vartheta } \| z _ { 2 } \| ^ { 2 } + \frac { 1 } { 2 \epsilon } \zeta ^ { \top } \dot { \beta } \dot { \beta } \zeta & \quad \text{remain} \\ & \quad + z _ { 2 } ^ { \top } \left ( z _ { 3 } + \alpha _ { 2 } - \dot { \alpha } _ { 1 } \right ). & \quad \text{(24)} \quad \text{that $\tau_{n}$}$$

Therefore, we design α 2 as

$$\alpha _ { 2 } = - c _ { 2 } z _ { 2 } + \dot { \alpha } _ { 1 }, \ c _ { 2 } > 0. \quad \quad ( 2 5 ) \ \eta _ { k } / \beta _ {, }$$

By replacing α 2 from (25) in (24), one has

$$\dot { V } _ { 2 } \leq - c _ { 1 } ^ { \prime \prime } \| z _ { 1 t } \| ^ { 2 } - c _ { 2 } ^ { \prime } \| z _ { 2 } \| ^ { 2 } + z _ { 2 } ^ { \top } z _ { 3 } + \frac { 1 } { 2 \epsilon } \zeta ^ { \top } \dot { \beta } \dot { \beta } \zeta. \quad ( 2 6 ) \quad \text{where} \ \, \jmath$$

Step q ( q = 3 , . . . , n -1) : Define the quadratic function V q as V q = V q -1 + 1 2 z ⊤ q z q , then its derivative is

$$\cdot \alpha \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \\ \dot { V } _ { q } & \leq - c _ { 1 } ^ { \prime \prime } \| z _ { 1 t } \| ^ { 2 } - c _ { 2 } ^ { \prime } \| z _ { 2 } \| ^ { 2 } - \sum _ { j = 3 } ^ { q - 1 } c _ { j } z _ { j } ^ { \top } z _ { j } + z _ { q - 1 } ^ { \top } z _ { q } & \quad \begin{cases} \dot { z } _ { 1 } \\ \dot { z } _ { q } \\ \dot { \tilde { \theta } } \end{cases} \\ \text{Following} \, the ideas in the above stems. we design $\alpha$ as$$

Following the ideas in the above steps, we design α q as

$$\alpha _ { q } = - c _ { q } z _ { q } - z _ { q - 1 } + \dot { \alpha } _ { q - 1 }, \ c _ { q } > 0. \quad \ ( 2 8 ) \ \ z _ { n } ( t ).$$

Substituting the virtual control α q into (27), one has

$$\dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \delta _ { \zeta } = \\ \dot { V } _ { q } \leq & - c _ { 1 } ^ { \prime \prime } \| z _ { 1 t } \| ^ { 2 } - c _ { 2 } ^ { \prime } \| z _ { 2 } \| ^ { 2 } - \sum _ { j = 3 } ^ { q - 1 } c _ { j } z _ { j } ^ { \top } z _ { j } \\ & + z _ { q } ^ { \top } z _ { q + 1 } + \frac { 1 } { 2 \epsilon } \zeta ^ { \top } \dot { \beta } \dot { \beta } \zeta. & ( 2 9 ) \begin{matrix} & \text{one do} \\ & \text{one do} \\ & \text{h} ( t, Z _ {. } \\ \\ \text{ten n} \text{: The derivative of $z$} \text{, is commuted as} \end{matrix}$$

Step n: The derivative of z n is computed as

$$\dot { z } _ { n } = \dot { x } _ { n } - \dot { \alpha } _ { n - 1 } = u + \varphi ^ { \top } \theta - \dot { \alpha } _ { n - 1 }. \quad \ \ ( 3 0 ) \quad \text{(Some)}$$

Define the Lyapunov function candidate V as

$$V = V _ { n - 1 } + \frac { 1 } { 2 } z _ { n } ^ { \top } z _ { n } + \frac { 1 } { 2 } \tilde { \theta } ^ { \top } \Gamma ^ { - 1 } \tilde { \theta }, \quad \quad ( 3 1 ) \quad _ { - - }$$

where ˜ = θ θ -ˆ θ is the estimate error of θ with ˆ θ being the estimate of θ , and Γ = diag [Γ ] i ∈ R N × N is a positive definite matrix with Γ i &gt; 0 . Then the derivative of V is

$$\dot { V } & \leq - c _ { 1 } ^ { \prime \prime } \| z _ { 1 t } \| ^ { 2 } - c _ { 2 } ^ { \prime } \| z _ { 2 } \| ^ { 2 } - \sum _ { j = 3 } ^ { n - 1 } c _ { j } z _ { j } ^ { \top } z _ { j } + z _ { n - 1 } ^ { \top } z _ { n } & \quad \text{where } \rho \\ & \quad + \frac { 1 } { 2 \epsilon } \zeta ^ { \top } \dot { \beta } \dot { \beta } \zeta + z _ { n } ^ { \top } \left ( u + \varphi ^ { \top } \theta - \dot { \alpha } _ { n - 1 } \right ) - \tilde { \theta } ^ { \top } \Gamma ^ { - 1 } \dot { \hat { \theta } }. \quad ( 3 2 ) \quad \dot { V } \leq \cdot$$

The actual controller and the adaptive law are designed as:

$$u = - c _ { n } z _ { n } - z _ { n - 1 } + \dot { \alpha } _ { n - 1 } - \varphi ^ { \top } \hat { \theta }, \ c _ { n } > 0, \quad ( 3 3 )$$

˙

ˆ = Γ

θ

φz .

n

(34)

By applying the actual control and adaptive law into (32), we have

$$\Sigma _ { \quad } \dot { V } & \leq - c _ { 1 } ^ { \prime \prime } \| z _ { 1 t } \| ^ { 2 } - c _ { 2 } ^ { \prime } \| z _ { 2 } \| ^ { 2 } - \sum _ { j = 3 } ^ { n - 1 } c _ { j } z _ { j } ^ { \top } z _ { j } + \frac { 1 } { 2 \epsilon } \zeta ^ { \top } \dot { \beta } \dot { \beta } \zeta. \quad ( 3 5 )$$

## B. Stability Analysis

Here we present our main result in the following theorem.

Theorem 1: Consider the uncertain nonlinear MASs described by (1) under Assumption 1. If the developed control (33) and adaptive law (34) are employed, then target formation can be achieved, ensuring the satisfaction of multiple performance requirements, as well as the boundedness of all closed-loop signals at all time.

Proof. We carry out the proof of Theorem 1 by three phases. Firstly, we address the existence and uniqueness of a maximal solution ζ k ( ) t over the open set Ω ζ k := { ζ k ( ) : t ζ k ( ) t ∈ ( -δ k , δ k ) } for a time interval [0 , τ max ) with τ max being a positive constant. Secondly, we prove that all internal signals are bounded and ζ k ( ) t remains strictly within a subset of Ω ζ k for t ∈ [0 , τ max ) under the control law (33) and the adaptive control law (34). Finally, we show that τ max can be extended to ∞ and subsequently ζ k ( ) t ∈ -( δ k , δ k ) is ensured over time, thus completing the proof.

Phase I. Note that the modulated error ζ k ( ) t is given by ζ k = η /β k k . Differentiating ζ k with respect to time yields:

$$\dot { \zeta } = h _ { 0 } ( t, \zeta ) = \overline { \beta } \left ( \xi E ^ { \top } ( z _ { 2 } - c _ { 1 } E _ { \odot } W z _ { 1 } ) + \dot { \beta } \zeta \right ), \quad ( 3 6 )$$

where β = diag [ 1 β 1 , . . . , 1 β N ] . Differentiating z q and ˜ θ , employing (12), (25), (28) as well as the fact that z q = x q -α q -1 , and invoking (33), (34), we arrive at

$$\begin{cases} \dot { z } _ { 1 } = h _ { 1 } ( t, z _ { 2 }, \zeta ) = W E ^ { \top } ( z _ { 2 } - c _ { 1 } E _ { \odot } W z _ { 1 } ) + D, \\ \dot { z } _ { q } = h _ { q } ( t, z _ { q + 1 } ) = z _ { q + 1 } - c _ { q } z _ { q }, \ q = 2, \dots, n - 1, \\ \dot { z } _ { n } = h _ { n } ( t, z _ { n }, \tilde { \theta } ) = - c _ { n } z _ { n } - z _ { n - 1 } + \varphi ^ { \top } \tilde { \theta }, \\ \dot { \tilde { \theta } } = h _ { n + 1 } ( t, z _ { n } ) = - \Gamma \varphi z _ { n }. \end{cases} \quad \cdots \quad \cdots \quad \cdots$$

Thus, the closed-loop dynamic system of Z t ( ) = [ ζ ( ) t , z 1 ( ) t , . . . , z n ( ) t , θ ˜ ( t )] ⊤ may be written in compact form as:

$$\dot { Z } = h ( t, Z ) = \left [ h _ { 0 } ^ { \top } ( t, \zeta ), h _ { 1 } ^ { \top } ( t, z _ { 2 }, \zeta ), \cdots, h _ { n + 1 } ^ { \top } ( t, z _ { n } ) \right ] ^ { \top }. \ \ ( 3 8 )$$

Let us also define the open set Ω Z = Ω ζ × R n +1 , where Ω ζ = Ω ζ 1 × Ω ζ 2 × · · · × Ω ζ m . Note that ζ k (0) ∈ ( δ k , δ k ) , one observes Z (0) = [ ζ (0) , z 1 (0) , . . . , z n (0) , θ ˜ (0)] ⊤ ∈ Ω Z . Since h t, Z ( ) is continuous in t and locally Lipschitz in Z over the set Ω Z , the hypotheses of Theorem 54 in [26] (see p. 476) are satisfied. Consequently, the ˙ Z -system has a unique and maximal solution over the time interval [0 , τ max ) , ensuing that ζ k ( ) t ∈ Ω ζ k for all t ∈ [0 , τ max ) , i.e.,

$$\zeta _ { k } ( t ) \in ( - \underline { \delta } _ { k }, \overline { \delta } _ { k } ), \ k \in I, \ \forall t \in [ 0, \tau _ { \max } ). \quad \ ( 3 9 )$$

Phase II. According to (39), one has

$$\| \zeta \| ^ { 2 } < \sum _ { k = 1 } ^ { m } \max \{ \underline { \delta } _ { k } ^ { 2 }, \bar { \delta } _ { k } ^ { 2 } \} = \rho, \text{ \quad \quad } ( 4 0 )$$

where ρ is a positive constant. Then it follows from (35) that

$$( 3 2 ) \quad \dot { V } \leq - c _ { 1 } ^ { \prime \prime } \| z _ { 1 t } \| ^ { 2 } - c _ { 2 } ^ { \prime } \| z _ { 2 } \| ^ { 2 } - \sum _ { j = 3 } ^ { n - 1 } c _ { j } z _ { j } ^ { \top } z _ { j } + \frac { \rho } { 2 \epsilon ^ { 2 } } \lambda _ { \max } ( \dot { \beta } ^ { 2 } ) \ ( 4 1 )$$

for ∀ t ∈ [0 , τ max ) . Integrating both sides of (41) in the interval [0 , τ max ) yields that

$$V ( t ) + & \int _ { 0 } ^ { \tau _ { \max } } \left ( c _ { 1 } ^ { \prime \prime } z _ { 1 t } ^ { \top } z _ { 1 t } + c _ { 2 } ^ { \prime } z _ { 2 } ^ { \top } z _ { 2 } + \sum _ { j = 3 } ^ { n } c _ { j } z _ { j } ^ { \top } z _ { j } \right ) d t & \int _ { e _ { 2 } } ^ { \tau _ { \max } } \\ & \leq V ( 0 ) + \frac { \rho } { 2 \epsilon ^ { 2 } } \int _ { 0 } ^ { \tau _ { \max } } \lambda _ { \max } ( \dot { \beta } ^ { 2 } ) d t. & ( 4 2 ) \\ \varsigma \dots \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots.$$

Given that ˙ β k is smooth and lim t →∞ ˙ β k = 0 , it follows that lim t →∞ λ max ( ˙ β 2 ) = 0 , implying that ∫ τ max 0 λ max ( ˙ β 2 ) dt is bounded. Based on the definition of V in (31), we establish that z 1 t and z q ( q = 2 , . . . , n ) and ˜ θ are bounded and integrable for all t ∈ [0 , τ max ) . Consequently, the parameter estimate ˆ θ remains bounded over the interval [0 , τ max ) . Furthermore, the boundedness of z q and ˜ e ensures that x i,q is also bounded, leading to the boundedness of φ i ( x i ) and ˙ α q ( q = 1 , . . . , n -1 ) within the same interval. As β k and p ∗ are continuously differentiable, ˙ z q ( q = 1 , . . . , n ) is also bounded. Then from (33) and (34), we conclude that u and ˙ ˆ θ are bounded for all t ∈ [0 , τ max ) . Therefore, all signals in the closedloop systems are bounded over the interval [0 , τ max ) . Additionally, due to the boundedness of s k , according to the Lemma 2, we have

$$- \underline { \delta } _ { k } < - \underline { \delta } _ { 1 k } \leq \zeta _ { k } ( t ) \leq \overline { \delta } _ { 1 k } < \overline { \delta } _ { k }, k \in I, t \in [ 0, \tau _ { \max } ). \ \ ( 4 3 ) \ \text{phase} ^ { \cdot }$$

Phase III. Note that from (43), ζ k ( ) t ∈ Ω ′ ζ k := { ζ k ( ) t ∈ [ δ 1 k , δ 1 k ] } for ∀ t ∈ [0 , τ max ) is a nonempty compact subset of Ω ζ k . Hence, assuming τ max ≤ ∞ and since Ω ′ ζ k ⊂ Ω ζ k , Proposition C.3.6 in [26] (see p. 481) dictates the existence of a time instant t ′ ∈ [0 , τ max ) such that ζ k ( t ′ ) / ∈ Ω ′ ζ k , which is a clear contradiction. Therefore, τ max = ∞ . Thus, all closed loop signals remain bounded and moreover ζ k ∈ Ω ′ ζ k ⊂ Ω ζ k , ∀ t &gt; 0 . Based on Lemma 2, (2) is guaranteed for ∀ t ≥ 0 . In addition, by applying the Barbalat's lemma, we have lim t →∞ z 1 t = 0 . Note that z 1 t = 0 implies S = 0 since R is full rank, which further implies ζ = 0 and η = 0 according to the nonlinear mapping as shown in (3) and (4), and hence ˜ e k = 0 for k ∈ I . The proof is completed. ■

Remark 3: It is worth noting that although most existing works (e.g., [12]-[15]) can guarantee the asymmetric prescribed performance for the virtual consensus error, it remains unclear that whether the relative error satisfies the performance constraint. Moreover, if the users would like to achieve other different performances, one has to redesign the controller and reanalyze the stability of the closedloop system, making the control implementation more complex and less flexible. Some efforts (e.g., [16]-[21]) have been developed to achieve the prescribed performance for the relative errors, but it still only guarantees one special constraint form. To solve the above issue, in this paper, we introduce a unified performance function (as defined in Definition 1) and a series of nonlinear mappings (as shown in (3) and (4)), so that the prescribed performance problem of relative errors is converted into the stabilization problem of a new strict-feedback-like system. Under such a framework, multiple kinds of prescribed performance behaviors can be guaranteed only by solely tuning the design parameters a priori, without the need for control redesign and stability reanalysis, making the design more user-friendly and the implementation less demanding. Especially for the global performance requirement case, the tedious and explosive verification process for the initial constraint in the existing methods is completely eliminated.

Remark 4: In this work, due to the considerations of multiple performance constraints for relative errors and the digraph setup, a coupled nonlinear interaction term WE ⊤ that contains global graphic information among agents and an additional coupling term Jβζ ˙ are inevitably embedded into the error dynamics (7), therefore, a mismatched term appearing in (24) may exist in each step of the

Fig. 1. The sensing graph for a group of 5 mobile robots: (a) directed spanning tree; (b) directed cycle; and (c) undirected graph.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00323v1-figure-b124638d-e8a0-4745-a245-66827320053c.jpg)

backstepping design procedure, which poses an additional challenge for the stability analysis and control design. To overcome these difficulties, inspired by [22], we propose an edge Laplacian-based adaptive control method, which exhibits the following appealing characteristics. Firstly, by skillful designing the first virtual controller α 1 , the nonlinear interconnections term WE ⊤ can be treated as a quadratic function W L W ⊤ e as shown in (13), then with the aid of the reduced-order representation and the eigenvalue analysis, we utilize the adaptive parameter estimate technique to establish the asymptotic stability of the formation manifold. Secondly, here we need to emphasize that multiple performance behaviors can be achieved without changing the control structure because the performance index has been completely decoupled from the structure of the dynamic systems and the underlying graphs. Furthermore, to handle the negative effects caused by the mismatched terms z ⊤ 1 WE z ⊤ 2 and z ⊤ 1 Jβζ ˙ , the robust technique is used for further analysis as shown in (15) and (16). As such, the stability of the closed-loop system can be ensured via utilizing the Lyapunov stability.

## V. SIMULATION RESULTS

To illustrate the effectiveness of the proposed formation controller, we consider a group of 5 mobile robots described by [7] interconnected over the directed spanning-tree/directed cycle/undirected graph showed in Fig. 1, moving in the plane (i.e., d = 2 ) with the following dynamics:

$$\begin{cases} \dot { x } _ { i } ( t ) = v _ { i } ( t ), \ i = 1, \dots, 5, \\ \dot { v } _ { i } ( t ) = u _ { i } ( t ) - \varphi _ { i } ^ { \top } ( x _ { i }, v _ { i } ) \theta _ { i }, \end{cases}$$

where x i = [ x i, 1 , x i, 2 ] ⊤ ∈ R 2 and v i = [ v i, 1 , v i, 2 ] ⊤ ∈ R 2 are, respectively, the position and velocity of the mobile robots, φ ⊤ i ( x , v i i ) θ i is an unknown friction satisfying φ i = diag [tanh(10 v i, 1 ) -tanh(100 v i, 1 ) , tanh(10 v i, 2 ) -tanh(100 v i, 2 )] and θ i = [0 25 0 25] . , . ⊤ .

The desired formation corresponds to a pentagon described by p d i ( ) t = [cos(2( i -1) π/ 5) , sin(2( i -1) π/ 5)] ⊤ for i = 1 , . . . , 5 , and the desired relative position vector is denoted as e d k = [ e d kx , e d ky ] ⊤ for k = 1 , . . . , 5 . For Fig. 1(a) (directed spanning tree case), e d k is set to e d 1 = p d 2 -p d 1 , e d 2 = p d 3 -p d 2 , e d 3 = p d 4 -p d 3 , and e d 4 = p d 5 -p d 4 . For Fig. 1(b) (directed cycle case), e d k is set to e d 1 = p d 2 -p d 1 , e d 2 = p d 3 -p d 2 , e d 3 = p d 4 -p d 3 , e d 4 = p d 5 -p d 4 , and e d 5 = p d 5 -p d 1 ; For Fig. 1(c) (undirected graph case), e d k is set to e d 1 = p d 2 -p d 1 , e d 2 = p d 3 -p d 2 , e d 3 = p d 4 -p d 3 , e d 4 = p d 5 -p d 4 , and e d 5 = p d 1 -p d 5 . In the following, to show the validity and the advantages of the proposed algorithm, we will verify the asymmetric and global performance behaviors just by choosing different design parameters under a fixed controller. To simplify notations, we choose the identical performance function for each edge, implying that each agent has the same communication/sensing capabilities, namely,

$$\overset { \prime \, \mu } { \dots } \, \ a \, \ a$$

$$\cdot \text{ of the } \text{ for } k = 1, \dots, 5, \text{ where } \beta ( t ) = ( \beta _ { 0 } - \beta _ { f } ) \exp ( - \lambda t ) + \beta _ { f }.$$

Fig. 2. Asymmetric performance behavior. (a) directed spanning tree. (b) directed cycle, and (c) undirected graph cases.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00323v1-figure-5690cd96-f132-4b27-b7c7-c45f603f8d51.jpg)

## A. Asymmetric Performance Behavior

To verify that the proposed control is able to achieve the asymmetric prescribed tracking performance. In this part, for the directed spanning tree case, the agents' initial positions are set as x 1 1 , (0) = [1 , 1 2] . ⊤ , x 2 1 , (0) = [0 6 . , 2 3] . ⊤ , x 3 1 , (0) = [ -0 8 . , 2 1] . ⊤ , x 4 1 , (0) = [ -1 4 . , 0 5] . ⊤ , and x 5 1 , (0) = [ -0 2 . , 0 35] . ⊤ . For the directed cycle and undirected graph case, the agents' initial positions are set as x 1 1 , (0) = [1 2 2] . , ⊤ , x 2 1 , (0) = [0 5 3] . , ⊤ , x 3 1 , (0) = [ -0 8 . , 2 5] . ⊤ , x 4 1 , (0) = [ -0 9 . , 1 5] . ⊤ , and x 5 1 , (0) = [0 4 1] . , ⊤ ; the initial velocities for each agent are set as v i (0) = zeros (2 1) , . To fulfill the initial constraint, the key design parameters δ , δ and β 0 are selected as δ = 0 3 . and δ = β 0 = 0 8 . such that each initial relative error falls within the initial prescribed boundary, i.e., -0 25 . = P -( δβ 0 ) &lt; e ˜ kq (0) &lt; P ( δβ 0 ) = 0 83 . , q ∈ { x, y } and k = 1 , . . . , 5 . The other parameters are chosen as: λ = 1 , c 1 = 1 , c 2 = 5 , Γ i = I 2 , and β f = 0 1 . . By implementing adaptive control law (33) and (34), the simulation results are shown in Fig. 2. From Fig. 2(top) it is observed that the five agents eventually achieve and maintain the target formation shape during their motions. The responses of edge errors are depicted in the middle and bottom of Fig. 2, respectively, from which it can be seen that in all cases-(a) directed spanning tree, (b) directed cycle, and (c) undirected graph-the edge errors asymptotically converge to zero within the prescribed performance bounds over time.

It should be noted that in order to guarantee the asymmetric performance behavior of the relative errors, it is crucial to carefully choose the key design parameters δ , δ and β 0 to meet the initial constraints (i.e., P -( δβ 0 ) &lt; e ˜ kq (0) &lt; P ( δβ 0 ) ). This situation also exists in [12]-[21]. However, this parameter selection process is exceedingly tedious and time-consuming, for example, in an N -agent MAS with m edges, it necessitates the design and execution of a performance boundary selection algorithm at least m times to identify feasible parameters. To circumvent such challenges, this work proposes a solution by considering global performance characteristics, as discussed below.

## B. Global Performance Behavior

In this part, the agents' initial positions are set as x 1 1 , (0) = [1 6 . , -0 2] . ⊤ , x 2 1 , (0) = [1 2 . , -0 1] . ⊤ , x 3 1 , (0) = [1 75 . , 0] ⊤ , x 4 1 , (0) = [1 4 . , 0 15] . ⊤ , and x 5 1 , (0) = [1 5 . , -0 1] . ⊤ , and the initial velocities are also set as v i (0) = zeros (2 1) , . The key design parameters δ , δ , and β 0 are selected as δ = δ = β 0 = 1 , i.e., -∞ = P -( δβ 0 ) &lt; e ˜ kq (0) &lt; P ( δβ 0 ) = ∞ , q ∈ { x, y } , which shows that there is no any constraint on the initial states. The other parameters are chosen as: λ = 1 , c 1 = 2 , c 2 = 5 , Γ i = I 2 , and β f = 0 1 . . By implementing adaptive control law (33) and (34), the simulation results are shown in Fig. 3. From Fig. 3(top) it can be seen that the five agents eventually achieve and maintain the target formation shape during their motions. Figs. 3 middle and bottom show that the edge errors asymptotically converge to zero within the

Fig. 3. Global performance behavior. (a) directed spanning tree. (b) directed cycle, and (c) undirected graph cases.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00323v1-figure-04744fba-3ea1-44dd-b3c6-6ed0a0f48c42.jpg)

prescribed boundary, which indicates that the performance behaviors are independent of initial conditions, illustrating the effectiveness of the proposed approach.

- [2] W. Ren and E. Atkins, 'Distributed multi-vehicle coordinated control via local information exchange,' International Journal of Robust and Nonlinear Control: IFAC-Affiliated Journal , vol. 17, no. 10-11, pp. 10021033, 2007.

## VI. CONCLUSION

In this paper, a novel adaptive formation control algorithm has been proposed for a group of agents described by uncertain highorder nonlinear dynamic models. The features that distinguish our method from the existing ones include: 1) the relative error constraints between the agents with physical meaning are taken into account, which makes the developed algorithm advantageous for practical applications; 2) it is able to fulfill various performance requirements merely by selecting design parameters a priori, making the control redesign and stability reanalysis not required; and 3) the dynamic model and underlying graph considered in this work are more general such that most of existing results are the special cases of the proposed control. The simulation results have shown the efficacy of our approach. Future research may consider multiple time-varying performance constraints and derive a formation control algorithm in the local reference frame.

## REFERENCES

- [1] R. Olfati-Saber and R. M. Murray, 'Consensus problems in networks of agents with switching topology and time-delays,' IEEE Transactions on automatic control , vol. 49, no. 9, pp. 1520-1533, 2004.
- [3] M. Zhong and C. G. Cassandras, 'Distributed coverage control and data collection with mobile sensor networks,' IEEE Transactions on Automatic Control , vol. 56, no. 10, pp. 2445-2455, 2011.
- [4] K.-K. Oh, M.-C. Park, and H.-S. Ahn, 'A survey of multi-agent formation control,' Automatica , vol. 53, pp. 424-440, 2015.
- [5] Y. Yang, Y. Xiao, and T. Li, 'A survey of autonomous underwater vehicle formation: Performance, formation control, and communication capability,' IEEE Communications Surveys &amp; Tutorials , vol. 23, no. 2, pp. 815-841, 2021.
- [6] C. Wei, X. Wu, B. Xiao, J. Wu, and C. Zhang, 'Adaptive leaderfollowing performance guaranteed formation control for multiple spacecraft with collision avoidance and connectivity assurance,' Aerospace Science and Technology , vol. 120, p. 107266, 2022.
- [7] C. Sun, G. Hu, L. Xie, and M. Egerstedt, 'Robust finite-time connectivity preserving coordination of second-order multi-agent systems,' Automatica , vol. 89, pp. 21-27, 2018.
- [8] E. Restrepo, A. Lor´ ıa, I. Sarras, and J. Marzat, 'Robust consensus of high-order systems under output constraints: Application to rendezvous of underactuated uavs,' IEEE Transactions on Automatic Control , vol. 68, no. 1, pp. 329-342, 2022.
- [9] X. Min, S. Baldi, W. Yu, and J. Cao, 'Low-complexity control with funnel performance for uncertain nonlinear multi-agent systems,' IEEE Transactions on Automatic Control , 2023.
- [10] C. K. Verginis, 'Asymptotic consensus of unknown nonlinear multiagent systems with prescribed transient response,' in 2023 American Control Conference (ACC) . IEEE, 2023, pp. 1992-1997.
- [11] J. G. Lee, S. Trenn, and H. Shim, 'Synchronization with prescribed

transient behavior: Heterogeneous multi-agent systems under funnel coupling,' Automatica , vol. 141, p. 110276, 2022.

- [12] C. P. Bechlioulis and G. A. Rovithakis, 'Decentralized robust synchronization of unknown high order nonlinear multi-agent systems with prescribed transient and steady state performance,' IEEE Transactions on Automatic Control , vol. 62, no. 1, pp. 123-134, 2016.
- [13] C. K. Verginis, C. P. Bechlioulis, D. V. Dimarogonas, and K. J. Kyriakopoulos, 'Robust distributed control protocols for large vehicular platoons with prescribed transient and steady-state performance,' IEEE Transactions on Control Systems Technology , vol. 26, no. 1, pp. 299304, 2017.
- [14] W. Cheng, K. Zhang, and B. Jiang, 'Fixed-time fault-tolerant formation control for a cooperative heterogeneous multiagent system with prescribed performance,' IEEE Transactions on Systems, Man, and Cybernetics: Systems , vol. 53, no. 1, pp. 462-474, 2022.
- [15] W. Wang, D. Wang, Z. Peng, and T. Li, 'Prescribed performance consensus of uncertain nonlinear strict-feedback systems with unknown control directions,' IEEE Transactions on Systems, Man, and Cybernetics: Systems , vol. 46, no. 9, pp. 1279-1286, 2015.
- [16] Y. Karayiannidis, D. V. Dimarogonas, and D. Kragic, 'Multi-agent average consensus control with prescribed performance guarantees,' in 2012 IEEE 51st IEEE Conference on Decision and Control (CDC) . IEEE, 2012, pp. 2219-2225.
- [17] L. Macellari, Y. Karayiannidis, and D. V. Dimarogonas, 'Multi-agent second order average consensus with prescribed transient behavior,' IEEE Transactions on Automatic Control , vol. 62, no. 10, pp. 52825288, 2016.
- [18] C. P. Bechlioulis and K. J. Kyriakopoulos, 'Robust model-free formation control with prescribed performance and connectivity maintenance for nonlinear multi-agent systems,' in 53rd IEEE Conference on Decision and Control . IEEE, 2014, pp. 4509-4514.
- [19] F. Chen and D. V. Dimarogonas, 'Leader-follower formation control with prescribed performance guarantees,' IEEE Transactions on Control of Network Systems , vol. 8, no. 1, pp. 450-461, 2020.
- [20] C. J. Stamouli, C. P. Bechlioulis, and K. J. Kyriakopoulos, 'Robust dynamic average consensus with prescribed transient and steady state performance,' Automatica , vol. 144, p. 110503, 2022.
- [21] Y. Huang, Z. Meng, and D. V. Dimarogonas, 'Prescribed performance formation control for second-order multi-agent systems with connectivity and collision constraints,' Automatica , vol. 160, p. 111412, 2024.
- [22] D. Zelazo, A. Rahmani, and M. Mesbahi, 'Agreement via the edge laplacian,' in 2007 46th IEEE Conference on Decision and Control . IEEE, 2007, pp. 2309-2314.
- [23] E. Restrepo, A. Lor´ ıa, I. Sarras, and J. Marzat, 'Edge-based strict lyapunov functions for consensus with connectivity preservation over directed graphs,' Automatica , vol. 132, p. 109812, 2021.
- [24] K. Zhao, Y. Song, C. P. Chen, and L. Chen, 'Adaptive asymptotic tracking with global performance for nonlinear systems with unknown control directions,' IEEE Transactions on Automatic Control , vol. 67, no. 3, pp. 1566-1573, 2021.
- [25] M. Krstic, P. V. Kokotovic, and I. Kanellakopoulos, Nonlinear and adaptive control design . John Wiley &amp; Sons, Inc., 1995.
- [26] E. D. Sontag, Mathematical control theory: deterministic finite dimensional systems . Springer-Verlag New York, Inc., 1998.
Kun Li Kai Zhao Yongduan Song Lihua Xie
A Novel Edge Laplacian-based Approach for Adaptive Formation Control of Uncertain Multi-agent Systems with Unified Relative Error Performance
For most existing prescribed performance formation control methods,
performance requirements are not directly imposed on the relative states
between agents but on the consensus error, which lacks a clear physical
interpretation of their solution. In this paper, we propose a novel adaptive
prescribed performance formation control strategy, capable of guaranteeing
prescribed performance on the relative errors, for uncertain high-order
multi-agent systems under a class of directed graphs. Due to the consideration
of performance constraints for relative errors, a coupled nonlinear interaction
term that contains global graphic information among agents is involved in the
error dynamics, leading to a fully distributed control design more difficult
and challenging. Here by proposing a series of nonlinear mappings and utilizing
the edge Laplacian along with Lyapunov stability theory, the presented
formation control scheme exhibits the following appealing features when
compared to existing results: 1) different performance requirements can be
guaranteed in a unified way by solely tuning the design parameters a priori,
without the need for control redesign and stability reanalysis under the
proposed fixed control protocol, making the design more user-friendly and the
implementation less demanding; 2) the complex and burdensome verification
process for the initial constraint, often encountered in existing prescribed
performance controls, is completely obviated if the performance requirements
are global; and 3) nonlinear interaction is completely decoupled and the
asymptotic stability of the formation manifold is ensured via using the
adaptive parameter estimate technique. Finally, simulations of various
performance behaviors are performed to show the efficiency of the theoretical
results.