Graphensignal-Wiener-Filterung im linearen kanonischen Bereich: Theorie und Methodendesign

Grundlegende Informationen

• Papier-ID: 2510.10512
• Titel: Graph Signal Wiener Filtering in the Linear Canonical Domain: Theory and Method Design
• Autoren: Xiaopeng Cheng, Zhichao Zhang
• Klassifizierung: eess.SP (Signalverarbeitung)
• Veröffentlichungsdatum: 12. Oktober 2025
• Papierlink: https://arxiv.org/abs/2510.10512

Zusammenfassung

Filtermethoden, die auf der graphischen linearen kanonischen Transformation (GLCT) basieren, optimieren typischerweise Transformationsparameter und Filter separat, was zu hohen Rechenkosten und begrenzter Stabilität führt. Um dieses Problem zu lösen, wird in diesem Papier ein trainierbares gemeinsames Optimierungsframework vorgeschlagen, das GLCT-Parameter und Wiener-Filterung in einen End-to-End-Lernprozess integriert und eine synergistische Optimierung zwischen Transformationsbereichskonstruktion und Filteroperationen erreicht. Die Methode eliminiert nicht nur die aufwändige Gittersuche, die von traditionellen Strategien erforderlich ist, sondern verbessert auch erheblich die Flexibilität und Trainingsstabilität des Filtersystems. Experimentelle Ergebnisse auf echten Graphdaten zeigen, dass die vorgeschlagene Methode bestehende Methoden bei Rauschunterdrückungsaufgaben übertrifft und bessere Rauschunterdrückungsleistung, höhere Robustheit und niedrigere Rechenkomplexität bietet.

Forschungshintergrund und Motivation

Problemdefinition

In unregelmäßigen Strukturen wie sozialen Netzwerken, Verkehrssystemen und biologischen Molekülnetzwerken liegen Daten häufig auf nicht-euklidischen Gittern vor, wodurch klassische Signalverarbeitungsmethoden nicht mehr anwendbar sind. Die Graphensignalverarbeitung (GSP) entstand, um unregelmäßige Strukturdaten als Graphen zu modellieren, wobei Knoten Datenentitäten darstellen, Kanten ihre Beziehungen kodieren und Signalwerte an Knoten angebracht sind.

Kernherausforderungen

1. Rauschstörung: Graphensignale werden während der Erfassung, Übertragung und Speicherung unvermeidlich durch Rauschen gestört
2. Adaptivität der Filtertheorie: Klassische lineare Filterung basiert auf euklidischen Raumeigenschaften und lässt sich nicht direkt auf Graphstrukturen übertragen, die nicht-euklidische Räume darstellen
3. Komplexität der Parameteroptimierung: Bestehende GLCT-Filtermethoden optimieren typischerweise Transformationsparameter und Filter separat, was zu hohen Rechenkosten und begrenzter Stabilität führt

Forschungsmotivation

Die Einschränkungen bestehender Methoden manifestieren sich hauptsächlich in:

• Niedrige Recheneffizienz: Gittersuche-Strategien erfordern die Aufzählung zahlreicher Kandidatenparameterkombinationen
• Unkoordinierte Optimierung: Die Trennung von Transformationsbereichswahl und Filterdesign führt zu suboptimalen Ergebnissen
• Schlechte Skalierbarkeit: Gittersuche in hochdimensionalen Parameterräumen ist mit kombinatorischer Explosion konfrontiert

Kernbeiträge

1. Neue GLCT-Definition: Schlägt die auf Laplace-Eigenbasis basierende CM-CC-CM-GLCT vor, füllt Lücken in der bestehenden CM-CC-CM-GLCT und organisiert das CDDHFs-GLCT- und CM-CC-CM-GLCT-Framework
2. Differenzierbarkeittheorie: Beweist die Differenzierbarkeit der GLCT-Kernmodule unter gewichteten Adjazenzmatrizen und Laplace-Matrizen und bietet theoretische Unterstützung für die End-to-End-Optimierung von Transformationsparametern und Filterkoeffizienten
3. Gemeinsames Optimierungsframework: Konstruiert das GLCT-GWF-Framework, das die End-to-End-Optimierung von GLCT-Parametern und Filterkoeffizienten realisiert und seine Effektivität und Robustheit bei echten Graphensignal-Rauschunterdrückungsaufgaben validiert

Methodendetails

Aufgabendefinition

Gegeben das Beobachtungsmodell: $\tilde{f} = Gf + n$, wobei $G$ eine bekannte Störungsmatrix ist, $f$ ein glattes Signal und $n$ ein additiver Rauschterm. Das Ziel ist es, eine optimale Filtermethode im transformierten Spektralbereich zu entwerfen, um das ursprüngliche Signal $f$ mit minimalem mittleren quadratischen Fehler (MSE) wiederherzustellen.

Kernkomponenten der Technik

1. Graphische lineare kanonische Transformation (GLCT)

Die GLCT wird durch eine 2×2-Matrix $M = (a,b;c,d)$ bestimmt, wobei $ad-bc=1$. Diese Matrix kann zerlegt werden als: $\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}$

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

2. Lap-CM-CC-CM-GLCT-Definition

Die in diesem Papier vorgeschlagene auf der Laplace-Matrix basierende CM-CC-CM-GLCT-Definition ist: $\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. Gemeinsames Optimierungsziel

Der zweistufige Optimierungsprozess traditioneller Methoden:

1. Fixieren Sie Transformationsparameter $(a,b,d)$ und lösen Sie den optimalen Filter $H^*$
2. Gittersuche $(a,b,d)$ zur Minimierung von MSE

Die in diesem Papier vorgeschlagene gemeinsame Optimierung: $\min_{a,b,d,H} E\{\|F^{M^{-1}}_L HF^M_L \tilde{f} - f\|^2_2\}$

Algorithmus-Framework

Algorithmus 1: Traditionelle Gittersuchmethode

Eingabe: Graphensignal f, Zielsignal f̃, Parametergitter A,B,D
Ausgabe: Optimale Parameter (a*,b*,d*), optimaler Filter H*
1. Spektrale Basis vorberechnen (Eigendekomposition)
2. for a ∈ A, b ∈ B, d ∈ D:
   - GLCT-Operatoren F^M und F^{M^{-1}} konstruieren
   - Wiener-Hopf-Gleichung lösen: h = T^{-1}q
   - Verlust MSE(H,a,b,d) evaluieren
   - Optimale Lösung aktualisieren

Algorithmus 2: Adam-Methode zur gemeinsamen Optimierung

Eingabe: Graphensignal f, Zielsignal f̃, Lernrate ε
Ausgabe: Gelernter Filter und GLCT-Parameter
1. Parameter (a₀,b₀,d₀) und Filter H₀ initialisieren
2. while Stoppbedingung nicht erfüllt:
   - Gradienten ∇_{H,a,b,d}MSE berechnen
   - Parameter aktualisieren: H ← H - ε∇_H MSE
   - GLCT-Parameter aktualisieren: a,b,d ← a,b,d - ε∇_{a,b,d}MSE

Technische Innovationspunkte

1. Differenzierbarkeitsgarantie: Beweist die Differenzierbarkeit der Verlustfunktion bezüglich Transformationsparametern und Filterkoeffizienten, wodurch End-to-End-Optimierung möglich wird
2. Optimierung der Rechenkomplexität:
   • Gittersuche-Komplexität: $O(n_a n_b n_d N^4)$
   • Adam-Komplexität der gemeinsamen Optimierung: $O(KN^2)$
3. Theoretische Eigenschaften: Die neu vorgeschlagene Lap-CM-CC-CM-GLCT erfüllt wichtige Eigenschaften wie Linearität, Nullrotation, Additivität, Invertierbarkeit und Unitarität

Experimentelle Einrichtung

Datensätze

Synthetische Daten

1. 5-nn-Graph: 15 Knoten, jeder Knoten verbunden mit 5 nächsten Nachbarn
2. Schweizer Rolle-Graph: 30 Knoten Mannigfaltigkeitsstruktur
3. Sensorgraph: 20-Knoten-Sensornetzwerk

Echte Daten

1. Meeresoberflächentemperatur (SST): 50 Knoten, k∈{2,6,10}
2. PM2.5: Luftqualitätsdaten, 50 Knoten
3. COVID-19: Epidemiologische Ausbreitungsdaten, 50 Knoten

Bewertungsmetriken

Der mittlere quadratische Fehler (MSE) wird als Hauptbewertungsmetrik zur Bewertung der Rauschunterdrückungsleistung verwendet.

Vergleichsmethoden

• GFRFT_W: Graphische Bruchfourier-Transformation basierend auf gewichteter Adjazenzmatrix
• GFRFT_L: Graphische Bruchfourier-Transformation basierend auf Laplace-Matrix
• Verschiedene GLCT-Varianten: wAdj-CDDHFs-GLCT, Lap-CDDHFs-GLCT usw.

Implementierungsdetails

• Gittersuche: Parameterbereich 0,2, Schrittweite 0,1
• Adam-Optimierung: Lernrate 0,005, 5000 Iterationen
• Rauscheinstellung: Gaußsches Rauschen, Standardabweichung s∈{0,5,1,0,1,5} (synthetische Daten), s∈{0,5,0,6,0,7} (echte Daten)

Experimentelle Ergebnisse

Hauptergebnisse

Ergebnisse synthetischer Daten

Experimente auf drei synthetischen Graphen zeigen, dass GLCT-Methoden GFRFT-Methoden konsistent übertreffen:

Ergebnisse echter Daten

Auf dem SST-Datensatz erreicht wAdj-CDDHFs-GLCT unter der Einstellung k=2, s=0,5 den niedrigsten MSE-Wert von 1,442, was einer Verbesserung von etwa 25% gegenüber traditionellen GFRFT-Methoden entspricht.

Leistungsanalyse

1. Rauschunterdrückungsleistung: GLCT-Methoden zeigen unter allen Testbedingungen überlegene Rauschunterdrückungsleistung
2. Robustheit: Mit zunehmender Rauschintensität verschlechtert sich die Leistung der GLCT-Methode eleganter
3. Recheneffizienz: Adam-Optimierung zur gemeinsamen Optimierung reduziert die Rechenkomplexität erheblich

Ablationsstudien

Durch den Vergleich verschiedener GLCT-Varianten wird die Bedeutung jeder Komponente validiert:

• CDDHFs-GLCT zeigt in den meisten Fällen die beste Leistung
• CM-CC-CM-GLCT hat Vorteile in bestimmten spezifischen Szenarien
• Laplace-Basis und Adjazenzmatrix-Basis haben jeweils anwendbare Szenarien

Verwandte Arbeiten

Grundlagen der Graphensignalverarbeitung

• Graphische Fourier-Transformation (GFT): Erweitert die klassische Fourier-Transformation auf Graphstrukturdaten
• Graphische Bruchfourier-Transformation (GFRFT): Führt Bruchordnungsparameter ein und interpoliert zwischen Identitätstransformation und vollständiger GFT

Entwicklung der linearen kanonischen Transformation

• Klassische LCT: Verallgemeinert Rotation zu affiner Transformation und erweitert FT und FRFT
• Graphische Domänenerweiterung: Zhang et al. definieren GLCT basierend auf CDDHFs, Li et al. schlagen CM-CC-CM-Implementierung vor

Wiener-Filterung-Erweiterung

Traditionelle Graphen-Wiener-Filterung arbeitet in fester Transformationsdomäne; dieses Papier erweitert sie auf lernbare Transformationsbereichswahl.

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

1. Das gemeinsame Optimierungsframework löst effektiv die Probleme der Rechenkomplexität und Stabilität traditioneller Methoden
2. Die reichhaltige Parametrisierung von GLCT kombiniert mit gemeinsamen Lernverfahren erreicht bessere Signal-Rausch-Trennung
3. End-to-End-Optimierung vermeidet suboptimale Ergebnisse, die durch mehrstufige Optimierung entstehen können

Einschränkungen

1. Nicht-konvexe Optimierung: Das gemeinsame Optimierungsproblem ist von Natur aus nicht-konvex und kann lokale Optima aufweisen
2. Empfindlichkeit der Parameterinitialisierung: Die Leistung der Adam-Optimierung kann von der Initialisierungsstrategie abhängen
3. Theoretische Konvergenzgarantie: Mangel an strenger globaler Konvergenztheorieanalyse

Zukünftige Richtungen

1. Entwicklung stärkerer theoretischer Konvergenzgarantien
2. Erkundung adaptiver Parameterinitialisierungsstrategien
3. Erweiterung auf zeitvariable Graphen und dynamische Graphensignalverarbeitung

Tiefgreifende Bewertung

Stärken

1. Solide theoretische Beiträge: Bietet vollständige Differenzierbarkeitsnachweise und Eigenschaftsanalyse
2. Starke Methodennovation: Erste Realisierung der gemeinsamen Optimierung von GLCT-Parametern und Filtern
3. Umfassende experimentelle Validierung: Validiert die Methodeneffektivität unter mehreren Datensätzen und Einstellungen
4. Signifikante Verbesserung der Recheneffizienz: Komplexitätsverbesserung von $O(N^4)$ auf $O(N^2)$

Mängel

1. Unzureichende Konvergenzanalyse: Mangel an tiefgreifender theoretischer Analyse der Konvergenz nicht-konvexer Optimierung
2. Parameterempfindlichkeit: Begrenzte Analyse der Empfindlichkeit gegenüber Lernrate und Initialisierung
3. Eingeschränkte Anwendungsszenarien: Konzentriert sich hauptsächlich auf Rauschunterdrückungsaufgaben; Anwendbarkeit auf andere Graphensignalverarbeitungsaufgaben bleibt zu überprüfen

Einfluss

1. Akademischer Wert: Bietet ein neues Optimierungsparadigma für das Graphensignalverarbeitungsfeld
2. Praktischer Wert: Signifikant reduzierte Rechenkomplexität ermöglicht großflächige Anwendungen
3. Reproduzierbarkeit: Bietet detaillierte Algorithmusbeschreibungen und experimentelle Einrichtungen

Anwendungsszenarien

• Großflächige Graphensignal-Rauschunterdrückungsaufgaben
• Graphensignalanwendungen, die Echtzeitverarbeitung erfordern
• Eingebettete Systeme mit strengeren Anforderungen an Recheneffizienz
• Mehrskalige Graphenstrukturanalyse

Literaturverzeichnis

Das Papier zitiert 49 verwandte Referenzen, die wichtige Arbeiten in den Kernbereichen Graphensignalverarbeitung, lineare kanonische Transformation und Wiener-Filterung abdecken und eine solide theoretische Grundlage für die Forschung bieten.

Gesamtbewertung: Dieses Papier leistet wichtige Beiträge zum Graphensignalverarbeitungsfeld, löst effektiv das Rechenkomplexitätsproblem traditioneller Methoden durch ein gemeinsames Optimierungsframework, bietet umfassende theoretische Analysen und umfassende experimentelle Validierung und hat hohen akademischen und praktischen Wert.