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

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

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

Содержание

Введение

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

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

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

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

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

Литература

Введение

Первый натиск и первые крики решают дело.

Тит Ливий. Римский историк, 59 г.д.н.э.-17 г.н.э.

Роль крика в драке - существенный вклад в теорию конфликтов. Но имеет ли он такое же значение при фильтрации данных - не очевидно, Лично я предпочитаю воплем завершать этот процесс, а не начинать.

Эдуард Павелко. Новосибирский геофизик Уральской школы, ХХ в.

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

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

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

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

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

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

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

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

y(kDt) = TL{x(kDt)}.

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

a_m y(kDt-mDt) =b_n x(kDt-nDt), (1.1)

где k = 0, 1, 2, …- порядковый номер отсчетов, Dt - интервал дискретизации сигнала, a_m и b_n - вещественные или комплексные коэффициенты. Положим a₀ = 1, что всегда может быть выполнено соответствующей нормировкой уравнения (1.1), и, принимая в дальнейшем Dt = 1, т.е. переходя к числовой нумерации цифровых последовательностей значений сигналов, приведем его к виду:

y(k) = b_n x(k-n) -a_m y(k-m). (1.2)

При k < n и m проведение фильтрации возможно только при задании начальных условий для точек x(-k), k = 1, 2, … , N, и y(-k), k = 1, 2, … , M. Как правило, в качестве начальных условий принимаются либо нулевые значения, либо выполняется продление отсчетов входных сигналов или его тренда по отрицательным значениям аргумента.

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

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

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

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

· Один фильтр может обрабатывать несколько входных каналов или сигналов.

· Входные и выходные данные можно сохранять для последующего использования.

· Точность цифровых фильтров ограничена только разрядностью отсчетов.

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

Нерекурсивные фильтры. При нулевых значениях коэффициентов a_m уравнение (1.2) переходит в уравнение линейной дискретной свертки функции x(k) с оператором b_n:

y(k) = b_n x(k-n). (1.3)

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

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

y(k) =b_n x(k-n). (1.4)

При N' = N фильтр называется двусторонним симметричным. Симметричные фильтры, в отличие от односторонних фильтров, не изменяют фазы обрабатываемого сигнала.

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

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

Представим, что на одной полоске бумаги выписаны по порядку сверху вниз значения данных x(k) ? s_k (см. рис. 1.1). На второй полоске бумаги находятся записанные в обратном порядке значения коэффициентов фильтра b_n ? h_n (обозначение h для коэффициентов операторов НЦФ является общепринятым). Для вычисления y_k ? y(k) располагаем вторую полоску против первой таким образом, чтобы значение h₀ совпало со значением s_k, перемножаем все значения h_n с расположенными против них значениями s_k-n, и суммируем все результаты перемножения. Результат суммирования является выходным значением сигнала y_k. Сдвигаем окно фильтра - полоску коэффициентов h_k, на один отсчет последовательности s_k вниз (или массив s_k сдвигаем на отсчет вверх) и вычисляем аналогично следующее значение выходного сигнала, и т.д.

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

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

y(k) = b(n) _* x(k-n) b(n) b x(k-n).

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

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

y_k = 0.2(x_k-2+x_k-1+x_k+x_k+1+x_k+2).

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

y_k =b_n x_k-n, b_n = 0,2. (1.5)

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

y_k = (-3x_k-2+12x_k-1+17x_k+12x_k+1-3x_k+2)/35. (1.6)

Это также НЦФ с коэффициентами: b₀ = 17/35, b₁ = b_-1 = 12/35, b₂ = b_-2 = -3/35.

Пример. Уравнение НЦФ: yk =bn xk-n, bn = 0,2. Начальные условия - нулевые.

Входной сигнал - скачок функции (ступень): xk = {0,0,0,0,0,0,10,10,10,10,…}.

Выходной сигнал: yk = {0,0,0,0,2,4, 6, 8,10,10,10,10,…}.

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

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

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

Повторите фильтрацию фильтром МНК на основе кубического уравнения. Сравните результаты фильтрации с результатами первого НЦФ (приведены на рис. 1.2(В)).

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

· Дистрибутивность: h(n) b [a(k)+b(k)] = h(n) b a(k)+h(n) b--b(k).

· Коммутативность: h(n) b a(k) b b(k) = a(k) b--b(k) b h(n).

· Ассоциативность: [a(k) b b(k)] b h(n) = h(n) b a(k) b b(k).

Фильтрация однозначно определяет выходной сигнал y(k) для установленного значения входного сигнала s(k) при известном значении импульсного отклика фильтра h(n).

Рекурсивные фильтры. Фильтры, которые описываются полным разностным уравнением (1.2)

y(k) = b_n x(k-n) -a_m y(k-m),

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

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

Пример. Уравнение РЦФ: yk = boxk+a1yk-1, при bo = a1 = 0.5, y-1 = 0.

Входной сигнал: xk = {0,0,1,0,0,0,0,0,0,0,1,1,1,1,1....}

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

уо = 0,5xo + 0,5y-1 = 0; y1 = 0,5x1 + 0,5yo =0; y2 = 0,5x2 + 0,5y1 = 0.5; y3 = 0,5x3 + 0,5y2 = 0.25;

y4 = 0,5x4 + 0,5y3 = 0.125; y5 = 0,5x5 + 0,5y4 = 0.0625; y6 = 0,5x6 + 0,5y5 = 0.03125; и т.д.

Выходной сигнал: yk = {0, 0, 0.5, 0.25, 0.125, 0.0625, 0.03125, 0.015625,...}

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

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

При положительной обратной связи (сумма коэффициентов a_m больше 1) фильтр становится неустойчивым (идет «в разнос» как на рис. 1.5)

Пример. Уравнение РЦФ: yk = boxk - a1yk-1, при bo =0.5, a1=1.1, y-1= 0

Входной сигнал: xk = {0, 10, 0, 0, 0,....}.

Выходной сигнал: yk = {0,0,5,-5.5,6.05,-6.655,7.321,-8.053,8.858,-9.744,10.718,-11.79,… и т.д.}

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

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

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

yk = (xk+xk-1)/2 + yk-1, (2.1.7)

т.е. здесь мы имеем РЦФ с коэффициентами: bo = b1 = 0.5, a1 = 1.

Пример. Уравнение РЦФ: yk=(xk+xk-1)/2+yk-1, начальные условия - нулевые.

Входной сигнал: xk={0,0,2,2,4,0,0,0,4,4,4,0,0,0,5,0,0,0,....}

Выполните фильтрацию. Контроль: yk= {0,0,0,1,3,6,8,8,8,10,14,18,20,20,20,22.5,25,25,25...}

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

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

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

y(k) = TL[d(0)] = b_n b--d(k-n) = h(k) ? b_n. (2.1)

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

y(k) = b_n d(k-n) -a_m y(k-m) = h_k. (2.1')

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

Если произвольный сигнал на входе фильтра представить в виде линейной комбинации взвешенных импульсов Кронекера

x(k) =d(n) x(k-n),

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

y(k) = h(n) (d(n) x(k-n)) h(n) x(k-n).

Для нерекурсивных фильтров пределы суммирования в последнем выражении устанавливаются непосредственно по длине импульсного отклика h(n).

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

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

Пример. Уравнение РЦФ: yk = xk + 0.5yk-1.

Входной сигнал: xk= do= {1,0,0,0,...}.

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

yo = xo+0.5 y-1 = 1+0 = 1 = ho. y1 = x1+0.5 yo = 0+0.5 = 0.5 = h1. y2 = x2+0.5 y1 = 0+0.25 = 0.25 = h2.

y3 = x3+0.5 y2 = 0.125 = h3. y4 = x4+0.5 y3 = 0.0625 = h4, и т.д.

Импульсный отклик фильтра: hk = (0.5)k, k = 0, 1, 2....

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

g(k) =h(n) u(k-n) =h(n).

h(k) = g(k) - g(k-1).

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

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

Z-преобразование. Удобным методом решения разностных уравнений линейных систем является z-преобразование.

Напоминание.

yn - yn zn = Y(z).

Аргумент z используется в двух вариантах: zn или z-n. В данном курсе по умолчанию используем аргумент zn.

Применяя z-преобразование к обеим частям равенства (1.1), c учетом сдвига функций (y(k-m) z^m Y(z)), получаем:

Y(z)a_mz^m = X(z)b_nzⁿ, (3.1)

где X(z),Y(z)- соответствующие z-образы входного и выходного сигнала. Отсюда, полагая a_o = 1, получаем в общей форме уравнение передаточной функции системы в z-области:

H(z) = Y(z)/X(z) =b_nzⁿ(1+a_mz^m). (3.2)

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

H(z) =b_nzⁿ. (3.3)

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

Y(z) = H(z)·X(z).

Y(z)·(1+a_m z^m) = X(z)b_n zⁿ

Y(z) = X(z)b_n zⁿ - Y(z)a_m z^m. (3.4)

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

y(k) =b_n x(k-n) -a_m y(k-m). (3.5)

При подаче на вход фильтра импульса Кронекера d_о, имеющего z-образ d(z) = zⁿ = 1, сигнал на выходе фильтра будет представлять собой импульсную реакцию фильтра y(k) ? h(k), при этом:

H(z) = Y(z)/d(z) = Y(z) = TZ[y(k)] =h(k) z^k, (3.6)

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

h(k) H(z). (3.7)

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

Примеры.

1. Передаточная функция РЦФ: H(z) = (1-z5)/(1-z).

Прямым делением числителя на знаменатель получаем: H(z) = 1+z+z2+z3+z4.

H(z) h(n) = {1,1,1,1,1}. Фильтр РЦФ является КИХ-фильтром.

2. Передаточная функция: H(z) = 1/(1-2z).

Методом обратного z-преобразования: h(n) = 2n. Фильтр РЦФ является БИХ-фильтром.

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

|h(n)| < . (3.8)

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

Пример.

Передаточная функция фильтра рис. 1.4: H(z) = b0/(1-a1z). При а1= 0.5 полюс знаменателя: zр= 2. |zр|>1. Фильтр устойчив.

Передаточная функция фильтра рис. 1.5: H(z) = b0/(1+a1z). При а1= 1.1 полюс знаменателя: zр= -0.909. |zр| < 1. Фильтр неустойчив, что и подтверждает пример фильтрации.

Передаточная функция фильтра рис. 1.6: H(z) = 0.5(1+z)/(1-z). Полюс знаменателя: zр= 1. В принципе, фильтр неустойчив, но эта неустойчивость проявляется только при k = ?. Импульсный отклик фильтра h(n) = {0.5,1,1,1, ….}, сумма которого равна ? только при n = ?, т.е. при интегрировании бесконечно больших массивов. При интегрировании конечных массивов результат всегда конечен.

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

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

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

Общие понятия. От z-образов передаточных функций подстановкой z = exp(-jwDt) в уравнение (3.2) можно перейти к Фурье-образам функций, т.е. к частотным спектрам сигналов и частотной характеристике фильтров, а точнее - к функциям их спектральных плотностей.

Можно применить и способ получения частотных характеристик непосредственно из разностного уравнения системы обработки данных. Так как цифровая фильтрация относится к числу линейных, то, принимая для сигнала на входе фильтра выражение x(kDt) = B(w) exp(jwkDt), мы вправе ожидать на выходе фильтра сигнал y(kDt) = A(w) exp(jwkDt). Подставляя эти выражения в разностное уравнение фильтра (1.1), получаем:

a_m A(w) exp(jwkDt-jwmDt) =b_n B(w) exp(jwkDt-jwnDt).

A(w) exp(jwkDt) a_m exp(-jwmDt) = B(w) exp(jwkDt)b_n exp(-jwnDt).

A(w)a_m exp(-jwmDt) = B(w)b_n exp(-jwnDt). (4.1)

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

H(w) = A(w)/B(w) =b_n exp(-jwnDt)[1+a_m exp(-jwmDt)]. (4.2)

Полученная частотная характеристика повторяет функцию (3.2) при z = exp(-jwDt), что и следовало ожидать. Аналогично z-преобразованию (3.7), частотная характеристика фильтра представляет собой Фурье-образ его импульсной реакции, и наоборот. При Dt = 1:

H(w) =h(n) exp(-jwn), (4.3)

h(n) = (1/2w)H(w) exp(jwn) dw. (4.4)

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

A(w) = |H(w)| =

j(w) = arctg(-Im H(w)/Re H(w)).

Выбор знака фазового угла ориентирован на каузальные системы с отрицательным временным запаздыванием сигналов. Допустим, что система осуществляет только сдвиг сигнала x(t) вправо по временной оси, т е. y(t) = x(t-t). Для преобразования Фурье функции y(t) имеем:

Y(f) =y(t) exp(-j2pft) dt =x(t-t) exp(-j2pft) dt =

= exp(-j2pft)x(t) exp(-j2pft) dt = exp(-j2pft) X(f).

Отсюда:

H(f) = Y(f)/X(f) = exp(-j2pft), |H(f)| = 1, j_h(f) = -2pft.

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

Y(f)= H(f) X(f)= |H(f)|exp(jj_h(f)) |X(f)|exp(jj_x(f))= |H(f)| |X(f)| exp{j--[j_h(f)+j_x(f)]},

|Y(f)| = |H(f)| |X(f)|, j_y(f) = j_h(f)+j_x(f).

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

На рис. 4.1-4.3 приведены частотные характеристики фильтров (модули и аргументы спектральных плотностей), которые были рассмотрены выше в примерах и на рис. 1.4 - 1.6. Графики приведены в границах главных диапазонов спектров, и получены непосредственной подстановкой z=exp(-jwDt) при Dt=1 в уравнения передаточных функций H(z).

Рис. 4.1. Спектр не имеет особых точек.



Рис. 4.2. Спектр имеет особые точки на границах диапазонов.

Рис. 4.3. Спектр интегрирующего фильтра. Особая точка на нулевой частоте.

При обработке ограниченных массивов амплитуда центрального пика равна количеству точек массива.

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

1. Частотные характеристики являются непрерывными функциями частоты.

2. При дискретизации данных по интервалам Dt функция H(w) является периодической. Период функции H(w) равен частоте дискретизации входных данных F = 1/Dt. Первый низкочастотный период (по аргументу w от -p/Dt до p/Dt, по f от -1/2Dt до 1/2Dt) называется главным частотным диапазоном. Граничные частоты главного частотного диапазона соответствуют частоте Найквиста w_N, w_N = p/Dt. Частота Найквиста определяет предельную частоту данных, которую способен обрабатывать фильтр.

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

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

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

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

sin(wt-j) = sin w(t-t_p), щt-j = щ(t-t_p).

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

t_p = j--/щ. (4.5')

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

t_p(w)= j(w)--/щ. (4.5)

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

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

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

s(t) = cos щ₁t + cos щ₂t.

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

s(t) = 2 cos 0.5(щ₁+щ₂)t · cos 0.5(щ₁-щ₂)t.

Эта запись показывает, что сумму двух гармоник с частотами щ₁ и щ₂ можно рассматривать, как амплитудную модуляцию гармоники с частотой (щ₁+щ₂)/2 гармоникой с частотой (щ₁-щ₂)/2.

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

s(t) = cos (щ₁t-j₁) + cos (щ₂t-j₂).

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

s(t) = 2 cos[0.5((щ₁+щ₂)t-(j₁+j₂))] · cos[0.5((щ₁-щ₂)t-(j₁-j₂))].

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

cos[0.5((щ₁-щ₂)t-(j₁-j₂))] = cos[0.5(щ₁-щ₂)·(t-t_g)].

Отсюда:

(щ₁-щ₂)·t_g = j₁-j₂.

t_g = (j₁-j₂)/(щ₁-щ₂) = Dj--/Dщ. (4.6)

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

t_g(щ)= dj(w)/dщ. (4.7)

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

Ln H(щ) = ln |H(щ)| + j·j(щ), j(щ) = Im [Ln H(щ)].

t_g(щ)= dj--/dщ = Im{d[Ln(H(щ))]/dщ} = Im{dH(щ)/[H(щ)dщ]}. (4.8)

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

t_g(k·щ) ? (2/щ) Im{(H_k+1-H_k) / (H_k+1+H_k)}. (4.9)

Различают фильтры с линейной и нелинейной фазовой характеристикой.

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

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

j(w)--=--aw,-- (4.10)

j(w)--=--b-----aw,-- (4.11)

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

h(n) = h(N-n-1), n = 0, 1, 2, …, (N-1)/2, N - нечетное;

n = 0, 1, 2, …, (N/2)-1, N - четное.

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

a--=--(N-1)/2.

Пример.

Импульсный отклик фильтра задан параметрами: N=7, h(0)=h(6), h(1)=h(5), h(2)=h(4), h(3).

Передаточная функция фильтра: H(z) = h(k) zk. Подставляем z=exp(-jwDt) при Dt = 1 и получаем частотную характеристику фильтра в главном диапазоне (-p, p):

H(w) = h(0)+h(1)exp(-jw)+h(2)exp(-2jw)+ h(3)exp(-3jw)+h(4)exp(-4jw)+

+h(5)exp(-5jw)+h(6)exp(-6jw) = exp(-3jw) {h(0)[exp(3jw)+exp(-3jw)] +

+h(1)[exp(2jw)+exp(-2jw)] + h(2)[exp(jw)+exp(-jw)] + h(3)} =

= exp(-3jw) {2h(0) cos 3jw + 2h(1) cos 2jw + 2h(2) cos jw + h(3)}.

Изменяя обозначения и переходя к индексации относительно центра симметрии a(0) = h(3), a(n) = 2h(3-n), n=1, 2, 3, записываем в компактной форме:

H(w) =a(n) cos(njw) exp(-3jw) = |H(w)|exp(jj(w)), j(w) = -3w, a=3(N-1)/2.

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

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

h(n) = -h(N-n-1),

a = (N-1)/2, b = p/2.

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

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

Примем для входного сигнала x(t)X(f) и выходного сигнала y(t)Y(f) за основу выражение преобразования в частотной области

Y(f) = H(f) X(f). (4.12)

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

M{X*(f) Y(f)} = M{X*(f) H(f) X(f)} = H(f) M{X*(f) X(f)}.

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

W_xy = H(f) W_x h(t) b B_x(t) = B_xy(t).

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

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

W_y(f) = |H(f)|² W_x(f) B_h(t) b B_x(t).

Области применения НЦФ и РЦФ обычно обуславливаются видом их передаточных функций.

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

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

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

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

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

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

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

Пример структурной схемы фильтра с передаточной функцией H(z) = (1+b₁z)/(1+a₁z) и графа, ей соответствующего, приведен на рисунке 5.2. С каждым i - узлом графа связано значение сигнала x_i(k) или его образа X_i(z), которые определяются суммой всех сигналов или z-образов входящих в узел ветвей. В каждой ij - ветви (из узла i в узел j) происходит преобразование сигнала в соответствии с передаточной функцией ветви, например задержка сигнала или умножение на коэффициент.

Рис. 5.2. Граф фильтра.

Соединения фильтров. Различают следующие соединения фильтров.

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

Рис. 5.3.

H(z) = H₁(z)H₂(z)...H_N(z).

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

Рис. 5.4.

H(z) = H₁(z)+H₂(z)+...+H_N(z).

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

H(z) = H₁(z)/(1H₁(z)H₂(z)).

Рис. 2.5.5.

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

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

y_k =b_nx_k-n -a_my_k-m,

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

H(z) =b_nzⁿ /(1+a_mz^m).

Рис. 5.6.

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

H(z) = Y(z)/X(z) = H₁(z)H₂(z),

H₁(z) = V(z)/X(z) = 1/(1+a_mz^m),

H₂(z) = Y(z)/V(z) =b_nzⁿ.

Отсюда: v(k) = x(k) -a_mv(k-m), (5.1)

y(k) =b_nv(k-n). (5.2)

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

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



Рис. 5.7.

H(z) =H_i(z).

H_i(z) - составляющие функции вида (1-r_iz)/(1-p_iz) при представлении H(z) в факторизованной форме, где r_i и p_i - нули и полюсы функции H(z). В качестве функций H_i(z) обычно используются передаточные функции биквадратных блоков - фильтров второго порядка:

H_i(z) = (b_0i + b_1i z + b_2i z²) / (1 + a_1i z + a_2i z²).

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

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



Рис. 5.8.

рекурсивный цифровой фильтр сигнал

Литература

1. Баскаков С.И. Радиотехнические цепи и сигналы Учебник для вузов. - М. Высшая школа, 1988.- 448 с.

2. Гольденберг Л.М. и др. Цифровая обработка сигналов: Справочник. - М.: Радио и связь, 1985.- 312 с.

3. Гольденберг Л.М. и др. Цифровая обработка сигналов: Учебное пособие для вузов. - М.: Радио и связь, 1990.- 256 с.

4. Клаербоут Д.Ф. Теоретические основы обработки геофизической информации с приложением к разведке нефти. - М.: Недра, 1981. - 304 с.

5. Рабинер Л., Гоулд Б. Теория и применение цифровой обработки сигналов. - М.: Мир, 1978. - 848 с.

6. Хемминг Р.В. Цифровые фильтры. - М.: Недра, 1987. - 221 с.

7. Айфичер Э., Джервис Б. Цифровая обработка сигналов. Практический подход. / М., "Вильямс", 2004, 992 с.

Размещено на Allbest.ru

...