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

книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации

.pdf
Скачиваний:
49
Добавлен:
25.10.2023
Размер:
14.41 Mб
Скачать

Для фильтра нечетного порядка величина М одно­ значно определяется через коэффициент 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)

330

Предварительно из 2я полюсов следует отобрать ппо люсов, удовлетворяющих неравенству

l m ^ g ^ ^ > 0 ,

(4.83)

которое вытекает из условия устойчивости фильтра В 2- плоскости

| 2"^ | < 1.

Для перехода от низкочастотного к высокочастотным полосовым и режекторным фильтрам можно воспользо­ ваться преобразованиями [139]. Для фильтра верхних частот

*

=

- Ш

<4-84>

где

 

 

 

k = — cos (°'Т'

-

J cos ("° Т * + ^ .

 

Пример. Рассчитаем цифровой эллиптический фильтр верхних частот второго порядка. Исходные данные для расчета: относитель­ ная частота среза 0,6л, относительная граничная частота режекции 0,08я, коэффициент пульсации в полосе пропускания є = 1 (3 д Б ) . Составим характеристическую функцию для соответствующего филь­ тра нижних частот и подставим ее в уравнение (4.75):

' •

f

соГд

где параметры характеристической функции определяются по фор­

мулам (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 .

331

Приравнивая нулю знаменатель и решая полученное уравнение, на­ ходим полюсы фильтра в тоіі же плоскости:

( 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,

где Лі = О і , А22,

Bi=—Ь\,

Вг=—Ьг-

4.3.5. Особенности других методов синтеза цифровых фильтров. К .настоящему времени наряду с рассмотрен­ ными выше -разработано] много -различных методов син­ теза цифровых фильтров, так что один и тот же цифро­ вой фильтр может быть синтезирован различными мето­ дами. Кратко охарактеризуем некоторые из них.

332

Метод

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

частот

на г-плоскости [139]

і состоит в

следующем.

Задана передаточная

функция

К (г)

цифрового фильтра

нижних

частот.

Необходимо

найти

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

g(z),

такое,

чтобы

K[g(z)]

стала

передаточной функцией одного из следующих типов циф­ рового фильтра: верхних частот, полосового или режекторното. Кроме того, это преобразование должно сохра­ нить вид частотной характеристики прототипа.

Как отмечается в [130], нахождение аналога непре­ рывного фильтра в дискретной области при помощи z- преобразований не всегда является оптимальным реше­ нием данной проблемы, так как не все методы синтеза приводят к цифровым фильтрам, имеющим ту ж е сте­ пень устойчивости, как и их аналоговый прототип. Изве­ стно, что цифровые фильтры, синтезированные в частот­ ной области, .не обладают удовлетворительными времен­ ными характеристиками [131], однако во многих практических ситуациях временная характеристика так ж е важна, как и частотная. Поэтому синтез цифровых фильтров во временной области связан с выравниванием временных характеристик аналогового и цифрового филь­ тров в момент выборки при подаче на входы фильтров •одного я того же сигнала. Процедура синтеза цифрового 'фильтра сводится к минимизации среднеквадратической ошибки между выходными выборками аналогового и цифрового фильтров [131].

Необходимо также остановиться на вопросах синтеза ^нерекурсивных фильтров, отклик которых определяется выражением (4.34). Как правило, нерекурсивный фильтр не имеет заметных преимуществ перед рекурсивным: для его реализации требуется больше элементов задержки и умножителей, но в то ж е время не возникает проблем устойчивости. Однако нерекурсивные фильтры имеют особенности, благодаря которым они находят применение в радиолокации. В частности, имеется возможность пост­ роения параллельного спектроанализатора (гребенки фильтров) на базе одних и тех же элементов задержки, используя матрицу весовых резисторов [132], что упро­ щает реализацию.

Нерекурсивные фильтры находят применение в

Виде

алгоритмов

Э Ц В М при машинном анализе

случайных

процессов, так как определение требуемой

импульсной

переходной

характеристики осуществляется

более

про­

стыми методами, чем определение передаточной функции

333

рекурсивного фильтра. Использование алгоритмов бы­ строго преобразования Фурье (БПФ) [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]. В отличие от

334

метода -весовой функции -импульсная переходная харак­ теристика не обрывается, а происходит сложение отбро­ шенных ее компонент с оставшимися (aliasing), что при­ водит « искажению частотной характеристики. Нерекур­ сивный фильтр, синтезированный по методу частотной выборки, является составным и состоит из последова­ тельно соединенных гребенчатого фильтра с передаточ­ ной функцией K(z)=l—zm, содержащей т нулей и гре­ бенку элементарных фильтров (параллельное соединение

тфильтров первого порядка).

4.3.6.Метод анализа «нулей» процесса. Широкое распростране­ ние получили цифровые фильтры, основанные на анализе фазовых свойств или «нулей» процессов во временной области. К их числу относится цифровой фазовый фильтр (см. п. 2.5.3) и цифровые филь­

тры, основанные на корреляторах совпадения полярностей (140— 142]. По сравнению с классическими цифровыми фильтрами подобные

/7орог счета

ult) І Детектор „нулей."

Л . * J

_5_

S R.

5 R

T 3

 

1 ^F

J?

1

1

 

Счетчик I

" П V - * • Счетчик iff

 

Уст.„іТ

 

An

1 L

 

Счетчик //

R

 

S

 

Уст.,,0"

 

alt)

 

Дет.,,0"

 

6")

Рис. 4.8. Щелевой фильтр (а) и цифровой аналог щелевого

фильтра ( б ) .

335

устройства значительно проще в технической реализации, однако им присущи недостатки, отмечавшиеся в п. 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 % .

 

 

336

 

От Нестабільностей

величин Т и т/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.Вводные замечания. Д о сих пор мы рассматри­ вали вопросы синтеза и анализа цифровых фильтров на основе теории линейных разностных уравнений с постоян­ ными коэффициентами. Тем самым предполагалось, что входящие в уравнения постоянные параметры непрерыв­ ны и принимают заданные значения. При реализациях ЦФ эти параметры являются дискретными из-за огра­ ниченности разрядной сетки Э Ц В У или Э Ц В М . Входные сигналы также являются дискретными, а при умноже­ нии возникают погрешности, связанные с округлением произведений. Таким образом, шумы квантования, ошиб­ ки округления и погрешности задания коэффициентов разностного уравнения составляют три основных источ­ ника ошибок в цифровых фильтрах.

338

Анализ и учет названных ошибок возможны только в том случае, если их величина значительно меньше ди­ намического диапазона входных воздействий и значений параметров, т. е. при большом числе разрядов аналогоцифрового преобразования и соответствующей разряд­ ной сетке ЭЦВУ . В этом случае эффекты квантования входных сигналов сводятся к аддитивным шумам кванто­ вания, воздействующим на вход фильтра, эффекты округ­ ления аналогичны шумам квантования, хотя источники

находятся не на входе фильтра, а на выходах

умножите­

лей. Учет погрешностей задания

параметров

фильтра

в нерекурсивных фильтрах связан

с оценкой

изменения

его импульсной переходной характеристики. В рекурсив­ ных ж е фильтрах в этом случае полюсы передаточной функции могут оказаться вне единичного круга на плос­ кости 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

Соседние файлы в папке книги из ГПНТБ