#### Abstract

The consensus of the multiagent system with directed topology and a leader is investigated, in which the leader is dynamic. Based on Laplace transform method, the accurate upper error bound between the leader and the followers can be obtained. It is also proved that all agents of the system will aggregate and eventually form a cohesive cluster following the leader if the leader is globally reachable. Finally, some simulation examples are given to illustrate the theoretical results.

#### 1. Introduction

In recent years, the consensus problem of multiagent systems has become a hot topic due to their broad applications, such as cooperative unmanned air vehicles, automated highway systems, air traffic control, and autonomous underwater vehicles [1–22].

Various algorithms and models about multiagent systems have been discussed based on different tasks or interests. The leader-following system is one of the most interesting topics in the motion control of the multiagent systems. Vicsek et al. [1] proposed a simple model about autonomous agents moving with a constant identical speed and tending to the average direction of its neighbors. It is demonstrated numerically that all the agents will move in the same direction at the same speed eventually. Jadbabaie et al. [2] gave the theoretical explanation for the numerical results of the Vicsek model by algebraic graph theory. Based on the Vicsek model, some researchers designed the leader-following model, where the leader agent can be regarded as the control input used to control the other agents in such system. Shi et al. [5] regarded the reference signal as a virtual leader for guiding the agent group to move at the desired constant velocity. Hong et al. [7] discussed the leader-following system with variable coupling topology. Olfati-Saber [3] used virtual leaders to accomplish obstacle avoidance.

The coupling topology plays an important role in the studies of multiagent systems. Because of the complexity in the consensus analysis with directed topology, most researchers focus on the undirected topology or balance topology to simplify the problem [7–22]. In this paper, we will discuss the consensus problem with directed topology and give a weak condition for reaching network consensus. According to Laplace transform, the accurate upper error bound between the leader and followers can be calculated and some assumptions needed in Lyapunov approach can be simplified. Finally, some simulation examples are provided to illustrate the theoretical results.

This paper is organized as follows. Section 2 proposes the model formulation and some preliminaries of graph theory. Section 3 gives the analysis of the convergence of the model. Section 4 gives simulations of the theoretical results. Finally, the main contribution of this paper is summarized in Section 5.

#### 2. Model and Preliminaries

##### 2.1. Model

Consider a multiagent system of agents, where the leader is labeled with and the followers are labeled with . The motion of the system is described by where represents the state of agent ; is the neighbor set of agent . is the coupling matrix, where is the weight parameter with and . is the term performing attractive effect from leader agent to agent . Define . means agent can measure the information of the leader with strength , while means there is no information from the leader to agent . is the velocity of the leader, which is a bounded function.

##### 2.2. Preliminaries

To discuss the coordinated control among the agents, graph theory is a very effective tool. Regard the agent as a node and the connection link between any two agents as an edge; the coupling topology is conveniently described by a directed graph. Let be a weighted digraph of order with the set of nodes and set of arcs . is adjacency matrix of graph . An arc of is denoted by , which starts from and ends at . 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 with its diagonal elements , . Then the Laplacian of the weighted digraph (or matrix ) is defined as In order to study the leader-following problem, a multiagent system with directed topology is considered, which is consisting of agents and one leader. In , if there exists at least one path from every node in to the leader, we say that the leader is globally reachable in .

It is easy to see that has a zero eigenvalue corresponding to the right eigenvector** 1 **.

Lemma 1 (see [9]). *A digraph has a globally reachable node if and only if 0 is a simple eigenvalue of (i.e., ).*

Lemma 2. *The nonzero eigenvalues of have positive real parts.*

*Proof. *The proof is similar to that of [3]; it is omitted here.

Lemma 3. *The eigenvalues of have positive real parts if and only if the leader is a globally reachable node in digraph .*

*Proof. *If node is globally reachable, the Laplacian matrix of can be written as
According to Lemma 1, then .

Notice that the block matrix of is
Then, we can have , where is the Laplacian matrix of .

#### 3. Coordinated Control Analysis

In this section, we focus on the coordinated control problem of model (1). Denote , , and ** 1**; then the model can be rewritten into matrix form
Based on system (5), we can derive the following theorem.

Theorem 4. *If leader is globally reachable in digraph , all the agents described by (1) will converge and form a cohesive cluster following the leader asymptotically. Moreover, the errors between the leader and followers will be included in a fixed bound.*

*Proof. *To solve this problem, we introduce the Laplace transform, succinctly denoted by in this paper. System (5) can be transformed as
Define as the error vector; then (6) can be written as
Applying to (7), then we have
That is,
For simplicity, denote and , where is the Laplace variable; then
According to* Cramer rule*, we can get solutions of (10):
where is the identity matrix of order ; is the determinant of matrix in which the th column has been replaced by ; and is the determinant of matrix . We can separate (11) into two parts as follows:
where is the determinant of matrix in which the th column has been replaced by ; is the determinant of matrix in which the th column has been replaced by . Because **1**, determinant can be transformed as
where is the determinant of matrix in which the th column has been replaced by** 1**. Since , , and are polynomials about , for simplicity, we denote , , and , respectively. Then (12) can be rewritten as
According to Lemma 3, the eigenvalues of have positive real parts, denoted as
with multiplicity as . So we have

If the fractional expressions in (14) are reducible, we should make reduction of the fraction expressions. In the following discussion, we assume that the fractional expressions in (14) are irreducible. Applying* Heaviside’s Method*, (14) can be expanded into
where
Applying* Inverse Laplace transform* to (17), we obtain
According to integral property, it is easy to have
Since as , we have , then we can get the following inequality:
as where is a function close to , is the upper bound of , and . It shows that the upper bound of is influenced by and the coupling matrix . Define and choose ; we have
as . Then, any agent of the system will enter into a region with the leader as the center bounded by within sufficient time.

This completes the proof.

*Remark 1. *In Theorem 4, the error between the leader and any agent is bounded. If the velocity of leader is zero, we can get that the error will tend to zero, which implies that all the agents will asymptotically approach the same states of the leader. If the velocity of the leader is nonzero and bounded, it is easy to see that the states and velocities of all the followers will keep in a bounded region following the leader.

#### 4. Numerical Simulations

In order to verify the above theoretical analysis, we present some numerical simulations to illustrate the systems. These simulations are performed with ten followers and one leader, and the initial positions of the agents are chosen randomly. The coupling matrix is given based on certain conditions.

Figures 1, 2, and 3 show the errors between leader and followers about system (1), under different coupling topology and control. It is easy to see that the states of the error are contained in a bounded region. Figures 4 and 5 depict the trajectories of the agents in 2-dimensional Euclidian space. If leader agent in the coupling topology is globally reachable, then the followers in the system will reach synchronization with the leader within sufficient time. Figures 6 and 7 show the trajectories of the agents in 3-dimensional Euclidian space, in which the states of the followers are influenced by the state of leader, and all agents in the system can achieve consensus eventually.

#### 5. Conclusion

In this paper, we have investigated the coordinated control of leader-following multiagent systems. It is proved that the agents of the systems will aggregate and form a cluster following the leader asymptotically. Meanwhile, we have studied the coupling topology among the agents in the general case. The systems considered in this paper can better reflect the collective behavior in practice. It is clear that the ideas and approaches about graph theory and linear dynamical system theory will play an important role in the analysis of the multiagent system. Finally, simulations give an effective demonstration of the leader-following systems.

#### 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 no. 61304049, the Science and Technology Development Plan Project of Beijing Education Commission (no. KM201310009011), and the Plan training project of excellent young teacher of North China University of Technology.