图论代数
广义代数连通性
定义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:可以用在信息传播上,如果大家都转发和点赞文章,就能促进大家的共同进步!
2621

被折叠的 条评论
为什么被折叠?



