Abstract

This paper investigates the consensus problem of multiagent systems with directed topologies. Different from the literatures, a new method, the Laplace transform, to study the consensus of multiagent systems with directed topology and communication time delay is proposed. The accurate state of the consensus center and the upper bound of the communication delay to make the agents reach consensus are given. It is proved that all the agents could aggregate and eventually form a cohesive cluster in finite time under certain conditions, and the consensus center is only determined by the initial states and the communication configuration among the agents. Finally, simulations are given to illustrate the theoretical results.

1. Introduction

In recent years, there has been an increasing interest in the study of consensus and rendezvous problems in the multiagent systems [116]. This is partly due to the wide applications in cooperative control of unmanned air vehicles, formation control of mobile robots, design of sensor networks, flocking of ants and birds, distributed decision making, and so on [510]. Vicsek et al. [1] proposed a discrete model of autonomous agents. Each agent of such model was moving at a constant identical velocity, and the direction was updated via a local rule based on the average of the directions of its neighbors. At the same time, they gave some numerical simulations to describe the dynamic behavior of the model. Jadbabaie et al. [3] gave a theoretical explanation for the numerical results of Vicsek’s model. In [5], Olfati-Saber and Murray analyzed the consensus problem of multiagent systems with time delay and obtained the accurate bound of the time delay with undirected topology. In [7, 9], the authors discussed the consensus of second-order multiagent systems with time delay by Lyapunov approach. However, due to certain limitations, the bound of delay in those papers is not specific.

The analysis of the coupling topology plays an important part in discussing the consensus problems. Jadbabaie et al. [3] discussed the communication information by applying an undirected graph to model the coupling topology among the agents. Olfati-Saber and Murray [5] investigated the average consensus problem with directed topology, where the coupling topology is undirected and needs to satisfy the balance condition, which is a strong condition for the consensus. Ren and Beard [11] introduced the definition of spanning tree to depict the coupling topology. It was shown that consensus can be achieved asymptotically if the directed interaction graph contains a spanning tree as the system evolves. A useful Lemma about the Laplacian matrix was given in [11], and Laplacian matrix has a simple zero eigenvalue if and only if the directed graph has a spanning tree. Lin et al. [12] introduced the definition of globally reachable node to describe the coupling topology and gave a similar lemma about the Laplacian matrix; that is, the digraph has a globally reachable node if and only if 0 is a simple eigenvalue of Laplacian matrix.

In this paper, we discuss the consensus problem of the multiagent system mentioned in [5] with directed topology and time delay. Compared with the previous references, the main contribution of this paper is to study the coupling topology in a more general case. According to the Laplace transform, we can not only have the specific consensus center of the model but also obtain the accurate upper bound of the communication delay value to achieve consensus. By the Laplace transform, we can eliminate some assumptions used in the Lyapunov approach and give further simplification of the assumptions.

This paper is organized as follows. Section 2 presents some preliminaries of graph theory. Section 3 proposes the model and gives the analysis for the consensus of the model. In Section 3, some extensions of the studies are discussed. Section 4 gives simulations to verify the theoretical results. Finally, we summarize our main contribution in Section 5.

2. Some Preliminaries

To discuss the coupling topology of the communication configuration of the agents, graph theory is a very effective tool. If each agent is regarded as a node, then the coupling topology is conveniently described by a directed graph. In particular, in the definition of the directed graph, self-loops are excluded. Let be a weighted digraph of order with the set of nodes and set of arcs , and is the adjacency matrix of graph . An arc of is denoted by , which is from to . The set of neighbors of node is denoted by . A digraph is strongly connected if there exists a path between any two distinct nodes. For a node , if there exists at least a path from every other node in to node , we say that node is globally reachable.

A diagonal matrix is a degree matrix of , and its diagonal elements for . The Laplacian of the weighted digraph (or matrix ) is denoted as It is obvious that the digraph , the adjacency matrix , and the Laplacian matrix are peer-to-peer. Some basic properties of the Laplacian matrix need to be introduced in the following.

Lemma 1 (see [12]). The digraph has a globally reachable node if and only if the Laplacian matrix has a simple zero eigenvalue (with eigenvector ).

Lemma 2 (see [5]). The nonzero eigenvalues of are of positive real part.

Lemma 3. For arbitrary row of the Laplacian matrix , the cofactor of any elements is equal.

Proof. For Laplacian without loss of generality, the first row of is chosen to discuss. The cofactors of element and , , are, respectively, denoted by For , adding all the other columns to the th column, we can have Since is arbitrary, we can obtain . Similarly, we can have .
This completes the proof.

3. Model Formulation and Analysis

Consider a multiagent system consisting of agents in -dimensional Euclidian space whose motion is governed by the following delay differential equations: where represents the state of agent , is the neighbor set of agent , is the coupling weight matrix, and is weight parameter with , if agent has information with agent ; otherwise, ; the time delay is a constant.

For simplicity of discussion, we take in model (5). Then, model (5) can be rewritten into a matrix form: where and is the Laplacian matrix.

For system (6), we can have the following main results by Laplace transform.

Theorem 4. For system (6), if there exists at least one globally reachable node in and the time delay parameter satisfies then the states of all agents of the system will asymptotically converge to a constant value; that is, where are the nonzero eigenvalues of and and are the argument and modulus of , respectively.

Proof. Applying Laplace transform on (6), succinctly denoted as , we can obtain that is, where is the Laplace variable. For simplicity, is denoted as .
According to Cramer rule, we can get the solutions of (10): where is the determinant of in which the th column has been replaced by ; is the determinant of matrix .
By calculating, we can obtain That is fractional expression about variable . If there are common factors between and , we can reduce them and make them irreducible. In the following, assume that and are irreducible. For simplicity, we denote and .
According to Lemma 1, the zero eigenvalue of is simple and the eigenvalues of are denoted as where are the multiple eigenvalues corresponding with the multiplicity as , respectively. Then, it follows that By Heaviside’s method, (12) can be expanded into where Equation (15) is the transfer function of ; the Laplace reverse transform of the first term is a constant value . In order to make asymptotically stable, the real part of the denominators’ roots must be negative, which implies that the roots of transcendental equation must have negative real part. Then, the roots of will be analyzed in the following. To solve (17), we let , and , where , and is the plural unit. Then, which can be expanded as that is, Separating the real and imaginary parts, then we can have Solving the above equations, we can get where is the argument of and .
According to (22), we can know that
Substituting those into the first equation of (22), we have that is, Let be equal to the critical value; that is, ; then In order to get the smallest positive value of , we have , and , since and the sign of is decided by and . So, we can get the critical value of the time delay, denoted as , and In the following, we will show that the roots of (17) have negative real part when .
Since the complex eigenvalues of are conjugated, their imaginary signs arguments are opposite. According to (25), we define the function as For arbitrary , let , and if , then For function , there exists a positive root , and it decreases monotonically in the interval .
Taking there exists only one solution in the interval , denoted as .
Then From (29) and (31), there exists a root of between and . Moreover, if , the same result holds.
For arbitrary , then Because is a monotonically decreasing function in the interval , and , which implies that they have no positive roots of (28) but have infinite negative roots of (28).
Define for arbitrary ; the roots of the denominator of each term of (15) must have negative real part.
Take arbitrary term of (15), simply denoted as where is the coefficient, is the degree, and . It has infinite singular points which have negative real part. The Laplace inverse transform of the function asymptotically tends to ; that is, for the function, when and . So, we can obtain when ; that is, as .
According to Lemma 1, we can get This completes the proof.

Remark 5. If matrix is symmetrical, then the imaginary parts of the eigenvalues are equal to zero; that is, . The critical value of the time delay is .

Furthermore, the multiagent system is described by We can also have the same result as follows.

Theorem 6. For model (38), if there exists at least one globally reachable node in digraph and the time delay parameter satisfies then all the agents of the system can achieve asymptotically consensus.

In model (5), we regarded the feedback gain value as 1. If the feedback gain value is , then the model can be described as which can be rewritten as the matrix form It can be easy to get that the consensus center is not changed; then, we have the following theorem.

Theorem 7. For model (40), if there exists at least one reserve tree in the digraph associated with matrix and the time delay parameter satisfies then all the agents of the system will asymptotically converge upon a fixed point , where

Remark 8. The feedback gain can be used to adjust the consensus velocity of the agents and the critical value of delay.

4. Numerical Simulations

In this section, we will give some numerical simulations to illustrate the theoretical results.

Consider a multiagent system with ten agents, where the initial states of agents are chosen randomly. The coupling matrices are given as follows, respectively:

Matrix is the associated Laplacian matrix with only one zero eigenvalue, and the other eigenvalues are all equal to 1. And the critical value of delay is .

For matrix , the associated digraph is strongly connected, and the eigenvalues of are;; .

The critical value of delay .

Figures 1 and 2 show the simulation results, where the coupling matrix is . We can find that, from Figure 1, the states of all agents of model (5) will asymptotically converge to a constant value when (), respectively; however, from Figure 2, the states of all agents of the system will diverge when (), respectively. Figure 3 shows, for model (38), the states of all agents of the system with the coupling matrix , and the feedback gain will asymptotically converge to a constant value when (), respectively.

5. Conclusion

In this paper, we have considered the consensus of multiagent systems with directed topology and communication time delay. We have proved that the system aggregates and forms a cluster in finite time if the time delays are smaller than the critical value. The methods and results of this paper can be extended to discuss the leader-follower second-order multiagent system with time delay.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant (nos. 61304049, 61174116, and 51308005) and Science and Technology Development Plan Project of Beijing Education Commission (no. KM201310009011).