有向图上的多智能体同步(下)
如何在有向图上实现同步
图论代数
广义代数连通性
定义1 对于具有拉普拉斯矩阵 L 的强连通网络,广义代数连通性定义为 a ( L ) = min  x T ξ = 0 , x ≠ 0 x T L ^ x x T Ξ x a(L) = \min_{\substack{x^{\mathrm{T}}\xi = 0, x \neq 0}} \frac{x^{\mathrm{T}}\widehat{L}x}{x^{\mathrm{T}}\Xi x} a(L)=xTξ=0,x=0minxTΞxxTL
                   x其中, L ^ = 1 2 ( Ξ L + L T Ξ ) \widehat{L} = \frac{1}{2}(\Xi L + L^{\mathrm{T}}\Xi) L
           =21(ΞL+LTΞ), Ξ = d i a g ( ξ 1 , … , ξ n ) \Xi = \mathrm{diag}(\xi_1, \dots, \xi_n) Ξ=diag(ξ1,…,ξn), ξ = ( ξ 1 , ξ 2 , … , ξ n ) T > 0 \xi = (\xi_1, \xi_2, \dots, \xi_n)^{\mathrm{T}} > 0 ξ=(ξ1,ξ2,…,ξn)T>0,且满足  ξ T L = 0 \xi^{\mathrm{T}}L = 0 ξTL=0, ∑ i = 1 n ξ i = 1 \sum_{i=1}^n \xi_i = 1 ∑i=1nξi=1。注意,若  Ξ = η I N \Xi = \eta I_N Ξ=ηIN( η \eta η 为常数, I n I_n In 为  n n n 阶单位矩阵)且网络是无向的,则  a ( L ) = λ 2 ( L ) a(L) = \lambda_2(L) a(L)=λ2(L)(即无向网络的第二小拉普拉斯特征值,传统代数连通性)。
引理 4:假设矩阵  L ^ \widehat{L} L
            是对称的、不可约的,且满足  ∑ j = 1 N L ^ i j = 0 \sum_{j=1}^N \widehat{L}_{ij} = 0 ∑j=1NL
            ij=0(行和为 0),同时  L ^ i j ≥ 0 \widehat{L}_{ij} \geq 0 L
            ij≥0( i ≠ j i \neq j i=j, i , j = 1 , 2 , … , n i,j = 1,2,\dots,n i,j=1,2,…,n)。令 a ^ ( L ^ ) = min  x T ξ = 0 , x ≠ 0 x T L ^ x x T x \widehat{a}(\widehat{L}) = \min_{\substack{x^{\mathrm{T}}\xi = 0, x \neq 0}} \frac{x^{\mathrm{T}}\widehat{L}x}{x^{\mathrm{T}}x} a
            (L
            )=xTξ=0,x=0minxTxxTL
                   x则有  λ 2 ( L ^ ) ≥ a ^ ( L ^ ) ≥ 0 \lambda_2(\widehat{L}) \geq \widehat{a}(\widehat{L}) \geq 0 λ2(L
           )≥a
           (L
           )≥0。此外, a ^ ( L ^ ) = 0 \widehat{a}(\widehat{L}) = 0 a
           (L
           )=0 当且仅当  ξ \xi ξ 与  L ^ \widehat{L} L
            对应于特征值 0 的左特征向量正交;并且,若  ξ \xi ξ 是  L ^ \widehat{L} L
            对应于特征值 0 的左特征向量,则  a ^ ( L ^ ) = λ 2 ( L ^ ) \widehat{a}(\widehat{L}) = \lambda_2(\widehat{L}) a
           (L
           )=λ2(L
           )。
证明 由于  L ^ \widehat{L} L
            对称且行和为零,其最小特征值为  λ 1 ( L ^ ) = 0 \lambda_1(\widehat{L}) = 0 λ1(L
           )=0,对应的特征向量为  1 n \mathbf{1}_n 1n。由 Courant–Fischer 极小极大原理,第二小特征值可表示为
 λ 2 ( L ^ ) = min  x ≠ 0 x ⊤ 1 N = 0 x ⊤ L ^ x x ⊤ x . \lambda_2(\widehat{L}) = \min_{\substack{x \neq 0 \\ x^\top \mathbf{1}_N = 0}} \frac{x^\top \widehat{L} x}{x^\top x}. λ2(L
            )=x=0x⊤1N=0minx⊤xx⊤L
                   x.
  注意到约束条件  x ⊤ ξ = 0 x^\top \xi = 0 x⊤ξ=0 定义了一个超平面,而  x ⊤ 1 n = 0 x^\top \mathbf{1}_n = 0 x⊤1n=0 定义了另一个超平面。一般情况下,这两个超平面不同,因此在更受限的集合  { x ∣ x ⊤ ξ = 0 , x ≠ 0 } \{x \mid x^\top \xi = 0, x \neq 0\} {x∣x⊤ξ=0,x=0} 上取最小值,所得结果不会小于在  { x ∣ x ⊤ 1 N = 0 , x ≠ 0 } \{x \mid x^\top \mathbf{1}_N = 0, x \neq 0\} {x∣x⊤1N=0,x=0} 上的最小值。故有
 λ 2 ( L ^ ) ≥ a ^ ( L ^ ) ≥ 0. \lambda_2(\hat{L}) \geq \hat{a}(\hat{L}) \geq 0. λ2(L^)≥a^(L^)≥0.
进一步,若 ξ \xi ξ 与 1 n \mathbf{1}_n 1n 正交(即 ξ ⊤ 1 n = 0 \xi^\top \mathbf{1}_n = 0 ξ⊤1n=0),则存在非零向量 x x x 同时满足 x ⊤ ξ = 0 x^\top \xi = 0 x⊤ξ=0 和 x ∈ s p a n { 1 n } x \in \mathrm{span}\{\mathbf{1}_n\} x∈span{1n},此时 x ⊤ L ^ x = 0 x^\top \widehat{L} x = 0 x⊤L x=0,从而 a ^ ( L ^ ) = 0 \widehat{a}(\widehat{L}) = 0 a (L )=0。反之,若 a ^ ( L ^ ) = 0 \hat{a}(\hat{L}) = 0 a^(L^)=0,则存在非零 x x x 满足 x ⊤ ξ = 0 x^\top \xi = 0 x⊤ξ=0 且 x ⊤ L ^ x = 0 x^\top \widehat{L} x = 0 x⊤L x=0,即 x x x 属于 L ^ \widehat{L} L 的零空间,故 x ∝ 1 N x \propto \mathbf{1}_N x∝1N,从而 ξ ⊤ 1 N = 0 \xi^\top \mathbf{1}_N = 0 ξ⊤1N=0。
特别地,若 ξ = 1 n / n \xi = \mathbf{1}_n / n ξ=1n/n(即 ξ \xi ξ 是 L ^ \widehat{L} L 对应于零特征值的归一化左特征向量),则约束 x ⊤ ξ = 0 x^\top \xi = 0 x⊤ξ=0 等价于 x ⊤ 1 n = 0 x^\top \mathbf{1}_n = 0 x⊤1n=0,此时极小化问题与 λ 2 ( L ^ ) \lambda_2(\widehat{L}) λ2(L ) 的定义完全一致,故 a ^ ( L ^ ) = λ 2 ( L ^ ) \widehat{a}(\widehat{L}) = \lambda_2(\hat{L}) a (L )=λ2(L^)。 □ \quad\square □
引理 5 若 Laplacian 矩阵 L L L 不可约(即对应的有向图强连通),则广义代数连通度 a ( L ) > 0 a(L) > 0 a(L)>0。
证明 当  L L L 不可约时,存在正向量  ξ = ( ξ 1 , … , ξ n ) ⊤ > 0 \xi = (\xi_1, \dots, \xi_n)^\top > 0 ξ=(ξ1,…,ξn)⊤>0 满足  ξ ⊤ L = 0 \xi^\top L = 0 ξ⊤L=0 且  ξ ⊤ 1 n = 1 \xi^\top \mathbf{1}_n = 1 ξ⊤1n=1。定义对角矩阵  Ξ = d i a g ( ξ 1 , … , ξ N ) \Xi = \mathrm{diag}(\xi_1, \dots, \xi_N) Ξ=diag(ξ1,…,ξN),并令
 L ^ = 1 2 ( Ξ L + L ⊤ Ξ ) . \widehat{L} = \frac{1}{2}(\Xi L + L^\top \Xi). L
            =21(ΞL+L⊤Ξ).
  易见  L ^ \widehat{L} L
            为对称矩阵,且满足  L ^ 1 n = 0 \widehat{L} \mathbf{1}_n = 0 L
           1n=0,即其行和与列和均为零。此外,由于  L L L 不可约且  ξ > 0 \xi > 0 ξ>0,可证  L ^ \hat{L} L^ 亦不可约(否则将导致  L L L 可约,矛盾)。
根据定义1,广义代数连通度为
 a ( L ) = min  x ⊤ ξ = 0 x ≠ 0 x ⊤ L ^ x x ⊤ Ξ x . a(L) = \min_{\substack{x^\top \xi = 0 \\ x \neq 0}} \frac{x^\top \widehat{L} x}{x^\top \Xi x}. a(L)=x⊤ξ=0x=0minx⊤Ξxx⊤L
                   x.
  注意到  Ξ \Xi Ξ 正定,故对任意非零  x x x 满足  x ⊤ ξ = 0 x^\top \xi = 0 x⊤ξ=0,分母  x ⊤ Ξ x > 0 x^\top \Xi x > 0 x⊤Ξx>0。又因  L ^ \hat{L} L^ 对称半正定且零空间为  s p a n { 1 N } \mathrm{span}\{\mathbf{1}_N\} span{1N},而约束  x ⊤ ξ = 0 x^\top \xi = 0 x⊤ξ=0 与  1 n \mathbf{1}_n 1n 不正交(因  ξ ⊤ 1 n = 1 ≠ 0 \xi^\top \mathbf{1}_n = 1 \neq 0 ξ⊤1n=1=0),故在该约束下  x ∉ k e r ( L ^ ) x \notin \mathrm{ker}(\hat{L}) x∈/ker(L^),从而  x ⊤ L ^ x > 0 x^\top \hat{L} x > 0 x⊤L^x>0。
进一步,利用引理4中的辅助量  a ^ ( L ^ ) = min  x ⊤ ξ = 0 , x ≠ 0 x ⊤ L ^ x x ⊤ x \widehat{a}(\widehat{L}) = \min_{x^\top \xi = 0, x \neq 0} \frac{x^\top \widehat{L} x}{x^\top x} a
           (L
           )=minx⊤ξ=0,x=0x⊤xx⊤L
                   x,可得
 a ( L ) = min  x ⊤ ξ = 0 , x ≠ 0 x ⊤ L ^ x x ⊤ Ξ x ≥ 1 max  i ξ i min  x ⊤ ξ = 0 , x ≠ 0 x ⊤ L ^ x x ⊤ x = a ^ ( L ^ ) max  i ξ i . a(L) = \min_{x^\top \xi = 0, x \neq 0} \frac{x^\top \widehat{L} x}{x^\top \Xi x} \geq \frac{1}{\max_i \xi_i} \min_{x^\top \xi = 0, x \neq 0} \frac{x^\top \widehat{L} x}{x^\top x} = \frac{\widehat{a}(\widehat{L})}{\max_i \xi_i}. a(L)=x⊤ξ=0,x=0minx⊤Ξxx⊤L
                   x≥maxiξi1x⊤ξ=0,x=0minx⊤xx⊤L
                   x=maxiξia
                   (L
                   ).
  由于  L ^ \hat{L} L^ 不可约且  ξ ⊤ 1 N = 1 ≠ 0 \xi^\top \mathbf{1}_N = 1 \neq 0 ξ⊤1N=1=0,由引理6知  a ^ ( L ^ ) > 0 \hat{a}(\hat{L}) > 0 a^(L^)>0,故  a ( L ) > 0 a(L) > 0 a(L)>0。  □ \quad \square □
多智能体同步
在多智能体协同控制中,一致性(consensus)是最基本且核心的问题之一:每个智能体仅能获取局部邻居信息,却希望全局状态最终趋于一致。例如,一群无人机,彼此只能“单向”接收邻居的位置信息(比如 A 能看到 B,但 B 看不到 A),它们还能不能最终排成整齐队形?答案是肯定的,只要通信图是“强连通”的,就能同步。
模型
考虑  N N N 个智能体,每个状态为  x i ( t ) ∈ R x_i(t) \in \mathbb{R} xi(t)∈R,遵循如下规则:
 x ˙ i ( t ) = ∑ j = 1 N a i j ( x j ( t ) − x i ( t ) ) , \dot{x}_i(t) = \sum_{j=1}^N a_{ij} \big( x_j(t) - x_i(t) \big), x˙i(t)=j=1∑Naij(xj(t)−xi(t)),
  其中  a i j ≥ 0 a_{ij} \geq 0 aij≥0 表示智能体  i i i 对  j j j 的“信任权重”。若  a i j > 0 a_{ij} > 0 aij>0,说明存在一条从  j j j 指向  i i i 的有向边( i i i 能接收  j j j 的信息)。
这个模型的物理意义非常直观:每个智能体把自己的状态向邻居的加权平均靠拢。它是一阶共识(first-order consensus)的标准形式,也是更复杂编队、协同控制的基础。
令 Laplacian 矩阵  L = D − A L = D - A L=D−A( D D D 为出度对角阵),系统可写为:
 x ˙ = − L x . \dot{x} = -L x. x˙=−Lx.
同步的定义
我们说系统达成同步(consensus),若
 lim  t → ∞ ∣ x i ( t ) − x j ( t ) ∣ = 0 , ∀ i , j . \lim_{t \to \infty} |x_i(t) - x_j(t)| = 0, \quad \forall i,j. t→∞lim∣xi(t)−xj(t)∣=0,∀i,j.
在无向图中,极限是算术平均  1 N ∑ x i ( 0 ) \frac{1}{N}\sum x_i(0) N1∑xi(0)。但在有向图中,由于信息流不对称,极限是一个加权平均:
 x ˉ = ξ ⊤ x ( 0 ) , \bar{x} = \xi^\top x(0), xˉ=ξ⊤x(0),
  其中  ξ > 0 \xi > 0 ξ>0 是  L L L 的左零特征向量,满足  ξ ⊤ L = 0 \xi^\top L = 0 ξ⊤L=0 且  ξ ⊤ 1 N = 1 \xi^\top \mathbf{1}_N = 1 ξ⊤1N=1。这个  ξ i \xi_i ξi 可理解为智能体  i i i 在网络中的“影响力权重”——越能被他人间接影响,权重越高。
结论与证明
定理一 当通信图是强连通(strongly connected)时,即任意两个节点间都存在有向路径,则系统是指数同步的。
证明 注意上一节的定义1,广义代数连通度:
 a ( L ) = min  x ⊤ ξ = 0 x ≠ 0 x ⊤ L ^ x x ⊤ Ξ x 。 a(L) = \min_{\substack{x^\top \xi = 0 \\ x \ne 0}} \frac{x^\top \hat{L} x}{x^\top \Xi x}。 a(L)=x⊤ξ=0x=0minx⊤Ξxx⊤L^x。
  同时,引理5证明,只要图强连通,就有  a ( L ) > 0 a(L) > 0 a(L)>0。 此时,定义误差  x ^ = x − 1 N x ˉ \hat{x} = x - \mathbf{1}_N \bar{x} x^=x−1Nxˉ,则  ξ ⊤ x ^ = 0 \xi^\top \hat{x} = 0 ξ⊤x^=0,且  x ^ ˙ = − L x ^ \dot{\hat{x}} = -L \hat{x} x^˙=−Lx^。构造 Lyapunov 函数:
 V ( t ) = x ^ ⊤ Ξ x ^ ≥ 0 , V(t) = \hat{x}^\top \Xi \hat{x} \geq 0, V(t)=x^⊤Ξx^≥0,
  其导数为:
 V ˙ = − 2 x ^ ⊤ L ^ x ^ ≤ − 2 a ( L ) x ^ ⊤ Ξ x ^ = − 2 a ( L ) V ( t ) . \dot{V} = -2 \hat{x}^\top \hat{L} \hat{x} \leq -2 a(L) \hat{x}^\top \Xi \hat{x} = -2 a(L) V(t). V˙=−2x^⊤L^x^≤−2a(L)x^⊤Ξx^=−2a(L)V(t).
于是 V ( t ) ≤ V ( 0 ) e − 2 a ( L ) t → 0 V(t) \leq V(0) e^{-2a(L)t} \to 0 V(t)≤V(0)e−2a(L)t→0,即指数同步!
注 当图无向时, L L L 对称, ξ = 1 N 1 N \xi = \frac{1}{N}\mathbf{1}_N ξ=N11N, Ξ = 1 N I \Xi = \frac{1}{N}I Ξ=N1I, L ^ = L \hat{L} = L L^=L,此时
 a ( L ) = λ 2 ( L ) , a(L) = \lambda_2(L), a(L)=λ2(L),
  完全退化为经典结果。
潜在应用
这项关于有向图同步的理论,可以直接回应了现实世界中大量非对称通信场景的控制需求。 同步与否,不取决于通信是否“双向”,而取决于信息能否“遍历全网”。我们可以将模型用在例如无人机编队、机器人控制、网络传输以及网络算法的收敛和鲁棒性上。其实也有点像团队协作,每个人做好自己的局部沟通,整个团队就会实现全局目标。
PS:可以用在信息传播上,如果大家都转发和点赞文章,就能促进大家的共同进步!
更多推荐
 

所有评论(0)