Для фильтра нечетного порядка величина М одно значно определяется через коэффициент N.
Молено показать, что, исходя из оптимальности функ ций (4.76) и (4.77), существует связь нулей и полюсов характеристическоп функции:
aiPi = a 2 p 2 = . . . =amP,„ = ,v0 ) |
(4.80) |
где в случае цифрового фильтра
Так, для фильтра |
второго порядка, |
решая совместно |
(4.78) и (4.80), |
получаем |
|
|
|
|
(4.81) |
Используя |
тот |
факт, что функция |
Р ( ы Т д ) — 1 для |
фильтра третьего порядка имеет два корня в полосе про пускания, находим дополнительное уравнение
решая которое совместно с (4.79) и |
(4.80), определим |
нули н |
полюсы характеристической |
функции фильтра |
третьего |
порядка. Нули и полюсы |
характеристических |
функций цифровых эллиптических фильтров высших по
рядков вычисляются |
по |
аналогии |
с |
определением |
соот |
ветствующих величин для непрерывных фильтров. |
|
|
Подставляя найденную |
характеристическую функцию |
в |
уравнение |
(4.75), |
молено |
опргделить |
нули ^tg - ^ p j |
и |
полюсы |
( t g ш 9 д ) |
квадоата |
амплитудно-частотной |
характеристики в комплексной плоскости tg" |
Пере |
ход к полюсам и нулям |
в |
плоскости |
z |
осуществляется |
при помощи |
соотношения |
[139] |
|
|
|
|
|
* = ( |
l |
+ < ( « |
^ ) |
) / |
( l - « |
( t |
g |
^ ) ) . |
(4.82) |
Предварительно из 2я полюсов следует отобрать ппо люсов, удовлетворяющих неравенству
которое вытекает из условия устойчивости фильтра В 2- плоскости
| 2"^ | < 1.
Для перехода от низкочастотного к высокочастотным полосовым и режекторным фильтрам можно воспользо ваться преобразованиями [139]. Для фильтра верхних частот
* |
= |
- Ш |
<4-84> |
где |
|
|
|
k = — cos (°'Т' |
- |
J cos ("° Т * + ^ . |
|
Пример. Рассчитаем цифровой эллиптический фильтр верхних частот второго порядка. Исходные данные для расчета: относитель ная частота среза 0,6л, относительная граничная частота режекции 0,08я, коэффициент пульсации в полосе пропускания є = 1 (3 д Б ) . Составим характеристическую функцию для соответствующего филь тра нижних частот и подставим ее в уравнение (4.75):
где параметры характеристической функции определяются по фор
мулам (4.81) |
при |
заданных •соо7'д=0,6я и |
шС Г П = 0 , О 8 л : а = 0 , 7 0 8 , |
Р = |
1-5,4, М=472,9. |
|
|
|
Корни числителя |
в найденном выражении |
определяют нули филь- |
тра |
в плоскости |
tg—g—• |
|
( t g ^ ) B = ± p t £ ^ = ± l , 9 4 5 .
Приравнивая нулю знаменатель и решая полученное уравнение, на ходим полюсы фильтра в тоіі же плоскости:
( t g ^ ) = ± a t g ^ |
- j / " |
T^v+M±rm=m-. |
+ a - ^ 2 t g — L e |
= + (0,0407 + /0,0983). |
Нули и полюсы низкочастотного фильтра z " и z" вычислим в со ответствии с (4.82), предварительно выбрав полюсы в плоскости
tg — , — с положительной мнимой частью. Для перехода к фильтру
верхних частот воспользуемся преобразованием (4.84). Окончательно получим
2" = - f q f l F = е ± Ш ' ' 3 ' |
= 0,9841 + «0,1774, |
2° = - JTf^T = 0 , 7 0 2 е ± Ш 8 ° 5 |
8 ' = - 0,3400 + (0,6141. |
От передаточной функции найденного фильтра
к, л _ AM — ( 2 - - п ) ( г - г - * ) _
А( г ) ~ В (г) ( г — 2 п ) ( г — г"*)
|
_ |
1 — ( 2 R e z " ) z - ' + z - 2 |
|
|
1 — (2 Re 2 n ) z - I + | z n | 2 z " 2 |
можно перейти к разностному уравнению |
|
|
2 |
2 |
v (kTA) = |
£ |
atu f(ft - |
0 7-„] + S bp [ (k - і) Гд], |
|
/=0 |
і= 1 |
коэффициенты которого равны |
a 0 = f l i = l , а і = — 2 R e z H = — 1 , 9 6 8 2 , b\ = |
= 2Rez»= — 0,6400, |
b2=—0,4928. |
|
Всоответствии с разностным уравнением структурная схема
цифрового |
фильтра принимает |
вид, аналогичный изображенному на |
рис. 4.4,6, |
где Лі = О і , А2=а2, |
Bi=—Ь\, |
Вг=—Ьг- |
4.3.5. Особенности других методов синтеза цифровых фильтров. К .настоящему времени наряду с рассмотрен ными выше -разработано] много -различных методов син теза цифровых фильтров, так что один и тот же цифро вой фильтр может быть синтезирован различными мето дами. Кратко охарактеризуем некоторые из них.
|
|
|
|
|
|
|
|
Метод |
преобразования |
частот |
на г-плоскости [139] |
і состоит в |
следующем. |
Задана передаточная |
функция |
К (г) |
цифрового фильтра |
нижних |
частот. |
Необходимо |
найти |
преобразование |
g(z), |
такое, |
чтобы |
K[g(z)] |
стала |
передаточной функцией одного из следующих типов циф рового фильтра: верхних частот, полосового или режекторното. Кроме того, это преобразование должно сохра нить вид частотной характеристики прототипа.
Как отмечается в [130], нахождение аналога непре рывного фильтра в дискретной области при помощи z- преобразований не всегда является оптимальным реше нием данной проблемы, так как не все методы синтеза приводят к цифровым фильтрам, имеющим ту ж е сте пень устойчивости, как и их аналоговый прототип. Изве стно, что цифровые фильтры, синтезированные в частот ной области, .не обладают удовлетворительными времен ными характеристиками [131], однако во многих практических ситуациях временная характеристика так ж е важна, как и частотная. Поэтому синтез цифровых фильтров во временной области связан с выравниванием временных характеристик аналогового и цифрового филь тров в момент выборки при подаче на входы фильтров •одного я того же сигнала. Процедура синтеза цифрового 'фильтра сводится к минимизации среднеквадратической ошибки между выходными выборками аналогового и цифрового фильтров [131].
Необходимо также остановиться на вопросах синтеза ^нерекурсивных фильтров, отклик которых определяется выражением (4.34). Как правило, нерекурсивный фильтр не имеет заметных преимуществ перед рекурсивным: для его реализации требуется больше элементов задержки и умножителей, но в то ж е время не возникает проблем устойчивости. Однако нерекурсивные фильтры имеют особенности, благодаря которым они находят применение в радиолокации. В частности, имеется возможность пост роения параллельного спектроанализатора (гребенки фильтров) на базе одних и тех же элементов задержки, используя матрицу весовых резисторов [132], что упро щает реализацию.
Нерекурсивные фильтры находят применение в |
Виде |
алгоритмов |
Э Ц В М при машинном анализе |
случайных |
процессов, так как определение требуемой |
импульсной |
переходной |
характеристики осуществляется |
более |
про |
стыми методами, чем определение передаточной функции
рекурсивного фильтра. Использование алгоритмов бы строго преобразования Фурье (БПФ) [122] ускоряет процедуру вычисления свертки (4.34).
Нерекурсивные фильтры имеют конечную импульсную переходную характеристику, а их передаточные функции содержат только нули <и не имеют полюсов на г-плоско- сти. Конечность импульсной переходной характеристики приводит к неравномерности частотной характеристики. Эту неравномерность ослабляют введением специальной весовой функции w(kTjj) с передаточной функцией w(z). Тогда передаточная функция сглаженного фильтра запи
шется |
в соответствии с теоремой 'Комплексной |
свертки |
(4.13а) |
как |
z-преобразоваиие произведения |
/г(йТд )Х |
Xw(kTx), |
т. |
е. |
|
i*l=i
Выбирая параметры весовой функции, можно ослабить неравномерность частотной характеристики расширяя по лосы пропускания. Наиболее часто используются реали зуемые аппроксимации весовой функции Дольф — Чебы шева [133]. Синтез нерекурсивного фильтра методом ве совой функции сводится к определению коэффициентов h(kTs), k = 0, 1, ... , уравнения (4.34) и зависит от того, как задан его прототип. Если задана передаточная функ ция аналогового фильтра K(S), то с помощью обратного преобразования Лапласа определяется соответствующая импульсная переходная характеристика h (t):
|
|
|
h(t) |
= |
~§K(S)e"dS. |
|
|
Далее берутся |
равноотстоящие отсчеты от h(t) |
:}і(кТр)у |
£ = 0 , 1 , 2 , . . . |
|
|
|
|
|
|
Если |
же |
задана |
передаточная функция |
цифрового |
фильтра, |
то |
соответствующая |
импульсная |
переходная |
характеристика |
вычисляется в |
соответствии |
с (4.12а). |
Далее, |
предварительные коэффициенты, |
полученные |
одним из двух методов, обрываются и оставшиеся умно жаются на весовую функцию. Весовая функция Хэмминга имеет частотную характеристику
да ((D) == 0,08 - j - 0,92 cos2 ~
Другим методом синтеза нерекурсивных фильтров является метод частотной выборки [134]. В отличие от
метода -весовой функции -импульсная переходная харак теристика не обрывается, а происходит сложение отбро шенных ее компонент с оставшимися (aliasing), что при водит « искажению частотной характеристики. Нерекур сивный фильтр, синтезированный по методу частотной выборки, является составным и состоит из последова тельно соединенных гребенчатого фильтра с передаточ ной функцией K(z)=l—zm, содержащей т нулей и гре бенку элементарных фильтров (параллельное соединение
тфильтров первого порядка).
4.3.6.Метод анализа «нулей» процесса. Широкое распростране ние получили цифровые фильтры, основанные на анализе фазовых свойств или «нулей» процессов во временной области. К их числу относится цифровой фазовый фильтр (см. п. 2.5.3) и цифровые филь
тры, основанные на корреляторах совпадения полярностей (140— 142]. По сравнению с классическими цифровыми фильтрами подобные
/7орог счета
ult) І Детектор „нулей."
Л . * J
_5_
S R.
5 R
T 3
Счетчик I |
" П V - * • Счетчик iff |
|
Уст.„іТ |
|
An
1 L |
|
Счетчик // |
R |
|
S |
|
Уст.,,0" |
|
alt) |
|
Дет.,,0" |
|
6") |
Рис. 4.8. Щелевой фильтр (а) и цифровой аналог щелевого
фильтра ( б ) .
устройства значительно проще в технической реализации, однако им присущи недостатки, отмечавшиеся в п. 2.5.3.
|
|
|
|
|
|
|
Рассмотрим ЦФ, изображенный на рис. 4.8,0, |
с помощью |
которо |
го может осуществляться полосовая |
фильтрация. Такой |
фильтр |
имеет |
гребенчатую |
частотную |
зависимость с |
полосами! |
прозрач |
ности |
на частотах |
fj £ = /cf«^-/.'/'/', |
k—\, 2 |
где Т — задержка |
в эле |
менте. Крутизна фронтов частотной зависимости теоретически пря
моугольна и определяется |
крутизной |
фронтов импульсов |
формиро |
вателя, |
стабильностью срабатывания |
схемы И и стабильностью сра |
батывания счетчика от импульсов малой длительности. |
|
С |
помощью детектора |
«нулей» и |
формирователя F |
образуются |
прямоугольные импульсы длительностью т/2, соответствующие мо
менту пересечения входным напряжением нулевого |
уровня |
снизу |
вверх. После электронной задержки |
(например, с помощью кипп-реле) |
на время Т задержанные и незадержанные |
импульсы |
селектируются |
с, помощью схемы И. Счетчик регистрирует |
число совпавших импуль |
сов за время наблюдения Т„. Если |
это число |
превышает |
определен |
ную |
величину, |
то выносится |
решение о |
наличии сигнала |
с частотой |
/о. Отметим, что свойства фильтра будут иными, если |
задержка |
осу |
ществляется |
с помощью |
личин задержки. |
|
|
|
|
|
|
|
|
Огибающая «зубьев» |
частотной характеристики с ростом / убы |
вает |
по закону |
k=\, |
2, |
.. Ширина полосы |
пропускания каждого |
зуба |
равна |
(при т < 7 " ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
к |
|
к |
|
ki |
|
г |
|
|
|
Д/Іі = |
/щах |
' m l i i |
~ |
Т |
і |
Т |
Т |
(Т |
х) ^ |
^ Т*~' |
(^-^) |
а добротность |
не зависит |
от к: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Qh4h!Afu |
= Tlx. |
|
' |
|
|
(4.86) |
При дополнении фильтра схемой НЕ, обозначенной |
на |
рис. 4.8,а |
штриховыми линиями, фильтр превращается в режекторный и |
вме |
сто полос прозрачности появляются полосы затухания. |
|
|
Если данный фильтр используется как одпополосный, то фильт |
рацию гармоник частоты /о можно |
осуществить |
дополнительным |
входным фильтром с полосой AF, выбираемой из условия |
Afi^AF< |
<2/о—A/i/2, |
и |
центральной |
частотой /0 , |
при |
этом |
к |
форме |
его |
амплитудно-частотной характеристики не |
предъявляется |
высоких |
требований. |
Необходимо, |
чтобы |
его |
фазо-частотная |
характеристика |
была линейна в пределах полосы Д/і. |
|
|
|
|
|
|
|
• Рассмотренный |
фильтр |
имеет |
одну |
особенность, |
ограничиваю |
щую |
его применение. |
При наличии в спектре |
входного |
процесса ко |
лебаний с частотами, не попадающими в полосу i4.fi, но имеющими
амплитуды, |
большие, |
чем полезный |
сигнал, статистика |
«нулей» |
входного процесса не |
будет соответствовать |
статистике |
«нулей» |
фильтруемого |
сигнала. |
Происходит |
подавление |
слабого |
сигнала |
сильным. С помощью дополнительного полосового фильтра, уста навливаемого на входе данного фильтра, этот эффект можно осла бить.
Из-за аппаратурных погрешностей (нестабильности Т и т/2) частотная характеристика фильтра ухудшается. Величина нестабиль
ности |
должна быть такой, чтобы AT<xjn, |
где п>\, |
ATIT<i/nQ. |
Так, |
при / о = 3 0 0 кГц, Г = 3 , 3 3 мке, т / 2 = 5 |
ис, іД/, = 1 кГц, |
0=330, |
я = 2 , Д 7 7 Г < 0 , 1 5 % . |
|
|
|
От Нестабільностей |
величин Т и т/2 можно |
избавиться, форми |
руя |
|
соответствующие |
импульсы |
с |
помощью |
цифровых |
схем |
(рис. 4.8,6). Счетчик I, |
заменяющий |
схему |
электронной |
задержки |
D, |
выполняет |
функции |
.преобразования кода во временной интер |
вал, счетчик II используется для формирования импульса т. Теперь |
Т=пТ3, |
х=тТз, |
|
где Т3—период |
|
следования |
импульсов |
генератора |
меток; |
/1=1, 2,..., |
т = 1 , |
2 |
/ ; > т . |
Последовательности |
импуль |
сов |
с |
выходов |
1 |
и 2 генератора меток сдвинуты |
относительно |
друг |
друга на величину Г3 /2. |
Входной |
процесс |
поступает на детектор |
«нулей». Импульс |
«нуля» |
обрасы- |
|
|
|
|
|
|
состоящий из триггера Тг и схемы Иг. Импульсы с выхода 2 гене ратора меток, следующие с перио
дом Тз, заполняют |
счетчик II и от |
крывают ключ 3 |
(Т 3 и Из) . Сиг |
нал |
переполнения |
счетчика |
I I за |
крывает ключ |
сбрасывает счет |
чик |
I и открывает ключ / |
(Ті и |
И і ) , |
после чего импульсы |
с выхо |
да 1 генератора меток заполняют счетчик I. Если период следова ния входного колебания совпадает ао временем заполнения счетчика I, то на выходе схемы Из появ ляется единица. С помощью счет чика I I I производится накопление
~ |
і |
|
0,75 |
Aft |
V |
|
0,5 |
' У |
0,25 |
|
\ |
|
|
' тіл. |
|
f |
Рис. 4.9. Частотная характери стика.
единиц. Если за время Тп накапливается N^kB единиц, где k0 — порог счета, то считается, что присутствует выделяемый сигнал.
Введением в счетчик I I I кода |
Ani=Miu—fto+'l, |
где М ш — емкость |
счетчика |
I I I , обеспечивается |
установка заданного порога счета ка- |
Точно так же |
производится |
установка |
длительностей |
Г = / г Г 3 и |
х=1пТя |
путем |
введения в |
счегчики I и |
II |
кодов Ai=Mi—п + 1, |
і4ц = Л І ц — m + l . |
|
|
|
|
|
Если |
время |
прихода сигнала неизвестно, |
то вместо |
счетчика I I I |
можно использовать устройство, реализующее оператор экспоненци ального сглаживания [4], либо регистр сдвига, вход и выход которого соединены соответственно с суммирующим и вычитающим входами реверсивного счетчика (см. рис. 2.12).
Вследствие случайного фазового сдвига между частотами гене ратора меток и входного колебания появление единицы на выходе схемы Из в некоторых областях частот носит случайный характер. Плотность вероятности этого события в функции частоты wi(f) изо бражена на рис. 4.9. При наличии на выходе ЦФ накопителя (счет чик I I I ) зависимость Wy(f) на рис. 4.9 можно интерпретировать как частотную характеристику. Легко убедиться, что
1 |
|
|
|
f,m m , — T3(n + |
m) |
T3 |
( , Z - 1 + 772)' |
|
|
|
(4.87) |
fmax, — 7"3 |
(/i—1) |
• fmaxa — |
T3n |
22—1410 |
|
|
337 |
Из (4.87) можно .получить соотношения для следующих .пара метров ЦФ.
— центральная частота
|
|
. |
|
Апах3 |
+ |
|
/ІІІ:П„ |
|
|
2» — 1 -f- т |
|
|
|
|
|
' |
0 |
= |
|
2 |
|
|
" |
|
2 7 > ( м _ 1 |
- j - |
in) |
' |
— |
полоса |
пропускания |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, r |
= f |
|
_ |
|
г |
|
|
|
|
'JLzJ |
|
|
H I ) ' |
|
|
|
|
|
'mix, |
|
'rniii, |
|
Г3і\ |
(n — |
I -f- |
|
— |
полная |
полоса пропускання |
|
|
|
|
|
|
|
|
|
|
. |
„ |
г |
|
|
f |
|
|
|
|
|
m + |
1 |
|
|
|
|
|
a |
l |
— 'max, |
|
|
'mill, |
- |
|
T3 (ll —!)(/) |
+ |
Ill) ' |
— коэффициент перекрытия |
|
|
|
|
|
|
|
|
|
|
|
|
|
С ~ |
'maxj/'min, |
' -Ь |
/г |
» |
|
|
|
— коэффициент прямоугольное™ |
|
|
|
|
|
|
|
|
|
г |
_ |
, „ |
f |
Jf, |
- |
я (" |
— |
1 + |
/и) (я + |
ш) |
• |
|
|
t\ — ц |
|
|
^ |
_ |
|
(Я |
_ |
(/ я |
+ |
/ ; ) |
В заключение отметим, что такие ЦФ могут применяться для измерения частоты монохроматических и квазпмонохроматических сигналов и выделения их из шумов с гладким спектром. Основное достоинство подобных ЦФ — простота технической реализации, стабильность частотной зависимости, простота перестройки /о и Д/і и возможность получения высоких значений коэффициента .прямо- \толыюстн.
4.4. ЭФФЕКТЫ КВАНТОВАНИЯ
ВЦ И Ф Р О В Ы Х Ф И Л Ь Т Р А Х
4.4.1.Вводные замечания. Д о сих пор мы рассматри вали вопросы синтеза и анализа цифровых фильтров на основе теории линейных разностных уравнений с постоян ными коэффициентами. Тем самым предполагалось, что входящие в уравнения постоянные параметры непрерыв ны и принимают заданные значения. При реализациях ЦФ эти параметры являются дискретными из-за огра ниченности разрядной сетки Э Ц В У или Э Ц В М . Входные сигналы также являются дискретными, а при умноже нии возникают погрешности, связанные с округлением произведений. Таким образом, шумы квантования, ошиб ки округления и погрешности задания коэффициентов разностного уравнения составляют три основных источ ника ошибок в цифровых фильтрах.
Анализ и учет названных ошибок возможны только в том случае, если их величина значительно меньше ди намического диапазона входных воздействий и значений параметров, т. е. при большом числе разрядов аналогоцифрового преобразования и соответствующей разряд ной сетке ЭЦВУ . В этом случае эффекты квантования входных сигналов сводятся к аддитивным шумам кванто вания, воздействующим на вход фильтра, эффекты округ ления аналогичны шумам квантования, хотя источники
находятся не на входе фильтра, а на выходах |
умножите |
лей. Учет погрешностей задания |
параметров |
фильтра |
в нерекурсивных фильтрах связан |
с оценкой |
изменения |
его импульсной переходной характеристики. В рекурсив ных ж е фильтрах в этом случае полюсы передаточной функции могут оказаться вне единичного круга на плос кости 2, что сделает фильтр неустойчивым.
Как показывают исследования [143], прямая форма (см. рис. 4.2) и каноническая форма (см. рис. 4.3) весьма чувствительны в смысле устойчивости к точности задания коэффициентов разностного уравнения и на практике они почти не используются. Преимущественное применение находят последовательная (см. рис. 4.4) и параллельная (см. рис. 4.5) канонические формы, менее чувствительные
кточности задания их параметров.
4.4.2.Шумы квантования. При квантовании входного воздействия в аналого-цифровом преобразователе с ха рактеристикой вида рис. В.1 квантованный сигнал u(t)
можно |
представить |
как |
разность |
исходного |
сигнала |
u(t) |
и ошибки округления n(t) — ш у м а |
квантования |
|
|
|
|
U(t)=u(t)—n(t). |
|
|
|
|
|
(4.88) |
Как показывают |
исследования [4], при &и/а< |
\ (сг2 — дис |
персия |
входного |
воздействия) среднее |
значение |
шумов |
квантования равно /?іі[/г]=0, а дисперсия — |
|
|
|
|
|
|
|
|
оЩ=Аи*/\2, |
|
|
|
|
(4.89) |
поскольку |
шум равномерно |
распределен |
в |
интервале |
квантования 1А-и. |
|
|
|
|
|
|
|
|
|
|
|
Коэффициент корреляции шумов квантования пропор |
ционален |
коэффициенту |
корреляции |
|
входного |
сигнала |
Ru |
и отношению А]и/а. |
Д а ж е |
при А«/сг=1 |
шумы |
кванто |
вания практически иекоррелированы, |
если |
|
Ru<0,9.-Кро |
ме того, при Дм/ог<1 |
взаимная корреляция-сигнала |
u(t) |
и шума |
n(t) практически |
отсутствует. |
Поэтому |
учет |
21* |
|
|
|
|
|
|
|
|
|
|
|
|
339 |