图信号处理——拉普拉斯矩阵

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

G(V,E)G(V,E)G(V,E) 表示图,VVV 为图的节点集,EEE 为图的边集。

A=[aij]A = [a_{ij}]A=[aij] 表示图的邻接矩阵。图的拉普拉斯矩阵为:L=D−AL = D-AL=DA,其中 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}L1=0,暗示含有特征值 0;
  • 对称矩阵,有N个实特征值和特征向量;
  • 对于任意图信号xxx(Lx)i=∑j∈Ni(xi−xj)(Lx)_i = \sum_{j\in N_i} (x_i - x_j)(Lx)i=jNi(xixj),即 LxLxLx 的第 iii 个元素的含义为节点 iii 的信号与其邻居节点的信号的差距之和,因此称其为 variation
  • 对于任意向量 x∈RNx\in\mathbb{R}^NxRNxTLx=∑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=ixijNi(xixj)=ijNixi(xixj)=i,j(xixj)20 。说明 LLL 是半正定矩阵。xTLxx^TLxxTLx 称为 total variation,常作为 Graph Laplacian Regularizer 用于正则化.
  • 拉普拉斯矩阵有几种归一化方式:L^=I−D−1A\hat{L} = I-D^{-1}AL^=ID1A或者Lˉ=I−D−12AD−12\bar{L} = I-D^{-\frac{1}{2}}AD^{-\frac{1}{2}}Lˉ=ID21AD21其中 L^\hat{L}L^ 保证每一行元素之和为 1,Lˉ\bar{L}Lˉ 为对称矩阵。

图信号滤波

假设向量 x⃗∈RN\vec{x}\in\mathbb{R}^NxRN 为原始图信号,其中 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=jNixj,即新的图信号 yyy 表示邻居节点的原始信号之和。易得 D−1AD^{-1}AD1A 表示所有邻居节点的原始信号的均值。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

颹蕭蕭

白嫖?

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值