International Journal of Differential Equations
Volume 2016 (2016), Article ID 9192127, 9 pages
http://dx.doi.org/10.1155/2016/9192127

Research Article

Static Consensus in Passifiable Linear Networks

Ibragim A. Junussov

Department of Control Engineering, Czech Technical University in Prague, Karlovo namesti 13, 121 35 Praha 2, Czech Republic

Received 29 December 2015; Revised 21 March 2016; Accepted 27 March 2016

Academic Editor: Jinde Cao

Copyright © 2016 Ibragim A. Junussov. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Sufficient conditions of consensus (synchronization) in networks described by digraphs and consisting of identical deterministic SIMO systems are derived. Identical and nonidentical control gains (positive arc weights) are considered. Connection between admissible digraphs and nonsmooth hypersurfaces (sufficient gain boundary) is established. Necessary and sufficient conditions for static consensus by output feedback in networks consisting of certain class of double integrators are rediscovered. Scalability for circle digraph in terms of gain magnitudes is studied. Examples and results of numerical simulations are presented.

1. Introduction

Control of multiagent systems has attracted significant interest in last decade since it has a great technical importance [1–5] and relates to biological systems [6].

In consensus problems agents communicate via decentralized controllers using relative measurements with a final goal to achieve common behaviour (synchronization) which can evolve in time. Many approaches have been developed for different problem settings.

Laplace matrix, its spectrum, and eigenspace play crucial role in description and analysis of consensus problems. Laplace matrix has broad applications, for example, [7]. Not all possible digraph topologies can provide consensus over dynamical networks. Admissible digraph topologies and connection with algebraic properties of Laplace matrix have been found in [8]. Tree structure analysis and properties of Laplace matrix spectrum of digraphs are also studied by these authors. Work [9] contains examples of out-forests as well as useful graph theoretical concepts and can be recommended as an entry reading to the research of these authors on algebraic digraph theory and consensus problems.

Concept of synchronization region in complex plane for networks consisting of linear dynamical systems is introduced in [10]. In [11] this concept is used for analysis of synchronization with leader. Problem is solved using Linear Quadratic Regulator approach in cases when full state is available for measurement and when it is not. In the last case observers are constructed.

Analysis of consensus with scalar coupling strengths [10, 11] is fruitful in a sense that conditions on gains (which depend on connection topology and single agent properties) give more insight to problem. A lot of works on topic consider dynamic couplings; however, for certain type of connections it might happen that tunable parameters will exceed upper bound on possible control gains, that is, not meet physical limitations. So, necessary and sufficient conditions on consensus achievement for different connection types in terms of coupling strengths are needed.

Celebrated Kalman-Yakubovich-Popov Lemma (Positive Real Lemma) establishes important connection between passivity (positive-realness) of transfer function and matrix relations on its minimal state-space realization ; see [12, 13]. Positive Real Lemma is a basis for Passification Method [14, 15] (“Feedback Kalman-Yakubovich Lemma”) which answers question when a linear system can be made passive, that is, strictly positive real (SPR) by static output feedback. Powerful idea of rendering system into passive by feedback has been also studied for nonlinear systems, for example, [13, 16, 17].

In consensus-type problems considering SPR agents with stable (Hurwitz) matrix leads to a synchronous behaviour when all states go to zero. The latter is undesirable in essence, since such behaviour can be reached by local control without communication. So, instead of SPR systems it is possible to consider passifiable systems, with an opportunity that a study is extendable to nonlinear systems. Also, Passification Method allows avoiding constructing observers for reaching full-state consensus by output feedback. Observers implementation increases dimension of overall phase space and complexity of a dynamic network. However, finding a passification matrix (vector ) is long standing open problem; still set of passifiable systems is nonempty.

In this paper Passification Method is used to synthesize a decentralized control law and to derive sufficient conditions of full-state synchronization by relative output feedbacks in networks described by digraphs with Linear Time Invariant dynamical nodes in continuous time. Assumptions made on network topology are minimal. Synchronous behaviour is described, including case of nonidentical gains. It is determined that boundary of sufficient gain region geometrically is a hypersurface in corresponding gain space. For certain three-node network this geometrical observation connects algebraic properties of Laplace matrix with constructed hypersurface. Namely, Jordan block appears in a direction of a cusp (nonsmooth) extremal point of the hypersurface.

Necessary and sufficient conditions for static consensus by output feedback in networks consisting of certain class of double integrators have been rediscovered. Conditions are given in terms of Laplace matrix spectrum.

Scalability in a circle digraphs in terms of gain (coupling strength) is studied. It is shown that common gain in large cycle digraphs consisting of double integrators should grow not slower than quadratically in number of agents.

Results of numerical simulations in 3- and 20-node double-integrator networks are presented.

Neighbouring papers, which influenced this work, are cited before Conclusions.

2. Theoretical Study

2.1. Preliminaries and Notations

Notations, some terms of graph theory, and Passification Lemma are listed in this section.

2.1.1. Notations

Notation stands for Euclidian norm. For two symmetric matrices , inequality means that matrix is positive definite. Notation stands for vector Identity matrix of size is denoted by Vector is vector of size and consisting of ones. Vector is defined similarly. Matrix is square matrix whose th element on main diagonal is ; other entries are zeroes. Notation stands for Kronecker product of matrices. Definition and properties of Kronecker product, including eigenvalues property, can be found in [18, 19]. Direct sum of matrices [20] is denoted by .

2.1.2. Terms of Graph Theory

A pair , where is set of vertices and is set of arcs (ordered pairs) is called digraph (directed graph). Let have elements, It is assumed hereafter that graphs do not have self-loops, that is, for any vertex, .

Digraph is called directed tree if all its vertices except one (called root) have exactly one parent. Let us agree that in any vertex is parent or neighbour. Directed spanning tree of a digraph is a directed tree formed of all digraph vertices and some of its arcs such that there exists path from any vertex to the root vertex in this tree. Existence of directed spanning tree has connection to principal achievement of synchronization in consensus-like problems.

A digraph is called weighted if to any pair of vertices number is assigned such that

A digraph in which all nonzero weights are equal to will be referred to as unit weighted.

An adjacency matrix is matrix whose th, th entry is equal to .

Laplace matrix of digraph is defined as follows:

Matrix always has zero eigenvalue with corresponding right eigenvector By construction and Gershgorin Circle Theorem all nonzero eigenvalues of have positive real parts. Let us denote by left eigenvector of which is corresponding to zero eigenvalue and scaled such that . It is known that vector describes synchronous behaviour if reached.

Suppose that a digraph has directed spanning tree. A set of digraph vertices is called Leading Set (“basic bicomponent” in terms of [9]) if subdigraph constructed of them is strongly connected and no vertex in this set has neighbours in the remaining part of digraph. Nonzero components of and only them correspond to vertices of Leading Set. Definition of basic bicomponent is wider and applicable for digraphs with no directed spanning trees.

For illustration, by [21], there are 16 different types of digraphs which can be constructed on 3 nodes. 12 of them have directed spanning tree; among these, 5 digraphs have Leading Set with 3 nodes, 2 digraphs have Leading Set with 2 nodes, and 5 digraphs have Leading Set with 1 node.

2.1.3. Passification Lemma

Problem of linear system passification is a problem of finding static linear output feedback which is making initial system passive. It was solved in [14, 15] for nonsquare SIMO and MIMO systems including case of complex parameters. Brief outline of SIMO systems passification is given below.

Let , , be real matrices of sizes , , accordingly. Denote by . Let vector . If numerator of function is Hurwitz with degree and has positive coefficients then function is called hyper-minimum-phase.

Lemma 1 (Passification Lemma [14, 15]). The following statements are equivalent.

(1)There exists vector such that function is hyper-minimum-phase.(2)Number is positive and for any there exists real matrix satisfying the following matrix relations:

2.2. Problem Statement and Assumptions

Consider a network consisting of agents modelled as linear dynamical systems:

where , is state vector, is output or measurements vector, is input or control, and , , are real matrices of appropriate size. By associating agents with vertices of unit weighted digraph and introducing set of arcs one can describe information flow in the network. For let us introduce notation for relative outputs

where is a set of th agents neighbours.

Problem is to design controllers which use relative outputs and ensure achievement of the state synchronization (consensus) of all agents:

In the case of synchronization achievement asymptotical behaviour of all agents will be described by the same time-dependant consensus vector which is denoted hereafter by :

Let us make the following assumption about dynamics of a single agent.

(A1)There exists vector such that transfer function is hyper-minimum-phase.

Now let us make an assumption on graph topology.

(A2)Digraph has at least one directed spanning tree.

Zero eigenvalue of Laplace matrix has unit multiplicity iff this assumption holds [8].

2.3. Static Identical Control

Denote

where are eigenvalues of . Under assumption (A2) zero eigenvalue is simple. By properties of other eigenvalues lie in open right half of complex plane, so is positive number.

Suppose that assumption (A1) holds with known vector . Consider the following static consensus controller with gain which is the same for all agents:

where relative output has been defined in the previous section. Denote .

Theorem 2. Let assumptions (A1) and (A2) hold. Then for all such that

controller (9) ensures achievement of goal (6) in dynamical network (4); asymptotical behaviour in the case of goal achievement is described by the following consensus vector:

Proof. Closed loop system (4) and (9) can be rewritten in the following form:

Consider nonsingular matrix (real or complex) such that

where or . All eigenvalues of have positive real parts. By considering first (zero) columns of matrices and we can accept that first column of is and first row of is .
Let us apply coordinate transformation and rewrite (12):

where , or . Note that zero solution of (15) is globally asymptotically stable iff goal (6) is achieved.
For simplicity let , , and be real till the end of proof. For any fixed satisfying (10) there exists such that

Eigenvalues of matrix have positive real parts. Therefore, according to [22], there exists real matrix such that the following Lyapunov inequality holds:

We can rewrite the last inequality:

By assumption (A1) there exists such that (3) is true with , since . Considering following Lyapunov function,

and derivating it along the nonzero trajectories of (15), we obtain

where

Matrix relations (3) have been used here. Last inequality concludes the proof.

Assumptions of this theorem are relaxed in comparison with Theorem  2 from [23]. Proof of Theorem 2 also provides the following auxiliary result.

Lemma 3. Let assumption (A2) hold. Controller (9) ensures achievement of goal (6) in dynamical network (4) if and only if all eigenvalues of matrix

have negative real parts. In the case of goal achievement asymptotical behaviour is described by (11).

2.4. Nonidentical Control and Gain Region

Let the initial digraph be unit weighted. Let us fix Laplace matrix and consider static control with nonidentical gains :

Without loss of generality we can assume that network does not have a leader (formally: cardinality of Leading Set is more than 1), since in leader case we can reduce the following consideration of synchronization gain region to lower dimension .

Let us denote by and by point which is projection of point on unit sphere

where scalar common gain is radius vector magnitude of point . Points , lie in orthant

Denote . Laplace matrices and correspond to the same digraphs which differ only in arc weights. Equation for closed loop system (4) and (23) can be rewritten as follows:

By repeating proof of Theorem 2 we can formulate the following result.

Theorem 4. Let assumptions (A1) and (A2) hold. Then for all such that

controller (23) ensures achievement of goal (6) in dynamical network (4); asymptotical behaviour in the case of goal achievement is described by the following consensus vector:

Let us denote by region in orthant such that for any point in this region control (23) ensures achievement of the goal (6) in network (4), (23). Consider the following region:

which is subset of . Let us consider closed part of unit sphere

Point on determines ray (half-line) in with initial point at the origin. According to Theorem 4, by moving along this ray from origin, that is, increasing , we will reach . Consider map

which is continuous as a composition of continuous maps ([20], continuous dependence of matrix eigenvalues on parameters). Image of this map is a subset of boundary ; therefore, by continuity of map , boundary is a hypersurface in . Further, let us consider induced map

Domain is compact, so we can apply Weierstrass Extreme Value Theorem and arrive at the following lemma.

Lemma 5. Map is continuous. Map is continuous and has minimum and maximum.

Generally, hypersurface is not smooth in all its points. Alternatively, part of a simplex of appropriate dimension can be taken instead of sphere part to serve as the domain for maps and .

Pairwise ratios of nonidentical gains and common gain define homogeneous coordinates in orthant. Common gain coefficient relates to reachability of consensus and to speed of convergence but it does not influence consensus vector. Also, consensus vector can be changed only by gain ratios variation within Leading Set of agents; see [24].

3. Double-Integrator Networks

3.1. Agents Description

Suppose that each agent in a network is modelled as follows:

For transfer function is hyper-minimum-phase. It can be shown that number .

First and second components of can describe (or can be interpreted as) velocity and position. Single system (32) can be viewed as double integrator with transfer function and proportionally differential (PD) control applied to it.

Since , static consensus controller (9) has the following form:

3.2. Necessary and Sufficient Conditions on Consensus

Let us denote by Laplace matrix of unit weighted cycle digraph which is consisting of nodes with exactly arcs:

Eigenvalues of are evenly located at circle in complex plane [25]:

Theorem 6. Controller (33) ensures achievement of goal (6) in dynamical network consisting of double integrators (32) connected in directed cycle if and only if

Proof. Let us diagonalize . Matrix from Lemma 3 in our case is block diagonal:

where

So, matrix is stable iff matrices are stable for all . Characteristic polynomial of is

Let , , . Taking in account (35) we can obtain

Now let argument of run on imaginary axis and let us decompose this polynomial on real and imaginary parts:

where and

According to Hermite-Biehler Theorem, polynomial is stable iff both of the following conditions are satisfied:

(i)roots of and are interlacing;(ii)Wronskian is positive

     for at least one value of argument

Wronskian is positive for . Root interlacing property is equivalently transformable to

Right parts of these inequalities reach maximum when (also when ) and this concludes proof.

Therefore, for a large increasing number of agents gain should grow as :

It is possible to conclude that consensus in large cycle digraphs is hard to achieve, at least for agents (32), since arbitrary high gains are not physically realizable.

On the other hand, it is worth noting that cycle digraph is the graph with minimal number of edges which is delivering average consensus among all its nodes; it is strongly connected.

Remark 7 (see [24]). Minimality in edges number causes simple relations on nonidentical gains and left eigenvector components for agents in form (4):

In other words, all pairs lie on the same hyperbola.

The following result can be obtained by repeating proof of Theorem 6.

Theorem 8. Consider network consisting of agents (32). Let a digraph, describing information flow, contain directed spanning tree. Let Laplace matrix of the digraph have real spectrum. Controller (33) ensures achievement of goal (6) in dynamical network consisting of double integrators (32) if and only if

Proof. For diagonalizable Laplace matrix with real spectrum statement is following from the well-known fact that polynomial (39) with real coefficients is stable iff its coefficients are positive. For nondiagonalizable Laplace matrix let us transform it to Jordan form. Expansion of matrix determinant shows that only determinants across main diagonal are forming (factorizing) characteristic polynomial of .

Note that undirected graphs (i.e., digraphs with symmetric ) have real spectrum and some class of digraphs have real spectrum too, for example, directed path graphs [26].

We can formulate similar result for general digraphs.

Theorem 9. Consider network consisting of agents (32). Let digraph , describing information flow, contain at least one directed spanning tree. Let all nonzero eigenvalues of Laplace matrix be denoted by . Controller (33) ensures achievement of goal (6) in dynamical network consisting of double integrators (32) if and only if

Proof. Using similar argumentation as in proofs of Theorems 8 and 6 we arrive at study of polynomial (39) stability with

Wronskian property does hold for . For root interlacing property is equivalent to trigonometric inequality

or

4. Examples and Numerical Simulations Results

4.1. Three-Node Digraph and Gain Region

Consider digraph shown in Figure 1 with dynamic nodes described in Section 3.1. By Lemma 5 distance from origin to reaches minimum. Let us draw . First, let . Let be small, and . Let be unit weighted. Eigenvalues of matrix are real: . Using Theorem 8 we conclude that . Any determines angle

and radius vector Note that pair is polar coordinates of boundary . We can conclude that minimum on is realized on a point for which and . Boundary is presented in Figure 2. For all matrix is similar to

Figure 1: Digraph of 3 nodes.

Figure 2: Gain region boundary .

Let us consider two cases:

(1), identical gains ;(2), nonidentical gains .

By Theorem 4 common gain is as follows: , . Identical gains are chosen such that .

Let us choose the following initial conditions:

Denote by sum of error norms or disagreement measure: error in the first case and error in the second case. Results of 25 sec. simulations are shown in Figure 3.

Figure 3: Performance with identical (; black line) and nonidentical (; blue line) gains.

Note that consensus vector (27) does not change for all since subsystem is leader.

4.2. Twenty-Node Digraph and Nonidentical Control

Let us consider digraph shown in Figure 4 consisting of 20 agents described in Section 3.1. This dodecahedron-like digraph has Leading Set consisting of dynamic nodes which are connected in directed circle. Let us choose and . According to Theorem 6 gain should be chosen For faster convergence let us choose . Simulations show that can be chosen considerably less than . Let us choose , and let agents have different initial conditions. Results of numerical simulation show that such nonidentical gain choice provides achievement of consensus. All trajectories of 20 agents on the same phase plane are shown in Figure 5.

Figure 4: Digraph of 20 agents.

Figure 5: Trajectories of 20 agents on the same phase plane. Trajectories of colored blue, colored red, and colored green.

Numerical simulations also show that by choosing and applying Theorem 9 for resulting Laplace matrix one can obtain lower bound approximation .

5. Reference Remarks

Assumptions of Theorem 2 are relaxed in comparison with Theorem  2 from [23]. Proof of these results uses coordinate transformation as in [27]. Lemma 3 partially succeeds Theorem  3 from [25]. Theorem 9 is a trigonometric form of Theorem  1 from [28] with different proof, which is potentially extendable to higher orders of state space.

6. Conclusions

By means of Passification Method sufficient conditions of consensus with identical and nonidentical gains are derived. Synchronous behaviour (consensus vector) is described; it can be affected by nonidentical gains (nonidentity in actuation) within Leading Set. Gain asymptote in growing cycle digraphs which have lowest communication cost for reaching average consensus and consisting of double integrators is studied.

It is rediscovered that cycle digraph connection with nonidentical actuation of nodes causes nonidentical impact on synchronous behaviour. Reachability of synchronization corresponds to positive scalar—common gain. By constructing boundary of sufficient gain region in 3-node digraph it is found that Jordan block of Laplace matrix (which affects transient dynamics) appears in a direction of extremal point. Comparison of dynamics is a subject of a future study. Geometrical interpretations which might be useful in theory and applications were developed.

Competing Interests

The author declares that he has no competing interests.

Acknowledgments

The work was supported by the Ministry of Education, Youth and Sports of the Czech Republic within Project no. CZ.1.07/2.3.00/30.0034 (support for improving R&D teams and the development of intersectoral mobility at CTU in Prague). The author would like to thank A. L. Fradkov for problem statement and Z. Hurak, A. Selivanov, and I. Herman for useful discussions.

References

F. Bullo, J. Cortés, and S. Martínez, Distributed Control of Robotic Networks, Princeton University Press, 2009.

Y. Cao, W. Yu, W. Ren, and G. Chen, “An overview of recent progress in the study of distributed multi-agent coordination,” IEEE Transactions on Industrial Informatics, vol. 9, no. 1, pp. 427–438, 2013.

M. Mesbahi and M. Egerstedt, Graph Theoretic Methods in Multiagent Networks, Princeton University Press, 2010.

R. Olfati-Saber, J. A. Fax, and R. M. Murray, “Consensus and cooperation in networked multi-agent systems,” Proceedings of the IEEE, vol. 95, no. 1, pp. 215–233, 2007.

W. Ren, R. W. Beard, and E. M. Atkins, “A survey of consensus problems in multi-agent coordination,” in Proceedings of the American Control Conference (ACC '05), pp. 1859–1864, Portland, Ore, USA, June 2005.

C. W. Reynolds, “Flocks, herds, and schools: a distributed behavioral model,” Computer Graphics, vol. 21, no. 4, pp. 25–34, 1987.

I. Gutman and D. Vidovic, “The largest eigenvalues of adjacency and Laplacian matrices, and ionization potentials of alkanes,” Indian Journal of Chemistry—Section A, vol. 41, no. 5, pp. 893–896, 2002.

R. P. Agaev and P. Y. Chebotarev, “The matrix of maximal outgoing forests of a digraph and its applications,” Automation and Remote Control, vol. 61, no. 9, pp. 15–43, 2000.

P. Chebotarev and R. Agaev, “The forest consensus theorem and asymptotic properties of coordination protocols,” in Proceedings of the 4th IFAC Workshop on Distributed Estimation and Control in Networked Systems, pp. 95–101, Koblenz, Germany, September 2013.

Z. Li, Z. Duan, G. Chen, and L. Huang, “Consensus of multiagent systems and synchronization of complex networks: a unified viewpoint,” IEEE Transactions on Circuits and Systems. I. Regular Papers, vol. 57, no. 1, pp. 213–224, 2010.

H. Zhang, F. L. Lewis, and A. Das, “Optimal design for synchronization of cooperative systems: state feedback, observer and output feedback,” IEEE Transactions on Automatic Control, vol. 56, no. 8, pp. 1948–1952, 2011.

B. Brogliato, B. Maschke, R. Lozano, and O. Egeland, Dissipative Systems Analysis and Control—Theory and Applications, Springer, New York, NY, USA, 2007.

P. Kokotović and M. Arcak, “Constructive nonlinear control: a historical perspective,” Automatica, vol. 37, no. 5, pp. 637–662, 2001.

A. L. Fradkov, “Quadratic Lyapunov functions in the adaptive stability problem of a linear dynamic target,” Siberian Mathematical Journal, vol. 17, no. 2, pp. 341–348, 1976.

A. L. Fradkov, “Passification of non-square linear systems and feedback Yakubovich-Kalman-Popov Lemma,” European Journal of Control, vol. 9, no. 6, pp. 577–586, 2003.

C. I. Byrnes, A. Isidori, and J. C. Willems, “Passivity, feedback equivalence, and the global stabilization of minimum phase nonlinear systems,” IEEE Transactions on Automatic Control, vol. 36, no. 11, pp. 1228–1240, 1991.

P. V. Kokotovic and H. J. Sussmann, “A positive real condition for global stabilization of nonlinear systems,” Systems & Control Letters, vol. 13, no. 2, pp. 125–133, 1989.

R. Bellman, Introduction to Matrix Analysis, McGraw-Hill, New York, NY, USA, 1960.

M. Marcus and H. Minc, A Survey Of Matrix Theory And Matrix Inequalities, Allyn and Bacon, Boston, Mass, USA, 1964.

R. A. Horn and C. R. Johnson, Matrix Analysis, Cambridge University Press, Cambridge, UK, 1990.

F. Harary, Graph Theory, Addison-Wesley, Reading, Mass, USA, 1969.

S. Boyd, L. Ghaoui, E. Feron, and V. Balakrishnan, Linear Matrix Inequalities in System and Control Theory, vol. 15 of Studies in Applied Mathematics, SIAM, 1994.

A. Fradkov and I. Junussov, “Synchronization of linear object networks by output feedback,” in Proceedings of the 50th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC '11), pp. 8188–8192, Orlando, Fla, USA, December 2011.

R. P. Agaev and P. Y. Chebotarev, “A cyclic representation of discrete coordination procedures,” Automation and Remote Control, vol. 73, no. 1, pp. 161–166, 2012.

J. A. Fax and R. M. Murray, “Information flow and cooperative control of vehicle formations,” IEEE Transactions on Automatic Control, vol. 49, no. 9, pp. 1465–1476, 2004.

I. Herman, D. Martinec, Z. Hurak, and M. Shebek, “Scaling in bidirectional platoons with dynamic controllers and proportional asymmetry,” IEEE Transactions on Automatic Control, http://arxiv.org/abs/1410.3943

C. Yoshioka and T. Namerikawa, “Observer-based consensus control strategy for multi-agent system with communication time delay,” in Proceedings of the 17th IEEE International Conference on Control Applications (CCA '08), pp. 1037–1042, San Antonio, Tex, USA, September 2008.

W. Yu, G. Chen, and M. Cao, “Some necessary and sufficient conditions for second-order consensus in multi-agent dynamical systems,” Automatica, vol. 46, no. 6, pp. 1089–1095, 2010.