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

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

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

хи, компоненты которого определяются из (2.216) и (2.217).

2.7.2. Цифровые программные обнаружители. На прак­ тике находят применение обнаружители, функционирую­ щие в соответствии с алгоритмом: начало пачки фикси­ руется по наличию к единиц в т позициях (или k еди­ ниц в к позициях), конец пачки — по наличию / нулей подряд. При этом критерий начала пачки является одно­ временно и критерием ее обнаружения, а конец пачки фиксируется с целью измерения азимута обнаруженной цели [3]. Такие обнаружители называются цифровыми программными и обозначаются «к/т—/» или «к/к—Ь. В состав цифрового программного обнаружителя входит цифровой автомат, реализующий критерий обнаружения пачки эхо-сигналов и цифровой накопитель, подсчиты­ вающий позиции внутри обнаруженной пачки и осу­ ществляющий установку в нулевое состояние автомата при выполнении критерия конца пачки.

Анализ цифровых программных обнаружителей в ус­ ловиях воздействия некоррелированных помех рассмо­ трен в [3]. Представляет интерес оценка эффективности таких устройств при воздействии коррелированных по­ мех. Составим переходные матрицы цепей Маркова, опи­ сывающих функционирование цифрового программного обнаружителя «к/к—/», когда на его вход поступает по­ следовательность нулей и единиц, соответствующая це­ пи С2 с переходной матрицей (2.213). Такая матрица имеет вид

 

 

Л0

Л, Аг

 

 

 

 

 

 

 

 

А

>оо

Pot

о

о

о

о

о • .

О

О

 

 

А А о 0 Ри 0

 

0 0 0 .

О О

к — 1

 

 

 

 

 

 

 

 

 

 

 

 

 

piQ

0 0 ри

0 0 о

о

о

 

р =

 

0 0 0 0 1 0 0

о

о

 

 

 

 

 

 

 

 

 

 

о

о

 

 

А1п

qv

0

0

0

0

0

0

о

 

 

 

1 0 0 0 0 0 0

о

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.220)

где

qi =

p10p'00

'—вероятность

выполнения критерия кон­

ца

пачки,

 

рі=1ці.

 

 

 

 

 

200

Для определения вектора финальных вероятностей следует использовать методику п. 2.7.1. Вероятность лож­ ной тревоги определяется как /г-я компонента финально­ го вектора:

Р = Рь = РоР,Л1.

k =

T77[,

(2.2.21)

где

 

 

 

Ро = Л о / 4 ~ ' {Р.оРоО- 1 +

РоЛҐ

(1 —

Р\\) +

+Р о Л п [ 1 + ( / - 2 ) р , 0 ^ - 1 - ( 1 - р 1 „ ^ 1 ) " - ^ - 3 ' ] } - 1 .

(2.222)

Дальнейший анализ аналогичен приведенному в п. 2.7.1. Если входная последовательность нулей и единиц — бер-

нуллиева,

т. е.

poo=Pio — q, poi — pu — p, то

переходная

матрица

(2.220)

и формулы

(2.221), (2.222)

совпадают

с известными соотношениями,

соответствующими работе

программного обнаружителя в условиях некоррелирован­ ных помех [3].

Приведем результаты оценки эффективности цифро­

вого программного

обнаружителя

«k/k—I»

в

условиях

воздействия

некоррелированных

помех

 

по

методике

п. 2.3.5. При отсутствии

флуктуации ЭОП

оптимальное

значение

k

принималось

равным £ = l , 5 l / « o , 5

(см.

табл. 2.1), 1 = 2. В табл. 2.4 представлены

потери в поро­

говой мощности по сравнению с весовым

обнаружителем

(2.57) при £> = 0,5,

F = 1 0 - 6 ,

я = 1 0 , 30, 50.

Из

таблицы

следует,

что эффективность

обнаружителя

падает

с ро­

стом числа импульсов в пачке п. При

 

Т А Б Л И Ц А 2.4

шумоподобных флуктуациях ЭОП эти

 

 

 

 

потери

еще

значительнее.

Указан­

 

п

Потерн, д Б

ное обстоятельство

ограничивает

при­

 

 

 

 

менение

программных

обнаружите­

 

10

 

1,7

лей.

 

 

 

 

 

 

 

 

2.7.3. Алгоритмы

автозахвата

тра­

 

30

 

3,0

 

50

 

5,0

екторий и сброса слежения. Вторичная

 

 

 

 

 

 

(траєкторная)

обработка радиолокаци­

 

 

 

 

онной информации связана с объединением в траектории (трассы) результатов точечных замеров координат (от­ меток) целей на фоне помеховых отметок. Все операции вторичной обработки выполняются алгоритмически в спе­ циализированных Э Ц В М . К числу задач вторичной об­ работки относятся: автоматическое обнаружение цели

201

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

На рис. 2.43 изображена упрощенная структурная схема, отображающая взаимосвязь основных операций, выполняемых при автоматическом обнаружении и COnpO-

tf/w У/Я?

Буферная

Блок

Блок

Блок

 

апробиро­

селекции,

обнаруже­

0(x,y,tk)

память

вания

отметок

ния и сброса

 

 

 

 

траектории

 

 

 

Строб

 

 

 

Начальный

Блок

Блок

 

 

empoff

 

 

 

экстрапо­

сглажива­

 

 

 

ляции

ния

 

 

 

Блок обна -

 

 

 

 

ружения

 

 

 

 

маневра

 

Рис. 2.43. Структурная схема автоматического обнаружения и сопро­ вождения одной цели.

вождении одной цели. Из устройства первичной обработ­ ки (УПО) в буферную память Э Ц В М вторичной обра-

ботки поступают отметки 0(х, у, 4 ) , обнаруженные РЛС совместно с УПО. Каждая отметка содержит: вектор ко­ ординат обнаруженной цели (илипомехи) х= (хи Xz, .. .)т, время обнаружения и /г=1, 2, ... , вектор некоординат­ ных признаков у={Уи Уг, • - • ) т (например, информацию о мощности принятого сигнала, поляризации и т. д . ) . Считанные из буферной памяти отметки поступают в бло­ ки стробирования, где производится отбор отметок, по­ падающих в заданную область (строб) пространства координат. Стробирование может быть как физическим, так и математическим. Координаты и размеры строба определяются в процессе автосопровождения.

В блоке селекции из множества отметок, попавших

202

в строб, отбирается одна (например, ближайшая к цен­ тру строба или наиболее мощная). При этом использу­ ются координаты отметки и некоордииатные признаки. В условиях большой плотности отметок селекция не про­ изводится, а траектория продолжается по каждой от­ метке, находящейся в стробе, а также по экстраполиро­ ванным отметкам. При этом число наблюдаемых траек­ торий (истинных и ложных) резко увеличивается, что предъявляет повышенные требования к объему памяти и быстродействию ЭЦВМ .

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

Координаты отобранной отметки поступают в блок экстраполяции, где вычисляются упрежденные на время Го (период обзора), координаты центра строба, а также в блок сглаживания, на выходе которого получаются сглаженные координаты сопровождаемой цели. В блоке обнаружения и сброса в соответствии с заданными ста­ тистическими критериями решаются задачи автозахвата и сброса автосопровождения. .Входной информацией это­ го блока является последовательность нулей и единиц из блока селекции. Единица соответствует наличию отметки в стробе, нуль — пустому стробу. Информация 'блока об­ наружения маневра используется для переключения алгоритмов экстраполяции и сглаживания. Следует от­ метить, что разделение задач вторичной обработки на автосопровождение и автозахват является условным и делается с целью удобства анализа.

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

203

накладывает ограничения на положение цели в прост­ ранстве. Область неопределенности в пространстве ко­

ординат,

определяемую

разрешающей

способностью

РЛС,

будем

называть

элементом

пространства.

Для

обнаружения

возможных сигналов

каждый

уча­

сток

пространства

периодически

обследуется. Период

обзора Г 0

принимается за единицу времени. -

 

 

Данные, поступающие в систему вторичной обработ­

ки, являются либо началом

новой

траектории, либо

про­

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

дой из координат на ограниченном участке

наблюдения

в виде полинома степени р:

 

 

 

 

 

 

 

р

 

х (t) =

х0+vt

+ \

f + . . . =

JJ X { t t ,

(2.223)

где x(t)—координата;

х\— коэффициенты

полинома;

Ха — начальное

положение;

xt=v—скорость;

х2/2 =

= g — ускорение и т. д. — параметры

движения.

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

Придерживаясь терминологии и обозначений [3], выражение для помехи в общем виде можно записать так:

n(t)=KO(i),P (t),n(t)]+m(t), (2.224)

где rrn(t) —случайные колебания цели вокруг заданной траектории, обусловленные влиянием случайных возму-

204

щений среды,

а также ошибками системы управления

полетом

цели;

0(>Ё)—ошибки измерения координат це­

ли

РЛС

совместно с ошибками первичной обработки,

p(t)

— ложные

отметки, обусловленные внутренними шу­

мами приемной аппаратуры, а также искусственными и естественными помехами; пЦ)—пропуски отметок от цели, обусловленные флуктуациями отражающей поверх­ ности цели и влиянием собственных шумов, искусствен­ ных и естественных помех.

Составляющая m(t) представляет собой стационар­ ный случайный процесс- с нормальным законом распре­ деления вероятностей. Математическое ожидание распре­ деления равно нулю. Корреляционная функция процесса может быть найдена экспериментально для каждого ти­ па цели и чаще всего она является экспонентной. Ошиб­ ки измерения координат имеют систематическую и слу­ чайную составляющие. Случайные ошибки измерения также подчинены нормальному закону распределения вероятностей. Для наземных Р Л С с периодом обзора в несколько секунд случайные ошибки измерения коор­ динат можно считать некоррелированными в соседних обзорах.

Если помехи вызваны внутренними шумами прием­ ника Р Л С , то ложные отметки равномерно распределе­ ны во времени и в пространстве. Их характеристикой может быть либо средняя плотность отметок в секунду в заданном объеме пространства, либо вероятность лож­ ной тревоги в элементе пространства. Степень статисти­

ческой

связи ложных

отметок в данном цикле

обзора

в большой степени зависит

от способа первичной

обра­

ботки.

 

 

 

 

 

 

 

При

автоматическом

обнаружении траекторий исполь­

зуются

критерии,

основанные как на

классической про­

цедуре

Неймана — Пирсона,

так и на

последовательной

процедуре Вальда .(главным

образом весовые

алгоритмы

с усечением). В

зависимости от критериев

автозахвата

и сброса, а также плотности шумовых отметок из блока автозахвата будут передаваться на сопровождение обна­ руженные ложные траектории. Будем считать, что веро­ ятность ложной тревоги на выходе УПО i f = c o n s t і(т. е. среднее число ложных отметок за обзор /Vj^const), среднее число отметок от целей за обзор Ыцconst и общее число отметок за обзор меньше числа элементар­ ных участков в зоне обзора (-/Ул+М») <Ng где Д/8 для

205

двухкоординатнои Р Л С определяется в соответствии

свыражением

у'max

' о

целевые и шумовые траектории ие пересекаются, за вре­ мя наблюдения размеры стробов постоянны.

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

Из блока

 

 

Необнаружение

селекции

Критерий,

Накопле­

Решение Обнаружение

 

 

начала

ние

 

 

 

 

Подтверждение

 

 

Нано/гме-

Решение

 

 

ние

 

 

Неподтверждение

 

 

 

Рис. 2.44. Структурная схема алгоритмов автозахвата и подтвержде­ ния траектории.

алгоритм автозахвата включает в себя: а) критерии фиксации начала траектории, б) алгоритм накопления за ограниченное время наблюдения после фиксации на­ чала, в) критерий вынесения решения (обнаружения или необнаружения). В свою очередь, алгоритм подтвержде­ ния сопровождения обнаруженной траектории состоит из алгоритма накопления одиночных обнаружений вдвижущемся окне и критерия вынесения решения (подтверж­ дения или неподтверждения траектории). Сказанное ил­ люстрируется рис. '2.44, на котором изображена струк­ турная схема алгоритмов автозахвата и подтверждения (сброса) траектории. Фиксация начала траектории мо­ жет производиться: а) при появлении отметки, не при­ надлежащей ни одной траектории, б) при появлении от­ метки значительной амплитуды с вероятностью Q, в) при появлении двух отметок в двух смежных ииклах обзора.

206

В последнем

случае

облегчается задача

экстраполяции

на следующие циклы

обзора.

 

 

 

Накопление после фиксации начала траектории осу­

ществляется

за ограниченное время наблюдения

(3 —

10 циклов обзора). При наличии статистической

связи

между отметками от цели

или помехи

целесообразно

использовать

оптимальный

весовой алгоритм накопле­

ния, рассматриваемый ниже. Находят применение также следующие критерии: весовой вариант усеченной после­ довательной процедуры (см. п. 2 . 6 . 2), цифровые нако­ пители позиций (см. п. 2.7.1), цифровые программные обнаружители (см. п. 2 . 7 . 2) . В соответствии с критерием накопления функционирует блок решения, который срав­ нивает накопленную сумму с одним либо с двумя порога­ ми. Алгоритм подтверждения сопровождения наиболее просто реализуется в виде обнаружителя движущегося окна. Траектория подтверждается при наличии /г единиц на п позициях. Возможно также использование крите­ рия — «/ пропусков подряд».

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

Вероятность появления отметок в стробе

автозахвата,

состоящем из М независимых элементов,

определяется

из (2. 22):

 

 

при отсутствии цели

 

 

P^FM^MF,

при М Л < 1,

(2. 225)

при наличии цели

 

 

/>«'£>+(1—D)FM

при M F ! « 1 .

(2.226)

Переходные «матрицы накопителей, используемых при автозахвате, рассматривались выше. Так, например, для цифровых накопителей позиций и' цифрового програм­ много обнаружителя переходные матрицы определяются из (2.215) и (2.220) соответственно. Начальные векто­ ры в зависимости от критерия фиксации начала траек­ тории для переходных матриц (2.215) и (2.220) будут следующими:

207

а) при появлении отметки, не принадлежащей ни од­ ной из траекторий

 

 

?о =

(А . Рі<

Р»

А . ) т =

(1 -

Р.

Р. О, 0,..

,

0)',

 

 

где р

определяется из

('2. 225),

(2.

226);

 

 

 

 

 

 

 

б)

при появлении

отметки-значительной

 

амплитуды

 

 

 

 

 

J e

=

( l — Q . Q . 0 ,

0

 

0)';

 

 

 

 

 

 

в)

при появлении

двух отметок

подряд

 

 

 

 

 

 

 

 

 

 

Я

= ( 1 - Р 2 . 0 , / Л 0

 

 

0)т.

 

 

 

 

 

АІ

Таким

образом,

траектория

возникает

состоянии

или Л2 ,

сбрасывается в процессе ее обнаружения

при

достижении

состояния

Л о, обнаружение

происходит

при

 

 

 

 

 

 

tOOff ООО//

достижении

состояния,

соот­

 

 

 

 

 

 

ветствующего критерию

об­

 

АВтоэах.

 

/ /

наружения.

 

 

 

 

 

пра­

 

 

 

 

 

 

 

Расчет

вероятностей

 

1

*

 

 

^0/1

. 00///

 

вильного

и ложного

обнару.

 

 

р\

 

Р2?г/

 

жения

(автозахвата) для ци­

 

 

 

t

/

J

 

 

фровых накопителя

позиций

 

/ff/

 

Off

/

Сорос

и программного обнаружив­

 

 

ая рассматривался в ли. 2.7.1

 

 

Р5

 

* р

у

 

 

 

І

 

f

/

 

 

 

 

 

и

2.7.2

соответственно,

 

и мы

 

 

 

 

 

 

не

 

будет

 

останавливаться

 

 

 

(о/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

на этом вопросе. В ряде слу­

 

 

Р1

 

 

 

 

 

 

I

 

 

 

 

 

 

 

чаев

 

получение

 

указанных

 

 

 

°ч

 

 

 

 

характеристик можно

упро­

 

 

 

 

 

 

 

стить,

использовав

методи­

 

0

 

/

2

 

3

 

 

ку, связанную -с плоскостью

 

 

 

Число

нулей."

 

случайных

блужданий

 

(ом.

Рис. 2.45. Схема расчета харак­

п. 2.6.2).

Рассмотрим

в

ка­

честве

примера

расчет

ха­

теристик

автозахвата

для

кри­

рактеристик

автозахвата

 

 

терия

«3/5—3».

 

 

 

 

 

 

цифрового

программного об­

 

 

 

 

 

 

 

 

 

наружителя, фиксирующего

начало траектории

по

 

двум

единицам подряд, автозахват траектории — по крите­ рию «3/5—3», т. е. обнаружение происходит при наличии трех единиц в пяти позициях и сброс автозахвата — при

наличии трех

нулей

подряд. Схема расчета приведена

на рис. 2. 45.

На первом шаге возможно появление еди­

ницы с вероятностью

р и нуля с вероятностью q=\—р.

208

;В последнем случае сразу осуществляется пересечение порога «сброс». Заполнение следующей строки схемы производится- в соответствии с результатами предыду­ щей. Стрелками указывают возможные направления пе­ реходов из одного состояния в другое. Из поглощающих состояний «автозахват» и «сброс» переходов быть не может.

Нетрудно убедиться, что для выбранного критерия вероятность автозахвата за три шага равна p^J = р*, за четыре—р*ю = р1*} +p*q и т. д.

Для расчета среднего числа шагов до автозахвата используем формулу математического ожидания диск­ ретной случайной величины (2.202). При определении характеристик правильного автозахвата ра з=-Ьаз сле­ дует подставить значения р из (2.226), а при определе­ нии характеристик ложного' /?аз=^лз из (2.225). При наличии корреляции между обзорами методика сохра­ няется, изменяются лишь вероятности появления той или иной комбинации из единиц и нулей.

Перейдем теперь к вопросу построения оптимально­ го весового алгоритма автозахвата, когда последователь­ ность обнаружений и необнаружений образует однород­ ную односвязную цепь Маркова [63, 101]. Начало траек­

тории будем

фиксировать

по двум

единицам

подряд.

В соответствии

с алгоритмом

 

(2. 204)

имеем

 

 

 

 

 

' l n A = U 7 0 + S

 

£

 

(k) Wa?

,

 

 

(2.227)

 

 

 

 

 

 

ft=l

oB

 

 

 

 

 

 

где W0 = I n ( p n

s N I P U N ) —

начальный

вес;

 

^ ( ^ — и н д и ­

катор

перехода

из

состояния

 

Аа

в

состояние

А^ [см.

(1.11)]; ар=Ы(ра^/ра^),

 

 

а,

р = 0 Т Т ,

pa?SN

 

и p a ? N -

переходные

вероятности

цепей

Маркова,

соответствую­

щих последовательностям

единиц и

нулей

для

 

сигнала

с помехой и помехи соответственно. Переходные

матри­

цы таких цепей имеют вид

(2.

213).

 

 

 

 

 

Структурная

схема .алгоритма

|(2. 227),

выполненная

для пояснения

ее

функционирования

на

элементах вы­

числительной

техники,

приведена

на

рис.

2. 46.

Входная

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

единиц

и нулей

из блока

селекции

отметок

подается

на

'вход

 

двухразрядного

регистра

14—1410

209

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