上接有向图上的多智能体同步(上)

图论代数

广义代数连通性

定义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 ij0 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=0x1N=0minxxxL x.
注意到约束条件 x ⊤ ξ = 0 x^\top \xi = 0 xξ=0 定义了一个超平面,而 x ⊤ 1 n = 0 x^\top \mathbf{1}_n = 0 x1n=0 定义了另一个超平面。一般情况下,这两个超平面不同,因此在更受限的集合 { x ∣ x ⊤ ξ = 0 , x ≠ 0 } \{x \mid x^\top \xi = 0, x \neq 0\} {xxξ=0,x=0} 上取最小值,所得结果不会小于在 { x ∣ x ⊤ 1 N = 0 , x ≠ 0 } \{x \mid x^\top \mathbf{1}_N = 0, x \neq 0\} {xx1N=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\} xspan{1n},此时 x ⊤ L ^ x = 0 x^\top \widehat{L} x = 0 xL 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 xL x=0,即 x x x 属于 L ^ \widehat{L} L 的零空间,故 x ∝ 1 N x \propto \mathbf{1}_N x1N,从而 ξ ⊤ 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 x1n=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ΞxxL 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 xL^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=0xxxL 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ΞxxL xmaxiξi1xξ=0,x=0minxxxL 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=1Naij(xj(t)xi(t)),
其中 a i j ≥ 0 a_{ij} \geq 0 aij0 表示智能体 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=DA 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. tlimxi(t)xj(t)=0,i,j.

在无向图中,极限是算术平均 1 N ∑ x i ( 0 ) \frac{1}{N}\sum x_i(0) N1xi(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ΞxxL^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^=x1Nxˉ,则 ξ ⊤ 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)e2a(L)t0,即指数同步

当图无向时, 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:可以用在信息传播上,如果大家都转发和点赞文章,就能促进大家的共同进步!

Logo

有“AI”的1024 = 2048,欢迎大家加入2048 AI社区

更多推荐