<figure><img src="image_1.png" /></figure>

HAL
open science

# Robustness Analysis for Multi-Agent Consensus Systems with Application to DC Motor Synchronization

Daniel Olivares, Gerardo Romero, Jose Guerrero, Rogelio Lozano

► To cite this version:

Daniel Olivares, Gerardo Romero, Jose Guerrero, Rogelio Lozano. Robustness Analysis for Multi-Agent Consensus Systems with Application to DC Motor Synchronization. Applied Sciences, 2020, 10 (18), pp.6521. 10.3390/app10186521. hal-02965278

HAL Id: hal-02965278

https://hal.science/hal-02965278

Submitted on 2 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Article

# Robustness Analysis for Multi-Agent Consensus Systems with Application to DC Motor Synchronization

Daniel Olivares <fnref n="1" />,†,✉, Gerardo Romero <fnref n="2" />,*,†,✉, Jose A. Guerrero <fnref n="3" />,†,✉ and Rogelio Lozano <fnref n="4" />,<fnref n="5" />,†,✉

Received: 1 August 2020; Accepted: 15 September 2020; Published: 18 September 2020

<figure><img src="image_4.png" /></figure>

check for
updates

**Abstract:** DC motor speed synchronization is a critical problem in industrial and robotic applications. To tackle this problem, we propose to use a multi-agent consensus-based control scheme that guarantees the convergence of the DC motor speeds to either fixed or time-varying reference. A detailed robustness analysis considering parametric uncertainty and time delay in the multi-agent system is performed to guarantee the consensus on the speed of DC motors in actual practice. The results obtained concerning the robustness analysis allowed us to implement experimental tests on a three-motor system using a wireless communication system to achieve satisfactory performance.

**Keywords:** multi-agent systems; consensus; synchronization; DC motors; time delay systems; robust control

## 1. Introduction

Multi-agent coordination control has been a research field of interest during the last few decades. A wide variety of multi-agent control methods can be found in the literature. These control approaches are classified as (a) virtual structure [1], where each agent is considered a component of a bigger structure; (b) behavioral control [2], which deals with conflicting tasks such as obstacle avoidance and move to goal, for example; and (c) leader-follower [3,4], where one of the agents works as a group leader and the others are designated as followers. The consensus is one of the most studied strategies to control multi-agent systems [4–7], and its applications include power systems [8], robotics [3,9], flight formation control [10], to mention a few.

Arguably, the most common actuator in robotics and industrial applications is the DC motor. It is well known that motor synchronization plays an important role in the quality of manufactured products. For example, in rolling mill machines, the conveyor belt's tension must lie in an allowable tension to guarantee the thickness of the rolled steel [11]. As discussed in [12], many manufacturing processes require coordinated motion of two or more motion axes to achieve the desired performance, e.g., CNC machines, paper machines. Poor motor synchronization in such industrial applications would result in low-quality products that might be considered unusable in the worst case. The problem

<footnote n="1">¹ Department of Electrical and Electronics Engineering, Instituto Tecnológico de Nuevo Laredo, Tecnológico Nacional de México, Nuevo Laredo, Tamaulipas C.P. 88000, Mexico; daniel.oc@nlaredo.tecnm.mx</footnote>

<footnote n="2">² Electronics Department at U.A.M. Reynosa-Rodhe, Universidad Autónoma de Tamaulipas, Reynosa, Tamaulipas C.P. 88779, Mexico</footnote>

<footnote n="3">³ Institut Pascal, Université Clermont Auvergne, CNRS, SIGMA Clermont, F-63000 Clermont-Ferrand, France; jguerrer@ieee.org</footnote>

<footnote n="4">⁴ Sorbonne Universites, UTC, CNRS, UMR, 7253 Heudiasyc, France; rlozano@hds.utc.fr</footnote>

<footnote n="5">⁵ UMI LAFMIA 3175 CNRS-CINVESTAV, Ciudad de México C.P. 07360, Mexico</footnote>

<footnote>* Correspondence: gromero@docentes.uat.edu.mx</footnote>

<footnote>† These authors contributed equally to this work.</footnote>
of designing control algorithms for motor synchronization has been considered in different research works. In [13], the authors presented some results on the synchronization of two motors' speed using the leader-follower approach, where leader speed was used as a reference for the follower system, and the error was controlled by an LQR controller. A control scheme based on an adaptive sliding mode control technique was employed into the synchronization of the control system for four motors in [14]. In [15], a synchronization control strategy of multiple induction motors was proposed. The approach is based on an adjacent cross-coupling control structure, and then it can use a total sliding mode control method. A stator flux regulation method is used in [16] to apply an adaptive time-delay feedback control and the direct torque control to synchronize the chaotic speed of induction motors in a master-slave scheme. Three common synchronization methods for chaotic systems, namely an SMC control, a PI control, and an adaptive controller, were used in [17] to synchronize two brushless DC motor models in a master-slave configuration and were compared by using well-known performance indices.

Recently, the theory and applications of robust control in the presence of uncertainty and time delays have been addressed in the literature using a variety of control approaches [18–21], for example. The problem of the robustness of a multi-agent consensus control in the presence of uncertainty and time delays has been studied using various approaches [22–25]. In [22], a robust controller using algebraic Ricatti equations is presented. The agents are assumed to have linear dynamics, while uncertainty is treated as additive perturbations, which are assumed to be stable and $H_\infty$ bounded. Sliding mode control is proposed in [23] to deal with uncertainty as additive perturbations in a leader-follower multi-agent system. In a similar way than [22], it is assumed that additive perturbation is stable and $H_\infty$ bounded. In [24], a stochastic Lyapunov function to solve the consensus-seeking problem in networks where the agents can only obtain noisy measurements of its neighbor's state is proposed. In these previously mentioned papers, dynamic uncertainty is considered, which is introduced to the system as an additive disturbance, and, in other cases, it is additive noise in contrast to this paper where uncertainty is considered in each of the parameters. A parametric uncertainty-based approach is proposed in [25]; here, the multi-agent system is transformed into a controllable form where the state matrix's coefficients are assumed to be uncertain. Then, the maximal input's time delay for which the multi-agent system stability is guaranteed is determined using the value set method. Here, an assumption is made in the multi-agent system to achieve the particular structure that allows for using the value set tool in contrast with our paper where that assumption is not made.

In this paper, we present a multi-agent consensus-based control scheme to synchronize the speed of DC motors. Unlike the methods discussed above, it can control multiple motors. In addition, in the multi-agent scheme, the motors are synchronized first to reach the consensus and then follow the reference, unlike the parallel or master-slave scheme. In this work, the leader's control is based on the center of the state's mass and the motors' speed; thus, the response depends on all speeds, making the transient smoother. Previous to experimentation, an in-depth analysis was carried out in the system's robustness to guarantee the proper function of the multi-agent system; that is, the tolerance to phenomena not considered in the mathematical model. First, uncertainty limits that may be withstood by the closed-loop multi-agent system are calculated in the model parameters. Then, time delay is added to the system, and the maximum value it can withstand is computed; both calculations are performed without losing the stability property. The above provides the confidence to apply the control strategy for synchronizing the speed of the three DC motors, which is perfectly validated with the experimental results presented at the end of this work.

This paper's organization is as follows: the modeling of the information flow among agents using algebraic graph theory is presented in Section 2. The DC motor dynamic model and a control law based on multi-agent consensus control developed for DC motor speed synchronization are presented in Section 3. In Section 4, the robustness analysis of the multi-agent system in the presence of uncertainty and time delays is presented. In Section 5, experimental results are discussed. Finally, conclusions and future work are presented in Section 6.
## 2. Multi-Agent Consensus System Preliminaries

As discussed in [6], information exchange among agents is a necessary condition for multi-agent consensus. Such information exchange is often modelled using graph theory [6,26,27]. Then, in this section, some preliminary results on graph theory and forced multi-agent consensus are presented.

### 2.1. Graph Theory

A graph $\mathcal{G}$ is a pair $(\mathcal{N}, \mathcal{E})$ which consists of a set of nodes $\mathcal{N}$ and a set of edges $\mathcal{E}$ on $\mathcal{N}$. Every link between a pair of nodes $(n_i, n_j)$ is called an edge, $e \in \mathcal{E}$. A graph is said to be bidirectional if nodes $n_i$ and $n_j$ can get information from each other, i.e., $(n_i, n_j) \in \mathcal{E} \Leftrightarrow (n_j, n_i) \in \mathcal{E}$. A graph is called digraph if the information between nodes only flows from node $n_i$ to $n_j$. Thus, information exchange between agents can be modeled as a digraph but also as an undirected graph. A graph is said to be balanced if its in-degree (the number of communication links arriving at the node) is equal to its out-degree (the number of communication links leaving the node). A graph is connected if, for every pair $(n_i, n_j)$ of distinct vertices, there is a path from $n_i$ to $n_j$. A connected graph allows the communication between all agents through the network. This leads us to the fact that a graph $\mathcal{G}$ has a spanning tree if and only if $\mathcal{G}$ is connected [27]. We assume that every communication link $(n_i, n_j)$ has a gain (weight) $g_{ij}$. A matrix representation of the information flow $\mathcal{G}$ among agents can be obtained using Laplacian matrix. The definition of the Laplacian matrix is $\mathcal{L} = [l_{ij}]$ where $l_{ii} = \sum_{j \in \mathcal{N}_i} g_{ij}$ and $l_{ij} = -g_{ij}$ for $i \neq j$ and $\mathcal{N}_i$ is the set of neighbors of the $i^{th}$ agent.

### 2.2. Multi-Agent Forced Consensus Control

Let us consider the following multi-agent system:

$$ \dot{\xi}_i = \bar{u}_i \quad (1) $$

where $\xi$ is the state and $\bar{u}_i$ is a forced consensus control strategy given by

$$ \bar{u}_i = \sum_{j \in \mathcal{N}_i} -k_p (\xi_i - \xi_j) + b_i u_i \quad (2) $$

where $k_p$ is a positive gain and $b_i$ is defined as follows:

$$ b_i = \begin{cases} 0 & \text{if the } i^{th} \text{ agent is follower} \\ 1 & \text{if the } i^{th} \text{ agent is leader} \end{cases} $$

Then, system (1) can be rewritten:

$$ \dot{\xi} = -\mathcal{L}\xi + Bu \quad (3) $$

$$ y = C\xi \quad (4) $$

where $C = B^T$ are the output and input vectors and $\mathcal{L}$ is the Laplacian matrix of $\mathcal{G}$, the multi-agent information graph. Notice that, without loss of generality, agent 1 can be designated as the leader agent with access to the reference as shown in Figure 1. The leader control law considered here is based on the center of mass (average) of the states [4]

$$ u_l = N k_{cm} \left( \xi_{cm}^d - \xi_{cm} \right) \quad (5) $$

where $k_{cm}$ is a positive gain, $N$ is the number of agents in the multi-agent system, and $\xi_{cm}^d$ and $\xi_{cm}$ are the desired state center of mass and the calculated system state center of mass, respectively.
## 2.3. Information Flow Topology

The information flow between agents is defined by the employed topology. In [4], it is stated that cyclic and chain topology, shown in Figure 1, can be combined to model all kinds of complex multi-agent information flow topologies. It is important to recall that the forced consensus of multi-agent systems of arbitrary length communicating over networks with cyclic and chain topologies are controllable and observable. Details on controllability and observability properties of leader-based multi-agent systems can be found in [4].

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Information flow configuration: (a) cyclic topology; (b) chain topology. The gray circle represents the leader agent and white circles are followers.</figcaption></figure>

The state space representation of a multi-agent system with cyclic topology is the following:

$$ \begin{bmatrix} \dot{\xi}_1 \\ \dot{\xi}_2 \\ \vdots \\ \dot{\xi}_n \end{bmatrix} = \begin{bmatrix} -1 & 1 & 0 & 0 \\ -1 & 1 & & \\ & & \ddots & \\ 1 & & & -1 \end{bmatrix} \begin{bmatrix} \xi_1 \\ \xi_2 \\ \vdots \\ \xi_n \end{bmatrix} + \begin{bmatrix} 1 \\ 0 \\ \vdots \\ 0 \end{bmatrix} u \quad (6) $$

$$ y = \begin{bmatrix} 1 & 0 & \cdots & 0 \end{bmatrix} \begin{bmatrix} \xi_1 \\ \xi_2 \\ \vdots \\ \xi_n \end{bmatrix} $$

The state space representation of a multi-agent system with the chain topology is the following:

$$ \begin{bmatrix} \dot{\xi}_1 \\ \dot{\xi}_2 \\ \vdots \\ \dot{\xi}_n \end{bmatrix} = \begin{bmatrix} -1 & 1 & & & \\ 1 & -2 & 1 & & \\ & \ddots & \ddots & \ddots & \\ & & 1 & -2 & 1 \\ & & & 1 & -1 \end{bmatrix} \begin{bmatrix} \xi_1 \\ \xi_2 \\ \vdots \\ \xi_n \end{bmatrix} + \begin{bmatrix} 1 \\ 0 \\ \vdots \\ 0 \end{bmatrix} u \quad (7) $$

$$ y = \begin{bmatrix} 1 & 0 & \cdots & 0 \end{bmatrix} \begin{bmatrix} \xi_1 \\ \xi_2 \\ \vdots \\ \xi_n \end{bmatrix} $$
For simplicity, Equations (6) and (7) can be rewritten as

$$ 
\begin{bmatrix} \dot{\xi}_1 \\ \dot{\xi}_2 \\ \vdots \\ \dot{\xi}_n \end{bmatrix} = -\mathcal{L} \begin{bmatrix} \xi_1 \\ \xi_2 \\ \vdots \\ \xi_n \end{bmatrix} + B u \\
y = C \begin{bmatrix} \xi_1 \\ \xi_2 \\ \vdots \\ \xi_n \end{bmatrix} \qquad (8)
$$

### 3. DC Motor Speed Consensus Control

The DC motor model employed in this work was presented in [13]

$$ \left( \frac{d}{dt} \omega_l(t) \right) J_{eq} + B_{eq,v} \omega_l(t) = A_m V_m(t) \qquad (9) $$

with

$$ B_{eq,v} = \frac{\eta_g K_g^2 \eta_m k_t k_m + B_{eq} R_m}{R_m}, \quad A_m = \frac{\eta_g K_g \eta_m k_t}{R_m} $$

where $J_{eq} = \eta_g K_g^2 J_m + J_l$ and $B_{eq} = \eta_g K_g^2 B_m + B_l$. $\omega_l$ is the motor angular velocity. Motor parameters employed in (9) are: $R_m$ is the motor armature resistance, $k_m$ is the motor back-emf constant, $k_t$ is the motor current-torque constant, $K_g$ is the gear ratio, $\eta_m$ is the motor efficiency, $\eta_g$ is the gearbox efficiency, $J_m$ is the motor shaft moment of inertia, $J_l$ is the load moment of inertia, $B_m$ is the viscous friction on the motor shaft, and $B_l$ is the viscous friction on the load shaft. The final model is obtained

from (9)

$$ \dot{x}(t) = \alpha x(t) + \beta u(t) \qquad (10) $$

where $x(t) = \omega_l(t)$, $\alpha = -\frac{B_{eq,v}}{J_{eq}}$, $\beta = \frac{A_m}{J_{eq}}$, and $u(t) = V_m(t)$, which is designed to compensate the dynamics as follows:

$$ u(t) = \beta^{-1}(-\alpha x(t) + \bar{u}(t)) \qquad (11) $$

where $\bar{u}(t)$ is the control law obtained in the following sections.

Although multi-agent consensus control can be used to synchronize any number of DC motors, for simplicity, we illustrate the application of the consensus control on a three DC-motor system communicating over a network with cyclic topology.

#### 3.1. Three Agents Cyclic Topology

The particular case of a three agents system with a cyclic topology, as in Figure 1a, is analyzed below. Thus, for this case, Equation (1) simplifies to

$$ \begin{bmatrix} \dot{\xi}_1 \\ \dot{\xi}_2 \\ \dot{\xi}_3 \end{bmatrix} = \begin{bmatrix} \bar{u}_1 \\ \bar{u}_2 \\ \bar{u}_3 \end{bmatrix} \qquad (12) $$

Now, according to the desired reference characteristics, we consider the situations of fixed reference and time variant reference.
## 3.2. Fixed Reference

Similarly to [28], based on Figure 1a, we use

$$ \bar{u}_1 = k_p (\xi_2 - \xi_1) + u_1 \quad (13a) $$

$$ \bar{u}_2 = k_p (\xi_3 - \xi_2) \quad (13b) $$

$$ \bar{u}_3 = k_p (\xi_1 - \xi_3) \quad (13c) $$

so the resulting state space representation is

$$ \begin{bmatrix} \dot{\xi}_1 \\ \dot{\xi}_2 \\ \dot{\xi}_3 \end{bmatrix} = k_p \begin{bmatrix} -1 & 1 & 0 \\ 0 & -1 & 1 \\ 1 & 0 & -1 \end{bmatrix} \begin{bmatrix} \xi_1 \\ \xi_2 \\ \xi_3 \end{bmatrix} + \begin{bmatrix} u_1 \\ 0 \\ 0 \end{bmatrix} \quad (14) $$

The resulting matrix $L_r$ has an eigenvector $\omega^T = [1 1 1]$. Pre multiplying (14) by $\omega^T$,

$$ \dot{\xi}_1 + \dot{\xi}_2 + \dot{\xi}_3 = u_1 \quad (15) $$

The center of mass of the states is

$$ \xi_{cm} = \frac{\xi_1 + \xi_2 + \xi_3}{3} \quad (16) $$

Then, we take the leader control law (5) with $n = 3$

$$ u_1 = 3k_{cm} (\xi_{cm}^d - \xi_{cm}) \quad (17) $$

where $\xi_{cm}^d$ is the desired center of mass value and $k_{cm}$ is a constant. By substituting (17) in (15), we obtain

$$ \dot{\xi}_{cm} = -k_{cm} \xi_{cm} + k_{cm} \xi_{cm}^d \quad (18) $$

If $\xi_{cm}^d$ is a constant and $k_{cm} = 1$, then $(\xi_{cm} - \xi_{cm}^d) \to 0$ as $t \to \infty$. By substituting (17) in (14), we obtain

$$ \begin{bmatrix} \dot{\xi}_1 \\ \dot{\xi}_2 \\ \dot{\xi}_3 \end{bmatrix} = \begin{bmatrix} -2 & 0 & -1 \\ 0 & -1 & 1 \\ 1 & 0 & -1 \end{bmatrix} \begin{bmatrix} \xi_1 \\ \xi_2 \\ \xi_3 \end{bmatrix} + \begin{bmatrix} 3 \\ 0 \\ 0 \end{bmatrix} \xi_{cm}^d \quad (19) $$

From (19), it is possible to see that $\xi_1, \xi_2$ and $\xi_3$ converges to $\xi_{cm}^d$.

**Remark 1.** It is worth to recall that since the leader control is based on the center of mass of the agent's states, the general case for DC motor synchronization can be achieved by using control law (5).

## 3.3. Time Varying Reference

When a time varying reference is considered, it is necessary to include additional controls in all the agents besides the leader

$$ \dot{\xi}_1 = k_p (\xi_2 - \xi_1) + u_1 \quad (20a) $$

$$ \dot{\xi}_2 = k_p (\xi_3 - \xi_2) + u_2 \quad (20b) $$

$$ \dot{\xi}_3 = k_p (\xi_1 - \xi_3) + u_3 \quad (20c) $$
where

$$u_1 = k_l (\xi^d - \xi_2) + \dot{\xi}^d \qquad (21a)$$

$$u_2 = \dot{\xi}_3 \qquad (21b)$$

$$u_3 = \dot{\xi}_1 \qquad (21c)$$

by substituting (21) in (20), we obtain

$$\dot{\xi}_1 = k_p (\xi_2 - \xi_1) + k_l (\xi^d - \xi_2) + \dot{\xi}^d \qquad (22a)$$

$$\dot{\xi}_2 = k_p (\xi_3 - \xi_2) + \dot{\xi}_3 \qquad (22b)$$

$$\dot{\xi}_3 = k_p (\xi_1 - \xi_3) + \dot{\xi}_1 \qquad (22c)$$

again, by adding (22a), (22b), and (22c), we obtain

$$\dot{\xi}_2 - \dot{\xi}^d = k_l (\xi^d - \xi_2) \qquad (23)$$

Equation (23) shows that $(\xi_2 - \xi^d) \to 0$ as $t \to \infty$. From Equation (22b),

$$\dot{\xi}_2 - \dot{\xi}_3 = k_p (\xi_3 - \xi_2) \qquad (24)$$

then $(\xi_2 - \xi_3) \to 0$ as $t \to \infty$. Finally, from (22c),

$$\dot{\xi}_3 - \dot{\xi}_1 = k_p (\xi_1 - \xi_3) \qquad (25)$$

so we can see that $(\xi_3 - \xi_1) \to 0$ as $t \to \infty$.

To generalize this procedure for N motors, modify (20) as follows:

$$\dot{\xi}_1 = k_p (\xi_2 - \xi_1) + u_1 \qquad (26a)$$

$$\dot{\xi}_2 = k_p (\xi_3 - \xi_2) + u_2 \qquad (26b)$$

$\vdots$

$$\dot{\xi}_N = k_p (\xi_1 - \xi_N) + u_N \qquad (26c)$$

where $u_1$ is the same as in (21a), $u_k = \dot{\xi}_{k+1}$ for $k = 2, \dots, N-1$ and $u_N = \xi_1$, considering a cyclic case.

## 4. Robust Consensus Analysis

Equation (14) is the mathematical model of the consensus system for the particular case of three DC motors considering a cyclic topology; however, it is worth noting that the method presented in this section is valid for any topology that could be considered. In real life, it does not exist a perfect mathematical model that represents the dynamic performance of a physical system, so different areas of control theory, such as robust control, seek to compensate for modeling errors. Robust control considers the uncertainty in the mathematical model, taking it into account in the analysis stage so that the results are closer to reality than in the case when uncertainty is not considered. The state of the multi-agent system (14) represents the information exchange of the variables that must be synchronized, which, in our case, are the speeds of the DC motors. For the above, it is essential to use sensor devices frequently affected by the measurement noise. This phenomenon is modeled considering parametric uncertainty in the coefficients of the Laplacian matrix, as shown below, see [29]

$$\dot{\xi} = \mathcal{L}(q)\xi + Bu_1 \qquad (27)$$
where

$$ \mathcal{L}(\mathbf{q}) = \begin{bmatrix} q_{11} & q_{12} & 0 \\ 0 & q_{22} & q_{23} \\ q_{31} & 0 & q_{33} \end{bmatrix} \qquad B = \begin{bmatrix} 1 \\ 0 \\ 0 \end{bmatrix} \qquad \mathbf{q} = \begin{bmatrix} q_{11} \\ q_{12} \\ \vdots \\ q_{33} \end{bmatrix} \in Q $$

and Q is defined as

$$ Q = \left\{ \mathbf{q} = \begin{bmatrix} q_{11} & \cdots & q_{33} \end{bmatrix}^T : q_{ij} \bar{\le} q_{ij} \le q_{ij}^{+} \right\} \qquad (28) $$

the sets defined as in (28) are generally called boxes because of the geometric shape they have. On the other hand, manipulating Equation (27), the following equation can be obtained:

$$ \mathcal{L}(\mathbf{q}) = \mathcal{L}(\mathbf{q}^0) + \mathcal{L}(\mathbf{r}) \qquad (29) $$

the vectors $\mathbf{q}^0$ and $\mathbf{r}$ are the nominal and uncertain parameters of the multi-agent system, respectively. Their structure is as follows:

$$ \mathbf{q}^0 = \begin{bmatrix} q_{11}^0 \\ q_{12}^0 \\ \vdots \\ q_{nn}^0 \end{bmatrix} \qquad \mathbf{r} = \begin{bmatrix} r_{11} \\ r_{12} \\ \vdots \\ r_{nn} \end{bmatrix} \quad \text{and} \quad r_{ij} \in [q_{ij} - q_{ij}^0, q_{ij}^0 - q_{ij}^0] $$

the new way of representing the parametric uncertainty in (29), allows for describing the system (27) as follows:

$$ \dot{\xi} = \mathcal{L}(\mathbf{q}^0)\xi + Bu_1 + \mathcal{L}(\mathbf{r})\xi \qquad (30) $$

As we can see in the above equation, the uncertainty was only considered in the multi-agent system because our main objective was to establish the consensus in this system and achieve synchronization in the DC motors' speed. The uncertainty in the DC motors' mathematical model was treated independently through the control law applied to each motor (11).

### 4.1. Robust Consensus for the Multi-Agent System

In this section, the Lyapunov method is going to use to find the robustness margins of the multi-agent system. We will start by defining the following Lyapunov function taking into account the representation of the multi-agent system shown in (30):

$$ V(\xi) = \xi^T P \xi \qquad (31) $$

where, as stated in the Lyapunov method, the matrix P is a real symmetric positive definite matrix. Taking V(ξ) as in (31), the derivative of V along the trajectories of the multi-agent system is given by

$$ \begin{aligned} \dot{V}(\xi) &= \xi^T P \dot{\xi} + \xi^T P \xi \\ &= \xi^T P [\mathcal{L}(\mathbf{q}^0)\xi + Bu_1 + \mathcal{L}(\mathbf{r})\xi] \\ &\quad + [\mathcal{L}(\mathbf{q}^0)\xi + Bu_1 + \mathcal{L}(\mathbf{r})\xi]^T P \xi \end{aligned} $$
Now, considering the control law as defined in (17), it is possible to obtain the following equation:

$$ 
\begin{aligned}
\dot{V}(\xi) &= \xi^T P [\mathcal{L}(\mathbf{q}^0)\xi + B u_1 + \mathcal{L}(\mathbf{r})\xi] \\
&\quad + [\mathcal{L}(\mathbf{q}^0)\xi + B u_1 + \mathcal{L}(\mathbf{r})\xi]^T P \xi \\
&= \xi^T [P (\mathcal{L}(\mathbf{q}^0) - Bk) + (\mathcal{L}(\mathbf{q}^0) - Bk)^T P] \xi \\
&\quad + 2\xi^T P \mathcal{L}(\mathbf{r})\xi.
\end{aligned}
\quad (32) $$

The term in Equation (32) represents an inner product so the Cauchy–Schwarz inequality can be applied to obtain the following expression, see [30]

$$ 
\begin{aligned}
2\xi^T P \mathcal{L}(\mathbf{r})\xi &= 2\lambda(P)\xi^T \mathcal{L}(\mathbf{r})\xi \\
&\le 2\lambda_{\max}(P)\xi^T \mathcal{L}(\mathbf{r})\xi \\
&\le 2\lambda_{\max}(P)\gamma\|\xi\|^2
\end{aligned}
$$

Here, $\gamma$ is the upper bound of uncertainty $\|\mathcal{L}(\mathbf{r})\| \le \gamma$ and $\lambda(P)$ are the eigenvalues of the matrix $P$. When we replace the previous terms in (32), the following condition is obtained:

$$ 
\begin{aligned}
\dot{V}(\xi) &\le \xi^T [P (\mathcal{L}(\mathbf{q}^0) - Bk) + (\mathcal{L}(\mathbf{q}^0) - Bk)^T P] \xi \\
&\quad + 2\lambda_{\max}(P)\gamma \|\xi\|^2 \\
&\le -\xi^T Q \xi + 2\lambda_{\max}(P)\gamma \|\xi\|^2
\end{aligned}
$$

where $Q$ is a real symmetric positive definite matrix and satisfies the following equation:

$$ P (\mathcal{L}(\mathbf{q}^0) - Bk) + (\mathcal{L}(\mathbf{q}^0) - Bk)^T P = -Q \quad (33) $$

Finally,

$$ 
\begin{aligned}
\dot{V}(\xi) &\le -\xi^T Q \xi + 2\lambda_{\max}(P)\gamma \|\xi\|^2 \\
&\le (-\lambda_{\min}(Q) + 2\lambda_{\max}(P)\gamma) \|\xi\|^2
\end{aligned}
$$

Therefore, the condition to ensure the stability property of the uncertain multi-agent system (30) is the following:

$$ \gamma \le \frac{\lambda_{\min}(Q)}{2\lambda_{\max}(P)} \quad (34) $$

In the previous equation, it can be seen that the robustness margin is a function of the $Q$ and $P$ matrices that solve the Equation (33). Now, considering the closed-loop multi-agent system from (19)

$$ \mathcal{L}(\mathbf{q}^0) - Bk = \begin{bmatrix} -2 & 0 & -1 \\ 0 & -1 & 1 \\ 1 & 0 & -1 \end{bmatrix} \quad (35) $$

and defining $Q = I$, the solution of Equation (33) is as follows:

$$ P = \begin{bmatrix} 0.3095 & 0.0714 & 0.1190 \\ 0.0714 & 0.5000 & 0.2143 \\ 0.1190 & 0.2143 & 0.5952 \end{bmatrix} \quad (36) $$
Therefore, it can be noted that the uncertainty tolerated by the closed-loop multi-agent system must be $||\mathcal{L}(\mathbf{r})|| < 0.6208$. This implies that the equivalent uncertainty supported by the system is equal to 70% in each parameter of the Laplacian matrix, that is,

$$ \mathcal{L}(\mathbf{q}) = \begin{bmatrix} q_{11} & q_{12} & 0 \\ 0 & q_{22} & q_{23} \\ q_{31} & 0 & q_{33} \end{bmatrix} \qquad (37) $$

where $q_{11} = q_{22} = q_{33} \in [-1.35, -0.65]$ and $q_{12} = q_{23} = q_{31} \in [0.65, 1.35]$. It is possible to appreciate this robustness margin in the experimental tests that will be presented later.

## 4.2. A Practical Example of Computing the Time Delay Margin in Communication between Agents

A fundamental aspect to achieve the consensus of multi-agent systems is communication between agents, which is why it is important to consider all the possible phenomena that occur in this process. An important point is the time delay that exists when one agent sends messages to another; see [5,6]; this can be represented by

$$ \dot{\xi}(t) = A_0\xi(t) + A_r\xi(t - \tau) + Bu_1(t) \qquad (38) $$

where

$$ A_0 = \begin{bmatrix} -1 & 0 & 0 \\ 0 & -1 & 0 \\ 0 & 0 & -1 \end{bmatrix} \qquad A_r = \begin{bmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 1.25 & 0 & 0 \end{bmatrix} \qquad B = \begin{bmatrix} 3 \\ 0 \\ 0 \end{bmatrix} $$

The structure of the multi-agent system previously presented is due to the topology selected in Section 3.1, and it is important to note that the element of the matrix $Ar(3, 1)$ is equal to 1.25 because uncertainty was considered in that parameter. However, it could have been considered in all the matrix terms (37) without affecting the methodology. The calculations would be similar to those made for this particular case of the DC motor synchronization. Now, considering the control law from (17), the following closed-loop multi-agent system is obtained:

$$ \dot{\xi}(t) = A_0\xi(t) + A_1\xi(t - \tau) + B\xi_{cm}^d \qquad (39) $$

where

$$ A_1 = A_r - Bk = \begin{bmatrix} -1 & 0 & -1 \\ 0 & 0 & 1 \\ 1.25 & 0 & 0 \end{bmatrix} $$

Once the structure of the system is well defined, the fundamental problem will be to calculate the maximum time delay that the multi-agent system supports without the stability property being lost. It is well known that the property of stability of (39) is determined by the roots of the following characteristic equation; see [31]

$$ \det(sI - A_0 - A_1 e^{-\tau s}) \qquad (40) $$

Once this equation is developed, it represents functions called quasipolynomials that have the following form:

$$ \det(sI - A_0 - A_1 e^{-\tau s}) = p(s, e^{-\tau s}) \qquad (41) $$

$$ p(s, e^{-\tau s}) = \sum_{k=0}^{3} p_k(s) e^{-k\tau s} \qquad (42) $$

$$ p(s, e^{-\tau s}) = p_0(s) + p_1(s)e^{-\tau s} + p_2(s)e^{-2\tau s} + p_3(s)e^{-3\tau s} \qquad (43) $$
where

$$p_0(s) = (s+1)^3$$

$$p_1(s) = (s+1)^2$$

$$p_2(s) = 1.25(s+1)$$

$$p_3(s) = 0$$

The multi-agent system (39) is stable if and only if the following condition is satisfied:

$$\det(sI - A_0 - A_1e^{-\tau s}) \neq 0 \quad \forall s \in \mathbb{C}_+ \quad (44)$$

where $\mathbb{C}_+$ denotes the set of complex numbers with real part greater or equal to zero. Like an extension of the stability concept, the quasipolynomials that satisfy the previous condition are called stable quasipolynomials. Due to difficulty in determining the roots of quasipolinomial functions, we will define the following associated polynomials:

$$p(s, T) = \sum_{k=0}^{n} p_k(s)(1 - T_s)^{2k}(1 + T_s)^{2n-2k} \quad (45)$$

$$p(s, e^{-\tau s}) = \sum_{k=0}^{n} p_k(s)e^{-k\tau s} \quad (46)$$

These associate polynomials have an important relation between their roots: suppose that $s^0 = j\omega$ for some value of $\omega^0 \ge 0$; then, $s^0 = j\omega^0$ is a root of the characteristic equation $p(s, e^{-\tau s})$ for some value of $\tau \ge 0$ if and only if $s^0 = j\omega^0$ is a root of $p(s, T)$ for some value of $T \ge 0$. This relation that exists between $p(s, e^{-\tau s})$ and $p(s, T)$ is obtained from a transformation that is carried out in the exponential function that represents the time delay; see [32], while the relation between $T$ and $\tau$ is determined as follows (see [33]):

$$T = \frac{\tan(\frac{\tau\omega_i}{4})}{\omega_i} \quad (47)$$

where $\omega_i \in W$; the set of W is defined as follows:

$$W = \{0 < \omega < \omega^* : p(j\omega^*, T) = 0 \text{ for some value of } T > 0\} \quad (48)$$

in the case of the multi-agent system (39), the associated polynomial is the following:

$$p(s, T) = \sum_{k=0}^{7} a_k(T)s^k \quad (49)$$

where

$$a_0(T) = 3.25$$

$$a_1(T) = 6.25 - T$$

$$a_2(T) = 11.5T^2 + 7T + 4$$

$$a_3(T) = -T^3 + 21.5T^2 + 12T + 1$$

$$a_4(T) = 3.25T^4 + 7T^3 + 16T^2 + 4T$$

$$a_5(T) = 6.25T^4 + 12T^3 + 6T^2$$

$$a_6(T) = 4T^4 + 4T^3$$

$$a_7(T) = T^4$$
Now, define the polynomial associated Hurwitz matrix as follows:

$$p(s, T) = a_7(T)s^7 + a_6(T)s^6 + \dots + a_1(T)s + a_0(T) \quad (50)$$

$$H[p(s, T)] = \begin{bmatrix} a_6(T) & a_4(T) & a_2(T) & a_0(T) & 0 & 0 & 0 \\ a_7(T) & a_5(T) & a_3(T) & a_1(T) & 0 & 0 & 0 \\ 0 & a_6(T) & a_4(T) & a_2(T) & a_0(T) & 0 & 0 \\ 0 & a_7(T) & a_5(T) & a_3(T) & a_1(T) & 0 & 0 \\ 0 & 0 & a_6(T) & a_4(T) & a_2(T) & a_0(T) & 0 \\ 0 & 0 & a_7(T) & a_5(T) & a_3(T) & a_1(T) & 0 \\ 0 & 0 & 0 & a_6(T) & a_4(T) & a_2(T) & a_0(T) \end{bmatrix} \quad (51)$$

and let $p(s, T)$ the polynomial associated with $p(s, e^{-\tau s})$ and consider $H[p(s, T)]$ as the Hurwitz matrix in the functioning of $T$ of the associate polynomial $p(s, T)$; if $p(s, e^{-\tau s})$ is stable for $\tau = 0$, then $p(s, e^{-\tau s})$ is stable for all $\tau \in [0, \infty)$ if and only if

$$\det\{H[p(s, T)] \neq 0\} \quad \forall T > 0 \quad (52)$$

This result shows stability conditions for any time delay; when this happens, it is said that the system is stable independent of the delay. However, our goal is to determine the maximum delay that the multi-agent system can support without losing the stability property. To do this, it is necessary to compute the values of $T$ where the determinant of the Hurwitz matrix $H[p(s, T)]$ becomes zero and using (47) and (48) compute the time delays for which the system loses the stability property. Carrying out the corresponding computation, we get the values of $T$ that make the determinant of the Hurwitz matrix (51) zero; these are $T_1 = 0.828427$ and $T_2 = 2.774851$ and their corresponding $\omega_1 = \omega_2 = 0.5$. These, in turn, are used to compute $\tau_1 = 3.1416$ and $\tau_2 = 7.5702$; thus, the maximum delay supported by the multi-agent system is $\tau_{max} = 3.1416$ s. The robustness analysis performed in this section guarantees that adequate performance is obtained in experimental testing of the DC motor synchronization.

## 5. Experimental Results

The agent speed consensus control addressed in Section 2 was implemented by using a Bluetooth wireless communication system and an embedded microcontroller-based control platform, as shown in Figure 2. This platform consists of three identical modules whose elements are: a DC motor managed by a dual-core ARM Cortex-A9 embedded microcontroller system, with the ability to be programmed via WiFi, a Bluetooth transmitter, and a receiver (115,200 baud), as well as an encoder to measure motor speed, with a sampling frequency of 25 samples per second. Each motor has a load consisting of a wheel. The control algorithm was programmed using LabVIEW®. Before the tests, the WiFi communication was used to configure the controllers parameters, but, during the tests, it was only used to send the reference value from the computer to the leader controller. The gains $k_p = 0.07$ and $k_{cm} = 1$ were chosen to have an underdamped response. Note that while controllers 2 and 3 need one speed state information via Bluetooth, the controller 1, the leader, needs to have all speeds to calculate the center of mass. Approximate motor parameter values are shown in Table 1. It is important to note that, in this system, there may be several sources of uncertainty, namely motor parameters, measurement noise, and transmission delays. A schematic diagram of the multi-agent system is shown in Figure 2, and the experimental implementation is presented in Figure 3.
<figure><img src="image_6.png" /><figcaption>Figure 2. Experimental agent interconnection diagram.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Figure 3. Three-motor synchronization experimental platform.</figcaption></figure>

<table><caption>Table 1. Motor parameters and values.</caption><tr><td>R<sub>m</sub></td><td>3.3 Ω</td><td>B<sub>l</sub></td><td>4.41 × 10<sup>-6</sup> Nm/(rad/s)</td></tr><tr><td>k<sub>m</sub></td><td>7.68 × 10<sup>-3</sup> V/(rad/s)</td><td>k<sub>t</sub></td><td>6.9 × 10<sup>-3</sup> Nm/A</td></tr><tr><td>k<sub>g</sub></td><td>53</td><td>η<sub>m</sub></td><td>0.69</td></tr><tr><td>η<sub>g</sub></td><td>0.90</td><td>J<sub>m</sub></td><td>3.90 × 10<sup>-7</sup> kgm<sup>2</sup></td></tr><tr><td>J<sub>l</sub></td><td>1.03 × 10<sup>-4</sup> kgm<sup>2</sup></td><td>B<sub>m</sub></td><td>1.07 × 10<sup>-9</sup> Nm/(rad/s)</td></tr></table>

The tests consisted of a step reference change and a sinusoidal reference change, $r(t) = 70 + 30\sin(0.2\pi t)$. The response of the DC motor speeds ($\omega_1(0) = 60.13$ rpm, $\omega_2(0) = 59.75$ rpm and $\omega_3(0) = 59.67$ rpm) and controls to a reference step change are shown in Figures 4 and 5, respectively. It is possible to note that the leader control $u_1$ has a bigger transitory because it is responsible for carrying the system to the reference; the other two controls carry the states to the center of mass. The behavior of the center of mass is plotted in Figure 6, where the center of mass is very close to the reference during stationary state and smoother than the speeds, acting as a low pass filter. This variation, of approximately 0.5 rpm, is mainly due to the encoder used to measure speed. Figures 7–9 show the response, control outputs of the motor system, and the center of mass with a time-varying reference ($\omega_1(0) = 0.40$ rpm, $\omega_2(0) = 1.10$ rpm and $\omega_3(0) = 1.14$ rpm). The synchronization of angular speeds is appreciated under reference step change and sinusoidal reference. Observe that small errors appear when using time-varying references, due to the propagation delay of the reference over the wireless communication system and the states, working as uncertainty, but is not enough to avoid synchronization.
<figure><img src="image_8.png" /><figcaption>Figure 4. Experimental speed consensus with a reference step change.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>Figure 5. Experimental control outputs with a reference step change.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>Figure 6. Experimental state center of mass with a reference step change.</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 7.** Experimental speed tracking using a sinusoidal reference.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 8.** Experimental control outputs for each agent following a sinusoidal reference.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 9.** Experimental center of mass of agents following a sinusoidal reference.</figcaption></figure>
## 6. Conclusions and Future Work

In this article, an application of a control methodology that uses a multi-agent consensus theory was presented. The convergence of the states to both a fixed and a time-varying reference was demonstrated on a three-agents system with cyclic information flow topology. This methodology was applied in order to synchronize the speed of three DC motors. Through a robustness analysis of the uncertainty of the parameters of DC motors and the delays in communication between agents, the consensus was ensured. In addition, experimental results on a three DC motor system with a wireless communication system were presented. Most of the reviewed works presented just two motors. It was implemented by using a Bluetooth wireless agent communication system and an embedded micro-controller based control platform. It was noticed that the proposed control approach allows us to synchronize the multiple DC motor system. The small agent consensus error concerning the desired reference is due, besides the motor dynamics, to wireless communication systems, and the information flow topology.

From the information above, we identify three different research lines that may be addressed as future work; the first is to perform an independent robustness analysis in each of the parameters of the *Laplacian matrix* prioritizing the parameters that may present most variation. The second research line is to consider in the analysis of the time delay system a higher number of parameters with uncertainty since this paper only considers one of them. Regarding the latter, it is clear that the maximum time delay that the close loop system may withstand without losing the stability property will be smaller than the one obtained when considering only one parameter. Finally, another future work is to perform experimental testing considering different topologies of information exchange such as chain topology or balances graph topology, adapting the experimental platform for this purpose.

**Author Contributions:** Conceptualization, D.O., G.R., J.A.G., and R.L.; methodology, G.R. and J.A.G.; software, D.O., G.R., and J.A.G.; validation, D.O.; formal analysis, D.O., G.R., J.A.G., and R.L.; investigation, D.O., G.R., J.A.G., and R.L.; resources, D.O., G.R., J.A.G., and R.L.; writing—original draft preparation, D.O., G.R., and J.A.G.; writing—review and editing, D.O., G.R., and J.A.G.; visualization, D.O., G.R., and J.A.G.; supervision, G.R. and J.A.G.; project administration, G.R.; funding acquisition, D.O. and G.R. All authors have read and agreed to the published version of the manuscript.

**Funding:** The work of Gerardo Romero was supported by Universidad Autónoma de Tamaulipas, under grants PROINNOVA-2018-250105, PROINNOVA-2018-250586, PROINNOVA-2018-250540, PROINNOVA-2018-250113, PROINNOVA-2018-250117, PROINNOVA-2018-250160, and PROINNOVA-2018-250255.

**Acknowledgments:** Daniel Olivares appreciates the support and use of facilities of the Tecnológico Nacional de México—Instituto Tecnológico de Nuevo Laredo.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Leonard, N.E.; Fiorelli, E. Virtual leaders, artificial potentials and coordinated control of groups. In Proceedings of the 40th IEEE Conference on Decision and Control (Cat. No.01CH37228), Orlando, FL, USA, 4–7 December 2001; pp. 2968–2973.

2. Arrichiello, F.; Chiaverini, S.; Fossen, T. Formation Control of Marine Surface Vessels Using the Null-Space-Based Behavioral Control. In *Lecture Notes in Control and Information Sciences: Group Coordination and Cooperative Control*; Springer: Berlin/Heidelberg, Germany, 2006; pp. 1–19.

3. Chen, X.; Serrani, A. ISS-Based Robust Leader/Follower Trailing Control. In *Lecture Notes in Control and Information Sciences: Group Coordination and Cooperative Control*; Springer: Berlin/Heidelberg, Germany, 2006; pp. 75–91.

4. Lozano, R.; Spong, M.W.; Guerrero, J.A.; Chopra, N. Controllability and observability of leader-based multi-agent systems. In Proceedings of the 47th IEEE Conference on Decision and Control, Cancun, Mexico, 9–11 December 2008; pp. 3713–3718.

5. Olfati-Saber, R.; Fax, J.A.; Murray, R.M. Consensus and Cooperation in Networked Multi-Agent Systems. *Proc. IEEE* **2007**, *95*, 215–233.
6. Ren, W.; Beard, R.W.; Atkins, E.M. Information Consensus in Multivehicle Cooperative Control. *IEEE Control Syst. Mag.* **2007**, *27*, 71–82.

7. Guerrero, J.A.; Olivares, D.; Romero, G. Forced Bipartite Consensus for Multi-Agent Systems. In Proceedings of the 44th Annual Conference of the IEEE Industrial Electronics Society (IECON), Washington, DC, USA, 21–23 October 2018; pp. 2335–2340.

8. Nagata, T.; Sasaki, H. A multi-agent approach to power system restoration. *IEEE Trans. Power Syst.* **2002**, *17*, 457–462. [CrossRef]

9. Balch, T.; Arkin, R.C. Behavior-based formation control for multirobot teams. *IEEE Trans. Robot. Autom.* **1998**, *14*, 926–939. [CrossRef]

10. Guerrero, J.A.; Lozano, R. Flight Formation of Multiple Mini Rotorcraft based on Nested Saturations. In Proceedings of the 2010 IEEE/RSJ International Conference on Intelligent Robots and Systems, Taipei, Taiwan, 18–22 October 2010; pp. 634–639.

11. Chuang, C.W.; Haung, C.L.; Lee, C.D.; Kao, C.C.; Fung, R.F. Synchronization and tension control of dual motor systems via MIMO discrete pseudo model following integral variable structure control. *Mech. Mach. Theory* **2009**, *44*, 499–510. [CrossRef]

12. Chiu, G.T.C.; Tomizuka, M. Coordinated Position Control of Multi-Axis Mechanical Systems. *J. Dyn. Syst. Meas.* **Control** **1998**, *12*, 389–393. [CrossRef]

13. Guzey, H.C.; Dumlu, A.; Guzey, N.; Alpay, A. Optimal Synchronizing Speed Control of Multiple DC Motors. In Proceedings of the 4th International Conference on Optimization and Applications, Mohammedia, Morocco, 26–27 April 2018; pp. 1–5.

14. Li, L.; Sun, L.; Zhang, S.; Yang, Q. Speed tracking and synchronization of multiple motors using ring coupling control and adaptive sliding mode control. *ISA Trans.* **2015**, *58*, 635–649. [CrossRef] [PubMed]

15. Zhao, D.; Li, C.; Ren, J. Speed Synchronization of Multiple Induction Motors with Total Sliding Mode Control. *Syst. Eng. Theory Pract.* **2009**, *29*, 110–117. [CrossRef]

16. Zhang, Z.; Chau, K.T.; Wang, Z. Chaotic Speed Synchronization Control of Multiple Induction Motors Using Stator Flux Regulation. *IEEE Trans. Magn.* **2012**, *48*, 4487–4490. [CrossRef]

17. Kose, E.; Muhurcu, A. The Control of Brushless DC Motor for Electric Vehicle by Using Chaotic Synchronization Method. *Stud. Inform. Control* **2018**, *27*, 403–412. [CrossRef]

18. Bhattacharyya, S.P. Robust control under parametric uncertainty: An overview and recent results. *Annu. Rev. Control* **2017**, *44*, 45–77. [CrossRef]

19. Matusu, R.; Senol, B.; Pekar, L. Value-Set-Based Approach to Robust Stability Analysis for Ellipsoidal Families of Fractional-Order Polynomials with Complicated Uncertainty Structure. *Appl. Sci.* **2019**, *9*, 5451. [CrossRef]

20. Sun, Y.-G.; Xie, S.; Xu, J.-Q.; Lin, G.-B. A Robust Levitation Control of Maglev Vehicles Subject to Time Delay and Disturbances: Design and Hardware Experimentation. *Appl. Sci.* **2020**, *10*, 1179. [CrossRef]

21. Wu, L.; Zhao, R.; Li, Y.; Chen, Y.-H. Optimal Design of Adaptive Robust Control for the Delta Robot with Uncertainty: Fuzzy Set-Based Approach. *Appl. Sci.* **2020**, *10*, 3472. [CrossRef]

22. Trentelman, H.L.; Takaba, K.; Monshizadeh, N. Robust Synchronization of Uncertain Linear Multi-Agent Systems. *IEEE Trans. Autom. Control* **2013**, *58*, 1511–1523. [CrossRef]

23. Wang, X.; Li, S. Distributed finite-time consensus algorithms for leader–follower second-order multi-agent systems with mismatched disturbances. In Proceedings of the IEEE American Control Conference, Boston, MA, USA, 6–8 July 2016; pp. 2814–2819.

24. Huang, M.; Manton, J.H. Coordination and Consensus of Networked Agents with Noisy Measurements: Stochastic Algorithms and Asymptotic Behavior. *SIAM J. Control Optim.* **2009**, *8*, 134–161. [CrossRef]

25. Guerrero, J.A.; Romero, G.; Lozano, R. Robust consensus tracking of leader-based multi-agent systems. In Proceedings of the 2010 American Control Conference, Baltimore, MD, USA, 30 June–2 July 2010; pp. 6299–6305.

26. Godsil, C.; Royle, G. *Graduate Texts in Mathematics: Algebraic Graph Theory*; Springer: Berlin, Germany, 2001.

27. Mesbahi, M.; Egerstedt, M. *Graph Theoretic Methods in Multiagent Networks*; Princeton University Press: Princeton, NJ, USA, 2010.

28. Guerrero, J.A.; Castillo, P.; Challal, Y. Trajectory Tracking for a Group of Mini Rotorcraft Flying in Formation. In Proceedings of the IFAC World Congress, Milano, Italy, 27 August–2 September 2011; pp. 6331–6336.
29. Romero, G.; Guerrero, J.A.; Reyes, L.D.; Lozano, R.; Olivares, D. Robustness Margin for Leader-based Multi-agent Consensus Systems in Presence of Parametric Uncertainty. In Proceedings of the 44th Annual Conference of the IEEE Industrial Electronics Society (IECON), Washington, DC, USA, 21–23 October 2018; pp. 2625–2630.

30. Luenberger, D.G. *Optimization by Vector Space Methods*; John Wiley & Sons, Inc.: New York, NY, USA, 1969.

31. Romero, G.; Diaz, I.; Perez, I.; Guerrero, J.A.; Lara, D.; Rivera, J. New results on robust stability for differential-difference systems with affine linear parametric uncertainty. *Int. J. Syst. Sci.* **2013**, *44*, 14–22.
[CrossRef]

32. Thowsen, A. An Analytic Stability Test for a Class of Time-Delay Systems. *IEEE Trans. Autom. Control* **1981**, *26*, 735–736. [CrossRef]

33. Thowsen, A. The Routh-Hurwitz Method for Stability Test for a Class of Time-Delay Systems. *Int. J. Control* **1981**, *33*, 991–995. [CrossRef]

<figure><img src="image_14.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).