Mathematical Problems in Engineering

Volume 2017, Article ID 7307834, 12 pages

https://doi.org/10.1155/2017/7307834

## Consensus Conditions for High-Order Multiagent Systems with Nonuniform Delays

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

Correspondence should be addressed to Kaiyu Qin; nc.ude.ctseu@niqyk

Received 20 May 2016; Revised 8 October 2016; Accepted 16 October 2016; Published 9 April 2017

Academic Editor: Guido Ala

Copyright © 2017 Mengji Shi et al. 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

Consensus of first-order and second-order multiagent systems has been wildly studied. However, the convergence of high-order (especially the third-order to the sixth-order) state variables is also ubiquitous in various fields. The paper handles consensus problems of high-order multiagent systems in the presence of multiple time delays. Obtained by a novel frequency domain approach which properly resolves the challenges associated with nonuniform time delays, the consensus conditions for the first-order and second-order systems are proven to be nonconservative, and those for the third-order to the sixth-order systems are provided in the form of simple inequalities. The method revealed in this article is applicable to arbitrary-order systems, and the results are less conservative than those based on Lyapunov approaches, because it roots in sufficient and necessary criteria of stabilities. Simulations are carried out to validate the theoretical results.

#### 1. Introduction

Consensus problems of multiagent systems have found many applications in the fields that hold great promise, including (but not limited to) biosciences, robotics, and computer sciences. Consensus is the agreement regarding a certain quality of interest on specific states of all the agents, which is widely demanded in the engineering applications. The research on consensus problems has lasted for decades. Various techniques are developed to solve consensus problems of numerous multiagent systems [1–23].

This paper addresses the consensus control problems of high-order multiagents systems with nonuniform time delays. One motivation for studying high-order systems is to achieve accurate control of complex motion: for example, when performing consensus motion that requires abrupt change of heading, a team of vehicles should maintain consistency of acceleration (as well as position and velocity) among them by controlling the third-order state (acceleration), while lower-order (first-order and second-order) consensus protocols are usually designed for more regular motion (e.g., rectilinear [16] and rotational [23] motion). Besides high-order dynamics [3–9], nonlinearities [1, 2, 24], time delays [8–21, 24, 25] and fuzziness [26, 27] also bring complexities to the control systems, which often lead to difficulties in stability analysis.

A novel frequency-domain-based method is developed to challenge the system complexities and derive the consensus conditions. Comparing to the universal stability analysis tool Lyapunov approaches, frequency domain methods are more possibly conducing to less conservative results as it roots in sufficient and necessary stability criteria. On the other hand, Lyapunov approaches applied in many literatures yield consensus conditions in the form of Linear Matrix Inequalities (LMIs) [7–9], while with frequency domain methods authors of [12–16, 22] as well as this note obtain consensus conditions in the form of inequalities which are more perspicuous and simple to calculate. However, frequency domain methods are limited to linear and time-invariant systems, and consequently most of the aforementioned articles [1–11, 24, 25] especially those coping with nonlinear and high-order systems have adopted Lyapunov approaches instead. Moreover, both high-order dynamics and nonuniform time delays give rise to dramatic increment of the systems’ dimensionality, which makes it a knotty problem. In existing literatures, methods based on the properties of nonnegative matrices [17, 18] and Nyquist stability theories [19–21] are introduced as alternative stability analysis tools.

The main idea of the proposed approach is to transform the high-order systems’ dynamics into high-degree polynomials with respect to hypothetically existing imaginary eigenvalues of the systems. By studying the monotonicity of the polynomials and their derivatives, consensus conditions can be figured out in the form of inequalities. The present work first brings out sufficient consensus conditions for the first-order to sixth-order nonuniformly delayed systems which are most likely to apply to practical engineering applications [28]: if all the delays are bounded by a given value and all the parameters agree to corresponding inequalities, the systems can achieve consensus and in addition, for the sake of nonconservativeness, the paper provides stronger conditions for the first-order and second-order systems by thorough derivation: the states converge when all the delays are bounded by a given value but diverge when all the delays exceed that value.

Literature [29] has proposed a high-order nonlinear consensus tracking algorithm with unmeasurable system states which applies to wide-range multiagent systems and proven the achievement of consensus by constructing Lyapunov functions, while time delays are not considered. In [9], Zhang et al. have solved average consensus problem of high-order multiagent system with time-varying delays and provided stability conditions in the form of LMIs via a Lyapunov-Krasovskii approach. The authors of [19] have derived necessary and sufficient consensus conditions for large-scale high-order linear multiagent systems with heterogeneous communication delays by using the generalized Nyquist criterion; however, the derived consensus conditions are set-valued graphical conditions, and inequality conditions are only derived for the first-order system. The previous work in [30] has studied consensus motion of delayed second-order multiagent systems by a nonconservative frequency-domain-based method, while this article will present stationary consensus conditions for more complex high-order multiagent systems.

The remainder of this note is organized as follows: Section 2 states the consensus problem with the help of graph theory; Section 3 presents the main results by demonstrating the stability analysis; Section 4 depicts the selected simulation experiments; Section 5 draws conclusions with future research directions.

#### 2. Problem Statement

This section starts with some definitions and results in graph theory.

Consider an -agent system. 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 neighbors of node is denoted by . The Laplacian corresponding to the graph is defined as , where and , . A path is a sequence of indexed edges , where . If there is a path between every pair of nodes in graph , the graph is said to be connected. The following lemma is given by [31].

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

Consider an th-order multiagent system consisting of agents. The dynamics of the th agent () iswhere is the th state variable of the th agent, , and is the control input. Let be the state vector of the th agent; we assume that the initial conditions are and , , for . The control input is said to solve the consensus problem asymptotically, if and only if for all .

In [8], a discrete-time control input was introduced as In this paper, we introduce an continuous-time consensus algorithm for system (1) with multiple time delays, and the input delays are supposed to occur. The protocol isfor any , where for ; denotes the edge weight, and is the time delay for the th agent to get the state information of the th agent. We assume that the system has different time delays, denoted by .

Let , and Under the control input given by (3), the network dynamics of the multiagent system becomeswith the initial condition , , where denotes the Laplacian of a subgraph associated with the delay . Clearly, . If all the time delays are equal to zero, system (5) could be rewritten as

This paper assumes that the graph is always connected and undirected.

#### 3. Main Results

The following lemma presents a sufficient condition for the stability of high-degree polynomials given by [32], which is helpful in the present work.

Lemma 2. *Consider a polynomial , where , , , with coefficients of determination defined as if all the coefficients of determination satisfy that , where , , and is the only real root of equationthen all the roots of have negative real parts.*

Let and suppose that the eigenvalues of are according to Lemma 1.

*Assumption 3. *Assume for that all () satisfy the following:where and is the largest eigenvalue of .

Under Assumption 3, the following lemma can be proven.

Lemma 4. *Matrix has a singleton zero eigenvalue and all other eigenvalues have negative real parts if Assumption 3 is satisfied.*

*Proof. *According to Lemma 1, there exists an orthogonal matrix , such that and then it follows thatThrough simple calculations, the eigenpolynomial of could be obtained; then we getTo simplify the following statements, let , and (13) could be written asFor the first-order system, since , it is evident that are the eigenvalues of the system; thus the lemma is proven. For the second-order system, (14) becomes and its roots areIf , the real part of (15) is , and when , (15) apparently is a pair of nonpositive real numbers, and the bigger one equals zero if and only if ; then the lemma is proven.

When , note that for , and is the largest eigenvalue of . It is apparent that ; according to Lemma 2, (10) is a group of more conservative condition, which ensures that all the roots of (14) have negative real parts except that there exists one singleton zero root for . Thus, matrix has a singleton zero eigenvalue and all its other eigenvalues have negative real parts.

*Remark 5. *Lemma 4 has shown that each of the nondelayed multiagent systems given by (6) has all its eigenvalues on the open LHP except one equals zero. That implies that system (6) is stable, and all the states of each agent will reach a common value. The existence of the only zero root indicates that only the first-order state variable of each agent reaches a value that is decided by the initial state and all the other high-order state variables return to zero at last; that is, if all , multiagent system (5) will reach consensus.

By analyzing the effect of nonuniform time delays on the stability of the systems, we will give a proof to the ensuing theorem.

Theorem 6. *Consider th-order () system given by (5) that satisfies Assumption 3, and the following inequalities (16), (17), and (18) are satisfied for , respectively:Define functions and () as follows:where is defined in (27). If all satisfy for th-order multiagent system, wherethen control input (3) can solve the consensus problem of system (5).*

*Proof. *Consider the network of high-order multiagents with nonuniform time delays. Let , where is the Laplace transform of , and 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 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 all of the first elements of the vector obtained by calculating the left part of (23) are equal to zero, we getfor all . Multiplied by (the conjugate transpose of ) on the left side of the left part of (23), and with (24) substituted, we obtainwhere , and Rewrite (25) asTake modulus of both sides of (27); then, we haveAnother necessary condition of (27) isAs , should be discussed in the positive interval . From (29), it is obvious thatDefinewhere and denote the imaginary part and real part of the complex number , respectively. is an access to .

Consider the first-order system; then, we have and . It is apparent that , and ; according to (28), we should only consider ; if we set all , then which contradicts (30). Then when all , the first-order system is impossible to have an imaginary eigenvalue which presents the first contact of the eigenvalues from the stable region to the unstable one. Hence the system is still stable then, and it can reach consensus, and the theorem is proven for .

Unlike , are not a fixed value. But we have found that the derivatives of listed below are negative values when Assumption 3 is applied: means that are monotonically decreasing with the growth of . Then it can be deduced that the arguments also decrease monotonically and continuously because the values of vary smoothly. Evidently, we havefor , which implies that the arguments start at when .

By investigating the locus of on a complex plane, we have found that the value of the argument first falls from , tending to ; as the trajectory of passes from the first quadrant to the fourth, does not really turn negative but jumps to instead and then it falls again, without second “jump” to perform. Similar phenomena occur in the other high-order systems when ; in each system, the argument performs a jump at most once and no jumps for some of these systems, such as the second-order system, because the real and imaginary parts of are always positive when . The trajectories of and on the complex plane are shown in Figure 1 to give evidence of the variation on their arguments. For those whose arguments perform jumps, we consider the value of each in two continuous intervals: and .

For , if then and , which we do not want. Thus, is the minimum to consider. Then holds when , which implies . Setting all , we have which contradicts (30); then by the same idea of the proof for the first-order system, the theorem for is proven.

Investigate the following parabola when : Apparently, its maximal value comes at the point when . According to the symmetry of parabolas, are the minima of . According to Assumption 3 that , which implies , we have Thus, . Then we have for , when .

Let be the first contact of its trajectory from the first quadrant to the fourth, as is shown in Figure 1. We analyze in two intervals: and . Obviously, is the smallest positive real root of .

Consider the interval . If , then would be smaller than the minimal possible value of to satisfy (27), because if .

Taking the fourth-order system into account, we have According to (16), one can obtainwhere . Because is monotonically increasing when , we have for . Then we know that holds if . The only situation that needs to be considered is , where is the smallest value of ; if we set all , it can be obtained thatwhich brings about the impossibility of (30) and provides the fourth-order system with consensus achievement.

For , , then Comparing and , we haveAccording to (17) and (18), , so that under Assumption 3. Thus (42) is also positive, which means . When , if we set all , for , (30) is impossible.

Consider the interval . Evidently, when , is monotonically increasing, According to Assumption 3, we have . Then Therefore, all would make (30) impossible when . As is described in Theorem 6, if all , the consensus problem of the fifth-order system is solved.

Likewise, for , according to (18) we havewhen , is a monotonically increasing function. Since is positive already, would be further greater than as the growth of has passed through where the situation should be ignored. Consequently, if all , the sixth-order system is stable, and the theorem is proven out.