用 G(V,E)G(V,E)G(V,E) 表示图,VVV 为图的节点集,EEE 为图的边集。
用 A=[aij]A = [a_{ij}]A=[aij] 表示图的邻接矩阵。图的拉普拉斯矩阵为:L=D−AL = D-AL=D−A,其中 D=[dij]D = [d_{ij}]D=[dij],为对角矩阵,对角元满足 dii=∑jaijd_{ii} = \sum_j a_{ij}dii=∑jaij,即 AAA 的第iii行元素之和。
拉普拉斯矩阵有几个重要性质:
- L⋅1⃗=0⃗L\cdot \vec{1} = \vec{0}L⋅1=0,暗示含有特征值 0;
- 对称矩阵,有N个实特征值和特征向量;
- 对于任意图信号xxx,(Lx)i=∑j∈Ni(xi−xj)(Lx)_i = \sum_{j\in N_i} (x_i - x_j)(Lx)i=∑j∈Ni(xi−xj),即 LxLxLx 的第 iii 个元素的含义为节点 iii 的信号与其邻居节点的信号的差距之和,因此称其为 variation;
- 对于任意向量 x∈RNx\in\mathbb{R}^Nx∈RN,xTLx=∑ixi∑j∈Ni(xi−xj)=∑i∑j∈Nixi(xi−xj)=∑i,j(xi−xj)2≥0x^TLx = \sum_i x_i\sum_{j\in N_i} (x_i - x_j)= \sum_i \sum_{j\in N_i} x_i(x_i - x_j)= \sum_{i,j}(x_i - x_j)^2\geq 0xTLx=∑ixi∑j∈Ni(xi−xj)=∑i∑j∈Nixi(xi−xj)=∑i,j(xi−xj)2≥0 。说明 LLL 是半正定矩阵。xTLxx^TLxxTLx 称为 total variation,常作为 Graph Laplacian Regularizer 用于正则化.
- 拉普拉斯矩阵有几种归一化方式:L^=I−D−1A\hat{L} = I-D^{-1}AL^=I−D−1A或者Lˉ=I−D−12AD−12\bar{L} = I-D^{-\frac{1}{2}}AD^{-\frac{1}{2}}Lˉ=I−D−21AD−21其中 L^\hat{L}L^ 保证每一行元素之和为 1,Lˉ\bar{L}Lˉ 为对称矩阵。
图信号滤波
假设向量 x⃗∈RN\vec{x}\in\mathbb{R}^Nx∈RN 为原始图信号,其中 NNN 为图的节点数,我们可以用 A,W,D,L,L1,L2A,W,D,L,L1,L2A,W,D,L,L1,L2 对原始图信号滤波。
具体来看:y=Axy = Axy=Ax,则 yi=∑j∈Nixjy_i = \sum_{j\in N_i} x_jyi=∑j∈Nixj,即新的图信号 yyy 表示邻居节点的原始信号之和。易得 D−1AD^{-1}AD−1A 表示所有邻居节点的原始信号的均值。
图拉普拉斯矩阵详解

本文深入探讨了图拉普拉斯矩阵的定义及其在图信号处理中的应用,包括邻接矩阵、对角矩阵的构成,以及拉普拉斯矩阵的性质如特征值、特征向量、半正定性和图信号滤波等。
174

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



