Abstract

A formation control problem for second-order multiagent systems with time-varying delays is considered. First, a leader-following consensus protocol is proposed for theoretical preparation. With the help of Lyapunov-Krasovskii functional, a sufficient condition under this protocol is derived for stability of the multiagent systems. Then, the protocol is extended to the formation control based on a multiple leaders’ architecture. It is shown that the agents will attain the expected formation. Finally, some simulations are provided to demonstrate the effectiveness of our theoretical results.

1. Introduction

Recent years have witnessed a rapidly growing interest in coordinated control of multiagent systems due to its broad applications in various disciplines [112]. As one of the important topics in this field, formation control has attracted great attention. Generally speaking, the main objective in formation control is to design appropriate protocol and algorithm such that the agents can achieve and preserve a predefined geometrical shape, such as a chain or a wedge. Potential applications of formation control include lots of cooperative tasks such as surveillance, exploration, search and rescue, transporting large objects, and control of arrays of satellites.

In the literature, researchers have proposed numerous approaches for the multiagent systems to achieve the anticipated formation, roughly categorized as behavior-based strategy [13], virtual structure method [14], and leader-following approach [1521], to name a few. For instance, Xiao et al. [16] developed a formation framework of multiple leaders and applied a class of nonlinear consensus protocols to the formation control. Under the proposed framework, all the agents in the first-order systems could reach the expected formation. An adaptive formation control approach, in the absence of the velocity information of the leader, was proposed in [17]. Besides, the authors in [18] investigated the leader-following formation control problems for nonlinear systems under fixed and switching topologies. The above works, however, did not take into account the effects of time delays.

Owing to the finite speed of information transmission and processing, time delays are inevitable in multiagent systems. In particular, one type of time delays is communication delays, whose effects on multiagent systems have been addressed by many researchers [39]. For formation control with time delays, Luo et al. [19] gave a sufficient condition of formation control of multiagent systems by using Lyapunov stability theory. Also, Rezaee and Abdollahi [20] provided a motion synchronization strategy with time delays. Note that the time delays in both papers were assumed to be constant. In reality, it is more practical when the time-varying delays are accounted for. Lu et al. [21] studied the formation control of second-order multiagent systems with time-varying delays, where the time delays existed only in the transmission of position information between neighbors. For second-order multiagent systems, it is worthwhile to mention that the information exchanged between neighbors may include velocity information as well as position information.

Motivated by the above analysis, we consider a leader-following formation control problem for second-order multiagent systems, with time-varying delays existing in the transmission of both velocity and position. Here, we adopt the formation framework proposed in [16]. More specifically, the formation information is divided into two parts: the global one and the local one, where the former determines the geometric pattern of the desired formation and the latter decides the relative information of agents with respect to their neighbors. Only a small number of agents called leaders have access to the global information, and the other agents called followers regulate their states according to the local information.

An outline of this paper is as follows. Section 2 provides some preliminary notions of graph theory and formulates the formation control problem. By utilizing Lyapunov-Krasovskii functional, Section 3 presents the main results under the proposed control protocol. Then several simulations are illustrated in Section 4. Finally, the conclusion is drawn in Section 5.

The following notations will be used throughout this paper. Given a matrix, the superscripts “” and “−1” stand for its transposition and inverse, respectively; and denote the set of all eigenvalues and the spectral norm of the matrix, respectively. Let be an identity matrix, , and represents an zero matrix. And denotes a diagonal matrix with diagonal elements being . For a complex number , , , and are its real part, imaginary part, and modulus, respectively. denotes the Kronecker product.

2. Preliminaries and Problem Formulation

Let be a weighted directed graph with the set of nodes , set of edges , and a weighted adjacency matrix with nonnegative elements . The node indexes belong to a finite index set . An edge in is denoted by , where if and only if . In this case, we say node is a neighbor of node and denote the neighbors of node by . Moreover, we assume for all . Let the Laplacian matrix associated with be defined as and .

A directed path in directed graph from to is a sequence of edges of the form , where for . A directed graph is called strongly connected if any two distinct nodes of the graph can be connected via a directed path that follows the edges of the graph. A directed tree is a directed graph, where every node, except one special node without any parent, which is called the root, has exactly one parent, and the root can be connected to any other nodes through paths. A spanning tree of a digraph is a directed tree formed by graph edges that connect all the nodes of the graph.

For the -agent system considered in this paper, suppose that there are    leaders and followers. For convenience, we use and to denote, respectively, the leader set and the follower set. Further, the interaction topology among agents is modeled by a direct graph .

The dynamics of the autonomous agents are given by where , , denote the position, velocity and control input of agent , respectively.

Using the formation framework in [16], the formation information is divided into two independent parts: the global one and the local one. The global information, which determines the geometric pattern of the expected formation, is represented by a time-dependent column vector with   . The local formation information is denoted by a time-independent nonnegative matrix with unit entry sum for each row    .

Denote and ; we have the following.

Definition 1. We say that system (1) solves the formation problem if there exists a -valued function with , , and such that , for and , for as . In particular, the formation problem is called a time-invariant formation (TIF) problem, a time-varying formation (TVF) problem, and a time-varying formation for trajectory tracking (TVFT) problem if , with and with , respectively.

The column vector , which specifies the state of the formation and may be dependent on initial states or may be an external input, used to guide the group of agents to track a prescribed trajectory. The vector defines the basic frame of the anticipated formation formed by the leaders, and the nonnegative matrix specifies the local-state restrictions of followers with respect to their leader neighbors. Since each row entry sum of equals 1, the followers should lie in the convex region covered by the leaders.

3. Main Results

3.1. Leader-Following Consensus

For the better understanding of the formation control problem, we first consider a leader-following consensus problem. The multiagent systems solve a leader-following consensus asymptotically if and for all . Under this circumstance, all the agents can obtain the acceleration , but only the root agent can obtain the difference between its state and the formation information, and hence it can pin the other leaders to attain the anticipated formation. Define a matrix as , where if the interaction topology has a spanning tree rooted at , and otherwise.

The time-varying delay is the time delay for information communicated between agents at time . Owing to the communication time delays, each agent cannot instantly get the information from others. With time-varying delays existing in the transmission of both velocity and position, we now provide the following consensus protocol: where the control parameter , is the adjacency matrix corresponding to the graph , if the interaction topology has a spanning tree rooted at and otherwise, and the time-varying delay is a continuously differentiable function with

Denote , ; then, system (1) with protocol (2) can be written in a matrix form: where is the Laplacian matrix associated with , , if the interaction topology has a spanning tree rooted at and otherwise.

Let , ; then, (4) can be rewritten as

Taking , (5) can be equally expressed as where To proceed further, we need the following lemmas.

Lemma 2 (see [22]). Given a complex-coefficient polynomial where , is Hurwitz stable if and only if and .

Lemma 3. The matrix is Hurwitz stable if and only if is positive stable and .

Proof. Let be an eigenvalue of . Then one has Clearly, the Hurwitz stability of matrix is equivalent to that of the polynomial: , where is the th eigenvalue of . It follows from Lemma 2 that is Hurwitz stable if and only if and , .

Lemma 4 (see [4]). The matrix is positive stable if has a spanning tree.

Now we give the main result of this subsection.

Theorem 5. For system (6), take and suppose that Then, the consensus is reached asymptotically if the graph has a spanning tree.

Proof. Since has a spanning tree, is positive stable based on Lemma 4. It follows from Lemma 3 that is Hurwitz stable. Therefore, there exists a positive definite matrix such that To analyze the convergence of system (6), we define a Lyapunov-Krasovskii functional where and .
Calculating along the solution of (6), we have
Due to the fact that and holds for any appropriate positive definite matrix , we can obtain that Similarly, Substituting (15) and (16) into (14) leads to Consequently, a sufficient condition for is From (18), we can obtain that Therefore, by Lyapunov-Krasovskii Theorem (see [23]), the error system (6) is uniformly asymptotically stable. Namely, the consensus is reached asymptotically.

Remark 6. It can be seen that many zoom techniques are applied during the derivation of , which result in a conservative estimation of .

3.2. Time-Varying Formation for Trajectory Tracking

In this subsection, we consider the case that the agents form a time-varying formation as they track the desired trajectory. More specifically, the desired trajectory of the formation, represented by , is assumed to be determined by the following equation: Also make the following assumptions about the multiagent system:(A1)the local interaction topology of the leaders has a spanning tree, and the leaders’ dynamic is unaffected by the followers;(A2)the root agent in the local interaction topology of the leaders is able to access the reference trajectory;(A3)in addition to the local information from their neighbors, the leaders can also obtain the global information ;(A4)in addition to the local information from their follower neighbors, each follower can also get the local formation directly or indirectly from the leaders.

Using the following control protocol: where , is the adjacency matrix corresponding to the graph , and if there is a spanning tree rooted at in the graph and otherwise, we have the following result.

Theorem 7. Suppose that the graph has a spanning tree. Take With protocol (21), the multiagent systems attain a time-varying formation for trajectory tracking (TVFT) under assumptions (A1)–(A4).

Proof. Let , for , and , for . Then we can rewrite protocol (21) as Denote and ; (23) can be expressed in a matrix form: It follows from Theorem 5 that and converge to zero asymptotically, and equally, , for and , for as .

Remark 8. Obviously, Theorem 7 still holds if , or and . In the form case, the multiagent systems attain a time-invariant formation (TIF). In the latter case, the multiagent systems attain a time-varying formation.

4. Simulations

In this section, to illustrate our theoretical results derived in the above section, we will provide several examples. Consider a multiagent system consisting of 25 agents moving in a plane (i.e., ), with the direct interaction topology described in Figure 1. The expected formation is a hexagram. Assume that the first 7 agents are leaders with the spanning tree rooted at agent 7. For simplicity, let if agent is a neighbor of agent and otherwise. With simple calculations, we can obtain that and . Take , , and .

Consider

Figure 2 shows the trajectories of some agents in time-invariant formation (TIF), where the initial states of the agents are randomly generated in a given bounded region. From Figure 2, we can see that the agents attain the vertexes of a hexagram. In other words, the expected time-invariant formation is reached. For the trajectory tracking formation, Figure 3 illustrates the formation states of agents. Since the expected trajectory of the formation satisfies the equation in (25), the formation in Figure 3 changes with time as it moves along a sinusoidal curve. Furthermore, the position and velocity errors of agents in TVFT with time delays existing only in the transmission of position are shown in Figures 4(a) and 5(a), respectively, and the position and velocity errors of agents in TVFT with time delays existing in the transmission of both position and velocity are shown in Figures 4(b) and 5(b), respectively. It can be seen that all of the errors converge to zero ultimately, while the errors in the latter figures converge to zero faster than those in the former figures.

5. Conclusion

This paper investigates the formation control problem for second-order multiagent systems with time-varying delays. We first consider a leader-following consensus problem. By employing Lyapunov-Krasovskii functional, we prove that the multiagent systems can reach consensus. Then, under a special multiple leaders’ framework, we apply the protocol to the formation control, and derive a sufficient condition for the system to achieve prescribed formation. Moreover, several numerical simulations are shown to verify the theoretical analysis.

Conflict of Interests

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

Acknowledgments

The authors would like to thank the anonymous reviewers and the Editor Su for their insightful comments and suggestions. This research is supported by the Chinese Universities Specialized Research Fund for the Doctoral Program (20110185110020), NSFC (61170309, 61104141, and 61301052), the Key Technology R&D Program of Sichuan Province (2012SZ0159), and 973 Program (2013CB329404).