Лайонс Р. Цифровая обработка сигналов. Второе издание. Пер. с англ. (2006) (1095937), страница 41
Текст из файла (страница 41)
Как подробно описано в разделе 5.2, мы выполняем дискретную фильтрацию с помощью КИХ-фильтра, вычисляя свертку входной последовательности х(п) с импульсной характеристикой фильтра л()!), а в случае КИХ-фильтра отсчеты его импульсной характеристики оказываются равными коэффициентам фильтра. Результат этой свертки представляет собой отфильтрованную последовательность во временной области, спектр которой модифицирован умножением на частотную характеристику фильтра Н(т). В разделе 13.10 описывается остроумная схема эффективной реализации КИХ-фильтра с использованием БПФ для выполнения свертки.
2 ! 4 5 8 1 Рис. 6.42. Соотношения между л(lг), к(lг), Н(гп) и Х(тл), показанными на рисун- ке 5,38, в процессе вычисления свертки ! Как мы увидим в главе 6, коэффициенты фильтра с бесконечной импульсной характеристикой (БИХ) не равны отсчетам импульсной характеристики. о о~ о о,' а о, + о о,' (Ч а а а + о о о о о + + о о сС О3 а 7 ю СМ СО ! а о о с~ о о а а О3 о + ~ Ю ' а а См о СМ <ч а о + о (О тй СЧ о а о а о а 'О о о~ СО СО Й ЬС Ф О о Б О. СО х х .й Х Ш о 4С )с Е Б (~ О) С2 Б Б Ф с о Б Е с х л Х Ф Э О.
Э с СС Б Х Ф СО со со МЪ СО з Ц Ю СО 1 Е ~~ В Е.' , ~)с,', оЕ: Е), Е, 9~< с ц, 4С С~ эС, Е. ~9к,', ~=п ,' 1 '*з„ , 'С~О '! ЗОС о о о о о о о о о мъ и') л + о О СМ СМ о о СО СО ~о о со о + <. о о СО о ! о о о + о а о о о о о 1 о о о о ~ С) 'Ф МЭ СО 5.9. Обобщенное описание диск етной све тки 213 Вследствие дуальности теоремы о свертке мы можем в наших рассуждениях о свертке и произведении как паре преобразований Фурье поменять местами временную и частотную области. Это значит, что подобно (5-31) мы можем написать ДПФ ЫА) ° х(й) Н(т) + Х(т) (5-32) ОЯПФ Следовательно, теорему о свертке можно сформулировать в более общем виде так: Свертка в одной области эквивалентна произведению в другой области.
На рисунке 5АЗ показано соотношение между умножением во временной области и сверткой в частотной области. Соотношение (5-32) представляет собой фундаментальное соотношение, используемое в процессе взвешивания данных с помощью окон для уменьшения утечки спектра, описанного в разделе 3.9. Временная область Па преобра Фур Частотная область Рис. 5.43. Соотношения теоремы о свертке для умножения во временной области 5.9.3. Применение теоремы о свертке Теорема о свертке полезна как инструмент качественного предсказания результатов различных операций в дискретных линейных инвариантных во времени системах.
Например, многие авторы используют теорему о свертке, чтобы показать, почему периодическая дискретизация непрерывных сигналов дает дискретные сигналы с периодическими по частоте спектрами. Рассмотрим действительный непрерывный сигнал, показанный на рисунке 5.44 (а), односторонний спектр которого имеет ширину В. Глава 5. Фнльт ысимп льснойха акте истнкойконечнойдлины 214 -В О В чаоаоаа кккка ка а1 (Ь) В ма 1 к а а а к г~~, ч 2К, ПК, 4 а Врача -21 -т О т 2т чаоаоаа Рис. 5.44. Использование теоремы о свертке для предсказания размножения спектра при периодической дискретизации 11оскольку это действительный сигнал, его спектр, конечно же, симметричен относительно 0 Гц.
(На рисунке 5А4 широкие стрелки, показывающие вправо, обозначают операции преобразования Фурье.) Дискретизация этого сигнала эквивалентна умножению его на периодическую последовательность импульсов, показанную на рисунке 5А4 (Ъ), отсчеты которой равны единице. Если частота дискретизации равна /, отсчетов в секунду, то период дискретизации г, = 1Д, секунд. Результатом этого перемножения является последовательность дискретных импульсов, показанная на рисунке 5-44 (с). Для предсказания того, какие эффекты в частотной области будут вызваны этим перемножением во временной области, мы можем использовать теорему о свертке. Согласно этой теореме, как мы теперь знаем, спектр произведения во временной области должен быть сверткой спектров сомножителей.
Хорошо, мы знаем, каков спектр непрерывного сигнала. А что можно сказать о спектре последовательности импульсов? Было показано, что спектр периодической последовательности импульсов, период которой равен г, секунд, представляет собой периодическую последовательность импульсов в частотной области, расстояние между которыми равно~; Гц, как показано на рисунке 5-44 (Ъ) [301 Теперь нам остается только вычислить свертку двух спектров.
В нашем случае это выполняется просто, поскольку оба спектра симметричны относительно частоты 0 Гц, и зеркальное отображение одного из них относительно 0 Гц будет излишним. Таким образом, мы просто смещаем один спектр относительно другого и 216 59. Обобщенноеописаниедиск етнойсае тки рисуем график их произведения. Свертка спектра исходного сигнала и импульсов в частотной области приводит к тому, что спектр сигнала повторяется через каждые/, Гц, как показано на рисунке 5А4 (с). Это напоминает нам о том, что ДПФ всегда периодично с периодом/, Гц. А вот еще один пример того, как теорема о свертке может прийти на помощь, когда мы пытаемся разобраться в операциях цифровой обработки сигналов. Автор однажды использовал эту теорему, чтобы разобраться в том, почему первый ноль спектра треугольного окна приходится на частоту, которая в два раза превышает частоту первого нуля прямоугольного окна.
Вопрос формулировался следующим образом: «Если прямоугольное окно длительностью Т имеет первый ноль спектра на частоте 1/Т Гц, почему треугольное окно такой же длительности Т имеет первый ноль спектра на частоте 2/Т Гц?» Мы можем ответить на этот вопрос, рассмотрев свертку во временной области. Посмотрите падве прямоугольные функции времени, показанные на рисунках 5А5 (а) и 5А5 (Ъ). Если длительность каждой из них равна Т секунд, их спектры имеют первые нули на частоте 1/Т Гц, как показано на графиках, приведенных на рисунках 5.45 (а) и 5А5 (Ъ).
Мы знаем, что модуль спектра представляет собой абсолютную величину классической функции яп(х)/х . Если мы выполним свертку этих двух прямоугольных функций длительностью Т, мы получим треугольную функцию, показанную на рисунке 5А5 (с). Здесь тоже зеркальное отображение одной из прямоугольных функций относительно нулевого времени не обязательно. Нам нужно просто сдвигать одну функцию относительно другой и вычислять площадь их перекрытия. Наибольшая площадь перекрытия будет при нулевом сдвиге.
Следовательно, свертка будет иметь максимум при нулевом сдвиге по времени, поскольку при этом имеется стопроцентное перекрытие. Если мы будем сдвигать одну из прямоугольных функций в любом направлении, свертка будет линейно уменьшаться до нуля. При сдвиге Т/2 секунд, прямоугольники перекрываются на 50 %. Свертка станет равной нулю при сдвиге на Т секунд, когда перекрытие прямоугольников прекратится. Заметьте, что длительность треугольника, полученного в результате свертки, равна 2Т, и это ключ к ответу на наш вопрос. Поскольку свертка во временной области эквивалентна умножению в частотной области, модуль преобразования Фурье нашей треугольной функции длительностью 2Травен произведению функции ~яп(х)/х(, показанной на рисунке 5А5 (а), на функцию ~яп(х)/х), показанную на рисунке 5.45 (Ь), или (яп(х)/х), как показано на рисунке 5А5 (с). Если треугольная функция длительностью 2Тимеет первый ноль на частоте 1/ТГц, то та же функция длительностью Т должна иметь частоту первого нуля, равную 2/Т Гц, как показано на рисунке 5.45 (б), и это как раз то, что требовалось показать.
Сравнение рисунков 5.45 (с) и 5.45 (д) демонстрирует фундаментальное свойство преобразования Фурье, которое состоит в том, что сжатие функции во времени приводит к растягиванию ее спектра в частотной области. 1 Мы познакомились с функцией яп(х)/х во время обсуждения функций окон в разделе 3.9 и более подробно описали ее в разделе 3.13.
216 Глава 5. Фильт ы с имл льсной ха акте истикой конечнойдлины -ттт О «т Чаотота (Ь) О тт2 Время -Вт О Вт "' ота (В (с) Время -Вт О Вт Чаотота О тт2 Время Рис. 8.45. Использование свертки для того, чтобы показать, что преобразование Фурье треугольной функции имеет первый ноль спектра на частоте, равной удвоенной частоте первого нуля преобразования Фурье прямоугольной функции Библиография 1.
ВЬупКЯ.Я. «Аг(арг(че ПК Ейейп8», ХЕЕЕАЯЯРМадагте, Арп1 1989. 2. ! апас)г(е, А. «АО1арГгве Ейегз ЕпаЫе ЯуВГеп2В Го Тгзс(г Чзг(аГ(опз», М(сгоглаоез О ЯР, ВерсегпЪег 1989. 3. Вп!!ос)г, Я, К. «Н(8Ь ЕгеО(пепсу Ас1арГгве Ейег», Мгсгоатаое /оигла1, Яергеп2Ьег 1990.
Библиот а ия 217 10. ВагЬ|еге, Р. «А МесЬос? ?ог Са1сп!абпб гЬе Свггепс Пзсг!Ьпс!оп о? СЬеЪузЬеч Аггауз», Ртосеес?1пбз ос йе 1КЕ, ч'о1. 40, ?аппагу 1952. 11. Соо?с, С. Е., апс? Вегп?е14, М. Кас?ат 518паЬ, Аеас?еш!с Ргезз, Ь?е|ч Уог?с, 1967, рр. 178-180. 12. Ка!зег, !. К «Р?8!са! ГПсегз», ш Еузгет Апа1уыз Ьу Ргфга? Сотригеп Ес?. Ьу Г. К Кпо апс? !. К Ка?зег,?оЬп чЧПеу апс? Бопз, Хе|ч г от?с, 1966, рр. 218-277. 13. ЖППашз, С. 5.
Резгдп?пд 0181га1 Р11гетз, Ргепске-НаП, Епб!е|чоос! СИ?з, Ь?еьч ?ег- зеу, 1986, р. 117. 14. Нагпз, К !. «Оп |Ье ??зе о?Ж!пс?о|чз ?ог Нагшоп!с Апа1уяз |ч!сЬ сЬе Р!зсгесе Го пбег Тгапз?огш», Ртосеес?1п8з о1'йе ?ЕЕЕ, сто!. 66, Ь?о. 1, !аппагу 1978. 15. МсС1еПап, !. Н., Раг1сз, Т. Ъ", апс1 КаЪ|пег, 1.. К.