图卷积神经网络之 - 图傅立叶变换的解释

L = np.array([
    [ 2, -1,  0, -1],
    [-1,  2, -1,  0],
    [ 0, -1,  2, -1],
    [-1,  0, -1,  2]
])

上述4个节点的简单无向图进行图傅立叶变换：

1. 我们计算了图拉普拉斯矩阵 $L$ 的特征值和特征向量。特征值是：

   $\lambda =\left[\begin{array}{cccc}0,& 2,& 2,& 4\end{array}\right]$

   特征向量 $U$（归一化后）是：

   $U=\left[\begin{array}{cccc}-0.5& 0& 0.707& 0.5\\ -0.5& -0.707& 0& -0.5\\ -0.5& 0& -0.707& 0.5\\ -0.5& 0.707& 0& -0.5\end{array}\right]$

2. 我们定义了一个信号 $x={\left[\begin{array}{cccc}1& 2& 3& 4\end{array}\right]}^T$，这个信号在图的节点上有不同的值。

3. 图傅立叶变换 $\hat{x}=U^{T}x$ 将这个信号转换到频域，结果是：

   $\hat{x}=\left[\begin{array}{cccc}-5,& 1.414,& -1.414,& -1\end{array}\right]$

4. 图傅立叶逆变换 $x=U\hat{x}$ 将变换后的信号 $\hat{x}$ 转换回空间域，结果是原始信号 $x$：

   $x={\left[\begin{array}{cccc}1,& 2,& 3,& 4\end{array}\right]}^T$

这表明图傅立叶变换及其逆变换可以无损地在图的空间域和频域之间转换信号。在这个过程中，特征向量定义了图的频域，而特征值类似于传统傅立叶变换中的频率。

from scipy.linalg import eigh

# 使用前面计算的拉普拉斯矩阵L
L = np.array([
    [ 2, -1,  0, -1],
    [-1,  2, -1,  0],
    [ 0, -1,  2, -1],
    [-1,  0, -1,  2]
])

# 计算特征值和特征向量
eigenvalues, eigenvectors = eigh(L)

# 假设有一个信号x，定义在图的节点上
x = np.array([1, 2, 3, 4])

# 图傅立叶变换
x_hat = eigenvectors.T @ x

# 图傅立叶逆变换
x_reconstructed = eigenvectors @ x_hat

# 显示特征值和特征向量，以及变换后的信号和重构的信号
eigenvalues, eigenvectors, x_hat, x_reconstructed

执行结果：

(array([2.66453526e-15, 2.00000000e+00, 2.00000000e+00, 4.00000000e+00]),

 array([[-5.00000000e-01,  0.00000000e+00,  7.07106781e-01, 5.00000000e-01],
        [-5.00000000e-01, -7.07106781e-01, -4.57558474e-16, -5.00000000e-01],
        [-5.00000000e-01,  2.87079964e-16, -7.07106781e-01, 5.00000000e-01],
        [-5.00000000e-01,  7.07106781e-01, -1.70478510e-16, -5.00000000e-01]]),

 array([-5.        ,  1.41421356, -1.41421356, -1.        ]),
 
 array([1., 2., 3., 4.]))