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

advertisement
2. Цифровые фильтры
2.1. Свойства Z-преобразования
Прямым Z-преобразованием дискретной последовательности xn, где n = 0,1, 2..,
называется функция комплексной переменной z, определяемая следующим соотношением

X(z)   x n z  n .
(2.1)
n 0
Функция X(z) определена для тех значений z, при которых ряд сходится.
Здесь и в дальнейшем последовательность отсчётов обозначается строчной, а ее Zпреобразование той же прописной буквой.
Примеры определения Z-преобразований трех простых последовательностей
приведены в таблице 2.1.
Таблица 2.1
Название
Последовательность
Z-преобразование
последовательности
последовательности
Единичный
отсчет
1 при n  0
xn  
0 при n  0
Единичный
скачок
1 при n  0
xn  
0 при n  0
Показательная
функция
x n  n при   1, n  0
X(z)=1 при z  0
X(z) 
X(z) 
1
при z  1
1  z 1
1
при z  
1   z 1
Рассмотрим основные свойства прямого Z-преобразования.
1.Линейность. Пусть последовательность yn представляет взвешенную сумму двух
последовательностей x1n и x2 n
y n  1x1 n   2 x 2 n ,
где 1 , 2  постоянные весовые коэффициенты.
Тогда Z-преобразование последовательности yn определяется следующим
соотношением
Y(z)  1X1 (z)  2X 2 (z) .
(2.2)
Таким
образом,
Z-преобразование
взвешенной
суммы
двух
последовательностей
равно
взвешенной
сумме
Z-преобразований
этих
последовательностей.
2.Сдвиг последовательностей. Пусть последовательность yn представляет собой
сдвинутую (задержанную) на m отсчетов последовательность xn (рисунок 2.1)
y n  x n m .
Тогда Z-преобразование Y(z) последовательности yn выражается через Zпреобразование X(z) последовательности xn следующим образом
Y(z)  z  m X(z) .
(2.3)
Таким
образом,
Z-преобразование
последовательности,
сдвинутой
относительно исходной на m отсчетов, равно Z-преобразованию исходной
последовательности, умноженной на z –m.
1
xn
n
0 1
2
3 4
5
6
7 8
yn
n
0
1
2
3 4
5
6 7
8
Рисунок 2.1 – Последовательность yn сдвинута относительно xn
на 2 отсчета (2 интервала дискретизации)
3.Дискретная свертка двух последовательностей. Дискретной сверткой двух
последовательностей xn и hn называется последовательность yn, определяемая следующим
соотношением
n
n
k 0
k 0
y n   h k x n k   x k h n k .
(2.4)
Z-преобразование Y(z) дискретной свертки yn двух последовательностей равно
произведению Z -преобразований H(z) и X(z) исходных последовательностей hn и xn
Y(z)  H(z) X(z) ,
(2.5)



n 0
n 0
n 0
где Y(z)   y n z  n , H(z)   h n z  n , X(z)   x nz  n .
2.2. Импульсная характеристика цифрового фильтра. Понятие о рекурсивных и
нерекурсивных цифровых фильтрах, БИХ- и КИХ-фильтрах
Цифровым фильтром дискретного сигнала называется линейная частотноизбирательная система, реализуемая на основе вычислительного устройства.
Пусть при действии на входе цифрового фильтра последовательности отсчетов xn
на его выходе действует последовательность yn .
Если n-ый отсчет выходного сигнала фильтра yn зависит только от отсчетов
входного сигнала в данный и предшествующие моменты дискретного времени xn, xn-1 ..и
т.д., то такой фильтр называется нерекурсивным.
Если n-ый отсчет выходного сигнала фильтра yn зависит не только от отсчетов
входного сигнала в данный и предшествующие моменты дискретного времени xn, xn-1 ..и
т.д., но и от отсчетов выходного сигнала в предшествующие моменты времени, то такой
фильтр называется рекурсивным.
Импульсной характеристикой цифрового фильтра называется выходной сигнал
фильтра при действии на его входе единичного отсчета и нулевых начальных условиях.
На рисунке 2.2 показаны входной сигнал фильтра в виде единичного отсчета xn и
реакция фильтра на этот сигнал – импульсная характеристика hn.
Фильтр с конечной импульсной характеристикой называется КИХ-фильтром
(КИХ-конечная импульсная характеристика). Фильтр с бесконечной импульсной
характеристикой называют БИХ-фильтром.
2
xn
1
hn 0 1
2
3
4
5
6
7
8
9 n
0
2
3
4
5
6
7
8
9 n
1
Рисунок 2.2 – Единичный отсчет xn и импульсная характеристика hn
2.3. Определение выходного сигнала фильтра по входному сигналу и импульсной
характеристике
Определение выходного сигнала цифрового фильтра по входному сигналу и
импульсной характеристике основано на определении импульсной характеристики и
принадлежности фильтра к линейным системам, для которых справедлив принцип
суперпозиции.
На рисунке 2.3 показаны три отсчета входного сигнала фильтра x0=3, x1=1, x2=2 и
его импульсная характеристика, представленная четырьмя отсчетами h0=2, h1= -1, h2=1,
h3=0.5.
Согласно принципу суперпозиции для определения выходного сигнала фильтра
можно определить сигналы на выходе фильтра при действии каждого отсчета,
изолированного от всех остальных, а затем сложить полученные сигналы.
Определим сигнал на выходе фильтра при действии на входе одного отсчета x0.
Если бы вместо x0 на входе действовал единичный отсчет, то выходным сигналом была
бы импульсная характеристика фильтра. Поскольку данный отсчет отличается от
единичного отсчёта в x0 раз, а фильтр является линейной системой, то и отсчеты
выходного сигнала будут в x0 раз больше отсчетов импульсной характеристики: h0x0, h1x0,
h2x0, h3x0.
Определим сигнал на выходе фильтра при действии на входе одного отсчета x1.
Если бы вместо x1 на входе действовал единичный отсчет, то выходным сигналом была
бы импульсная характеристика фильтра, смещенная на один отсчет вправо, а реакцией
фильтра на отсчет x1 будет последовательность отсчетов h0x1, h1x1, h2x1, h3x1, смещенная
относительно последовательности h0x0, h1x0, h2x0, h3x0 на один отсчет вправо.
Аналогичным образом определяется реакция фильтра на остальные отсчеты
сигнала. Из рисунка видно, что отсчеты выходного сигнала определяются следующими
соотношениями:
y0  h 0 x 0 ,
y1  h 0 x1  h1 x 0 ,
y 2  h 0 x 2  h1 x1  h 2 x 0 ,
y 3  h 0 x 3  h1 x 2  h 2 x1  h 3 x 0 ,
y 4  h 0 x 4  h1 x 3  h 2 x 2  h 3 x1  h 4 x 0 ,
y5  h 0 x 5  h1 x 4  h 2 x 3  h 3 x 2  h 4 x1  h 5 x 0 .
В выражении для y3 первое слагаемое равно нулю, т.к. x3 = 0.
В выражении для y4 два первых и последнее слагаемое равны нулю, т.к. x3 = 0, x4
= 0, h4 = 0.
3
В выражении для y5 три первых и два последних слагаемых равны нулю, т.к. x3 =
0, x4 = 0, x5 = 0, h4 = 0, h5 = 0.
xn
x0 =3
x2 = 2
x1 = 1
0
1
2
n
hn
h 2= 1
h 0= 2
0
1
h 1= -1
h 3= 0.5
2
3
n
h 0x0=6
-1 1
h 2x0=3
h 3x0=1.5
n
h 1x0=-3
h 0x1=2
h 2x1=1
h 3x1=0.5
n
h 1x1=-1
h 0x2=4
h 2x2=2
h 3x2=1
n
h 1x2=-2
yn
y2=6
y4 = 2.5
y0=6
y3=0.5
y5 = 1
n
y1= -1
Рисунок 2.3 – Определение выходного сигнала фильтра по входному
сигналу и импульсной характеристике
В общем случае n – ый отсчет выходного сигнала определяется следующими
соотношениями:
n
y n   h k x nk
(2.6)
k 0
4
n
y n   x k h nk .
или
(2.7)
k 0
Из них следует, что выходной сигнал фильтра представляет собой дискретную
свертку входного сигнала и импульсной характеристики.
На рисунке 2.4 представлена графическая интерпретация соотношения (2.6) при
конечной импульсной характеристике нерекурсивного цифрового фильтра, содержащей
N+1 отсчет. В дальнейшем такое графическое представление алгоритма фильтрации
будем называть схемой фильтра.
h0
yn
x
n
z -1
h1
z -1
h2
z -1
hN
x n-1
x n-2
x n-N
Рисунок 2.4 – Нерекурсивный цифровой фильтр
2.4. Системная функция цифрового фильтра. Формы программной реализации
фильтра
Системной функцией цифрового фильтра называется отношение Zпреобразования выходного сигнала фильтра к Z-преобразованию входного сигнала
Y(z)
.
H(z) 
X(z)
Воспользовавшись (2.6) и теоремой о дискретной свертке (раздел 2.1), выразим Zпреобразование Y(z) выходного сигнала фильтра yn через Z-преобразование X(z) входного
сигнала xn
Y(z) = H(z) X(z),



n 0
n 0
n 0
где Y(z)   y n z  n , H(z)   h n z  n , X(z)   x nz  n .
Из последних соотношений следует, что системная функция H(z) представляет
собой Z-преобразование импульсной характеристики цифрового фильтра.
Полюсом системной функции называется значение комплексной переменной z,
при котором системная функция H(z) стремится к бесконечности.
Нулем системной функции называется значение комплексной переменной z, при
котором системная функция H(z) равна нулю.
Рассмотрим формы программной реализации фильтра:
1. Прямая форма
На рисунке 2.5 представлен алгоритм функционирования цифрового фильтра при
прямой форме реализации. Прямая форма следует из определения фильтра как линейной
системы. Следовательно, n – ый отсчет выходного сигнала фильтра yn должен быть связан
линейными соотношениями с отсчетами входного сигнала в данный и предшествующие
моменты дискретного времени xn, xn-1, ..xn-N и отсчетами выходного сигнала в
предшествующие моменты времени yn-1, yn-2, .. yn-N. Соответствующие коэффициенты
пропорциональности B0, B1, .. BN, A1, A2, .. AN определяют свойства фильтра.
5
В0
xn
z -1
В1
yn
-A1
z -1
xn-1
yn-1
z-1
В2
-A2
z-1
xn-2
yn-2
z
xn-N
-1
ВN
-AN
z
-1
yn-N
Рисунок 2.5 – Прямая форма программной реализации фильтра
Согласно схеме
yn  B0 x n  B1x n 1  ..BM x n  N  A1yn 1  A 2 yn 2 ..  A N y n  N . (2.8)
Выразим Z - преобразование выходного сигнала Y(z) через Z-преобразование
входного сигнала
Y(z)  B0 X(z)  z 1B1X(z)  z 2B2X(z)  ..  z  N BN X(z) 
 z 1A1Y(z)  z 2 A 2 Y(z)  ..  z  N A N Y(z)
Из последнего соотношения получим
B  B1z 1  B2 z 2  ..  BN z  N
H(z)  0
.
(2.9)
1  A1z 1  A 2 z 2  ..  A N z  N
Таким образом, системная функция цифрового фильтра в общем случае
представляет собой дробно-рациональную функцию. Полином числителя описывает
нерекурсивную часть фильтра, а полином знаменателя – рекурсивную.
Чтобы найти нули системной функции, нужно полином числителя приравнять
нулю и найти корни полученного уравнения.
Чтобы найти полюсы системной функции, нужно полином знаменателя приравнять
нулю и найти корни полученного уравнения.
Отметим, что знаки перед коэффициентами A в выражении для системной функции
и в разностном уравнении (2.8) противоположны.
2.Каноническая форма.
Представим выражение (1.9) в виде произведения двух функций
H(z)  H A (z) H B (z) ,
1
где H A (z) 
,
1
1  A1z  A 2 z 2  ..A N z  N
(2.10)
H B (z)  B0  B1z 1  ..BN z  N .
Согласно (2.10) цифровой фильтр с системной функцией H(z) можно представить
в виде последовательного соединения двух фильтров с системными функциями HA(z) и
HB(z) (рисунок 2.6).
xn
yn
vn
HA(z)
HB(z)
Рисунок 2.6 – Представление фильтра с прямой формой реализации в виде
последовательного соединения двух фильтров
6
Действительно,
V(z)  X(z)H A (z),
.
Y(z)  V(z)H B (z)  X(z)H A (z)H B (z)  X(z)H(z)
Заменив укрупненный алгоритм рисунка 2.6 детальным, получим схему фильтра,
изображенную на рисунке 2.7.
xn
В0
vn
yn
-A1
z -1
z -1
В1
vn-1
-A2
z-1
z-1
В2
vn-2
-AN
z-1
z-1
ВN
vn-N
Рисунок 2.7 – Детальный алгоритм представления фильтра с прямой
реализацией в виде последовательного соединения двух фильтров
Из рисунка видно, что для хранения одних и тех же переменных используются две
линии задержки, поэтому одну из них можно удалить. При этом схема фильтра
преобразуется к виду, представленному на рисунке 2.8. Это и есть каноническая форма
программной реализации фильтра.
В0
vn
xn
yn
-A1
-A2
-AN
z -1
z-1
vn-1
В1
vn-2
В2
vn-N
ВN
z-1
Рисунок 2.8 – Каноническая форма программной реализации фильтра
Достоинством канонической формы является в два раза меньшее количество
элементов задержки, следовательно, ячеек памяти вычислительного устройства.
На рисунке 2.8 показана каноническая форма фильтра N-го порядка на одной
линии задержки, состоящей из N элементов. Однако обычно вместо структуры,
изображенной на рисунке 2.8, используется параллельное или последовательное
соединение звеньев второго порядка. Такое
представление фильтра связано с
возможностью представления системной функции (2.9) в виде произведения или суммы
системных функций с полиномами второго порядка в числителе и знаменателе
B0L  B1Lz 1  B2Lz 2
H(z)  
,
1
 A 2Lz 2
L 1 1  A1L z
Lmax
(2.11)
7
B0L  B1Lz 1  B2Lz 2
H(z)  
,
1
 A 2Lz 2
L 1 1  A1L z
Lmax
(2.12)
где L – порядковый номер звена, Lmax – максимальное значение номера звена
N
 2 , åñëè N  ÷åòí î å ÷èñëî ,
L max  
1  N  1 , åñëè N  í å÷åòí î å ÷èñëî .

2
При четном N фильтр состоит из N/2 звеньев второго порядка, при нечетном N
фильтр состоит из одного звена первого порядка и (N-1)/2 звеньев второго порядка.
Системная функция звена первого порядка отличается от системной функции звена
второго порядка тем, что коэффициенты B2 и A2 равны нулю.
Соотношению (2.11) соответствует схема рисунка 2.9а, а соотношению (2.12) –
схема рисунка 2.9б.
а)
H1(z)
H2(z)
HLmax(z)
б)
H1(z)
H2(z)
HLmax(z)
Рисунок 2.9- Последовательное (а) и параллельное (б) соединение
звеньев фильтра
Типовая схема звена второго порядка приведена на рисунке 2.10. На входе звена
показан масштабный коэффициент ML (как правило, меньше единицы),
предотвращающий появление в процессе вычислений значений сигналов фильтра,
выходящих за пределы разрядной сетки вычислительного устройства.
xn
ML
vn
-A1L
-A2L
z -1
z-1
B0L
yn
vn-1 В1L
vn-2 В
2L
Рисунок 2.9 – Типовое звено второго порядка
8
2.5.Частотная характеристика цифрового фильтра
Комплексным коэффициентом передачи фильтра K является отношение
комплексной амплитуды Y выходного сигнала фильтра к комплексной амплитуде
входного сигнала X
Y
K .
X
Частотной характеристикой цифрового фильтра K( j) называется зависимость
комплексного коэффициента передачи фильтра от частоты.
Амплитудно-частотной характеристикой (АЧХ) K( ) называется зависимость
модуля комплексного коэффициента передачи от частоты
K( )  K( j) .
Фазочастотной характеристикой (ФЧХ) называется зависимость аргумента
комплексного коэффициента передачи фильтра от частоты.
()  arg(K( j)) .
Для определения комплексного коэффициента передачи фильтра подадим на вход
фильтра с прямой формой реализации (рисунок 1.5) комплексный сигнал с единичной
амплитудой
jnT
x n  e Д  cos( nTД )  jsin( nTД ) .
Согласно определению комплексного коэффициента передачи комплексный
выходной сигнал должен быть равен
j nT
yn  K( j)e Д .
Из схемы рисунка 2.5 следует, что выходной комплексный сигнал фильтра
определяется следующим соотношением
j nT
j nT
j (n 1)TД
j( n  2)TД
j (n  N)TД
y n  K( j) e Д  B0 e Д  B1 e
 B2 e
 ..  B N e

.
j (n 1)TД
j (n  2)TД
j (n  N)TД
 A1K( j) e
 A 2 K( j) e
 ..  A N K( j) e
Из последнего соотношения получим
 jT
2 jTÄ
 N jTÄ
B0  B1e Ä  B2e
 ..  BN e
(2.13)
K( j) 
 jT
2 jTÄ
 N jTÄ
1  A1e Ä  A 2e
 ..  A Ne
Сравнивая последнее соотношение с выражением для системной функции
цифрового фильтра, можно сформулировать правило определения комплексного
коэффициента передачи при известной системной функции фильтра: для нахождения
комплексного коэффициента передачи нужно в выражении для системной функции
jT
заменить z на e Д :
z  e j ,
(2.14)
f
где   TД  2f N , f N 
- нормированная частота – отношение текущей частоты f к
FД
частоте дискретизации FД.
2.6.Устойчивость цифровых фильтров
Рассмотрим критерии устойчивости цифровых фильтров.
1.Критерий «ОВ-ОВ» («Ограниченный вход – ограниченный выход»)
Цифровой фильтр устойчив, если при ограниченном входном сигнале
выходной сигнал фильтра также ограничен.
9
Условие ограниченности входного сигнала определяется соотношением x n  R ,
где R   , а условием ограниченности выходного сигнала является yn   .
Непосредственное использование этого критерия весьма затруднительно, т.к.
требует определения значений отсчетов выходного сигнала при всех возможных
значениях отсчетов входного сигнала. Поэтому требуются критерии, позволяющие
оценить устойчивость фильтра на основании его характеристик.
2.Критерий оценки устойчивости по импульсной характеристике фильтра
В разделе 2.3 было доказано, что выходной сигнал фильтра представляет собой
дискретную свертку входного сигнала и импульсной характеристики фильтра
n
y n   h k x n k .
k 0
Абсолютное значения отсчетов выходного сигнала удовлетворяет неравенству
n
y n   h k x n k .
k 0
При x n  R справедливо неравенство
n

k 0
k 0
 h k x nk  R  h k .
Следовательно,

yn  R  hk .
k 0
Таким образом, чтобы обеспечить выполнение условия
выполнить условие

h
k 0
k
.
yn   , достаточно
(2.15)
Последнее соотношение определяет критерий устойчивости цифрового фильтра,
который формулируется так: цифровой фильтр устойчив, если сумма абсолютных
значений отсчетов его импульсной характеристики конечна.
Из этого критерия следует, что все фильтры с конечной импульсной
характеристикой абсолютно устойчивы.
В качестве примера воспользуемся критерием (2.15) для проверки устойчивости
фильтра, импульсная характеристика которого бесконечна и описывается соотношением

TД
k
hk  e  ,
где    – положительная константа, от которой зависит скорость убывания
отсчетов импульсной характеристики.
Учитывая, что h k  h k , получим

e

TД

k
k 0

1

1 e

TД
.

Tд
Так как e   1 , то фильтр устойчив.
3.Критерий оценки устойчивости по системной функции фильтра
В разделе 2.4 показано, что системная функция представляет собой Zпреобразование импульсной характеристики фильтра

H(z)   h n z  n .
n 0
Модуль системной функции удовлетворяет неравенству
10
n
1
H(z)   h n   .
n 0
z
При z  1 справедливо неравенство

n

1
h

  hn .

n 
z
n 0
n 0
 

При

h
n 0
n
  и при z  1 модуль системной функции H(z)   . Последнее
соотношение означает, что в устойчивом цифровом фильтре должны отсутствовать
полюсы системной функции в области комплексной переменной z, которая удовлетворяет
неравенству z  1 .
Следовательно, если полюсы существуют, то в устойчивом фильтре они должны
располагаться в области комплексной плоскости, для которой выполняется условие z  1 .
Поэтому критерий устойчивости, связанный с системной функцией фильтра,
формулируется следующим образом: цифровой фильтр устойчив, если полюсы
системной функции располагаются внутри круга единичного радиуса с центром в
начале координат ( z  1 ).
Оценим устойчивость фильтра, системная функция которого описывается
соотношением
1
H(z) 
,
1  A1z 1
где A1= - 0.5.
Приравняем знаменатель системной функции нулю и определим корень
полученного уравнения, который является координатой полюса
1  A1z 1  0,
z   A1 .
На рисунке 2.10 показан круг единичного радиуса и полюс системной функции,
располагающийся внутри круга. Следовательно, фильтр устойчив.
Im(z)
1
1
Re(z)
Рисунок 2.10 – Область устойчивости и полюс системной функции
2.7. Коэффициенты системной функции устойчивого звена второго порядка
Системная функция звена второго порядка определяется соотношением
B  B1z 1  B2z 2
H(z)  0
.
1  A1z 1  A 2z 2
Для определения полюсов системной функции приравняем знаменатель нулю и
найдем корни полученного квадратного уравнения
11
1  A1z 1  A 2z 2  0,
2
A
A 
z1,2   1   1   A 2
2
 2 
Фильтр реализуется в виде звеньев второго порядка
в случае комплексносопряженных корней, т.е. при
2
 A1 
(2.16)
   A2  0 .
 2 
В этом случае корни уравнения определяются следующим соотношением
2
A1
A 
 j A2   1  .
2
 2 
Из последнего соотношения находим
z1,2  A 2 .
z1,2  
Условием устойчивости звена является
z1,2  1 .
Поэтому коэффициент
удовлетворять условию
A2
устойчивого
звена
второго
порядка
должен
0  A2  1 .
(2.17)
Из неравенств (2.16) и (2.17) следует неравенство для коэффициента A1
(2.18)
A1  2 .
12
Download