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

LR_TsOS_2_BIKh

.pdf
Скачиваний:
0
Добавлен:
22.04.2024
Размер:
2.74 Mб
Скачать

Федеральное агентство связи ордена Трудового Красного Знамени

федеральное государственное бюджетное образовательное учреждение высшего образования

«Московский технический университет связи и информатики»

Кафедра радиотехнических систем

Практикум по дисциплине

ЦИФРОВАЯ ОБРАБОТКА СИГНАЛОВ

Лабораторная работа № 3

СИНТЕЗ ЦИФРОВЫХ БИХ-ФИЛЬТРОВ МЕТОДОМ БИЛИНЕЙНОГО Z -ПРЕОБРАЗОВАНИЯ

Москва, 2016 г.

УДК 621.391:519.27

План подготовки УМД 2015/2016 уч. года

Практикум по дисциплине

«ЦИФРОВАЯ ОБРАБОТКА СИГНАЛОВ»

СИНТЕЗ ЦИФРОВЫХ БИХ-ФИЛЬТРОВ МЕТОДОМ БИЛИНЕЙНОГО Z - ПРЕОБРАЗОВАНИЯ

В лабораторной работе № 3 изучается метод билинейного z - преобразования на примере синтеза цифровых ФНЧ и ФВЧ по аналоговым фильтрам прототипам Чебышёва и Баттерворта. Провести моделирование синтезированного цифрового фильтра в среде «Спектр-2».

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

Список лит. 11 назв., ил. 9, табл. 5.

Составители: Лобов Е.М., Смердова Е.О. Рецензент

Издание утверждено советом факультета Радио и Телевидения. Протокол № … от

_._.2016 г.

2

ЛАБОРАТОРНАЯ РАБОТА № 3

СИНТЕЗ ЦИФРОВЫХ БИХ-ФИЛЬТРОВ МЕТОДОМ БИЛИНЕЙНОГО Z -ПРЕОБРАЗОВАНИЯ

Цель работы:

Изучение метода билинейного z - преобразования на примере синтеза цифровых ФНЧ и ФВЧ по аналоговым фильтрам прототипам Чебышёва и Баттерворта. Моделирование синтезированного цифрового фильтра проводится в среде «Спектр-2».

Номер варианта выбирается студентом из следующей таблицы:

Условные

 

Список требований

Задаваемые значения

обозначения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

Частота

F

10000 200N

 

 

 

 

 

 

 

 

 

бр

 

 

 

 

 

д

 

 

 

дискретизации

д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F 1200 20N

бр

– Баттерворт

 

 

1

 

 

 

Граничная частота ПП

1

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

– Чебышев

 

 

 

 

 

 

 

F 1500 30N

бр

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

бр

 

 

 

 

 

 

 

 

Fз

3120 50Nбр

10

 

– Баттерворт

 

 

 

 

 

 

 

5

 

F

 

 

 

Граничная частота ПЗ

 

 

 

 

 

 

 

 

 

з

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fз

2800 40Nбр

 

– Чебышев

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Максимально

 

 

 

 

 

 

 

 

 

 

a

max

 

 

 

допустимое затухание

a

 

0.4455

 

 

 

 

 

 

 

 

 

 

 

max

 

 

 

 

 

 

 

 

 

 

 

 

 

в ПП (дБ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Минимально

 

 

 

 

 

 

 

 

 

 

amin

 

 

 

допустимое затухание

amin 40

 

 

 

 

 

 

 

 

 

 

 

 

в ПЗ (дБ)

 

 

 

 

 

 

 

 

 

Nb

1, 2,...40

,

Nb

- результат целочисленного деления номера бригады на 5 (с

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

отбрасыванием цифр после запятой у результата деления).

Тип функции фильтрации определяется по чётности/нечётности последней цифры в номере студенческого билета. Если последняя цифра номера нечётная, то требуется провести синтез фильтра Чебышёва, если чётная, то рассчитывается фильтр Баттерворта.

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

3

нечётная, то рассчитывается фильтр нижних частот, чётная – фильтр верхних частот.

1 Домашний расчёт.

1.1 Домашний расчёт состоит из следующих пунктов:

1.Формулировка требований к ЦФ.

2.Синтез передаточной функции АФП, включающий в себя: 2.1.оценку порядка ПФ АФП;

2.2.расчет полюсов и нулей, а также коэффициента усиления передаточной функции АФП;

2.3.запись аналитического выражения передаточной функции АФП.

3.Переход от передаточной функции АФП к передаточной функции ЦФ с помощью билинейного z -преобразования.

4.Запись аналитического выражения передаточной функции ЦФ в виде произведения передаточных функций биквадратных звеньев.

1.2 Пример расчёта домашнего задания для ФНЧ Чебышёва.

Спроектируем ФНЧ Чебышева I рода. Сформулируем требования к фильтру, - для ФНЧ это:

1.значения граничных частот полосы задержания и полосы пропускания ( Fз и F1 соответственно);

2.значение коэффициента максимального ослабления в полосе пропускания aр max , дБ;

3.значение коэффициента минимального ослабления в полосе

задержания aр min

, дБ;

4. частота дискретизации сигналов, обрабатываемых фильтром,

Fд

.

Примем значения граничных частот и коэффициентов

a

р max

 

,

a

р min

 

низкочастотного аналогового фильтра прототипа равными:

F 1900Гц, F 4940Гц,

 

1

з

 

 

a

p max

0.4455дБ, a

p min

40дБ,

 

 

 

Fд 17000 Гц.

Синтез АФП осуществляется в области нормированных частот. При этом граничная угловая частота полосы пропускания принимается равной 1 рад/с. Проведём соответствующее нормирование оси частот. Тогда

4

нормированные граничные частоты полосы пропускания и полосы задержания будут равны соответственно

 

2 F

1,

 

 

2 F

2.6.

 

1

 

з

 

 

 

 

 

 

 

1

2 F

 

з

 

2 F

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

Следующим шагом синтеза является определение порядка фильтра и

коэффициента неравномерности ослабления в

полосе пропускания

.

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

Н (j )

2

 

 

 

 

 

1

1

2

ch

2

nArch

 

 

,

(1)

где 2 F модуля квадрата, ch(x) - гиперболический косинус, Arch(x) - гиперболический арккосинус.

С учетом (1) квадрат модуля передаточной функции АФП можно записать в форме:

Н ( p)

 

2

H ( p)H ( p)

2

 

2

1

 

 

 

 

 

nArch( jp)

 

 

 

1

ch

 

 

 

 

 

,

(2)

где

H ( p)

- операторная передаточная функция АФП,

T ( ) ch n Arch

- полином Чебышёва, являющийся

p j .

функцией

фильтрации данного фильтра.

Коэффициент неравномерности рабочего ослабления найдём через коэффициент максимального подавления и значение граничной частоты полосы пропускания:

A

 

 

 

 

 

2

2

( )) a

 

,

10lg(1

T

 

p max

p

 

 

 

 

 

 

 

 

 

1

 

T ( ) T(1) 1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

10lg(1

2

) a

 

 

 

,

 

 

 

 

p max

 

 

 

 

 

 

 

 

 

 

 

 

 

10

0.1a

p max

1.

 

 

 

 

 

Подставив численное значение

 

aр max , получаем 0.32867 .

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

5

 

 

 

10

0.1a

p min

1

 

 

 

 

 

 

nArch

3

Arch

 

 

2

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

0.1ap min

1

 

 

 

Arch

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

.

Arch

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

Порядок фильтра может быть только целым числом,

A

2

2

(

)) a

 

10lg(1 T

 

p min

p

 

 

з

 

тогда справедливо:

,

 

2ch2

 

nArch

3

 

 

 

10lg 1

 

 

 

 

a

p min

 

,

 

 

 

 

 

10

0.1a

p min

1

ch nArch

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

0.1ap min

 

1

 

 

 

Arch

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Arch

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

(3)

где x – операция округления в сторону ближайшего целого значения,

превышающего x .

Проведя арифметические операции с учетом заданных требований к фильтру, получим n 4 . Теперь мы имеем все данные для нахождения корней знаменателя квадрата модуля АЧХ.

Корни находятся по формуле:

pk

где

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

Arsh

 

 

2k 1

 

 

Arsh

 

 

2k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sh

 

sin

 

 

 

jch

 

cos

 

 

 

, k 1,2,...,2n.

 

 

 

 

 

 

 

 

 

n

 

 

 

2n

 

 

n

 

 

 

2n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sh(x)

- гиперболический синус.

 

 

 

 

 

 

 

 

 

Из рассчитанных корней выбираем отвечающие соображениям устойчивости (система устойчива, если полюсы её передаточной функции лежат в левой полуплоскости её p -плоскости, т.е. имеют отрицательную вещественную часть).

С учетом изложенного, отобранные корни равны:

6

p -0.43804-j0.42352,

1

 

p

-0.18144-j1.0225,

2

 

p

-0.18144+j1.0225,

3

 

p

-0.43804+j0.42352.

4

 

(4)

В результате расчетов получены 2 пары комплексно-сопряжённых корней, которые мы относим к H ( p) . Последним шагом в формировании низкочастотного АФП Чебышёва является запись аналитического выражения его операторной передаточной функции, которая имеет вид:

Н ( p)

1

 

n 1

 

2

 

 

где

 

1

 

– коэффициент усиления.

2

n 1

 

 

 

 

 

 

В нашем случае H ( p) принимает вид:

1

n ( p

k 1

,

p

k

)

 

 

Н ( p)

0.38032

 

 

 

 

 

( p p )( p p )( p p )( p p )

1

2

3

4

где

1

0.38032 .

3

 

2

 

,

(5)

Дальнейшие действия синтеза аналоговых фильтров рассматриваться не будут, так как в них нет необходимости при получении передаточной функции цифрового БИХ-фильтра.

Для перехода от передаточной функции АФП к передаточной функции ЦФ используется билинейное z - преобразование, которое заключается в подстановке в передаточную функцию выражения

где

Fд

 

1 z

1

p 2F

 

 

 

1

д

1

z

 

 

– частота дискретизации сигналов,

,

обрабатываемых фильтром,

(6)

z

комплексный параметр Z-преобразования. Перед использованием билинейного

передаточную функцию Н ( p) (5) в удобном

Z-преобразования, преобразуем для этого виде:

Н ( p)

 

 

0.38032

 

 

 

 

 

 

 

 

( p2 p( p p ) p p )( p2

p( p p ) p p )

 

 

1

4

1

4

2

3

2

3

(7)

 

 

 

 

0.38032

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

( p2 0.87608 p 0.37124)( p2 0.36288 p 1.0784)

 

В (7) проведено перемножение множителей, содержащих комплексносопряжённые корни полинома знаменателя ПФ (5) (см. значение корней (4)).

7

Теперь можно проводить билинейное

z

-

предварительно нормировав частоту дискретизации (примем

 

 

F

 

 

Fs

 

д

1.424

,

2 F

 

 

 

 

 

 

1

 

 

тогда замена примет вид

преобразование,

Fд

17000Гц )

p 2F 1 z

s 1 z

В результате подстановки в

Н

1

1

z 1

 

2.848

 

 

1

1

z 1

( p) выражения

,

(8)

(8) получаем:

Н (z)

 

 

 

 

 

 

 

 

 

 

 

 

 

0.38032

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1 z

 

 

 

 

 

1 z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.848

 

 

 

 

0.87608

2.848

 

 

 

0.37124

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

1 z

 

 

 

 

 

1 z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1 z

 

 

 

 

 

 

1 z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.848

 

 

0.36288

2.848

 

1.0784

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1 z

 

 

 

 

 

1 z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Умножим числитель и знаменатель каждой дроби на

 

 

1 z 1

2

 

и получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.38032 1 z 1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Н (z)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.848 1 z 1

 

 

 

0.87608 2.848 1

z 1

 

 

 

1 z 1

 

 

0.37124 1 z 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 z 1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2.848 1 z 1

 

 

0.36288 2.848

1 z 1

 

 

1 z 1

 

1.0784

1 z 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Раскроем скобки

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2z 1 z 2

 

 

 

 

 

 

 

 

 

 

 

 

 

Н (z)

 

 

 

 

 

 

 

 

 

 

 

0.38032 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z 2

 

 

 

 

z 2

 

 

 

 

 

 

 

 

2z 1

z 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8.1111 1 2z 1

 

 

2.4951 1

 

 

0.37124 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2z 1

z 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.0335

 

 

 

 

 

 

 

 

 

z 2

 

 

 

 

 

 

 

 

 

 

 

 

 

8.1111 1 2z 1 z 2

 

1 z 2

 

 

 

1.0784 1 2z 1

 

 

 

 

 

 

 

и сложим слагаемые с одинаковыми степенями при z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.38032 1

 

2z 1 z 2

 

 

 

 

 

 

 

 

 

 

 

 

1 2z 1

z 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

15.4797z 1 5.9872z 2

 

 

 

 

 

 

 

 

 

 

 

 

8.1560z 2

 

 

Н (z)

10.9774

 

 

 

10.2230 14.0654z 1

 

.

Поделим числитель и знаменатель первой дроби на 10.9774 , а второй дроби - на 10.2230 , для того чтобы свободные члены в знаменателе стали равны единице, и получим

8

 

 

 

2z 1

z 2

 

 

 

 

2z 1 z 2

 

 

 

Н (z)

0.03464 1

 

0.09781 1

 

 

 

.

 

1.4101z 1

0.5454z 2

 

1.3759z 1

 

0.7978z 2

 

 

 

 

1.0

1.0

 

 

Тем самым

мы

Н (z) 0.0033888 * H

(z)H

(z)

1

2

 

привели

наше

выражение

к

виду

и получили

два

звена второго

порядка, где

0.0033888 0.03464 0.09781 – общий коэффициент усиления. Выпишем отдельно выражения для каждого звена:

 

 

 

 

1 2z

1

z

2

 

 

b

b z

1

b z

2

H

(z)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1 1.4101z

1

0.5454z

2

 

1 a z

1

a z

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

 

 

 

 

1 2z 1

z 2

 

 

b b z 1

b z 2

H

 

(z)

 

 

 

 

 

 

 

 

 

0

1

 

 

 

2

 

 

2

 

1.3759z 1

0.7978z 2

1 a z 1

 

 

z 2

 

1

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

(9)

(10)

Выпишем отдельно коэффициенты каждого звена в таблицу, без учета общего коэффициента усиления фильтра.

Таблица 1. Коэффициенты звеньев ФНЧ Чебышева

Звено №

b

b

b

a

a

a

 

0

1

2

0

1

2

1

1

2

1

1

-1.4101

0.54542

2

1

2

1

1

-1.3759

0.79781

Определим искажение частот, получившееся в

билинейного z

-преобразования. Частоты АЧХ АФП

соотношениями:

 

ходе проведения и ЦФ связаны

где

ˆ

/

Fд

ˆ 2arctg

 

2F

 

 

д

F

F

д

 

 

 

- нормированная частота

 

2arctg

F

F

 

 

 

 

 

 

д

tg

ˆ

,

 

2

 

 

 

 

АЧХ цифрового фильтра,

2

(11)

(12) f

угловая частота АЧХ цифрового фильтра.

Подставляя в выражение для нахождения частоты АЧХ цифрового фильтра, значения граничных частот полосы пропускания и полосы задержания аналогового фильтра, получим граничные частоты ЦФ:

f

 

 

 

Fд

 

arctg

F1

1827.2311 Гц,

 

 

 

 

1

 

 

 

F

 

 

 

 

 

 

 

 

 

д

 

(13)

 

 

 

 

Fд

arctg Fз

f

 

 

4003.8081 Гц.

з

 

 

 

 

 

Fд

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

 

На этом домашний расчёт БИХ-фильтра Чебышёва окончен.

1.3 Пример расчёта домашнего задания для ФНЧ Баттерворта.

Требования к фильтру и, соответственно, значения нормированных частот такие же, как и для ФНЧ Чебышёва (см. пункт 1.1 домашнего расчёта)

Квадрат модуля АЧХ фильтра Баттерворта имеет вид:

Н (j )

2

 

1

,

 

 

 

1 2 2n

где ( ) 2n - функция фильтрации фильтра Баттерворта. Квадрат модуля передаточной функции АФП можно записать виде:

где

H (

p)

Н ( p)

2

H ( p)H ( p)

 

1

 

2

( jp

 

 

 

 

 

 

1

- операторная передаточная функция АФП,

)

2

 

p

n

,

 

 

j

.

Коэффициент неравномерности рабочего ослабления в полосе пропускания и порядок фильтра для фильтра Баттерворта определяется по формулам:

 

 

 

 

 

10

0.1a

p max

1

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

0.1ap min

1

 

 

 

 

 

10

0.1a

p min

1

 

 

lg

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

lg

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

,

 

lg

 

 

 

 

 

 

 

2lg

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где x – операция округления в сторону ближайшего целого значения,

превышающего x .

Проведя несложный арифметический расчёт, учитывающий выдвинутые к фильтру требования, получим значения 0.32867 , n 6 .

Все данные для нахождения корней знаменателя квадрата модуля АЧХ получены.

Корни найдём по формуле для чётных значений порядка фильтра Баттерворта:

p

 

 

 

1

 

cos

 

2k 1

 

 

j sin

 

k

n

 

 

 

2n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из рассчитанных корней, выбираем те, устойчивости:

2k 1

 

 

 

2n

 

 

которые

,

k 1, 2,..., 2 n .

отвечают соображениям

10

Соседние файлы в предмете Цифровая обработка сигналов