Filtrage de Wiener des Signaux Graphiques dans le Domaine Canonique Linéaire : Théorie et Conception de Méthode

Informations Fondamentales

• ID de l'article : 2510.10512
• Titre : Graph Signal Wiener Filtering in the Linear Canonical Domain: Theory and Method Design
• Auteurs : Xiaopeng Cheng, Zhichao Zhang
• Classification : eess.SP (Traitement du Signal)
• Date de publication : 12 octobre 2025
• Lien de l'article : https://arxiv.org/abs/2510.10512

Résumé

Les méthodes de filtrage basées sur la transformation canonique linéaire graphique (GLCT) optimisent généralement séparément les paramètres de transformation et le filtre, ce qui entraîne des coûts de calcul élevés et une stabilité limitée. Pour résoudre ce problème, cet article propose un cadre d'optimisation conjointe entraînable qui intègre les paramètres GLCT et le filtrage de Wiener dans un processus d'apprentissage de bout en bout, réalisant une optimisation synergique entre la construction du domaine de transformation et l'opération de filtrage. Cette méthode non seulement élimine la recherche en grille fastidieuse requise par les stratégies traditionnelles, mais améliore également considérablement la flexibilité et la stabilité d'entraînement du système de filtrage. Les résultats expérimentaux sur des données graphiques réelles démontrent que la méthode proposée surpasse les méthodes existantes dans les tâches de débruitage, avec une meilleure performance de débruitage, une robustesse supérieure et une complexité de calcul réduite.

Contexte de Recherche et Motivation

Définition du Problème

Dans les structures irrégulières telles que les réseaux sociaux, les systèmes de transport et les réseaux de molécules biologiques, les données résident souvent sur des grilles non-euclidiennes, rendant les méthodes classiques de traitement du signal inapplicables. Le traitement des signaux graphiques (GSP) émerge pour modéliser les données de structures irrégulières en tant que graphes, où les nœuds représentent les entités de données, les arêtes codent leurs relations, et les valeurs de signal sont attachées aux nœuds.

Défis Fondamentaux

1. Interférence du bruit : Les signaux graphiques subissent inévitablement des perturbations de bruit lors de l'acquisition, la transmission et le stockage
2. Adaptabilité de la théorie de filtrage : Le filtrage linéaire classique basé sur les propriétés de l'espace euclidien est difficile à transférer directement aux structures graphiques représentant l'espace non-euclidien
3. Complexité de l'optimisation des paramètres : Les méthodes GLCT existantes optimisent généralement séparément les paramètres de transformation et le filtre, entraînant des coûts de calcul élevés et une stabilité limitée

Motivation de la Recherche

Les limitations des méthodes existantes se manifestent principalement par :

• Efficacité de calcul faible : La stratégie de recherche en grille nécessite d'énumérer de nombreuses combinaisons de paramètres candidats
• Optimisation non coordonnée : La séparation entre la sélection du domaine de transformation et la conception du filtrage entraîne des résultats sous-optimaux
• Mauvaise extensibilité : La recherche en grille dans l'espace de paramètres de haute dimension fait face à une explosion combinatoire

Contributions Fondamentales

1. Nouvelle définition GLCT : Propose la CM-CC-CM-GLCT basée sur la base de vecteurs propres laplaciens, comblant les lacunes de la CM-CC-CM-GLCT existante, et organisant le cadre CDDHFs-GLCT et CM-CC-CM-GLCT
2. Théorie de différentiabilité : Démontre la différentiabilité des modules fondamentaux GLCT sous les matrices d'adjacence pondérées et les matrices laplaciennes, fournissant un support théorique pour l'optimisation de bout en bout des paramètres de transformation et des coefficients de filtre
3. Cadre d'optimisation conjointe : Construit le cadre GLCT-GWF, réalisant l'optimisation conjointe de bout en bout des paramètres GLCT et des coefficients de filtre, validant son efficacité et sa robustesse dans les tâches réelles de débruitage de signaux graphiques

Explication Détaillée de la Méthode

Définition de la Tâche

Étant donné le modèle d'observation : $\tilde{f} = Gf + n$, où $G$ est une matrice de perturbation connue, $f$ est un signal lisse, et $n$ est un terme de bruit additif. L'objectif est de concevoir une méthode de filtrage optimale pour récupérer le signal original $f$ dans le domaine spectral de transformation avec une erreur quadratique moyenne (MSE) minimale.

Composants Techniques Fondamentaux

1. Transformation Canonique Linéaire Graphique (GLCT)

La GLCT est déterminée par une matrice 2×2 $M = (a,b;c,d)$ où $ad-bc=1$. Cette matrice peut être décomposée comme : $\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}$

où $\xi_1 = \frac{d-1}{b}$, $\xi_2 = -b$, $\xi_3 = \frac{a-1}{b}$.

2. Définition de Lap-CM-CC-CM-GLCT

La définition proposée de CM-CC-CM-GLCT basée sur la matrice laplacienne est : $\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. Objectif d'Optimisation Conjointe

Le processus d'optimisation en deux étapes des méthodes traditionnelles :

1. Fixer les paramètres de transformation $(a,b,d)$, résoudre le filtre optimal $H^*$
2. Recherche en grille $(a,b,d)$ pour minimiser la MSE

L'optimisation conjointe proposée : $\min_{a,b,d,H} E\{\|F^{M^{-1}}_L HF^M_L \tilde{f} - f\|^2_2\}$

Cadre Algorithmique

Algorithme 1 : Méthode Traditionnelle de Recherche en Grille

Entrée : signal graphique f, signal cible f̃, grille de paramètres A,B,D
Sortie : paramètres optimaux (a*,b*,d*), filtre optimal H*
1. Pré-calculer la base spectrale (décomposition en valeurs propres)
2. pour a ∈ A, b ∈ B, d ∈ D :
   - Construire l'opérateur GLCT F^M et F^{M^{-1}}
   - Résoudre l'équation de Wiener-Hopf : h = T^{-1}q
   - Évaluer la perte MSE(H,a,b,d)
   - Mettre à jour la solution optimale

Algorithme 2 : Méthode d'Optimisation Conjointe Adam

Entrée : signal graphique f, signal cible f̃, taux d'apprentissage ε
Sortie : filtre appris et paramètres GLCT
1. Initialiser les paramètres (a₀,b₀,d₀) et le filtre H₀
2. tant que la condition d'arrêt n'est pas satisfaite :
   - Calculer les gradients ∇_{H,a,b,d}MSE
   - Mettre à jour les paramètres : H ← H - ε∇_H MSE
   - Mettre à jour les paramètres GLCT : a,b,d ← a,b,d - ε∇_{a,b,d}MSE

Points d'Innovation Technique

1. Garantie de différentiabilité : Démontre la différentiabilité de la fonction de perte par rapport aux paramètres de transformation et aux coefficients de filtre, rendant possible l'optimisation de bout en bout
2. Optimisation de la complexité de calcul :
   • Complexité de la recherche en grille : $O(n_a n_b n_d N^4)$
   • Complexité de l'optimisation conjointe Adam : $O(KN^2)$
3. Propriétés théoriques : La Lap-CM-CC-CM-GLCT nouvellement proposée satisfait les propriétés importantes de linéarité, rotation nulle, additivité, inversibilité et unitarité

Configuration Expérimentale

Ensembles de Données

Données Synthétiques

1. Graphe 5-nn : 15 nœuds, chaque nœud connecté à 5 plus proches voisins
2. Graphe Swiss Roll : 30 nœuds avec structure de variété
3. Graphe de capteurs : 20 nœuds de réseau de capteurs

Données Réelles

1. Température de Surface Marine (SST) : 50 nœuds, k∈{2,6,10}
2. PM2.5 : Données de qualité de l'air, 50 nœuds
3. COVID-19 : Données de propagation épidémique, 50 nœuds

Métriques d'Évaluation

L'erreur quadratique moyenne (MSE) est utilisée comme métrique d'évaluation principale pour évaluer la performance de débruitage.

Méthodes de Comparaison

• GFRFT_W : Transformation de Fourier fractionnaire graphique basée sur la matrice d'adjacence pondérée
• GFRFT_L : Transformation de Fourier fractionnaire graphique basée sur la matrice laplacienne
• Variantes GLCT : wAdj-CDDHFs-GLCT, Lap-CDDHFs-GLCT, etc.

Détails d'Implémentation

• Recherche en grille : Plage de paramètres 0,2, pas 0.1
• Optimisation Adam : Taux d'apprentissage 0.005, 5000 itérations
• Configuration du bruit : Bruit gaussien, écart-type s∈{0.5,1.0,1.5} (données synthétiques), s∈{0.5,0.6,0.7} (données réelles)

Résultats Expérimentaux

Résultats Principaux

Résultats sur Données Synthétiques

Les expériences sur trois graphes synthétiques montrent que la méthode GLCT surpasse systématiquement la méthode GFRFT :

Résultats sur Données Réelles

Sur l'ensemble de données SST, wAdj-CDDHFs-GLCT atteint la valeur MSE la plus basse de 1.442 dans la configuration k=2, s=0.5, représentant une amélioration d'environ 25% par rapport aux méthodes GFRFT traditionnelles.

Analyse de Performance

1. Performance de débruitage : La méthode GLCT démontre une performance de débruitage supérieure dans toutes les conditions de test
2. Robustesse : Avec l'augmentation de l'intensité du bruit, la dégradation de performance de la méthode GLCT est plus gracieuse
3. Efficacité de calcul : L'optimisation conjointe Adam réduit considérablement la complexité de calcul

Expériences d'Ablation

En comparant différentes variantes GLCT, on valide l'importance de chaque composant :

• CDDHFs-GLCT affiche les meilleures performances dans la plupart des cas
• CM-CC-CM-GLCT présente des avantages dans certains scénarios spécifiques
• La base laplacienne et la base de matrice d'adjacence ont chacune leurs scénarios d'application

Travaux Connexes

Fondamentaux du Traitement des Signaux Graphiques

• Transformation de Fourier Graphique (GFT) : Extension de la transformée de Fourier classique aux données de structures graphiques
• Transformation de Fourier Fractionnaire Graphique (GFRFT) : Introduction de paramètres fractionnaires, interpolation entre la transformation identité et la GFT complète

Développement de la Transformation Canonique Linéaire

• LCT Classique : Généralisation de la rotation en transformation affine, extension de FT et FRFT
• Extension au Domaine Graphique : Zhang et al. définissent GLCT basée sur CDDHFs, Li et al. proposent la réalisation CM-CC-CM

Extension du Filtrage de Wiener

Le filtrage de Wiener graphique traditionnel opère dans un domaine de transformation fixe, cet article l'étend à la sélection de domaine de transformation apprenable.

Conclusion et Discussion

Conclusions Principales

1. Le cadre d'optimisation conjointe résout efficacement les problèmes de complexité de calcul et de stabilité des méthodes traditionnelles
2. La paramétrisation riche de GLCT combinée à l'apprentissage conjoint réalise une meilleure séparation signal-bruit
3. L'optimisation de bout en bout évite les résultats sous-optimaux qui pourraient résulter de l'optimisation par étapes

Limitations

1. Optimisation non-convexe : Le problème d'optimisation conjointe est intrinsèquement non-convexe, pouvant présenter des optima locaux
2. Sensibilité à l'initialisation des paramètres : La performance de l'optimisation Adam peut dépendre de la stratégie d'initialisation
3. Garanties théoriques de convergence : Absence d'analyse théorique rigoureuse de convergence globale

Directions Futures

1. Développer des garanties théoriques de convergence plus fortes
2. Explorer des stratégies d'initialisation de paramètres adaptatives
3. Étendre au traitement des signaux graphiques temporels et dynamiques

Évaluation Approfondie

Points Forts

1. Contributions théoriques solides : Fournit des preuves complètes de différentiabilité et une analyse des propriétés
2. Innovation méthodologique forte : Première réalisation de l'optimisation conjointe des paramètres GLCT et du filtre
3. Vérification expérimentale suffisante : Validation de l'efficacité de la méthode sur plusieurs ensembles de données et configurations
4. Amélioration significative de l'efficacité de calcul : Amélioration de la complexité de $O(N^4)$ à $O(N^2)$

Insuffisances

1. Analyse de convergence insuffisante : Manque d'analyse théorique approfondie de la convergence de l'optimisation non-convexe
2. Sensibilité des paramètres : Analyse limitée de la sensibilité au taux d'apprentissage et à l'initialisation
3. Limitation des scénarios d'application : Concentration principale sur les tâches de débruitage, applicabilité à d'autres tâches de traitement de signaux graphiques à vérifier

Impact

1. Valeur académique : Fournit un nouveau paradigme d'optimisation pour le domaine du traitement des signaux graphiques
2. Valeur pratique : La complexité de calcul considérablement réduite rend les applications à grande échelle possibles
3. Reproductibilité : Fournit une description détaillée des algorithmes et des configurations expérimentales

Scénarios d'Application

• Tâches de débruitage de signaux graphiques à grande échelle
• Applications de traitement de signaux graphiques nécessitant un traitement en temps réel
• Systèmes embarqués avec exigences strictes d'efficacité de calcul
• Analyse de structures graphiques multi-échelles

Références

L'article cite 49 références connexes, couvrant les travaux importants dans les domaines fondamentaux du traitement des signaux graphiques, de la transformation canonique linéaire et du filtrage de Wiener, fournissant une base théorique solide pour la recherche.

Évaluation Globale : Cet article apporte des contributions importantes au domaine du traitement des signaux graphiques, résolvant efficacement le problème de complexité de calcul des méthodes traditionnelles par le biais d'un cadre d'optimisation conjointe, avec une analyse théorique suffisante et une vérification expérimentale complète, possédant une valeur académique et pratique considérable.