Цифровые фильтры в каналах передачи одномерных сигналов

Материал из Национальной библиотеки им. Н. Э. Баумана
Последнее изменение этой страницы: 23:46, 17 ноября 2016.


Задачей обнаружения и распознавания сообщений в КПС является установление факта передачи скрытого сообщения и его идентификации или отдельных конкретных объектов по данным наблюдения и анализа сигналов, несущих сообщения. Поэтому разработка методов математической обработки и интерпретации анализируемых сигналов, а также установления взаимосвязи между физическими свойствами носителя сообщения, способов кодирования и модуляции имеет большое значение для решения указанной задачи.

Цифровая фильтрация распознаваемых сигналов является непременной процедурой при обнаружении и распознавании сообщений. Под фильтрацией будем понимать любое преобразование сигналов, при котором во входной последовательности обрабатываемых данных целенаправленно изменяются определенные соотношения (динамические или частотные) между различными параметрами сигналов. Как известно, преобразование динамики сигналов (и данных, которые несут эти сигналы) осуществляется в системах, являющихся частью каналов передачи сообщений. Системы, избирательно меняющие форму сигналов (амплитудно-частотную или фазово-частотную характеристику), устранение или уменьшение помех, извлечение из сигналов определенной информации, разделение сигналов на определенные составляющие, и т.п. называют фильтрами. Соответственно, фильтры с любым целевым назначением являются частным случаем систем преобразования сигналов.

К основным операциям фильтрации информации относят операции сглаживания, прогнозирования, дифференцирования, интегрирования и разделения сигналов, а также выделение информационных (полезных) сигналов и подавление шумов (помех). Основными методами цифровой фильтрации данных являются частотная селекция сигналов и оптимальная (адаптивная) фильтрация.

Здесь рассматриваются, в основном, методы линейной обработки данных (носителей этих данных - сигналов) линейными дискретными системами. Линейными называют системы, которые осуществляют преобразование линейных комбинаций входных сигналов в суперпозицию выходных сигналов. Принцип реализации линейных систем, аппаратный - в виде специальных микропроцессорных устройств, или алгоритмический - в виде программ на ЭВМ, существенного значения не имеет и определяет только их потенциальные возможности.

В общем случае термином цифровой фильтр называют аппаратную или программную реализацию математического алгоритма, входом которого является цифровой сигнал, а выходом – другой цифровой сигнал с определенным образом модифицированной формой и/или амплитудной и фазовой характеристикой. Классификация цифровых фильтров обычно базируется на функциональных признаках алгоритмов цифровой фильтрации, согласно которому ЦФ подразделяются на 4 группы: фильтры частотной селекции, оптимальные (квазиоптимальные), адаптивные и эвристические. Наиболее изученными и опробованными на практике являются ЦФ частотной селекции.

Цифровые фильтры

Общие понятия

В одномерной дискретной линейной системе связь между входом и выходом (входной и выходной дискретными последовательностями значений сигнала – отсчетами), задается линейным оператором преобразования :

Это выражение отображает краткую запись линейного разностного уравнения:

где - порядковый номер отсчетов,
- интервал дискретизации сигнала,
и - вещественные или, в общем случае, комплексные коэффициенты.

Положим , что всегда может быть выполнено соответствующей нормировкой уравнения , и, принимая в дальнейшем , т.е. переходя к числовой нумерации цифровых последовательностей значений сигналов, приведем его к виду:

Оператор, представленный правой частью данного уравнения, получил название цифрового фильтра (ЦФ), а выполняемая им операция - цифровой фильтрации данных (информации, сигналов). Если хотя бы один из коэффициентов или зависит от переменной , то фильтр называется параметрическим, т.е. с переменными параметрами. Ниже мы будем рассматривать фильтры с постоянными коэффициентами (инвариантными по аргументу).

Основные достоинства цифровых фильтров

Основные достоинства цифровых фильтров по сравнению с аналоговыми:

  • Цифровые фильтры могут иметь параметры, реализация которых невозможна в аналоговых фильтрах, например, линейную фазовую характеристику.
  • ЦФ не требуют периодического контроля и калибровки, т.к. их работоспособность не зависит от дестабилизирующих факторов внешней среды, например, температуры.
  • Один фильтр может обрабатывать несколько входных каналов или сигналов.
  • Входные и выходные данные можно сохранять для последующего использования.
  • Точность цифровых фильтров ограничена только используемой разрядностью отсчетов (длиной слов).
  • Фильтры могут использоваться при очень низких частотах и в большом диапазоне частот, для чего достаточно только изменять частоту дискретизации данных.

Нерекурсивные фильтры

При нулевых значениях коэффициентов уравнение переходит в уравнение линейной дискретной свертки функции с оператором :

Значения выходных отсчетов свертки для любого аргумента определяются текущим и "прошлыми" значениями входных отсчетов. Такой фильтр называется нерекурсивным цифровым фильтром (НЦФ). Интервал суммирования по получил название "окна" фильтра. Окно фильтра составляет отсчет, фильтр является односторонним каузальным, т.е. причинно обусловленным текущими и "прошлыми" значениями входного сигнала, и выходной сигнал не может опережать входного. Каузальный фильтр может быть реализован физически в реальном масштабе времени. При , а также при для фильтра , проведение фильтрации возможно только при задании начальных условий для точек и . Как правило, в качестве начальных условий принимаются нулевые значения, или продление первых отсчетов входных сигналов или его тренда назад по аргументу.

При обработке данных на ЭВМ ограничение по каузальности снимается. В программном распоряжении фильтра могут находиться как "прошлые", так и "будущие" значения входной последовательности отсчетов относительно текущей точки вычислений , при этом уравнение будет иметь вид:

При фильтр называется двусторонним симметричным. Симметричные фильтры, в отличие от односторонних фильтров, не изменяют фазы обрабатываемого сигнала. Так как реакция НЦФ на единичный входной импульс (а равно и на любой произвольный входной сигнал) всегда конечна и ограничена размером окна фильтра, такие фильтры называют также фильтрами с конечной импульсной характеристикой (КИХ-фильтры). Техника выполнения фильтрации не отличается от техники выполнения обычной дискретной свертки двух массивов данных.

Представим, что на одной полоске бумаги выписаны по порядку сверху вниз значения данных (см. рис. 1.1). На второй полоске бумаги находятся записанные в обратном порядке значения коэффициентов фильтра (обозначение для коэффициентов операторов НЦФ является общепринятым). Для вычисления располагаем вторую полоску против первой таким образом, чтобы значение совпало со значением , перемножаем все значения с расположенными против них значениями , и суммируем все результаты перемножения. Результат суммирования является выходным значением сигнала . Сдвигаем окно фильтра - полоску коэффициентов , на один отсчет последовательности вниз (или массив сдвигаем на отсчет вверх) и вычисляем аналогично следующее значение выходного сигнала, и т.д.

Рис. 1.1. Нерекурсивный ЦФ.

Описанный процесс является основной операцией цифровой фильтрации, и называется сверткой в вещественной области массива данных с функцией (оператором) фильтра (массивом коэффициентов фильтра). Для математического описания наряду с формулами применяется также символические формы записи фильтрации:

Сумма коэффициентов фильтра определяет коэффициент передачи (усиления) средних значений сигнала в окне фильтра и постоянной составляющей в целом по массиву данных (с учетом начальных и конечных условий). Как правило, сумма коэффициентов фильтра нормируется к 1. Имеется целый ряд методов обработки данных, достаточно давно и широко известных, которые по существу относятся к методам цифровой фильтрации, хотя и не называются таковыми. Например, методы сглаживания отсчетов в скользящем окне постоянной длительности. Так, для линейного сглаживания данных по пяти точкам с одинаковыми весовыми коэффициентами используется формула:

С позиций цифровой фильтрации это не что иное, как двусторонний симметричный нерекурсивный цифровой фильтр:

Аналогично, при сглаживании данных методом наименьших квадратов (МНК) на основе кубического уравнения:


Это также НЦФ с коэффициентами: .

Пример:

Уравнение НЦФ:

Начальные условия - нулевые.

Входной сигнал – скачок функции (ступень):

Выходной сигнал:

Результат фильтрации приведен на рис. 1.2(А). Проверьте результат (выполните фильтрацию, как это показано на рис. 1.1, с учетом четности фильтра).

Заметим: сумма коэффициентов сглаживающих НЦФ всегда должна быть равна 1, при этом сумма значений массива выходного сигнала равна сумме значений массива входного сигнала. Координатная детальность выходного сигнала ниже входного, резкие изменения входных сигналов "размазываются" по аргументу. Повторите фильтрацию фильтром МНК на основе кубического уравнения. Сравните результаты фильтрации с результатами первого НЦФ (приведены на рис. 1.2)).

Рис. 1.2. Сглаживание МНК в скользящем окне по пяти точкам.

Для операции фильтрации характерны следующие основные свойства:

  • Дистрибутивность: .
  • Коммутативность: .
  • Ассоциативность: .

Фильтрация однозначно определяет выходной сигнал для установленного значения входного сигнала при известном значении импульсного отклика фильтра .

Рекурсивные фильтры

Фильтры, которые описываются полным разностным уравнением , принято называть рекурсивными цифровыми фильтрами (РЦФ), так как в вычислении текущих выходных значений участвуют не только входные данные, но и значения выходных данных фильтрации, вычисленные в предшествующих циклах расчетов. С учетом последнего фактора рекурсивные фильтры называют также фильтрами с обратной связью, положительной или отрицательной в зависимости от знака суммы коэффициентов . По существу, полное окно рекурсивного фильтра состоит из двух составляющих: нерекурсивной части , ограниченной в работе текущими и "прошлыми" значениями входного сигнала (при реализации на ЭВМ возможно использование и “будущих” отсчетов сигнала) и рекурсивной части , которая работает только с "прошлыми" значениями выходного сигнала. Техника вычислений для РЦФ приведена на рис. 1.3.

Рис. 1.3. Рекурсивный ЦФ.
Пример:

Уравнение РЦФ: при

Входной сигнал:

Расчет выходного сигнала:

и т.д.

Выходной сигнал:

Рис. 1.4. Рекурсивная фильтрация.

Из примера можно видеть, что реакция РЦФ на конечный входной сигнал (например, на единичный импульс Кронекера в точке 2), в результате действующей положительной обратной связи, в принципе, может иметь бесконечную длительность (в данном случае с близкими к нулю, но не нулевыми значениями), в отличие от реакции НЦФ, которая всегда ограничена количеством членов (окном фильтра). Фильтры такого типа называют также фильтрами с бесконечной импульсной характеристикой (БИХ-фильтры).

Пример:

Уравнение РЦФ: при

Входной сигнал:

Выходной сигнал: и т.д.

Заметим: коэффициент обратной связи больше и выходной сигнал идет "в разнос".

Рис. 1.5. Неустойчивый рекурсивный фильтр.

Операции, относящиеся к рекурсивной фильтрации, также известны в обычной практике, например - интегрирование. При интегрировании по формуле трапеций:

т.е. здесь мы имеем РЦФ с коэффициентами: .

Пример:

Уравнение РЦФ: начальные условия - нулевые.

Входной сигнал:

Выполните фильтрацию. Контроль: и т.д.

Рис. 1.6. Интегрирующий рекурсивный фильтр.

Импульсная реакция фильтров

Функция отклика

Если на вход нерекурсивного фильтра подать единичный импульс (импульс Кронекера), расположенный в точке , то на выходе фильтра мы получим его реакцию на единичный входной сигнал (формула ), которая определяется весовыми коэффициентами оператора фильтра:

Для рекурсивных фильтров реакция на импульс Кронекера зависит как от коэффициентов фильтра, так и от коэффициентов обратной связи . С использованием формулы :

Функция , которая связывает вход и выход фильтра по реакции на единичный входной сигнал и однозначно определяется оператором преобразования фильтра, получила название импульсного отклика фильтра (функции отклика). Если произвольный сигнал на входе фильтра представить в виде линейной комбинации взвешенных импульсов Кронекера

то, с использованием функции отклика, сигнал на выходе фильтра можно рассматривать как суперпозицию запаздывающих импульсных реакций на входную последовательность взвешенных импульсов:

Для нерекурсивных фильтров пределы суммирования в последнем выражении устанавливаются непосредственно по длине импульсного отклика . Для рекурсивных фильтров длина импульсного отклика, в принципе, может быть бесконечной.

Определение импульсной реакции

Определение импульсной реакции на практике требуется, как правило, только для рекурсивных фильтров, так как импульсная реакция для НЦФ при известных значениях коэффициентов , как это следует из выражения , специального определения не требует: .

Если выражение для системы известно в общей форме , определение импульсной реакции производится подстановкой в уравнение системы импульса Кронекера с координатой при нулевых начальных условиях. В соответствии с выражением сигнал на выходе системы будет представлять собой импульсную реакцию системы.

Пример:

Уравнение РЦФ: начальные условия - нулевые.

Входной сигнал:

Расчет выходного сигнала при нулевых начальных условиях:

и т.д.

Импульсный отклик фильтра:

Определение импульсной реакции физической системы обычно производится подачей на вход системы ступенчатой функции (функции Хевисайда), которая равна при , и при :

Отсюда:

Функция получила название переходной характеристики системы (перехода из одного статического состояния в другое). Форму реакции фильтра на функцию Хевисайда можно видеть на рис. 1.4 (с точки и далее) в сопоставлении с реакцией на импульс Кронекера в точке .

Передаточные функции фильтров

Z-преобразование

Удобным методом решения разностных уравнений линейных систем является z-преобразование. Применяя z-преобразование к обеим частям равенства , c учетом сдвига функций , получаем:

где - соответствующие z-образы входного и выходного сигнала.

Отсюда, полагая , получаем в общей форме функцию связи выхода фильтра с его входом - уравнение передаточной функции системы в z-области:

Для НЦФ, при нулевых коэффициентах :

При проектировании фильтров исходной, как правило, является частотная передаточная функция фильтра , по которой вычисляется ее Z-образ и обратным переходом в пространство сигналов определяется алгоритм обработки данных. В общей форме для выходных сигналов фильтра:

После обратного Z-преобразования выражения :

При подаче на вход фильтра единичного импульса Кронекера , имеющего z-образ , сигнал на выходе фильтра будет представлять собой импульсную реакцию фильтра , при этом:

т.е. передаточная функция фильтра является z-образом его импульсной реакции. При обратном z-преобразовании передаточной функции соответственно получаем импульсную характеристику фильтра:

Если функция представлена конечным степенным полиномом, что, как правило, характерно для НЦФ, являющихся КИХ-фильтрами, то обратное z-преобразование осуществляется элементарно идентификацией коэффициентов по степеням z. Передаточная функция РЦФ также может быть представлена степенным полиномом прямым делением числителя на знаменатель правой части выражения , однако результат при этом может оказаться как конечным, так и бесконечным, т.е. система может иметь либо конечную, либо бесконечную импульсную характеристику. Практически используемые рекурсивные фильтры обычно имеют бесконечную импульсную характеристику (БИХ-фильтры) при конечном числе членов алгоритма фильтрации .

Пример:

1. Передаточная функция РЦФ: . Прямым делением числителя на знаменатель получаем: .. Фильтр РЦФ является КИХ-фильтром.

2. Передаточная функция: . Методом обратного z-преобразования: . Фильтр РЦФ является БИХ-фильтром.

Устойчивость фильтров

Фильтр называется устойчивым, если при любых начальных условиях реакция фильтра на любое ограниченное воздействие также ограничена. Критерием устойчивости фильтра является абсолютная сходимость отсчетов его импульсного отклика:

Анализ устойчивости может быть проведен по передаточной функции. В устойчивой системе значение должно быть конечным во всех точках z-плоскости, где , а, следовательно, передаточная функция не должна иметь особых точек (полюсов) на и внутри единичного круга на z-плоскости. Полюсы определяются корнями многочлена знаменателя передаточной функции .

Пример:

Передаточная функция фильтра рис. 1.4: . При полюс знаменателя: Фильтр устойчив.

Передаточная функция фильтра рис. 1.5: . При полюс знаменателя: . Фильтр неустойчив, что и подтверждает пример фильтрации.

Передаточная функция фильтра рис. 1.6: . Полюс знаменателя: . В принципе, фильтр неустойчив, но эта неустойчивость проявляется только при . Импульсный отклик фильтра , сумма которого равна только при , т.е. при интегрировании бесконечно больших массивов. При интегрировании конечных массивов результат всегда конечен.

Приведенный критерий устойчивости относится к несократимой дроби, т.к. в противном случае возможна компенсация полюса нулем передаточной функции, и следует проверить наличие однозначных нулей и полюсов. Проверка на устойчивость требуется только для рекурсивных цифровых фильтров(систем с обратной связью), не рекурсивные системы всегда устойчивы.

Частотные характеристики фильтров

Общие понятия

От z-образов сигналов и передаточных функций подстановкой в уравнение можно перейти к Фурье-образам функций, т.е. к частотным спектрам сигналов и частотной характеристике фильтров, а точнее – к функциям их спектральных плотностей. Можно применить и способ получения частотных характеристик непосредственно из разностного уравнения системы обработки данных. Так как цифровая фильтрация относится к числу линейных операций, то, принимая для сигнала на входе фильтра выражение , мы вправе ожидать на выходе фильтра сигнал . Подставляя эти выражения в разностное уравнение фильтра , получаем:

Передаточная частотная функция (частотная характеристика при ):

Нетрудно убедиться, что полученная частотная характеристика повторяет функцию при , что и следовало ожидать. Аналогично z-преобразованию , частотная характеристика фильтра представляет собой Фурье-образ его импульсной реакции, и наоборот. При :

В общем случае является комплексной функцией, модуль которой называется амплитудно-частотной характеристикой (АЧХ), а аргумент фазово-частотной характеристикой (ФЧХ).

Выбор знака фазового угла ориентирован на каузальные системы с отрицательным временным запаздыванием сигналов. Допустим, что система осуществляет только сдвиг сигнала вправо по временной оси, т е. . Для преобразования Фурье функции имеем

Отсюда:

Из последнего равенства следует, что фаза представляет собой прямую с отрицательным тангенсом угла наклона Соответственно, для всех каузальных фильтров, осуществляющих преобразование с определенной задержкой сигнала на выходе, при выполнении операции над частотными составляющими сигнала имеет место:

C учетом отрицательного знака фазовой характеристики каузальных фильтров это вызывает сдвиг в "минус" всех частотных составляющих сигнала и соответствующую задержку выходного сигнала относительно входного.

На рис. 2.1 - 2.3 приведены частотные характеристики фильтров (модули и аргументы спектральных плотностей), которые были рассмотрены выше в примерах и на рис. 1.4 – 1.6. Графики приведены в границах главных диапазонов спектров, и получены непосредственной подстановкой при в уравнения передаточных функций .

Рис. 2.1. Спектр не имеет особых точек.
Рис. 2.2. Спектр имеет особые точки на границах диапазонов.
Рис. 2.3. Спектр интегрирующего фильтра. Особая точка на нулевой частоте. При обработке ограниченных массивов амплитуда центрального пика равна количеству точек массива.

Основные свойства

Основные свойства частотных характеристик цифровых фильтров:

  1. Частотные характеристики являются непрерывными функциями частоты.
  2. При дискретизации данных по интервалам функция является периодической. Период функции равен частоте дискретизации входных данных . Первый низкочастотный период (по аргументу от до , по от до ) называется главным частотным диапазоном. Граничные частоты главного частотного диапазона соответствуют частоте Найквиста , . Частота Найквиста определяет предельную частоту данных, которую способен обрабатывать фильтр.
  3. Для фильтров с вещественными коэффициентами импульсной реакции функция АЧХ является четной, а функция ФЧХ - нечетной. С учетом этого частотные характеристики фильтров обычно задаются только на интервале положительных частот главного частотного диапазона. Значения функций на интервале отрицательных частот являются комплексно сопряженными со значениями на интервале положительных частот.

Как правило, при частотном анализе фильтров значение интервала дискретизации принимают за 1, что соответственно определяет задание частотных характеристик на интервале по частоте или по . При использовании быстрых преобразований Фурье (БПФ) вычисления спектров осуществляются в одностороннем варианте положительных частот в частотном интервале от 0 до (от 0 до 1 Гц), где комплексно сопряженная часть спектра главного диапазона (от до ) занимает интервал от до (для ускорения вычислений используется принцип периодичности дискретных спектров). Заметим, что при выполнении БПФ количество точек спектра равно количеству точек входной функции, а, следовательно, отсчет на частоте , комплексно сопряженный с отсчетом на частоте , отсутствует. При нумерации точек входной функции от до он принадлежит точке - начальной точке следующего периода, при этом шаг по частоте равен . Современное программное обеспечение БПФ допускает любое количество точек входной функции, при этом для нечетного значения N частоте p соответствует отсчет на точке , не имеющий сопряженного отсчета, а при четном значении N отсутствует отчет на частоте p (она располагается между отсчетами и ). Отсчетам с номерами k главного диапазона БПФ (за исключением точки k=0) соответствуют комплексно сопряженные отсчеты (за исключением точки при нечетном ).

Фазовая и групповая задержка

Задержка сигналов во времени относится к характерной особенности каузальных систем в целом, а, следовательно, рекурсивных и односторонних нерекурсивных фильтров. Фазовая задержка, это прямая характеристика временной задержки фильтром гармонических колебаний. При подаче на вход фильтра гармоники , сигнал на выходе каузального фильтра, без учета изменения его амплитуды, равен , при этом:

Отсюда, фазовая задержка на частоте равна:

При распространении в целом на спектральную передаточную функцию фильтра получаем:

Постоянство значения в определенном частотном диапазоне обеспечивает для всех гармоник сигнала такое же соотношение их фазовых характеристик, какое было на входе системы, т.е. не изменяет формы сигнала, если его спектр полностью сосредоточен в этом частотном диапазоне, и значения АЧХ в этом диапазоне также имеют постоянное значение. Это условие является определяющим, например, для систем передачи данных, для сглаживающих и полосовых частотных фильтров.

Что касается каузальных фильтров, то они, как правило, имеют в рабочем диапазоне определенную зависимость значения от частоты, которая характеризуется групповым временем задержки (ГВЗ). ГВЗ характеризует среднюю временную задержку составного сигнала. Допустим, что сигнал на входе фильтра представляет собой сумму двух гармоник с близкими частотами:

Тождественная тригонометрическая запись:

Эта запись показывает, что сумму двух гармоник с частотами и можно рассматривать, как амплитудную модуляцию гармоники с частотой гармоникой с частотой . При прохождении через фильтр каждая из гармоник и может получить различную задержку, при этом сигнал на выходе фильтра, без учета амплитудных изменений:

Тождественная запись:

Пульсацию колебаний выразим через групповую временную задержку :

Отсюда:

При распространении этого выражения на непрерывную частотную характеристику фильтра:

Для вычислений ГВЗ удобно использовать комплексный логарифм передаточной функции фильтра:

Приближение для дискретных спектральных функций:

Различают фильтры с линейной и нелинейной фазовой характеристикой. В фильтрах с нелинейной фазовой характеристикой частотные компоненты сигнала задерживаются на величину, не пропорциональную частоте, и тем самым в выходном сигнале изменяется гармоническая связь между его компонентами, что может быть недопустимо во многих случаях обработки сигналов (передача данных, обработка биосигналов, воспроизведение музыки и видео, и пр.).

Чтобы фильтр имел линейную фазовую характеристику необходимо и достаточно, если выполняется одно из следующих условий:

где и - константы.

Условие обеспечивает постоянные значения групповой и фазовой задержки. Оно выполняется, если импульсная характеристика фильтра имеет положительную симметрию:

где нечётное;
чётное.

При этом фазовая характеристика является функцией длины фильтра:

Пример:

Импульсный отклик фильтра задан параметрами:. Передаточная функция фильтра: . Подставляем при и получаем частотную характеристику фильтра в главном диапазоне :

Изменяя обозначения и переходя к индексации относительно центра симметрии записываем в компактной форме:

Частотная характеристика фильтра линейна.

Условие обеспечивает постоянную групповую задержку и выполняется при отрицательной симметрии импульсной характеристики фильтра:

Для того чтобы убедиться в последнем, достаточно рассмотреть пример, аналогичный вышеприведенному.

Корреляция входа и выхода фильтров

Корреляция входа и выхода фильтров может быть получена на основе следующих простых соображений. Примем для входного сигнала и выходного сигнала за основу выражение преобразования в частотной области:

Умножим обе части этого выражения на комплексно сопряженную функцию и найдем математические ожидания левой и правой части:

Но математические ожидания этих произведений спектров представляют собой спектры плотности мощности, и, при обратном преобразовании Фурье, зависимость взаимной корреляционной функции входного и выходного сигналов фильтра от корреляционной функции входного сигнала и функции импульсного отклика фильтра:

Это выражение в спектральной области может использоваться для практического определения частотных передаточных функций фильтров с неизвестной формой импульсных откликов.

Если математические ожидания взять от квадратов модулей левой и правой части исходного выражения , то в результате получим выражения:

Области применения НЦФ и РЦФ

Области применения НЦФ и РЦФ обычно обуславливаются видом их передаточных функций. В принципе, нерекурсивные цифровые фильтры универсальны и способны реализовать любые практические задачи обработки сигналов. Это и понятно, т.к. реакция РЦФ на единичный импульс Кронекера представляет собой импульсный отклик НЦФ, а, следовательно, задачи, решаемые РЦФ, могут выполняться и НЦФ, но при условии отсутствия ограничений по размерам окна. В первую очередь это касается реализации БИХ-фильтров с незатухающим или слабо затухающим импульсным откликом, например, интегрирующих или фильтров рекурсивной деконволюции. Ограничение по размерам окна является скорее не теоретическим (бесконечных операторов НЦФ не требуется, максимум – двойная длина входного сигнала для двусторонних НЦФ), а чисто практическим. Нет смысла применять НЦФ с огромными размерами операторов и тратить машинное время, если та же задача во много раз быстрее решается рекурсивным фильтром.

Существенным преимуществом НЦФ является их устойчивость, возможность выполнения в виде двусторонних симметричных фильтров, не изменяющих фазу выходных сигналов относительно входных, и реализации строго линейных фазовых характеристик.

С другой стороны, нерекурсивные фильтры могут быть преобразованы в рекурсивные фильтры, если есть возможность z-полином передаточной функции НЦФ выразить в виде отношения двух коротких z-полиномов РЦФ типа , что может дать существенное повышение производительности вычислений. Как правило, такая возможность имеется для сходящихся степенных рядов. Отношение двух z-полиномов позволяет реализовать короткие и очень эффективные фильтры с крутыми срезами на частотных характеристиках.

Структурные схемы цифровых фильтров

Структурные схемы

Алгоритмы цифровой фильтрации сигналов (цифровых фильтров) представляются в виде структурных схем, базовые элементы которых показаны на рисунке 5.1 вместе с примерами структурных схем фильтров. Как правило, структурные схемы соответствуют программной реализации фильтров на ЭВМ, но не определяют аппаратной реализации в специальных радиотехнических устройствах, которая может существенно отличаться от программной реализации.

Рис. 5.1. Структурные схемы цифровых фильтров
Рис. 5.2. Граф фильтра

Наряду со структурной схемой фильтр может быть представлен в виде графа, который отображает диаграмму прохождения сигналов, и состоит из направленных ветвей и узлов. Пример структурной схемы фильтра с передаточной функцией и графа, ей соответствующего, приведен на рисунке 5.2. С каждым - узлом графа связано значение сигнала или его образа , которые определяются суммой всех сигналов или z-образов входящих в узел ветвей. В каждой - ветви (из узла в узел ) происходит преобразование сигнала в соответствии с передаточной функцией ветви, например задержка сигнала или умножение на коэффициент.

Соединения фильтров

    Рис. 5.3. Последовательное соединение.
  1. Последовательное соединение (рис. 5.3). Выходной сигнал предшествующего фильтра является входным для последующего. Эквивалентная передаточная функция общей системы равна произведению передаточных функций фильтров, в нее входящих:



    Рис. 5.4. Параллельное соединение.
  2. Параллельное соединение (рис. 5.4). Сигнал подается на входы всех параллельно соединенных фильтров одновременно, выходные сигналы фильтров суммируются. Эквивалентная передаточная функция общей системы равна сумме передаточных функций фильтров, в нее входящих:



    Рис. 5.5. Соединение обратной связи.
  3. Соединение обратной связи (рис. 5.5). Выходной сигнал первого фильтра подается на выход системы и одновременно на вход фильтра обратной связи, выходной сигнал которого суммируется, со знаком плюс или минус в зависимости от вида связи (отрицательной или положительной), с входным сигналом системы. Эквивалентная передаточная функция системы:




Схемы реализации фильтров

По принципам структурной реализации фильтров различают следующие схемы:

    Рис. 5.6. Прямая форма.
  1. Прямая форма (рис. 5.6) реализуется непосредственно по разностному уравнению:

    или по передаточной функции

    Рис. 5.7. Прямая каноническая форма .


  2. Прямая каноническая форма содержит минимальное число элементов задержки. Передаточную функцию РЦФ можно представить в следующем виде:




    Отсюда:



    В разностных уравнениях осуществляется только задержка сигналов . Граф реализации РЦФ в прямой канонической форме приведен на рисунке 5.7.

  3. Каскадная (последовательная) форма соответствует представлению передаточной функции в виде произведения:

    - составляющие функции вида при представлении в факторизованной форме, где и - нули и полюсы функции . В качестве функций обычно используются передаточные функции биквадратных блоков - фильтров второго порядка:


  4. Параллельная форма используется много реже, и соответствует представлению передаточной функции в виде суммы биквадратных блоков или более простых функций.

Обращенные формы

В теории линейных направленных сигнальных графов существуют процедуры преобразования исходных графов с сохранением передаточных функций. Одна из таких процедур - обращение (транспозиция) графов, которая выполняется путем изменения направления всех ветвей цепи, при этом вход и выход графа также меняются местами. Для ряда систем такая транспозиция позволяет реализовать более эффективные алгоритмы обработки данных. Пример обращения графа прямой канонической формы рекурсивной системы (с перестроением графа на привычное расположение входа с левой стороны) приведен на рис. 2.5.8.

Литература

  1. Баскаков С.И. Радиотехнические цепи и сигналы: Учебник для вузов. - М.: Высшая школа, 1988.- 448 с.
  2. Гольденберг Л.М. и др. Цифровая обработка сигналов: Справочник. - М.: Радио и связь, 1985.- 312 с.
  3. Гольденберг Л.М. и др. Цифровая обработка сигналов: Учебное пособие для вузов. - М.: Радио и связь, 1990.- 256 с.
  4. Клаербоут Д.Ф. Теоретические основы обработки геофизической информации с приложением к разведке нефти. – М.: Недра, 1981. – 304 с.
  5. Рабинер Л., Гоулд Б. Теория и применение цифровой обработки сигналов. – М.: Мир, 1978. – 848 с.
  6. Хемминг Р.В. Цифровые фильтры. – М.: Недра, 1987. – 221 с.
  7. Айфичер Э., Джервис Б. Цифровая обработка сигналов. Практический подход. / М., "Вильямс", 2004, 992

См. также