2799-up_ch2
.pdf141
7.11 Синтез цифрового ПФ Баттерворта методом инвариантной импульсной характеристики
Разработать цифровой ПФ с параметрами:
частота дискретизации – 100 кГц;
граничные частоты полосы пропускания – 2,5 кГц и 10 кГц;
граничные частоты полосы заграждения – 1,2 кГц и 20 кГц;
затухание АЧХ в пределах полосы пропускания – 3 дБ;
затухание АЧХ в пределах полосы заграждения – 15 дБ.
1.Расчет нормированных цифровых частот ПФ
ˆï 1 ï 1 ä fï 1 fä 2,5100 0,025;ˆï 2 ï 2 ä fï 2 fä 10100 0,1;ˆç1 ç1 ä fç1 fä 1,2100 0,012;ˆç2 ç2 ä fç2 fä 20 /100 0,2.
2. Расчет деформированных частот цифрового ФНЧ-прототипа
п частота цифрового ФНЧ-прототипа (обратите внимание на верхний индекс, обозначающий прототип)
Границу полосы пропускания цифрового НЧ-прототипа определим как половину полосы пропускания проектируемого ПФ, для того чтобы минимизировать деформирование частот при применении преобразования Константинидиса.
|
|
ˆ |
ï |
|
ˆ |
ï |
2 |
ˆï |
1 |
|
0,1 0,025 |
0,0375. |
|
|
|
|
|
|
|||
|
|
ï |
|
|
|
2 |
|
|
2 |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Нормированные |
|
|
граничные |
частоты |
полосы |
|
пропускания |
||||||||||||
проектируемого |
|
фильтра равны соответственно |
равны |
ˆï |
1 |
0,025 |
|
и |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ˆп |
2 |
0,1. Пересчет граничных частот полосы заграждения ПФ ˆç |
и ˆç |
2 |
в |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
частоты цифрового НЧ-прототипа ˆçï1 и ˆçï 2 проведем в соответствии с формулами таблицы 6.8.
|
|
|
|
1 |
|
|
1 |
arctg |
|
sin2 ˆ |
, ïðè |
ños2 ˆ ; |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
ños2 ˆ |
|
|
||||||||||
ˆ ï |
|
1 |
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
sin2 ˆ |
|
ños2 ˆ . |
||||||||||
|
|
|
|
|
|
|
arctg |
|
|
, |
ïðè |
|||||||||||||||
|
|
|
2 |
|
ños2 ˆ |
|||||||||||||||||||||
tg( ˆ |
ï |
|
|
) tg[ ( ˆ |
ï |
|
ˆ |
ï |
|
)] |
|
|
|
|
||||||||||||
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
cos ˆ |
ï |
2 |
|
ˆ |
ï |
1 |
|
2 |
|
|
1 |
|
|
|
|
|
|||||||||
cos ˆ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
ï |
|
ˆ |
ï |
|
|
|
|
|
|
|
|
|
||||||||||||||
2 |
|
1 |
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
142
tg( 0,0375) |
tg[ (0,1 0,025)] 0,1184 0,2401 0,4932 |
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||
|
cos 0,1 0,025 |
|
0,9239 |
0,9501 |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
cos 0,1 0,025 |
0,9724 |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
Проверим выполнение условия ños2 ˆ ( ) для ˆç |
и ˆç |
2 |
. |
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|||
ños2 ˆç |
|
|
|
|
|
ños(2 0,012) 0,9972 0,9501 |
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
ños2 ˆç |
2 |
|
|
ños(2 0,2) 0,309 0,9501 |
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
ˆçï1 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
sin2 ˆç |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
arctg |
ños2 ˆç1 |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
1 |
arctg |
|
|
|
|
|
|
|
|
sin(2 0,012) |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
0,4932 [ños(2 0,012) |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
0,9501] |
|
|
|
|
|
|
|
||||||||||||||||||||||||
0,5 |
|
|
1 |
|
arctg3,2386 0,5 |
|
1,2713 |
0,0953 |
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
ˆ |
ï |
|
|
|
|
1 |
|
|
|
|
|
1 |
arctg |
|
|
|
|
|
sin2 ˆç |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
ç2 |
|
|
|
|
|
|
|
|
|
|
|
|
ños2 ˆç2 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
1 |
|
|
1 |
|
arctg |
|
|
|
|
|
|
|
sin(2 0,2) |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,4932 [ños(2 0,2) 0,9501] |
|
|
|
|
|
|
|
||||||||||||||||||||
0,5 |
|
1 |
arctg( 3,0079) 0,5 |
|
1,2498 |
0,1022 |
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
В условиях поставленной задачи пересчет граничных частот полосы |
|||||||||||||||||||||||||||||||||||||||||||||||||
заграждения ПФ ˆç |
и ˆ |
ç |
2 |
в частоты цифрового НЧ-прототипа ˆ |
ï |
и ˆ |
ï |
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ç |
|
ç |
2 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
можно провести в соответствии с таблицей 9.2.1.
ˆç1 0,012 3-я строка, 4-й столбец таблицы 9.2.1 ˆçï1 0,095.
ˆç2 0,2 22-я строка, 2-й столбец таблицы 9.2.1 ˆçï 2 0,1021.
Выберем из двух возможных значений ˆçï1 и ˆçï 2 то, которое
соответствует более жесткому требованию по порядку, то есть меньшее по модулю ˆçï 0,095.
Отношение граничных частот цифрового ФНЧ-прототипа и отношение граничных частот аналогового НЧ-прототипа совпадают, то есть
ˆ ç ˆ ï ˆçï ˆïï .
3. Расчет порядка аналогового ФНЧ-прототипа Баттерворта
|
|
|
lg |
|
100,1Aç |
1 |
lg |
100,1 15 1 |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
N |
Á |
|
|
100,1Aï |
1 |
|
|
100,1 3 1 |
|
|
0,7441 |
1,8432 |
||||||||||
|
|
|
|
|
|
|||||||||||||||||
|
|
|
ˆ |
ˆ |
|
|
0,095/ 0,0375 |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
lg |
0,4037 |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
lg |
ç |
ï |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
143
Примем порядок NБ =2.
4. Определение импульсной характеристики аналогового ФНЧ-
прототипа |
Баттерворта, нормированной относительно граничной частоты |
|||||||
|
|
|
ˆ |
|
|
|
|
|
полосы пропускания ï согласно таблице 9.3.9 |
|
|||||||
g t |
|
п |
Ae a t sin(w t) , где |
t |
п |
t , где |
п |
п . |
|
|
a |
|
|
п |
|||
A 1,41589356, a 0,70794678, |
wa 0,70794678. |
|
5. Денормирование импульсной характеристики Учитывая, что импульсная характеристика аналогового НЧ-прототипа
нормирована относительно граничной частоты полосы пропускания п , ее
нужно преобразовать следующим образом: g(t) п Ae a t п sin(wa t п )
6. Дискретизация импульсной характеристики g(n) Tд g t n Tд
|
ï |
T |
2 ˆ ï |
, поскольку |
п |
|
п . |
|
|
|||||||||
|
ä |
|
|
|
ï |
|
|
|
|
|
|
|
п |
|
|
|||
|
|
|
|
|
|
ˆ ï |
|
2 a ˆïï |
n |
|
ˆ ï |
|
, |
|||||
g(n) A 2 ï |
e |
|
|
|
sin 2 wa ï |
n |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где ˆ ï |
|
ˆï |
2 |
ˆï |
1 |
|
|
0,1 0,025 |
0,0375. |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
ï |
|
|
|
2 |
|
|
|
|
2 |
|
|
|
|
|
|
|
g(n) 0,3336 e0,1668 n sin 0,1668 n |
|
|
g(n) 0,3336 0,8464 n sin 0,1668 n
7. Определение системной функции цифрового ФНЧ-прототипа в
соответствии с таблицей 2.2 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
Ksin (z) |
|
|
|
|
z 1 [0,3336 e 0,1668 sin(0,1668) |
, |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
1 2z 1 e 0,1668 cos(0,1668) z 2 (e 0,1668)2 |
|
||||||||||||||||||||||
Kфнч (z) |
|
|
|
|
|
0,0469 z 1 |
|
|
|
|
|
|
|
||||||||||||
|
1,6693 z 1 0,7164 z 2 |
|
|
|
|
|
|
||||||||||||||||||
1 |
|
|
|
|
|
|
|||||||||||||||||||
8. Расчет параметров преобразования Константинидиса |
|
||||||||||||||||||||||||
|
cos[ ( ˆï |
2 |
ˆ |
ï |
1 |
)] |
|
|
cos[ (0,1 0,025)] |
0,9239 |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,9501 |
||||||
cos[ ( ˆï |
2 |
ˆ |
ï |
1 |
)] |
cos[ (0,1 0,025)] |
|
0,9274 |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
tg( ˆ ï ) |
|
|
|
|
|
|
|
|
tg(0,0375 ) |
0,1184 |
|
|
|
||||||||||
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
0,493 |
||||||
tg[ ( ˆï |
2 |
ˆï |
1 |
)] |
tg[ (0,1- 0,025)] |
|
0,2401 |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
144 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
1 |
|
|
2 |
|
z 1 |
z 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
0,3396 0,6275 z 1 z 2 |
|
|
|
|
|
|
||||||||||||||||||||||||||
|
z 1 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
1 |
|
2 |
z 1 |
|
|
1 |
z 2 |
1 |
0,6275 z 1 |
0,3396 z 2 |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
9. Определение системной функции цифрового ПФ Баттерворта |
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,0469 |
|
|
0,3396 0,6275 z-1 z 2 |
|
|
|
|
|
|
|
|
|||||||||||||||||||
Kпф (z) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
0,6275 z 1 0,3396 z 2 |
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
2 2 |
|||||||||
1 |
|
|
|
|
|
|
|
|
|
|
0,3396 0,6275z |
z |
|
|
|
0,3396 0,6275z |
z |
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
1 1,6693 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,7164 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
1 |
|
|
|
|
2 |
||||||||||||||
|
|
|
|
|
|
|
|
1 0,6275z |
0,3396z |
|
|
|
|
|
0,6275z |
0,3396z |
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
||||||||||||||||||||||||||||||
|
Kпф |
|
|
(z) |
0,0309 0,0377 z 1 0,1372 z 2 |
0,0377 z 3 0,0309 z 4 |
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
1 |
|
1 3,1826 z 1 3,9342 z 2 |
2,2578 z 3 |
0,5134 z 4 |
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Амплитудно-частотная |
|
характеристика |
|
Kïô |
1 |
( ˆ ) |
|
получена |
из |
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||
системной |
|
функции |
K |
пф |
(z) |
|
при |
|
z e j Tä |
e j2ˆ |
и |
представлена |
на |
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
рисунке 7.11.
Рисунок 7.11 – Амплитудно-частотные характеристики двух ПФ, рассчитанных в пунктах 7.11 и 7.12
7.12 Синтез цифрового ПФ Баттерворта методом инвариантной импульсной характеристики
Разработать цифровой ПФ с параметрами:
частота дискретизации – 100 кГц;
граничные частоты полосы пропускания – 21,25 кГц и 28,75 кГц;
145
граничные частоты полосы заграждения – 15 кГц и 40 кГц;
затухание АЧХ в пределах полосы пропускания – 3 дБ;
затухание АЧХ в пределах полосы заграждения – 17 дБ.
1.Расчет нормированных цифровых частот ПФ
ˆï 1 ï 1 ä fï 1 fä 21,25100 0,2125;ˆï 2 ï 2 ä fï 2 fä 28,75100 0,2875;ˆç1 ç1 ä fç1 fä 15100 0,15;
ˆç2 ç2 ä fç2 fä 40 /100 0,4.
2.Расчет деформированных частот цифрового ФНЧ-прототипа
Граница полосы пропускания цифрового НЧ-прототипа определим как половину полосы пропускания проектируемого ПФ, для того чтобы минимизировать деформирование частот при применении преобразования Константинидиса.
ˆ |
ï |
|
ˆ |
ï |
2 |
ˆï |
1 |
|
0,1 0,025 |
0,0375. |
|
ï |
|
|
|
2 |
|
2 |
|
||||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
||
За |
неимением |
пересчетных таблиц для случая проектирования |
цифрового ПФ с граничными частотами полосы пропускания, соответственно
равными ˆï |
1 |
|
|
0,2125 |
|
и |
ˆï |
2 |
|
0,2875 |
, |
|
расчет деформированных |
частот |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
цифрового ФНЧ-прототипа ведется по аналогии с примером 7.11. |
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
Пересчет граничных частот полосы заграждения ПФ |
ˆç |
|
и |
ˆç |
2 |
в |
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
частоты цифрового НЧ-прототипа ˆ ï |
и |
ˆ ï |
|
проведем в соответствии с |
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ç |
|
|
|
ç |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
формулами таблицы 6.8. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
1 |
|
|
1 |
arctg |
|
|
|
sin2 ˆ |
|
|
, åñëè |
ños2 ˆ ; |
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
ños2 ˆ |
|
|
|
|
|
|
|
|
|
||||||||||||||||
ˆ ï |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
sin2 ˆ |
|
|
|
|
ños2 ˆ . |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
arctg |
|
|
|
|
|
, |
åñëè |
|
|
|
|
|
||||||||||||||||||||||
|
|
2 |
|
ños2 ˆ |
|
|
|
|
|
|||||||||||||||||||||||||||||||
tg( ˆ |
ï |
|
|
) tg[ ( ˆ |
ï |
|
ˆ |
ï |
|
|
)] |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
cos ˆ |
ï |
2 |
|
ˆ |
ï |
1 |
|
2 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
ˆ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
cos |
ï |
|
ˆ |
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
2 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
tg( 0,0375) |
|
tg[ (0,2875 0,2125)] 0,1184 0,2401 0,4931 |
|
|
|
|||||||||||||||||||||||||||||||||||
|
cos 0,2875 0,2125 |
|
|
|
|
0 |
|
0 |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
cos 0,2875 0,2125 |
|
0,9724 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
Проверим выполнение условия ños2 ˆ ( ) для ˆç |
и ˆç |
2 |
. |
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
ños2 ˆç |
|
|
|
ños(2 0,15) 0,5878 0 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
ños2 ˆç |
2 |
|
ños(2 0,4) 0,809 0 |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
146
ˆçï1 |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|
sin2 ˆç |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
arctg |
ños2 ˆç1 |
|
|
|
|
|
||||||||||||||||||||
|
|
2 |
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
1 |
arctg |
|
|
|
sin(2 0,15) |
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
0,4931 [ños(2 0,15) 0] |
|
|
|
|
||||||||||||||
0,5 |
|
1 |
|
arctg2,7913 0,5 0,3905 0,1095 |
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
ˆ |
ï |
|
|
|
1 |
|
|
|
|
1 |
arctg |
|
|
sin2 ˆç |
2 |
|
|
|
|
|
|
|
||||||||||||||
ç2 |
|
|
|
|
|
|
|
|
ños2 ˆç2 |
|
|
|
|
|||||||||||||||||||||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
1 |
|
|
1 |
|
arctg |
|
|
|
|
|
sin(2 0,4) |
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
0,4931 [ños(2 0,4) |
|
|
|
|
|
|||||||||||||||||||||||||
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0] |
|
|
|
|
|
|||||||||||
0,5 |
|
1 |
arctg( 1,4734) 0,5 |
0,9745 0,1898 |
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Выберем |
|
|
из двух |
возможных |
значений ˆ |
ï |
и ˆ |
ï |
то, которое |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ç |
|
ç |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
соответствует более жесткому требованию по порядку, то есть меньшее по модулю зп 0,1095.
Отношение граничных частот цифрового ФНЧ-прототипа и отношение граничных частот аналогового НЧ-прототипа совпадают, то есть
ˆ ç ˆ ï ˆçï ˆïï .
3. Расчет порядка аналогового ФНЧ-прототипа Баттерворта
|
|
|
lg |
|
100,1Aç |
1 |
|
|
lg |
100,1 17 1 |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
N |
Á |
|
|
100,1Aï |
1 |
|
|
|
100,1 3 1 |
|
|
0,8467 |
1,8196 |
|||||||||
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
ˆ |
ˆ |
|
|
|
|
|
0,095/ 0,0375 |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
lg |
0,4653 |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
lg |
ç |
ï |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Примем порядок NБ 2 .
4. Определение импульсной характеристики аналогового ФНЧ-
прототипа Баттерворта, нормированной относительно граничной частоты |
|||||||
полосы пропускания ˆ ï согласно таблице 9.3.9 |
|
|
|||||
|
|
ï |
|
|
|
|
|
g t |
п |
Ae a t sin(w t) , где |
t |
п |
t , где |
п |
п . |
|
a |
|
|
п |
|||
A 1,41589356, a 0,70794678, |
wa 0,70794678. |
|
5. Денормирование импульсной характеристики Учитывая, что импульсная характеристика фильтра-прототипа
нормирована относительно граничной частоты полосы пропускания ˆïï , ее
следует преобразовать следующим образом g(t) п Ae a t п sin(wa t п )
6. Дискретизация импульсной характеристики g(n) Tд g t n Tд
147
|
ï |
T |
2 ˆ ï |
, поскольку |
|
п |
|
п . |
|
|
||||||||||
|
ä |
|
|
|
ï |
|
|
|
|
|
|
|
п |
|
|
|
||||
|
|
|
|
|
|
ˆ ï |
|
2 a ˆïï |
n |
|
|
|
|
ˆ ï |
|
, |
||||
g(n) A 2 ï |
e |
|
|
|
|
sin 2 wa ï |
n |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где ˆ ï |
|
ˆï |
2 |
ˆï |
1 |
|
|
0,2875 0,2125 |
0,0375. |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
ï |
|
|
|
2 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
g(n) 0,3336 e0,1668 n sin 0,1668 n |
|
|
g(n) 0,3336 0,8464 n sin 0,1668 n
7. Определение системной функции цифрового ФНЧ-прототипа в
соответствии с таблицей 2.2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
Ksin (z) |
|
|
|
|
|
|
|
z 1 [0,3336 e 0,1668 sin(0,1668) |
|
, |
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
1 2z 1 e 0,1668 cos(0,1668) z 2 (e 0,1668)2 |
|
|
||||||||||||||||||||||||||||||||||||||||||
Kфнч |
(z) |
|
|
|
|
|
|
|
0,0469 z 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
1 1,6693 z 1 0,7164 z 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
8. Расчет параметров преобразования Константинидиса |
|
|
|||||||||||||||||||||||||||||||||||||||||||||||
|
|
cos[ ( ˆï |
2 |
ˆï |
1 |
)] |
|
|
cos[ (0,2875 0,2125)] |
|
|
|
|
|
|
0 |
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
||||||||
cos[ ( ˆ |
|
|
|
ˆ |
1 |
|
)] |
|
cos[ (0,2875 0,2125)] |
0,9274 |
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
ï 2 |
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
tg( ˆ ï |
) |
|
|
|
|
|
|
|
|
|
tg(0,0375 ) |
0,1184 |
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,493 |
|
|
||||
tg[ ( ˆï |
2 |
|
ˆï |
1 |
)] |
|
tg[ (0,1- 0,025)] |
0,2401 |
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
1 |
|
|
2 |
|
z 1 z 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
0,3396 z 2 |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
z 1 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
1 |
|
|
1 |
0,3396 z 2 |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
1 |
z 1 |
|
z 2 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
9. Определение системной функции цифрового ПФ Баттерворта |
|||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,3396 z |
2 |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,0469 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 0,3396 z |
2 |
|
|
|
|
|
||||||||||||||
Kпф |
|
(z) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
0,3396 z |
2 |
|
|
|
|
|
|
|
|
0,3396 z 2 |
2 |
||||||||||||||||||||||||
|
|
|
|
|
|
|
1 1,6693 |
|
|
|
|
|
|
|
|
|
|
|
|
|
0,7164 |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
2 |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 0,3396 z |
|
|
|
|
|
|
|
1 |
0,3396 z |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
Kпф 2 |
(z) |
0,0309 0,1014 z 2 0,0309 z 4 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
1 3,9495 z |
2 |
0,5134 z |
4 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
148
Амплитудно-частотная |
характеристика |
|
Kïô |
2 |
( ˆ ) |
получена |
из |
||
системной функции K |
пф |
(z) |
при z e j Tä |
e j2ˆ |
и |
представлена |
на |
||
|
2 |
|
|
|
|
|
|
|
|
рисунке 7.11.
149
8СИНТЕЗ ТРАНСВЕРСАЛЬНЫХ ЦИФРОВЫХ ФИЛЬТРОВ
8.1Метод частотной выборки
Трансверсальный цифровой фильтр имеет импульсную характеристику конечной длины N, и его называют КИХ-фильтром.
Стандартный КИХ-фильтр характеризуется следующими уравнениями:
|
N 1 |
|
|
|
|
|
yn |
gk xn k , |
|
|
|
(8.1) |
|
|
k 0 |
|
|
|
|
|
K z |
N 1 |
gn z n , |
|
|
|
|
|
|
|
|
(8.2) |
||
|
n 0 |
|
|
|
|
|
|
|
N 1 |
|
|
|
|
|
( ˆ ) |
gn (e |
j2 ˆ |
) |
n |
(8.3) |
Kêèõ |
|
|
||||
|
|
n 0 |
|
|
|
|
где gn - коэффициенты импульсной характеристики ЦФ; |
|
|||||
K z - системная функция ЦФ, |
z e j Tä e j2 ˆ ; |
|
Kêèõ ( ˆ ) - комплексная передаточная функция КИХ-фильтра.
Отсчеты импульсной характеристики КИХ-фильтров одновременно являются и коэффициентами системной функции ЦФ, поэтому задача синтеза
сводится к получению импульсной характеристики ЦФ. |
|
|
|
g t и |
|||||
В аналоговых линейных системах импульсная характеристика |
|||||||||
комплексная передаточная |
функция |
K линейной |
электрической |
цепи |
|||||
|
|
|
|
|
|
|
|
|
|
(ЛЭЦ) связаны между собой двойным интегралом Фурье: |
|
|
|
||||||
|
|
1 |
|
j t |
|
|
j t |
|
|
K |
|
K e |
|
d |
e |
|
dt |
(8.4) |
|
|
0 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
g t |
|
|
|
|
|
|
Графическая иллюстрация импульсных и частотных характеристик идеального аналогового ФНЧ представлена на рисунке 8.1.
|
|
150 |
|
|
gи(t) |
|
|
|K(ω)| |
|
|
|
|
|
|
0 |
|
t |
0 |
ω |
а) |
|
|
б) |
|
g(t) |
|
|
|
|
|
|
|
argK(ω) |
|
0 |
tз |
t |
0 |
ω |
в) |
|
|
г) |
|
Рисунок 8.1 – б) АЧХ г) ФЧХ идеального ФНЧ и импульсные характеристики а) без учета ФЧХ и в) с учетом линейной ФЧХ
Переход от аналоговых линейных систем к дискретным сопровождается заменой аналоговых переменных t и двойного интеграла Фурье (8.4) на дискретные переменные n и m дискретных преобразований Фурье (ДПФ)
Для этого периодическую частотную характеристику идеального цифрового фильтра дискретизируют, разбивая интервал (0; д ) на N равных частей ( д N 1). Интервал взятия выборки на частотной оси 1 должен быть таким, чтобы передать все особенности АЧХ. При этом импульсная
характеристика НЧ-прототипа становится |
периодической |
(с периодом |
Tоg 2 1 ), дискретизированной с шагом |
дискретизации |
Tд 2 д и |
периодом Tоg , разбитым на N равных частей Tоg N Tд .
Графическая иллюстрация дискретизированных импульсных и частотных характеристик идеального цифрового ФНЧ представлена на рисунке 8.2.