Лайонс Р. Цифровая обработка сигналов. Второе издание. Пер. с англ. (2006) (1095938), страница 89
Текст из файла (страница 89)
Входной сигнал усредняющего филь 1-я 1И-2)-я 1»г 1) „ 2-я Выходной сигнал усредняющего фильтра Рис. 11.11. Устройство усреднения, изображенное как КИХ-фильтр Вопрос, на который мы должны ответить здесь, формулируется следующим образом: «Какую АЧХ имеет обобщенный Ю-точечный усредняющий фильтру» Мы можем вычислить АЧХ по формуле (6-28) при всех а(Й) = О, которая дает частотную характеристику КИХ-фильтра с Аг ответвлениями. В этом выражении мы должны присвоить всем коэффициентам от Ь(О) до Ь()т' — 1) значения 1/)т'и вычислить модуль Нагл(го) в диапазоне нормированных круговых частот О < ог < л. Этот диапазон соответствует диапазону реювных частот О < / < ~, /2 (где 1', — частота дискретизации в Гц).
Но проще будет вспомнить раздел 5.2, в котором указывалось, что мы можем вычислить частотную характеристику КИХ-фильтра как ДПФ последовательности коэффициентов фильтра. В этом случаем мы можем использовать программу М-точечного БПФ для преобразования последовательности Аг коэффициентов, равных 1/Аг.
Конечно, М должно быть больше Аг, чтобы стала заметной форма частотной характеристики вида яп(х)/х. Следуя этой процедуре и используя программу 128-точечного БПФ, получаем частотные характеристики У-точечного фильтра усреднения при разных значениях Н, приведенные на рисунке 11.12. Чтобы придать этим графикам болыпе наглядности, частотная ось на рисунке определена через 1, в отсчетах/секунду.
430 Глава 11. Ус еднение сигналов о.в о.в 0.4 о.г о о г /4 гув Рис. 11.12. Амплитудно-частотная характеристика И-точечного усредняющего фильтра в зависимости от И 11.5. Экспоненциальное усреднение Существует вид усреднения во временной области, который используется в некоторых приборах измерения мощности — он называется экспоненциальным усреднением 18 - 11]. Этот метод обеспечивает подавление шума путем умножения входного отсчета на константу и сложения этого произведения с предыдущим значением среднего, умноженным на дополнение константы до единицы.
На словах это кажется запутанным, но закон функционирования устройства экспоненциального усреднения представляет собой простое выражение у(л) = ах(л) + (1 — а) у(л — 1), (11-22) где у(л) — текущий выходной отсчет устройства усреднения, у(л-1) — предыдущий отсчет среднего и а — постоянный весовой коэффициент. Процесс, описываемый уравнением (11-22) реализуется так, как показано на рисунке 11.13.
Преимущество экспоненциального усреднения в том, что оно требует всего одну ячейку памяти для хранения значения у(л-1) до получения следующего входного отсчета х(л). Название этого устройства происходит от его импульсной характеристики. Предположим, что на вход устройства усреднения поступает длинная последовательность нулей и в момент времени 1-0 мы подаем на его вход один единичный отсчет. Затем входная последовательность снова представляет собой последовательность нулей. Выходная последовательность устройства усреднения при а - ОА показана на рисунке 11.14.
В момент г = О входной отсчет умножается на а, следовательно, выходной отсчет равен 0.4. На следующем такте входной отсчет равен нулю, а старое значение ОА умножается на (1 — 0.4), или 0.6, давая значение выходного отсчета 0.24. На следующем такте входной отсчет опять равен нулю, старое значение, равное 0.24, умножается на 0.6 и дает выходной отсчет, равный 0.144. Дальше выходная последовательность, или импульсная характеристика, устройства экспоненциально уменьшается вследствие последовательных умножений на 0.6. 431 11.5. Экспоненциальное с еднение у(л) 1-а Рис. 11.13. Устройство экспоненциального усреднения Выходной сигнал усредняющего фильтра 0.40 к 0.35 0.30 0.25 0.20 0.15 0.10 0.05 0 и 4м 0 1 2 3 4 5 5 7 8 9 Время Рис. 11.14.
Импульсная характеристика устройства экспоненциального усредне- ния приа = 0.4 Полезной особенностью устройства экспоненциального усреднения является то, что оно позволяет изменять степень подавления шума путем изменения значения весового коэффициента. Если а=1, входные отсчеты не ослабляются, предыдущие отсчеты среднего не учитываются, и никакого усреднения не происходит.
В этом случае изменения входной последовательности передаются непосредственно на выход. Когда а уменьшается по величине, входные отсчеты ослабляются, предыдущие выходные отсчеты начинают влиять на значение текущего выходного отсчета. Эти предыдущие отсчеты представляют собой сумму экспоненциально взвешенных входных отсчетов, и такое суммирование сглаживает зашумленный сигнал. Чем меньше а, тем сильнее подавляется шум. С другой стороны, чем меныпе а, тем медленнее устройство реагирует на изменения входного сигнала. Мы можем увидеть это, посмотрев на переходные характеристики устройства экспоненциального усреднения при разных значениях а, приведенные на рисунке 11.15'.
Итак, мы должны найти компромисс. Чем больше подавление шума, тем медленнее устройство усреднения будет реагировать на изменения входного сигнала. Мы видим по рисунку 11.15, что при уменьшении а, позволяющем лучше 1 Переходная характеристика представляет собой выходную последовательность устройства усреднения, когда на его вход подается последовательность, состоящая из нулей, за которой следует последовательность, состоящая из единиц.
432 Глава 11. Ус еднение сигналов подавлять шум, для установления выходного сигнала требуется больше времени. Некоторые производители измерительных приборов используют хитроумный прием для достижения компромисса между подавлением шума и временем установления выходного сигнала. В начале измерения они используют большое значение а, так что выходной сигнал сразу принимает ненулевое значение. Затем, с течением времени, значение а уменьшается для достижения требуемого подавления шума. гнал усредняющего Фильтра 1.0 0.9 0.8 2 и ° Г а 01 и' оа 0.8 0.5 0.4 0.З 0.2 оп и 0 г 4 В В 10 12 14 гв 1В Время Рис.
11.15. Выходной сигнал устройства экспоненциального усреднения при подаче на его вход ступенчатого сигнала в момент времени ~0 как функция а Как было показано в [9, 10], уменьшение дисперсии шума на выходе устройства усреднения, как функции весового коэффициента а описывается выражением дисперсия выходного шума/дисперсия входного шума = = 5ЛгЯ „= а,/(2 — а) . (11-23) Следовательно, уменьшение мощности шума в децибелах определяется выражением 5%Я Ив 10108101а/(2 — а)] .
(11-24) На рисунке 11.18 приведен график выражения (11-24), который иллюстрирует соотношение между уменьшением шума и временем реакции устройства усреднения. Чтобы продемонстрировать способность устройства экспоненциального усреднения уменьшать мощность шума, на рисунке 11.17 показан выходной сигнал при подаче на вход устройства косинусоиды, смешанной с шумом.
Весовой коэффициент а линейно уменьшается, начиная со значения 1.0, до значения 0.1 на 180-м отсчете входной последовательности. Обратите внимание на уменьшение шума при уменьшении а. Однако, амплитуда косинусоиды также уменьшается при маленьких значениях а. 11.5. Экспоненциальное с еднение Снижение мощности шума в дБ 15 Более быстрая реакция Более медленная реакция .. я гя инт ~и 1г оя 0 5 т е 5 4 з 'г 1 а Рис. 11.16. Уменьшение мощности шума на выходе устройства экспоненциального усреднения как функция весового коэффициента а а = .45 а = .1 1 45 90 135 150 225 Отсчеты (время) Рис.
11.17. Уменьшение выходного шума устройства экспоненциального усреднения при уменьшении а Читатель может обнаружить, что реализация устройства экспоненциального усреднения на рисунке 11.13 представляет собой БИХ-фильтр с одним ответвлением (12]. Действительно, это так и есть, и мы можем определить его частотную характеристику.
Для этого вспомним общее выражение для частотной характеристики БИХ- фильтра (6-28), приведенное в главе 6, которое мы повторяем здесь: Нггя(то) = [ ~~к~ Ь(Ь)сое(Ьто) — у,~к Ь(Ь)яп(Ьто)1/ (11-25) /[1 — ~т а®соз(Ьш) +у~~> аЯ)яп(Ьш)) 5=1 5=1 По сравнению с рисунком 6-18 мы должны модифицировать (11-25), положив Ж = 0 н М = 1, так что 434 Глава 11. Ус еднение сигналов Н (в) — [Ь(О) ° сох(Ов) — 1ЦО) ° яп(Ов)]/ /[1 — а(1) ° соз(1в) + 1а(1) ° з)п(1в)] . (11-26) Теперь, при о(О) = а и а(1) = 1 — а, частотная характеристика устройства экспоненциального усреднения приобретает вид комплексного выражения Н (в) = а/[1 — (1 — а) ° сов(в) + 1(1 — а) ° з)п(в)] .