Skip to content

图 Laplacian

Graph Laplacian · Laplacian matrix

用度矩阵减邻接矩阵得到的算子,以二次型衡量相邻顶点取值的不一致。

形式陈述

G=(V,E,w) 是有限无向图,每条边 {u,v} 有非负权重 wuv;无权图取 wuv=1。加权邻接矩阵 A 满足 Auv=wuv,度矩阵 D 是对角矩阵,Duu=du=vwuv。图的组合 Laplacian 定义为

L=DA.

对任意 xRV,对称性给出

xTLx={u,v}Ewuv(xuxv)20.

因此 L 是实对称半正定矩阵。每行元素之和为零,所以 L1=0;更一般地,Lx=0 当且仅当 x 在每个连通分量上为常数。由此再进入特征值分析,便知零特征值的重数恰等于 G 的连通分量数。

若图没有孤立点,常用归一化 Laplacian

L=D1/2LD1/2=ID1/2AD1/2.

有孤立点时需声明约定;本页令 Duu1/2=0,因而对应行列为零。

直觉

xu 看作顶点 u 上的温度或电势。每条边都惩罚两端数值的差,xTLx 汇总了全图的不协调能量。若能量为零,正权边两端必须取相同值;沿路径传播后,每个连通分量只能保持一个常数。Laplacian 因而把图的连通结构转写成矩阵的核。

Lx 在顶点 u 的分量是 duxuvwuvxv,也就是本点值与邻居加权值之间的差。它与连续空间中的 Laplace 算子承担相似角色:检测局部偏离平衡的程度,但图上没有坐标或微分,全部信息由边给出。

例子与边界

三顶点路径 P3 的 Laplacian 为

L=(110121011),

其特征值是 0,1,3。唯一零特征方向由常向量张成,反映图连通。完全图 Kn 的 Laplacian 是 nIJ,谱为 0 一次、n 重数 n1;高对称性把所有与常向量正交的振动方向赋予相同能量。若图由两个互不相连的连通图组成,分别取两个分量的示性向量便得到两个线性无关的零特征向量。

这些结论不能原样套到有向图:邻接矩阵不再对称,DA 可能不是对称半正定矩阵。允许负边权时,上式会成为带符号的平方项之和,二次型可能取负值,因而失去半正定性。组合 Laplacian、归一化 Laplacian 与随机游走 Laplacian ID1A 还具有不同的特征值尺度;比较谱隙时必须先固定所用规范。

推论与应用

由核的刻画可立即得到:有限无向图连通,当且仅当 0L 的单特征值。第二小特征值因此衡量离断开还有多远,并把连通性与割的瓶颈联系起来。Rayleigh 商

xTLxxTx

把寻找低能量顶点函数转为特征值问题,是谱划分和扩张估计的入口。

Laplacian 的任一主余子式还出现在矩阵树定理中,用行列式计数生成树;其伪逆编码电网络的有效电阻。随机游走、扩散、谱聚类和图信号平滑看似来自不同应用,实质上都在使用同一个“邻边差异能量”。

参考资料
  • Daniel A. Spielman, Spectral and Algebraic Graph Theory, Yale lecture notes, Chapters 2–4.
  • Fan R. K. Chung, Spectral Graph Theory, American Mathematical Society, 1997, Chapters 1–2.