Многомерное преобразование
В математическом анализе и приложениях многомерные преобразования используются для анализа частотного содержания сигналов в области двух или более измерений.
Многомерное преобразование Фурье
[ редактировать ]Одним из наиболее популярных многомерных преобразований является преобразование Фурье , которое преобразует сигнал из представления во временной/пространственной области в представление в частотной области. [1] Многомерное преобразование Фурье (FT) в дискретной области можно вычислить следующим образом:
где F означает многомерное преобразование Фурье, m означает многомерное измерение. Задайте f как многомерный сигнал дискретной области. Обратное многомерное преобразование Фурье имеет вид
Многомерное преобразование Фурье для сигналов в непрерывной области определяется следующим образом: [1]
Свойства преобразования Фурье
[ редактировать ]Применяются аналогичные свойства одномерного преобразования Фурье, но вместо того, чтобы входной параметр был единственной записью, это многомерный (MD) массив или вектор. Следовательно, это x(n 1 ,…,n M ) вместо x(n).
Линейность
[ редактировать ]если , и затем,
Сдвиг
[ редактировать ]если , затем
если , затем
Умножение
[ редактировать ]если , и
затем,
( MD-свертка в частотной области ) |
или,
( MD-свертка в частотной области ) |
Дифференциация
[ редактировать ]Если , затем
Транспонирование
[ редактировать ]Если , затем
Отражение
[ редактировать ]Если , затем
Комплексное сопряжение
[ редактировать ]Если , затем
Теорема Парсеваля (MD)
[ редактировать ]если , и затем,
если , затем
Особым случаем теоремы Парсеваля является ситуация, когда два многомерных сигнала одинаковы. В этом случае теорема описывает сохранение энергии сигнала, а член суммирования или интеграла представляет собой плотность энергии сигнала.
Разделимость
[ редактировать ]Сигнал или система называются разделимыми, если их можно выразить как произведение одномерных функций с различными независимыми переменными. Это явление позволяет вычислять преобразование FT как произведение одномерных FT вместо многомерных FT.
если , , ... , и если , затем
, так
МД БПФ
[ редактировать ]( БПФ Быстрое преобразование Фурье ) — это алгоритм вычисления дискретного преобразования Фурье (ДПФ) и обратного ему преобразования. БПФ вычисляет ДПФ и дает точно такой же результат, как и непосредственная оценка определения ДПФ; единственная разница в том, что БПФ работает намного быстрее. (При наличии ошибки округления многие алгоритмы БПФ также намного более точны, чем непосредственная оценка определения ДПФ). Существует множество различных алгоритмов БПФ, включающих широкий диапазон математики, от простой арифметики комплексных чисел до теории групп и чисел. теория. Подробнее см. в БПФ .
МД ДПФ
[ редактировать ]Многомерное дискретное преобразование Фурье (ДПФ) представляет собой дискретную версию ПФ в дискретной области, оценивающую его на равномерно распределенных частотах дискретизации. [2] ДПФ N 1 × N 2 × ... Н м определяется выражением:
для 0 ≤ K i ≤ N i - 1 , я знак равно 1, 2, ..., м .
Обратное многомерное уравнение ДПФ имеет вид
для 0 ≤ n 1 , n 2 , ... , n m ≤ N (1, 2, ... , m ) – 1 .
Многомерное дискретное косинусное преобразование
[ редактировать ]Дискретное косинусное преобразование (DCT) используется в широком спектре приложений, таких как сжатие данных , извлечение признаков , реконструкция изображений , многокадровое обнаружение и так далее. Многомерное DCT определяется следующим образом:
для k i знак равно 0, 1, ..., N i - 1 , я знак равно 1, 2, ..., р .
Многомерное преобразование Лапласа
[ редактировать ]Многомерное преобразование Лапласа полезно для решения краевых задач. Краевые задачи с двумя или более переменными, характеризуемые уравнениями в частных производных, могут быть решены путем прямого использования преобразования Лапласа. [3] Преобразование Лапласа для M-мерного случая определено [3] как
где F обозначает представление сигнала f(t) в s-области.
Определен частный случай (в двух измерениях) многомерного преобразования Лапласа функции f(x,y). [4] как
называется изображением и известен как оригинал . [ нужна ссылка ] Этот особый случай можно использовать для решения уравнений Телеграфиста . [ нужна ссылка ] }
Многомерное Z-преобразование [5]
[ редактировать ]Многомерное Z-преобразование используется для отображения многомерного сигнала в области дискретного времени в область Z. Это можно использовать для проверки стабильности фильтров. Уравнение многомерного Z-преобразования имеет вид
где F обозначает представление сигнала f(n) в z-области.
Особым случаем многомерного Z-преобразования является 2D Z-преобразование, которое задается как
Преобразование Фурье является частным случаем преобразования Z, оцениваемого по единичному кругу (в 1D) и единичному би-кругу (в 2D). то есть в
где z и w — векторы.
Область конвергенции
[ редактировать ]Точки ( z 1 , z 2 ), для которых находятся в РПЦ.
Пример:
Если последовательность имеет опору, как показано на рисунке 1.1а, то ее ROC показан на рисунке 1.1b. Отсюда следует, что | F ( z 1 , z 2 )| < ∞ .
лежит в РПЦ, то все точки которые удовлетворяют |z1|≥|z01| и |z2|≥|z02 лежат в ROC.
Следовательно, для рисунков 1.1a и 1.1b ROC будет
где L – наклон.
Двумерное Z-преобразование , аналогичное Z-преобразованию, используется при обработке многомерных сигналов для связи двумерного сигнала дискретного времени с комплексной частотной областью, в которой известна двумерная поверхность в четырехмерном пространстве, на которой лежит преобразование Фурье. как единичная поверхность или единичный бицикл.
Приложения
[ редактировать ]DCT и DFT часто используются при обработке сигналов. [6] и обработка изображений, а также они используются для эффективного решения уравнений в частных производных спектральными методами. ДПФ также можно использовать для выполнения других операций, таких как свертка или умножение больших целых чисел. ДПФ и ДКП широко используются во многих областях, ниже мы лишь набросаем несколько примеров.
Обработка изображений
[ редактировать ]DCT используется при JPEG сжатии изображений , MJPEG , MPEG , DV , Daala и Theora сжатии видео . Там вычисляются двумерные DCT-II из блоков N x N , а результаты квантоваются и энтропийно кодируются . В этом случае N обычно равно 8, и формула DCT-II применяется к каждой строке и столбцу блока. Результатом является массив коэффициентов преобразования 8x8, в котором элемент (0,0) (вверху слева) является компонентом постоянного тока (нулевой частоты), а записи с увеличивающимися значениями вертикального и горизонтального индекса представляют более высокие вертикальные и горизонтальные пространственные частоты, как показано на рисунке справа.
При обработке изображений можно также анализировать и описывать нетрадиционные криптографические методы, основанные на 2D DCT, для вставки невидимых двоичных водяных знаков в плоскость 2D изображения. [7] и в соответствии с различными ориентациями двумерное направленное гибридное преобразование DCT-DWT может применяться для шумоподавления ультразвуковых изображений. [8] 3-D DCT также можно использовать для преобразования видеоданных или данных трехмерного изображения в схемах внедрения водяных знаков в области преобразования. [9] [10]
Спектральный анализ
[ редактировать ]Когда ДПФ используется для спектрального анализа , последовательность { x n } обычно представляет собой конечный набор равномерно расположенных временных отсчетов некоторого сигнала x ( t ), где t представляет время. Преобразование непрерывного времени в выборки (дискретное время) изменяет базовое преобразование Фурье x называемый ( t ) в преобразование Фурье с дискретным временем (DTFT), что обычно влечет за собой тип искажения, сглаживанием . Выбор подходящей частоты дискретизации (см. Частота Найквиста ) является ключом к минимизации этого искажения. Аналогичным образом, преобразование очень длинной (или бесконечной) последовательности в управляемый размер влечет за собой тип искажения, называемого утечкой , которое проявляется как потеря детализации (так называемого разрешения) в DTFT. Выбор подходящей длины подпоследовательности является основным ключом к минимизации этого эффекта. Когда доступные данные (и время на их обработку) превышают объем, необходимый для достижения желаемого разрешения по частоте, стандартным методом является выполнение нескольких ДПФ, например, для создания спектрограммы . Если желаемым результатом является спектр мощности и в данных присутствует шум или случайность, усреднение компонентов величины нескольких ДПФ является полезной процедурой для уменьшения дисперсия спектра ( также называемая периодограммой в этом контексте ); двумя примерами таких методов являются метод Уэлча и метод Бартлетта ; Общий предмет оценки спектра мощности зашумленного сигнала называется спектральной оценкой .
Последним источником искажений (или, возможно, иллюзий ) является само ДПФ, поскольку это всего лишь дискретная выборка ДПФ, которая является функцией непрерывной частотной области. Это можно смягчить, увеличив разрешение ДПФ. Эта процедура проиллюстрирована в § Выборка DTFT .
- Эту процедуру иногда называют заполнением нулями . Это конкретная реализация, используемая в сочетании с алгоритмом быстрого преобразования Фурье (БПФ). Неэффективность выполнения умножения и сложения с «выборками» с нулевым значением более чем компенсируется присущей БПФ эффективностью.
- Как уже отмечалось, утечка накладывает ограничение на собственное разрешение DTFT. Таким образом, существует практический предел выгоды, которую можно получить от мелкозернистого ДПФ.
Уравнения в частных производных
[ редактировать ]Дискретные преобразования Фурье часто используются для решения уравнений в частных производных , где снова ДПФ используется в качестве аппроксимации ряда Фурье (который восстанавливается в пределе бесконечного N ). Преимущество этого подхода в том, что он разлагает сигнал в комплексную экспоненту e inx , которые являются собственными функциями дифференцирования: d / dx e inx = в е inx . Таким образом, в представлении Фурье дифференцирование простое — мы просто умножаем на . (Однако обратите внимание, что выбор n не является уникальным из-за наложения псевдонимов; чтобы метод сходился, тригонометрической интерполяции следует использовать выбор, аналогичный выбору в разделе выше.) Линейное дифференциальное уравнение с постоянными коэффициентами преобразуется в легко решаемое алгебраическое уравнение. Затем используется обратное ДПФ для преобразования результата обратно в обычное пространственное представление. Такой подход называется спектральным методом .
ДКП также широко используются при решении уравнений в частных производных спектральными методами, где разные варианты ДКП соответствуют немного отличающимся четным/нечетным граничным условиям на двух концах массива.
Преобразования Лапласа используются для решения уравнений в частных производных. Общая теория получения решений в этом методе развита теоремами о преобразовании Лапласа в n измерениях. [3]
Многомерное Z-преобразование также можно использовать для решения уравнений в частных производных. [11]
Обработка изображений для анализа поверхности произведений искусства с помощью БПФ
[ редактировать ]Очень важным фактором является то, что мы должны применять неразрушающий метод для получения информации о редких ценностях (с точки зрения HVS, сосредоточенной на всей колориметрической и пространственной информации) о произведениях искусства и нулевом нанесении им ущерба.Мы можем понять искусство, наблюдая за изменением цвета или измеряя изменение однородности поверхности. Поскольку все изображение будет очень огромным, мы используем окно двойного приподнятого косинуса, чтобы обрезать изображение: [12]
где N — размер изображения, а x , y — координаты центра изображения от 0 до N /2.Автор хотел вычислить равное значение пространственной частоты, например: [12]
где «БПФ» обозначает быстрое преобразование Фурье, а f — диапазон пространственных частот от 0 до N /2 – 1 . Предлагаемый подход к визуализации на основе БПФ представляет собой диагностическую технологию, обеспечивающую долгую жизнь и стабильность в культурном искусстве. Это простой и дешевый вариант, который можно использовать в музеи, не влияя на их повседневное использование. Но этот метод не позволяет количественно измерить скорость коррозии.
Приложение к моделированию слабо нелинейных цепей [13]
[ редактировать ]Обратное многомерное преобразование Лапласа можно применять для моделирования нелинейных цепей. Это делается путем формулирования схемы как пространства состояний и расширения обратного преобразования Лапласа на основе расширения функции Лагерра .
Метод Лагерра можно использовать для моделирования слабонелинейной схемы, а метод Лагерра может эффективно инвертировать многомерное преобразование Лапласа с высокой точностью.
Замечено, что высокая точность и значительное ускорение могут быть достигнуты при моделировании больших нелинейных цепей с использованием многомерных преобразований Лапласа.
См. также
[ редактировать ]- Дискретное косинусное преобразование
- Список преобразований, связанных с Фурье
- Список тем анализа Фурье
- Многомерная дискретная свертка
- 2D Z-преобразование
- Многомерная декомпозиция по эмпирическим модам
- Многомерная реконструкция сигнала
Ссылки
[ редактировать ]- ^ Jump up to: а б Смит, В. Справочник по быстрым преобразованиям Фурье в реальном времени: алгоритмы тестирования продуктов, Wiley_IEEE Press, издание 1, страницы 73–80, 1995 г.
- ^ Даджен и Мерсеро, Многомерная цифровая обработка сигналов, 2-е издание, 1995 г.
- ^ Jump up to: а б с Дебнат, Джояти; Дахия, РС (1 января 1989 г.). «Теоремы о многомерном преобразовании Лапласа для решения краевых задач» . Компьютеры и математика с приложениями . 18 (12): 1033–1056. дои : 10.1016/0898-1221(89)90031-X .
- ^ Операционное исчисление в двух переменных и его применение (1-е английское издание) - перевод DMG Wishart (Calcul opérationnel) .
- ^ «Народная книга» (PDF) .
- ^ Тан Сяо, Шао-хай Ху, Ян Сяо. Приложение 2-D DFT-DWT для обработки многомерных сигналов. Материалы ICSP2006, IEEE 2006 г.
- ^ Питер КУЛЛАИ, Павол САБАКАЙ, ЮзефХУСКАИ. Простые возможности применения 2D DCT в цифровом форматеКриптография монохромного изображения. Радиоэлектроника, 17-я Международная конференция, IEEE, 2007, стр. 1–6.
- ^ Синь-лин Вэнь, Ян Сяо. Двумерное направленное гибридное преобразование DCT-DWT и его применение для шумоподавления ультразвукового изображения. Обработка сигналов. ICSP 2008. 9-я Международная конференция, Страницы: 946–949.
- ^ Цзиньвэй Ван, Шиго Лянь, Чжунсюань Лю, Чжэнь Жэнь, Юэвей Дай, Хайла Ван. Схема водяных знаков изображения на основе 3-D DCT. Промышленная электроника и приложения, 1-я конференция IEEE, 2006 г., стр. 1–6.
- ^ Джин Ли, Монсеф Габбуж, Ярмо Такала, Хексин Чен. Алгоритм прямого изменения размера 3-D DCT в DCT для кодирования видео. Обработка и анализ изображений и сигналов, 2009. ISPA 2009. Материалы 6-го международного симпозиума.стр. 105–110
- ^ Грегор, Иржи (1998). «Кибернетика» (PDF) . Кибернетика . 24 .
- ^ Jump up to: а б Анджелини Э., Грассен С.; Пиантанида, М.; Корбеллини, С.; Феррарис, Ф.; Нери, А.; Парвис, М. Обработка изображений на основе БПФ для мониторинга культурного наследия. Конференция по приборостроению и измерительным технологиям (I2MTC), 2010 г., IEEE
- ^ Ван, Тинтин (2012). «Анализ слабонелинейных цепей на основе быстрого многомерного обратного преобразования Лапласа» . 17-я Азиатско-Тихоокеанская конференция по автоматизации проектирования . стр. 547–552. дои : 10.1109/ASPDAC.2012.6165013 . ISBN 978-1-4673-0772-7 . S2CID 15427178 .