Введение в теорию сигналов. Цифровая обработка сигналов. Учебное пособие
.pdf
|
1, |
|
|
|
|
при n0 |
2l0 ; |
|
V00 |
|
|
00 |
|
|
|
|
|
(s0 ) |
|
|
|
, 00 |
00 , при n0 |
2l0 1. |
||
|
|
|
|
|
|
|||
|
(s |
|
00 |
) |
||||
|
|
0 |
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
0i |
, 2 |
2 |
2 |
|
|
|
|||||
|
|
|
|
|
|
|||||||||
|
|
(s |
0i |
)2 2 |
0i |
|
|
0i |
0i |
|
|
|
||
V0i (s0 ) |
|
0 |
|
0i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
для фильтровБаттервортаиЧебышёва; |
||||||||||
|
|
2 |
s2 |
2 |
|
|
|
2 |
2 |
|
|
|
||
|
|
2 |
|
|
|
|
||||||||
|
|
0i |
|
0 |
0i |
|
|
0i |
0i |
|
|
|
||
|
|
|
|
|
, 0i |
|
|
2 |
для фильтровКауэра. |
|||||
|
(s |
|
)2 2 |
|
||||||||||
|
|
0 |
|
0i |
0i |
|
|
|
|
|
0i |
|
|
|
Для сравнения в табл. 2.1 приведены порядки фильтров полиномов |
||||||||||||||
Баттерворта, |
Кауэра и Чебышёва при 0,02; 0,01 в зависимости |
|||||||||||||
от параметра . |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 2.1 |
|
|
|
|
Зависимость порядка n0 от |
|
|||||||||
|
|
|
|
|
Порядок n0 |
полинома фильтра |
|
|||||||
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
||||||
v |
|
Баттерворта |
|
|
|
Чебышёва |
|
Кауэра |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1,05 |
|
|
|
127 |
|
|
|
|
|
22 |
|
|
8 |
|
1,10 |
|
|
|
65 |
|
|
|
|
|
16 |
|
|
7 |
|
1,20 |
|
|
|
34 |
|
|
|
|
|
11 |
|
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1,40 |
|
|
|
19 |
|
|
|
|
|
8 |
|
|
5 |
|
1,60 |
|
|
|
14 |
|
|
|
|
|
7 |
|
|
5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1,80 |
|
|
|
11 |
|
|
|
|
|
6 |
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2,00 |
|
|
|
|
9 |
|
|
|
|
|
6 |
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Следует отметить, что параметры передаточной функции фильтра Кауэра очень чувствительны к вычислительным погрешностям. В результате нули и полюса фильтра Кауэра необходимо вычислять с удвоенной точностью. В противном случае задача проектирования фильтра
81
становится некорректно поставленной, поскольку небольшие погрешности при вычислении полюсов и нулей приводят к существенным отклонениям в АЧХ. Более того, возникает проблема устойчивости рекурсивного фильтра Кауэра.
2.4.5. ПЕРЕХОД К ПЕРЕДАТОЧНОЙ ФУНКЦИИ НЕПРЕРЫВНОГО ФИЛЬТРА
Переход от передаточной функции W0 (s0 ) непрерывного НЧФ-П к передаточной функции W (s) требуемого непрерывного фильтра осуществляется на основе представления W0 (s0 ) в виде произведения простых сомножителей (2.98). В результате будем иметь
|
|
|
|
|
|
|
|
l |
|
|
|
|
|
|
|
|
|
W (s) V0 (s) Vi (s) , |
|
|
|
(2.99) |
|||||||
|
|
|
|
|
|
|
|
i 1 |
|
|
|
|
|
|
где 0kф , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1, |
|
|
|
|
|
|
n0 2l0 , |
|
|
|
|
||
|
b |
s b |
|
|
|
|
|
|
|
|
|
|
|
|
V (s) |
02 |
01 |
, |
|
|
n |
2l 1 для НЧФ, ВЧФ, |
|||||||
|
|
|
|
|||||||||||
0 |
|
s |
0 |
|
|
|
|
|
0 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
b |
s2 b |
|
s b |
|
|
|
|
|
|
|
|
||
|
|
03 |
02 |
01 |
, |
|
n0 2l0 1 для ПФ, РФ, |
|||||||
|
|
(s 0 ) |
2 |
2 |
|
|||||||||
|
|
|
|
0 |
|
|
|
|
|
|
|
|
||
l |
l |
для НЧФ, ВЧФ, |
V |
(s) bi3s |
2 |
bi2s bi1 . |
||||||||
|
0 |
|
|
|
|
|
|
|||||||
|
|
|
для ПФ, РФ, |
|
|
i |
(s |
)2 |
2 |
|||||
|
2l |
|
|
|
||||||||||
|
|
0 |
|
|
|
|
|
|
|
|
|
i |
|
i |
Параметры передаточной функции |
|
|
|
|
|
|
||||||||
|
|
W (s) i , i , bi3 , bi2 , bi1, i [0, l] |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
, |
ф при замене s0 следующим |
||||
вычисляются по исходным данным ф |
||||||||||||||
образом:
82
|
s |
|
|
|
|
|
|
дляНЧФ, |
|
|
|
|
|
|
|
|
|
ф |
|
|
|
|||||
|
|
|
|
|||||
|
|
|
|
|
||||
s2 |
|
|
|
|||||
|
|
|
|
|
|
ф |
ф |
дляПФ, |
|
( |
|
|
|
)s |
|||
|
|
ф |
|
|
ф |
(2.100) |
||
s0 |
|
|
|
|
|
|
|
|
( ф ф)s |
дляРФ, |
|||||||
|
|
|
|
|
||||
s2 |
|
|
||||||
|
|
|
|
|
|
ф |
ф |
|
|
ф |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
дляВЧФ. |
|
|
s |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Вкачествепримерарассмотримопределениепараметровпередаточной функции для НЧФ.
Осуществляется переход от к в виде ф, при этом:
* ф,п ( * 1) ф,
i i0 ф, vi vi0 ф, 0 00 ф,
b01 0 , b02 0.
Дополнительно для фильтров Кауэра вычисляются:
bi1 02i ф2 ф2 , bi2 0, bi3 02i ,
а для всех остальных фильтров
bi1 02i ф2 , b01 b02 0.
Аналогично определяются параметры W (s) и у ВЧФ. Для ПФ и РФ вычисления сложнее приблизительно в два раза.
83
2.4.6. ПЕРЕХОД К ПЕРЕДАТОЧНОЙ ФУНКЦИИ РЕКУРСИВНОГО ФИЛЬТРА
Переход к передаточной функции рекурсивного фильтра W (z) осуществляется по билинейному преобразованию, когда W (s) представ-
лена в виде произведения простых сомножителей. В соответствии с (2.99) имеем
l |
|
|
W (s) V0 (s) Vi (s) |
(2.101) |
|
i 1 |
|
|
и в результате преобразования получим |
|
|
l |
(u) , u z 1 , |
|
W (u) V 0 (u) V i |
(2.102) |
|
i 1
где t , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 b |
u b |
|
u2 |
|
|
|
|
V |
(u) |
i1 |
i2 |
|
|||
|
|
|
1 a u a |
u2 |
||||||
|
|
|
i |
|
i |
|||||
|
|
|
|
|
|
i1 |
i |
2 |
|
|
|
1, |
|
|
|
|
|
n0 2l0 ; |
|||
V |
|
|
1 b |
u |
, |
|
n 2l 1 |
|||
(u) |
|
01 |
|
|||||||
0 |
|
0 |
1 a u |
|
|
0 |
0 |
|
|
|
|
|
|
01 |
|
|
|
|
|
|
|
|
|
|
1 b |
u b |
u2 |
|
|
|
|
|
|
|
|
01 |
02 |
u2 , n0 2l0 1 |
|||||
|
0 |
1 a u |
a |
|||||||
|
|
|
01 |
02 |
|
|
|
|
|
|
Определение параметров i , |
a i1, |
a i2 , |
b i2 , |
|||||||
ляется при замене s zz 11 .
,
дляНЧФ, ВЧФ,
дляПФ, РФ.
b i1, i [0, l], осуществ-
2.4.7. РЕАЛИЗАЦИЯ РЕКУРСИВНЫХ ФИЛЬТРОВ
Реализовать рекурсивный фильтр можно по трем схемам:
1)прямой;
2)каскадной;
3)параллельной.
84
В случае прямой схемы передаточная функция рекурсивного фильтра должна быть представлена в виде суммы простых дробей или произведения простых сомножителей. В любом случае она приводится к виду
|
|
|
|
|
|
n |
b ui |
|
|
n |
|
, |
(2.103) |
W (u) W (u) , |
W (u) |
|
1 |
a ui |
|
||||||||
|
|
|
|
|
|
i |
|
|
i |
|
|
||
|
|
|
|
|
|
i 0 |
|
|
|
i 1 |
|
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2l, |
n0 2l0 , |
|
|
|
|
|
|
|
|
||
|
n |
|
n0 2l0 1 |
дляНЧФ, ВЧФ, |
|
|
|
||||||
|
2l 1, |
|
|
|
|||||||||
|
|
2l 2, |
n |
2l |
2 |
дляПФ, РФ. |
|
|
|
||||
|
|
|
0 |
0 |
|
|
|
|
|
|
|
|
|
Структурная схема фильтра представлена на рис. 2.18.
Рис. 2.18. Прямая форма реализации рекурсивного фильтра
Отсчеты выходной переменной определяются при этом как
|
|
n |
n |
|
(2.104) |
yk yk , |
yk bi x(k i) ai y(k i) . |
||||
|
|
i 0 |
i 1 |
|
|
Недостатком прямой схемы реализации является высокая чувствительность полюсов и нулей передаточной функции фильтра к погреш-
ностямзаданияеекоэффициентов ai и bi . Чувствительностьтемвыше,
чем больше порядок фильтра n. В итоге даже вычислительные погрешности могут привести к неустойчивости фильтра, реализуемого по прямой схеме.
При каскадной схеме реализации передаточная функция представляется в виде произведения простых сомножителей:
l
W (u) V 0 (u) V i (u) ,
1
85
а соответствующая схема показана на рис. 2.19.
y(l – 1)k |
|
ylk |
|
|
|
|
κ |
||
|
|
|
|
|
|
Рис. 2.19. Каскадная схема реализации рекурсивного фильтра
Отсчеты выходного сигнала определяются согласно выражению
yk ylk , |
(2.105) |
где ylk находится по рекуррентному i [1, l] соотношению |
|
y |
y |
|
b y |
|
|
b y |
|
|
a y |
a y |
, |
||||
ik |
i |
(i 1)k |
|
i1 (i 1)(k 1) |
i2 (i 1)(k 2) |
i1 i(k 1) |
i2 i(k 2) |
|
|||||||
а y0k |
определяется в виде |
|
|
|
|
|
|
|
|
|
|||||
|
x |
|
, |
|
|
|
|
|
|
|
n |
2l |
, |
|
|
|
k |
x |
|
|
|
|
|
|
|
0 |
0 |
|
|
|
|
|
|
|
b |
x |
a |
y |
, |
n |
2l |
1 для НЧФ, ВЧФ, |
|||||
|
|
0 k |
|
01 (k 1) |
|
01 0(k 1) |
|
0 |
0 |
|
|
|
|||
y0k |
|
|
b01x(k 1) |
|
|
|
a01 y0(k 1) a02 y0(k 2) , |
|
|||||||
|
|
|
|
|
|
||||||||||
|
0 xk |
b02 x(k 2) |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
n |
2l |
1 для ПФ, РФ. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
0 |
0 |
|
|
|
Каскадная схема в сравнении с прямой более устойчива к вычислительным погрешностям.
Параллельная схема требует представления передаточной функцией фильтра W (u) в виде суммы простых дробей:
W (u) |
b |
W |
(u) |
l |
W |
(u) |
|
, |
|
|
|
|
0 |
|
|
i |
|
|
|
|
|
|
|
|
i 1 |
|
|
|
|
а дискретные значения выходного сигнала вычисляются по выражению
y |
|
b x |
|
l |
y |
|
, |
||
k |
|
|
|
k |
|
ik |
|
||
|
|
|
|
|
|
i 0 |
|
|
|
где yik , i [1, l], находятся рекуррентно, согласно
yik i[xk ci1x(k 1) ci2 x(k 2) ] ai1 yi(k 1) ai2 yi(k 2) ,
86
причем
0, при n0 2l0 ,
0[xk c01xi(k 1) ] a01 y0(k 1) ,
y0k при n0 2l0 1 для НЧФ и ВЧФ,
0[xk c01x(k 1) c02 x(k 2) ] a01 y0(k 1) a02 y0(k 2) ,
при n0 2l0 1 для ПФиРФ.
Параллельная схема реализации изображена на рис. 2.20.
W l (u) |
ylk |
κ |
|
ylk
W l (u) ylk
Рис. 2.20. Параллельная форма реализации рекурсивного фильтра
Параллельная схема также более устойчива к вычислительным погрешностям.
87
2.4.8.СРАВНИТЕЛЬНЫЕ ХАРАКТЕРИСТИКИ РЕКУРСИВНЫХ ФИЛЬТРОВ
Вкачестве критерия для сравнения будем использовать эффективность фильтра (количество вычислительных операций, необходимых для вычисления одного отсчета выходного сигнала фильтра) при заданных качествах его характеристик. В результате получается, что эффективность рекурсивных фильтров будет определяться:
1) схемной реализацией рекурсивного фильтра;
2) видом функции, аппроксимирующей фильтр-прототип;
3) способом проектирования рекурсивного фильтра (используется стандартное или билинейное Z-преобразование).
Если обозначить через l количество параметров передаточной функции рекурсивного фильтра, то каскадная схема реализации требует
5(l 1) 1 операций умножения и 4(l 1) 2 операций сложения, а параллельная: 5(l 1) 2 операций умножения и 5(l 1) операций сложе-
ния. Кроме того, в литературе указывается, что каскадная схема дает меньшую ошибку фильтрации.
Прямая схема реализации фильтров не распространена из-за проблемы ее устойчивости.
С точки зрения вида аппроксимирующей функции самыми эффективными являютсяфильтрыЧебышёваиблизкиек нимфильтрыКауэра (особенно Кауэра, см. табл. 2.1). Фильтры Баттерворта используются при строгих требованиях к неравномерности в полосе пропускания АЧХ. Однако фильтры Кауэра теряют свое преимущество при их практической реализации из-за сложностей с вычислением (они требуют двойной точности).
При проектировании рекурсивного фильтра, как уже упоминалось выше, удобнее использовать билинейное Z-преобразование.
Таким образом, предпочтение следует отдать рекурсивным фильтрам Кауэра (Чебышёва) с использованием билинейного Z-преобразования, реализуемымпо каскаднойсхеме.
88
2.5.СРАВНЕНИЕ РЕКУРСИВНЫХ
ИНЕРЕКУРСИВНЫХ ФИЛЬТРОВ
Критерии сравнения фильтров следующие:
1)качество АЧХ фильтра;
2)линейность ФЧХ;
3)устойчивость фильтра;
4)эффективность.
Результаты сравнения отражены в табл. 2.2.
|
|
Таблица 2.2 |
|
|
|
Параметр |
Нерекурсивные |
Рекурсивные |
сравнения |
фильтры |
фильтры |
Качество АЧХ фильтра |
+ |
+ |
Линейность ФЧХ |
+ |
?– |
|
|
|
Устойчивость фильтра |
+ |
?– |
|
|
|
Эффективность |
– |
+ |
Требуемое качество АЧХ обеспечивают оба вида фильтров. Проблемы устойчивости у нерекурсивных фильтров не существует. Рекурсивные фильтры могут быть и устойчивыми, и неустойчивыми, что является их недостатком, который следует учитывать при проектировании. Количество вычислений в случае нерекурсивного фильтра определяется длительностью его ИХ. Если обозначить число параметров передаточной функции нерекурсивного фильтра через L, а через l – рекурсивного, то при одном и том же качестве АЧХ будет иметь место соотношение L 5(l 1) .
Таким образом, можно сделать вывод, что при отсутствии требований к линейности ФЧХ отдается предпочтение рекурсивным фильтрам при условии их устойчивости.
Нерекурсивные фильтры используются при отсутствии переходных и фазовых искажений и требованиях к линейности ФЧХ.
89
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1.Френкс Л. Теория сигналов / Л. Френкс. – М.: Советское радио, 1974,
344 с.
2.ГоноровскийИ.С. Радиотехническиецепиисигналы: учебникдлявузов / И.С. Гоноровский. – М.: Радио и связь, 1986. – 512 с.
3.Баскаков С.И. Радиотехнические цепи и сигналы / С.И. Баскаков. – М.: Высшая школа, 2000. – 462 с.
4.Анисимов А.С. Методы цифровой фильтрации: учебное пособие / А.С. Анисимов. – Новосибирск: Изд-во НГТУ, 1992. – 82 с.
5.КононовВ.Т. Цифроваяфильтрация сигналов: учебное пособие/ В.Т. Кононов, Д.С. Худяков, Г.П. Чикильдин. – Новосибирск: Изд-во НГТУ, 2008. – 64 с.
6.Рабинер Л. Теория и применение цифровой обработки сигналов: пер.
сангл. / Л. Рабинер, Б. Гоулд. – М.: Мир, 1978. – 848 с.
7.Голд Б. Цифровая обработка сигналов: пер. с англ. / Б. Голд, Ч. Рэйдер. – М.: Советское радио, 1973. – 368 с.
8.Каппелини В. Цифровые фильтры и их применение: пер. с англ. / В. Каппелини, А.Д. Константинидис, П. Эмилиани. – М.: Энергоатомиздат, 1983. – 360 с.
9.Лэм Г. Аналоговые и цифровые фильтры. Расчет и реализация: пер.
сангл. – М.: Мир, 1982. – 592 с.
10.Алгоритмическое и программное обеспечение нерекурсивной фильтрации: отчетоНИР(промежуточный) / Новосиб. электротехн. ин-т; руководитель А.С. Анисимов. – № ГР 81067432; Инв. № 02840082120. – Новосибирск, 1984. – 35 с.
11.Алгоритмическоеипрограммноеобеспечениенепрерывнойфильтрации: отчет о НИР (промежуточный) / Новосиб. электротехн. ин-т; руководитель А.С. Анисимов. – № ГР 81067432; инв. № 02840073921. – Новосибирск, 1984. – 81 с.
90
