Фильтрация графических сигналов Винера в линейной канонической области: теория и проектирование методов

Основная информация

• ID статьи: 2510.10512
• Название: Graph Signal Wiener Filtering in the Linear Canonical Domain: Theory and Method Design
• Авторы: Xiaopeng Cheng, Zhichao Zhang
• Классификация: eess.SP (Обработка сигналов)
• Дата публикации: 12 октября 2025 г.
• Ссылка на статью: https://arxiv.org/abs/2510.10512

Аннотация

Методы фильтрации, основанные на графическом линейном каноническом преобразовании (GLCT), обычно оптимизируют параметры преобразования и фильтр отдельно, что приводит к высокой вычислительной стоимости и ограниченной стабильности. Для решения этой проблемы в статье предлагается обучаемая структура совместной оптимизации, которая объединяет параметры GLCT и фильтрацию Винера в процесс сквозного обучения, достигая синергетической оптимизации между построением области преобразования и операциями фильтрации. Предложенный метод не только устраняет утомительный поиск по сетке, требуемый традиционными стратегиями, но и значительно повышает гибкость и стабильность обучения системы фильтрации. Результаты экспериментов на реальных графических данных показывают, что предложенный метод превосходит существующие методы в задачах подавления шума с лучшей производительностью подавления, повышенной робастностью и более низкой вычислительной сложностью.

Исследовательский контекст и мотивация

Определение проблемы

В социальных сетях, транспортных системах, сетях биологических молекул и других нерегулярных структурах данные часто располагаются на неевклидовых сетках, что делает классические методы обработки сигналов неприменимыми. Обработка графических сигналов (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 и коэффициентов фильтра, с проверкой её эффективности и робастности в задачах подавления шума графических сигналов

Описание методов

Определение задачи

Дана модель наблюдения: $\tilde{f} = Gf + n$, где $G$ — известная матрица возмущения, $f$ — гладкий сигнал, $n$ — член аддитивного шума. Цель состоит в разработке оптимального метода фильтрации для восстановления исходного сигнала $f$ в спектральной области преобразования с минимальной среднеквадратичной ошибкой (MSE).

Основные технические компоненты

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. Целевая функция совместной оптимизации

Двухэтапный процесс оптимизации традиционных методов:

1. Фиксирование параметров преобразования $(a,b,d)$, решение оптимального фильтра $H^*$
2. Поиск по сетке $(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. while условие остановки не выполнено:
   - Вычисление градиента ∇_{H,a,b,d}MSE
   - Обновление параметров: H ← H - ε∇_H MSE
   - Обновление параметров GLCT: a,b,d ← a,b,d - ε∇_{a,b,d}MSE

Технические инновации

1. Гарантия дифференцируемости: доказана дифференцируемость функции потерь относительно параметров преобразования и коэффициентов фильтра, что делает возможной сквозную оптимизацию
2. Оптимизация вычислительной сложности:
   • Сложность поиска по сетке: $O(n_a n_b n_d N^4)$
   • Сложность совместной оптимизации Adam: $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 последовательно превосходят методы GFRFT:

Результаты на реальных данных

На наборе данных SST wAdj-CDDHFs-GLCT достигает минимального значения MSE 1.442 при установке k=2, s=0.5, что представляет улучшение примерно на 25% по сравнению с традиционным методом GFRFT.

Анализ производительности

1. Производительность подавления шума: методы GLCT демонстрируют превосходную производительность подавления шума при всех условиях тестирования
2. Робастность: при увеличении интенсивности шума производительность методов GLCT деградирует более плавно
3. Вычислительная эффективность: совместная оптимизация Adam значительно снижает вычислительную сложность

Абляционные исследования

Путём сравнения различных вариантов GLCT проверена важность каждого компонента:

• CDDHFs-GLCT показывает лучшую производительность в большинстве случаев
• CM-CC-CM-GLCT имеет преимущества в некоторых специфических сценариях
• Лапласовский базис и базис матрицы смежности имеют свои области применения

Связанные работы

Основы обработки графических сигналов

• Графическое преобразование Фурье (GFT): расширение классического преобразования Фурье на данные графических структур
• Графическое дробное преобразование Фурье (GFRFT): введение дробного параметра, интерполирующего между тождественным преобразованием и полным GFT

Развитие линейного канонического преобразования

• Классическое LCT: обобщение вращения в аффинное преобразование, расширение FT и FRFT
• Расширение на графический домен: определение GLCT Zhang и соавторами на основе CDDHFs, реализация CM-CC-CM Li и соавторами

Расширение фильтрации Винера

Традиционная фильтрация Винера для графических сигналов работает в фиксированной области преобразования; в данной статье расширено на выбираемую область преобразования.

Выводы и обсуждение

Основные выводы

1. Структура совместной оптимизации эффективно решает проблемы вычислительной сложности и стабильности традиционных методов
2. Богатая параметризация GLCT в сочетании с совместным обучением достигает лучшего разделения сигнала и шума
3. Сквозная оптимизация избегает субоптимальных результатов, которые могут возникнуть при пошаговой оптимизации

Ограничения

1. Невыпуклая оптимизация: задача совместной оптимизации по своей природе невыпукла и может иметь локальные оптимумы
2. Чувствительность инициализации параметров: производительность оптимизации Adam может зависеть от стратегии инициализации
3. Отсутствие теоретических гарантий сходимости: недостаточно строгого анализа глобальной сходимости

Будущие направления

1. Разработка более сильных теоретических гарантий сходимости
2. Исследование адаптивных стратегий инициализации параметров
3. Расширение на обработку временных графов и динамических графических сигналов

Глубокая оценка

Преимущества

1. Прочный теоретический вклад: предоставлены полные доказательства дифференцируемости и анализ свойств
2. Сильная инновационность метода: первая реализация совместной оптимизации параметров GLCT и фильтра
3. Достаточная экспериментальная проверка: метод проверен на нескольких наборах данных и установках
4. Значительное улучшение вычислительной эффективности: улучшение сложности с $O(N^4)$ до $O(N^2)$

Недостатки

1. Недостаточный анализ сходимости: отсутствует глубокий теоретический анализ сходимости невыпуклой оптимизации
2. Чувствительность параметров: ограниченный анализ чувствительности к скорости обучения и инициализации
3. Ограничение сценариев применения: основное внимание уделяется задачам подавления шума; применимость к другим задачам обработки графических сигналов требует проверки

Влияние

1. Академическая ценность: предоставляет новую парадигму оптимизации для области обработки графических сигналов
2. Практическая ценность: значительно сниженная вычислительная сложность делает возможным крупномасштабное применение
3. Воспроизводимость: предоставлены подробные описания алгоритмов и экспериментальные установки

Сценарии применения

• Крупномасштабные задачи подавления шума графических сигналов
• Приложения, требующие обработки в реальном времени
• Встроенные системы со строгими требованиями к вычислительной эффективности
• Анализ графических структур на нескольких масштабах

Библиография

В статье цитируется 49 связанных работ, охватывающих важные работы в основных областях обработки графических сигналов, линейного канонического преобразования и фильтрации Винера, обеспечивающие прочную теоретическую основу для исследования.

Общая оценка: Данная статья вносит важный вклад в область обработки графических сигналов, эффективно решая проблему вычислительной сложности традиционных методов посредством структуры совместной оптимизации. Теоретический анализ полный, экспериментальная проверка всеобъемлющая, работа имеет высокую академическую и практическую ценность.