Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Введение в теорию сигналов. Цифровая обработка сигналов. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

Другой формой математической модели непрерывного фильтра может служить интегральное уравнение свертки

 

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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]