Лайонс Р. Цифровая обработка сигналов. Второе издание. Пер. с англ. (2006) (1095938), страница 27
Текст из файла (страница 27)
Это значит, что частотная характеристики отдельного бина при подаче на вход ДПФ действительной синусоиды описывается функцией з)пс, которая определяется выражениями (3-74) — (3-79). 3.17. Интерпретация ДПФ Теперь, когда мы кое-что узнали о ДПФ, самое время убедиться в том, что мы правильно понимаем, что в действительности представляет собой ДПФ, и постараться избежать обычных заблуждений относительно его поведения. В литературе по ЦОС мы найдем разделы, посвященные непрерывному преобразованию Фурье, рядам Фурье, дискретно-временному преобразованию Фурье, дискретному преобразованию Фурье и периодическим спектрам. Чтобы не запутаться во всех этих понятиях, требуется приложить определенные усилия, особенно когда вы читаете или слышите что-то вроде «ДПФ предполагает, что его входная последовательность периодична во времени».
(Вы удивляетесь, почему это так, поскольку ничего не стоит вычислить ДПФ апериодической последовательности.) Такое замечание в лучшем случае сбивает с толку, т. е. ДПФ не требует подобных предположений. Далее я постараюсь изложить свое понимание природы периодичности последовательностей во времени и по частоте. Рассмотрим непрерывный сигнал бесконечной длительности, содержащий единственный импульс конечной длительности, показанный на рисунке 3.46 (а).
Модуль непрерывного преобразования Фурье (НПФ) такого сигнала представляет собой непрерывную функцию частоты Хг(ге). Если этот импульс можно описать аналитически (математической формулой), то функцию Х1(ге) можно также найти аналитически, взяв интеграл Фурье. (Весьма вероятно, что вы делали это в 3.17. Инте п етация ПФ 133 качестве домашней или контрольной работы.) Непрерывная частота го измеряется в радианах в секунду.
Если мы возьмем НПФ бесконечного по длительности сигнала, представляющего собой периодическую последовательность импульсов, показанную на рисунке ЗА6 (Ь), мы получим линейчатый спектр, известный как ряд Фурье' Х~(го). Линии этого спектра имеют бесконечно малую ширину, а Х~(со) хорошо определен в промежутках между этими линиями, т. к. Х~(го) является непрерывной функцией частоты. (а) т-З т-1 ты т+З т-3 т-1 т+1 т+Э (с) т-З т-1 т+1 т+з т-3 т-1 т+1 анэ частота Х Рис. 3.46. Определение модуля т-го бина И-точечного ДПФ: (а) действительная последовательность х,(и) имеет й = и1-2 5 периодов на анализируемом интервале; (о) действительная последовательность х (и) на анализируемом интервале имеет )с = т-1.5 периодов; (с) действительная последовательность х (и) на анализируемом интервале имеет )с = ги периодов; (с)) частотная характеристика отдельного бина ДПФ ги = )с ! Точнее, спектр представляет собой коэффициенты ряда Фурье — (прим.
перев.). 134 Глава 3. Диск егное и еоб азование Ф ье й ческий НПФ (а) Время Частота Непрерывный и периодический Линейчатый спектр — НП (Ь) Частота Время Апериодическая дискретная последовательность бесконечной длины ° ° ° ° ° ° ° ° х(п) Непрерывный и периодический ДВПФ— (с) еееееееее ° - — . — - - ее«вен аеВремя -е, 0 х,(п) (-г,) ДПФ Дискретный дискретный и периодический Обратное йь. и периодический ° ° ° ° ° 'ек е'кека ее «В )(ПФ ' Х чп' х (и) Х,(пт) . 2 — НПФ-~ ° ° (б) ° ° к Время -/, О Г, Частота и, Частота (т,) Рис.
3.46. Сигналы и последовательности во временной области и модули их преобразований в частотной области (3-80) и-- Чтобы проиллюстрировать ДВПФ, предподожим, что у нас есть последовательность видах,(л) = (0.75)" при п ) О. Бе ДВПФ имеет вид м м Хо(пт) = „)' 0.75"е ) "=. „Я (0.75е х )" (3-81) и= — а Выражение (3-81) представляет собой геометрическую прогрессию (см. приложение В) сумма которой равна Рисунок 3.46 (Ь) показывает пример непрерывной периодической функции, спектр которой содержит ряд дискретных компонентов.
Вы можете считать ряд Фурье Х2(то) дискретизированной версией непрерывного спектра на рисунке 3.46 (а). Это соотношение между временным и частотным представлением сигнала х2(г) и Хз(от) показывает, что периодической функции в одной области соответствует дискретная по своей природе функция в другой области. Далее, рассмотрим бесконечную по длительности дискретную последовательность х(п), содержащую несколько ненулевых отсчетов, показанную на рисунке 3.46 (с).
Мы можем вычислить НПФ х(п), которое описывает спектр в виде непрерывной функции Хз(со). Этот непрерывный спектр называется дискретно-временным преобразованием Фурье (ДВПФ), которое определено как (см. стр. 48 в 16) ОО Х(ш) =')'х(п)а-) и 136 3.17. Инте п етация ПФ Х„(4О) = 1/(1 — 0.75е Хв) = еУ",'(е~"-0.75) .
(3-82) Спектр Х,(в) непрерывен и периодичен с периодом 2я, его модуль показан на рисунке 3.47. Это пример дискретной функции времени, имеющей периодический спектр. Для особо любознательного читателя мы можем проверить периодичность ДВПФ, используя целое число Й в следующем выражении: Х(„+ 2пь) = ~~~ х(п)е-1(и ~2"л)л = ) ( ) -ллл -12плл = (3-83) л = — л~ = ~~)' х(п)е л =Х(щ) л — ю потому что е 12"ьл = 1 для целых значений й. Хз(в) на рисунке 3.46 (с) также имеет период 2л, если представить его как функцию круговой частоты ил = 2тг/;, где частота/; представляет собой величину, обратную периоду отсчетов х(п). Непрерывная периодическая функция Хз(м) представляет собой то, что мы хотели бы иметь на практике, но мы не можем этого добиться.
Мы используем компьютеры, и как это ни грустно, не можем выполнять анализ непрерывных сигналов из-за дискретной природы самих компьютеров. Любые операции обработки данных выполняются над дискретными числами, хранящимися в памяти компьютера, и вследствие этого все сигналы во временной области и все частотные спектры представляют собой дискретные последовательности, полученные в результате дискретизации. Следовательно, НПФ и обратное НПФ последовательностей, с которыми мы работаем, будут периодическими.
= з, з )с 2 О -л О л 2л бл ю 4л Рис. 3.47. Модуль ДВПФ ~Х„(и) ~ Преобразования, показанные на рисунках ЗА6 (а) — (с) представляют собой результат математических выкладок на бумаге. В компьютере, использующем дискретные последовательности конечной длины, мы можем только аппроксимировать НПФ (ДВПФ) последовательности х(п) бесконечной длины, показанное на рисунке 3.46 (с). Эта аппроксимация называется дискретным преобразованием Фурье (ДПФ), и это единственный инструмент Фурье, который нам доступен. Вычисляя ДПФ последовательности хз(п), которая является конечной частью х(п), мы получаем дискретные периодические отсчеты Х~(т), показанные на рисунке 3.46 (д).
Обратите внимание на то, что Х~(т) есть дискретизированная 1ЗВ Глава 3. иск егноел еоб взование Ф ье версия непрерывного периодического спектра Хэ(а|). Однако спектр Хг(т) в точности равен НПФ периодической последовательности х~(п) на рисунке ЗАб (д). Так что, когда кто-нибудь говорит: «ДПФ предполагает, что входная последовательность периодична во времени», — он имеет ввиду, что ДПФ равно непрерывному преобразованию Фурье (которое называется ДВПФ) периодической дискретной последовательности. Итог всего сказанного заключается в следующем: если функция периодична, ее прямое/обратное ДВПФ будет дискретным; если функция дискретна, ее прямое/обратное ДВПФ будет периодическим.
Библиография 1. Вгасе|че11, К. «ТЬе РопНег Тгапз(огш», Ес(епсфс Атепсап, )ппе 1989. 2. Зсгпй, 1). А Сопсгхе НЫогу о/Майетаг(сз, Почег РцЫ!сас(опз 1пс., Хе|ч 'г'ог1с, 1967, р. 142 (неоднократно издавались русские переводы, например: Стройк Д. Я. Краткий очерк истории математики. М. Наука, 1990, 256 с.). 3. %'11!саша, С. 5. Рея8пт8 И8(га1 Р1(гегх. бесс!оп 8.6, Ргепсссе-На!1, Еп81е|чоод СИ(з, Хе|ч !егзеу, 1986, р.
122. 4. Ргеэз, ЪЪ'., ес а!. 7|/итепса! Кес(рез — Т7|е Агг о/5с(епгфс Сотриг(пй. СашЬПдйе Пшчегз(су Ргезз, 1989, р. 426. 5. ОесЫп11, Х. С., апд г'ачцх, П. «Зоше Хоче! ЪЪ'!вдовая апд а Сопсгэе Тпсопа1 Сошраг!хоп о(ЪЪ|1пдо|ч Рапп11ез,» 1ЕЕЕ Тгапя оп Асоизс 5реес7|, апг75(япа! Ргос., Ъ'о!. А85 Р-26, Хо. 6, ПесегпЪег 1978. (Кстати, на странице 505 этой статьи фраза «такой, что ЪЪг(Г) ) 0|| 1» означает, что функция ЪЪг(Г) никогда не бывает отрицательной. Символ Ъ| значит «для всех».) 6.
О'Роппе!Ц. «1ооЬ(пб ТЬгоц8Ь сЬе К!8Ьс ЪЪ'1пдочч 1шргочез брессга! Апа1уяз,» ЕТ)1ч', ХочешЪег 1984. 7. Ка(зег, !. К «О!ясса! Р!!сегз,» ш 5узгет Апа1уяз бу 1)18(га1 Сотригег. Е|1. Ъу Г. К Кпо апд !. К Касзег, 1оЬп ЪЪ|11еу апд бонз, Хе|и 'г'огас, 1966, рр. 218-277. 8. КаЪ|пег, 1.. К., апд Оо!д, В.
77|е Т7|еогуапг7Арр(ссагюпо/018(га(518па1Ртосея(п8. Ргепгссе-На!1, Еп81еч оод СИ(з, Хе|ч )егзеу, 1975, р. 88 (есть русский перевод: Рабинер Л., Голд Б. «Теория и применение цифровой обработки сигналов», Мс Мир, 1978, доступен по адресу Ьсср://йео81п.пагод.гц/агЬ(ч/дзр/дэрЗ.Ьсш.