图神经网络中的图傅里叶变换与图信号处理原理详解
字数 1621 2025-11-29 18:56:49
图神经网络中的图傅里叶变换与图信号处理原理详解
图傅里叶变换是图信号处理的核心工具,它将信号从图空间变换到谱域,从而便于分析图的频域特性。以下是逐步讲解:
1. 图信号与拉普拉斯矩阵
- 图信号定义:图信号是一个函数 \(f: V \to \mathbb{R}\),其中每个节点 \(v_i\) 对应一个标量值 \(f(i)\)。
- 拉普拉斯矩阵:
- 度矩阵 \(D\):对角矩阵,\(D_{ii} = \deg(v_i)\)。
- 邻接矩阵 \(A\):表示节点连接关系。
- 拉普拉斯矩阵 \(L = D - A\)。
- 归一化拉普拉斯矩阵 \(L_{\text{norm}} = I - D^{-1/2} A D^{-1/2}\)(常用形式)。
关键性质:
- \(L\) 是实对称半正定矩阵,其特征值 \(\lambda_k\) 非负,特征向量 \(u_k\) 构成正交基。
- 特征值 \(\lambda_k\) 表示频率:小特征值对应低频(平滑信号),大特征值对应高频(剧烈变化信号)。
2. 图傅里叶变换
- 传统傅里叶变换:将信号分解为正弦波(频域基)。
- 图傅里叶变换:将图信号投影到拉普拉斯矩阵的特征向量基上:
\[ \hat{f}(\lambda_k) = \langle f, u_k \rangle = \sum_{i=1}^n f(i) u_k(i), \]
其中 \(\hat{f}\) 是谱域信号,\(u_k\) 是第 \(k\) 个特征向量。
- 逆变换:
\[ f(i) = \sum_{k=1}^n \hat{f}(\lambda_k) u_k(i). \]
物理意义:
- 特征向量 \(u_k\) 是图的“振动模式”,特征值 \(\lambda_k\) 表示振动频率。
- 低频分量对应图中缓慢变化的信号(如社区结构),高频分量对应局部剧烈变化(如噪声)。
3. 图滤波与卷积
- 谱域滤波:在谱域对信号进行滤波操作:
\[ \hat{f}_{\text{out}}(\lambda_k) = g(\lambda_k) \cdot \hat{f}_{\text{in}}(\lambda_k), \]
其中 \(g(\cdot)\) 是滤波器函数(如低通、高通滤波器)。
- 空域卷积:通过逆变换将滤波后的信号映射回空域:
\[ f_{\text{out}} = U \cdot \text{diag}(g(\lambda_k)) \cdot U^T f_{\text{in}}. \]
挑战与改进:
- 直接计算特征分解复杂度高(\(O(n^3)\))。
- 切比雪夫多项式近似:用多项式逼近滤波器 \(g(\lambda)\),避免显式特征分解。
- GCN的简化:使用一阶近似 \(g(\lambda) \approx \theta(1 - \lambda)\),得到空域消息传递公式。
4. 应用与实例
- 信号去噪:通过低通滤波保留平滑信号,抑制高频噪声。
- 社区检测:低频特征向量(如Fiedler向量)用于划分社区。
- 图卷积网络(GCN):利用谱域滤波学习节点表示,如:
\[ H^{(l+1)} = \sigma\left( \tilde{D}^{-1/2} \tilde{A} \tilde{D}^{-1/2} H^{(l)} W^{(l)} \right), \]
其中 \(\tilde{A} = A + I\) 是带自环的邻接矩阵,\(\tilde{D}\) 是对应的度矩阵。
5. 总结
- 图傅里叶变换将图信号从空域转换到谱域,便于频域分析。
- 拉普拉斯矩阵的特征基提供了图的频域表示,低频与全局结构相关,高频与局部细节相关。
- 通过谱域滤波可实现图卷积,为GNN提供了理论基础。