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

線形正準領域におけるグラフ信号ウィーナフィルタリング：理論と方法設計

基本情報

論文ID: 2510.10512
タイトル: Graph Signal Wiener Filtering in the Linear Canonical Domain: Theory and Method Design
著者: Xiaopeng Cheng, Zhichao Zhang
分類: eess.SP（信号処理）
発表日: 2025年10月12日
論文リンク: https://arxiv.org/abs/2510.10512

要約

グラフ線形正準変換（GLCT）に基づくフィルタリング方法は、従来、変換パラメータとフィルタを個別に最適化していたため、計算コストが高く安定性が限定的でした。この問題を解決するため、本論文はGLCTパラメータとウィーナフィルタを端末間学習プロセスに統合した訓練可能な共同最適化フレームワークを提案します。これにより、変換領域の構築とフィルタリング操作間の協調最適化を実現します。本手法は従来の戦略に必要だった煩雑なグリッド探索を排除するだけでなく、フィルタリングシステムの柔軟性と訓練安定性を大幅に向上させます。実グラフデータ上の実験結果は、提案手法がノイズ除去タスクにおいて既存手法を上回り、優れたノイズ除去性能、高いロバスト性、低い計算複雑度を実現することを示しています。

研究背景と動機

問題定義

ソーシャルネットワーク、交通システム、生物分子ネットワークなどの不規則な構造では、データはしばしば非ユークリッド格子上に位置し、古典的信号処理手法は適用できなくなります。グラフ信号処理（GSP）はこの課題に対応し、不規則な構造データをグラフとしてモデル化します。ここでノードはデータエンティティを表し、エッジはそれらの関係を符号化し、信号値はノードに付加されます。

中核的課題

ノイズ干渉：グラフ信号は取得、伝送、保存プロセス中に不可避的にノイズの影響を受けます
フィルタリング理論の適応性：古典的線形フィルタリングはユークリッド空間の性質に基づいており、非ユークリッド空間を表現するグラフ構造への直接的な転用は困難です
パラメータ最適化の複雑性：既存のGLCTフィルタリング手法は通常、変換パラメータとフィルタを個別に最適化するため、計算コストが高く安定性が限定的です

研究動機

既存手法の限界は主に以下の点に表れます：

計算効率の低さ：グリッド探索戦略は多数の候補パラメータ組み合わせの列挙が必要です
最適化の不調和：変換領域の選択とフィルタ設計の分離は次善的な結果をもたらします
拡張性の欠如：高次元パラメータ空間におけるグリッド探索は組み合わせ爆発に直面します

中核的貢献

新しいGLCT定義：ラプラシアン固有基に基づくCM-CC-CM-GLCTを提案し、既存のCM-CC-CM-GLCTの空白を埋め、CDDHFs-GLCTおよびCM-CC-CM-GLCTフレームワークを組織化しました
微分可能性理論：重み付き隣接行列およびラプラシアン行列下でのGLCT中核モジュールの微分可能性を証明し、変換パラメータとフィルタ係数の端末間最適化に理論的支援を提供します
共同最適化フレームワーク：GLCT-GWFフレームワークを構築し、GLCTパラメータとフィルタ係数の端末間共同最適化を実現し、実グラフ信号ノイズ除去タスクでその有効性とロバスト性を検証しました

方法の詳細説明

タスク定義

観測モデル $\tilde{f} = Gf + n$ が与えられます。ここで $G$ は既知の摂動行列、 $f$ は平滑信号、 $n$ は加法ノイズ項です。目標は、変換スペクトル領域で最小二乗平均誤差（MSE）を最小化して元の信号 $f$ を復元する最適フィルタリング手法を設計することです。

中核技術コンポーネント

1. グラフ線形正準変換（GLCT）

GLCTは2×2行列 $M = (a,b;c,d)$ で決定されます。ここで $ad-bc=1$ です。この行列は以下のように分解できます： $\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}$

ここで $\xi_1 = \frac{d-1}{b}$ 、 $\xi_2 = -b$ 、 $\xi_3 = \frac{a-1}{b}$ です。

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

本論文で提案されたラプラシアン行列に基づくCM-CC-CM-GLCT定義は以下の通りです： $\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. 共同最適化目標

従来の手法の2段階最適化プロセス：

変換パラメータ $(a,b,d)$ を固定し、最適フィルタ $H^*$ を求解
$(a,b,d)$ をグリッド探索してMSEを最小化

本論文で提案される共同最適化： $\min_{a,b,d,H} E\{\|F^{M^{-1}}_L HF^M_L \tilde{f} - f\|^2_2\}$

アルゴリズムフレームワーク

アルゴリズム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}}を構築
   - ウィーナ・ホップ方程式を求解: h = T^{-1}q
   - 損失MSE(H,a,b,d)を評価
   - 最適解を更新

アルゴリズム2：Adam共同最適化手法

入力: グラフ信号f、目標信号f̃、学習率ε
出力: 学習されたフィルタおよびGLCTパラメータ
1. パラメータ(a₀,b₀,d₀)およびフィルタH₀を初期化
2. 停止条件を満たすまで:
   - 勾配∇_{H,a,b,d}MSEを計算
   - パラメータを更新: H ← H - ε∇_H MSE
   - GLCTパラメータを更新: a,b,d ← a,b,d - ε∇_{a,b,d}MSE