2025-11-18T07:58:12.738440

Graph Signal Wiener Filtering in the Linear Canonical Domain: Theory and Method Design

Cheng, Zhang
The graph linear canonical transform (GLCT)-based filtering methods often optimize transform parameters and filters separately, which results in high computational costs and limited stability. To address this issue, this paper proposes a trainable joint optimization framework that combines GLCT parameters and Wiener filtering into an end-to-end learning process, allowing for synergistic optimization between transform domain construction and filtering operations. The proposed method not only eliminates the cumbersome grid search required by traditional strategies but also significantly enhances the flexibility and training stability of the filtering system. Experimental results on real-world graph data show the proposed method outperforms existing methods in denoising tasks, featuring superior denoising performance, higher robustness and lower computational complexity.
academic

Graph Signal Wiener Filtering in the Linear Canonical Domain: Theory and Method Design

基本信息

  • 论文ID: 2510.10512
  • 标题: Graph Signal Wiener Filtering in the Linear Canonical Domain: Theory and Method Design
  • 作者: Xiaopeng Cheng, Zhichao Zhang
  • 分类: eess.SP (Signal Processing)
  • 发表时间: 2025年10月12日
  • 论文链接: https://arxiv.org/abs/2510.10512

摘要

基于图线性正则变换(GLCT)的滤波方法通常分别优化变换参数和滤波器,这导致了高计算成本和有限的稳定性。为解决这一问题,本文提出了一个可训练的联合优化框架,将GLCT参数和Wiener滤波结合到端到端学习过程中,实现变换域构建和滤波操作之间的协同优化。该方法不仅消除了传统策略所需的繁琐网格搜索,还显著增强了滤波系统的灵活性和训练稳定性。在真实图数据上的实验结果表明,所提方法在去噪任务中优于现有方法,具有更优的去噪性能、更高的鲁棒性和更低的计算复杂度。

研究背景与动机

问题定义

在社交网络、交通系统、生物分子网络等不规则结构中,数据往往位于非欧几里得网格上,使得经典信号处理方法不再适用。图信号处理(GSP)应运而生,将不规则结构数据建模为图,其中节点表示数据实体,边编码它们的关系,信号值附加到节点上。

核心挑战

  1. 噪声干扰:图信号在获取、传输和存储过程中不可避免地受到噪声干扰
  2. 滤波理论适应性:经典线性滤波基于欧几里得空间性质,难以直接迁移到表示非欧几里得空间的图结构
  3. 参数优化复杂性:现有GLCT滤波方法通常分别优化变换参数和滤波器,导致高计算成本和有限稳定性

研究动机

现有方法的局限性主要体现在:

  • 计算效率低:网格搜索策略需要枚举大量候选参数组合
  • 优化不协调:变换域选择和滤波设计的分离导致次优结果
  • 扩展性差:高维参数空间中的网格搜索面临组合爆炸问题

核心贡献

  1. 新的GLCT定义:提出了基于拉普拉斯特征基的CM-CC-CM-GLCT,填补了现有CM-CC-CM-GLCT的空白,并组织了CDDHFs-GLCT和CM-CC-CM-GLCT框架
  2. 可微性理论:证明了GLCT核心模块在加权邻接矩阵和拉普拉斯矩阵下的可微性,为变换参数和滤波器系数的端到端优化提供理论支撑
  3. 联合优化框架:构建了GLCT-GWF框架,实现GLCT参数和滤波器系数的端到端联合优化,在真实图信号去噪任务中验证了其有效性和鲁棒性

方法详解

任务定义

给定观测模型:f~=Gf+n\tilde{f} = Gf + n,其中GG是已知扰动矩阵,ff是光滑信号,nn是加性噪声项。目标是设计最优滤波方法在变换谱域中以最小均方误差(MSE)恢复原始信号ff

核心技术组件

1. 图线性正则变换(GLCT)

GLCT由2×2矩阵M=(a,b;c,d)M = (a,b;c,d)确定,其中adbc=1ad-bc=1。该矩阵可分解为: [abcd]=[10ξ11][1b01][10ξ31]\begin{bmatrix} a & b \\ c & d \end{bmatrix} = \begin{bmatrix} 1 & 0 \\ \xi_1 & 1 \end{bmatrix} \begin{bmatrix} 1 & b \\ 0 & 1 \end{bmatrix} \begin{bmatrix} 1 & 0 \\ \xi_3 & 1 \end{bmatrix}

其中ξ1=d1b\xi_1 = \frac{d-1}{b}ξ2=b\xi_2 = -bξ3=a1b\xi_3 = \frac{a-1}{b}

2. Lap-CM-CC-CM-GLCT定义

本文提出的基于拉普拉斯矩阵的CM-CC-CM-GLCT定义为: f^MIV=FLMf=CMLξ1ULCMLξ2UL1CMLξ3f\hat{f}^{IV}_M = F^M_L f = CM^{\xi_1}_L U_L CM^{\xi_2}_L U^{-1}_L CM^{\xi_3}_L f

3. 联合优化目标

传统方法的两步优化过程:

  1. 固定变换参数(a,b,d)(a,b,d),求解最优滤波器HH^*
  2. 网格搜索(a,b,d)(a,b,d)最小化MSE

本文提出的联合优化: mina,b,d,HE{FLM1HFLMf~f22}\min_{a,b,d,H} E\{\|F^{M^{-1}}_L HF^M_L \tilde{f} - f\|^2_2\}

算法框架

Algorithm 1: 传统网格搜索方法

输入: 图信号f,目标信号f̃,参数网格A,B,D
输出: 最优参数(a*,b*,d*),最优滤波器H*
1. 预计算谱基(特征分解)
2. for a ∈ A, b ∈ B, d ∈ D:
   - 构造GLCT算子F^M和F^{M^{-1}}
   - 求解Wiener-Hopf方程: h = T^{-1}q
   - 评估损失MSE(H,a,b,d)
   - 更新最优解

Algorithm 2: Adam联合优化方法

输入: 图信号f,目标信号f̃,学习率ε
输出: 学习的滤波器和GLCT参数
1. 初始化参数(a₀,b₀,d₀)和滤波器H₀
2. while 未满足停止条件:
   - 计算梯度∇_{H,a,b,d}MSE
   - 更新参数: H ← H - ε∇_H MSE
   - 更新GLCT参数: a,b,d ← a,b,d - ε∇_{a,b,d}MSE

技术创新点

  1. 可微性保证:证明了损失函数关于变换参数和滤波系数的可微性,使得端到端优化成为可能
  2. 计算复杂度优化
    • 网格搜索复杂度:O(nanbndN4)O(n_a n_b n_d N^4)
    • Adam联合优化复杂度:O(KN2)O(KN^2)
  3. 理论性质:新提出的Lap-CM-CC-CM-GLCT满足线性性、零旋转、可加性、可逆性和酉性等重要性质

实验设置

数据集

合成数据

  1. 5-nn图:15个节点,每个节点连接5个最近邻
  2. 瑞士卷图:30个节点的流形结构
  3. 传感器图:20个节点的传感器网络

真实数据

  1. 海表温度(SST):50个节点,k∈{2,6,10}
  2. PM2.5:空气质量数据,50个节点
  3. COVID-19:疫情传播数据,50个节点

评价指标

使用均方误差(MSE)作为主要评价指标,评估去噪性能。

对比方法

  • GFRFT_W:基于加权邻接矩阵的图分数傅里叶变换
  • GFRFT_L:基于拉普拉斯矩阵的图分数傅里叶变换
  • 各种GLCT变体:wAdj-CDDHFs-GLCT、Lap-CDDHFs-GLCT等

实现细节

  • 网格搜索:参数范围0,2,步长0.1
  • Adam优化:学习率0.005,迭代5000次
  • 噪声设置:高斯噪声,标准差s∈{0.5,1.0,1.5}(合成数据),s∈{0.5,0.6,0.7}(真实数据)

实验结果

主要结果

合成数据结果

在三个合成图上的实验表明GLCT方法consistently优于GFRFT方法:

方法5-nn (s=0.5)Swiss Roll (s=0.5)Sensor (s=0.5)
GFRFT_W2.6335.4133.383
GFRFT_L2.8465.2923.432
wAdj-CDDHFs-GLCT2.5375.1163.066

真实数据结果

在SST数据集上,wAdj-CDDHFs-GLCT在k=2, s=0.5设置下达到最低MSE值1.442,相比传统GFRFT方法提升约25%。

性能分析

  1. 去噪性能:GLCT方法在所有测试条件下均表现出优异的去噪性能
  2. 鲁棒性:随着噪声强度增加,GLCT方法的性能退化更加优雅
  3. 计算效率:Adam联合优化显著降低了计算复杂度

消融实验

通过对比不同GLCT变体,验证了各个组件的重要性:

  • CDDHFs-GLCT在多数情况下表现最佳
  • CM-CC-CM-GLCT在某些特定场景下具有优势
  • 拉普拉斯基和邻接矩阵基各有适用场景

相关工作

图信号处理基础

  • 图傅里叶变换(GFT):将经典傅里叶变换扩展到图结构数据
  • 图分数傅里叶变换(GFRFT):引入分数阶参数,在恒等变换和完整GFT之间插值

线性正则变换发展

  • 经典LCT:将旋转推广为仿射变换,扩展FT和FRFT
  • 图域扩展:Zhang等人基于CDDHFs定义GLCT,Li等人提出CM-CC-CM实现

Wiener滤波扩展

传统图Wiener滤波在固定变换域中操作,本文扩展到可学习的变换域选择。

结论与讨论

主要结论

  1. 联合优化框架有效解决了传统方法的计算复杂度和稳定性问题
  2. GLCT的丰富参数化结合联合学习实现了更好的信号-噪声分离
  3. 端到端优化避免了分步优化可能导致的次优结果

局限性

  1. 非凸优化:联合优化问题本质上是非凸的,可能存在局部最优
  2. 参数初始化敏感性:Adam优化的性能可能依赖于初始化策略
  3. 理论收敛保证:缺乏严格的全局收敛性理论分析

未来方向

  1. 发展更强的理论收敛保证
  2. 探索自适应参数初始化策略
  3. 扩展到时变图和动态图信号处理

深度评价

优点

  1. 理论贡献扎实:提供了完整的可微性证明和性质分析
  2. 方法创新性强:首次实现GLCT参数和滤波器的联合优化
  3. 实验验证充分:在多个数据集和设置下验证了方法有效性
  4. 计算效率显著提升:从O(N4)O(N^4)降至O(N2)O(N^2)的复杂度改进

不足

  1. 收敛性分析不足:缺乏对非凸优化收敛性的深入理论分析
  2. 参数敏感性:对学习率和初始化的敏感性分析有限
  3. 应用场景限制:主要关注去噪任务,其他图信号处理任务的适用性待验证

影响力

  1. 学术价值:为图信号处理领域提供了新的优化范式
  2. 实用价值:显著降低的计算复杂度使大规模应用成为可能
  3. 可复现性:提供了详细的算法描述和实验设置

适用场景

  • 大规模图信号去噪任务
  • 需要实时处理的图信号应用
  • 对计算效率有严格要求的嵌入式系统
  • 多尺度图结构分析

参考文献

论文引用了49篇相关文献,涵盖了图信号处理、线性正则变换、Wiener滤波等核心领域的重要工作,为研究提供了坚实的理论基础。


总体评价:本文在图信号处理领域做出了重要贡献,通过联合优化框架有效解决了传统方法的计算复杂度问题,理论分析充分,实验验证全面,具有较高的学术价值和实用价值。