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

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

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

2.4. РЕКУРСИВНЫЕ ФИЛЬТРЫ

2.4.1. ОБЩИЕ ПОЛОЖЕНИЯ

Проектирование рекурсивного фильтра сводится к синтезу его передаточной функции W (z) , при которой АЧХ фильтра удовлетворяет за-

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

прерывного фильтра, а синтез W (s) – через передаточную функцию W0 (s0 ) ( s0 j – нормированная частота) непрерывного низкочастотного фильтра-прототипа (НЧФ-П). Идеализированная и идеальная W0ид( j ) АЧХ непрерывного НЧФ-П показаны на рис. 2.15.

Рис. 2.15

Такимобразом, непрерывныйНЧФ-П, ккоторому предъявляютсяте же требования, что и к синтезируемому фильтру, один, а непрерывные фильтры могут быть различные.

На этапе перехода от передаточной функции W (s) непрерывного фильтра к передаточной функции W (z) требуемого рекурсивного

фильтра возникает ряд особенностей, связанных с тем, каким способом этот переход выполнен. Это можно осуществить, используя стандартное Z-преобразование или билинейное Z-преобразование.

71

2.4.2. СТАНДАРТНОЕ Z-ПРЕОБРАЗОВАНИЕ

 

Стандартное Z-преобразование имеет вид

 

 

z e ts ,

 

(2.75)

если s j , получаем

 

 

z e t e j t e t [cos( t) j sin ( t)]

 

e t [cos( t i 2 ) j sin ( t i2 )]

i 0, 1, 2, ...

(2.76)

Следовательно, левая полуплоскость и мнимая ось комплексной плоскости s преобразуются в круг и окружность единичного радиуса комплексной плоскости z (рис. 2.16), а из-за периодичности тригонометрических функций в (2.76) стандартное преобразование допускает неоднозначность перехода к плоскости z и порождает периодичность (с периодом 2 в, в t ) частотных характеристик рекурсивного

фильтра. В результате переход нужно выполнять при i 0 , чтобы не выйти за границу нужного диапазона частот:

 

 

 

 

.

(2.77)

 

m

в

 

t

 

Рис. 2.16

Если выйти за границы этого диапазона, будут происходить искажения частотных характеристик за счет отражения и наложения боковых полос ( i 1, 2, ... ). Обеспечить основной диапазон частот ( i 0 )

можно за счет соответствующего выбора шага дискретизации t .

72

Таким образом, те особенности, которые возникают при использовании стандартного Z-преобразования, позволяют выполнить переход для любого типа фильтра в следующем виде:

1)синтез передаточной функции W0 (s0 ) непрерывного НЧФ-П по заданным параметрам: , , ;

2)представление передаточной функции W0 (s0 ) в виде суммы про-

стых дробей;

3) переход по стандартному Z-преобразованию от простых дробей передаточной функции W0 (s0 ) непрерывного НЧФ-П к дискретным

с последующим восстановлением W 0 (z0 ) – дискретной передаточной функции НЧФ-П;

4) переходот W 0 (z0 ) к W (z) сучетомтребуемых ф , ф, kф , t .

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

1)синтез передаточной функции W0 (s0 ) непрерывного НЧФ-П по задаваемым параметрам: , , ;

2)представление передаточной функции W0 (s0 ) в виде суммы про-

стых дробей;

3)переход от передаточной функции W0 (s0 ) согласно заданнымф , ф, kф к непрерывному требуемому фильтру W (s) ;

4)переход от W (s) к передаточной функции W (z) с использованием стандартного Z-преобразования.

73

2.4.3. БИЛИНЕЙНОЕ Z-ПРЕОБРАЗОВАНИЕ

Билинейное Z-преобразование имеет вид

 

 

 

 

 

 

 

z

 

2

s t

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

s t

 

 

 

 

 

 

 

 

 

 

и если положить s j , то будем иметь

 

 

 

 

 

 

 

 

 

z

 

4 2 t2 2 t2

j

 

 

 

4 t

 

 

 

 

 

 

 

 

(2 t)2 2 t2

(2

t)2 2

t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8 t

 

 

 

 

 

 

 

 

4 t

 

 

 

 

 

1

 

 

 

 

 

 

 

 

exp jarctg

 

 

 

 

 

 

 

.

(2

t)

2

2

2

 

4

2

t

2

 

2

2

 

 

 

t

 

 

 

 

 

 

 

 

t

 

 

Рис. 2.17

(2.78)

(2.79)

Таким образом при билинейном Z-преобразовании (рис. 2.17) происходит трансформация комплексной полуплоскости s в круг единичного радиуса плоскости z . Однако в этом случае переход однозначный

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

отлюбого непрерывного фильтра, асампереходотнепрерывного фильтра (передаточная функция W (s) ) к дискретному фильтру (передаточ-

ная функция W (z ) ) осуществить при замене s следующим образом:

 

2 z 1

 

s

 

 

 

.

(2.80)

t

 

z

1

 

 

 

 

 

 

 

74

Переходя в частотную область

s j и z

e

j

t

(2.81)

 

 

 

 

 

 

 

и подставляя (2.81) в (2.80), после элементарных преобразований получаем соотношение

 

2

tg

 

 

t

,

(2.82)

 

 

 

 

t

 

 

 

2

 

 

устанавливающее связь между частотами непрерывного ( ) и рекурсивного ( ) фильтра. Из выражения (2.82) следует, что полуось[0, ) трансформируется в конечный отрезок [0, в / t] .

В соответствии с (2.82) связь между и нелинейная, поэтому

после перехода характерные частоты АЧХ рекурсивного отличаются от задаваемых частот АЧХ непрерывного фильтра. Например, если АЧХ

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

достаток билинейного Z-преобразования, который может быть устранен, еслинаэтапепроектированиянепрерывногофильтравместочастот

ф , ф использовать их скорректированные согласно (2.82) значения

 

2

 

 

t

 

2

 

 

t

 

ф

t

tg

ф

2

 

, ф

t

tg

ф

2

.

(2.83)

 

 

 

 

 

 

 

 

 

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

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

мент z, в результате чего получим передаточную функцию W (z)

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

75

Таким образом, проектирование рекурсивного фильтра с помощью билинейного Z-преобразования предполагает выполнение следующих операций:

1)синтез передаточной функции W0 (s0 ) непрерывного НЧФ-П по задаваемым параметрам , , ;

2)переход от передаточной функции W0 (s0 ) непрерывного НЧФ-П

кпередаточной функции W (s) требуемого непрерывного фильтра на

основании исходных данных ф , ф, kф , t ;

3)переход от передаточной функции W (s) непрерывного фильтра

кпередаточной функции W (z) требуемого рекурсивного фильтра со-

гласно (2.82).

Если сравнивать билинейное и стандартное Z-преобразование, то предпочтение следует отдать билинейному Z-преобразованию, и в дальнейшем именно это преобразование будет использоваться при проектировании рекурсивного фильтра.

2.4.4.СИНТЕЗ ПЕРЕДАТОЧНОЙ ФУНКЦИИ НЕПРЕРЫВНОГО Í×Ô-Ï

Воснове синтеза передаточной функции W0 (s0 ) лежит аппрокси-

мация квадрата идеальной АЧХ

 

W

( j )

 

2 непрерывного НЧФ-П

 

 

(рис. 2.15) в виде

 

 

0ид

 

 

 

 

 

 

 

 

 

 

 

|W ( j ) |2

 

1

 

.

(2.84)

 

 

 

0

 

1 F

2 ( )

 

 

 

 

 

 

 

 

Функция F( ) задает вид аппроксимации и соответственно вид фильтра, и тем самым его порядок n0 передаточной функции W0 (s0 )

и обеспечение таких исходных данных, как неравномерности и АЧХ в полосах пропускания и задерживания, ширину ( 1) переход-

ной полосы АЧХ. Поскольку порядок n0 определяет сложность реали-

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

же , и .

76

Минимальный порядок аппроксимации при одной и той же погрешности обеспечивают равноволновые функции, представителями которых являются, например, полиномы Чебышёва, ортогональные на интервале [ 1,1] . В этом смысле представляется целесообразным, с точки

зрения минимума параметра n0 , рассматривать фильтры Чебышёва на

основе чебышёвской аппроксимации и фильтры Кауэра, использующие эллиптическую аппроксимацию. Вызывают определенный интерес фильтры Баттерворта, которые обеспечивают минимальную неравномерность в полосе пропускания, и фильтры Бесселя, в наибольшей мере выполняющие требование линейности ФЧХ в полосе пропускания АЧХ.

Аппроксимирующие функции для указанных фильтров имеют вид: 1) фильтры Чебышёва

 

 

 

 

F ( ) Tn

( ) ,

 

 

(2.85)

 

 

 

 

 

 

0

 

 

 

 

где Tn

( ) – полиномы Чебышёва первого рода;

 

 

 

0

 

 

 

 

 

 

 

 

 

 

2) фильтры Кауэра

 

 

 

 

 

 

 

 

 

 

 

K(q)

 

1

 

,

(2.86)

 

 

 

 

F( ) sn n0

 

sn

[ ; g] c; q

 

 

 

 

K(g)

 

 

 

 

 

 

 

 

 

 

где q

 

2

 

, g 1 / 2 ,

 

 

 

 

 

 

1 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/2

[1 xsin2 ] 1/2 d – полныйэллиптическийинтегралпер-

K(x)

 

0

 

 

 

 

 

 

 

 

вого рода, sn[x, y] , sn 1[ ; g]

– эллиптические синус и обратный си-

нус, а условие равноволновости функций (2.86) записывается в форме

 

K (q)

 

K (q )

 

0,

n

нечетное,

 

 

 

 

 

0

 

 

n

 

 

1

, c

 

 

(2.87)

 

K (g )

n

четное;

0 K (g)

 

 

K(q),

 

 

 

 

1

 

 

0

 

 

 

 

 

q1 1 q, g1 1 g .

 

 

В случае выбора параметра

в виде

 

 

 

 

 

 

(2 ) / (1 )2

 

(2.88)

77

обеспечивается соотношение

 

W ( j1)

 

2 (1 )2

,

гарантирующее за-

 

 

даваемую неравномерность

 

0

 

 

 

 

 

АЧХ W0 ( j ) в полосе пропускания

[0, 1];

 

 

 

 

 

 

3) фильтры Баттерворта

 

 

 

 

 

 

F ( ) n0 .

 

(2.89)

Для фильтров Баттерворта и Чебышёва порядок n0 передаточной функции W0 (s0 ) определяется из условия обеспечения требуемой неравномерности АЧХ W0 ( j ) в полосе задерживания [ , ) , которое согласно (2.85) имеет вид

W0 ( j ) 2 1 / [1 F 2 ( )] 2.

4.Фильтры Бесселя обеспечивают максимальную близость ФЧХ

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

Алгоритм синтеза непрерывного НЧФ-П следующий.

1.Определяется порядок n0 передаточной функции W0 (s0 ) в виде

 

 

 

 

 

 

 

 

 

INT

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n0

 

 

 

INT

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

INT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

1

;

q 1 q ;

 

0

q

 

1

 

 

 

 

 

 

 

 

дляфильтровБаттерворта;

ln q0

ln

 

 

 

 

 

 

ln q

 

q2

1

 

 

0

 

0

 

 

 

 

 

 

 

 

 

дляфильтровЧебышёва; (2.90)

ln v

 

v2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K (q1)K(g)

 

дляфильтровКауэра,

K (q)K(g )

 

 

 

1

 

 

 

 

g1 1 g ; INT – целая часть числа.

78

2. Определяютсяполюсаи нули W0 (s0 ) , длячеговычисляетсяколичество пар комплексно-сопряженных полюсов:

l

INT

n0

 

;

(2.91)

0

 

 

2

 

 

 

 

 

 

 

 

 

 

 

P0(2i 1) 0i

j 0i ,

 

P0(2i) 0i

 

j 0i ,

i [1,

l0 ] ,

 

 

(2.92)

при n0 2l0 1 дополнительный вещественный полюс

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P00 00 ,

 

 

 

 

 

 

 

 

 

 

(2.93)

а для фильтра Кауэра также l0 пар мнимых нулей

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z0(2i 1) j 0i ,

 

z0(2i) j 0i ,

i [1, l0 ] .

 

 

 

 

 

(2.94)

 

Полюса имеют следующие параметры:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1) для фильтра Баттерворта

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0i

d0 cos(i0 ) , 0i

d0 sin (i0 ) , 00 d0 ,

 

 

 

(2.95)

причем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

2i 1

при n

2l

, i

 

i

 

при

n 2l

 

1 , d

0

n0

 

0

,

 

0

 

1

;

 

 

 

 

 

0

 

2n0

 

0

0

0

 

n0

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) для фильтра Чебышёва

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0i

 

2

d

 

 

 

 

 

 

sin(i0

),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0i

 

2

d

 

 

 

 

 

cos(i0

),

 

 

 

 

 

 

 

(2.96)

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

00

 

 

d

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

причем d n0

0

 

2 1,

i

 

2i 1

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

2n0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

79

1) для фильтра Кауэра

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0i

1 g sn2 i0 K(g); g cn i0 K(g); g sn v0 ; g1

cn v0

; g1

 

,

(2.97)

 

 

 

 

 

2

 

 

 

 

2

[v0 ; g1

]

 

 

 

 

 

 

 

 

1 1 g sn

 

i0 K(g); g sn

 

 

 

 

 

 

 

 

 

 

sn i0 K (g); g

1 g1sn2 (v0 ; g1)

 

 

 

 

 

 

 

 

sn(v ; g )

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

0

1

 

,

0i

1

1

g sn2 i K

(g); g sn2 v

0

; g

 

 

 

00

 

cn(v0 ; g1)

 

 

 

 

 

0

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0i

v

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sn i0 K(g); g

 

 

 

 

 

 

 

 

 

 

 

 

причем

 

 

2i 1

при

n

2l

,

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

i0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2i

 

 

 

 

 

 

 

 

 

 

 

 

i

 

при

n

2l

1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

n0

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

1 K (g)

 

 

 

2

 

 

cn[v; a]

1 sn

 

[v; a] ,

ν0

 

 

K (q)

ln

0

 

0

1

.

 

n0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, передаточная функция W0 (s0 ) непрерывного НЧФ-П

полностью определяется в соответствии c выражениями (2.90)–(2.97) и в общем случае может быть представлена в виде произведения элементарных сомножителей:

 

 

l0

 

 

 

W0 (s0 ) 0V00 (s0 ) V0i (s0 ) ,

(2.98)

 

 

i 1

 

где

 

 

 

 

1

при n0 2l0 1

 

0

 

длявсех типовфильтров и n0 2l0 для Баттерворта;

 

 

 

при n0 2l0 дляфильтров Чебышёва иКауэра.

 

 

1

 

80

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