Введение в теорию сигналов. Цифровая обработка сигналов. Учебное пособие
.pdfДругой формой математической модели непрерывного фильтра может служить интегральное уравнение свертки
|
t |
|
|
|
|
y(t) w( ) x(t )d . |
|
|
(2.2) |
|
0 |
|
|
|
Математической моделью рекурсивного фильтра является линейное |
||||
разностное уравнение n -го порядка |
|
|
|
|
n |
m 1 |
|
|
|
yk ai |
y(k (n 1 i)) bi x(k (n 1 i)) , |
m |
n , |
(2.3) |
i 1 |
i 1 |
|
|
|
представляющее собой дискретный аналог дифференциального уравнения (2.1), или соответствующая уравнению (2.3) передаточная функция
|
|
|
|
m 1 |
z(i 1) |
|
|
Y (z) |
|
|
bi |
||
W (z) |
|
|
i 1 |
|
, |
|
|
|
n |
|
|||
|
X (z) |
|
|
|
|
|
|
|
|
zn ai z(i 1) |
|||
|
|
|
|
i 1 |
|
|
где z – переменная Z-преобразования; z exp ( t s) , |
||||||
|
|
|
|
|
|
|
X (z) xk z k , Y (z) |
yk z k |
, W (z) wk z k , |
||||
k 0 |
|
k 0 |
|
|
k 0 |
|
а wk w(k t) – ИХ фильтра, имеющая бесконечную длительность.
Передаточнаяфункцияэтихфильтровсодержитполюсаинули, ивыход фильтра определяется помимо входа еще и предыдущим выходом.
Математическоймодельюнерекурсивного фильтраслужитдискретная форма интегрального уравнения свертки (2.2):
k
yk wi x(k i) , i 0
причем ИХ нерекурсивного фильтра является финитной, т. е. имеет ко-
нечную длительность – wk 0 при k 0 , N 1 , |
wk 0 при k |
N , |
и его передаточная функция записывается в виде |
|
|
41
N 1 |
|
W (z) wk z k . |
(2.4) |
k 0
Эти фильтры физически нереализуемы, так как их передаточная функция имеет порядок n 0 , содержит нули, но не содержит полюсов и выход фильтра определяется только входом.
Частотные характеристики непрерывных фильтров являются непериодическими функциями частоты и определяются на интервале [0, ) . Частотные характеристики дискретных фильтров пред-
ставляют собой периодические с периодом 2 в , в / t , функции
частоты , в результате чего они полностью определяются на интервале [0, в] .
В зависимости от характера обработки входного сигнала фильтры обладают различными АЧХ. При этом могут быть выделены четыре основных типа фильтров (низкочастотный (НЧФ), полосовой (ПФ), режекторный (РФ) и высокочастотный (ВЧФ)), а также представляющий значительный практический интерес дифференцирующий фильтр (ДФ). Идеализированные АЧХ Aи( ) таких фильтров в интервале [0, в]
приведены на рис. 2.2–2.6 соответственно.
На каждой АЧХ можно выделить следующие полосы: 1) полоса пропускания
для НЧФ и ДФ [0, ф] ;
для ПФ [ ф, ф] ;
для РФ [0, ф] и [ ф, в] ;
для ВЧФ [ ф, в] ; 2) переходная полоса
для НЧФ и ДФ [ ф, δ] ; для ПФ [ δ, ф] и [ ф, δ] ; для РФ [ ф, δ] и [ δ, ф] ; для ВЧФ [ δ, ф] ;
42
Aи(ω)
Kф
Kф
δKф |
Δωп |
|
|
|
|
ω |
|
|
|
|
|||
|
|
|
|
|
||
0 |
|
|
|
|
|
|
|
Ωф Ωδ ωв |
|
||||
Рис. 2.2. АЧХ НЧФ
Aи(ω) |
|
Kф |
Kф |
|
|
Δωп |
Δωп |
δKф
0 |
Ωф Ωδ |
ωδ ωф ωв ω |
Рис. 2.4. АЧХ РФ
Aи(ω)
Aи(ω) |
|
|
|
|
|
Kф |
|
Kф |
|
|
|
|
|
|
|
|
|
|
|
Δωп |
|
Δωп |
|
δKф |
|
|
|
|
|
|
0 ωδ ωф |
Ωф Ωδ ωв |
ω |
||
|
|
Рис. 2.3. АЧХ ПФ |
|
||
Aи(ω) |
|
|
|
|
|
Kф |
Kф |
|
|
|
|
|
|
|
|
|
|
|
|
Δωп |
|
|
|
δKф |
|
|
|
|
|
|
0 |
ωδ |
ωв |
ω |
|
|
Рис. 2.5. АЧХ ВЧФ |
|
|||
Δωп 
δ
|
|
ω |
0 |
Ωф Ωδ |
ωв |
Рис. 2.6. АЧХ ДФ
3) полоса задерживания
для НЧФ и ДФ [ δ, в] ;
43
для ПФ [0, δ] и [ δ, в] ;
для РФ [ , ] ; для ВЧФ [0, δ] .
Граничные частоты ф и ф характеризуют тип фильтра и кон-
кретный вид его идеальной АЧХ (при 0 , |
0 , п 0 ). |
|
Относительные неравномерности и |
характеризуют |
ши- |
рину Kф , и Kф зон допустимых изменений реальной АЧХ |
A( ) |
|
фильтравполосахпропусканияизадерживания, атемсамымсовместно с шириной п переходных полос определяют качество АЧХ фильтра.
С уменьшением , и п качество фильтра улучшается.
При проектировании фильтра удобно характеризовать ширину п переходных полос АЧХ через относительный параметр в виде
( 1) ф |
для НЧФ, |
|
|||
|
|
|
|
для ПФ, |
|
( 1)( ф ф) |
(2.5) |
||||
п |
|
/ )( ф ф) |
|
||
(1 |
1 |
для РФ, |
|
||
(1 1 |
/ |
) |
для ВЧФ, |
|
|
|
|
|
ф |
|
|
где 1 – относительная частота, характеризующая переходную полосу (в случае идеального фильтра 1 ).
В заключение отметим два характерных режима фильтрации:
1) режим реального времени, при котором для вычисления выходного сигнала фильтра в текущий k -й момент могут быть использованы значения входного сигнала только в k -й и предыдущие моменты;
2) режим обработки, при котором весь массив дискретных значений входного сигнала хранится в памяти ЦВМ, и для вычисления выходногосигналафильтравлюбоймоментмогутбытьиспользованызначения входного сигнала как в предыдущие, так и в последующие моменты.
44
2.2.2. ПОГРЕШНОСТИ ФИЛЬТРАЦИИ
Выходной сигнал реального фильтра в процессе его функционирования всегда отличен от идеального. Это происходит за счет наличия погрешностей фильтрации.
Составляющие погрешности
1. Динамическая, или переходная, составляющая обусловлена динамикой фильтра и может принимать достаточно большие значения по уровню на интервале времени t [0, Tw ] , Tw – длительность ИХ (в слу-
чае нерекурсивного фильтра) либо |
эффективная длительность ИХ |
(в случае устойчивого рекурсивного |
фильтра). При t Tw переходная |
погрешность отсутствует полностью у нерекурсивных фильтров и принимает малые значения у рекурсивных. Уровень этой погрешности соответствует уровню усечения w ИХ w(t) , определяющему эффектив-
ную длительность ИХ. Обычно w = 0,01…0,05 от max w(t) .
t [0, )
2. Шумовая составляющая погрешности возникает из-за составляющих спектра помехи вне полосы пропускания фильтра и зависит от неравномерности и ширины переходной полосы Δωп АЧХ фильтра.
3. Смешанная составляющая обусловлена составляющими спектра полезного сигнала вне полосы пропускания АЧХ и составляющими АЧХ шума | x( j ) | в полосе пропускания фильтра. Уровень этой со-
ставляющей погрешности зависит от частоты ф .
4.Амплитудная составляющая возникает из-за неидеальности АЧХ
вполосе пропускания и определяется неравномерностью .
5.Фазовая составляющая обусловлена неидеальностью ФЧХ в полосе пропускания АЧХ. ФЧХ нерекурсивных фильтров строго линейна:
( ) ф , поэтому если в выходной сигнал ввести запаздывание величиной ф , то фазовая составляющая погрешности может быть све-
дена к нулю. ФЧХ рекурсивных фильтров нелинейная, и потому в пределах полосы пропускания АЧХ она линеаризуется, и так же, как и в нерекурсивных фильтрах, вводится запаздывание, однако свести к нулю фазовую составляющую погрешности здесь не удается.
45
6. Дискретнаясоставляющаявозникаетиз-задискретногопредстав- ления обрабатываемых реализаций сигнала. Проявляется это в отражении и наложении спектров относительно частоты в
t при невер-
ном выборе шага t .
7. Вычислительная составляющая погрешности возникает из-за округления числовых значений в реализации сигнала в процессе обработки.
2.2.3. ИСХОДНЫЕ ДАННЫЕ ПРИ ПРОЕКТИРОВАНИИ ФИЛЬТРОВ
Как правило, при проектировании фильтра предъявляются требования к его АЧХ, а необходимая ФЧХ обеспечивается выбором класса фильтра.
Исходными данными, определяющими вид АЧХ фильтра и ее качество, являются следующие:
1)относительная неравномерность АЧХ в полосе пропускания;
2)относительная неравномерность δ АЧХ в полосе задерживания;
3)относительная частота , определяющая ширину Δωп переход-
ных полос АЧХ; 4) верхняя граничная частота ф полосы пропускания АЧХ в слу-
чае НЧФ, ПФ, РФ и ДФ; 5) нижняя граничная частота ωф полосы пропускания АЧХ в случае
ПФ, РФ и ВЧФ;
6)коэффициент усиления kф в полосе пропускания АЧХ;
7)шаг дискретизации по времени t .
2.3.НЕРЕКУРСИВНЫЕ ФИЛЬТРЫ
2.3.1. ХАРАКТЕРИСТИКИ ФИЛЬТРОВ
Нерекурсивные фильтры имеют конечную длительность ИХ и строго линейную ФЧХ. При проектировании таких фильтров отсутствует проблема обеспечения устойчивости.
Проектирование нерекурсивного фильтра сводится к синтезу N отсчетовегоИХ wk , k [0, N 1] , котораяобеспечиваеттребованиякАЧХ
46
фильтра. По свойствам ИХ можно выделить четыре вида нерекурсивных фильтров.
Фильтры вида 1 |
|
|
ИХ фильтра вида 1 |
содержит нечетное количество |
отсчетов |
N 2 L 1 , L (N 1) / 2 |
и является симметричной относительно от- |
|
счета L , т. е. w(2L k) wk , k [0, L 1] , wL 0 . На рис. 2.7 |
приведен |
|
пример ИХ фильтра вида 1 при N 11 . |
|
|
Рис. 2.7
Передаточная функция таких фильтров имеет вид
W (z) z L |
|
L 1 |
z(L k ) z (L k) |
|
, |
(2.6) |
|
wL |
wk |
|
|||||
|
|
k 0 |
|
|
|
|
|
АЧХ и ФЧХ: |
|
|
|
|
|
|
|
|
L |
|
|
|
|
|
|
A( ) | wL 2w(L k ) cos |
(k t ) | , |
|
|
(2.7) |
|||
|
k 1 |
|
|
|
|
|
|
|
( ) L t . |
|
|
|
(2.8) |
||
Фильтрывида1 могутиспользоватьсяприсинтезеНЧФ, ВЧФ, РФ, ПФ.
Фильтры вида 2
ИХ фильтра вида 2 содержит нечетное число отсчетов N 2L 1 , L (N 1) / 2 и является антисимметричной относительно отсчета L ,
47
т. е. w(2L k ) wk , k [0, L 1] , |
wL 0 . На рис. 2.8 представлен при- |
||
мер ИХ фильтра вида 2 при N 11 . |
|
|
|
Передаточная функция таких фильтров имеет вид |
|
|
|
L 1 |
|
, |
(2.9) |
W (z) z L wk z(L k) z (L k) |
|||
k 0 |
|
|
|
Рис. 2.8
АЧХ и ФЧХ:
|
|
|
|
|
|
|
|
L |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A( ) | 2w(L k ) sin (k t ) | , |
|
|
|
|
|
(2.10) |
|||||||||||||
|
|
|
|
|
|
|
|
k 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( ) |
|
L t . |
|
|
|
|
|
|
|
(2.11) |
|||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Фильтры вида 2 могут использоваться только для синтеза ДФ. |
||||||||||||||||||||||
Фильтры вида 3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
ИХ фильтра вида |
3 содержит четное число отсчетов |
N 2L0 , |
||||||||||||||||||||
L0 N / 2 |
и является |
симметричной |
относительно |
точки |
L0 , т. е. |
|||||||||||||||||
w(2L 1 k ) |
wk , k [0, L0 1] . На рис. 2.9 представлен пример ИХ филь- |
|||||||||||||||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
тра вида 3 при N 10 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Передаточная функция таких фильтров имеет вид |
|
|
|
|
||||||||||||||||||
|
|
|
|
|
1 |
|
L 1 |
|
|
|
|
1 |
|
|
|
|
1 |
|
|
|
||
|
|
L |
2 |
|
|
0 |
|
|
L |
2 |
k |
|
L |
2 |
k |
|
|
|
||||
|
W (z) z |
|
0 |
|
|
|
|
0 |
|
|
z |
0 |
|
|
(2.12) |
|||||||
|
|
|
|
|
|
|
wk z |
|
|
|
|
|
|
|
|
|
, |
|||||
|
|
|
|
|
|
|
k 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
48
|
|
|
Рис. 2.9 |
|
|
|
|
|
|
|
|||
АЧХ и ФЧХ: |
|
|
|
|
|
|
|
|
|
|
|
||
A( ) |
|
L0 |
|
|
|
|
|
1 |
|
|
|
, |
(2.13) |
|
|
|
|
|
|
||||||||
|
2 w(L k ) cos k |
2 |
|
t |
|
||||||||
|
|
k 1 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
t . |
|
|
|
|
(2.14) |
|
|
|
( ) L |
|
|
|
|
|
||||||
|
|
|
|
0 |
|
2 |
|
|
|
|
|
|
|
Фильтрывида 3 |
могут использоваться только при синтезе НЧФ, ВЧФ. |
Фильтры вида 4 |
|
ИХ фильтра вида 4 содержит четное число отсчетов N и является |
|
антисимметричной |
относительно отсчета L0 , т. е. w(2 L 1 k ) wk , |
|
0 |
k 0 , L0 1 .
На рис. 2.10 представлен пример ИХ фильтра вида 4 при N 10 .
Рис. 2.10
49
Передаточная функция таких фильтров имеет вид
|
|
|
|
1 |
L 1 |
|
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
||
|
L |
|
|
0 |
|
L |
2 |
k |
|
L |
2 |
k |
|
|
|
|
||||||
W (z) z |
|
0 |
|
2 |
wk |
|
0 |
|
|
z |
|
0 |
|
|
, |
(2.15) |
||||||
|
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
k 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
АЧХ и ФЧХ:
A( ) |
|
L0 |
|
|
|
|
|
|
1 |
|
|
|
, |
|
|
|
|
|
|
|
|
||||||||
|
2 w(L |
k ) sin k |
2 |
|
t |
|
||||||||
|
|
k 1 |
0 |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
( ) |
|
|
|
|
1 |
|
t . |
|
|
|
||||
2 |
L |
2 |
|
|
|
|
||||||||
|
|
|
|
0 |
|
|
|
|
|
|
|
|
||
(2.16)
(2.17)
Фильтры вида 4 используются только при синтезе ДФ. Проектирование фильтров может быть осуществлено тремя спосо-
бами:
1)на основе аппроксимации АЧХ фильтра гармоническим рядом Фурье;
2)на основе частотной выборки;
3)на основе оптимальной аппроксимации АЧХ фильтра.
2.3.2. ФИЛЬТРЫ НА ОСНОВЕ ГАРМОНИЧЕСКОГО РЯДА ФУРЬЕ
Этот метод применим для проектирования фильтров вида 1 и 2. Целью проектирования фильтра является получение его математической модели на основе требований к АЧХ. В данном случае необходимо синтезировать такую W (z) , которая полностью определялась
бы N отсчетами ИХ.
Передаточные функции W (z) фильтров вида 1 и 2 имеют вид
(2.6), (2.9).
Переходя к АЧХ, полагаем z e j t . После несложных преобразований получаем АЧХ фильтра вида 1
L
A(ω) W (ω) | | a0 ak cos(k tω) (2.18)
k 1
50
