图神经网络中的图傅里叶变换与图信号处理原理
字数 1966 2025-11-19 17:05:32
图神经网络中的图傅里叶变换与图信号处理原理
描述
图傅里叶变换(Graph Fourier Transform, GFT)是图信号处理的核心工具,它将传统傅里叶变换扩展到图结构数据上,用于分析图信号的频率成分。通过图的拉普拉斯矩阵特征分解,GFT将图信号从顶点域变换到谱域,为图卷积等操作提供理论基础。理解GFT是掌握谱域图神经网络(如GCN、ChebNet)的关键前提。
1. 图信号与顶点域
- 图信号定义:设图 \(G=(V, E)\) 有 \(n\) 个顶点,图信号是函数 \(f: V \to \mathbb{R}\),可表示为向量 \(\mathbf{f} = [f(v_1), f(v_2), ..., f(v_n)]^T\)。例如,社交网络中每个用户的年龄构成一个图信号。
- 顶点域操作:在顶点域中,信号处理依赖邻接关系。例如,平滑操作可通过邻域平均实现:\(\mathbf{f}'(v_i) = \frac{1}{d_i} \sum_{v_j \in N(v_i)} f(v_j)\),其中 \(d_i\) 是顶点度。但顶点域操作缺乏全局频率视角。
2. 图的拉普拉斯矩阵
- 定义:拉普拉斯矩阵 \(L = D - A\),其中 \(D\) 是度矩阵(对角矩阵),\(A\) 是邻接矩阵。
- 归一化拉普拉斯:常用对称归一化形式 \(L_{sym} = D^{-1/2} L D^{-1/2} = I - D^{-1/2} A D^{-1/2}\),其特征值范围在 [0, 2] 内,利于数值稳定。
- 性质:\(L\) 是实对称半正定矩阵,有正交特征向量组 \(\mathbf{u}_1, ..., \mathbf{u}_n\) 和实非负特征值 \(\lambda_1 \leq ... \leq \lambda_n\)。特征值反映图的“频率”,小特征值对应低频(平滑信号),大特征值对应高频(剧烈变化信号)。
3. 图傅里叶变换推导
- 传统傅里叶变换类比:传统傅里叶变换将信号分解为正弦波(拉普拉斯算子的特征函数)。图拉普拉斯矩阵 \(L\) 是图上的拉普拉斯算子,其特征向量 \(\mathbf{u}_k\) 是图的“振动模式”。
- GFT定义:
- 正变换:图信号 \(\mathbf{f}\) 的GFT为 \(\hat{\mathbf{f}} = U^T \mathbf{f}\),其中 \(U = [\mathbf{u}_1, ..., \mathbf{u}_n]\) 是 \(L\) 的特征向量矩阵。分量 \(\hat{f}(\lambda_k) = \mathbf{u}_k^T \mathbf{f}\) 表示信号在频率 \(\lambda_k\) 上的强度。
- 逆变换:\(\mathbf{f} = U \hat{\mathbf{f}}\),即用特征向量线性组合重建原信号。
- 计算示例:对链状图(一维网格),\(L\) 的特征向量是离散余弦函数,GFT退化为离散余弦变换(DCT)。
4. 图滤波与卷积定理
- 顶点域卷积困难:图结构不规则,无法直接定义平移操作,故顶点域卷积无明确定义。
- 谱域卷积:借鉴卷积定理(时域卷积等价于频域乘积),定义图卷积为:
\[ \mathbf{f} * g = U \left( \hat{g}(\Lambda) \circ (U^T \mathbf{f}) \right) \]
其中 \(\hat{g}(\Lambda) = \text{diag}(g(\lambda_1), ..., g(\lambda_n))\) 是谱域滤波器函数。滤波操作在谱域是对角矩阵乘法,计算复杂度高(需特征分解)。
5. 应用与优化
- 图信号去噪:通过低通滤波(保留小特征值对应分量)平滑信号。例如,最小化 \(\|\mathbf{f} - \mathbf{y}\|^2 + \alpha \mathbf{f}^T L \mathbf{f}\),其中 \(\mathbf{y}\) 是噪声信号,正则项 \(\mathbf{f}^T L \mathbf{f}\) 惩罚高频成分。
- 图神经网络简化:为避免昂贵特征分解,ChebNet用切比雪夫多项式逼近滤波器 \(g(\lambda)\),GCN进一步简化为一阶近似,实现高效邻域聚合。
总结
图傅里叶变换通过拉普拉斯矩阵特征分解将图信号映射到谱域,奠定了谱域图卷积的理论基础。尽管直接应用计算成本高,但其思想推动了ChebNet、GCN等高效图神经网络的发展。