Event-Triggered Bounded Consensus Tracking for Second-Order Nonlinear Multi-Agent Systems with Uncertainties

Abstract

This paper is concerned with event-triggered bounded consensus tracking for a class of second-order nonlinear multi-agent systems with uncertainties (MASs). Remarkably, the considered MASs allow multiple uncertainties, including unknown control coefficients, parameterized unknown nonlinearities, uncertain external disturbances, and the leader’s control input being unknown. In this context, a new estimate-based adaptive control protocol with a triggering mechanism is proposed. We rule out Zeno behavior by testifying that the lower bound on the interval between two consecutive events is positive. It is shown that under the designed protocol, all signals caused by the closed-loop systems are bounded globally uniformly and tracking errors ultimately converge to a bounded set. The effectiveness of the devised control protocol is demonstrated through a simulation example.

1. Introduction

Nowadays, consensus tracking has been extensively investigated for MASs due to its widespread applications in many fields such as cooperative control of mobile robots and spacecraft formation flying control [1,2,3,4,5,6]. In MASs, each agent is usually equipped with embedded microprocessors, which have limited energy and computing resources. We remark that event-triggered control (ETC) makes it so the information transmitted or the controller is updated only when necessary for systems, which thus can effectively reduce resource consumption. As an effective tool for reducing resource consumption, event-triggered control has been proposed [7,8]. This inspires the study of the event-triggered consensus tracking problem of MASs.

In the past few decades, abundant work on consensus tracking of MASs with ETC has been reported [9,10,11,12,13,14,15,16,17,18,19]. Specifically, single-integrator agents and double-integrator agents were investigated in [9,10], respectively. However, the presented triggering mechanism required each agent to continuously monitor the states of neighboring agents. This requirement is removed in [11] by adopting a triggering threshold that is state-independent, while consensus tracking of general linear MASs was realized. As a further improvement, works [12,13] considered the consensus tracking problems by ETC for nonlinear MASs.

We remark that all the above works do not consider MASs with uncertainties, which are unavoidable when modeling real plants. Specifically, work [20] proposed a hybrid system approach to address the ETC problem for linear systems with uncertainties. Then, in [21], consensus tracking was realized for nonlinear MASs whose the parameters in nonlinearities have a known upper bound. As an extension, work [22] permitted completely unknown parameters in nonlinearities. Work [23] considered MASs with parameterized unknown nonlinearities, but unknown control coefficients were not taken into account. Work [24] considered uncertain nonlinear MASs with unknown control coefficients, and an estimate-based ETC protocol was proposed to achieve bounded consensus tracking. Then, uncertain nonlinear MASs subjected to unknown external disturbances were investigated in [25], and output consensus tracking was realized. However, in reality, second-order nonlinear MASs with uncertainties are more practical for modeling real systems such as torque motors and jet engines, which are tuned to achieve the desired motion directly through acceleration rather than velocity [26]. Hence, it is exigent to investigate consensus-tracking-based ETC for nonlinear MASs with uncertainties. In this context, how to design an ETC protocol to realize bounded consensus tracking and how to develop an appropriate compensation mechanism to counteract the uncertainties in systems deserve our efforts. A comparison of the features of the investigated MASs in the existing literature is depicted in

Table 1. Comparison of the features of the investigated MASs in the existing literature.

Work	Control Coefficient	Lipschitz Constant	External Disturbance	Leader Input
[21]	No	Known	No	Known
[22]	No	Unknown	No	Zero
[22]	No	Unknown	No	Zero
[23]	Yes	Unknown	No	Unknown
this paper	Yes	Unknown	Yes	Unknown

.

This paper is devoted to event-triggered bounded consensus tracking for second-order nonlinear MASs with uncertainties. Compared with the relevant literature, the MASs in this paper permit multiple uncertainties, including unknown control coefficients and parameterized unknown nonlinearities. In contrast with work [27], the considered nonlinear MASs with uncertainties permit unknown external disturbances, and event-triggered consensus tracking is investigated in this paper. Moreover, the leader’s control input is also allowed to be unknown, and the leader’s information is broadcast to only a few of the agents. To counteract the uncertainties and to realize bounded consensus tracking, a novel estimate-based adaptive control protocol with triggering mechanism is developed. Further, through Lyapunov analysis, it is proved that the devised adaptive control protocol can make certain that all signals caused by the closed-loop systems are bounded globally uniformly, and with the passage of time, tracking errors converge to a bounded set.

The remaining context of this paper consists of the following sections. The related preliminaries, which include the problem statement, graph theory, and notation, are formulated in Section 2. Section 3 proposes the triggering mechanism and a distributed adaptive control protocol with the triggering mechanism. The main results are summarized in Section 4. Simulations of nonlinear MASs with uncertainties are given to exemplify the effectiveness of the proposed control protocol are in Section 5. Concluding remarks are provided in Section 6.

2. Preliminaries

2.1. Graph Theory and Notation

Represent a graph by $\mathcal{G}=(\mathcal{V},\mathcal{E})$, which consists of a nonempty finite set of nodes and a set of edges $\mathcal{E}\subseteq \mathcal{V}\times \mathcal{V}$, where a set of nodes is denoted by $\mathcal{V}$ = (${\mathit{v}}_1$, …,${\mathit{v}}_{\mathit{N}}$). The edge (${\mathit{v}}_{\mathit{i}}$,${\mathit{v}}_{\mathit{j}}$) denotes that node ${\mathit{v}}_{\mathit{i}}$ can transmit information to node ${\mathit{v}}_{\mathit{j}}$, where ${\mathit{v}}_{\mathit{i}}$ is ${\mathit{v}}_{\mathit{j}}$’s neighbor. Graph $\mathcal{G}$ is called undirected if (${\mathit{v}}_{\mathit{i}}$,${\mathit{v}}_{\mathit{j}}$) $\in \mathcal{E}$ implies (${\mathit{v}}_{\mathit{j}}$,${\mathit{v}}_{\mathit{i}}$) $\in \mathcal{E}$. For a directed graph, an edge (${\mathit{v}}_{\mathit{i}}$,${\mathit{v}}_{\mathit{j}}$)$\in \mathcal{E}$ indicates that node ${\mathit{v}}_{\mathit{i}}$ can transmit information to node ${\mathit{v}}_{\mathit{j}}$, but node ${\mathit{v}}_{\mathit{j}}$ cannot transmit information to node ${\mathit{v}}_{\mathit{i}}$. A directed graph contains a spanning tree: the graph includes a root node with no parent and the root node has directed paths to every other node.

The adjacency matrix of a graph $\mathcal{G}$ is denoted by $\mathcal{A}$ = [$a_{ij}$] ∈ ${\mathbf{R}}^{N\times N}$. In $\mathcal{A}$, $a_{ij}$ = 1 if (${\mathit{v}}_{\mathit{i}}$,${\mathit{v}}_{\mathit{j}}$)$\in \mathcal{E}$, and $a_{ij}$ = 0 otherwise. This paper assumes that there are no self-loops; thus, $a_{ii}$ = 0. Correspondingly, the degree matrix is denoted by $\mathrm{\Delta }$ = diag{${\mathrm{\Delta }}_1,\cdots ,{\mathrm{\Delta }}_N$} with ${\mathrm{\Delta }}_i$ = ${\sum }_{j=1}^N$$a_{ij}$. Then, the Laplacian matrix associated with $\mathcal{G}$ can be defined as $\mathcal{L}$ = $\mathrm{\Delta }$−$\mathcal{A}$, where $\mathcal{L}$ = [$l_{ij}$] ∈ ${\mathbf{R}}^{N\times N}$ is defined as $l_{ii}$ = ${\sum }_{j=1,i\ne j}^N$$a_{ij}$, and $l_{ij}$ = $-a_{ij}$, $i\ne j$.

In this paper, we use ${\mu }_i$ = 1 to denote that the leader’s trajectory information is available for the i-th follower; otherwise, ${\mu }_i$ is set as 0. Let $\mathcal{P}=\mathrm{diag}\{{\mu }_1,\cdots ,{\mu }_N$}, and define W = $\mathcal{L}$ + $\mathcal{P}$.

Standard notations are presented as follows.

Q > 0	the matrix Q is positive definite
$Q^T$	the transposition of matrix Q
$\Vert Q\Vert$	the 2-norm of a matrix Q
${\lambda }_{\max }(·)$	the largest eigenvalues of the matrix
${\lambda }_{\min }(·)$	the smallest eigenvalues of the matrix
$1_N$	an N-dimensional column vector with all entries being 1
$D^{+}$	the upper-right-hand derivative

Lemma 1.

If a directed graph contains a spanning tree wherein the root node is the leader, then the matrix W = $\mathcal{L}$ + $\mathcal{P}$ is positive definite.

Lemma 2 (Schur Complement).

For a symmetric matrix Q = $\left[\begin{array}{cc}{Q}_{11}& Q_{12}\\ Q_{12}^T& Q_{22}\end{array}\right]$, the following statements are equivalent:

(1) 

$Q>0$,

(2) 

$Q_{11}>$0, $Q_{22}-Q_{12}^T{Q}_{11}^{-1}{Q}_{12}>$0,

(3) 

$Q_{22}>$0, $Q_{11}-Q_{12}{Q}_{22}^{-1}{Q}_{12}^T>$0.

2.2. Problem Statement

Consider an uncertain MAS with N followers and one leader. The i-th follower’s dynamics, $i=1,\cdots ,N$, are described as follows:

$$\left\{\begin{array}{cc}& {\dot{x}}_i(t)=v_i(t),\hfill \\ & {\dot{v}}_i(t)=b_i{u}_i(t)+{\theta }_i{f}_i(x_i,v_i,t)+d_i(t),i=1,\cdots ,N,\hfill \end{array}\right.$$

(1)

where $x_i(t)$∈ $\mathbf{R}$ and $v_i(t)$∈ $\mathbf{R}$ represent the i-th follower’s states; $u_i$∈ $\mathbf{R}$ is the i-th follower’s control input; $b_i$∈ $\mathbf{R}$ is a non-zero constant that is unknown, called the i-th follower’s control coefficient; ${\theta }_i$∈ $\mathbf{R}$ is an unknown constant; $f_i(x_i,v_i,t):\mathbf{R}\times \mathbf{R}\times {\mathbf{R}}_{+}\to \mathbf{R}$, is a known nonlinear function and is locally Lipschitz in ($x_i,v_i$) on $\mathbf{R}\times \mathbf{R}$ and continuous in t on ${\mathbf{R}}_{+}$; $d_i(t)$∈ $\mathbf{R}$ represents an external disturbance that is unknown and time-varying.

The leader’s dynamics are described by:

$$\left\{\begin{array}{cc}& {\dot{x}}_0(t)=v_0(t),\hfill \\ & {\dot{v}}_0(t)=u_0(t),\hfill \end{array}\right.$$

(2)

where $x_0(t)$∈ $\mathbf{R}$ and $v_0(t)$∈ $\mathbf{R}$ represent the states of the leader; $u_0(t)$ is the leader’s control input.

This paper aims to design a distributed adaptive event-triggered control protocol for each follower to achieve the following objectives:

(1)

All signals that are caused by closed-loop systems are bounded globally uniformly on [0, $+\infty$);

(2)

There is no Zeno behavior;

(3)

The tracking errors, ${\overline{x}}_i(t)$ = $x_i(t)$−$x_0(t)$, ${\overline{v}}_i(t)$ = $v_i(t)$−$v_0(t)$, converge to a small adjustable bounded set.

In order to realize the control objectives, the following assumptions are imposed on systems (1) and (2)

Assumption 1.

A directed graph $\mathcal{G}$ contains a spanning tree whose leader is the root node, and the communication topology between all followers is undirected.

Assumption 2.

The sign of $b_i$ is known.

Assumption 3.

External disturbance $d_i$(t) is bounded: that is, $|d_i(t)|\le \overline{d^{*}}$, where $\overline{d^{*}}$ is an unknown positive constant.

Assumption 4.

The leader’s control input $u_0(t)$ is bounded: i.e., $|u_0(t)|\le \rho$, where ρ is an unknown positive constant.

3. Event-Based Distributed Adaptive Control Protocol Design

3.1. Triggering Mechanism

Event times for the j-th agent are denoted as: $t_0^j,t_1^j,\cdots ,t_k^j,\cdots$ with 0 $\le t_0^j<t_1^j<t_2^j<\cdots <t_k^j<t_{k+1}^j<+\infty$ for k∈ ${\mathbf{Z}}^{+}$, j = 0, 1,⋯, N. The variable $t_k^j$ denotes the k-th event time for the j-th agent. In this paper, ${\check{x}}_j$ and ${\check{v}}_j$ are adopted to denote the sampling states for the j-th agent, and ${\check{x}}_j(t)$ = $x_j(t_k^j)$, ${\check{v}}_j(t)$ = $v_j(t_k^j)$, $t\in [t_k^j,t_{k+1}^j)$, which can be broadcast to neighbors at $t=t_k^j$. Meanwhile, the j-th agent updates the controller at $t=t_k^j$.

The event times are defined as

$$t_{k+1}^j=\inf \left\{t>t_k^j|p_j(t)\ge 0\right\},j=0,1,\cdots ,N,$$

where $t_1^j=0$, and the triggering mechanism is designed as

$$p_j(t)=|x_j(t)-{\check{x}}_j(t)|+|v_j(t)-{\check{v}}_j(t)|-q_j,$$

(3)

where $q_j$ is a positive constant to be designed. The triggering mechanism (3) means that: for agent j, once an event is triggered (i.e., $p_j(t_k^j)=0$), information transmission and control updating occur; after $t_k^j$, the triggering mechanism (3) is continuously evaluated until $p_j(t_{k+1}^j)=0$, which means that the ($k+1$)-th event is triggered. During any two consecutive events, if the neighbors of the j-th agent are triggered, the triggering mechanism (3) will not be updated. Since the states of the neighboring agents are not used in the triggering mechanism (3) for the j-th agent, the event times of the j-th agent are determined by themselves. It is noted that the continuous monitoring of the neighbors’ states is avoided since the triggering mechanism (3) only needs the current states of the j-th agent and the sampling states of the j-th agent.

3.2. Adaptive Control Protocol

In this section, we devise the adaptive control protocol for each follower based on the above triggering mechanism, and we analyze the existence and uniqueness of solutions that are caused by closed-loop systems.

Firstly, according to ${\overline{x}}_i(t)$ = $x_i(t)$−$x_0(t)$, ${\overline{v}}_i(t)$ = $v_i(t)$−$v_0(t)$, we deduce that

$$\left\{\begin{array}{cc}& {\dot{\overline{x}}}_i(t)={\overline{v}}_i(t),\hfill \\ & {\dot{\overline{v}}}_i(t)=b_i{u}_i(t)+{\theta }_i{f}_i(x_i,v_i,t)+d_i(t)-u_0(t),i=1,\cdots ,N.\hfill \end{array}\right.$$

(4)

For System (4), the following adaptive control protocol is devised for the i-th follower

$$\left\{\begin{array}{cc}& u_i={\hat{\beta }}_i{\alpha }_i,\hfill \\ & {\alpha }_i=-{\hat{\theta }}_i{f}_i-l_0({ϵ}_{i1}+{ϵ}_{i2})+{\hat{u}}_{i,0}-{\hat{d}}_i,\hfill \end{array}\right.$$

(5)

where ${\hat{\beta }}_i$ and ${\hat{\theta }}_i$ are the estimates of ${\beta }_i$ = $\frac{1}{b_i}$ and ${\theta }_i$, respectively; ${\hat{d}}_i$ is the estimate of $d_i$; ${\hat{u}}_{i,0}$ is the estimate of $u_0$ by the i-th follower; $l_0$ is a positive parameter to be designed; ${ϵ}_{i1}$ and ${ϵ}_{i2}$ are relative errors on $t\in [t_k,t_{k+1})$. The triggering mechanism is integrated into Control Protocol (5) by relative errors ${ϵ}_{i1}$ and ${ϵ}_{i2}$. In (5), an estimate-based method is used to estimate unknown parameters.

In (5), the estimates’ dynamics for the i-th follower are designed as

$$\left\{\begin{array}{cc}& {\dot{\hat{\theta }}}_i=r_{i1}(e_{i1}+e_{i2})f_i-r_{i1}{k}_{i1}({\hat{\theta }}_i-{\theta }_{i,0}),\hfill \\ & {\dot{\hat{\beta }}}_i=-r_{i2}(e_{i1}+e_{i2})\mathrm{sign}(b_i){\alpha }_i-r_{i2}{k}_{i2}({\hat{\beta }}_i-{\beta }_{i,0}),\hfill \\ & {\dot{\hat{u}}}_{i,0}=-r_{i3}(e_{i1}+e_{i2})-r_{i3}{k}_{i3}({\hat{u}}_{i,0}-u_{0,0}),\hfill \\ & {\dot{\hat{d}}}_i=r_{i4}(e_{i1}+e_{i2})-r_{i4}{k}_{i4}({\hat{d}}_i-d_{i,0}),\hfill \end{array}\right.$$

(6)

where $\mathrm{sign}(·)$ is the sign function; $r_{il}$ and $k_{il}$, l = 1, 2, 3, 4, are positive constants; ${\theta }_{i,0}$, ${\beta }_{i,0}$, $u_{0,0}$, $d_{i,0}$ are constants; $e_{i1}$ and $e_{i2}$ are relative errors.

In (6) and (7), the involved relative errors $e_{i1}$, $e_{i2}$, ${ϵ}_{i1}$, and ${ϵ}_{i2}$ are defined as

$$\left\{\begin{array}{cc}& e_{i1}(t)={\displaystyle \sum _{j=1}^N}{a}_{ij}(x_i(t)-x_j(t))+{\mu }_i(x_i(t)-x_0(t)),\hfill \\ & e_{i2}(t)={\displaystyle \sum _{j=1}^N}{a}_{ij}(v_i(t)-v_j(t))+{\mu }_i(v_i(t)-v_0(t)),\hfill \\ & {ϵ}_{i1}(t)={\displaystyle \sum _{j=1}^N}{a}_{ij}(x_i(t)-{\check{x}}_j(t))+{\mu }_i(x_i(t)-{\check{x}}_0(t)),\hfill \\ & {ϵ}_{i2}(t)={\displaystyle \sum _{j=1}^N}{a}_{ij}(v_i(t)-{\check{v}}_j(t))+{\mu }_i(v_i(t)-{\check{v}}_0(t)).\hfill \end{array}\right.$$

(7)

Next, we analyze the existence and uniqueness of solutions that are caused by closed-loop systems, consisting of (3)–(6). Note that the right-hand sides of differential Equations (4)–(6) are locally Lipschitz in (${\overline{x}}_i,{\overline{v}}_i$, ${\hat{\theta }}_i,{\hat{\beta }}_i,{\hat{u}}_{i,0},{\hat{d}}_i$) on $\mathbf{R}\times \mathbf{R}\times \mathbf{R}\times \mathbf{R}\times \mathbf{R}\times \mathbf{R}$ and continuous in t on ${\mathbf{R}}_{+}$. Therefore, by Peano’s existence theorem and the extension theorem [28], a closed-loop system has a unique solution on its maximal existence interval [0, $t_f$) with $0<$$t_f$≤$+\infty$ for any initial conditions ($x_i(0),v_i(0)$, ${\hat{\theta }}_i(0),{\hat{\beta }}_i(0),{\hat{u}}_{i,0}(0),{\hat{d}}_i(0)$, $x_0$(0), $v_0$(0)) on $\mathbf{R}\times \mathbf{R}\times \mathbf{R}\times \mathbf{R}\times \mathbf{R}\times \mathbf{R}\times \mathbf{R}\times \mathbf{R}$.

4. Main Results

In this section, we firstly construct a Lyapunov function. By Lemma 2, we prove the constructed Lyapunov function is positive definite. Then, we prove that all closed-loop system signals are bounded and exclude Zeno behavior on $[0,t_f)$ with $t_f<+\infty$. Furthermore, for all signals that cause the closed-loop system to be bounded, we deduce that $t_f=+\infty$. Finally, we prove that the devised adaptive event-triggered control protocol can realize bounded consensus tracking.

Theorem 1.

Consider Systems $(1)$ and $(2)$ under Assumptions 1–4. The designed distributed adaptive control protocol $(5)$ with event-triggering mechanism $(3)$ can guarantee: $(i)$ all closed-loop system signals $(i.e.,x_i,v_i,u_i,{\hat{\theta }}_i,{\hat{\beta }}_i,{\hat{u}}_{i,0},{\hat{d}}_i)$ are bounded on $[0,+\infty )$; $(ii)$ Zeno behavior is excluded; $(iii)$ for each follower, the tracking errors ${\overline{x}}_i(t)$ and ${\overline{v}}_i(t)$ converge to an adjustable bounded set.

Proof. 

Firstly, we prove that all signals that are caused by closed-loop systems are bounded on [0, $t_f$) with $t_f<+\infty$. According to the estimation errors for unknown parameters and the unknown leader’s control input, a Lyapunov function is constructed as

$$V=\frac{1}{2}{\delta }^T\mathrm{\Theta }\delta +\frac{1}{2}\sum _{i=1}^N\frac{1}{r_{i1}}{{\tilde{\theta }}_i}^2+\frac{1}{2}\sum _{i=1}^N\frac{|b_i|}{r_{i2}}{\tilde{\beta }}_i^2+\frac{1}{2}\sum _{i=1}^N\frac{1}{r_{i3}}{\tilde{u}}_{i,0}^2+\frac{1}{2}\sum _{i=1}^N\frac{1}{r_{i4}}{\tilde{d}}_i^2,$$

(8)

where $\mathrm{\Theta }$ = $\left(\begin{array}{cc}2W& W\\ W& W\end{array}\right)$, $\delta (t)$=$[{\delta }_1(t),\cdots ,{\delta }_N(t)]$ with ${\delta }_i$ = ${[{\overline{x}}_i,{\overline{v}}_i]}^T$, ${\tilde{\theta }}_i$ = ${\theta }_i$−${\hat{\theta }}_i$, ${\tilde{\beta }}_i$ = ${\beta }_i$−${\hat{\beta }}_i$, ${\tilde{u}}_{i,0}$ = $u_0$−${\hat{u}}_{i,0}$, ${\tilde{d}}_i$ = $d_i$−${\hat{d}}_i$. Under Assumption 1, W is symmetric, which means that the matrix $\mathrm{\Theta }$ is symmetric, and we obtain that W is positive definite by Lemma 1. In Lemma 2, $\mathrm{\Theta }$ can be regarded as Q; $2W$ can be regarded as $Q_{11}$; W can be regarded as $Q_{12}$ and $Q_{22}$. Then, we can deduce that $\mathrm{\Theta }$ > 0 by noting $2W$ > 0, $W^T{(2W)}^{-1}W>0$. Thus, we obtain that V is positive definite.

With this in hand, the derivative of V regarding time can be computed as

$$\dot{V}={\delta }^T\mathrm{\Theta }\dot{\delta }-\sum _{i=1}^N\frac{1}{r_{i1}}{\tilde{\theta }}_i{\dot{\hat{\theta }}}_i-\sum _{i=1}^N\frac{|b_i|}{r_{i2}}{\tilde{\beta }}_i{\dot{\hat{\beta }}}_i-\sum _{i=1}^N\frac{1}{r_{i3}}{\tilde{u}}_{i,0}{\dot{\hat{u}}}_{i,0}-\sum _{i=1}^N\frac{1}{r_{i4}}{\tilde{d}}_i{\dot{\hat{d}}}_i.$$

(9)

From (4), it follows that

$$\dot{\delta }(t)=\left(\begin{array}{cc}0& I_N\\ 0& 0\end{array}\right)\delta (t)+\left(\begin{array}{c}0\\ Bu(t)\end{array}\right)+\left(\begin{array}{c}0\\ \theta F(x,v,t)\end{array}\right)+\left(\begin{array}{c}0\\ D(t)\end{array}\right)-\left(\begin{array}{c}0\\ 1_N{u}_0(t)\end{array}\right),$$

(10)

where $B=\mathrm{diag}\{b_1,\cdots ,b_N\},u(t)={[u_1,\cdots ,u_N]}^T,\theta =\mathrm{diag}\{{\theta }_1,\cdots ,{\theta }_N\}$, $F(x,v,t)={[f_1(x_1,v_1,t),\cdots ,f_N(x_N,v_N,t)]}^T$, $D(t)={[d_1(t),\cdots ,d_N(t)]}^T.$

Based on (5) and (10), (9) can be further computed as

$$\begin{array}{cc}\hfill \dot{V}=& {\delta }^T\left(\begin{array}{cc}0& 2W\\ 0& W\end{array}\right)\delta +\sum _{i=1}^N\frac{1}{r_{i1}}{\tilde{\theta }}_i\left(r_{i1}(e_{i1}+e_{i2})f_i-{\dot{\hat{\theta }}}_i\right)+\sum _{i=1}^N\frac{1}{r_{i3}}{\tilde{u}}_{i,0}(-r_{i3}(e_{i1}+e_{i2})\hfill \\ \hfill & -{\dot{\hat{u}}}_{i,0})+\sum _{i=1}^N\frac{1}{r_{i4}}{\tilde{d}}_i\left(r_{i4}(e_{i1}+e_{i2})-{\dot{\hat{d}}}_i\right)-\sum _{i=1}^N{l}_0(e_{i1}+e_{i2})({ϵ}_{i1}+{ϵ}_{i2})\hfill \\ \hfill & +\sum _{i=1}^N\frac{|b_i|}{r_{i2}}{\tilde{\beta }}_i\left(-r_{i2}(e_{i1}+e_{i2})\mathrm{sign}(b_i){\alpha }_i-{\dot{\hat{\beta }}}_i\right).\hfill \end{array}$$

Noting that ${\delta }^T\left(\begin{array}{cc}0& W\\ W& W\end{array}\right)\delta$ = ${\delta }^T\left(\begin{array}{cc}0& 2W\\ 0& W\end{array}\right)\delta$, we have

$$\begin{array}{cc}\hfill \dot{V}=& {\delta }^T\left(\begin{array}{cc}0& W\\ W& W\end{array}\right)\delta +\sum _{i=1}^N\frac{1}{r_{i1}}{\tilde{\theta }}_i\left(r_{i1}(e_{i1}+e_{i2})f_i-{\dot{\hat{\theta }}}_i\right)+\sum _{i=1}^N\frac{1}{r_{i3}}{\tilde{u}}_{i,0}(-r_{i3}(e_{i1}+e_{i2})\hfill \\ \hfill & -{\dot{\hat{u}}}_{i,0})+\sum _{i=1}^N\frac{1}{r_{i4}}{\tilde{d}}_i\left(r_{i4}(e_{i1}+e_{i2})-{\dot{\hat{d}}}_i\right)-\sum _{i=1}^N{l}_0(e_{i1}+e_{i2})({ϵ}_{i1}+{ϵ}_{i2})\hfill \\ \hfill & +\sum _{i=1}^N\frac{|b_i|}{r_{i2}}{\tilde{\beta }}_i\left(-r_{i2}(e_{i1}+e_{i2})\mathrm{sign}(b_i){\alpha }_i-{\dot{\hat{\beta }}}_i\right).\hfill \end{array}$$

(11)

Substituting (6) into (11) yields

$$\begin{array}{cc}\hfill \dot{V}=& {\delta }^T\left(\begin{array}{cc}0& W\\ W& W\end{array}\right)\delta +\sum _{i=1}^N{k}_{i1}{\tilde{\theta }}_i({\hat{\theta }}_i-{\theta }_{i,0})+\sum _{i=1}^N{k}_{i2}\left|b_i\right|{\tilde{\beta }}_i\left({\hat{\beta }}_i-{\beta }_{i,0}\right)+\sum _{i=1}^N{k}_{i4}{\tilde{d}}_i({\hat{d}}_i-d_{i,0})\hfill \\ \hfill & +\sum _{i=1}^N{k}_{i3}{\tilde{u}}_{i,0}({\hat{u}}_{i,0}-u_{0,0})-\sum _{i=1}^N{l}_0(e_{i1}+e_{i2})({ϵ}_{i1}+{ϵ}_{i2}).\hfill \end{array}$$

(12)

By [24], the following inequalities are acquired:

$$\left\{\begin{array}{cc}& {\tilde{\theta }}_i({\hat{\theta }}_i-{\theta }_{i,0})\le -\frac{1}{2}{\tilde{\theta }}_i^2+\frac{1}{2}{({\theta }_i-{\theta }_{i,0})}^2,\hfill \\ & {\tilde{\beta }}_i({\hat{\beta }}_i-{\beta }_{i,0})\le -\frac{1}{2}{\tilde{\beta }}_i^2+\frac{1}{2}{({\beta }_i-{\beta }_{i,0})}^2,\hfill \\ & {\tilde{d}}_i({\hat{d}}_i-d_{i,0})\le -\frac{1}{2}{\tilde{d}}_i^2+\frac{1}{2}{(d_i-d_{i,0})}^2,\hfill \\ & {\tilde{u}}_{i,0}({\hat{u}}_{i,0}-u_{0,0})\le -\frac{1}{2}{\tilde{u}}_{i,0}^2+\frac{1}{2}{(u_{i,0}-u_{0,0})}^2.\hfill \end{array}\right.$$

(13)

Substituting (13) into (12), it can be deduced that

$$\begin{array}{cc}\hfill \dot{V}\le & {\delta }^T\left(\begin{array}{cc}0& W\\ W& W\end{array}\right)\delta -\sum _{i=1}^N\frac{k_{i1}}{2}{\tilde{\theta }}_i^2-\sum _{i=1}^N\frac{k_{i2}}{2}\left|b_i\right|{\tilde{\beta }}_i^2-\sum _{i=1}^N\frac{k_{i3}}{2}{\tilde{u}}_{i,0}^2-\sum _{i=1}^N\frac{k_{i4}}{2}{\tilde{d}}_i^2\hfill \\ \hfill & +\sum _{i=1}^N\frac{k_{i1}}{2}{({\theta }_i-{\theta }_{i,0})}^2+\sum _{i=1}^N\frac{k_{i2}}{2}\left|b_i\right|{({\beta }_i-{\beta }_{i,0})}^2+\sum _{i=1}^N\frac{k_{i3}}{2}{(u_0-u_{0,0})}^2\hfill \\ \hfill & +\sum _{i=1}^N\frac{k_{i4}}{2}{(d_i-d_{i,0})}^2-\sum _{i=1}^N{l}_0(e_{i1}+e_{i2})({ϵ}_{i1}+{ϵ}_{i2}).\hfill \end{array}$$

(14)

where ${ϵ}_{i1}+{ϵ}_{i2}=e_{i1}+e_{i2}+{\displaystyle \sum _{j=1}^N}{a}_{ij}(x_j\left(t\right)-{\check{x}}_j\left(t\right))+{\mu }_i(x_0\left(t\right)-{\check{x}}_0\left(t\right))+{\displaystyle \sum _{j=1}^N}{a}_{ij}(v_j\left(t\right)-{\check{v}}_j\left(t\right))+{\mu }_i(v_0\left(t\right)-{\check{v}}_0\left(t\right))$.

Moreover, by Young’s inequality, we have

$$\begin{array}{c}\hfill l_{0}q|e_{i1}+e_{i2}|·|{\Delta }_i+{\mu }_i|\le \frac{l_0}{2}{(e_{i1}+e_{i2})}^2+\frac{l_0}{2}{q}^2{({\Delta }_i+{\mu }_i)}^2,\end{array}$$

where q = max{$q_j$, $j=0,1,\cdots ,N$}.

By (7), note ${\displaystyle \sum _{i=1}^N}{(e_{i1}+e_{i2})}^2={\delta }^T{\left(\begin{array}{cc}W& 0\\ 0& W\end{array}\right)}^T\left(\begin{array}{cc}W& 0\\ 0& W\end{array}\right)\delta$. Then, from (14), it follows that

$$\begin{array}{cc}\hfill \dot{V}\le & -{\delta }^T\left(\frac{l_0}{2}{\left(\begin{array}{cc}W& 0\\ 0& W\end{array}\right)}^T\left(\begin{array}{cc}W& 0\\ 0& W\end{array}\right)-\left(\begin{array}{cc}0& W\\ W& W\end{array}\right)\right)\delta -\sum _{i=1}^N(\frac{k_{i1}}{2}{\tilde{\theta }}_i^2+\frac{k_{i2}}{2}\left|b_i\right|{\tilde{\beta }}_i^2\hfill \\ \hfill & +\frac{k_{i3}}{2}{\tilde{u}}_{i,0}^2+\frac{k_{i4}}{2}{\tilde{d}}_i^2)+\sum _{i=1}^N(\frac{l_0}{2}{q}^2{({\mathrm{\Delta }}_i+{\mu }_i)}^2+\frac{k_{i1}}{2}{({\theta }_i-{\theta }_{i,0})}^2\hfill \\ \hfill & +\frac{k_{i2}}{2}\left|b_i\right|{({\beta }_i-{\beta }_{i,0})}^2+\frac{k_{i3}}{2}{(u_0-u_{0,0})}^2+\frac{k_{i4}}{2}{(d_i-d_{i,0})}^2).\hfill \end{array}$$

(15)

Let $\Gamma$ = ${\left(\begin{array}{cc}W& 0\\ 0& W\end{array}\right)}^T$$\left(\begin{array}{cc}W& 0\\ 0& W\end{array}\right)$, $\prod$ = $\left(\begin{array}{cc}0& W\\ W& W\end{array}\right)$. Based on (15), we obtain

$$\begin{array}{cc}\hfill \dot{V}\le & -\left(\frac{l_0}{2}{\lambda }_{\min }(\Gamma )-{\lambda }_{\max }(\prod )\right){\delta }^T\delta -\sum _{i=1}^N\omega \left(\frac{{\tilde{\theta }}_i^2}{2r_{i1}}+\frac{|b_i|{\tilde{\beta }}_i^2}{2r_{i2}}+\frac{{\tilde{u}}_{i,0}^2}{2r_{i3}}+\frac{{\tilde{d}}_i^2}{2r_{i4}}\right)\hfill \\ \hfill & +\sum _{i=1}^N(\frac{l_0}{2}{m}^2{({\mathrm{\Delta }}_i+{\mu }_i)}^2+\frac{k_{i1}}{2}{({\theta }_i-{\theta }_{i,0})}^2+\frac{k_{i2}}{2}\left|b_i\right|{({\beta }_i-{\beta }_{i,0})}^2+\frac{k_{i3}}{2}{(u_0-u_{0,0})}^2\hfill \\ \hfill & +\frac{k_{i4}}{2}{(d_i-d_{i,0})}^2),\hfill \end{array}$$

(16)

where $\omega$ = min$\{k_{i1}{r}_{i1},k_{i2}{r}_{i2},k_{i3}{r}_{i3},k_{i4}{r}_{i4}\}$.

Furthermore, from (8), it follows that

$$\begin{array}{c}\hfill V\le \frac{1}{2}{\lambda }_{\max }(\mathrm{\Theta }){\delta }^T\delta +\sum _{i=1}^N\left(\frac{{\tilde{\theta }}_i^2}{2r_{i1}}+\frac{|b_i|{\tilde{\beta }}_i^2}{2r_{i2}}+\frac{{\tilde{u}}_{i,0}^2}{2r_{i3}}+\frac{{\tilde{d}}_i^2}{2r_{i4}}\right).\end{array}$$

(17)

Combining (16) with (17), we obtain

$$\begin{array}{c}\hfill \dot{V}\le -(\frac{l_0}{2}{\lambda }_{\min }(\Gamma )-{\lambda }_{\max }(\prod )-\frac{\omega }{2}{\lambda }_{\max }(\mathrm{\Theta })){\delta }^T\delta -\omega V+M^{*},\end{array}$$

where $M^{*}$ = ${\displaystyle \sum _{i=1}^N}$($\frac{l_0}{2}{q}^2$${({\mathrm{\Delta }}_i+{\mu }_i)}^2$ + $\frac{k_{i1}}{2}$${({\theta }_i-{\theta }_{i,0})}^2$ + $\frac{k_{i2}}{2}\left|b_i\right|$${({\beta }_i-{\beta }_{i,0})}^2$ + $\frac{k_{i3}}{2}$${(u_0-u_{0,0})}^2$ + $\frac{k_{i4}}{2}$${(d_i-d_{i,0})}^2$). Then, by choosing $l_0$>${\lambda }_{\min }^{-1}$($\Gamma$)($\omega {\lambda }_{\max }(\mathrm{\Theta })+2{\lambda }_{\max }(\prod )$), we deduce that

$$\begin{array}{c}\hfill \dot{V}\le -\omega V+M^{*}.\end{array}$$

(18)

By solving (18), we derive

$$\begin{array}{c}\hfill V(t)\le V(0)e^{-\omega t}+\frac{M^{*}}{\omega }(1-e^{-\omega t}),\end{array}$$

(19)

which shows that $V(t)$ is uniformly bounded on [0, $t_f$). Thus, the signals ${\delta }_i$, ${\hat{\theta }}_i$, ${\hat{\beta }}_i$, ${\hat{u}}_{i,0}$, and ${\hat{d}}_i$ are bounded on [0, $t_f$). Under Assumption 4, it is concluded that $x_0$ and $v_0$ are bounded on [0, $t_f$). Therefore, $x_i$ and $v_i$ are bounded by the boundedness of ${\delta }_i$, $x_0$, $v_0$ on [0, $t_f$). Furthermore, $u_i$ is bounded on [0, $t_f$) from (5).

Now, we prove that there is no Zeno behavior on [0, $t_f$): that is, the lower bound on the interval between two consecutive events {$t_{k+1}^j-$$t_k^j$} is positive.

Define ${\phi }_k^j(t)$ = $|x_j(t)$−${\check{x}}_j(t)|$ + $|v_j(t)$−${\check{v}}_j(t)|$ for t∈ $[t_k^j,t_{k+1}^j)$. First, the derivative of the triggering mechanism ${\phi }_k^j(t)$ regarding time needs to be acquired. Then, the maximum velocity $h_j$ of the change in ${\phi }_k^j(t)$ from 0 to $q_j$ can be acquired. We exclude Zeno behavior by testing if $\frac{q_j}{h_j}$ is positive. The derivative of ${\phi }_k^j(t)$ regarding time is computed as

$$\begin{array}{cc}\hfill {\dot{\phi }}_k^j(t)& =D^{+}|x_j(t)|+D^{+}\left|v_j(t)\right|\hfill \\ \hfill & =|v_j(t)|+|b_j{u}_j(t)+f_j(x_j,v_j,t){\theta }_j+d_j(t)|,j=0,1,\cdots ,N\hfill \end{array}$$

and ${\dot{\phi }}_k^0(t)$≤∣${\dot{x}}_0$(t)∣+∣${\dot{v}}_0$(t)∣. By the boundedness of $u_j$, $v_j$, ${\dot{x}}_0$, and ${\dot{v}}_0$ on [0, $t_f$), we acquire a positive number $h_j$ such that ${\dot{\phi }}_k^j(t)$≤$h_j$. Moreover, we get that ${\int }_{t_{k+1}^j}^{t_k^j}{\dot{\phi }}_k^j(t)\mathrm{d}t\le {\int }_{t_{k+1}^j}^{t_k^j}{h}_j\mathrm{d}t=h_j(t_{k+1}^j-t_k^j)$. From (3), it follows that ${\phi }_{k+1}^j(t_{k+1}^j)-{\phi }_k^j(t_k^j)>q_j$. Thus, we arrive at $t_{k+1}^j$−$t_k^j$≥$\frac{q_j}{h_j}$. Zeno behavior is excluded on [0, $t_f$). Furthermore, all signals that cause closed-loop systems are bounded; we deduce that $t_f$ = +∞.

Finally, that the tracking errors ${\overline{x}}_i(t)$ and ${\overline{v}}_i(t)$ converge to a small adjustable bounded set was testified. From (9), it follows that

$$V(t)\ge \frac{1}{2}{\delta }^T(t)\left(\begin{array}{cc}2W& W\\ W& W\end{array}\right)\delta (t)\ge \frac{1}{2}{\lambda }_{\min }(\mathrm{\Theta }){\Vert \delta (t)\Vert }^2,$$

then,

$${\Vert \delta (t)\Vert }^2\le \frac{2}{{\lambda }_{\min }(\mathrm{\Theta })}V(t).$$

(20)

Substituting (19) into (20), we obtain

$${\Vert \delta (t)\Vert }^2\le \frac{2}{{\lambda }_{\min }(\mathrm{\Theta })}\left[V(0)e^{-\omega t}+\frac{M^{*}}{\omega }(1-e^{-\omega t})\right],$$

which implies that ${\Vert \delta (t)\Vert }^2$ is bounded and converges to a compact set $E_r$ = {$\delta (t)$ | ${\Vert \delta (t)\Vert }^2$ ≤ $\frac{2}{\omega {\lambda }_{\min }(\mathrm{\Theta })}$ ($M^{*}$+$\sigma$)} for t ≥ (1/$\omega$) ln(∣$\omega$$V(0)$−$M^{*}$∣/$\sigma$) with $\sigma$ an arbitrarily small positive constant. Note ${\Vert \delta (t)\Vert }^2={\overline{x}}_1^2+\cdots +{\overline{x}}_N^2+{\overline{v}}_1^2+\cdots +{\overline{v}}_N^2$. Therefore, tracking errors ${\overline{x}}_i$(t) and ${\overline{v}}_i$(t) also converge to $E_r$ for t≥ (1/$\omega$) ln(∣$\omega$$V(0)$−$M^{*}$∣/$\sigma$). It is noted that the compact set $E_r$ can be made as small as desired by increasing $r_{i1},r_{i2},r_{i3},r_{i4}$ while fixing all the remaining design parameters. Thus, the compact set $E_r$ can be adjusted by modifying the relevant parameters. □

5. Simulation of Nonlinear MASs with Uncertainties

In this section, we provide a simulation example in detail to verify the theoretical results obtained. Consider an uncertain MAS with four followers and one leader. The i-th follower’s dynamics, $i=1,\cdots ,4$, are described as (1), where $b_1$ = $b_2$ = $b_3$ = 1, $b_4$ = 1.1, ${\theta }_1$ = 1, ${\theta }_2$ = 2, ${\theta }_3$ = 3, ${\theta }_4$ = 4, $f_1(x_1,v_1,t)=0.1x_1^2+0.1v_1^2$, $f_2(x_2,v_2,t)=0.2x_2^2+0.1v_2^3$, $f_3(x_3,v_3,t)=0.2x_3^2+0.1v_3^2$, $f_4(x_4,v_4,t)=0.1x_4^3+0.1v_4^2$, $d_1(t)=\sin (t),d_2(t)=1.1\sin (t),d_3(t)=1.2\sin (t),d_4(t)=2\sin (t)$. The leader’s dynamics are described by (2), where $u_0(t)=-\sin (t)$.

Note that the communication topology for systems (21) and (22) is shown in Figure 1. The four followers’ initial conditions and the leader’s conditions are given as follows: $x_1(0)$ = 0, $v_1(0)$ = 1, $x_2(0)$ = 1, $v_2(0)$ = 2, $x_3(0)$ = 2, $v_3(0)$ = 3, $x_4(0)$ = 4, $v_4(0)$ = 2, $x_0(0)$ = 0, $v_0(0)$ = 1. The initial conditions ${\hat{\theta }}_i(0)$, ${\hat{\beta }}_i(0)$, ${\hat{u}}_{i,0}(0)$, ${\hat{d}}_i(0)$, $i=1,\cdots ,4$ are set to zero. The corresponding parameters in (6) and (7) are chosen as $l_0$ = 170, $k_{i1}$ = $k_{i2}$ = $k_{i3}$ = $k_{i4}$ = 5, $r_{i1}$ = $r_{i2}$ = $r_{i3}$ = $r_{i4}$ = 0.005, ${\theta }_{i,0}$ = ${\beta }_{i,0}$ = $u_{0,0}$ = $d_{i,0}$ = 0. The designed parameters in (3) are taken as $m_1=0.95,m_2=1,m_3=2,m_4=2.65,m_0=1.1$. From Figure 2, we can see that the evolution of tracking errors, ${\overline{x}}_i=x_i-x_0$, ${\overline{v}}_i=v_i-v_0$, converges to a bounded set. Figure 3 illustrates that the evolution of the control protocols for the four followers is bounded. The four followers’ triggering times and the leader’s triggering times are exhibited in Figure 4.

In light of the work [27], a distributed adaptive control protocol without ETC is devised to achieve consensus tracking for systems (1) and (2). The trajectories of tracking errors, ${\overline{x}}_i$ and ${\overline{v}}_i$, are shown in Figure 5.

The simulation experiments are run to make

Table 2. Comparison between the consensus control protocol with triggering mechanism (CWT) and the consensus control protocol without triggering mechanism (CWNT).

Protocol	ACT	IAU
CWT	5 s	205.4932
CWNT	8.5 s	622.5718

, which shows the average consensus time (ACT) and the integral of the absolute value of the input (IAU). ACT means the average time it takes for followers to track the leader’s trajectory. IAU can evaluate the energy consumption of control signals [29], where larger values of IAU mean greater energy consumption for the control signals.

From Table 2, we acquire that the devised consensus protocol in this paper can make all followers track the leader faster. Further, when the triggering mechanism is considered in the design of the control protocol, energy resources can be better conserved.

6. Conclusions

This paper has addressed bounded consensus tracking for nonlinear MASs with uncertainties by ETC. The considered MASs permit multiple uncertainties, including an unknown control coefficient, parameterized unknown nonlinearities, unknown external disturbances, and unknown control input of the leader. To compensate for uncertainties, a new estimate-based adaptive control protocol has been proposed. Integrating a triggering mechanism into the adaptive control protocol, the bounded consensus tracking of MASs can be achieved with fewer communication resources. Then, by choosing an appropriate Lyapunov function candidate, it has been proved that the designed control protocol can guarantee that there is no Zeno behavior and that tracking errors ultimately converge to a bounded set. For future study, extending the results to finite-time consensus will be interesting work.