arxiv: v1 [cs.sy] 4 Nov 2015

Transcription

1 Novel Distributed Robust Adaptive Consensus Protocols for Linear Multi-agent Systems with Directed Graphs and ExternalDisturbances arxiv:5.033v [cs.sy] 4 Nov 05 YuezuLv a,zhongkuili a,zhishengduan a,gangfeng b a State Key Laboratory for Turbulence and Complex Systems, Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing 0087, China Abstract b Department of Mechanical and Biomedical Engineering, City University of Hong Kong, Hong Kong, China This paper addresses the distributed consensus protocol design problem for linear multi-agent systems with directed graphs and external unmatched disturbances. A novel distributed adaptive consensus protocol is proposed to achieve leader-follower consensus for any directed graph containing a directed spanning tree with the leader as the root node. It is noted that the adaptive protocol might suffer from a problem of undesirable parameter drift phenomenon when bounded external disturbances exist. To deal with this issue, a distributed robust adaptive consensus protocol is designed to guarantee the ultimate boundedness of both the consensus error and the adaptive coupling weights in the presence of external disturbances. Both adaptive protocols are fully distributed, relying on only the agent dynamics and the relative states of neighboring agents. Key words: Multi-agent systems, cooperative control, consensus, distributed control, adaptive control, robustness. Introduction In recent years, the consensus problem of multi-agent systems has been an emerging research topic in the field of control, due to its wide applications in many areas such as satellite formation flying, cooperative unmanned systems, and distributed reconfigurable sensor networks []. There has been remarkable progress in achieving consensus for different scenarios; see[,,3,4,5,6] and the references therein. For the consensus problem, the critical task is to design distributed consensus protocols based on local information, i.e., local state or output information of each agent and its neighbors. In this paper, we consider the consensus problem of multi-agent systems with general linear time in- This work was supported by the National Natural Science Foundation of China under grants 60453, 3300, 6503, a Foundation for the Author of National Excellent Doctoral Dissertation of PR China, and the State Key Laboratory of Complex System Intelligent Control and Decision. Corresponding author Zhongkui Li. addresses: yzlv@pku.edu.cn (Yuezu Lv), zhongkli@pku.edu.cn (Zhongkui Li), duanzs@pku.edu.cn (Zhisheng Duan), megfeng@cityu.edu.hk (Gang Feng). variant dynamics. Previous works [7,8,9,0,,] have presented various static and dynamic consensus protocols, which are proposed in a distributed fashion, using only the local information of each agent and its neighbors. However, those consensus protocols involves some design issues. To be specific, the design of those consensus protocols generally requires the knowledge of some eigenvalue information of the Laplacian matrix associated with the communication graph, that is, the smallest nonzero eigvenvalue of the Laplacian matrix for undirected graphs and the smallest real part of the nonzero eigenvalues of the Laplacian matrix for directed graphs. Note that the nonzero eigenvalue information of the Laplacian matrix is global information in the sense that each agent has to know the entire communication graph to compute it. Therefore, although these consensus protocols are proposed and can be implemented in a distributed fashion, they cannot be designed by each agent in a distributed fashion. In other words, those consensus protocols in [7,8,9,0,,] are not fully distributed. To remove the limitation of requiring global information of the communication graph, distributed adaptive consensus protocols are reported in[3,4,5,6], which, Preprint submitted to Automatica 5 November 05

2 depending on only local information of each agent and its neighbors, are fully distributed. It is worth noting that the adaptive protocols in [3,4,5,6] are applicable to only undirected communication graphs or leaderfollower graphs where the subgraphs among followers are undirected. Due to the asymmetry of the Laplacian matrices, it is much more difficult to design distributed adaptive consensus protocols for general directed communication graphs. By including monotonically increasing functions to provide additional design flexibility, a distributed adaptive consensus protocol is derived in[7] to achieve consensus for general leader-follower directed graphs containing a directed spanning tree. The robustness of the distributed adaptive protocols with respect to uncertainties or external disturbances is an important issue. The adaptive protocol in[7] can only be modified to be applicable to external disturbances satisfying the matching condition; see [9]. To the best of our knowledge, how to design distributed robust adaptive consensus protocols for the case with directed graphs and general unmatched disturbances is still open. In this paper, we aim to design distributed robust adaptive consensus protocols for linear multi-agent systems with directed communication graphs. A novel distributed adaptive protocol is presented and shown to achieve leader-follower consensus for directed communication graphs containing a directed spanning tree with leader as the root node. This novel adaptive protocol is fully distributed, relying on only the agent dynamics and the relative state information of neighboring agent. In the presence of external disturbances, it is pointed out the adaptive protocol might suffer from a problem of the parameter drift phenomenon [8]. Therefore, the adaptive protocol is not robust in the presence of external disturbances. To deal with this instability issue, a robust adaptive protocol is presented, which can guarantee the ultimate boundedness of both the consensus error and the adaptive coupling weights. The existence condition of the proposed adaptive protocols are also discussed. Compared to the previous works [7] and [9], the contribution of this paper is at least two-fold. First, the adaptive protocol proposed in this paper replaces the multiplicative functions in the adaptive protocol in [7] by novel additive functions. In this case, a simple quadratic-like Lyapunov function, rather than the complicated integral-like one in [7], can be used to derive the result. Second, in contrast to the adaptive protocol in [9] which works only for the case with disturbances satisfying the restrictive matching condition, the robust adaptive consensus protocol given in this paper is applicable to the case of general bounded disturbances. It should be mentioned that the methods used to derive the results in this paper arequite different fromthose in [7] and [9]. The rest of this paper is organized as follows. The mathematical preliminaries used in this paper is summarized in Section. The distributed adaptive consensus protocol is designed in Section 3 for general directed leader-follower graphs. A novel robust adaptive consensus protocol is presented in Section 4 to deal with external disturbances. Simulation results are presented in Section 5. Section 6 concludes this paper. Mathematical Preliminaries Let R n m be the set of n m real matrices and the superscript T donates transpose for real matrices. I N represents the identity matrix of dimension N and I denotes the identity matrix of an appropriate dimension. donates a column vector with all entries equal to. diag(a,,a N ) represents a diagonal matrix with elements a i, i =,,N, on its diagonal while λ min (A) donates the minimal eigenvalue of A. The matrix inequality A > B means A and B are symmetric matrices and A B is positive definite. A B represents the Kronecker product of matrices A and B. A nonsingular M-matrix A = [a ij ] means that a ij < 0, i j, and all eigenvalues of A have positive real parts. A directed graph G consists of a node set V and an edgesete V V,inwhichanedgeisrepresentedbyan orderedpair of distinct nodes. For an edge (v i,v j ), node v i is called the parent node, node v j the child node, and v i is a neighbor of v j. A path from node v i to node v il is a sequence of ordered edges of the form (v ik, v ik+ ), k =,,l. A directed graph contains a directed spanning tree if there exists a node called the root such that the node has directed paths to all other nodes in the graph. Suppose there are N nodes in the directed graph G. TheadjacencymatrixA = [a ij ] R N N ofg is defined bya ij = if (v i,v j ) E and 0 otherwise.the Laplacian matrix L = [l ij ] R N N is defined as l ii = N j= a ij and l ij = a ij, i j. Lemma ([0]) Zero is an eigenvalue of L with as a right eigenvector and all nonzero eigenvalues have positive real parts. Furthermore, zero is a simple eigenvalue of L if and only if G has a directed spanning tree. Lemma ([,]) Consider a nonsingular M- matrix L. There exists a diagonal matrix G so that G diag(g,,g N ) > 0 and GL+L T G > 0. Lemma 3 ([3]) If a and b are nonnegative real numbers and p and q are positive real numbers such that p + q =, then ab ap p + bq q, and the equality holds if and only if a p = b q. Lemma 4 ([4]) For a system ẋ = f(x,t), where f( ) islocallylipschitzinxandpiecewisecontinuousint,assume that there exists a continuously differentiable func-

3 tion V(x,t) such that along any trajectory of the system, α ( x ) V(x,t) α ( x ), V(x,t) α 3 ( x )+ǫ, where ǫ > 0 is a constant, α and α are class K functions, and α 3 is a class K function. Then, the solution x(t) of ẋ = f(x,t) is uniformly ultimately bounded. 3 Distributed Adaptive Consensus Protocol Design ConsideragroupofN+identicalagentswithgeneral linear dynamics, consisting of N followers and a leader. The dynamics of the i-th agent are described by ẋ i = Ax i +Bu i, i = 0,,N, () wherex i R n is the state, u i R p is the controlinput, A and B are constant matrices with compatible dimensions. Without loss of generality, let the agent in() indexed by 0 be the leader whose control input is assumed to be zero, i.e., u 0 = 0, and the other agents be the followers. The communication graph G among the N + agents is assumed to satisfy the following assumption. Assumption The graph G contains a directed spanning tree with the leader as the root node. Under Assumption, the Laplacian matrix [ L associated with G can be partitioned as L =. ] 0 0 N L L In light of Lemma and the definition of M-matrix, it is easy to verify that L R N N is a nonsingular M- matrix. The objective of this paper is to design distributed consensus protocols such that the N agents in() achieve leader-followerconsensusinthesenseoflim t x i (t) x 0 (t) = 0, i =,,N. Based on the relative states of neighboring agents, we propose a distributed adaptive consensus protocol to each follower as u i = (c i +ρ i )Kξ i, ċ i = ξ T i Γξ i, i =,,N, () where ξ i N j=0 a ij(x i x j ), c i (t) denotes the timevarying coupling weight associated with the i-th follower with c i (0) 0, K R p n and Γ R n n are the feedback gain matrices, and ρ i are smooth functions to be determined. Let ξ [ξ T,,ξT N ]T and x [x T,,xT N ]T. Then, we get ξ = (L I N )(x x 0 ). (3) Since the graph G satisfies Assumption, it follows from Lemma that L is a nonsingular M-matrix and that the leader-follower consensus problem is solved if and only if ξ asymptotically converges to zero. Hereafter, we refer to ξ as the consensus error. Substituting () into () yields ξ = [I N A+L (C +ρ) BK]ξ, ċ i = ξ T i Γξ i, i =,,N, where C diag(c,,c N ) and ρ diag(ρ,,ρ N ). (4) The following theorem provides a result on the design of the adaptive consensus protocol (). Theorem Suppose that the communication graph G satisfies Assumption. Then, the leader-follower consensus problem of the agents in () can be solved under the adaptive protocol () with K = B T P, Γ = P BB T P, and ρ i = ξ T i P ξ i, where P > 0 is a solution to the following linear matrix inequality (LMI): PA T +AP BB T < 0. (5) Moreover, each coupling weight c i converges to some finite steady-state value. Proof Consider the Lyapunov function candidate: V = g i(c i +ρ i )ρ i + λ 0 c i, (6) where G diag(g,,g N ) is a positive definite matrix such that GL + L T G > 0, λ 0 denotes the smallest eigenvalue of GL +L T G, and c i c i α, where α is a positive constant to be determined later. It follows fromassumption and Lemma that L is a nonsingular M-matrix. Thus we know from Lemma that such a positive definite matrix G does exist. Since c i (0) > 0, it follows from ċ i 0 that c i (t) > 0 for any t > 0. Then, it is not difficult to see that V is positive definite. The time derivative of V along the trajectory of (4) 3

4 is given by V = [g i (c i +ρ i )ξi T P ξi +g i ρ i ċ i ] N +λ 0 (c i α)ċ i = ξ T [(C +ρ)g (P A+A T P ) (C +ρ)(gl +L T G)(C +ρ) Γ]ξ +ξ T (ρg Γ)ξ +ξ T [λ 0 (C αi) Γ]ξ ξ T [(C +ρ)g (P A+A T P ) λ 0 (C +ρ) Γ]ξ +ξ T (ρg Γ)ξ +ξ T [λ 0 (C αi) Γ]ξ. By using Lemma 3, we can get that and ξ T (ρg Γ)ξ ξ T ( λ 0 ρ Γ)ξ +ξ T ( λ 0 G Γ)ξ, (7) (8) ξ T (λ 0 C Γ)ξ ξ T ( λ 0 C Γ)ξ +ξ T ( λ 0 I Γ)ξ. (9) Substituting (8) and (9) into (7) yields V ξ T [(C +ρ)g (P A+A T P ) ( λ 0 (C +ρ) +λ 0 αi λ 0 λ 0 G ) Γ]ξ. (0) Choose α ˆα+ λ0 +max,,n g i, where ˆα > 0 will λ 0 be determined later. Then, it follows from (0) that V ξ T [(C +ρ)g (P A+A T P ) ( λ 0 (C +ρ) +λ 0ˆαI) Γ]ξ ξ T [(C +ρ)g (P A+A T P ) ˆαλ 0 (C +ρ) Γ]ξ. () Let ξ = ( (C +ρ)g P )ξ and choose ˆα to be sufficiently large such that ˆαλ 0 G I. Then we can get from () that V ξ T (I N (AP +PA T BB T )) ξ 0, () where the last inequality comes directly from the LMI (5). Therefore, we can get that V (t) is bounded and so is each c i. Noting that ċ i 0, we can know that each couplingweightc i convergestosomefinitevalue.noting that V 0 is equivalent to ξ 0 and thereby ξ 0. By LaSalle s Invariance principle[5], it follows that the consensus error ξ asymptotically converges to zero. That is, the consensus problem is solved. Remark Contrary to the consensus protocols in [7,8,0,] which use the nonzero eigenvalues of the Laplacian matrix, the design of the proposed adaptive protocol () relies on only the agent dynamics and the relative states of neighbors, which can be conducted by each agent in a fully distributed way. As shown in [7], a necessary and sufficient condition for the existence of the solution P > 0 to the LMI (5) is that (A,B) is stabilizable. Therefore, the existence condition of an adaptive protocol () satisfying Theorem is that (A, B) is stablizable. Remark In contrast to the distributed adaptive protocols in [3,4,5,6] which are applicable to only undirected graphs, the proposed adaptive protocol () works for the case with general directed graphs satisfying Assumption. It is worth mentioning that similar distributed adaptive protocols were designed in the previous works [7] and [9] for directed graphs satisfying Assumption. In comparison to the adaptive protocols in [7] and [9], the novel adaptive protocol () has two distinct features. First, different from the adaptive protocol in [7] which uses multiplicative functions to provide additional design flexibility, the adaptive protocol () introduces additive functions ρ i instead. As a consequence, a simple quadratic-like Lyapunov function as in (6), instead of the complicated integral-like Lyapunov function in [7], can be used to show Theorem. Second, contrary to the adaptive protocol in [9] which can only deal with external disturbances satisfying the restrictive matching condition, the proposed adaptive protocol () can be modified to be applicable to general bounded external disturbances, which will be detailed in the following section. 4 Distributed Robust Adaptive Consensus Protocols Theorem in the previous section shows that the adaptive protocol() is applicable to any directed graph satisfying Assumption for the case without external disturbances. In many circumstances, the agents might be subject to various external disturbances, for which case it is necessary and interesting to investigate whether the adaptive protocol () is robust. The dynamics of the i-th agent are described by ẋ i = Ax i +Bu i +ω i, i = 0,,N, (3) whereω i R n denotesexternaldisturbancesassociated with the i-th agent, which satisfy the following assumption. Assumption There exist positive constants υ i such that ω i υ i, i =,,N, and Bu 0 +ω 0 v 0. 4

5 Note that due to the existence of disturbances ω i in (3), the relative states will not converge to zero any more but rather can only be expected to converge into some small neighborhood of the origin. Since the derivatives of the adaptive gains c i in () are of nonnegative quadratic forms in terms of the relative states, in this caseitiseasytoseefrom()thatc i willkeepgrowingto infinity, which is called the parameter drift phenomenon in the classic adaptive control literature [8]. Therefore, the adaptive protocol () is not robust in the presence of external disturbances. In the following, we aim to make modification on () to propose a distributed robust adaptive protocol which can guarantee the ultimate boundedness of the consensus error and adaptive weights for the agents in (3). We propose a new robust distributed adaptive consensus protocol as follows: u i = (d i +ρ i )Kξ i, d i = ϕ i (d i ) +ξ T i Γξ i, i =,,N, (4) where d i (t) denotes the time-varying coupling weight associated with the i-th follower with d i (0), ϕ i, i =,,N, are small positive constants, and the rest of the variables are defined as in (). Substituting (4) into (3), it follows that ξ = [I N A+L (D +ρ) BK]ξ +(L I n )ω d i = ϕ i (d i ) +ξi T Γξ (5) i, i =,,N, where C diag(c,,c N ), ω [ω T (Bu 0 + ω 0 ) T,,ω T N (Bu 0 + ω 0 ) T ] T, and the rest of the variables are defined as in (4). In light of Assumption, we have that ω N (v i +v 0 ). (6) Note that d i (0) and d i 0 when d i = in (5). Then,itisnotdifficulttoseethatd i (t) foranyt > 0. The following theorem presents a result on design of the robust adaptive consensus protocol (4). Theorem Suppose that Assumptions and hold. Then, both the consensus error ξ and the coupling weights d i, i =,,N, in (5) are uniformly ultimately bounded under the adaptive protocol (4) with K = B T Q, Γ = Q BB T Q, and ρ i = ξ T i Q ξ i, where Q > 0 is a solution to the LMI: AQ+QA T +εq BB T < 0, (7) where ε >. The upper bound of the consensus error ξ will be given in the proof. Proof Consider the Lyapunov function candidate: V = g i(d i +ρ i )ρ i + λ 0 d i, (8) where d i d i α, where α is a positive constant to be determined later, and the rest of the variables are defined as in (6). Since g i > 0, d i (t) for any t > 0, and ρ i 0, it can be similarly shown as in the proof of Theorem that V is positive definite. By following similar steps in deriving Theorem, we can obtain the time derivative of V along (5) as V ξ T [(D +ρ)g (Q A+A T Q Γ)]ξ +ξ T [(D +ρ)gl Q ]ω ξ T [ϕ(d I) G Q ]ξ λ 0 ϕ i (d i α)(d i ), (9) where α is chosen to be sufficiently large as in the proof of Theorem and ϕ diag(ϕ,...,ϕ N ). By choosing α > and using Lemma 3, we can get that (d i α)(d i ) = (d i ) 3 +(α )(d i ) = (d i ) 3 +[( 3 4 ) 3 (di ) ][( 4 3 ) 3 (α )] (d i ) 3 + (d i ) (α )3 = (d i ) (α )3. Note that ξ T [(D +ρ)gl Q ]ω (0) = ξ T [(D I) ϕg Q ]( ϕ GL Q )ω +ξ T ( G Q )( GL Q )ω +ξ T ( ρg Q )( ρgl Q )ω ξ T [(D I) ϕg Q ]ξ + ( ϕ GL Q )ω + ξt (G Q )ξ + ( GL Q )ω + ξt (ρg Q )ξ + ( ρgl Q )ω, () 5

6 where we have used Lemma 3 several times to get the last inequality, and ( ρgl Q )ω ρ G I n (G 4 L Q )ω ρ G I n 4 + (G 4 L Q )ω 4 ξt (ρg Q )ξ + (G 4 L Q )ω 4, () where we have used matrix norm properties to get the first inequality, and Lemma 3 to get the second inequality, and to get the last inequality we have used the fact that ρ G I n 4 = max ( ρ i gi ) 4,,N ( ρ i gi ) 4 = ξ T (ρg Q )ξ. Substituting (0), (), and () into (9) yields V ξ T [(D +ρ)g (Q A+A T Q Γ)]ξ +ξ T [(ρ+ λ 0 )G Q ]ξ ϕ i(d i ) 3 +Π W(ξ) ξt (G Q )ξ where we have used the fact that D I, Π = and λ 0 ϕ i(d i ) 3 +Π, (3) 6λ 0 7 ϕ i(α ) 3 + ( ϕ GL Q )ω + ( GL Q )ω + (G 4 L Q )ω 4. W(ξ) ξ T [(D +ρ)g (Q A+A T Q +Q Γ)]ξ 0. Note that for any positive δ, we have the following assertion: δλ 0 (d i α) δλ 0 (d i ) + δλ 0 (α ) = [( 3λ 0ϕ i 4 ) 3 (di ) ][( λ 0 δλ 9ϕ ) 0 3 δ]+ i (d i α) λ 0 ϕ i(d i ) 3 + λ 0δ 3 7ϕ + δλ 0 i (α ), (4) where we have used the fact that d i > and α > to get the first inequality, and Lemma 3 to get the last inequality. From (3) and (4), we can obtain that where and V δv +δv +W(ξ) ξt (G Q )ξ λ 0 ϕ i(d i ) 3 +Π δv + W(ξ) ξt (G Q )ξ +Π, Π = [ λ 0δ 3 7ϕ i + δλ 0 (α ) ]+Π, W = ξ T [(D +ρ)g (Q A+A T Q +(+δ)q Q BB T Q )]ξ, (5) By choosing δ such that ε +δ, we can obtain that W(ξ) 0. Then, it follows from (5) that V δv ξt (G Q )ξ +Π. (6) In light of Lemma 4, we can conclude from (6) that both the consensus errorξ and the adaptivegains d i are uniformly ultimately bounded. Further, from (6), we can get that V δv if ξ Therefore, ξ converges to the set D = ξ : Π ξ λ min (Q ) min Π λ min(q ) min g i.,,n,,n g i with a convergence rate faster than e δt. (7) Remark 3 As shown in Proposition in [8], there exists a Q > 0 satisfying (7) if and only if (A,B) is controllable. Thus, a sufficient condition for the existence of (4) satisfying Theorem is that (A, B) is controllable, which, compared to the existence condition of () satisfying Theorem, is more stringent. Theorem shows that the modified adaptive protocol (4) does ensure the ultimate boundedness of both the consensus error ξ and the adaptive gains d i. That is, the adaptive protocol (4) is robust in the presence of external disturbances. The upper bound of the consensus error ξ as given in (7) depends on the dynamics of each agent, the communication graph, the upper bounds of the disturbances, and the parameters ϕ i. We should choose appropriately small ϕ i to get an acceptable upper bound of ξ. 6

7 Remark 4 Compared to the robust adaptive protocol in [9] which are only applicable to the case with matching disturbances, the adaptive protocol(4) works for general external disturbances. This is a favorable consequence of introducing novel additive functions ρ i, rather than multiplicative ones as in [9], into (4). It is worth noting that the procedures in showing Theorem is quite different from those in [9]. xi x Simulation Consider a network of second-order integrators, described by (), with x i = [ xi x i ] [ ] [ ] 0 0, A =, B = t (s) Fig.. The consensus errors x i x 0, i =,, ci.4 The communication graph is given as in Fig., which clearly satisfies Assumption t (s) Fig. 3. The adaptive coupling weights c i. 3 5 Fig.. The leader-follower directed communication graph. Solving the LMI (5) by using the LMI toolbox [ of Matlab gives ] a feasible solution matrix P = Then, the feedback gain matrices of () are given by K = [ ] , Γ = 6 4 [ ] Let c i (0) =, i =,,6. The consensus errors x i x 0, i =,,6, of the second-order integrators under the adaptive protocol () are depicted in Fig. and the adaptive coupling weights c i are shown in Fig. 3. Further, consider the case where the second-order integrators are perturbed by external disturbances. For illustration, the disturbances associated [ with ] 0.sin(t) the agents are assumed to be ω 0 =, 0.3 sin(4t) [ ] [ ] 0.sin(3.5t) 0.5cos(4t) ω =, ω =, ω 3 = 0.3 cos(.5t) 0. sin(5t) [ ] [ ] [ 0.3sin(x3 ) 0.3e t 0.sin(4t),ω 4 =,ω 5 = 0.6 sin(3t) 0.5 cos(3t) 0.5 cos(3t) ω 6 = 0.3sin(5t), and the control input of the leader 0.4 x 6 + is assumed to be u 0 = e 0.t. Solving the LMI (7) [ ] with ε = gives Q = and then [ ] [ ] K = , Γ =. In (4), let ϕ i = 0. and d i (0) =.5, i =,,6. The consensus errors x i x 0, i =,,6, under the robust adaptive protocol (4) are depicted in Fig.4 and the coupling weights d i are shown in Fig. 5, both of which are obviously bounded. ], 7

8 xi x Fig. 4. The consensus errors under the robust adaptive protocol (4). di t (s) t (s) Fig. 5. The coupling weights d i in the presence of disturbances. 6 Conclusion In this paper, we have presented novel distributed adaptive consensus protocols for linear multi-agent systems with external disturbances and directed graphs containing a directed spanning tree with the leader as the root. The adaptive consensus protocols, depending on only the agent dynamics an the relative state information of neighboring agents, can be designed and implemented in a fully distributed way. One contribution of this paper is that the new distributed adaptive protocol is robust in the presence of general bounded external disturbances. An interesting topic for future investigation is to design fully distributed adaptive protocols for nonlinear multi-agent systems or the case with local output information of each agent and its neighbors. References [] W. Ren, R. Beard, and E. Atkins (007). Information consensus in multivehicle cooperative control. IEEE Control Systems Magazine vol. 7 no. (7 8). [] Y. Hong, G. Chen, and L. Bushnell (008). Distributed observers design for leader-following control of multi-agent networks. Automatica vol. 44 no. 3 ( ). [3] R. Olfati-Saber and R. Murray (004). Consensus problems in networks of agents with switching topology and timedelays. IEEE Transactions on Automatic Control vol. 49 no. 9 (50 533). [4] T. Li, M. Fu, L. Xie, and J. Zhang (0). Distributed consensus with limited communication data rate. IEEE Transactions on Automatic Control vol. 56 no. (79 9). [5] Z. Li, Z. Duan, and G. Chen (0). On H and H performance regions of multi-agent systems. Automatica vol. 47 no. 4 ( ). [6] M. Guo and D. V. Dimarogonas (03). Consensus with quantized relative state measurements. Automatica vol. 49 no. 8 (53 537). [7] Z. Li, Z. Duan, G. Chen, and L. Huang (00). Consensus of multiagent systems and synchronization of complex networks: A unified viewpoint. IEEE Transactions on Circuits and Systems I: Regular Papers vol. 57 no. (3 4). [8] Z. Li, Z. Duan, and G. Chen (0). Dynamic consensus of linear multi-agent systems. IET Control Theory and Applications vol. 5 no. (9 8). [9] S. Tuna (009). Conditions for synchronizability in arrays of coupled linear systems. IEEE Transactions on Automatic Control vol. 54 no. 0 (46 40). [0] J. Seo, H. Shim, and J. Back (009). Consensus of high-order linear systems using dynamic output feedback compensator: Low gain approach. Automatica vol. 45 no. ( ). [] H. Zhang, F. Lewis, and A. Das (0). Optimal design for synchronization of cooperative systems: State feedback, observer, and output feedback. IEEE Transactions on Automatic Control vol. 56 no. 8 (948 95). [] C. Ma and J. Zhang (00). Necessary and sufficient conditions for consensusability of linear multi-sgent systems. IEEE Transactions on Automatic Control vol. 55 no. 5(63 68). [3] Z. Li, W. Ren, X. Liu, and L. Xie (03). Distributed consensus of linear multi-agent systems with adaptive dynamic protocols. Automatica vol. 49 no. 7 ( ). [4] Z. Li, W. Ren, X. Liu, and M. Fu (03). Consensus of multi-agent systems with general linear and Lipschitz nonlinear dynamics using distributed adaptive protocols. IEEE Transactions on Automatic Control vol. 58 no. 7(786 79). [5] H. Su, G. Chen, X. Wang, and Z. Lin (0). Adaptive second-order consensus of networked mobile agents with nonlinear dynamics. Automatica vol. 47 no. ( ). [6] W. Yu, W. Ren, W. X. Zheng, G. Chen, and J. Lv (03). Distributed control gains design for consensus in multi-agent systems with second-order nonlinear dynamics. Automatica vol. 49 no. 7 (07 5). [7] Z. Li, G. Wen, Z. Duan, and W. Ren (04). Designing fully distributed consensus protocols for linear multi-agent systems with directed communication graphs. IEEE Transactions on Automatic Control, contionally accepted. [8] P. A. Ioannou and J. Sun (996). Robust Adaptive Control. New York NY: Prentice-Hall, Inc.. [9] Z. Li, Z. Duan (04). Distributed robust adaptive consensus protocols for linear multi-agent systems with directed graphs and external disturbances. The 04 Chinese Control Conference, in press. [0] W. Ren and R. Beard (005). Consensus seeking in multiagent systems under dynamically changing interaction topologies. IEEE Transactions on Automatic Control vol. 50 no. 5 (655 66). 8

9 [] Z. Qu (009). Cooperative Control of Dynamical Systems: Applications to Autonomous Vehicles. London, UK: Spinger- Verlag. [] H.Zhang, F.L.Lewis,and Z.Qu(0). Lyapunov, adaptive, and optimal design techniques for cooperative systems on directed communication graphs. IEEE Transactions on Industrial Electronics vol. 59 no. 7 ( ). [3] D. S. Bernstein (009). Matrix Mathematics: Theory, Facts, and Formulas. Princeton University Press. [4] M. Corless and G. Leitmann (98). Continuous state feedback guaranteeing uniform ultimate boundedness for uncertain dynamic systems. IEEE Transactions on Automatic Control vol. 6 no. 5 (39 44). [5] M. Krstic, I. Kanellakopoulos, and P. Kokotovic (995). Nonlinear and Adaptive Control Design. New York: John Wiley Sons. 9