Mathematical Problems in Engineering

Volume 2016 (2016), Article ID 3567682, 10 pages

http://dx.doi.org/10.1155/2016/3567682

## Distributed Control for Multiagent Consensus Motions with Nonuniform Time Delays

School of Aeronautics and Astronautics, University of Electronic Science and Technology of China, Chengdu 611731, China

Received 23 February 2016; Accepted 19 May 2016

Academic Editor: Zhan Shu

Copyright © 2016 Mengji Shi and Kaiyu Qin. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This paper solves control problems of agents achieving consensus motions in presence of nonuniform time delays by obtaining the maximal tolerable delay value. Two types of consensus motions are considered: the rectilinear motion and the rotational motion. Unlike former results, this paper has remarkably reduced conservativeness of the consensus conditions provided in such form: for each system, if all the nonuniform time delays are bounded by the maximal tolerable delay value which is referred to as “delay margin,” the system will achieve consensus motion; otherwise, if all the delays exceed the delay margin, the system will be unstable. When discussing the system which is intended to achieve rotational consensus motion, an expanded system whose state variables are real numbers (those of the original system are complex numbers) is introduced, and corresponding consensus condition is given also in the form of delay margin. Numerical examples are provided to illustrate the results.

#### 1. Introduction

As one of the fundamental problems in the study of the cooperative control for multiagent systems, consensus problem has received a fair amount of recent attention, because of its applications in formation control [1], target tracking [2], and sensor networks [3]. When researching consensus problems, time delay is hardly a matter of rarity, and numerous examples [4–16] have attested to its prevalence. It then raises a question: are the consensus conditions for the nondelayed systems still applicable in presence of time delays?

The treatments of systems with time delays have been widely studied, and researchers have developed diverse applicable approaches to analyzing the system stability under the occurrence of time delays, such as Lyapunov approach, the approach based on the properties of nonnegative matrices, and frequency domain approach. Lyapunov approach is a powerful tool of time domain analysis, with which one can easily handle nonlinearities and time-varying uncertainties. The consensus problems of delayed systems in [4, 5, 7, 8, 10, 12, 13, 16] have been worked out by establishing Lyapunov functions (functionals). Approaches based on the properties of nonnegative matrices can be available to the cases of dynamically changing topologies [15], but they are commonly limited to discrete-time systems. Linear delay systems are discussed in greater length with frequency domain approaches with which [6, 9, 11] have studied the stability of multiagent systems with time delays, and necessary and sufficient conditions for these systems to reach consensus have been derived. Usually, consensus conditions derived by Lyapunov-based approaches are more likely to be conservative than the ones obtained by frequency domain approaches.

Taking nonuniform time delays into account, the paper investigates the consensus problems of multiagent systems in terms of motion. We consider two types of consensus motions: one is rectilinear motion, and the other is rotational motion. In the two systems, the agents aim to collectively move at a uniform velocity and angular speed, respectively, depending on their initial states. These multiagent models have found various applications, such as formation control for unmanned aerial vehicles (UAVs), spacecraft docking, and circular mobile sensor networks. The second order control algorithms are extended from the previous works [11, 17], but nonuniform time delays are considered.

An efficient frequency domain approach is employed to derived sufficient conditions for these consensus problems. The results are presented in the form of “delay margins.” Each delay margin is a derived value that if all the time delays in the corresponding system are bounded by this value the system will reach consensus, but if all the delays exceed this value the agents will not achieve collective motion. Consequently, one cannot find out a greater value that if all the delays are bounded by that value the system still reaches consensus always. In [18], the author discussed the stabilities of second-order multiagent systems with nonuniform delays using Lyapunov methods, and the results are presented in the form of Linear Matrix Inequalities (LMIs), which are different from our work. The previous publications [11, 19] have studied multiagent systems with unstable modes, but nonuniform delays have not been considered. Reference [17] has proposed the rotating consensus algorithm, without consideration of delay effects. The rotating consensus algorithm we adopt enables the agents to achieve rotational consensus motion in the complex plane. However, in many cases, complex numbers cannot be processed by the processor. We expand the original system by separating the real part and imaginary part of each state variable. The new system has real state variables and coefficient matrices. We finally provide consensus conditions for the expanded system in the form of delay margin.

The remainder of this paper is organized as follows. Section 2 presents some mathematical preliminaries for the later proof, and with their help, we define the consensus problems and then protocols are given for the second-order multiagent systems with multiple time delays. Section 3 shows the stability analysis of the systems and discusses what will happen to the system to achieve rectilinear consensus motion if the input delays are absent; this section also provides an expansion to the system to reach rotating consensus, and this additional result allows us to simulate the algorithms in an easier way. Section 4 gives numerical simulations. Finally in Section 5, some conclusions are drawn.

#### 2. Model

Before presenting the model, we first provide a brief narrative on some definition and results in graph theory, which are used extensively in the rest of this paper.

Consider a multiagent system consisting of agents. The communication network topology among them is represented by an undirected graph , which consists of a set of nodes , a set of edges , and a weighted adjacency matrix , where , and (for is undirected). if and only if there exists an edge between the th and th nodes, which implies that they can get information from each other. The set of neighbours of node is denoted by . The Laplacian associated with the graph is defined as , where and . A path is a sequence of ordered edges , where . If there is a path between every pair of nodes in graph , the graph is said to be connected. The lemma below is given by [20].

Lemma 1. *If the undirected graph is connected, then its Laplacian has one singleton zero eigenvalue (with eigenvector ), and the rest of eigenvalues of are all positive.*

Consider -agent systems with second-order dynamics:where , denote the position and velocity of the th agent, respectively, and is the control input. Denoting the state vector , we assume that the initial condition is , for .

Two goals have been set for the system to achieve: rectilinear consensus motion and rotational consensus motion. Clearly, the control input decides the behavior of the agents. We will propose two algorithms, respectively, for the two goals; before that, we use mathematical symbols to represent what will happen to the agents when the goals are achieved.

For the first goal, rectilinear consensus motion, if is satisfied for all , the control input is said to solve the consensus problem asymptotically [11]. System (1) under this control input is referred to as “system ” hereinafter. In fact, in some cases (very extreme cases), the agents will not move collectively but all stay at a point as a final state; it depends on the initial conditions.

For the second goal, rotational consensus motion, the consensus is reached if for all , where is the constant angular velocity in which all agents move around a common point together [17]. We assume to simplify the future discussion. System (1) under a control input to satisfy (2) is called “system ” hereunder.

The state variables of agents in system are considered in the real domain for convenience of calculation; that is, in . However, the agents in system are intended to achieve rotational consensus motion in the complex plane; thus the state variables are complex values for system .

Taking nonuniform time delays into account, we modify the corresponding protocols proposed in [11, 17] and give the control inputs of our systems and asfor any , where denotes the edge weight, and is the time delay for the th agent to get the state information of the th agent. We also consider the system to have different time delays, denoted by .

We use matrix notation to write the models of both systems in compact form. To unify the form, we define where for all . and are the state vectors of system and , respectively. By defining coefficient matrices system (1) under the control inputs (3) can be rewritten asfor , with the initial condition , where denotes the Laplacian of a subgraph associated with the delay . Clearly, .

In this paper, the graph is assumed to be connected and undirected. Then according to Lemma 1, we suppose that the eigenvalues of are .

#### 3. Main Results

In this section, we derive consensus conditions for systems , with a frequency domain approach. The main idea of this method is to firstly prove the consensus achievement of each system when all the time delays are equal to zero and then to calculate the delay margin to ensure that if all the delays are under the bound of it, the stability of the system will not be affected.

##### 3.1. Consensus Conditions for the Original Systems

Firstly, we investigate the stability of in absence of time delays. Let and then given by (6) () with zero time delays can be written as

Lemma 2. *The matrix has a double eigenvalue at zero with a single Jordan block, and all other eigenvalues have negative real parts.*

*Proof. *According to Lemma 1, there exists an orthogonal matrix , such that and then it follows that Through simple calculations, we haveThe roots of each () areSince when , if , the real parts of (12) are , and when , (12) apparently are negative real numbers. Thus the polynomial (11) has a double root at zero, and the real part of all other roots is negative; then the lemma is proven.

*Remark 3. *According to Lemma 2, system (8) is unstable in the sense of traditional control due to the existence of the zero Jordan block in the system matrix. However, one can easily obtain that the vectors and are the eigenvectors associated with the zero eigenvalue. Therefore the nondelayed system (8) finally converges to the space spanned by the two vectors; that is, , for all . This has illustrated that the consensus is reached.

Reference [17] has proven that system will achieve rotational consensus motion when all the time delays are absent. Then we summarize the results in the following lemma.

Lemma 4. *Consider a network of second-order agents with fixed topology that is connected. Both of the multiagent systems () given by (6) will reach consensus if all the delays .*

The next step is to find out the conditions ensuring stable consensus motions of both systems.

Theorem 5. *Consider the second-order systems , that, given by (6), the consensus problems are solved, if in each system all the delays are smaller than a corresponding value , where *

*Proof. *Consider the multiagent systems (6) with nonuniform time delays. Take as the Laplace transform of ; we have , where for . According to the foregoing discussions, to study the stability of the delayed system, we only need to investigate the values of that guarantee the existence of nonzero roots of on imaginary axis, which represents the crossing of the characteristic roots from the stable region to the unstable one. The roots of characteristic polynomials such as are hereinafter referred to as “the eigenvalues of the system.”

Suppose that, for each , is an imaginary root of , and is a corresponding eigenvector, where . Then we haveNote that all the complex roots of each appeared in conjugated pairs; we only need to study the situation that . Since the elements of the vector obtained by calculating the left part of (15) equal zero, we haveMultiplying by (the conjugate transpose of ) on the left side of the left part of (15), and with (16) substituted, it yieldswhere Equation (17) can be rewritten asTake modulus of both sides of (19) to obtainfor . After some calculation, we find that the inequality from (20) holds true if and only if . This is a necessary condition for our hypothesis to hold.

The idea now is to analyze the argument of each , to find out how affects the validity of our hypothesis that the imaginary eigenvalues exist. Calculate the arguments of ; we get where . Define . Then the following proves that, for each , is the minimal value of when . For system , Let . Calculate derivative of : and the derivative of is Thus when . Consequently, when , which means decreases as grows. Therefore we have , when . For system , when , one can easily calculate that always holds; and when , is monotonically decreasing from ; moreover, when , (we omit the proof which is similar to the one given above). Because the network of the system is connected, holds true; it follows that . Therefore, one cannot find a smaller value of than when .

Note that are also the arguments of the left parts of (19); if for each system we set all , we have for , which conflicts with the idea that is the minimal value of when ; and this rejects our hypothesis. Consequently, each of the two systems has no nonzero imaginary eigenvalue when all the delays . Following the continuity of function, the zero crossing of eigenvalues from the left half plane to the right half plane has been prevented, and the bounded delays will not bring about instabilities to the systems. According to Lemma 4, the system can reach consensus with bounded delays, and the theorem is proven.

In fact, for each system, if all , then there must exist an imaginary eigenvalue and corresponding eigenvector that satisfy . That is, according to (19),holds true.

*Remark 6. * are termed the delay margins. If we draw trajectories of the delays that each value grows from 0 to , then the touch of leads to the first contact of the eigenvalues from the stable region to the unstable one. When all the delays exceed , the system will be unstable, because there must exist at least one eigenvalue in the unstable region (the RHP). The instabilities get worse with increasing .

*Remark 7. *Reference [9] has provided a general method of obtaining delay margins for delayed systems (including nonuniform delays) to solve consensus problems, and the method proposed in this paper is inspired by that idea. However, only first-order systems have been considered in [9], where the arguments of are bounded by the simple . In our paper, the arguments are bounded by variables changing with ; thus we have to find possible maximal delay values related to to cause the contradictions. Reference [21] has discussed the consensus problem of multiagent systems with delays using frequency approaches, but the agent dynamics are only considered in the form of single integrator, and the methods can not be directly used in second-order systems. In the literatures [11, 22], the authors studied the consensus problems of second-order systems with uniform delays; however, the situations of multiple time delays were not discussed. In [18], the stability analyses of second-order multiagent systems with nonuniform delays were performed using Lyapunov methods, and the results were presented in the form of Linear Matrix Inequalities (LMIs).

##### 3.2. Consensus Phenomenon of System without Input Delays

In many cases, the consensus is reached when all agents converge on a point and stay. References [11, 19] have proposed consensus algorithms enabling the vehicles to drive together at a uniform velocity after their convergence. These classical algorithms are then extensively applied in the subsequent works, such as the generalized use in [23] to solve the tracking problem. Therefore exploring such basic algorithms is very useful and the results can be generalized for many applications. In this paper, the idea of comes from the algorithm proposed in [11]; but we consider nonuniform time delays. In our systems, each agent gets delayed information from its neighbors, and besides these communication delays, input delays also occur. Consequently, one may come up with a question that what if the input delays are absent while communication delays still occur. In this situation, the control input of (1) for becomes

Suppose that there is a time when all the agents are moving forward together at the same point; then the value of depends on the value of and how fast the agents move. If this happens, (26), which determines the acceleration of each agent, will not equal zero, and it is difficult for and () to reach the same value due to the nonuniform delays. Then, in the next period of time, the agents have to separate from each other because of their different accelerations.

Simulations (see in Section 4) have shown a reasonable consensus phenomenon when applying (26) to the system: when appropriate delays are set to the system, the agents move to get together at first, and then slow down until they stop at the same point, which means that the agents will not finally achieve a collective motion when the input delays are absent while communication delays still occur.

##### 3.3. Conditions for an Expanded System to Achieve Rotational Consensus Motion

System given by (6) () has provided a solution that the agents are intended to perform a rotational consensus motion in the complex plane. However, in many practical cases, the processors of the systems do not accept complex numbers directly; as an example, it is difficult to draw a simulation model in the MATLAB for such a system involving complex state variables. Commonly, we “expand” the system by replacing each complex state variable with two real variables: the real part and the imaginary part of the original variable, respectively. The new system is equipped with real coefficient matrices.

Denote the new state vector and coefficient matrices: where and mean the real and imaginary part of the term enclosed, respectively, and

We get the new system expanded from the system given by (6) ():Equation (29) is another form of , and its state variables are real values. Agents in system (29) behave like those in given by (6) () when all the time delays are absent. Therefore, in the light of the real state variables, we usually simulate the expanded system instead of the original one. However, there is a notable difference between the two systems. The amount of eigenvalues of the new system (29) doubles that of (6) (). We introduce another related system for reference:where and are the the conjugate matrices of and , respectively. The subsequent lemma illustrates the relationships among the three systems (6) (), (30), and (29).

Lemma 8. *There are eigenvalues of system (29), which consist of all the eigenvalues of (6) () and all the eigenvalues of (30).*

*Proof. *Suppose that is one eigenvalue of (6) (), and the corresponding eigenvector is . Then we haveSubstituting into the characteristic polynomial of system (29), and multiplying by , it yields which implies that is an eigenvalue of system (29), and a corresponding eigenvector. It follows that the eigenvalues of (6) () also are the eigenvalues of (29). By the same way, one can prove that each eigenvalue of system (30), whose corresponding eigenvector is here denoted as , is also an eigenvalue of (29), and can be a corresponding eigenvector. As the eigenvalues fit the system (29) in amount, and (29) has no other eigenvalues, the lemma is proven.

Lemma 9. *The expanded multiagent system (29) with real state variables will reach consensus if all the time delays are absent; that is, .*

*Proof. *When all the time delays equal zero, the eigenvalues of system (30) are the conjugate numbers of the eigenvalues of (6) (). As each conjugated pair has the same real parts, the expanded system (29) has the same stability properties as (6) () when all time delays are absent, and the agents will achieve rotational consensus motion.

Theorem 10. *Considering the multiagent system (29) expanded from (6) (), it will reach consensus if all the delays , where *

*Proof. *Theorem 5 has proven that is the delay margin ensuring that the eigenvalues of (6) () are in the stable region. We can also prove that the eigenvalues of system (30) are in the stable region if all the delays are bounded by , and the system will be unstable if all the delays exceed , by the same way as in the proof of Theorem 5. Thus the proof is omitted.

The foregoing discussions have helped us to understand that when all the delays , the eigenvalues of both systems (6) () and (30) are in the stable region. According to Lemmas 8 and 9, system (30) can reach consensus with nonuniform time delays bounded by .

*Remark 11. *The theorem above has provided us with a delay bound ensuring the consensus achievement. Note that both and are delay margins of the corresponding systems. We set all the time delays of the expanded system (29) greater than ; when these delays are put into the systems (6) () and (30), at least one of the systems will have eigenvalues in the unstable region; this will bring about the failure of the expanded system (29) to achieve consensus motion. Hence is also the delay margin for the expanded system.

#### 4. Numerical Examples

In this section, some simulations are carried out to evidence the theoretical results obtained by the previous analysis.

Consider that a multiagent system consists of 4 agents, whose communication topology is described in graph shown in Figure 1. The 1st-2nd and 3rd-4th agent pairs are connected and the communication delay between each of them is ; the 2nd-3rd and 4th-1st agent pairs are connected in presence of communication delay . Suppose the weight of each edge on is 1; then .