[论文解读] Unitary Shift Operators on a Graph

Bruno Scalzo Dees, Ljubiša Stanković+3 more|arXiv (Cornell University)|Sep 12, 2019

advanced mathematical theories参考文献 12被引用 5

一句话总结

本文提出一种幺正图移位算子（GSO），在正向与反向图移位过程中保持信号能量不变，从而实现具有正交基底和自然频域解释的严格图离散傅里叶变换（GDFT）。该方法推导出图微分算子，并利用GSO定义一致的图解析信号与图离散希尔伯特变换（GDHT），并通过具有能量保持和等距性质的有向图上的直观示例加以验证。

ABSTRACT

A unitary shift operator (GSO) for signals on a graph is introduced, which exhibits the desired property of energy preservation over both backward and forward graph shifts. For rigour, the graph differential operator is also derived in an analytical form. The commutativity relation of the shift operator with the Fourier transform is next explored in conjunction with the proposed GSO to introduce a graph discrete Fourier transform (GDFT) which, unlike existing approaches, ensures the orthogonality of GDFT bases and admits a natural frequency-domain interpretation. The proposed GDFT is shown to allow for a coherent definition of the graph discrete Hilbert transform (GDHT) and the graph analytic signal. The advantages of the proposed GSO are demonstrated through illustrative examples.

研究动机与目标

• 解决图信号处理（GSP）中缺乏能量保持、等距图移位算子的问题，确保移位过程中信号范数不变。
• 开发具有正交基底和清晰频域解释的图傅里叶变换（GDFT），不同于现有非正交方法。
• 基于所提出的幺正GSO，实现图离散希尔伯特变换（GDHT）与图解析信号的统一定义。
• 推导与幺正移位框架一致的图微分算子。
• 通过有向图结构上的直观示例，展示该框架的实际应用价值。

提出的方法

• 提出幺正移位算子 S，满足 S⁻¹ = Sᵀ 且 ∥Sx∥ = ∥x∥，确保能量保持与可逆性。
• 通过将邻接矩阵的特征值映射到单位圆上，构建GSO，确保等距性的同时保持邻域连通性。
• 通过解析方式推导图微分算子为 ∇ϕ = ω，实现在图上的局部相位与频率估计。
• 通过幺正GSO的特征分解定义图离散傅里叶变换（GDFT），获得正交谐波基底。
• 利用GDFT系数定义图解析信号，其中 xₐ(n) = a(n)e^{jϕ(n)}，分别表示局部幅值 a(n) 与相位 ϕ(n)。
• 通过非负GDFT系数的频域处理建立GDHT，实现局部幅值、相位与频率等特征提取。

实验结果

研究问题

• RQ1 能否定义一种幺正图移位算子，使其在正向与反向移位过程中均保持信号能量，确保等距性？
• RQ2 所提出的GSO是否能实现具有正交基底和自然频域解释的图离散傅里叶变换（GDFT）？
• RQ3 能否基于所提出的GSO框架，严格定义图离散希尔伯特变换（GDHT）与图解析信号？
• RQ4 所推导的图微分算子与图信号中的局部相位和频率有何关系？
• RQ5 该框架在揭示不规则图域上隐藏信号特征（如局部幅值、相位与频率）方面具有何种实际效用？

主要发现

• 所提出的幺正GSO确保能量保持，满足 ∥Sx∥ = ∥x∥ 与 ∥S⁻¹x∥ = ∥x∥，正向与反向移位均满足等距性。
• 基于GSO特征分解的GDFT产生正交基底，实现类似经典傅里叶分析的连贯频域解释。
• 图解析信号通过GDFT系数定义，明确给出局部幅值 a(n) 与相位 ϕ(n) 的表达式，以及局部频率 ω(n) = dϕ(n)/dn。
• 图微分算子被解析推导为 ∇ϕ = ω，实现从相位梯度中估计局部频率。
• GDHT成功提取局部特征：幅值、相位与频率，如图3所示，非明显信号特征得以显现。
• 对于含8个顶点的有向图，GSO的特征值为 ω ∈ {0.238, 1.244, 2.209, 2.817} rad/s，对应GDFT中的不同谐波模式。