Лайонс Р. Цифровая обработка сигналов. Второе издание. Пер. с англ. (2006) (1095938), страница 34
Текст из файла (страница 34)
Таким образом, мы можем переписать уравнение для выходного отсчета усредняющего фильтра (5-4) в более общем виде у(п) = Ь(4)х(п — 4) + Ь(3)х(п — 3) + Ь(2)х(п — 2) + Ь(1)х(п — 1) + Ь(0)х(п) = 4 =ХЬ(Ь)х( -Ь), (5-5) к-0 где все коэффициенты Ь(0) — Ь(4) равны 1/5. Выражение (5-5) представляет собой компактный способ описания структуры фильтра, показанной на рисунке 5.4 и процесса вычисления, иллюстрируемого рисунком 5.5. 170 )лааа5. Фильт ысимп льснойха акте истикойконечнойдлины Лорита в«одиой отому, «ГО) Сумма порви«пити проитввдвиий П д (а) Вычислить у(4) Выходной сигнал у(л) Ствмиоиврмтв у а и а а а О 2 4 Сумма втори«пыи проптввдтпий а« д д д 0 д ', П д (Ь) Вычислить у(О) тт — ь Выходной сигнал у(л) апа а ° а О г (с) а а и а а О 2 4 д 0 0 Д Д П д д О 2 4 Литва вмаиой отп мт «Гву д д д а а д д и П д д Вычислить у(а) О 2 4 Выходной сигнал у(л) (е) а а а ° а Рис.
6.6. Свертка в усредняющем фильтре: (а) первые пять входных отсчетов, выровненные относительно неподвижных коэффициентов фильтра, индекс и = 4; (Ь) сдвиг входных отсчетов вправо, индекс и = 5; (с) индекс и = 6; (с)) индексп=7;(е) индексп=8 Сделаем еще одни шаг вперед и запишем п-й отсчет КИХ-фильтра с и ответвлениями: М вЂ” 1 у(п) =,~„)г()г)х(п-)г), )г-а (5-6) Итак, мы получили то, что хотели. Выражение (5-6) и есть печально известное уравнение свертки в приложении к цифровым КИХ-фильтрам. Начинающие в области цифровой обработки сигналов часто испытываю трудности в усвоении понятия свертки. Эти трудности совсем не обязательны. Уравнение 00 д д д 0 д д Вычислить у(О) и су Выходной сигнал у(л) Вычислить у(т) Ч:> Выходной силтад у(л) 5.2. Све тка вКИХ- ильт ах 171 (5-6) представляет собой просто некоторое количество умножений с последующим суммированием произведений.
В действительности процесс достаточно прост. Мы просто инвертируем порядок следования отсчетов входной последовательности и сдвигаем инвертированную последовательность относительно коэффициентов фильтра, как показано на рисунке 5.5. Для каждого нового входного отсчета мы суммируем произведения и получаем выходной отсчет. Задержимся на минуту и введем новое понятие, которое важно хорошо усвоить — импульсную характеристику. Импульсная характеристика фильтра — это именно то, что значит ее название, это выходная последовательность фильтра во временной области при подаче на вход фильтра единственного отсчета, равного единице (единичного импульса), которому предшествуют и за которым следуют нулевые отсчеты.
Рисунок 5.6 иллюстрирует это понятие точно таким же образом, как рисунок 5.5 иллюстрирует получение выходной последовательности фильтра. Левая часть рисунка 5.6 показывает соответствие коэффициентов фильтра, черных квадратиков, отсчетам входного импульса, отмеченным белыми (пустыми) квадратиками. На рисунках 5.6 (а) — 5.6 (е) мы сдвигаем входные отсчеты вправо и на каждом шаге вычисляем выходной отсчет фильтра с помощью выражения (5-4).
Выходные отсчеты в правой части рисунка 5.6 представляют собой импульсную характеристику фильтра. Обратите внимание на важную особенность: импульсная характеристика КИХ-фильтра идентична последовательности пяти коэффициентов. По этой причине термины коэффициенты КИХ-фильтра и импульсная характеристика являются синонимами.
Следовательно, когда кто-то говорит об импульсной характеристике некоторого КИХ-фильтра, он также говорит и о его коэффициентах. Возвращаясь к усредняющему фильтру, вспомним, что все коэффициенты, или отсчеты импульсной характеристики, Ь(0) — Ь(4) равны 1/5. Оказывается, качество нашего фильтра можно улучшить, если использовать неравные значения коэффициентов. Под качеством фильтра мы понимаем то, насколько хорошо фильтр пропускает требуемые сигналы и насколько хорошо он подавляет нежелательные сигналы. Мы оцениваем качество фильтра по форме его частотной характеристики, которую мы получаем посредством использования свойства свертки линейных систем. Для описания этого понятия повторим (5-6), используя сокращенную запись у(п) = Ь(/г) « х(п), (5-7) где символ «обозначает свертку.
(Формула 5-7 читается так: «у от п равно свертке Ь от Ь их от пм) Свойство свертки в приложении к КИХ-фильтрам формулируется следующим образом: ДПФ свертки импульсной характеристики фильтра и входной последовательности равно произведению спектра входной последовательности и ДПФ импульсной характеристики. Идея, которую мы стараемся донести до вас, состоит в том, что, если две последовательности во временной области Ь(Ь) и х(п) имеют ДПФ Н(т) и Х(т) соответственно, то ДПФ последовательности у(п) = Ь(Ь)» х(п) будет равно Н(т) ° Х(т). Формулируя более компактно, мы представляем это утверждение в форме выражения дпо у(п) = Ь(Ь) «х(п) Н(т) Х(т) (5-8) ОИПФ 172 Глава 5.
Филвт ы с имп льснойха акте нстикой конечной длины Сумма п«аамх пвм В О й Вычислить у(4) ы> Импульсная характеристика, Ип) Импуп а й 1 О «хайяма аппп (д) В и В В В Супцпа«пр«в юзффуммвпм О 2 4 Сумма« самки ум паамэкеапммй Вычислить у(Б) Импульсная характеристика, у(п) 1О (Ь) е В ° В О 2 4 Вычислить у(6) В и Импульсная характеристика, у(п) (с) В В В М В О 2 4 1О Вычислить у(у) В В В В Импульсная характеристика, у(п) ° В М В В О 2 4 Вычислить у(8) В в В ° В Импульамвя характеристика, у(п) (е) Рис. 5.6.
Свертка коэффициентов фильтра и входного импульса для получения импульсной характеристики фильтра: (а) единичный отсчет импульса выровнен с первым коэффициентом фильтра, индекс и = 4; (Ь) импульс сдвинутвправои индексп= 5;(с) индексп=б;(с() индексп=7;(е) индекс п=8 Здесь ОДПФ обозначает обратное ДПФ, и (5-8) показывает, что две последовательности, образованные в результате выполнения операций Ь()у) «2(л) и Н(т) Х(т), связаны парой преобразований Фурье. Таким образом, вычисляя ДПФ от свертки )т(к)»х(л), мы получаем произведение Н(т) «Х(т), т. е. спектр выходного сигнала нашего фильтра У(гл).
Аналогично, мы можем получить свертку уу(к) «х(л), вычисляя обратное ДПФ произведения Н(лу) «Х(т). Из (5-8) можно сделать важное заключение о том, что свертка во временной области эквивалентна умножению в частотной области. Чтобы помочь вам понять этот принцип, на рисунке 5.7 наглядно представлено соотношение между сверткой во временной области и произведением в частотной области.
Операция свертки для линейных систем обсуждается подробнее в разделе 5.9. Мы настоятельно рекомендуем начинающему читателю просмотреть этот материал, чтобы получить общее представление о том, как и когда операция сверки может быть использована для анализа цифровых фильтров. Выражение (5-8) и соотношения на рисунке 5.7 подсказывают нам, что необходимо сделать, чтобы получить частотную характеристику некоторого КИХ-фильтра. Произведение Х(т) Н(т) — это ДПФ выходного сигнала фильтра.
Поскольку 5.2. Све ткавКИХ- нльт ех 173 Х(т) является ДПФ входной последовательности, отсюда следует, что частотная характеристика фильтра есть Н(т), ДПФ импульсной характеристики фильтра а(тг) . Возвращаясь опять к нашей задаче, мы можем определить частотную характеристику усредняющего фильтра, вычислив ДПФ последовательности коэффициентов фильтра (импульсной характеристики) в (5-4). Если мы возьмем пять коэффициентов гт(1г), равных 1/5, и дополним их 59 нулями, мы получим последовательностгп изображенную на рисунке 5.8 (а). Выполнив 64-точеное ДПФ этой последовательности и пронормировав модули ДПФ, мы получаем амплитудно-частотную характеристику. фильтра ~Н(т) ~, показанную на рисунке 5.8 (Ъ), и фазо-частотную характеристику, показанную на рисунке 5.8 (с)'. Мы видим, что Н(т) — наша старая знакомая, функция яп(х)гтх из раздела 3.13. вх си аиаппаа впасть пр астстпаа бпасть Рис.
Б 7. Соотношения свертки в приложении к цифровым КИХ-фильтрам Теперь установим связь отсчетов дискретной частотной характеристики, показанной на рисунках 5.8 (Ь) и 5.8 (с), с физической размерностью частоты дискретизации у',. Из раздела 3.5 и нашего опыта манипуляций с ДПФ мы знаем, что отсчет ш = Ж,'2 = 32 в этом случае соответствует частоте заворота, равной половине частоты дискретизации, /;,'2. Учитывая это, мы можем преобразовать частотную 1 Здесь мы говорим «импульсная характеристика» вместо «коэффицентьт», потому что эта концепция применима также и к БИХ-фильтрам. Частотная характеристика БИХ- фильтра также равна ДПФ импульсной характеристики.