Guaranteed Cost Tracking for Uncertain Coupled Multiagent Systems
Using Consensus over a Directed Graph^{†}^{†}thanks: This work was supported by the Australian Research Council under projects DP0987369 and DP120102152.
To be presented at the 2013 Australian Control Conference,
Perth, Australia
Abstract
This paper considers the leaderfollower control problem for a linear multiagent system with directed communication topology and linear nonidentical uncertain coupling subject to integral quadratic constraints (IQCs). A consensustype control protocol is proposed based on each agent’s states relative to its neighbors and leader’s state relative to agents which observe the leader. A sufficient condition is obtained by overbounding the cost function. Based on this sufficient condition, a computational algorithm is introduced to minimize the proposed guaranteed bound on tracking performance, which yields a suboptimal bound on the system consensus control and tracking performance. The effectiveness of the proposed method is demonstrated using a simulation example.
1 Introduction
In recent years, theoretical studies of distributed coordination and control for multiagent systems have attracted much attention in the literature, with broad applications in various areas including unmanned air vehicles (UAVs), formation control, flocking, distributed sensor networks, etc. [1]. As a result, much progress has been made in the study of cooperative control of multiagent systems [2, 3, 4].
Efforts have recently been made to consider the leaderfollowing consensus problem. For example, the leaderfollowing consensus problem for higher order multiagent systems is presented for both fixed and switching topologies in [5]. In [6], distributed observers are designed for the system of secondorder agents where an active leader to be followed moves with an unknown velocity, and the interaction topology has a switching nature. The consensusbased approach to observerbased synchronization of multiagent systems to the leader has been explored in [7, 8].
A common feature of the above literature on leaderfollowing consensusbased control problems is that interactions between agents are not considered. However, in many physical systems, interactions between agents are inevitable and must be taken into account. Examples of systems with a dynamical interaction between subsystems include power systems and spacecraft control systems [9]. This necessitates considering systems of interconnected agents.
In this paper, the leaderfollower control problem for multiagent systems coupled via linear unmodelled dynamics is considered. Coupling among the agents is regarded as an uncertainty and is described in term of time domain integral quadratic constraints (IQCs) [10]. The IQC modeling is a well established technique to describe uncertain interactions between subsystems in a large scale system [11, 12, 13].
The motivation of this paper is to extend our previous work [14] as follows. Firstly, this paper considers the multiagent system with directed topology rather than undirected topology, which poses additional difficulty compared with [14], due to the Laplacian matrix of directed graphs being in general asymmetric. Therefore a different technique is used in this paper to obtain a sufficient condition for leader follower tracking which does not involve coordinate transformation; the latter was used in [14] and required the Laplacian matrix of the graph to be symmetric. Furthermore, we consider a more general, compared to [14], class of systems with nonidentical time varying uncertain coupling. In this paper, we also propose a different LQR based cost function which describes the cost on the tracking error between the leader and all of the followers. In contrast, in [14] a consensus based cost function is considered, which penalizes the system input, the state error between the agent and its neighbours, as well as the tracking error between the leader and selected agents which observe the leader. Furthermore, the graph topology of the control protocol does not need to be the same as the topology of interconnections between the agents. Even though both communication topologies are represented as directed graphs, these graphs can be different: the agents are coupled over one directed graph, but the control protocol for the system uses another directed graph.
The main contribution of this paper is to propose a sufficient condition for the design of a guaranteed performance leaderfollower control protocol for multiagent systems with directed interconnection topology and a quite general linear uncertain coupling subject to IQCs. The sufficient condition is obtained by using a direct overbounding technique and involves checking feasibility of parameterized linear matrix inequalities (LMIs). The computational algorithm is introduced to minimize the proposed guaranteed bound by choosing local tuning parameters and guarantee a suboptimal bound on the system tracking performance.
The remainder of the paper proceeds as follows. In Section 2 of the paper, we set up the leader follower control problem for a multiagent system with directed topology and nonidentical linear uncertain coupling and give some preliminaries. The main results are given in Section III. In section IV, the computational algorithms are introduced. Section V gives an example which illustrates the theory presented in the paper. Finally, the conclusions are given in Section VI.
2 Problem Formulation and Preliminaries
2.1 Graph theory
Consider a directed graph , where is a finite nonempty node set and is an edge set of ordered pairs of nodes. The edge in the edge set of an directed graph means that the node can obtain information from node . Node is called a neighbor of node if . The set of neighbors of node is defined as . is a simple graph if it has no selfloops or repeated edges. If there is a directed path between any two nodes of the graph , then the graph is strongly connected. The adjacency matrix of the directed graph is defined as if , and otherwise. The indegree matrix is a diagonal matrix, whose diagonal elements are for . Also, let be the outdegree of node . The Laplacian matrix of the graph is defined as
2.2 Problem Formulation
Consider a system consisting of agents and a leader. All agents are assumed to be linear dynamical agents, coupled with their neighbors via, in general nonidentical, linear uncertain coupling. The connection between agents is described by a directed graph , with the node set , an edge set and a corresponding adjacency matrix . The dynamics of the th agent are described as
(1) 
where the summation is over the set of neighbors of node in the graph . The notation describes a linear uncertain operator mapping functions , into . Also, is the state, is the control input. We note that the last term in (1) reflects a relative, time varying nature of interactions between agents.
Let be the space of functions such that .
Assumption 1
All the mappings satisfy the following assumptions:

.

is linear in the second argument, i.e., if , then .

Given a matrix , there exists a sequence , such that for every , the following IQC holds
(2)
The sequence is assumed to be the same for all . The class of operators that satisfy these assumptions will be denoted by . We note that matrices are assumed to be fixed.
In addition to the system (1), suppose a leader is given. The dynamics of the leader, labeled , is expressed as
(3) 
where is its state. The control communication topology between agents is described by a directed graph , with the same node set , but possibly different edge set and a corresponding adjacency matrix . The Laplacian matrix of the graph is denoted as . We assume throughout the paper that the leader node can be observed from a subset of nodes of the graph . If the leader is observed by the node , we extend the graph by adding the edge with weighting gain , otherwise let . We refer to node with as a pinned or controlled node. The diagonal matrix is commonly referred to as the pinning matrix. The system is assumed to have at least one agent connected to the leader, hence .
Define error vectors as , . Then dynamics of the synchronization errors satisfy the equation
(4) 
In this paper we are concerned with finding a control protocol for each node of the form
(5) 
where is the feedback gain matrix to be found, and is the set of neighbors of node in the graph . As a measure of system performance, we will use the quadratic cost function,
(6) 
where and are given weighting matrices.
Remark 1
In [14] we considered a different cost function,
Each addend in this cost function penalizes the th system input, the disagreement between the th and the th system states, where is a neighbor of , as well as the tracking errors between the leader and the pinned agents which observe the leader. In contrast, the cost function (6) in this paper describes the cost on the tracking error between the leader and all of the followers and system input.
Taking linearity of the operator into account, the synchronization error dynamic (4) can be represented as
(7) 
The problem in this paper is to find a control protocol (5) which solves the leader following consensus control problem as follows:
3 The Main Result
In this section, the main result of this paper is presented which is a sufficient condition for the system (1) to be able to track the leader with guaranteed tracking performance.
First we present the following result of [15] and some notation.
Assumption 2
The digraph contains a spanning tree and the root node obtains information from the leader node, i.e., .
Let be the maximum eigenvalue of and . Also let . According to Lemma 1, is a positive definite and symmetric matrix. Let be an orthogonal matrix such that
(10) 
and denote . For node of the graph , introduce matrices , , where are the elements of the neighbourhood set , and are the nodes with the property ; and are the indegree and the outdegree of node , respectively, in the graph . Also, introduce the matrix .
Theorem 1
Proof: Using the Schur complement, the LMIs (11) can be transformed into the following Riccati inequality
(13) 
After pre and postmultiplying (3) by and multiplying (3) by , then substituting into the Riccati inequality (3), we obtain
(14) 
Define and consider the following Lyapunov function candidate for the subsystems (7):
(15) 
Then
(16) 
Note the following inequality:
(17) 
where .
(19) 
Using the following identity,
one has
(20) 
The expression on the right hand side of the above inequality is independent of . Letting leads to
(23) 
(24) 
4 The Computational Algorithm
In this section, we provide an algorithm to calculate a suboptimal control gain . According to Theorem 1, the upper bound on tracking performance is given by the right hand side of (12). Hence, one can achieve a suboptimal guaranteed performance by optimizing this upper bound over the feasibility set of the LMIs (11):
(26) 
where the infimum is taken over the feasibility set of the LMIs (11), .
As in [12], the optimization problem (26) can be shown to be equivalent to the minimization of subject to the constraints
(27) 
By the Schur complement, (27) is equivalent to the LMI
(28) 
where
This leads us to introduce the following optimization problem in the variables and : Find
(29) 
where the infimum is with respect to and subject to (11) and (28).
Theorem 2
.
Proof: The proof of this theorem is similar to the proof of Theorem 15 in [12].
Based on the foregoing discussion, we propose an algorithm for the design of the suboptimal protocol (5) based on Theorems 1 and 2:
It must be noted that the above algorithms require the knowledge of initial conditions of the leader and agents. In practice, however, the initial state of the leader may not be known. It is possible to avoid using in these algorithms using the approaches outlined in [12].
5 Example
To illustrate the proposed method, consider a system consisting of three identical pendulums coupled by two springdamper systems. Each pendulum is subject to an input as shown in Fig. 1. The dynamics of the coupled system are governed by the following equations
(30)  
where is the length of the pendulum, and are the positions of the springdamper, is the gravitational acceleration constant, is the mass of each pendulum, and are the spring constant and damping coefficient for the leftmost springdamper pair, while and are the spring constant and damping coefficient for the rightmost springdamper pair. The position of the springdamper can change along the full length of the pendulums and is considered to be uncertain, that is , .
In addition to the three pendulums, consider the leader pendulum which is identical to those given. Its dynamics are described by the equation
(31) 
Choosing the state vectors as , , and , equations (5) and (31) can be written in the form of (1), (3) , where , , and , , , ,
The agents in this example are coupled according to the undirected graph shown in Fig. 2, which is treated here as a special case of directed graph with symmetric adjacency matrix. On the other hand, the control communication topology of the system is assumed to be a linear directed graph shown in Fig. 3. According to this graph, only agent observes the leader. The Laplacian matrix of the graph consisting of nodes 1, 2 and 3 and the pinning matrix are
To illustrate the design based on Theorems 1 and 2, the LMI problem in Theorem 2 was solved numerically, and then the trajectories of the coupled pendulum system with the obtained protocol were simulated. To this end, the parameters of the coupled pendulum system were chosen to be , , , , , , , , . In the cost function, we let and . Using the computational algorithm based on Theorem 2, the problem (26) was found to be feasible and yielded the gain matrix . The performance bound was minimized by , with parameters , , , , , , , . The simulation results for this protocol are shown in Fig. 4. Also, using the controller obtained by means of the computational algorithm proposed in Theorem 2, we directly computed the performance cost (6) for the system to be , while the theoretically predicted bound is .
6 Conclusions
The consensus control for leadertracking problem with guaranteed tracking performance for nonidentical uncertain coupled linear systems connected over a directed graph has been discussed in this paper. A sufficient condition was proposed by using the direct overbounding of the performance cost. According to the simulation results, the proposed computational algorithm based on Theorems 1, 2, which solve N coupled LMIs, guarantees a suboptimal performance.
References
 [1] R. OlfatiSaber and R. M. Murray, Consensus problems in networks of agents with switching topology and timedelays, IEEE Trans. Autom. Contr., vol. 49 (9), pp. 15201533, Sep. 2004.
 [2] R. OlfatiSaber, J. A. Fax, and R. M. Murray, Consensus and cooperation in networked multiagent systems, Proc. IEEE, vol. 95 (1), pp. 215233, Jan. 2007.
 [3] W. Ren, R. W. Beard, and E. M. Atkins, Information consensus in multivehicle cooperative control: Collective group behavior through local interaction, IEEE Control Syst. Mag., vol. 27 (2), pp. 7182, Apr. 2007.
 [4] A. Arenas, A. DiazGuilera, J. Kurths, Y.Moreno, and C. Zhou, Synchronization in complex networks, Phys. Rep., vol. 468 (3), pp. 93153, 2008.
 [5] W. Ni and D. Cheng, Leaderfollowing consensus of multiagent systems under fixed and switching inteconnection topologies, Syst. & Contr. Letters, 59 (34):209217, 2010.
 [6] Y. Hong, G. Chen, and L.Bushnell, Distributed observers design for leaderfollowing control of multiagent networks, Automatica, 44 (3):846850, 2008.
 [7] V. Ugrinovskii. Distributed robust filtering with