
книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации
.pdfхи, компоненты которого определяются из (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 |