Телекомунікації, радіолокація, радіонавігація та електроакустика
УДК 621.396.96
АЛГОРИТМ АДАПТИВНОГО МНОГОАЛЬТЕРНАТИВНОГО ПОСЛЕДОВАТЕЛЬНОГО ОБНАРУЖЕНИЯ ТРАЕКТОРИИ ЦЕЛИ С ИСПОЛЬЗОВАНИЕМ РЕШАЮЩИХ СТАТИСТИК ОТМЕТОК И
ОТБРАСЫВАНИЕМ НЕУДАЧНЫХ ГИПОТЕЗ1
Неуймин А. С. ; Мешков С. И. ; Жук С. Я. , д.т.н., профессор
1 Национальный технический университет Украины
«Киевский политехнический институт», г. Киев, Украина, dragon_sasha@meta.ua
2Военно-дипломатическая академия имени Евгения Березняка,
г. Киев, Украина
AN ADAPTIVE MULTIALTERNATIVE SEQUENTIAL TARGET TRACK DETECTION ALGORITHM USING THE DECISION STATISTICS OF PLOTS AND DISCARDING THE FAILED HYPOTHESES
1 2 1
Neuimin O. S. ; Mieshkov S. I. ; Zhuk S. Ya. , Doctor of Engineering, Professor
1 National Technical University of Ukraine “Kiev Polytechnic Institute ”, Kyiv, Ukraine,
2
Eugene BereznyakMilitary-Diplomatic Academy, Kyiv, Ukraine
Введение
Широкое распространение для оптимального обнаружения траектории цели находит последовательный критерий Вальда [1, 2]. Методы на его основе позволяют получить существенный выигрыш во времени при заданных требованиях к вероятностям принимаемых решений по сравнению с оптимальным методом принятия решения при фиксированном числе обзоров по критерию Неймана-Пирсона. В [2] на основе критерия Вальда получены алгоритмы последовательного обнаружения траектории цели c использованием решающих статистик (РС) отметок в стробе сопровождения. В данных алгоритмах полагается, что отношение сигнал/шум (ОСШ) известно. Такое предположение часто не выполняется на практике, поскольку класс обнаруживаемой цели часто является неизвестным. При этом эффективная площадь рассеивания (ЭПР) различных классов целей может различаться на порядок. Поэтому важное значение имеет разработка адаптивного алгоритма обнаружения траектории цели, имеющего возможность наряду с решением главной задачи обнаружения траектории цели, выполнить и оценивание ОСШ, что позволит также распознавать класс воздушной цели по величине ЭПР.
Особенностью рассматриваемой задачи является сложный нелинейный характер зависимости между измерением и неизвестным параметром, а
1 Електронний варіант статті: http://radap.kpi.ua/radiotechnique/article/view/1053
Вісник Національного технічного університету України «КПІ» 82 Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №61
Телекомунікації, радіолокація, радіонавігація та електроакустика
также незначительное число наблюдений. Поэтому для синтеза адаптивных алгоритмов целесообразно применить подход, при котором область возможных значений ОСШ дискретизируется, что эквивалентно выдвижению гипотез относительно его значений [3].
В [4, 5] на основе последовательного критерия простого дополнения получено решение задачи многоальтернативной проверки статистических гипотез с использованием верхних и нижних порогов, которые определяются на основе заданных условных вероятностей распознавания и априорных вероятностей гипотез. Использование нижних порогов позволяет реализовать процедуру отбрасывания неудачных гипотез.
В работе на основе последовательного критерия простого дополнения разработан алгоритм адаптивного многоальтернативного последовательного обнаружения траектории цели с использованием РС отметок и отбрасыванием неудачных гипотез. В данном алгоритме оценивается неизвестное ОСШ, что позволяет распознавать класс воздушной цели по величине ЭПР.
Постановка задачи
В оптимальном приемнике первичной обработки сигналов на каждом обзоре выполняется обнаружение цели путем сравнения РС в элементах разрешения с входным порогом Нвх, выбор которого проводится в соответствии с требуемой вероятностью ложной тревоги FBX, а также измерение ее координат. В случае обнаружения формируется отметка, а также сохраняется значение полученной РС.
На каждом обзоре по всем отметкам, не попавшим в стробы сопровождения, завязываются траектории. В дальнейшем для выделения траектории цели используется алгоритм сопровождения, с помощью которых выполняется экстраполяция координат цели, стробирование и идентификация отметок [2]. В общем случае в стробе сопровождения на к -м обзоре содержится Мк отметок с РС zk, m = 1, Mk, которые объединяются в вектор Zk.
В результате дискретизации области возможных значений ОСШ будем полагать, что ОСШ принимает L значений ql, l = 1, L. При этом имеют место L+1 простые гипотезы: Н0 — об отсутствии траектории цели и Hl, l = 1, L — о наличии траектории цели с ОСШ ql, l = 1, L.
Функции правдоподобия (ФП) f (Zk,Мк | Н) простых гипотез Н, l = 1, L при наличии цели определяются, как и в двухальтернативной задаче [2]. При этом необходимо рассмотреть два случая: Мкф0 — в стробе находятся отметки, Мк = 0 — в стробе отметок нет. ФП
Вісник Національного технічного університету України «КПІ»
Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №61 83
Телекомунікації, радіолокація, радіонавігація та електроакустика
f (Zk, Mk | Ht) при Mk Ф 0 определяется выражением
M,
f(Zk,Mk | H,) = Pm, (1 -Овх^стр)П fN(zm) +
m=1
Mk D Р Mk . __
+TPMk-idbXthlfiw) П pn(4), і=i,=,
k-1 M
m=1 iMk j=1, j Ф m
где DBX — вероятность превышения РС целевой отметки входного порога HBX; Рстр — вероятность попадания отметки цели в строб; Рщ — вероятность появления Mk ложных отметок на k -м шаге в стробе размером FCTp, которая определяется по закону Пуассона; fS(zm),
fу(zm) — усеченные законы распределения РС, превысивших входной
н
порог вх, при условии, что отметка является целевой и ложной соответственно.
Для гипотезы H0 ФП f (Zk,Mk | H0) определяется по формуле
Mk
f(Zk,Mk і hо)=pMt пmzm).
m =1
Случай отсутствия в стробе сопровождения отметок является вырожденным и известно лишь, что Mk = 0. При этом ФП f (Zk, Mk 1 Hl ) 5 ^ = 0 L имеет вид
f (Zk, MklHt) = P(Mk = 0|Hi)
(1 - АхРстр )Р0, і = 1, L,
Р
0’
і = 0.
где Р0 — вероятность отсутствия в стробе ложных отметок.
В соответствии с постановкой задачи многоальтернативной проверки гипотез [4, 5] необходимо при условных вероятностях ошибочного распознавания гипотез и вероятностях правильного распознавания, не ниже заданных, получить последовательное решающее правило, которое позволит по наблюдениям Z k, Mk, k = 1,2,3... принять одно из решений H, i = 0, L.
Разработка алгоритма адаптивного многоальтернативного последовательного обнаружения траектории цели
В соответствии с последовательным критерием простого дополнения, задача многоальтернативного распознавания L+1 простых гипотез сводится к L+1 двухальтернативным задачам проверки простой H гипотезы против сложной альтернативы 0j. Сложная гипотеза 0j является объединени-
84
Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №61
Телекомунікації, радіолокація, радіонавігація та електроакустика
ем простых гипотез Н t, / = 0,Л, , 0/ = (J Ні . Геометрическая ин-
/=0,/W
терпретация сложной гипотезы представлена на рис. 1.
0 Для решения L+1 двухальтерна-
1 тивных задач распознавания на каж-
дом k-м обзоре определяются L+1 отношений правдоподобия (ОП)
Л, (k)
f (Zk, Mk\H,)
Рис. 1.
f (Zk, Mk|0,)
, = 0, L,
где f(Zk, Mk | Hl) — ФП про-
стой гипотезы Hj, которая определяется с помощью выражения f (Zk, Mk\Hj) = f (Zk, Mk\H,), f (Zk - \ Mk -1 |Hj).
f (Zk,Mk | 0,) — ФП сложной гипотезы 0,; Zk — набор измерений, полученных вплоть до k-го обзора.
В соответствии с [4] ФП f (Zk,Mk | 0,) можно представить в виде
Z f (Zk,Mk | Hi)Pi
f (Zk,Mk 10,) =
L
Z Pi
i=1, i
Z Уі vf (Zk, Zk I Z),
i =1, i
где уi| ,, i = 1, L, i ^, — условная вероятность гипотезы Hj при условии,
что имеет место сложная альтернатива 0,, которая определяется по формуле
L
Уі | , = Pi1 Z Pi.
i=1, і
Для принятия решения ОП Л, (k) сравниваются с нижними порогами
Bl, I = 0, L [4]
Л,) =
f (Zk, Mk!Ht) f (Zk, Mk|0,)
< B,
, = 0, L.
(1)
При выполнении условия (1) при , = j принимается решение в пользу сложной гипотезы 0 ., которое эквивалентно решению об отклонении простой гипотезы H■. Таким образом, решение в пользу простой гипотезы H может быть принято после того, как отклоняются остальные L простых
Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №61
85
Телекомунікації, радіолокація, радіонавігація та електроакустика
гипотез H, / = 0, L, / ф і, для которых одновременно выполняются условия (1). Геометрическая интерпретация процедуры принятия решения с отклонением неудачных гипотез представлена на рис. 2.
Рис. 2.
В [4, 5] показано, что нижние пороги В/, / = 0,L определяются по формуле:
в 1 - P(H/1Н/)
1 1 -Р(Н/ і©/)
(2)
где P(Н/ і ©/) — вероятность ошибочного распознавания гипотезы Н/ при условии, что имеет место сложная альтернатива ©/, которая вычисляется по формуле
L
Р(Н,\®,) = Z Yi|/P(Н,\НІ).
і=1, і ф /
Возможность отклонения неудачных простых гипотез позволяет реализовать процедуру их отбрасывания.
Алгоритм многоальтернативного последовательного обнаружения траектории цели по нижним порогам с использованием решающих статистик отметок и отбрасыванием неудачных гипотез можно представить в следующем виде. После получения очередных данных на
k -м обзоре, происходит вычисление и сравнение Л/ (k) с нижними порогами В/, / = 0,L.
Л1 (k)
f( Zk, Mk\H/)
L------------------B
Z y, V f (Zk, Mk\Hi)
і=1, і ф/
(3)
Если все Л/ (k) > В /, 1 = 0, L, то проводится следующий обзор. В случае если выполняется одно или несколько условий Л/ (k) — В/, производится отбрасывание соответствующих простых гипотез Н . Если осталась одна
86
Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №61
Телекомунікації, радіолокація, радіонавігація та електроакустика
простая гипотеза Нь то выполняется прекращение наблюдений и принимается решение в ее пользу. В качестве оценки ОСШ q принимается значение соответствующее принятой гипотезе. На основе полученной оценки ОСШ, с использованием уравнения радиолокации, может быть определена ЭПР цели
а =
(4n)3R4kTL-2 -—-—
PuxG\2
что позволяет распознавать ее класс.
Если сохранилось 1 < N < L простых гипотез Hj, проводится коррекция
постановки задачи, которая включает следующие этапы:
1) Перенормировка априорных вероятностей оставшихся N гипотез
N
p'j=pj / s Pi ■
i=1
2) Уточнение заданных условных вероятностей распознавания P3(H | H), i, l = 0,N для оставшихся гипотез. При этом условные вероятности правильного распознавания Рз (Hl | Ht), l = 0, N задаются такие же, как и в исходной задаче, а условные вероятности ошибочного распознавания Рз (H | H), i, l = 0, N, і ф l определяются с использованием свойства нормировки условных вероятностей.
3) Вычисление на основе полученных p \, i = 0, N и Рз (H | H), i, l = 0, N значений новых порогов Bl, l = 0, N по формуле (2) для следующего шага.
Применение полученного алгоритма позволяет сократить вычислительные затраты, поскольку необходимо рассчитывать ФП и выполнять проверку условия в соответствии с (3) только для оставшихся гипотез.
Анализ эффективности разработанного алгоритма
Анализ эффективности разработанного алгоритма выполним на примере обнаружения траектории цели по данным обзорной РЛС, измеряющей дальность ги и радиальную скорость ги цели. При отсутствии и при наличии цели нормированные РС подчиняются усеченным законам распре-f у (zm) f у (zm)
деления JNK k J и JS v k J [1, 2], соответственно
fN (zm)
1 1
F 2
А вх
exp
—m
zk
, fS (zm)
1
1
2(1 + q2/2)
exp
—m
zk
2(1 + q1 / 2)
Для выделения траектории цели применяется простейший алгоритм сопровождения [2], базирующийся на описании движения цели с помощью модели второго порядка. При отсутствии в стробе отметок для продолже-
Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №61
87
Телекомунікації, радіолокація, радіонавігація та електроакустика
ния сопровождения используются экстраполированные характеристики параметров движения цели на текущем обзоре. Для продолжения траектории выбирается отметка с максимальной РС в стробе сопровождения [2].
Темп поступления данных полагался T = 1с. Координаты ложных отметок в стробе имеют равномерный закон распределения. Вероятность ложной тревоги в элементе разрешения полагалась ^вх = 0.3, среднее число ложных отметок Мл = 2. Вероятность попадания в стробы отметок цели
Рстр = 0.994. СКО измерения РЛС устанавливались равными ar = 300 м, Gf. = 10 м/с. СКО шума возбуждения, характеризующего случайное ускорение цели, полагалось a a = 2м/с2.
Полагается, что ОСШ может принимать значения q1 = 2.5, q2 = 7, q3 = 16. Таким образом, L =3 и имеет место четыре простые гипотезы H, l = 0,3. Моделирование проводилось по 104 испытаниям.
Условные вероятности распознавания Рз (H11H), i, l = 0,3 заданы в виде матрицы
P,
0.99 0.05
0.0005 0.95 0.003 0 0.0065 0
0.05 0.05 00 0.95 0 0 0.95
(4)
10-
3-1.00
п
о
а
о
со
ю
о
о
ч
о
Я
Я
о
о
я
я:
о
а
U
8-
6-
4-
2-
1 - Цели нет
2 - Цель с ОСШ q=2.5
3 - Цель с ОСШ q=7
4 - Цель с ОСШ q = 16
2 3
Ситуация
Рис. 3
же заданных (4).
9
7
5
3
1
1
4
На рис.3 представлены среднее число обзоров и вероятности правильного обнаружения ложной траектории и траектории цели с ОСШ qr = 2.5, qr = 7, qr = 16, полученные разработанным алгоритмом адаптивного многоальтернативного обнаружения траектории цели. Как следует из рис.3 с увеличением ОСШ цели время обнаружения траектории цели уменьшается до трех периодов обзора. При этом алгоритм обеспечивает показатели распознавания гипотез не ни-
Выводы
Разработанный алгоритм адаптивного многоальтернативного последо-
88
Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №61
Телекомунікації, радіолокація, радіонавігація та електроакустика
вательного обнаружения траектории цели с использованием РС отметок и отбрасыванием неудачных гипотез, наряду с решением главной задачи обнаружения траектории цели, выполняет также оценивание ОСШ, что позволяет распознавать класс воздушной цели по ее ЭПР.
Отбрасывание неудачных гипотез позволяет сократить вычислительные затраты, поскольку необходимо рассчитывать ФП и выполнять проверку условия (3) только для оставшихся гипотез.
Для рассмотренного примера с увеличением ОСШ цели время обнаружения траектории цели уменьшается до трех периодов обзора. При этом алгоритм обеспечивает показатели распознавания гипотез не ниже заданных.
Перечень источников
1. Keuk V. Sequential track extraction / Van Keuk G. // IEEE Transactions on Aerospace and Electronic Systems. — 1998. — V.34. — № 4. — P. 1135-1148. doi:10.1109/7.722699
2. Неуймин А.С. Последовательное обнаружение траектории цели с использованием решающих статистик отметок / А.С. Неуймин, С.Я. Жук // Изв. вузов. Радиоэлектроника.— 2014.— № 6.— С. 35-46.
http://radio.kpi.ua/article/view/S0021347014060041
3. Стратонович Р.Л. Принципы адаптивного приема / Р.Л. Стратонович. — М. : Сов. радио, 1973. — 143 с.
4. Жук С. Я. Многоальтернативное последовательное решающее правило с отбрасыванием неудачных гипотез / С.Я. Жук, В.И. Ковалев // Проблемы управления и информатики.— 2000.— № 4.— С. 88-96.
5. Жук С. Я. Совместная фильтрация состояния и распознавание типа структуры динамической системы с отбрасыванием неудачных гипотез / С.Я. Жук, В.И. Ковалев // Известия вузов. Радиоэлектроника.— 2001.— № 7.— С.16-26.
References
1. Keuk V. Sequential track extraction / Van Keuk G. // IEEE Transactions on Aerospace and Electronic Systems. — 1998. — V.34. — № 4. — P. 1135-1148. doi:10.1109/7.722699
2. Neuimin O. S., Zhuk S. Ya. Sequential detection of target trajectory using the decision statistics of pips. Radioelectronics and Communications Systems, 2014, vol. 57, no. 8, pp. 352-361. doi: 10.3103/S0735272714080032
3. Stratonovich R.L. Principy adaptivnogo priema [Adaptive receiving principles]. Moscow, Sov. radio Publ., 1973, 143 p.
4. Zhuk S.Ya., Kovalev V.I. The multiple-choice sequential decision rule with rejection of unfortunate hypotheses. Journal of Automation and Information Sciences, 2000, vol. 32, no. 10, pp. 52-58. DOI: 10.1615/JAutomatInfScien.v32.i10.70
5. Zhuk S.Ya., Kovalev V.I. State joint filtering and recognition of the structure type of a dynamic system with rejection of unsuccessful hypotheses Izvestiya Vysshikh Uchebnykh Zavedenij. Radioelektronika, 2001, vol. 44, no. 7, pp. 16-26 (in Russian).
Неуймін О. С., Мєшков С. І., Жук С. Я. Алгоритм адаптивного багатоальтерна-тивного послідовного виявлення траєкторії цілі з використанням вирішальних статистик відміток і відкиданням невдалих гіпотез. На основі послідовного критерію простого доповнення розроблений алгоритм адаптивного багатоальтернативного послідовного виявлення траєкторії цілі з використанням вирішальних статистик від-
Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №61
89
Телекомунікації, радіолокація, радіонавігація та електроакустика
міток і відкиданням невдалих гіпотез. На основі отриманої оцінки ВСШ, може бути визначена ЕПР цілі, що дозволяє розпізнавати її клас. Аналіз отриманого алгоритму виконаний за допомогою статистичного моделювання на прикладі виявлення траєкторії цілі за даними оглядової РЛС, що вимірює дальність і радіальну швидкість цілі.
Ключові слова: виявлення траєкторії цілі, вирішальні статистики, невідоме ВСШ.
Неуймин А. С., Мешков С. И, Жук С. Я. Алгоритм адаптивного многоальтернативного последовательного обнаружения траектории цели с использованием решающих статистик отметок и отбрасыванием неудачных гипотез. На основе последовательного критерия простого дополнения разработан алгоритм адаптивного многоальтернативного последовательного обнаружения траектории цели с использованием решающих статистик отметок и отбрасыванием неудачных гипотез. На основе полученной оценки ОСШ, может быть определено ЭПР цели, что позволяет распознавать ее класс. Анализ полученного алгоритма выполнен с помощью статистического моделирования на примере обнаружения траектории цели по данным обзорной РЛС, измеряющей дальность и радиальную скорость цели.
Ключевые слова: обнаружение траектории цели, решающие статистики, неизвестное ОСШ.
Neuimin O. S., Mieshkov S. I, Zhuk S. Ya. An adaptive multialternative sequential target track detection algorithm using the decision statistics of plots and discarding the failed hypotheses.
Introduction. It is shown that the synthesis of an adaptive target track detection algorithm having the opportunity along with the decision of the main task to perform SNR estimation is of great practical importance.
Problem statement. In accordance with the problem statement of multialternative hypotheses testing it is necessary to obtain a sequential decision rule at conditional probability of erroneous recognition hypotheses and correct recognition probabilities, not less than specified that will allow for observation to take one of the hypotheses.
Development of an multialternative sequential target track detection algorithm. On the basis of the simple additions sequential test an adaptive multialternative sequential target track detection algorithm using the decision statistics of plots and discarding the failed hypotheses is developed.
Effectiveness Analysis of Algorithm. Analysis of the adaptive algorithm is carried out as an example of target trajectory detection due to surveillance radar which measures range and range rate using the statistical modeling. The decision statistics related to noise- and target-originate measurements are described by an exponential distribution and Swerling 1 distribution respectively.
Conclusions. On the basis of the SNR estimates can be determined radar cross-section of the target that can recognize its class. Discarding the failed hypotheses can reduce the computational cost. Designed adaptive algorithm provides hypotheses recognition performance is not less than specified.
Keywords: target track detection, decision statistics, unknown SNR.
90
Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №61