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

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

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

Передаточные функции, флуктуационные и динамичен ские ошибки фильтров с эффективной конечной памятью, рассмотренных в § 4.6, приведены в приложении 2. Во­ просы анализа цифровых фильтров рассмотрены в п. 4.2.3.

3.3.Ц И Ф Р О В Ы Е Э К В И В А Л Е Н Т Ы

АН А Л О Г О В Ы Х И З М Е Р И Т Е Л Е Й

3.3.1.Введение. Синтез цифрового эквивалента ана­

логового измерителя состоит в следующем. Для непре-

—V —У

рывной плотности wn(u\X) синтезируется решающая

функция у (и). Она может соответствовать либо решению уравнений максимума правдоподобия (3.1), (3.2), либо уравнению оптимального дискриминатора (3.7), (3.8), либо оператору формирования апостериорной плотности

измеряемых параметров % на данном шаге п с учетом всей предшествующей информации [104]

 

 

00

оо

 

_^

^

 

wn{l\u)

= c

J

. . . j w

 

.%i)w[u\Xlt

Z

Xn)X

 

 

— 00

— 0 0

 

 

 

 

 

 

 

Х Д

- Д і - , ,

(3.43a)

с последующим

формировании

оценки

 

 

 

 

 

 

 

 

(3.436)

 

 

 

 

Ф

 

 

 

где до (А!

 

А„) — априорная

плотность

вероятности

из-

меняющихся

измеряемых

параметров; w(u\ A l f . . . , An ) —

функция

правдоподобия;

с — нормирующая постоянная.

-*•

По аналоговому прототипу у (и) синтезируется циф­ ровой эквивалент, т. е. производится замена

- V

Т(ы)=>Т(«).

В устройствах обработки радиолокационной информа­ ции преимущественное применение находят измерители двух типов. В Р Л С обзора используются многоканаль­ ные неследящие измерители. Многоканальная структура обусловлена необходимостью измерения координат боль­ шого числа целей. В Р Л С слежения используются одноканальные следящие измерители.

240

Чаще всего многоканальный измеритель является устройством совместного обнаружения — измерения и строится в соответствии с уравнениями (3.9). Многока­ нальная система может быть построена различными спо­ собами (см. п. 2.2.3). При параллельном анализе всей

области Ф определения параметров X, разбитой на М каналов, возможно использование М измерителей или сочетание одного измерителя и адресного или ассоциа­ тивного запоминающего устройства. При последователь­ ном анализе области Ф используется один измеритель. Часто применяются комбинированные последовательнопараллельные измерители [98]. В дальнейшем мы будем рассматривать преимущественно один канал многока­ нального неследящего измерителя. Для следящих изме­ рителей мы ограничимся рассмотрением цифровых ди­ скриминаторов, поскольку цифровые сглаживающие це­ пи могут быть построены по оптимальным линейным фильтрам в соответствии с (3.43), а вопросы синтеза по : следних рассмотрены в § 4.5, 4.6.

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

3.3.2. Цифровые многоканальные неследящие дально­ меры и дискриминаторы следящих дальномеров. Анало­ говыми прототипами цифрового многоканального даль­ номера в импульсной РЛ С обзора являются индикатор кругового обзора (ИКО) [63] и рециркулятор [50]. Прин­ ципы организации многоканальной структуры цифровых дальномеров рассмотрены в [98]. В каждом канале даль­ номера производится обнаружение бинарно-квантован­ ных сигналов. Обнаружитель каждого канала может быть любым: детектор движущегося окна, весовой би­ нарный, программный, накопитель позиций, либо весо­ вой вариант последовательной процедуры.

Остановимся на вопросе точности измерения дально­ сти такими дальномерами. Дальность цели определяется по номеру канала, в котором произошло обнаружение цели (выполнен критерий обнаружения). С учетом того,

16—1410

241

что распределение Истинной дальности цели относительно импульсов дальности равновероятно, среднеквадратнческая ошибка оценки времени запаздывания равна

а [ Г ] ^ 7 д / 2 ] / " 3

(при Г д ^ т и ) ,

(3.44)

где Тп длительность зондирующего импульса.

В первом приближении

о [Г] не зависит от отношения

сигнал/шум q°-. Для пояснения этого факта

зафиксируем

q2 и будем уменьшать интервал дискретности по дально­ сти Гд. Тогда при Г д < т ц вследствие флуктуации времени запаздывания некоторые импульсы пачки будут попадать в смежные дальномерные каналы, что приведет к сниже­ нию вероятности правильного обнаружения цели. До тех пор пока можно пренебречь вероятностями ложных об­ наружений в смежных каналах, точность измерения вре­ мени запаздывания будет соответствовать (3.44). Нали­ чие ложных обнаружений данной цели в смежных кана­ лах увеличивает ошибку измерения времени запаздыва­ ния. Как правило, в подобных системах не предъявляет­ ся высоких требований к точности дальнометрии, по­

скольку ошибка

в определении местоположения

цели

в большей степени зависит от точности измерения

ази­

мута >8ц и угла

места 0Ц .

 

Перейдем теперь к вопросу построения цифровых ди­ скриминаторов следящих дальномеров. Уравнение опти­ мального дискриминатора имеет вид (3.8). На практике

находят применение

дискриминаторы типа

 

 

 

X =

K-~L^\l)

 

= Xu-F{s)^

 

(3.45)

где Ь(и\Х) —информационный

эквивалент

условного ко­

эффициента

правдоподобия А(и\Х) в общем

случае не­

оптимальный

[см.

(3.9)].

Образование

производной

в (3.45) заменяется формированием конечной

разности:

F[B) ^

L(u*\T-^)-L{u\T

+ ^) _

 

( 3

 

 

 

 

->

 

 

 

Чаще всего статистика Ь(и\Т)

является

корреляцион­

ным интегралом

 

 

 

 

 

 

 

L(u\T)=

j

v(t—T)«оц(t)dt,

 

(3.47)

242

где uon(t)—опорное напряжение, согласованное по фор­ ме с излучаемым сигналом. Структурная схема дискри­ минатора, соответствующего уравнениям (3.46). (3.47), хорошо известна, она приведена на рис. 3.7. Функцию умножителей выполняют стробируемые усилители (СУ).

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

аналого-цифрового

преобразователя

«временной

интер­

в а л — код», не

включенного в

замкнутый контур

следя­

щей системы.

Если

сигнальный

вход

будет аналоговым,

то и вход опорного сигнала также должен быть

аналого-

 

 

СУ

ґ

 

 

 

 

 

v(t-7)

J

 

+

eft)

 

 

 

E

 

СУ

" L

V, (t-T-1)

Рис. 3.7. Структурная схема дискриминатора задержки.

вым, что потребует включения цифро-аналогового преоб­ разователя «код — временной интервал» в цепь обратной связи. Перемножители и интеграторы в этом случае вы­

полняются

на цифровых элементах.

 

 

На рис. 3.8 представлена структурная схема цифро­

вого аналога дискриминатора рис. 3.7.

Перемножение

принятого импульса v(tТ)

и двух опорных

Vi(tТ—та),

vz(tТ+Ти)

производится

в двух схемах

И с

заполнени­

ем временных интервалов импульсами масштабной ча­ стоты. Период Тл=Хц/п, где м ^ 8 . Функции интеграторов и вычитающего устройства выполняет реверсивный счет­

чик,

формирующий

сигнал

рассогласования. В

на­

чале

каждого зондирования

счетчик устанавливается

на

0. Для уменьшения

ошибок

квантования в дискримина-

16*

243

торе следует уменьшать период следования масштабных импульсов 7"д, но это требует увеличения быстродействия реверсивного счетчика. Ширина дискриминационной ха­ рактеристики равна 2тц .

 

Более простои в реализации является схема, пред­

ставленная на рис. 3.9. Короткий

импульс v\(tТ—ти/2)

длительностью т/=т„//г с выхода преобразователя

«код-—

временной

интервал» подается на

вход линии

задержки

т п

с .п

отводами. Импульс

цели

с

выхода

приемника

v(t—Т)

длительностью

Тп поступает

на

входы

схем

И. В зависимости от величины временного

рассогласова­

ния

между

отраженным

и

опорным

импульсом

будет

 

 

 

 

 

г

 

 

 

 

 

 

 

 

 

 

_П_

 

 

 

 

 

 

 

t>(t-T)

 

Рев.

счетчик:

Уст. 0

 

 

 

 

 

 

 

 

 

 

 

г

ЧУ

 

 

Efk)

 

 

 

 

 

 

 

I

v,(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 3.8. Структурная схема цифрового аналога дискриминатора задержки.

сформирован код, запоминаемый в регистре. В частно­ сти, при отсутствии рассогласования в регистре будет код (1, 1, . .. ., 1, 1). Код рассогласования, записанный в ре­ гистре, с помощью преобразователя преобразуется в двоичный код. В отличие от схемы рис. 3.8 здесь не тре­ буется быстродействующий реверсивный счетчик. Шири­ на дискриминационной характеристики равна ти -

3.3.3. Цифровые многоканальные неследящие частото­ меры и дискриминаторы следящих частотомеров. В им- пульсно-допплеровских Р Л С и непрерывных Р Л С для измерения допплеровской частоты используются набо­ ры (гребенки) полосовых фильтров [63], перекрывающих диапазон допплеровских частот /вшах—fomin- Число филь­

тров M=(fDmax—fuminJ/Aif, где Af — полоса пропускания

244

фильтра. Обычно М составляет несколько сотен. Чаще всего применяются электромеханические фильтры, осно­ ванные на магнитострикционном эффекте .[115]. Построе­ ние цифрового эквивалента подобной гребенки в принци­ пе возможно на основе использования как рекурсивных, так и нерекурсивных цифровых фильтров. Разрядность

v,tt-r-z/2)_

1

и

|

 

oft-Г)

1

1

 

 

1 *

* *

 

 

Регистр

 

Уст О

 

 

 

Преобра

зователь

е(к)

 

кода

 

 

Рис. 3.9. Структурная схема цифрового дискриминатора задержки.

аналого-цифрового преобразования должна быть не ме­ нее 8—12 разрядов.

Пусть каждый фильтр соответствует элементарному звену второго порядка последовательной канонической формы с передаточной функцией /^ - (z - 1 ) (4.21) рис. 4.4,6. Техническая реализация такого фильтра требует наличия

-j(k) і elf)

v,(k)

Рис. ЗЛО. Цифровой дискриминатор системы частотной автопод­ стройки.

четырех умножителей, двух регистров сдвига и двух трехвходовых сумматоров. Д а ж е при выполнении опера­ ций в последовательном коде объем оборудования зна­ чителен. Цифровые гребенки фильтров могут быть кон-

245

курентно-способными по отношению к современным ана­ логовым устройствам только при использовании больших интегральных схем.

Рассмотрим принципы построения цифрового дискри­ минатора системы частотной автоподстройки [116]. Его структурная схема приведена на рис. 3.10. Оба входа дискриминатора — сигнальный и опорного напряжения— цифровые. Дискриминатор содержит два перестраивае­ мых фильтра с передаточными функциями Ki(z) и /Сг(Х),

два квадратора, сумматор и устройство

усреднения 2,

с выхода которого и снимается сигнал

рассогласования

Рис. 3.11. Полюсы и нули частотного дискриминатора (а) и дискри­ минационная характеристика цифрового частотного дискримина­ тора (б) .

e(f)

=Д1/(£). Передаточные функции фильтров таковы, что

Ki(z)

Кг(г) — K(z),

где

K(z)—передаточная

функция

восьмиполюсного

цифрового

фильтра

Баттерворта

[5].

Оба

фильтра

настроены

на

близкие

частоты так,

что

Ki(z)

имеет

резонансную

частоту ниже

опорной

ю0 ш

a Кг(г)—выше

 

опорной

ш0п- Расположение полюсов и

нулей на 2-плоскости изображено на рис. 3.11,а. Четыре нуля расположены так, чтобы сделать фильтр с большей

крутизной

склонов частотной

характеристики.

Передаточные функции обоих фильтров имеют вид

К

Ы

 

 

K{\-2z-*

+ z-*)

4

, v

'

(1 — Л . г - і

B , z - 2 ) ( 1 —

A2z-i—B2z~*y

K

Ы

_

/ С ( 1 - 2 г - 2 + 2 - * )

 

2 v

' ~ ~ (1 — A3z-і

B 3 z - 2 ) (1 -

Л ^ - 1 — B 4 z - 2 ) '

246

где

 

 

 

 

 

 

 

i4rf =

2exp(Pf 7"B )cos(Qi 7, „);

 

 

 

 

В І =

exp ( - 2Р,;7д);

P, = -

(Дш/2) sin (ir/8);

 

Р 0 =

- ( Д ш / 2 )

cos (тс/8); Л = Р 4 ,

^ =

Л ;

 

Q , = c o o a + P 2 ;

 

P 2 = (DO U +

P I ;

 

 

Q3 =

t o ou Л ;

 

Q.i=«>0 u — pi\

 

 

До>ширина

дискриминационной

характеристики;

/ ( — постоянная,

не зависящая

от won (при шоп ^Дсо), ее

величина выбирается из условия, чтобы

K\(z)Kz(z)

обес­

печивало усиление на частоте со0п, равное единице.

 

Коэффициенты ВІ

зависят

только

от Дсо ширины

дискриминационной характеристики и вычисляются за­ ранее. Центральную частоту дискриминационной харак­

теристики

можно

подстраивать, изменяя коэффициенты

А і, Аг,

Л з,

Л4 , для чего требуется вычислять четыре ко­

синуса

и четыре

произведения.

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

З.ЗЛ. Цифровой фазовый дискриминатор. В аналого­ вых измерителях функции фазового дискриминатора обычно выполняет фазовый детектор. Для построения цифрового эквивалентна фазового детектора имеется, по крайней мере, две возможности. Во-первых, можно проквантовать выходное напряжение аналогового фазового детектора. Подобная схема используется в цифровых измерителях, когда частоты входного колебания и опор­ ного напряжения достаточно высоки. В этом случае ис­ пользуются цифровые сглаживающие цепи и в цепи об­

ратной связи включается

преобразователь код — часто­

т а — синтезатор

частот.

Во-вторых, можно построить

цифровой аналог фазового дискриминатора, преобразо­

вав в цифровую форму отсчеты

входного колебания.

Опорное колебание также должно

быть

представлено

цифровым кодом.

 

 

 

Дискриминационная характеристика

фазового

детек­

тора имеет вид

 

 

 

Я д [ ' в ( 0 ї = ( 1 / 2 ) и ( * ) О і ( * ) c o s f o t a W - q p o n t f ) ] ,

(3.48)

247

-где Овх, £>опамплитуды входного и опорного напряже­ ний; (рвх, фоп фазы входного и опорного напряжений. Соответственно дискриминационная характеристика циф­ рового фазового дискриминатора

^д[е(/г)]= (1/2т с о з { ф в х ( й ) - ф о п ( / г ) ] .

(3.49)

Структурная схема, соответствующая уравнению (3.49), приведена на рис. 3.12. В качестве усредняющего устрой­ ства 2 удобно использовать простой в реализации опера­ тор экспоненциального сглаживания. Подобные дискри­ минаторы могут использоваться на не очень высоких частотах входного и опорного напряжений. Ограничения

vjt,y>sJt)]

£(*)

X

Рис. 3.12. Структурная схема цифрового фазового дискриминатора.

связаны

с быстродействием преобразователя напряже­

н и е — код, умножителя и усредняющего устройства.

3.3.5.

Цифровые многоканальные неследящие угломе­

ры. В импульсных РЛС обзора измерение азимута про­ изводится цифровыми устройствами первичной обработ­ ки. Вся проблема углометрии может быть разделена на ряд задач: а) синтез алгоритма углометрии и определе­

ние

его параметров, б) оценка точности углометрии,

в)

техническая реализация алгоритма.

В плане синтеза данная задача соответствует оценке параметра, постоянного на интервале наблюдения и она тесно связана с задачами обнаружения некогерентных сигналов (см. § 2.3). Статистические характеристики шу­ мов и эхо-сигналов приведены в п. 2.3.2.

В соответствии с (3.9) вместо максимизации по из­ меряемому параметру функции правдоподобия или ее информационного эквивалента можно максимизировать коэффициент правдоподобия или его информационный эквивалент. Поэтому для определения алгоритма опти­ мальной углометрии воспользуемся соотношением (2.58). С учетом (3.9) имеем

Ыр)=

=Е'^(р)-,>с. (3.50)

/=о

1=0

248

За оценку принимается значение к, обеспечивающее max Ук($). В этом случае принимается р=&.

Фиксация момента времени, когда сумма (3.50) до­ стигает максимума максиморума, неудобна, поскольку

Считывание

Рис. 3.13. Преобразователь угла поворота в код.

неизвестна его величина. Удобнее определять момент прохождения через нуль выражения

M =

2 £ G . ( R

) ^

V k . ^ W ' ^

V k ^ Q

(3.51)

 

/ = 0

 

/=0

 

 

при

выполнении

условия

(2.58).

 

 

Входящая в формулу (3.51) функция <?л(|3) опреде­

ляется в соответствии с (2.41). Таким образом,

алгоритм

оптимальной углометрии

(определение

центра пачки не-

флуктуирующих сигналов) при малом отношении сиг­ нал/шум состоит в весовом скользящем суммировании выходных напряжений квадратичного детектора и опре­ делении момента времени, соответствующего максимуму суммы (3.50), либо в фиксации момента прохождения взвешенной суммой через нулевой уровень. Импульс от­ счета азимута подается в схему преобразователя угла поворота в цифровой код — рис. 3.13. Данный преобразо­ ватель основан на методе считывания И' работает сле­ дующим образом. На оси вращения антенны закреплен диск (кодовая маска), на котором фотоспособом нанесе­ на двоичная кодовая комбинация прозрачных и непро­ зрачных участков. Каждому дискретному углу поворота

?49

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