共识问题即一致性问题,是多自主体网络的一个最基本的分布式协调控制问题.在过去几十年里,多自主体的共识问题在许多领域都引起了极大关注[1-5].其中,出现了一类线性切换多自主体系统一致性问题,包括无领导者一致性问题[6]和领导者跟随一致性问题[7-8].广义系统有动态系统的自然表示,比一般的线性系统有更广泛的背景[9-10].文献[11]和[12]分别利用状态反馈和输出反馈来设计控制协议,给出了广义多自主体系统达到一致的充分必要条件.但都只考虑了固定拓扑情形下的一致性问题,对于更一般的动态拓扑(切换拓扑)没有研究,给出的条件虽然是充分必要条件,但是证明过程却是从系统达到一致性的条件出发.
本文将文献[8]中的一般线性系统推广到广义系统,并讨论在切换拓扑下的一类广义多自主体系统的一致性问题.与文献[12]相比,这里讨论的拓扑图是动态图, 而且分两种情形(无领导者和领导者跟随)来研究其一致性.通过代数图论[13]和广义系统理论[14]得到结论:要解决广义多自主体系统的两个一致性问题,只需要相应的慢子系统达到一致性.
1 预备知识对于给定的向量或矩阵X,‖ X ‖表示X的欧几里得模.向量1N表示所有元素都是1的列向量,span{X}表示由X的列向量张成的线性子空间,|a|表示不超过实数a的最大的整数,A1/2表示正定矩阵A的二次方根.⊗代表Kronecker积,满足如下性质:(A⊗B)T=AT⊗BT, (A⊗B)(C⊗D)=(AC)⊗ (BD), (A+B)⊗ C=A⊗C+B⊗C, A⊗(B+C)=A⊗B+A⊗C.
通常,一个多自主体系统中每个自主体之间的信息交换可以通过有向图或无向图来描述[13].
2 主要结果本文将考虑式如
$ \mathit{\boldsymbol{E}}{{\dot x}_i} = \mathit{\boldsymbol{A}}{\mathit{\boldsymbol{x}}_i} + \mathit{\boldsymbol{B}}{\mathit{\boldsymbol{u}}_i} $ | (1) |
广义多自主体系统的稳定性,其中:xi∈Rn,ui∈Rm分别是第i个自主体的状态和输入;E, A∈Rn×n,B∈Rn×m是常数矩阵;(E, A)是正则无脉冲的,且rank E=r≤n.
对于系统(1),定义动态图Gσ(t)=(V, εσ(t)),其中V={1, …,N},且(j, i)∈εσ(t),当且仅当控制ui在时刻t利用(xj-xi)作为反馈.令Aσ(t)=[aij(t)]N×N是动态图Gσ(t)的邻接权矩阵,则可定义状态反馈拓扑
$ {\mathit{\boldsymbol{u}}_i} = \mathit{\boldsymbol{K}}\sum\limits_{j = 1}^i {{a_{ij}}\left( t \right)\left( {{\mathit{\boldsymbol{x}}_j} - {\mathit{\boldsymbol{x}}_i}} \right)} ,i = 1, \cdots ,N, $ | (2) |
这里K ∈Rm×n是增益矩阵.
定义1[8] 无领导者一致性问题.给定系统(1)和一个动态图Gσ(t),找到一个状态反馈拓扑(2)的反馈增益矩阵K,使得对于i, j=1, …, N,当t→∞时,xi(t)-xj(t)→ 0.
对于以上描述的无领导者的一致性问题,每一个子系统解的稳态行为是无足轻重的.还有一个一致性问题称为领导者跟随一致性问题,而这个问题就要求每一个子系统的解都要渐近趋近于信号x0(t).假设信号x0(t)由线性系统
$ \mathit{\boldsymbol{E}}{{\mathit{\boldsymbol{\dot x}}}_0} = \mathit{\boldsymbol{A}}{\mathit{\boldsymbol{x}}_0} $ | (3) |
产生,其中x0∈Rn.
下面分别称系统(1)和系统(3)为跟随者系统和领导者系统.联合系统(1)和系统(3),定义另外一个动态图Gσ(t)=(V, εσ(t)), 这里V={0, 1, …, N}.显然G是G的子图,因为可以从图G移去V中的结点0,和在t时刻εσ(t)中所有属于结点0的边得到.令Δσ(t)是一个N×N的非负对角矩阵,其中第i个对角元是ai0(t),这里如果(0, i)∈ε,则ai0(t)>0,否则ai0(t)=0.考虑状态反馈拓扑
$ {\mathit{\boldsymbol{u}}_i} = \mathit{\boldsymbol{K}}\sum\limits_{j = 0}^N {{a_{ij}}\left( t \right)\left( {{\mathit{\boldsymbol{x}}_j} - {\mathit{\boldsymbol{x}}_i}} \right)} ,i = 1, \cdots ,N. $ | (4) |
定义2[8] 领导者跟随一致性问题.给定领导者系统(3)、跟随者系统(1)和一个动态图Gσ(t),找到一个状态反馈拓扑(4)的反馈增益矩阵K,使得对于i=1, …, N,当t→∞时,xi(t)-x0(t)→0.
为了解决以上两个一致性问题, 我们将系统(1)进行正则性分解.由于(E, A)是正则无脉冲的,由文献[14]可知,存在可逆矩阵P, Q∈Rn×n,使得
$ \mathit{\boldsymbol{PEQ}} = \left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{I}}_r}}&0\\ 0&0 \end{array}} \right],\mathit{\boldsymbol{PAQ}} = \left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{A}}_1}}&0\\ 0&{{\mathit{\boldsymbol{I}}_{n - r}}} \end{array}} \right],\mathit{\boldsymbol{PB}} = \left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{B}}_1}}\\ {{\mathit{\boldsymbol{B}}_2}} \end{array}} \right], $ | (5) |
进行坐标变换,
$ {\mathit{\boldsymbol{Q}}^{ - 1}}{\mathit{\boldsymbol{x}}_i} = {\left[ {\begin{array}{*{20}{c}} {\mathit{\boldsymbol{\bar x}}_{1i}^{\rm{T}}}&{\mathit{\boldsymbol{\bar x}}_{2i}^{\rm{T}}} \end{array}} \right]^{\rm{T}}},i = 1, \cdots ,N, $ | (6) |
其中:x1iT∈Rr, x2iT∈Rn-r.由式(5)得到无领导者系统(1)等价于
$ {{\mathit{\boldsymbol{\dot {\bar x}}}}_{1i}} = {\mathit{\boldsymbol{A}}_1}{{\mathit{\boldsymbol{\bar x}}}_{1i}} + {\mathit{\boldsymbol{B}}_1}{\mathit{\boldsymbol{u}}_i}, $ | (7) |
$ {\bf{0}} = {{\mathit{\boldsymbol{\bar x}}}_{2i}} + {\mathit{\boldsymbol{B}}_2}{\mathit{\boldsymbol{u}}_i},i = 1, \cdots ,N, $ | (8) |
而领导者系统(3)等价于
$ {{\mathit{\boldsymbol{\dot {\bar x}}}}_{10}} = {\mathit{\boldsymbol{A}}_1}{{\mathit{\boldsymbol{\bar x}}}_{10}}, $ | (9) |
$ {\bf{0}} = {{\mathit{\boldsymbol{\bar x}}}_{20}}. $ | (10) |
这种分解通常称为快、慢子系统分解,式(7)和(9)为慢子系统,式(8)和(10)为快子系统.通过这种分解,证明无领导者系统和领导者跟随系统的一致性可分别由相应的慢子系统的一致性来得到.
定理1 如果无领导者慢子系统(7)达到一致性, 则无领导者系统(1)也达到一致性.
证明 假设存在一个增益矩阵K1∈Rm×r,使得
$ {\mathit{\boldsymbol{u}}_i} = {\mathit{\boldsymbol{K}}_1}\sum\limits_{j = 1}^N {{a_{ij}}\left( t \right)\left( {{{\mathit{\boldsymbol{\bar x}}}_{1j}} - {{\mathit{\boldsymbol{\bar x}}}_{1i}}} \right)} ,i = 1, \cdots ,N, $ | (11) |
解决无领导者系统(7)的一致性问题,由定义1,我们有
由式(2)可得,系统(1)的一致性问题也得以解决.
定理2 如果无领导者慢子系统(7)和领导者慢子系统(9)的领导者跟随一致性问题得以解决,则系统(3)和系统(1)的领导者跟随一致性问题也得以解决.
证明 假设存在一个增益矩阵K1∈Rm×r,使得ui=
取K=[K1 0]Q-1,所以系统(3)和(1)的领导者跟随一致性问题也得以解决.
要解决无领导者系统(1)的一致性问题,只需讨论无领导者慢子系统(7)的一致性问题,而系统(7)是一个一般的线性多自主体系统,研究其一致性问题就简单得多.这里考虑关于慢子系统(7)的线性切换系统
$ \mathit{\boldsymbol{\dot {\bar x}}}\left( t \right) = \left( {{\mathit{\boldsymbol{I}}_N} \otimes {\mathit{\boldsymbol{A}}_1} - {\mathit{\boldsymbol{F}}_{\sigma \left( t \right)}} \otimes \left( {{\mathit{\boldsymbol{B}}_1}\mathit{\boldsymbol{B}}_1^{\rm{T}}\mathit{\boldsymbol{X}}} \right)} \right)\mathit{\boldsymbol{\bar x}}\left( t \right),\sigma \left( t \right) \in P, $ | (12) |
其中:X是正定矩阵;x(t)∈Rr; A1、B1如式(5)中定义一样;IN∈RN×N是单位矩阵;σ(t):[0, +∞)→P={1, 2, …, ρ}, ρ≥1,是右连续的分段常值的切换信号,切换瞬时{ti:i=0, 1, …}满足对任意i≥1和正常数τ,都有ti-ti-1≥τ,且对所有的t≥0;Fσ(t)∈RN×N是半正定的矩阵.
假设1 (A1, B1)能控, 令X是一个正定矩阵,满足不等式
$ \mathit{\boldsymbol{A}}_1^{\rm{T}}\mathit{\boldsymbol{X}} + \mathit{\boldsymbol{X}}{\mathit{\boldsymbol{A}}_1} < 0. $ | (13) |
假设2 动态图Gσ(t)是无向图,∀t≥0.
假设3 存在{i:i=0, 1, …}的子序列{ik},tik+1-tik<v, v>0,使得连接图G([tik, tik+1))是连通的.
在假设2下,图Gσ(t)的拉普拉斯矩阵Lσ(t)是对称半正定的,∀t≥0.若一个动态图满足假设3,就称图在[0, ∞)上是一致连通的,或者称在[tik, tik+1)上是共连通的.
注意到
引理1[8] 考虑系统(12),在假设1下,X是满足式(13)的正定矩阵.σ(t)是驻留时间为τ的分段常值切换信号,对任意的t≥0,Fσ(t)是对称半正定的矩阵,则
1) 如果系统(12)的解x(t)满足性质:如果存在{i:i=0, 1, …}的一个子序列ik,tik+1-tik<v, v>0,使得x(tik)与矩阵(
2) 如果存在{i:i=0, 1, …}的一个子序列ik,tik+1-tik<v, v>0,使得矩阵(
定理3 如果假设1~3成立,则系统(1)在控制协议(2)下达到一致性.
证明 设式(11)的增益矩阵为K1=B1TX,系统(7)的第i个自主体的闭环系统为
$ {{\mathit{\boldsymbol{\dot {\bar x}}}}_i}\left( t \right) = {\mathit{\boldsymbol{A}}_1}{{\mathit{\boldsymbol{\bar x}}}_i} + {\mathit{\boldsymbol{B}}_1}\mathit{\boldsymbol{B}}_1^{\rm{T}}\mathit{\boldsymbol{X}}\sum\limits_{j = 1}^N {{a_{ij}}\left( t \right)\left( {{{\mathit{\boldsymbol{\bar x}}}_j} - {{\mathit{\boldsymbol{\bar x}}}_i}} \right)} . $ | (14) |
令xc(t)=(x1(t)+x2(t)+…+xN(t))/N,xc(t)称在t时刻所有自主体的中心.图是无向图,得
$ {{\dot x}_c}\left( t \right) = \left( {\sum\limits_{i = 1}^N {{{\mathit{\boldsymbol{\dot {\bar x}}}}_i}\left( t \right)} } \right)/N = {\mathit{\boldsymbol{A}}_1}\left( {\left( {\sum\limits_{i = 1}^N {{{\mathit{\boldsymbol{\bar x}}}_i}\left( t \right)} } \right)/N} \right) = {\mathit{\boldsymbol{A}}_1}{\mathit{\boldsymbol{x}}_c}\left( t \right), $ | (15) |
将xi(t)分解为
$ {{\mathit{\boldsymbol{\bar x}}}_i}\left( t \right) = {\mathit{\boldsymbol{x}}_c}\left( t \right) + {\mathit{\boldsymbol{\omega }}_i}\left( t \right),i = 1, \cdots ,N, $ | (16) |
则
$ \mathit{\boldsymbol{\bar x}}\left( t \right) = {{\bf{1}}_N} \otimes {\mathit{\boldsymbol{x}}_c}\left( t \right) + \mathit{\boldsymbol{\omega }}\left( t \right), $ | (17) |
其中x(t)=[x1(t)T, x2(t)T, …, xN(t)T]T,且ω(t)=[ω1(t)T, ω2(t)T, …, ωN(t)T]T.由式(14)、(15)和(17)得向量ω(t)满足
例1 考虑无领导者广义多自主体系统(1),N=4,系统矩阵为
$ \mathit{\boldsymbol{E}} = \left[ {\begin{array}{*{20}{c}} 1&0&0&1\\ 0&1&0&0\\ 0&0&1&0\\ 0&0&0&0 \end{array}} \right],\mathit{\boldsymbol{A}} = \left[ {\begin{array}{*{20}{c}} { - 1}&1&1&{ - 1}\\ { - 1}&{ - 1}&{ - 1}&{ - 1}\\ 1&1&{ - 1}&1\\ 0&0&0&{ - 1} \end{array}} \right],\mathit{\boldsymbol{B}} = \left[ {\begin{array}{*{20}{c}} {0.2}\\ {0.2}\\ {0.2}\\ {0.2} \end{array}} \right]. $ |
按(5)式将系统矩阵进行正则性分解,得到PEQ=
$ {\mathit{\boldsymbol{A}}_1} = \left[ {\begin{array}{*{20}{c}} { - 1}&1&1\\ { - 1}&{ - 1}&{ - 1}\\ 1&1&{ - 1} \end{array}} \right];{\mathit{\boldsymbol{B}}_1} = \left[ {\begin{array}{*{20}{c}} {0.2}\\ {0.2}\\ {0.2} \end{array}} \right];{\mathit{\boldsymbol{B}}_2} = 0.2;\mathit{\boldsymbol{P}} = {\mathit{\boldsymbol{I}}_4};\mathit{\boldsymbol{Q}} = \left[ {\begin{array}{*{20}{c}} 1&0&0&1\\ 0&1&0&0\\ 0&0&1&0\\ 0&0&0&0 \end{array}} \right]. $ |
切换动态图Gσ(t),由分段常值的切换信号定义
$ \sigma \left( t \right) = \left\{ {\begin{array}{*{20}{c}} \begin{array}{l} 1,\\ 2,\\ 3,\\ 4, \end{array}&\begin{array}{l} 若\;2s \le t < 2s + 0.5,\\ 若\;2s + 0.5 \le t < 2s + 1,\\ 若\;2s + 1 \le t < 2s + 1.5,\\ 若\;2s + 1.5 \le t < 2s + 2, \end{array} \end{array}} \right. $ |
其中:s=0, 1, 2, …;切换拓扑图Gi=(Vi, εi, Ai),i=1, 2, 3, 4,ε1={(1, 2), (2, 1)},ε2={(2, 3), (3, 2)},ε3={(3, 4), (4, 3)},ε4={(1, 4), (4, 1)}.没有一个图Gi是连通的,但是对于任意的s=0, 1, …, ∪i=2 s2 s+1.5,Gi是连通的.因此假设2和假设3成立,可以证明假设1也成立.解不等式A1TX+XA1<0,得到正定解
$ \mathit{\boldsymbol{X}} = \left[ {\begin{array}{*{20}{c}} {0.309\;0}&{0.075\;7}&{ - 0.200\;8}\\ {0.075\;7}&{0.302\;9}&{0.075\;7}\\ { - 0.200\;8}&{0.075\;7}&{0.309\;0} \end{array}} \right] \times {10^{ - 11}}, $ |
考虑慢子系统的线性切换(12),由定理1,控制增益为K1=B1TX=[0.367 7 0.908 7 0.367 7]×10-12,因此K=[K1 0]Q-1=[0.367 7 0.908 7 0.367 7 0.367 7]×10-12.
最后得到广义系统的状态与中心的误差趋近于零,即系统最终达到了一致性.
4 结论本文讨论了一类具有切换拓扑的广义多自主体系统,利用代数拓扑和广义系统理论知识,将广义多自主体系统进行快慢子系统分解,通过研究其慢子系统的性质,即设计了状态反馈控制协议,得到慢子系统的一致性,从而得到广义多自主体系统的一致性问题,包括无领导者一致性问题和领导者跟随一致性问题.
[1] |
LIU Q, JIANG D Q, SHI N Z, et al. Asymptotic behavior of stochastic multi-group epidemic models with distributed delays[J]. Physica a: statistical mechanics and its applications, 2017, 467: 527-541. DOI:10.1016/j.physa.2016.10.034 ( ![]() |
[2] |
XIAO L, BOYD S. Fast linear iterations for distributed averaging[J]. Systems and control letters, 2004, 53(1): 65-78. DOI:10.1016/j.sysconle.2004.02.022 ( ![]() |
[3] |
LI T, ZHANG J. Consensus conditions of multi-agent system with time-varying topologies and stochastic communication noises[J]. IEEE transactions on automatic control, 2010, 55(9): 2043-2057. DOI:10.1109/TAC.2010.2042982 ( ![]() |
[4] |
MA C, ZHANG J. Necessary and sufficient conditions for consensus ability of linear multi-agent systems[J]. IEEE transactions on automatic control, 2010, 55(5): 1263-1268. DOI:10.1109/TAC.2010.2042764 ( ![]() |
[5] |
QIN J, ZHENG W, GAO H. Consensus of multiple second-order vehicles with a time-varying reference signal under directed topology[J]. Automatica, 2011, 47(9): 1983-1991. DOI:10.1016/j.automatica.2011.05.014 ( ![]() |
[6] |
NI W, CHENG D. Leader-following consensus of multi-agent systems under fixed and switching topologies[J]. Systems and control letters, 2010, 59(3): 209-217. ( ![]() |
[7] |
OLFATI-SABER R, MURRAY R M. Consensus problems in networks of agents with switching topology and time-delays[J]. IEEE transactions on automatic control, 2004, 49(9): 1520-1533. DOI:10.1109/TAC.2004.834113 ( ![]() |
[8] |
SU Y, HUANG J. Stability of a class of linear switching systems with applications to two consensus problems[J]. IEEE transactions on automatic control, 2012, 57(6): 1420-1430. DOI:10.1109/TAC.2011.2176391 ( ![]() |
[9] |
DUAN G R. Analysis and design of descriptor linear systems[M]. New York: Springer-Verlag, 2010.
( ![]() |
[10] |
LIN Z Y, FRANCIS B, MAGGIORE M. Necessary and sufficient graphical conditions for formation control of unicycles[J]. IEEE transactions on automatic control, 2005, 50(1): 121-127. DOI:10.1109/TAC.2004.841121 ( ![]() |
[11] |
HONG Y, HU J, GAO L. Tracking control for multi-agent consensus with an active leader and variable topology[J]. Automatica, 2006, 42(7): 1177-1182. DOI:10.1016/j.automatica.2006.02.013 ( ![]() |
[12] |
YANG X R, LIU G P. Necessary and sufficient consensus conditions of descriptor multi-agent systems[J]. IEEE transactions on circuits and systems-Ⅰ, 2012, 59(11): 2669-2677. DOI:10.1109/TCSI.2012.2190663 ( ![]() |
[13] |
GODSIL C, ROYLE G. Algebraic graph theory[M]. Berlin: Springer-Verlag, 2001.
( ![]() |
[14] |
YANG D M, ZHANG Q L, YAO B. Singular systems[M]. Beijing: Science Press, 2004.
( ![]() |