图论与网络基础

1. 图的表示

G=(V,E,A)G=(V,E,A)G=(V,E,A)表示一个图三元组,其中V={v1,...,vn}V=\{v_1,...,v_n\}V={v1,...,vn}表示节点的集合,E⊆V×VE⊆V×VEV×V是有向边的集合,A=(aij)n×nA=(a_{ij})_{n\times n}A=(aij)n×n 是加权邻接矩阵,用来量化节点之间的关系。aij>0a_{ij}>0aij>0表示节点vjv_jvj接受到viv_ivi的信息的大小,也就是有一个边从viv_ivi出发链接到vjv_jvj(也有部分文献是反过来写的)。

2. 有向图和无向图

如果对于图中任意两个节点viv_ivivjv_jvj ,若存在连接关系就有aij=aji>0a_{ij} = a_{ji} > 0aij=aji>0 ,否则aij=aji=0a_{ij} = a_{ji} = 0aij=aji=0, (i≠ji \neq ji=ji,j=1,2,…,ni, j = 1, 2, \ldots, ni,j=1,2,,n),这样的图就是无向图;若从节点viv_ivivjv_jvj 存在连接时aij>0a_{ij} > 0aij>0 ,否则aij=0a_{ij} = 0aij=0i≠ji \neq ji=ji,j=1,2,…,ni, j = 1, 2, \ldots, ni,j=1,2,,n),则是有向图。可以看出,无向图其实是有向图的一种特殊情况。

3. 连通性

在图论中,从节点vjv_jvjviv_ivi 的有向(或无向)路径是有向(或无向)网络中由不同节点vikv_{ik}vikk=1,2,…,lk = 1, 2, \ldots, lk=1,2,,l)组成的边序列(vi,vi1),(vi1,vi2),…,(vil,vj)(v_i, v_{i1}), (v_{i1}, v_{i2}), \ldots, (v_{il}, v_j)(vi,vi1),(vi1,vi2),,(vil,vj)。如果一个有向(或无向)网络GGG中,任意两个不同节点viv_ivivjv_jvj 之间都存在从viv_ivivjv_jvj 的有向(或无向)路径,那么这个网络就是强连通(或连通)的。

4. 拉普拉斯矩阵

在图论和机器学习(尤其是图神经网络)的研究中,拉普拉斯矩阵(Laplacian Matrix)是连接图结构与代数运算的核心工具。对于矩阵AAA,其队形的拉普拉斯矩阵 L 是 n×nn \times nn×n的实矩阵,记为 L=(lij)∈Rn×nL = (l_{ij}) \in \mathbb{R}^{n \times n}L=(lij)Rn×n,它的元素定义遵循两个规则:当 i≠ji \neq ji=j(即非对角元素)时,lij=−aijl_{ij} = -a_{ij}lij=aij,直接取邻接矩阵对应元素的负值;当 i=ji = ji=j(即对角元素)时,lii=∑j=1,j≠inaijl_{ii} = \sum_{j=1, j \neq i}^n a_{ij}lii=j=1,j=inaij,也就是邻接矩阵中第 i 行所有非对角元素的和,这个值其实就是节点 i 的 “度”(Degree)。从这个定义能直接推导出拉普拉斯矩阵的一个关键性质:它的每一行元素之和都为 0,即 ∑j=1Nlij=0\sum_{j=1}^N l_{ij} = 0j=1Nlij=0
举个简单的例子帮助理解:假设一个包含 3 个节点的无向图,节点 1 与 2互相 相连、节点 1 与 3 互相相连,邻接矩阵 A=(011100100)A = \begin{pmatrix}0&1&1\\1&0&0\\1&0&0\end{pmatrix}A= 011100100 按照定义计算,拉普拉斯矩阵 L=(2−1−1−110−101)L = \begin{pmatrix}2&-1&-1\\-1&1&0\\-1&0&1\end{pmatrix}L= 211110101 可以看到每一行元素相加都等于 0,完全符合定义要求。

5. 不可约矩阵

一个 n×nn \times nn×n 的实矩阵 A(通常讨论非负矩阵),若不存在置换矩阵 P,使得 PTAPP^TAPPTAP 化为如下分块上三角形式:PTAP=(A11A120A22)P^TAP = \begin{pmatrix} A_{11} & A_{12} \\ 0 & A_{22} \end{pmatrix}PTAP=(A110A12A22)其中 A11A_{11}A11A22A_{22}A22 是阶数非零的方阵,则称 A 为不可约矩阵;反之则为可约矩阵。

6.树

若一个有向网络忽略边的方向后,其底层结构是一棵树,则该网络称为有向树(directed tree)。有向根树(directed rooted tree) 是一种有向网络,它至少存在一个根节点r,且满足:对于除 rrr 之外的任意节点 vvv,都存在从rrrvvv的唯一有向路径。网络 GGG 的有向生成树(directed spanning tree) 是一棵有向根树,它包 GGG 中的所有节点和部分边。
不难看出,有向生成树是强链接图(强连通图)的 “必要非充分条件”,即强链接图一定包含有向生成树,但存在有向生成树的图不一定是强链接图。同时,若有向图 GGG 是强连通图(强链接),则 GGG 中必然存在以任意节点为根的有向生成树(出树或入树)。

图论代数

引理 1 一个矩阵A是不可约的,当且仅当其对应的图G是强链接的。

R. A. Brualdi and H. J. Ryser, Combinatorial Matrix Theory. Cambridge, U.K.: Cambridge Univ. Press, 1991.

引理 2 矩阵A是不可约的,其对应的拉普拉斯矩阵LLL也是不可约的。
引理 3 若拉普拉斯矩阵LLL是不可约的,则有以下性质成立:

  • L⋅1n=0L \cdot \mathbf{1}_n = \mathbf{0}L1n=0(其中1n\mathbf{1}_n1n是元素全为 1 的N维列向量);
  • 存在正向量ξ=(ξ1,ξ2,…,ξn)T\boldsymbol{\xi} = (\xi_1, \xi_2, \ldots, \xi_n)^Tξ=(ξ1,ξ2,,ξn)T,使得ξT⋅L=0\boldsymbol{\xi}^T \cdot L = \mathbf{0}ξTL=0
  • 存在正定对角矩阵Ξ=diag(ξ1,ξ2,…,ξn)\Xi = \text{diag}(\xi_1, \xi_2, \ldots, \xi_n)Ξ=diag(ξ1,ξ2,,ξn),使得L~=12(ΞL+LTΞ)\widetilde{L} = \frac{1}{2} \left( \Xi L + L^T \Xi \right)L =21(ΞL+LTΞ)为对称矩阵,且对于所有i=1,2,…,Ni = 1, 2, \ldots, Ni=1,2,,N,均满足:∑j=1nL~ij=∑j=1nL~ji=0\sum_{j=1}^n \widetilde{L}_{ij} = \sum_{j=1}^n \widetilde{L}_{ji} = 0j=1nL ij=j=1nL ji=0

W. Yu, G. Chen, M. Cao, and J. Kurths, “Second-order consensus for multiagent systemsWith directed topologies and nonlinear dynamics,” IEEE Transactions on Systems Man & Cybernetics Part B: Cybernetics, vol. 40, no. 3, pp. 881-891, June 2010.

这里给个例子给大家直观体现引理3:
(1)邻接矩阵A(非对称)A=(0230000104005000)A = \begin{pmatrix} 0 & 2 & 3 & 0 \\ % 节点1的出边:2→2,3→3 0 & 0 & 0 & 1 \\ % 节点2的出边:1→4 0 & 4 & 0 & 0 \\ % 节点3的出边:4→2 5 & 0 & 0 & 0 % 节点4的出边:5→1 \end{pmatrix}A= 0005204030000100
(2)不可约拉普拉斯矩阵L对角元lii=∑j≠iaijl_{ii} = \sum_{j \neq i} a_{ij}lii=j=iaij(节点出度和):l11=2+3=5l_{11}=2+3=5l11=2+3=5l22=1l_{22}=1l22=1l33=4l_{33}=4l33=4l44=5l_{44}=5l44=5非对角元lij=−aijl_{ij} = -a_{ij}lij=aiji≠ji \neq ji=jL=(5−2−30010−10−440−5005)L = \begin{pmatrix} 5 & -2 & -3 & 0 \\ % 节点1:出度5,指向2(-2)、3(-3) 0 & 1 & 0 & -1 \\ % 节点2:出度1,指向4(-1) 0 & -4 & 4 & 0 \\ % 节点3:出度4,指向2(-4) -5 & 0 & 0 & 5 % 节点4:出度5,指向1(-5) \end{pmatrix}L= 5005214030400105
性质 1:验证不可约拉普拉斯矩阵的性质性质 1:L⋅14=0L \cdot \mathbf{1}_4 = \mathbf{0}L14=0, 14=(1,1,1,1)T\mathbf{1}_4 = (1,1,1,1)^T14=(1,1,1,1)T,计算:L⋅14=(5−2−3+00+1+0−10−4+4+0−5+0+0+5)=(0000)=0L \cdot \mathbf{1}_4 = \begin{pmatrix} 5-2-3+0 \\ 0+1+0-1 \\ 0-4+4+0 \\ -5+0+0+5 \end{pmatrix} = \begin{pmatrix}0 \\ 0 \\ 0 \\ 0\end{pmatrix} = \mathbf{0}L14= 523+00+1+0104+4+05+0+0+5 = 0000 =0成立。
性质 2:存在正向量ξ=(ξ1,ξ2,ξ3,ξ4)T\boldsymbol{\xi} = (\xi_1,\xi_2,\xi_3,\xi_4)^Tξ=(ξ1,ξ2,ξ3,ξ4)T,使得ξTL=0\boldsymbol{\xi}^T L = \mathbf{0}ξTL=0解方程组ξTL=(0,0,0,0)\boldsymbol{\xi}^T L = (0,0,0,0)ξTL=(0,0,0,0){5ξ1−5ξ4=0−2ξ1+ξ2−4ξ3=0−3ξ1+4ξ3=0−ξ2+5ξ4=0\begin{cases} 5\xi_1 - 5\xi_4 = 0 \\ -2\xi_1 + \xi_2 - 4\xi_3 = 0 \\ -3\xi_1 + 4\xi_3 = 0 \\ -\xi_2 + 5\xi_4 = 0 \end{cases} 5ξ15ξ4=02ξ1+ξ24ξ3=03ξ1+4ξ3=0ξ2+5ξ4=0
ξ1=4/31\xi_1=4/31ξ1=4/31(任意正数),解得:正向量ξ=(4/31,20/31,3/31,4/31)T\boldsymbol{\xi} = (4/31,20/31,3/31,4/31)^Tξ=(4/31,20/31,3/31,4/31)T,验证:1/31⋅(4,20,3,4)⋅L=(0,0,0,0)1/31 \cdot (4,20,3,4) \cdot L = (0,0,0,0)1/31(4,20,3,4)L=(0,0,0,0)成立。
性质 3:正定对角矩阵Ξ=1/31⋅diag(4,20,3,4)\Xi = 1/31\cdot \text{diag}(4,20,3,4)Ξ=1/31diag(4,20,3,4)对角元均为正数,故正定。L~=12(ΞL+LTΞ)\widetilde{L} = \frac{1}{2}(\Xi L + L^T \Xi)L =21(ΞL+LTΞ)是对称矩阵,且行 / 列和为 0(1)计算ΞL\Xi LΞLLTΞL^T \XiLTΞ, ΞL=131⋅(20−8−1200200−200−12120−200020),LTΞ=131⋅(2000−20−820−120−1201200−20020)\Xi L =\frac{1}{31}\cdot \begin{pmatrix} 20 & -8 & -12 & 0 \\ 0 & 20 & 0 & -20 \\ 0 & -12 & 12 & 0 \\ -20 & 0 & 0 & 20 \end{pmatrix}, \quad L^T \Xi =\frac{1}{31}\cdot \begin{pmatrix} 20 & 0 & 0 & -20 \\ -8 & 20 & -12 & 0 \\ -12 & 0 & 12 & 0 \\ 0 & -20 & 0 & 20 \end{pmatrix}ΞL=311 200020820120120120020020 ,LTΞ=311 208120020020012120200020
(2)对称矩阵L~\widetilde{L}L L~=12(ΞL+LTΞ)=162(20−4−6−10−420−6−10−6−6120−10−10020)\widetilde{L} = \frac{1}{2}(\Xi L + L^T \Xi) = \frac{1}{62}\begin{pmatrix} 20 & -4 & -6 & -10 \\ -4 & 20 & -6 & -10 \\ -6 & -6 & 12 & 0 \\ -10 & -10 & 0 & 20 \end{pmatrix}L =21(ΞL+LTΞ)=621 204610420610661201010020
对称成立且行 / 列和为 0。

未完待续

Logo

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

更多推荐