Преобразование Фурье дискретного времени
Преобразования Фурье |
---|
В математике преобразование Фурье с дискретным временем ( DTFT ) — это форма анализа Фурье , применимая к последовательности дискретных значений.
DTFT часто используется для анализа выборок непрерывной функции. Термин дискретное время относится к тому факту, что преобразование работает с дискретными данными, часто с выборками, интервал которых имеет единицы времени. Из равномерно расположенных выборок он создает функцию частоты, которая представляет собой периодическое суммирование непрерывного преобразования Фурье исходной непрерывной функции. При определенных теоретических условиях, описываемых теоремой выборки , исходная непрерывная функция может быть полностью восстановлена из DTFT и, следовательно, из исходных дискретных выборок. DTFT само по себе является непрерывной функцией частоты, но его дискретные выборки можно легко вычислить с помощью дискретного преобразования Фурье (DFT) (см. § Выборка DTFT ), которое на сегодняшний день является наиболее распространенным методом современного анализа Фурье.
Оба преобразования обратимы. Обратное DTFT — это исходная выборочная последовательность данных. Обратное ДПФ представляет собой периодическое суммирование исходной последовательности. ( БПФ Быстрое преобразование Фурье ) — это алгоритм вычисления одного цикла ДПФ, а его обратное преобразование дает один цикл обратного ДПФ.
Введение [ править ]
с Фурье преобразованием Связь
Начнем с общего определения преобразования Фурье интеграла :
Это сводится к суммированию (см. Преобразование Фурье § Численное интегрирование ряда упорядоченных пар ), когда заменяется дискретной последовательностью его выборок, для целых значений также заменяется на уход :
который представляет собой ряд Фурье по частоте с периодичностью Нижний индекс отличает его от и от формы угловой частоты DTFT. Т.е., когда частотная переменная, имеет нормированные единицы радианы /выборка , периодичность равна и ряд Фурье : [1] : стр.147
( Уравнение 1 ) |

Полезность DTFT коренится в формуле суммирования Пуассона , которая говорит нам, что периодическая функция, представленная рядом Фурье, представляет собой периодическое суммирование преобразования Фурье : [а] [А]
( Уравнение 2 ) |
Целое число имеет единицы циклов/выборки и частота дискретизации, ( выборок/сек ). Так содержит точные копии которые сдвинуты на кратное число герц и объединяются сложением. Для достаточно большого тот срок можно наблюдать в регионе с небольшими искажениями ( алиасингами ) от других терминов или без них. На рис.1 изображен пример, где недостаточно велик, чтобы предотвратить наложение псевдонимов.
Мы также отмечаем, что представляет собой преобразование Фурье Таким образом, альтернативное определение DTFT : [Б]
( Уравнение 3 ) |
Модулированная гребенчатая функция Дирака — это математическая абстракция, которую иногда называют импульсной выборкой . [3]
Обратное преобразование [ править ]
Операция, которая восстанавливает последовательность дискретных данных из функции DTFT, называется обратным DTFT . Например, обратное непрерывное преобразование Фурье обеих частей уравнения 3 дает последовательность в форме модулированной гребенчатой функции Дирака :
Однако, отметив, что является периодическим, вся необходимая информация содержится в любом интервале длины И в уравнении 1 , и в уравнении 2 суммирование по представляют собой ряд Фурье с коэффициентами Стандартные формулы для коэффициентов Фурье также являются обратными преобразованиями :
( Уравнение 4 ) |
Периодические данные [ править ]
Когда последовательность входных данных является -периодическое, уравнение 2 можно вычислительно свести к дискретному преобразованию Фурье (ДПФ), потому что :
- Вся доступная информация содержится внутри образцы.
- сходится к нулю везде, кроме целых чисел, кратных известные как гармонические частоты. На этих частотах DTFT расходится с разной частотно-зависимой скоростью. И эти скорости определяются ДПФ одного цикла последовательность.
- DTFT является периодическим, поэтому максимальное количество уникальных амплитуд гармоник равно
ДПФ одного цикла последовательность такая :
И можно выразить через обратное преобразование :
Обратное ДПФ иногда называют дискретным рядом Фурье (ДФС). [1] : стр. 542
Из-за -периодичность обеих функций это можно упростить до :
который удовлетворяет требованию обратного преобразования :
Выборка DTFT [ править ]
Когда DTFT является непрерывным, обычной практикой является вычисление произвольного количества выборок. одного цикла периодической функции : [1] : стр. 557–559 и 703.
где представляет собой периодическое суммирование :
- (см. Дискретный ряд Фурье )
The последовательность является обратным ДПФ. Таким образом, наша выборка DTFT приводит к тому, что обратное преобразование становится периодическим. Массив значения известны как периодограмма , а параметр называется NFFT в одноименной функции Matlab. [4]
Чтобы оценить один цикл численно нам требуется конечная длина последовательность. Например, длинная последовательность может быть усечена оконной функцией длины в результате были выявлены три случая, заслуживающие особого упоминания. Для простоты обозначений рассмотрим значения ниже, чтобы представить значения, измененные оконной функцией.
Случай: прореживание частоты. для некоторого целого числа (обычно 6 или 8)
Цикл сводится к суммированию сегменты длины Затем ДПФ носит разные названия, например :
- многофазное ДПФ [9] [10]
- блок многофазных фильтров [12]
- Многоблочное управление окнами и временные псевдонимы . [13]
Напомним, что прореживание выборочных данных в одном домене (по времени или частоте) приводит к перекрытию (иногда называемому псевдонимом ) в другом, и наоборот. По сравнению с -длина ДПФ, суммирование/перекрытие приводит к уменьшению частоты, [1] : стр.558 оставляя только образцы DTFT, наименее подверженные утечке спектра . Обычно это является приоритетом при реализации банка фильтров БПФ (канализатора). С обычной оконной функцией длины гребешковые потери были бы неприемлемы. Таким образом, многоблочные окна создаются с использованием инструментов проектирования КИХ-фильтров . [14] [15] Их частотный профиль плоский в самой высокой точке и быстро падает в средней точке между остальными выборками DTFT. Чем больше значение параметра тем лучше потенциальная производительность.
Случай:
Когда симметричный, -длина оконной функции ( ) усекается на 1 коэффициент, его называют периодическим или ДПФ-четным . Это обычная практика, но усечение незначительно влияет на DTFT (утечку спектра). Охарактеризовать этот эффект представляет по крайней мере академический интерес. Ан ДПФ усеченного окна по длине создает выборки частоты с интервалом вместо Образцы имеют реальную стоимость, [16] : стр.52 но их значения не совсем соответствуют DTFT симметричного окна. Периодическое суммирование, вместе с -длинное ДПФ, также может использоваться для выборки ДВПФ с интервалами Эти выборки также имеют действительные значения и точно соответствуют DTFT (пример: File:Sampling the Discrete-time Fourier Transform.svg ). Чтобы использовать полное симметричное окно для спектрального анализа на интервал, можно было бы объединить и выборки данных (дополнительно, поскольку симметричное окно взвешивает их одинаково), а затем применить усеченное симметричное окно и -длина ДПФ.


Случай: Частотная интерполяция.
В этом случае ДПФ упрощается до более привычной формы :
Чтобы воспользоваться преимуществами алгоритма быстрого преобразования Фурье для вычисления ДПФ, суммирование обычно выполняется по всем термины, хотя из них нули. Следовательно, случай часто называют заполнением нулями .
Спектральная утечка, которая увеличивается по мере уменьшается, вредно для некоторых важных показателей производительности, таких как разрешение нескольких частотных компонентов и количество шума, измеряемого каждой выборкой DTFT. Но эти вещи не всегда имеют значение, например, когда последовательность представляет собой бесшумную синусоиду (или константу), сформированную оконной функцией. Тогда обычной практикой является использование заполнения нулями для графического отображения и сравнения подробных закономерностей утечки оконных функций. Чтобы проиллюстрировать это для прямоугольного окна, рассмотрим последовательность:
- и
На рисунках 2 и 3 представлены графики величин двух ДПФ разного размера, как указано на их этикетках. В обоих случаях доминирующая составляющая находится на частоте сигнала: . На рис. 2 также видна картина спектральной утечки прямоугольное окно. Иллюзия на рис. 3 является результатом выборки DTFT только при пересечении нуля. Вместо DTFT последовательности конечной длины он создает впечатление бесконечно длинной синусоидальной последовательности. Факторами, способствующими возникновению иллюзии, являются использование прямоугольного окна и выбор частоты (1/8 = 8/64) ровно с 8 (целыми числами) циклами на 64 выборки. Окно Ханна даст аналогичный результат, за исключением того, что пик будет расширен до 3 выборок (см. Окно Ханна с четным ДПФ ).
Свертка [ править ]
Теорема свертки для последовательностей :
Важным частным случаем является круговая свертка последовательностей x и y, определяемая формулой где представляет собой периодическое суммирование. Дискретно-частотная природа означает, что произведение с непрерывной функцией также дискретно, что приводит к значительному упрощению обратного преобразования :
Для последовательностей x и y , ненулевая длительность которых меньше или равна N , окончательное упрощение :
Значение этого результата объясняется в разделах « Алгоритмы круговой свертки» и «Быстрая свертка» .
Свойства симметрии [ править ]
Когда действительная и мнимая части сложной функции разлагаются на четные и нечетные части , получается четыре компонента, обозначенные ниже индексами RE, RO, IE и IO. Между четырьмя компонентами комплексной функции времени и четырьмя компонентами ее комплексного частотного преобразования существует взаимно однозначное соответствие : [17] : стр.291
Отсюда вытекают различные зависимости, например :
- Преобразование действительной функции это четная симметричная функция И наоборот, четно-симметричное преобразование подразумевает вещественную временную область.
- Преобразование мнимой функции нечетная симметричная функция и обратное верно.
- Преобразование четно-симметричной функции это действительная функция и обратное верно.
- Преобразование нечетно-симметричной функции — мнимая функция и обратное верно.
Связь с Z-преобразованием [ править ]
представляет собой ряд Фурье , который также можно выразить через двустороннее Z-преобразование . Т.е .:
где Обозначение отличает Z-преобразование от преобразования Фурье. Следовательно, мы также можем выразить часть Z-преобразования через преобразование Фурье :
Обратите внимание, что при изменении параметра T члены оставаться постоянной разлукой друг от друга, а их ширина увеличивается или уменьшается. Члены X 1/ T ( f ) остаются постоянной ширины, а их расстояние 1/ T масштабируется вверх или вниз.
дискретных преобразований Таблица Фурье
Некоторые распространенные пары преобразований показаны в таблице ниже. Применяются следующие обозначения :
- — действительное число, представляющее непрерывную угловую частоту (в радианах на выборку). ( измеряется в циклах/сек, и находится в секундах/выборку.) Во всех случаях в таблице DTFT является 2π-периодическим (в ).
- обозначает функцию, определенную на .
- обозначает функцию, определенную на , и ноль в другом месте. Затем:
- это дельта-функция Дирака
- нормализованная функция sinc
- это функция треугольника
- n — целое число, представляющее область дискретного времени (в выборках).
- - дискретная ступенчатая функция
- это дельта Кронекера
Временной интервал х [ н ] | Частотная область Икс 2 п ( ω ) | Примечания | Ссылка |
---|---|---|---|
[17] : стр.305 | |||
целое число | |||
нечетное М | целое число | ||
The Этот термин следует интерпретировать как распределение в смысле главного значения Коши вокруг его полюсов в точке . | |||
[17] : стр.305 | |||
-π < а < π | действительное число | ||
действительное число с | |||
действительное число с | |||
целое число и нечетное целое число | |||
действительные числа с | |||
действительное число , | |||
он работает как дифференцирующий фильтр | |||
действительные числа с | |||
Преобразование Гильберта | |||
![]() | действительные числа сложный |
Свойства [ править ]
В этой таблице показаны некоторые математические операции во временной области и соответствующие эффекты в частотной области.
- это дискретная свертка двух последовательностей
- является комплексно-сопряженным числом x [ n ] .
Свойство | Временной интервал х [ н ] | Частотная область | Примечания | Ссылка |
---|---|---|---|---|
Линейность | комплексные числа | [17] : стр.294 | ||
Обращение времени/Обращение частоты | [17] : стр.297 | |||
Сопряжение времени | [17] : стр.291 | |||
Обращение времени и сопряжение | [17] : стр.291 | |||
Реальная часть времени | [17] : стр.291 | |||
Мнимая часть времени | [17] : стр.291 | |||
Действительная часть по частоте | [17] : стр.291 | |||
Мнимая часть частоты | [17] : стр.291 | |||
Сдвиг во времени/Модуляция по частоте | целое число к | [17] : стр.296 | ||
Сдвиг частоты/модуляция во времени | действительное число | [17] : стр.300 | ||
Децимация | [И] | целое число | ||
Расширение времени | целое число | [1] : стр.172 | ||
Производная по частоте | [17] : стр.303 | |||
Интегрирование по частоте | ||||
Разница во времени | ||||
Суммирование по времени | ||||
Свертка по времени/Умножение по частоте | [17] : стр.297 | |||
Умножение по времени/Свертка по частоте | Периодическая свертка | [17] : стр.302 | ||
Взаимная корреляция | ||||
Теорема Парсеваля | [17] : стр.302 |
См. также [ править ]
Примечания [ править ]
- ^ Когда зависимость от T не важна, обычно ее заменяют на Тогда f имеет единицы измерения ( циклы/выборка ), называемые нормализованной частотой .
- ^ Фактически уравнение 2 часто оправдывается следующим образом : [1] : стр.143
- ^ WOLA не следует путать с методом кусочной свертки с перекрытием .
- ^ Пример WOLA: Файл: пример канализатора WOLA.png
- ^ Это выражение получается следующим образом: [1] : стр. 168
Цитаты страниц [ править ]
- ^ Оппенгейм и Шафер, [1] стр. 147 (4.20), стр. 694 (10.1), а также Прандони и Веттерли, [2] стр. 255, (9.33), где: поэтому Также и
- ^ Оппенгейм и Шафер, [1] стр. 551 (8.35), а также Прандони и Веттерли, [2] стр. 82, (4.43). С определениями : и это выражение отличается от ссылки в раз потому что они потеряли его при переходе от 3-го шага к 4-му. В частности, DTFT в § Таблица дискретных преобразований Фурье имеет фактор, который в ссылках опущен.
- ^ Оппенгейм и Шафер, [1] стр. 60, (2.169), а также Прандони и Веттерли, [2] стр. 122, (5.21)
Ссылки [ править ]
- ^ Перейти обратно: а б с д и ж г час я дж к Оппенгейм, Алан В .; Шафер, Рональд В .; Бак, Джон Р. (1999). «4.2, 8.4». Дискретная обработка сигналов (2-е изд.). Река Аппер-Сэдл, Нью-Джерси: Прентис-Холл. ISBN 0-13-754920-2 .
образцы преобразования Фурье апериодической последовательности x[n] можно рассматривать как коэффициенты DFS периодической последовательности, полученные путем суммирования периодических реплик x[n].
- ^ Перейти обратно: а б с Прандони, Паоло; Веттерли, Мартин (2008). Обработка сигналов для связи (PDF) (1-е изд.). Бока-Ратон, Флорида: CRC Press. стр. 72, 76. ISBN. 978-1-4200-7046-0 . Проверено 4 октября 2020 г.
коэффициенты DFS для периодизованного сигнала представляют собой дискретный набор значений для его DTFT.
- ^ Рао, Р. (2008). Сигналы и системы . Prentice-Hall Of India Pvt. Ограничено. ISBN 9788120338593 .
- ^ «Оценка спектральной плотности мощности периодограммы - периодограмма MATLAB» .
- ^ Гумас, Чарльз Константин (июль 1997 г.). «Окно-предполагаемое БПФ обеспечивает широкий динамический диапазон и разрешение» . Новости индивидуальной инженерии и приборостроения : 58–64. Архивировано из оригинала 10 февраля 2001 г.
{{cite journal}}
: CS1 maint: bot: исходный статус URL неизвестен ( ссылка ) - ^ Крошер, RE; Рабинер, Л.Р. (1983). «7,2». Многоскоростная цифровая обработка сигналов . Энглвуд Клиффс, Нью-Джерси: Прентис-Холл. стр. 313–326. ISBN 0136051626 .
- ^ Ван, Хун; Лу, Юсинь; Ван, Сюэган (16 октября 2006 г.). «Канальный приемник с банком фильтров WOLA». 2006 Международная конференция CIE по радиолокации . Шанхай, Китай: IEEE. стр. 1–3. дои : 10.1109/ICR.2006.343463 . ISBN 0-7803-9582-4 . S2CID 42688070 .
- ^ Лайонс, Ричард Г. (июнь 2008 г.). «Хитрости DSP: создание практичного анализатора спектра» . ЭЭ Таймс . Проверено 20 февраля 2020 г. Однако обратите внимание, что он содержит структуру взвешенного перекрытия-добавления с меткой ссылки , которая неправильно переходит к методу Overlap-add .
- ^ Перейти обратно: а б Лиллингтон, Джон (март 2003 г.). «Сравнение архитектур широкополосного каналирования» (PDF) . Даллас: Международная конференция по обработке сигналов. п. 4 (рис. 7). S2CID 31525301 . Архивировано из оригинала (PDF) 8 марта 2019 г. Проверено 6 сентября 2020 г.
Метод «Перекрытие и добавление веса», или WOLA, или его подмножество, «Многофазное ДПФ», становится все более популярным и, безусловно, очень эффективным там, где требуются большие, высококачественные банки фильтров.
- ^ Перейти обратно: а б Лиллингтон, Джон. «Обзор методов набора фильтров — радиочастотные и цифровые» (PDF) . Armms.org . Остров Уайт, Великобритания: Libra Design Associates Ltd., с. 11 . Проверено 6 сентября 2020 г.
К счастью, существует гораздо более элегантное решение, как показано на рисунке 20 ниже, известное как многофазное или WOLA (взвешивание, перекрытие и сложение) БПФ.
- ^ Хохгюртель, Стефан (2013). «Эффективные реализации широкополосных БПФ-спектрометров высокого разрешения и их применение для исследования линий центра галактики APEX» (PDF) . hss.ulb.uni-bonn.de . Бонн: Рейнский Боннский университет имени Фридриха Вильгельма. стр. 26–27 . Проверено 6 сентября 2020 г.
Чтобы выполнить M-кратный WOLA для N-точечного ДПФ, M·N реальных входных выборок a j сначала умножаются на оконную функцию w j того же размера.
- ^ Ченнамангалам, Джаянт (18 октября 2016 г.). «Техника банка многофазных фильтров» . Группа компаний КАСПЕР . Проверено 30 октября 2016 г.
- ^ Даль, Джейсон Ф. (6 февраля 2003 г.). Методы сглаживания времени для оценки спектра (доктор философии). Университет Бригама Янга . Проверено 31 октября 2016 г.
- ^ Линь, Юань-Пей; Вайдьянатан, П.П. (июнь 1998 г.). «Подход окна Кайзера для проектирования прототипов фильтров косинусно-модулированных наборов фильтров» (PDF) . Письма об обработке сигналов IEEE . 5 (6): 132–134. Бибкод : 1998ISPL....5..132L . дои : 10.1109/97.681427 . S2CID 18159105 . Проверено 16 марта 2017 г.
- ^ Харрис, Фредерик Дж. (24 мая 2004 г.). «9». Многоскоростная обработка сигналов для систем связи . Река Аппер-Сэддл, Нью-Джерси: PTR Prentice Hall. стр. 226–253. ISBN 0131465112 .
- ^ Харрис, Фредрик Дж. (январь 1978 г.). «Об использовании Windows для гармонического анализа с дискретным преобразованием Фурье» (PDF) . Труды IEEE . 66 (1): 51–83. Бибкод : 1978IEEP..66...51H . CiteSeerX 10.1.1.649.9880 . дои : 10.1109/PROC.1978.10837 . S2CID 426548 .
- ^ Перейти обратно: а б с д и ж г час я дж к л м н тот п д р Проакис, Джон Г.; Манолакис, Дмитрий Г. (1996). Цифровая обработка сигналов: принципы, алгоритмы и приложения (3-е изд.). Нью-Джерси: Прентис-Холл Интернэшнл. Бибкод : 1996dspp.book.....P . ISBN 9780133942897 . sAcfAQAAIAAJ.
- ^ Рабинер, Лоуренс Р .; Голд, Бернард (1975). Теория и применение цифровой обработки сигналов . Энглвуд Клиффс, Нью-Джерси: Prentice-Hall, Inc., с. 59 (2,163). ISBN 978-0139141010 .
Дальнейшее чтение [ править ]
- Порат, Боаз (1996). Курс цифровой обработки сигналов . Джон Уайли и сыновья. стр. 27–29 и 104–105. ISBN 0-471-14961-6 .
- Зиберт, Уильям М. (1986). Цепи, сигналы и системы . Серия Массачусетского технологического института по электротехнике и информатике. Кембридж, Массачусетс: MIT Press. ISBN 0262690950 .
- Лайонс, Ричард Г. (2010). Понимание цифровой обработки сигналов (3-е изд.). Прентис Холл. ISBN 978-0137027415 .