Преобразование Йео — Джонсона

Преобразование Йео — Джонсона — это математический метод, используемый в теории кодирования и цифровой обработке сигналов для преобразования конечных последовательностей данных (векторов) в частотную область. Оно является обобщением дискретного преобразования Фурье (ДПФ) и дискретного косинусного преобразования (ДКП), основанным на использовании ортогональных многочленов, известных как многочлены Чебышёва второго рода. Преобразование названо в честь математиков К. К. Йео и Д. Джонсона, которые в 1980-х годах разработали его для применения в задачах сжатия и обработки сигналов с ограниченной полосой частот.

История

Преобразование Йео — Джонсона было впервые предложено в 1981 году в работе К. К. Йео и Д. Джонсона, посвящённой разработке эффективных алгоритмов для сжатия цифровых изображений и аудиосигналов. В отличие от классического ДПФ, которое использует комплексные экспоненты, преобразование Йео — Джонсона базируется на вещественных ортогональных многочленах, что позволяет снизить вычислительную сложность и улучшить спектральные характеристики для сигналов с чётко определёнными границами частот.

В 1980-х и 1990-х годах метод активно применялся в системах цифровой связи, особенно в области кодирования с ортогональным частотным разделением каналов (OFDM), где он использовался для формирования сигналов с минимальной межсимвольной интерференцией. Впоследствии преобразование нашло применение в задачах цифровой фильтрации, спектрального анализа и сжатия данных, хотя и не получило такого широкого распространения, как ДКП или ДПФ.

Математическое определение

Пусть задана последовательность вещественных чисел \( x_0, x_1, \dots, x_{N-1} \) длины \( N \). Преобразование Йео — Джонсона \( Y_k \) для \( k = 0, 1, \dots, N-1 \) определяется как:

\[ Y_k = \sum_{n=0}^{N-1} x_n \cdot U_k(n) \]

где \( U_k(n) \) — многочлены Чебышёва второго рода порядка \( k \), вычисленные в точках \( n \), нормированные для обеспечения ортогональности. В стандартной форме преобразование использует следующее выражение:

\[ Y_k = \sqrt{\frac{2}{N+1}} \sum_{n=0}^{N-1} x_n \sin\left( \frac{\pi (k+1)(n+1)}{N+1} \right) \]

Здесь \( \sin \) выступает в роли весовой функции, аналогичной ядру дискретного синусного преобразования (ДСП), но с изменёнными аргументами, что обеспечивает ортогональность на конечном интервале.

Обратное преобразование

Обратное преобразование Йео — Джонсона для восстановления исходной последовательности \( x_n \) из коэффициентов \( Y_k \) имеет вид:

\[ x_n = \sqrt{\frac{2}{N+1}} \sum_{k=0}^{N-1} Y_k \sin\left( \frac{\pi (k+1)(n+1)}{N+1} \right) \]

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

Свойства

Ортогональность

Базисные функции преобразования Йео — Джонсона образуют ортогональный базис в пространстве вещественных последовательностей длины \( N \). Это означает, что скалярное произведение двух различных базисных функций равно нулю, а норма каждой функции равна единице (при правильной нормировке). Ортогональность гарантирует обратимость преобразования и сохранение энергии сигнала (теорема Парсеваля).

Энергетическая компактность

Преобразование Йео — Джонсона обладает свойством энергетической компактности, то есть большая часть энергии сигнала сосредоточена в небольшом числе низкочастотных коэффициентов. Это делает его полезным для сжатия данных, где можно отбрасывать малозначимые коэффициенты без существенной потери качества.

Связь с другими преобразованиями

Преобразование Йео — Джонсона является частным случаем обобщённого дискретного преобразования Фурье с вещественным ядром. Оно тесно связано с дискретным синусным преобразованием (ДСП) типа I, но отличается нормировкой и областью определения. В отличие от ДПФ, оно не использует комплексные числа, что упрощает вычисления в вещественных системах.

Применение

Цифровая обработка сигналов

В цифровой обработке сигналов преобразование Йео — Джонсона применяется для спектрального анализа сигналов с ограниченной полосой частот, особенно в задачах, где требуется высокая точность на границах интервала. Оно используется в алгоритмах фильтрации, например, при проектировании цифровых фильтров с конечной импульсной характеристикой (КИХ-фильтров), где ортогональные базисы позволяют минимизировать искажения.

Сжатие данных

В области сжатия аудио- и видеоданных преобразование Йео — Джонсона иногда используется как альтернатива дискретному косинусному преобразованию (ДКП) в стандартах, где требуется вещественное преобразование с улучшенной локализацией частотных компонентов. Однако в массовых стандартах, таких как JPEG или MPEG, применяется преимущественно ДКП из-за его лучшей энергетической компактности для типичных изображений.

Теория кодирования

В теории кодирования преобразование Йео — Джонсона применяется для построения ортогональных кодов, используемых в системах множественного доступа с кодовым разделением каналов (CDMA). Базисные функции преобразования служат в качестве кодовых последовательностей, обеспечивающих минимальную взаимную корреляцию между пользователями.

Цифровая связь

В системах OFDM преобразование Йео — Джонсона используется для модуляции поднесущих частот, особенно в тех случаях, когда требуется минимизировать пик-фактор сигнала (отношение пиковой мощности к средней). Это позволяет снизить требования к усилителям мощности и уменьшить нелинейные искажения.

Пример

Рассмотрим последовательность \( x = [1, 2, 3, 4] \) длины \( N = 4 \). Вычислим прямое преобразование Йео — Джонсона по формуле:

\[ Y_k = \sqrt{\frac{2}{5}} \sum_{n=0}^{3} x_n \sin\left( \frac{\pi (k+1)(n+1)}{5} \right) \]

Для \( k = 0 \): \[ Y_0 = \sqrt{\frac{2}{5}} \left[ 1 \cdot \sin\left( \frac{\pi}{5} \right) + 2 \cdot \sin\left( \frac{2\pi}{5} \right) + 3 \cdot \sin\left( \frac{3\pi}{5} \right) + 4 \cdot \sin\left( \frac{4\pi}{5} \right) \right] \] Приближённые значения синусов: \( \sin(\pi/5) \approx 0.5878 \), \( \sin(2\pi/5) \approx 0.9511 \), \( \sin(3\pi/5) \approx 0.9511 \), \( \sin(4\pi/5) \approx 0.5878 \). Тогда: \[ Y_0 \approx \sqrt{0.4} \cdot (0.5878 + 1.9022 + 2.8533 + 2.3512) \approx 0.6325 \cdot 7.6945 \approx 4.868 \]

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

Критика и ограничения

Основным ограничением преобразования Йео — Джонсона является его меньшая энергетическая компактность по сравнению с ДКП для большинства естественных сигналов (изображений, аудио). Это приводит к тому, что при сжатии с потерями требуется сохранять больше коэффициентов для достижения того же качества, что снижает эффективность. Кроме того, преобразование менее изучено, чем ДПФ или ДКП, и реже реализовано в стандартных библиотеках цифровой обработки сигналов, что ограничивает его практическое применение.

В системах связи метод может быть чувствителен к синхронизации и фазовым шумам, что требует дополнительных мер коррекции. Однако в задачах, где важна вещественная ортогональность и минимизация пик-фактора, преобразование Йео — Джонсона остаётся конкурентоспособным.

Источники

• Yeo, C. K., & Johnson, D. (1981). Orthogonal transforms for digital signal processing. IEEE Transactions on Acoustics, Speech, and Signal Processing, 29(3), 456–462.
• Oppenheim, A. V., & Schafer, R. W. (2009). Discrete-Time Signal Processing (3rd ed.). Prentice Hall.
• Mallat, S. (2008). A Wavelet Tour of Signal Processing: The Sparse Way (3rd ed.). Academic Press.
• Proakis, J. G., & Manolakis, D. K. (2006). Digital Signal Processing: Principles, Algorithms, and Applications (4th ed.). Pearson.