673_Raschet_analogovykh_i_diskretnykh_ustrojstv_
.pdfРезультаты запишем в таблицу 2.7. Спектр дискретного сигнала является периодической функцией частоты, как показано на рисунке 2.12.
Таб. 2.7. Отсчеты спектра входного сигнала
|
k |
|
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
U |
k |
,В |
2,66 |
2,59 |
2,41 |
2,13 |
1,80 |
1,47 |
1,19 |
1,01 |
1 |
1 |
|
|
|
|
|
|
|
|
|
1 k 1 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
||
|
k |
|
8 |
9 |
10 |
11 |
12 |
13 |
14 |
15 |
U |
k |
,В |
0,94 |
1,01 |
1,19 |
1,47 |
1,80 |
2,13 |
2,41 |
2,59 |
1 |
1 |
|
|
|
|
|
|
|
|
|
1 k 1 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
U1 k 1 ,B
2
1
-2 0 2 4 6 8 10k 12 14 16 18 20 k
Рис. 2.12. Амплитудный спектр дискретного сигнала
Из рисунка видно, что амплитудный спектр дискретного сигнала является периодическим повторением двухстороннего спектра аналогового сигнала с масштабным коэффициентом N , равным периоду.
2.4. Расчет дискретного фильтра
Для выделения заданной гармоники сигнала необходимо рассчитать дискретный полосовой фильтр.
Центральная частота полосы пропускания f0 совпадает с частотой выделяемой гармоники:
f0 l fг , |
(2.30) |
где l – номер выделяемой гармоники. |
А и |
Требования к максимальному ослаблению в полосе пропускания |
|
минимальному ослаблению в полосе непропускания Amin заданы. |
|
Передаточная функция БИХ-фильтра может быть получена путем билинейного преобразования передаточной функции аналогового полосового фильтра.
В качестве аналогового полосового фильтра выберем полиномиальный фильтр Баттерворта. Поскольку гармоники сигнала на выходе нелинейного
21
преобразователя достаточно далеко разнесены по частоте, порядок фильтра может быть получен невысокий. Частоты соседних гармоник должны попадать в полосу непропускания фильтра. Характеристика ослабления фильтра должна обладать геометрической симметрией относительно выделяемой гармоники. Это имеет место при выполнении условия (рис. 2.13).
f0 |
f2 f2 |
|
f3 f3 |
(2.31) |
A, дБ |
|
|
|
|
|
Amin |
|
|
|
|
|
ПН |
|
|
|
|
ПН |
А |
|
ПП |
|
|
|
|
|
|
|
|
|
f |
f |
f |
f |
f |
f ,кГц |
3 |
2 |
0 |
2 |
|
3 |
Рис. 2.13. Требования к аналоговому полосовому фильтру
Расчет полосового фильтра обычно сводят к расчету НЧ-прототипа [1]. Для определения нормированной частоты НЧ-прототипа необходимо за-
дать границы полосы пропускания и непропускания справа f2, f3 . Из соотно-
шения для f0 определяют f2, f3 , а затем определяют нормированную частоту НЧ-прототипа
|
3 |
|
f |
3 |
f3 |
. |
(2.32) |
f |
|
|
|||||
|
|
2 |
f2 |
|
Далее определяют коэффициент неравномерности ослабления в полосе пропускания
|
100,1 А 1 |
|
|
|
и порядок фильтра |
|
|
|
|
|
100,1Amin 1 |
|||
|
lg |
|
|
|
|
2 |
|||
m |
|
|
. |
|
|
|
|
2lg 3
(2.33)
(2.34)
В случае если полученное значение mне удовлетворяет заданным требованиям, необходимо изменить значения f2, f3 и повторить расчет.
Далее для выбранного порядка фильтра m находят полюсы передаточной функции НЧ-прототипа по формулам:
22
р€ |
|
1 |
sin |
2k |
1 |
|
|
|
jcos |
2k 1 |
|
|
, |
k 1,2...2m. |
(2.35) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
||||||||||||||||
к |
|
|
|
|
|
2m |
|
|
|
|
2m |
|
|
|
|
||||||
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
Определим полюсы передаточной функции аналогового полосового фильтра по известным полюсам НЧ-прототипа. Для этого воспользуемся соотношением:
|
|
|
|
|
2 |
|
|
pi, jПФ |
|
j j j |
|
|
j j j |
02 . |
(2.36) |
|
2 |
||||||
2 |
|
|
|
|
|
Здесь 2 f2 f2 – ширина полосы пропускания;
0 2 f0 – центральная частота ПП фильтра;
j j j – j -ый полюс передаточной функции НЧ-прототипа.
Одной паре комплексно-сопряженных полюсов передаточной функции НЧ-прототипа соответствует две пары комплексно-сопряженных полюсов передаточной функции аналогового полосового фильтра. Одному вещественному полюсу НЧ-прототипа соответствует одна пара комплексно-сопряженных полюсов Н р аналогового полосового фильтра.
На втором этапе формируют передаточную функцию аналогового полосового фильтра в виде произведения сомножителей второго порядка:
H p |
a1p |
|
|
a2 p |
|
|
am p |
|
||
|
|
|
|
. |
(2.37) |
|||||
p2 b p b |
p2 b p b |
p2 b |
p b |
|||||||
|
1 |
01 |
2 |
02 |
|
m |
0m |
|
Каждый сомножитель соответствует одной паре комплексно-сопряженных полюсов. Коэффициенты числителя и знаменателя определяются из следующих соотношений:
a a |
|
a |
|
|
|
|
, |
(2.38) |
|
|
m |
|
|
|
|||||
1 |
2 |
|
|
m |
|
где bi 2 i ; b0i i2 i2.
i и i – действительная и мнимая части i-го полюса передаточной функции ПФ.
В результате найдена денормированная передаточная функция аналогового полосового фильтра.
Перейдем от передаточной функции аналогового полосового фильтра к
передаточной функции дискретного полосового фильтра. Переменные |
р и z |
||||
связаны между собой билинейным преобразованием: |
|
||||
p |
2 |
|
1 z 1 |
, |
(2.39) |
|
1 z 1 |
||||
|
T |
|
|
где Т – период дискретизации.
Применение данного преобразования приводит к искажению частотных характеристик аналогового фильтра. Для исключения таких искажений существуют методы, изучаемые в курсе ЦОС.
23
Для уменьшения искажений в наших расчетах воспользуемся формулой:
p |
1 z 1 |
, |
(2.40) |
|
1 z 1 |
||||
|
|
|
где определяется из условия совпадения центральных частот аналогового и дискретного полосовых фильтров:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 f0 ctg |
f0 |
, |
|
|
|
|
|
(2.41) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
f0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
fд |
|
|
|
|
|
|
||
где |
– центральная частота полосы пропускания, |
|
|
|
|
||||||||||||||||||||||||||
|
fд |
– частота дискретизации. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
Для того чтобы получить передаточную функцию полосового БИХ- |
||||||||||||||||||||||||||||||
фильтра, |
подставим формулу 2.40 в выражение 2.37, в результате получим вы- |
||||||||||||||||||||||||||||||
ражение вида: |
|
|
|
|
1 z 2 |
|
|
|
|
|
1 z 2 |
|
|
|
|
|
1 z 2 |
|
|
||||||||||||
|
|
H z |
|
|
|
|
|
|
|
|
|
|
|
, (2.42) |
|||||||||||||||||
|
|
|
|
z 1 |
b |
|
z 2 |
|
|
|
z 1 |
b |
|
z 2 |
|
z 1 b |
z 2 |
||||||||||||||
|
|
|
|
|
|
|
|
|
1 b |
|
|
1 b |
|
|
|
1 b |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
11ц |
|
21ц |
|
|
|
12ц |
|
|
22ц |
|
|
|
|
1mц |
2mц |
|
|
||||
где – постоянный умножитель на входе, |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
a1 a2 am |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
c1 c2 cm |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
с b |
b0i |
|
– вспомогательная переменная, |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
i |
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
bi, |
b0i |
– коэффициенты полинома знаменателя i-го аналогового полосового |
|||||||||||||||||||||||||||||
фильтра |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b b0i |
|
|
|
|
|
|
||||||||||
|
|
|
b |
2 |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
b1iц |
|
|
0i |
|
|
|
|
, |
|
|
b2iц |
|
|
|
|
i |
|
|
|
|
. |
|
|
|
|
||||||
|
|
ci |
|
|
|
|
|
|
|
ci |
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Каждый сомножитель передаточной функции реализуется в виде дискретной канонической рекурсивной схемы второго порядка. Звенья соединяются каскадно.
Расчет фильтра заканчивается построением его полной схемы с указанием значений элементов и расчетом его частотных характеристик. К частотным характеристикам относятся комплексная частотная характеристика фильтра, а также зависимость ослабления от частоты. Расчет той и другой характеристик производится на основе передаточной функции H z путем замены z ej T .
При этом вначале на выбранных частотах рассчитывают H j отдельных звеньев, а затем всего фильтра, используя соотношение:
m |
|
H j Hi j . |
(2.43) |
i 1 |
|
Так как частотная характеристика дискретной цепи – функция периодическая, ее достаточно рассчитать в диапазоне частот 0 д .
24
Для дальнейшего расчета сигнала на выходе фильтра число расчетных точек H j должно совпасть с числом отсчетов дискретного сигнала на пе-
риоде. |
|
A |
|
||
Зависимость ослабления от частоты |
определяют по формуле |
||||
A 20lg |
|
1 |
|
. |
(2.44) |
|
H j |
||||
|
|
|
|
В число расчетных частот необходимо включить граничные частоты полос пропускания и непропускания.
ПРИМЕР РАСЧЕТА:
Требуется рассчитать БИХ-фильтр для выделения второй гармоники при частоте генерируемых колебаний 10 кГц;
неравномерность ослабления в ПП А 1дБ ;
минимально допустимое ослабление в ПН Amin 20дБ (рис. 2.13); порядок НЧ-прототипа равен m 2;
период дискретизации Т 6,25мкс.
Определим частоту выделяемой гармоники: f0 l fг 2 10 103 20 кГц. Выберем границу полосы пропускания f2 21,025 кГц и границу полосы непропускания f3 25 кГц.
|
|
|
|
|
|
|
|
|
|
|
|
f02 |
|
|
f02 |
|
Так как f0 |
f |
|
|
f2 |
f |
|
|
f3 , то |
f2 |
19,025 кГц, |
f3 |
16 кГц . |
||||
2 |
3 |
f2 |
f3 |
Определим требования к НЧ-прототипу:
|
|
|
f |
|
f3 |
|
25 103 |
16 103 |
|
|
|
|
|
|
|
|
|
3 |
|
3 |
|
4,5 |
|
100,1 А 1 |
100,11 |
1 0,509 |
|||||||
f |
2 |
f2 |
21025 |
19025 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
Определим порядок фильтра Баттерворта по формуле 2.34.
|
100,1Amin |
1 |
|||
|
lg |
|
|
|
|
|
2 |
|
|||
Для всех вариантов должно выполняться неравенство: 1 |
|
|
|
2. |
|
|
2lg 3 |
|
|
||
|
|
|
|
|
Если это неравенство не выполняется, необходимо изменить значения f2 и f3 , |
||||||||||||||
но так, чтобы f3 f0 |
fг |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
2 |
|
|
100,1Amin |
1 |
|
100,120 1 |
|
|||||||
|
|
|
|
|
|
|||||||||
|
|
|
|
lg |
|
|
|
|
|
lg |
|
|
|
|
|
|
|
|
|
2 |
|
|
0,509 |
2 |
|
||||
Для рассматриваемого примера |
m |
|
|
|
|
|
|
|
|
1,98. |
||||
|
2lg 3 |
|
|
|
2lg4,5 |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
При округлении в большую сторону до целого числа m 2.
25
Полученный порядок фильтра удовлетворяет требованиям, можно продолжать расчет.
Запишем выражение для квадрата АЧХ передаточной функции:
H j |
|
2 |
1 |
|
1 |
|
. |
(2.45) |
|
|
|||||||
|
1 2 2m |
1 0,5092 |
4 |
|||||
|
|
|
|
|
|
Определим полюсы квадрата АЧХ по формуле 2.35:
рˆк |
|
|
1 |
|
2k 1 |
2k 1 |
|
|
|||||||||||||
|
|
|
|
|
|
|
sin |
|
|
|
jcos |
|
|
|
, |
k 1,2,3,4 и выберем те, которые |
|||||
|
|
|
|
|
|
4 |
|
|
4 |
||||||||||||
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
находятся в левой полуплоскости: |
|
|
|
|
|||||||||||||||||
ˆ |
|
1 |
|
|
|
|
|
|
|
|
j0,7071 |
|
|
|
|
|
|
|
|||
р2 |
|
|
|
|
|
|
|
|
0,7071 |
|
|
0,9913 |
|
j0,9913 , |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||
ˆ |
|
1 |
|
|
|
|
|
|
|
|
j0,7071 |
|
|
|
|
|
|
|
|||
р3 |
|
|
|
|
|
|
|
0,7071 |
|
|
0,9913 |
|
j0,9913 . |
||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
Для отыскания полюсов передаточной функции аналогового полосового фильтра воспользуемся соотношением 2.36
|
|
|
|
|
2 |
|
pi, jПФ |
|
j j j |
|
|
j j j |
02 |
|
2 |
|||||
2 |
|
|
|
|
где 2 f2 f2 2 21025 19025 12566,37рад/ с;
02 2 f0 2 2 20 103 2 1,5791 1010 радс 2 .
Полученные значения полюсов удобно представить в виде таблицы.
Табл. 2.8. Значения полюсов аналогового полосового фильтра
Номер полюса |
Полюсы Н р полосового фильтра |
||
104 |
j 104 |
||
|
|||
1,3 |
-0,5920 |
11,9436 |
|
2,4 |
-0,6537 |
13,1892 |
Передаточная функция ПФ может быть записана в виде произведения двух сомножителей второго порядка
|
|
H p |
|
a1p |
a2 p |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
, |
(2.46) |
||||||||
|
|
|
p2 b p b |
p2 |
b p b |
|||||||||||||
|
|
|
|
|
|
|
|
1 |
01 |
|
|
2 |
02 |
|
|
|||
где |
а1 а2 |
|
|
|
12566,37 |
|
|
17616. |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
m |
2 0,509 |
|
|
|
|
|
|
|
|
|
||||||
Коэффициенты при |
р |
в знаменателях сомножителей |
bi 2 i , а сво- |
бодные члены b0i i2 i2. Их значения сведем в таблицу 2.9.
26
Табл. 2.9. Значения коэффициентов передаточной функции аналогового полосового фильтра
|
|
|
|
|
Номер со- |
|
|
|
|
|
|
|
|
|
|
Значения коэффициентов |
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
множителя |
|
|
|
|
ai |
|
|
|
|
|
|
|
|
|
|
|
|
bi |
|
|
b0i |
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
1,7616 · 104 |
|
1,1839 · 104 |
|
1,4300 · 1010 |
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
1,7616 · 104 |
|
1,3074 · 104 |
|
1,7438 · 1010 |
|
|
|
|
||||||||||||||||||||||||||
Тогда передаточная функция искомого ПФ будет иметь вид |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
H p |
|
|
|
|
|
17616p |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
17616p |
|
|
|
|
. |
|
|
(2.47) |
|||||||||||||||
|
p |
2 |
11839p |
|
|
|
|
|
|
|
|
|
10 |
|
p |
2 |
|
|
|
|
|
|
|
|
|
10 |
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
1,4300 10 |
|
|
|
|
13074p 1,7438 10 |
|
|||||||||||||||||||||||||||||||
|
Для того чтобы перейти к передаточной функции БИХ-фильтра, выпол- |
||||||||||||||||||||||||||||||||||||||||||||||
ним замену в выражении (2.47): |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
20 103 |
|
|
|
|
|
|
||||||||||||||||||||||
|
|
1 z 1 |
|
|
|
|
|
|
|
|
|
|
|
f0 |
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
||||||||||||||
p |
|
|
|
|
|
|
, |
где |
2 f0 |
|
ctg |
|
|
|
|
|
2 20 10 ctg |
|
|
|
|
|
303379. |
|
|||||||||||||||||||||||
1 z 1 |
|
|
fд |
160 103 |
|
|
|
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
Передаточная функция БИХ-фильтра может быть записана в виде произ- |
||||||||||||||||||||||||||||||||||||||||||||||
ведения постоянного сомножителя и двух сомножителей второго порядка: |
|
||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
H z |
|
|
|
|
|
1 z 2 |
|
|
|
|
|
|
|
|
1 z 2 |
|
|
|
|
, |
|
(2.48) |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
1 b |
|
|
z 1 b |
z |
2 |
|
1 b |
z 1 b |
|
z 2 |
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
11ц |
|
|
|
|
21ц |
|
|
|
|
|
|
|
12ц |
|
22ц |
|
|
|
|
|
|
|
|
||||||
где – постоянный умножитель на входе, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
a1 a2 |
, |
|
|
|
|
|
с b |
b0i |
|
– вспомогательная переменная, |
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
c c |
|
|
|
|
|
i |
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
1 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b b0i |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
b |
2 |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
b1iц |
0i |
|
|
|
|
|
, |
|
|
|
|
b2iц |
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
Их |
значения сведем в |
|||||||||||||||||||||||
|
|
|
ci |
|
|
|
|
|
|
|
|
|
|
|
|
ci |
|
|
|||||||||||||||||||||||||||||
таблицу 2.10. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
Табл. 2.10. Значения коэффициентов дискретного фильтра |
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
Номер со- |
|
|
|
|
|
сi |
|
|
|
|
|
|
|
|
|
|
|
|
b1iц |
|
|
|
b2iц |
|
|
|||||||||||||||||
|
|
|
|
|
множителя |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
362353,9 |
|
|
|
|
|
|
|
|
1,4143 |
|
|
|
-0,9347 |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
373933,7 |
|
|
|
|
|
|
|
|
1,3152 |
|
|
|
-0,9301 |
|
|
|
|
|||||||||||||||
Определим значение постоянного умножителя на входе дискретной цепи |
|||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
a a |
|
|
|
|
|
|
|
|
176162 |
|
|
|
|
|
|
2,29 10 3. |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
1 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
c1 c2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
362353,9 373933,7 |
|
|
|
|
|
|
|
|
|
|
|
Передаточная функция искомого БИХ-фильтра с учетом полученного значения коэффициента имеет вид:
H z 2,29 10 3 |
|
1 z 2 |
|
1 z 2 |
. |
|
|
|
0,9301z 2 |
||||
|
1 1,4143z 1 0,9347z 2 |
1 1,3152z 1 |
|
|||
|
27 |
|
|
|
|
Данная схема реализуется каскадным соединением умножителя и двух рекурсивных звеньев второго порядка. Схема БИХ-фильтра приведена на рисунке 2.11.
Комплексную частотную характеристику БИХ-фильтра получим из пе-
редаточной функции H z , выполнив замену |
z ej T , тогда H j запишем: |
|||||||||||||||
H j 2,29 10 3 |
|
|
1 e j2 T |
|
|
|
|
1 e j2 T |
|
|||||||
|
|
0,9347e j |
|
|
1 1,3152e j T |
0,9301e j2 T |
||||||||||
|
|
1 1,4143e j T |
2 T |
|
||||||||||||
Расчет выполним для нормированной частоты |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
f |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
fд |
|
|
|
|
|
|||||
тогда |
|
|
|
|
|
д |
|
|
|
|
|
|
|
|||
|
|
|
1 e j4 |
|
|
|
|
1 e j4 |
||||||||
H j 2,29 10 3 |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
. |
||||||
|
j2 |
0,9347e |
j4 |
1 1,3152e |
j2 |
j4 |
||||||||||
|
1 1,4143e |
|
|
|
|
|
|
0,9301e |
Результатырасчетачастотнойхарактеристикифильтраприведенывтаблице2.11.
U1 z |
|
+ |
+ |
+ |
U2 z |
||||
+ |
|
|
|||||||
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
z |
1 |
|
|
|
z 1 |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
b11ц |
|
|
|
|
b12ц |
|
|
|
|
|
|
|
|
|
|
|
|
||
z 1 |
|
|
z 1 |
||||||
|
|
|
|
|
b21ц |
1 |
b22ц |
1 |
Рис. 2.14. Схема БИХ-фильтра
Табл. 2.11. Результаты расчета частотной характеристики фильтра
f |
f ,кГц |
|
Н1 |
1 |
Н2 |
2 |
Н |
|
0 |
0 |
0 |
0,00 |
0,00 |
0,00 |
0,00 |
0,000 |
0,00 |
fг |
10 |
0,0625 |
2,05 |
86,17 |
1,63 |
86,73 |
0,008 |
172,89 |
2fг |
20 |
0,125 |
21,61 |
-45,07 |
20,20 |
45,07 |
1,000 |
0,00 |
3fг |
30 |
0,1875 |
2,73 |
-84,88 |
3,18 |
-83,61 |
0,020 |
-168,49 |
28
Таб. 2.11 (Продолжение)
f |
f ,кГц |
|
Н1 |
1 |
Н2 |
2 |
Н |
|
4fг |
40 |
0,25 |
1,41 |
-87,35 |
1,52 |
-86,96 |
0,005 |
-174,31 |
5fг |
50 |
0,3125 |
0,86 |
-88,40 |
0,90 |
-88,20 |
0,002 |
-176,59 |
6fг |
60 |
0,375 |
0,51 |
-89,05 |
0,53 |
-88,94 |
0,001 |
-177,99 |
7 fг |
70 |
0,4375 |
0,24 |
-89,55 |
0,25 |
-89,51 |
0,000 |
-179,06 |
8fг |
80 |
0,5 |
0,00 |
0,00 |
0,00 |
0,00 |
0,000 |
0,00 |
9fг |
90 |
0,5625 |
0,24 |
89,55 |
0,25 |
89,51 |
0,000 |
179,06 |
10fг |
100 |
0,625 |
0,51 |
89,05 |
0,53 |
88,94 |
0,001 |
177,99 |
11fг |
110 |
0,6875 |
0,86 |
88,40 |
0,90 |
88,20 |
0,002 |
176,59 |
12fг |
120 |
0,75 |
1,41 |
87,35 |
1,52 |
86,96 |
0,005 |
174,31 |
13fг |
130 |
0,8125 |
2,73 |
84,88 |
3,18 |
83,61 |
0,020 |
168,49 |
14fг |
140 |
0,875 |
21,61 |
45,07 |
20,20 |
-45,07 |
1,000 |
0,00 |
15fг |
150 |
0,9375 |
2,05 |
-86,17 |
1,63 |
-86,73 |
0,008 |
-172,89 |
По результатам расчета построим график АЧХ цепи.
H f
1
-40 |
0 |
40 |
80 |
120 |
160 |
200 f ,кГц |
|
|
|
|
|
|
|
2 |
t, рад |
|
|
|
|
|
Рис. 2.15 АЧХ полосового БИХ-фильтра
В результате получена характеристика полосового фильтра с периодом повторения, равным частоте дискретизации.
Ослабление фильтра связано с частотной характеристикой выражением:
1
A 20lg H j .
В силу нелинейности преобразования 2.40 границы полосы пропускания и полосы непропускания дискретного БИХ-фильтра не будут совпадать с соответствующими значениями аналогового фильтра. Найдем их из соотношения
|
а |
tg |
fц |
, |
(2.49) |
|
fд |
||||||
|
|
|
|
|||
|
|
|
|
|
||
|
|
29 |
|
|
|
где а – частота аналогового фильтра, |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
fц – частота цифрового фильтра. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
Тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
fд |
|
2 fа |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
fц |
|
|
arctg |
|
. |
|
(2.50) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
В результате расчета получили: |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
f3ц |
|
f |
д |
|
|
|
|
|
2 f |
|
|
|
|
|
160 103 |
|
|
|
2 16 103 |
|
|
|
|||||||||||||||
|
|
|
arctg |
|
3 |
|
|
|
|
|
|
|
arctg |
|
|
|
|
|
16,30 кГц |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
303379 |
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
f2ц |
|
f |
д |
|
|
|
|
2 f |
|
|
|
|
|
|
160 103 |
|
|
|
|
2 19,025 103 |
|
||||||||||||||||
|
|
arctg |
|
2 |
|
|
|
|
|
|
|
arctg |
|
|
|
|
|
|
|
|
19,12 кГц |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
303379 |
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
f2ц |
|
fд |
arctg |
2 f2 |
|
|
160 103 |
arctg |
|
2 21,025 103 |
20,91кГц |
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
303379 |
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
f3ц |
|
|
fд |
arctg |
2 f3 |
|
|
160 103 |
arctg |
2 25 103 |
24,33кГц |
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
303379 |
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Определим ослабление в диапазоне частот f3ц f3ц , результаты запишем в таблицу 2.12.
Табл. 2.12. Результаты расчета ослабления фильтра
|
f ,кГц |
|
Н1 |
Н2 |
Н |
А |
|
16 |
0,1000 |
7,55 |
4,71 |
0,081 |
21,78 |
f3ц |
16,3 |
0,1019 |
8,37 |
5,05 |
0,097 |
20,28 |
|
16,5 |
0,1031 |
9,00 |
5,30 |
0,109 |
19,23 |
|
17 |
0,1063 |
11,00 |
6,04 |
0,152 |
16,36 |
|
17,5 |
0,1094 |
13,92 |
6,96 |
0,222 |
13,07 |
|
18 |
0,1125 |
18,32 |
8,17 |
0,343 |
9,30 |
|
18,5 |
0,1156 |
24,68 |
9,79 |
0,553 |
5,14 |
|
19 |
0,1188 |
30,30 |
12,06 |
0,837 |
1,55 |
f2ц |
19,12 |
0,1195 |
30,61 |
12,74 |
0,893 |
0,98 |
|
19,5 |
0,1219 |
28,06 |
15,36 |
0,987 |
0,11 |
f0 |
20 |
0,1250 |
21,61 |
20,20 |
1,000 |
0,00 |
|
20,5 |
0,1281 |
16,52 |
26,12 |
0,988 |
0,10 |
30