Consensus of the Multiagent System with a Dynamic Leader Based on Directed Topology Using Laplace Transform

The consensus of themultiagent systemwith 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.


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 .
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][8][9][10][11][12][13][14][15][16][17][18][19][20][21][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.

Model and Preliminaries
2.1.Model.Consider a multiagent system of  + 1 agents, where the leader is labeled with 0 and the followers are labeled with 1, 2, . . ., .The motion of the system is described by where   ∈ R  represents the state of agent ; N  is the neighbor set of agent . = [  ] ∈ R × is the coupling matrix, where   is the weight parameter with   ≥ 0 and   = 0.   ( 0 −   ) is the term performing attractive effect from leader agent 0 to agent .Define  = diag{ 1 ,  2 , . . .,   } ∈ R × .  > 0 means agent  can measure the information of the leader with strength   , while   = 0 means there is no information from the leader to agent .V() is the velocity of the leader, which is a bounded function.

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 G = (V, E, ) be a weighted digraph of order  with the set of nodes V = {1, 2, 3, . . ., } and set of arcs E ⊆ V × V.  = [  ] is adjacency matrix of graph G.An arc of G is denoted by (, ), which starts from  and ends at .The set of neighbors of node  is denoted by N  = { ∈ V : (, ) ∈ E}.A digraph G 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 G to node , we say that node  is globally reachable.
A diagonal matrix  = diag{ 1 ,  2 , . . .,   } ∈ R × is a degree matrix of G with its diagonal elements   = ∑ ∈N    ,  = 1, 2, . . ., .Then the Laplacian of the weighted digraph G(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 G, if there exists at least one path from every node in G to the leader, we say that the leader is globally reachable in G.

Lemma 2. The nonzero eigenvalues of 𝐿 have positive real parts.
Proof.The proof is similar to that of [3]; it is omitted here.
Notice that the block matrix of  is Then, we can have rank( + ) = , where  is the Laplacian matrix of G.
Theorem 4. If leader 0 is globally reachable in digraph G, 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 L(⋅) in this paper.System (5) can be transformed as Define − 0 =  as the error vector; then (6) can be written as Applying L to (7), then we have That is, For simplicity, denote L() = Y() and L() = V(), 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 (0) − V(); and Δ[  + ( + )] is the determinant of matrix   + ( + ).We can separate (11) into two parts as follows: where Δ[  + ( + )]  1 is the determinant of matrix   + ( + ) in which the th column has been replaced by (0); Δ[  + ( + )]  2 is the determinant of matrix   + ( + ) in which the th column has been replaced by V().Because 2 can be transformed as where 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) where Applying Inverse Laplace transform to (17), we obtain According to integral property, it is easy to have Since as  → +∞, we have  (−) → 0, then we can get the following inequality: 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 0 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.

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.

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.

15 Figure 1 :
Figure 1: Errors between leader and followers about system (1), in which 0 is globally reachable in G and V() = 1.3.

Figure 6 :Figure 7 :
Figure 6: Trajectories of the agents about system (1), in which 0 is globally reachable in G.