Свойства линейных дискретных фильтров

Содержание

Слайд 2

Там же было установлено, что устойчивость ЛДФ определяется условием, накладываемым на

Там же было установлено, что устойчивость ЛДФ определяется условием, накладываемым

на его импульсную характеристику. Была доказана следующая теорема.

Теорема 1. Необходимым и достаточным условием устойчивости ЛДФ, является выполнение условия.

(1)

Здесь C - некоторая константа.

Теперь мы посмотрим на проблему устойчивости линейного фильтра с другой точки зрения. Для выяснения устойчивости ЛДФ мы обратимся к передаточной функции H(z).

Слайд 3

Теорема 2. Для того чтобы ЛДФ был устойчив, необходимо и достаточно,

Теорема 2. Для того чтобы ЛДФ был устойчив, необходимо и достаточно,

чтобы все полюсы zn передаточной функции H(z) лежали в комплексной плоскости внутри единичного круга |z | < 1.
Теорему 2 примем без доказательства.
На следующих двух рисунках изображены полюсы передаточных функций двух ЛДФ.
Слайд 4

На первом рисунке показаны полюсы передаточной функции H1(z) первого фильтра. Полюсов

На первом рисунке показаны полюсы передаточной функции H1(z) первого фильтра.

Полюсов три и все они находятся внутри единичного круга | z | < 1. Поэтому первый фильтр устойчив.
Слайд 5

На втором рисунке показаны полюсы передаточной функции H2(z) второго фильтра. Полюсов

На втором рисунке показаны полюсы передаточной функции H2(z) второго фильтра.

Полюсов тоже три, два из них z1, z2 лежат внутри единичного круга | z | < 1, а один полюс z3 вне единичного круга | z | ≥ 1. Значит, второй фильтр неустойчив.
Слайд 6

Пример 1. На прошлой лекции мы рассматривали рекурсивный фильтр, который описывался

Пример 1. На прошлой лекции мы рассматривали рекурсивный фильтр, который описывался

следующим разностным уравнением.

(2)

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

(3)

Слайд 7

Единственная особая точка передаточной функции – это полюс первого порядка .

Единственная особая точка передаточной функции – это полюс первого порядка

. Согласно теореме 2, фильтр будет устойчив, если особая точка будет лежать внутри единичного круга комплексной плоскости. Поэтому фильтр будет устойчив, если выполняется неравенство.

(4)

Слайд 8

Частотная характеристика ЛДФ Определение. Частотная характеристика линейного дискретного фильтра – это

Частотная характеристика ЛДФ

Определение. Частотная характеристика линейного дискретного фильтра – это

комплексная функция K(ω) действительного переменного ω, и определяется соотношением.

(5)

Здесь ω - это циклическая частота ω = 2πf, а Δt – шаг дискретизации. Таким образом, чтобы найти частотную характеристику фильтра, надо в переходной функции сделать замену.

(6)

Слайд 9

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

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

характеристики h(n) фильтра.

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

(7)

Заменяя в сумме (7) комплексную переменную z с помощью замены (6), получаем формулу для нахождения частотной характеристики фильтра.

Слайд 10

Поэтому частотную характеристику можно представить в виде ряда. (8) Если частотную

Поэтому частотную характеристику можно представить в виде ряда.

(8)

Если

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

(9)

Слайд 11

Вспомним, как определяется спектр дискретного сигнала. Для аналогового сигнала s(t) выбирается

Вспомним, как определяется спектр дискретного сигнала. Для аналогового сигнала s(t)

выбирается шаг дискретизации Δt , и дискретный сигнал sn находят по формуле

(10)

Затем спектр дискретного сигнала находят по формуле

(11)

Здесь F - частота Найквиста.

(12)

Слайд 12

Вспомним, что условие физической реализуемости ЛДФ накладывает на импульсную характеристику условие.

Вспомним, что условие физической реализуемости ЛДФ накладывает на импульсную характеристику

условие.

(13)

Используя это условие, нижний предел суммирования в формуле (9) для частотной характеристики можно формально отодвинуть до - ∞ (минус бесконечности). Кроме того, в этой сумме заменим шаг дискретизации через частоту Найквиста по формуле (12). В результате получим.

(14)

Слайд 13

Мы видим, что формула (14) для частотной характеристики ЛДФ, и формула

Мы видим, что формула (14) для частотной характеристики ЛДФ, и

формула (11) для спектра дискретного сигнала имеют одинаковую структуру. Поэтому, если импульсную характеристику h(n) рассматривать как дискретный сигнал, то частотная характеристика H( f ) будет равна спектру этого дискретного сигнала HD( f ) , умноженному на удвоенную частоту Найквиста.

(15)

Так как частотная характеристика является спектром дискретного сигнала, то все свойства спектра дискретного сигнала, применимы к частотной характеристики ЛДФ.

Слайд 14

Так, например, частотная характеристика является периодической функцией частоты с периодом 2F

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

2F .

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

(16)

Слайд 15

КИХ и БИХ фильтры Определение. Конечной импульсной характеристикой ЛДФ будем называть

КИХ и БИХ фильтры

Определение. Конечной импульсной характеристикой ЛДФ будем называть

импульсную характеристику, имеющую конечное число элементов отличных от нуля.

Другими словами, всегда существует такое число M , что для любого номера m большего M все элементы импульсной характеристики h(m) равны нулю.

Фильтры с конечной импульсной характеристикой будем называть КИХ - фильтрами (английский термин – finite impulse response, FIR).

Слайд 16

Бесконечной импульсной характеристикой ЛДФ будем называть импульсную характеристику, имеющую бесчисленное число

Бесконечной импульсной характеристикой ЛДФ будем называть импульсную характеристику, имеющую бесчисленное число

элементов отличных от нуля.

Другими словами, какое бы большое число M мы не взяли, обязательно найдется такой номер m , для которого элемент импульсной характеристики h(m) будет отличен от нуля.
Фильтры с бесконечной импульсной характеристикой будем называть БИХ - фильтрами (английский термин – infinite impulse response, IIR).

На рисунках показаны импульсные характеристики КИХ и БИХ фильтров.

Слайд 17

Слайд 18

Слайд 19

Рекурсивные и нерекурсивные фильтры и их связь с КИХ и БИХ

Рекурсивные и нерекурсивные фильтры и их связь с КИХ и БИХ

фильтрами

Напомним, что основное разностное уравнение ЛДФ, имеет вид.

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

(18)

(17)

Слайд 20

В зависимости от того, равны нулю все или не все коэффициенты

В зависимости от того, равны нулю все или не все

коэффициенты an , фильтры разделяют на рекурсивные и нерекурсивные.
Если все коэффициенты an равны нулю, то мы получаем уравнение нерекурсивного фильтра.

(19)

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

(20)

Слайд 21

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

Найдем импульсную характеристику h(n) нерекурсивного фильтра. Вспомним, что импульсную характеристику

называют так же реакцией системы на единичный импульс. Это означает, что если на вход ЛДФ подать единичный импульс в виде символа Кронекера, то на выходе получился сигнал, совпадающий с импульсной характеристикой h(n).
Возьмем уравнение (19) нерекурсивного фильтра и поставим вместо входящего сигнала символ Кронекера.

В результате получим импульсную характеристику.

(21)

(22)

Слайд 22

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

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

равна коэффициентам bk основного разностного уравнения фильтра.

(23)

Формула (23) показывает, что полученная импульсная характеристика имеет конечное число отличных от нуля элементов, поэтому она является конечной импульсной характеристикой.
Значит все нерекурсивные фильтры являются КИХ – фильтрами.

Слайд 23

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

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

фильтра (20) перепишем в следующем виде.

(24)

Из формулы (24) мы видим, что у переходной функции имеется одна особая точка z = 0. Эта особая точка для членов суммы (24) является полюсом порядка m , где m - номер члена в сумме (24). Например, для члена

это будет полюс второго порядка m = 2 .

Слайд 24

Для анализа устойчивости фильтра, важно, что эта особая точка лежит в

Для анализа устойчивости фильтра, важно, что эта особая точка лежит

в комплексной плоскости в центре единичного круга | z | < 1. Отсюда по теореме 2 следует, что фильтр устойчив.
Значит все КИХ – фильтры являются устойчивыми фильтрами.

Пример 2. Нерекурсивный фильтр описывается следующим разностным уравнением.

(25)

Слайд 25

Найти импульсную характеристику фильтра h(n), переходную функцию фильтра H(z) и частотную

Найти импульсную характеристику фильтра h(n), переходную функцию фильтра H(z) и

частотную характеристику фильтра K(ω).
По формуле (23) находим импульсную характеристику.

(26)

По формуле (24) находим переходную функцию.

(27)

Слайд 26

Далее делая в формуле (27) для переходной функции подстановку получаем частотную

Далее делая в формуле (27) для переходной функции подстановку

получаем частотную характеристику.

(28)


Если хотя бы один коэффициент отличен от нуля, то фильтр называется рекурсивным фильтром. Рекурсивный фильтр – это фильтр с обратными связями.

Слайд 27

Пример 3. Рекурсивный фильтр описывается следующим разностным уравнением. (29) Найти импульсную

Пример 3. Рекурсивный фильтр описывается следующим разностным уравнением.

(29)

Найти импульсную

характеристику фильтра h(n), переходную функцию фильтра H(z) и частотную характеристику фильтра K(ω). Этот фильтр характеризуется следующими коэффициентами.

На рисунке показана структурная схема этого фильтра.

(30)

Слайд 28

Слайд 29

Этот фильтр был нами разобран на прошлой лекции. Поэтому воспользуемся полученными

Этот фильтр был нами разобран на прошлой лекции. Поэтому воспользуемся

полученными результатами.
Используя коэффициенты (30), по общей формуле для переходной функции

(31)

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

Слайд 30

Далее делая в формуле (32) для переходной функции подстановку (33) получаем частотную характеристику. (32)

Далее делая в формуле (32) для переходной функции подстановку

(33)

получаем частотную

характеристику.

(32)

Слайд 31

Затем, используя явный вид переходной функции (32) , с помощью теоремы

Затем, используя явный вид переходной функции (32) , с помощью

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

(34)

Формула (34) показывает, что полученная импульсная характеристика имеет бесконечное число отличных от нуля элементов, поэтому она является бесконечной импульсной характеристикой.
Таким образом, рассматриваемый рекурсивный фильтр является БИХ – фильтром.

Слайд 32

Исследуем этот фильтр на устойчивость. Из формулы (32) видно, что переходная

Исследуем этот фильтр на устойчивость. Из формулы (32) видно, что

переходная функция H(z) имеет одну особую точку . Эта особая точка не попадает внутрь единичного круга | z | < 1 комплексной плоскости. Эта точка находится на краю этого круга как это видно из рисунка.

Отсюда по теореме 2 следует, что рассматриваемый БИХ - фильтр неустойчив.

Слайд 33

Пример 4. На рисунке показана структурная схема рекурсивного фильтра.

Пример 4. На рисунке показана структурная схема рекурсивного фильтра.

Слайд 34

Написать основное разностное уравнение этого фильтра. Найти импульсную характеристику фильтра h(n),

Написать основное разностное уравнение этого фильтра. Найти импульсную характеристику фильтра

h(n), переходную функцию фильтра H(z) и частотную характеристику фильтра K(ω). Из структурной схемы видно, что фильтр характеризуется следующими коэффициентами.

(35)

Подставляя коэффициенты (35) в общее разностное уравнение ЛДФ (17) получаем.

(36)

Слайд 35

По общей формуле (31) для переходной функции, подставляя туда коэффициенты (35)

По общей формуле (31) для переходной функции, подставляя туда коэффициенты

(35) находим переходную функцию рассматриваемого фильтра.

(38)

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

(37)

Слайд 36

Выпишем несколько значений функции f(z) для некоторых n . (39) Из

Выпишем несколько значений функции f(z) для некоторых n .

(39)

Из

соотношений (39) видно, что особая точка у функции f(z) будет только в одном случае, когда n = 0 . По общей теории вычетов это означает, что все элементы импульсной характеристики равны нулю кроме одного.

(40)

Слайд 37

Если n = 0 , то функция f(z) имеет полюс первого

Если n = 0 , то функция f(z) имеет полюс

первого порядка в точке . Теперь используя теорему вычетов, находим h(z).

(41)

Мы видим что, рассматриваемый рекурсивный фильтр имеет частотную характеристику, у которой отличен от нуля только один элемент. Поэтому частотная характеристика является конечной импульсной характеристикой.
Таким образом, рассматриваемый рекурсивный фильтр является КИХ – фильтром.

Слайд 38

Как мы показали выше, все КИХ – фильтры являются устойчивыми фильтрами.

Как мы показали выше, все КИХ – фильтры являются устойчивыми

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

(42)

Слайд 39

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

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

Подведем промежуточный итог, касающийся

свойств рекурсивных и нерекурсивных ЛДФ
Слайд 40

Слайд 41

Аналоговые фильтры Как мы говорили ранее, в основе многих методов проектирования

Аналоговые фильтры

Как мы говорили ранее, в основе многих методов

проектирования дискретных фильтров лежат методы проектирования аналоговых фильтров. Поэтому в этом разделе мы познакомимся с некоторыми известными аналоговыми фильтрами.
Аналоговый фильтр относится к линейным системам с постоянными параметрами (ЛПП). Вспомним, что ЛПП (а значит и фильтр) описывается функцией времени h(t) - импульсной характеристикой. Эта функция определяется соотношениями.
Слайд 42

(43) Здесь x(t) - входящий сигнал, а y(t) - выходящий сигнал.

(43)

Здесь x(t) - входящий сигнал, а y(t) - выходящий сигнал.

Комплексный

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

(44)

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

(45)

Слайд 43

Классификация фильтров по их АЧХ 1. Фильтры нижних частот (ФНЧ, английский

Классификация фильтров по их АЧХ

1. Фильтры нижних частот (ФНЧ, английский

термин – low-pass filter), пропускающие частоты, меньшие некоторой частоты среза ω0 .
На рисунке показана АЧХ идеального фильтра нижних частот.
Слайд 44

2. Фильтры верхних частот (ФВЧ, английский термин – high-pass filter), пропускающие

2. Фильтры верхних частот (ФВЧ, английский термин – high-pass filter), пропускающие

частоты, больше некоторой частоты среза ω0 .
На рисунке показана АЧХ идеального фильтра верхних частот.
Слайд 45

3. Полосовые фильтры (ПФ, английский термин – band-pass filter), пропускающие частоты

3. Полосовые фильтры (ПФ, английский термин – band-pass filter), пропускающие частоты

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

и шириной полосы пропускания.

На рисунке показана АЧХ идеального полосового фильтра.

Слайд 46

Слайд 47

3. Режекторные фильтры (РФ, английский термин – band-stop filter). Имеются другие

3. Режекторные фильтры (РФ, английский термин – band-stop filter). Имеются другие

названия таких фильтров – заграждающий фильтр, фильтр-пробка, полосно-задерживающий фильтр. Режекторные фильтры – это фильтры пропускающие на выход все частоты, кроме лежащих в некотором диапазоне [ω1, ω2] . Такие фильтры характеризуются средней частотой ωc и шириной полосы задержки Δω .
На рисунке показана АЧХ идеального режекторного фильтра.
Слайд 48

Слайд 49

Передаточная функция аналогового фильтра Передаточная функция H(s) аналогового фильтра является изображением

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

Передаточная функция H(s) аналогового фильтра является изображением Лапласа

импульсной характеристики h(t).

(46)

Здесь s - комплексное число.

Слайд 50

Основным уравнением аналогового фильтра является следующее дифференциальное уравнение. (47) В этом

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

(47)

В

этом дифференциальном уравнении ai , bi - постоянные коэффициенты, которые определяют свойства фильтра.
Слайд 51

Если эти коэффициенты заданы, то передаточная функция определятся следующей формулой. (48)

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


(48)

Комплексный коэффициент передачи фильтра и то передаточная функция H(s) связаны простым соотношением.

(49)

Слайд 52

Поэтому АЧХ и ФЧХ аналогового фильтра легко найти, если мы знаем

Поэтому АЧХ и ФЧХ аналогового фильтра легко найти, если мы

знаем переходную функцию.
Разложив числитель и знаменатель функции передачи (48) на элементарные множители, мы получаем функцию передачи в следующем виде.

(50)

Здесь k = bm /an - коэффициент усиления, zi – нули функции передачи, pi - полюсы функции передачи. Поэтому для определения свойств аналогового фильтра вместо коэффициентов ai ,bi основного дифференциального уравнения, можно использовать нули и полюса zi , pi функции передачи.

Слайд 53

Фильтры Баттерворта Фильтры Баттерворта это фильтры нижних частот. Передаточная функция H(s)

Фильтры Баттерворта

Фильтры Баттерворта это фильтры нижних частот. Передаточная функция

H(s) определяется n полюсами, которые задаются формулой.

(51)

Здесь ω0 - некоторая заданная частота, называемая частотой среза. Полюсы (51) лежат в комплексной плоскости на окружности радиуса ω0 , потому что для них выполняется условие.

(52)

Слайд 54

Число n определяет порядок фильтра Баттерворта. На рисунке показано расположение полюсов

Число n определяет порядок фильтра Баттерворта. На рисунке показано расположение

полюсов для фильтра Баттерворта 5-го порядка в комплексной s - плоскости. Для простоты частота среза взята раной единице ω0 = 1 .
Слайд 55

Передаточная функция H(s) для фильтра Баттерворта конструируется из полюсов следующим образом.

Передаточная функция H(s) для фильтра Баттерворта конструируется из полюсов следующим

образом.

(53)

Здесь k0 - нормировочный множитель.

Используя связь (49) между передаточной функцией H(s) и комплексным коэффициентом передачи K(ω), находим комплексный коэффициент передачи.

(54)

Слайд 56

Теперь из формулы (54) можно найти АЧХ и ФЧХ фильтра Баттерворта

Теперь из формулы (54) можно найти АЧХ и ФЧХ фильтра

Баттерворта n - го порядка. Простые, но громоздкие вычисления позволяют получить для АЧХ простую аналитическую формулу. Если нормировочный множитель положить равным

то АЧХ фильтра Баттерворта будет выражаться формулой.

(55)

Слайд 57

На рисунке показано АЧХ фильтра Баттерворта 5-го порядка Для простоты частота

На рисунке показано АЧХ фильтра Баттерворта 5-го порядка Для простоты частота

среза взята раной единице ω0 = 1 .
Слайд 58

На рисунке показано ФЧХ фильтра Баттерворта 5-го порядка Для простоты частота

На рисунке показано ФЧХ фильтра Баттерворта 5-го порядка Для простоты

частота среза взята раной единице ω0 = 1 .

Для ФЧХ такой простой формулы как для АЧХ не существует. Поэтому надо брать формулу (54) и по ней находить ФЧХ.

Слайд 59

Фильтр Чебышева первого рода Фильтр Чебышева первого рода является фильтром нижних

Фильтр Чебышева первого рода

Фильтр Чебышева первого рода является фильтром

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

(56)

Слайд 60

Здесь ω0 - частота среза, Tn (x)- полином Чебышева n -

Здесь ω0 - частота среза, Tn (x)- полином Чебышева n

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

(57)

На рисунке показано расположение полюсов для фильтра Чебышева первого рода 5-го порядка в комплексной s - плоскости. Уровень пульсаций в полосе пропускания взят 0.5 дБ. Для простоты частота среза взята раной единице ω0 = 1.

Слайд 61

На следующих двух рисунках показаны АЧХ ФЧХ для фильтра Чебышева первого

На следующих двух рисунках показаны АЧХ ФЧХ для фильтра Чебышева

первого рода 5-го порядка. Уровень пульсаций в полосе пропускания взят 0.5 дБ. Для простоты частота среза взята раной единице ω0 = 1 .
Слайд 62