Математические методы идентификации пеленгов беспилотных летательных аппаратов в группе

Размещено на http://www.allbest.ru/

Казанский национальный исследовательский технический университет

им. А.Н. Туполева-КАИ, Казань

Математические методы идентификации пеленгов беспилотных летательных аппаратов в группе

А.В. Лукин2

И.С. Гайнутдинов2

Е.П. Денисенко3

Революционное применение беспилотных летательных аппаратов (БПЛА) в различных отраслях промышленности, услуг, связи, решении антитеррористических задач и др., выдвигают на первое место задачи их траекторных измерений и идентификации пеленгов, особенно при движении в «рое» группе [1-3], при меняющихся погодных условиях [4-8] и использовании различных оптических и оптико-электронных средств [9-13].

Задача по идентификации пеленгов БПЛА в группе может быть сформулирована следующим образом.

В некоторой области пространства размещены произвольным образом n однотипных объектов. Имеются три измерительных пункта (ИП) с координатами средств измерения (СИ) X1,Y1,Z1, X2,Y2,Z2 и X3,Y3,Z3. Из этих пунктов в один и тот же момент времени измерены пеленги объектов: (ij, ij), где i = 1, 2, 3, j = 1, 2, …, n, а затем пронумерованы в произвольном порядке независимо друг от друга. Требуется идентифицировать пеленги, полученные разными СИ.

Поскольку в задаче предполагается наличие угловых координат объектов, то идентификация пеленгов может производиться с использованием либо только горизонтальных (азимутальных) углов, полученных с разных ИП, либо не только горизонтальных, но и вертикальных углов (углов места). Первый вариант назовем идентификацией горизонтальных пеленгов (ИГП), второй - идентификацией пространственных пеленгов (ИПП).

Идентификация горизонтальных пеленгов

В ситуации, когда группа БПЛА наблюдается и регистрируется под углами места, близкими или равными нулю, при идентификации пеленгов могут быть использованы только азимутальные углы.

Пусть в какой-то момент времени пеленгатор первого измерительного пункта (ИП-1) зарегистрировал l объектов, азимуты пеленгов которых равны 11, 12, …, 1l. Пеленгаторы ИП-2 и ИП-3 зарегистрировали в это же время соответственно m и n объектов (l m n), азимуты которых равны 21, 22, …, 2m и 31, 32, …, 3n (угловые координаты приведены к единой нормальной земной системе).

Уравнение прямой, проходящей в горизонтальной плоскости через точку стояния пеленгатора ИП-1 в направлении 1i (i номер объекта), имеет вид

x=x1+tg(/21i)(zz1),

где x1, z1 координаты точки стояния пеленгатора ИП-1, м.

Аналогично записывается уравнение прямой с использованием j-го объекта, зарегистрированного пеленгатором ИП-2:

x=x2+tg(/22j)(zz2),

где x2, z2 координаты точки стояния пеленгатора ИП-2, м.

Решая совместно уравнения (1) и (2), получим координаты точки пересечения проекций на плоскость XOZ направлений на i-й объект с ИП-1 и j-й объект с ИП-2:

xij=;

zij=.

Находим азимутальное направление на точку (xij, zij) с третьего ИП:

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

3=

где x3, z3 координаты точки стояния пеленгатора ИП-3, м.

Затем, осуществляя перебор измерений пеленгатора ИП-3 3k (k=1, 2, …, n), и проверяя выполнение неравенства принадлежности азимутов 1i, 2j, 3k (соответственно с ИП-1, ИП-2 и ИП-3) одному и тому же объекту измерения. В случае выполнения неравенства (5) считаем, что триада пеленгов (1i, 1i)_(2j, 2j)_(3k, 3k) прошла идентификацию по азимутальным углам.

3k3<,

При нормальном законе распределения в качестве критерия, согласно известному правилу «трех сигм», можно принять выражение:

,

где радикал - СКО разности двух независимых азимутов, - СКО k-го азимута, измеренного пеленгатором ИП-3, а второе слагаемое - дисперсия расчетного азимута , значение которой может быть получено численными методами из выражения:



Однако на практике нашел применение и показал хорошие результаты более гибкий - адаптивный критерий, величина которого настраивается в зависимости от конкретных значений погрешностей всей совокупности пеленгов.

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

Первоначально идентификация проводится с некоторой априорно заданной величиной критерия 0. Если результат оказался успешным, то значение критерия для следующего шага уменьшается вдвое: i+1 = i /2.

В противном случае значение критерия идентификации увеличивается. Увеличение критерия происходит по-разному в зависимости от того, была хотя бы один раз на предыдущих шагах успешная идентификация или нет. Если успешная идентификация была, то значение критерия рассчитывается по формуле i+1=i+(yi)/2, где y значение критерия последней успешной идентификации пеленгов.

Если же успешной идентификации не было, то значение критерия увеличивается вдвое i+1 = 2i.

Процесс подбора критерия идентификации проводится до тех пор, пока не выполнится условие: i+1i , в котором - задаваемое минимальное изменение критерия идентификации.

Идентификация пространственных пеленгов

Алгоритм идентификации пространственных пеленгов основан на следующей аксиоматической посылке: линии визирования (пеленги) одного и того же объекта с двух точек с точностью до погрешностей визирования принадлежат одной плоскости.

Пусть направления на объект с двух ИП с координатами X1,Y1,Z1 и X2,Y2,Z2 определяются горизонтальными и вертикальными углами 1,1 и 2,2 (координаты и углы даны в единой нормальной земной системе координат). Уравнение плоскости P12, содержащей точку (X2,Y2,Z2) и направление на объект из точки (X1,Y1,Z1) имеет вид

где l1 =cos1cos1, m1=sin1, n1 =cos1sin1.

Уравнение прямой, проходящей через точку (X2,Y2,Z2) и объект,

.

где l2 =cos2cos2, m2=sin2, n2 =cos2sin2.

Угол между прямой (7) и плоскостью (6), имеющий место из-за случайных погрешностей измерений, вычисляется по формуле

,

в которой коэффициенты общего уравнения плоскости P12 (A1X+B1Y+C1Z+D1=0) в соответствии с (7) определяются выражениями

A1=(y2y1)n1(z2z1)m1,

B1=(z2z1)l1(x2x1)n1,

C1=(x2x1)m1(y2y1)l1.

При известных оценках СКО измерений (1, 1, 2, 2) можно оценить дисперсию угла

,

и с априорно заданной вероятностью утверждать, что

K(q)2.

Для нормального закона распределения в соответствии с «правилом трех сигм» можно принять K(q)=3, при этом q=0,9972. Если неравенство (6.11) не будет выполняться, то следует предположить либо грубую ошибку, хотя бы в одной из величин 1,1, 2,2, либо принадлежность угловых координат (попарно) к различным объектам.

Таким образом, неравенство (11) может служить критерием идентификации. Оно используется следующим образом. Выберем любую пару угловых координат, зарегистрированную пеленгатором ИП-1 и определяющую единичный радиус-вектор (пеленг) одного из объектов. Через этот вектор и точку стояния пеленгатора ИП-2 проведем плоскость P12. Вычислим по формуле (8) углы между плоскостью P12 и пеленгами (единичными радиус-векторами) объектов, зарегистрированными с ИП-2. При проверке этих углов по критерию (11) возможен один из следующих двух исходов.

Неравенство (11) не выполняется ни для одного из углов. Вывод: координаты выбранного объекта пеленгатором ИП-2 не зарегистрированы; необходимо перейти к идентификации измерений другого объекта.

Неравенство (11) выполняется хотя бы для одного из углов. Вывод: выполнено необходимое условие идентификации.

В плоскости Р12 может оказаться и несколько объектов. Проведем плоскость Р13 (через выбранный ранее единичный радиус-вектор в точке стояния пеленгатора ИП-1 и точку стояния пеленгатора ИП-3). Вычислим углы между плоскостью Р13 и единичными радиус-векторами объектов с началом в точке стояния пеленгатора ИП-3. В результате применения критерия (11) возможны исходы, аналогичные уже рассмотренным. Разумеется перспективным является второй исход выполнение еще одного необходимого условия.

Итак, имеется хотя бы один пеленг объекта, зарегистрированный пеленгатором ИП-2, в плоскости Р12 и хотя бы один пеленгатором ИП-3 в плоскости Р13. пространственный пеленг засечка статистический

В частном случае, когда в плоскостях Р12 и Р13 имеется только по одному пеленгу, проведем плоскость Р23 (плоскость Р32), проходящую через пеленг с ИП-2 (ИП-3) и точку стояния пеленгатора ИП-3 (ИП-2). Вычислим угол между плоскостью Р23 (Р32) и пеленгом с ИП-3 (ИП-2). Если неравенство (11) не выполняется, то переходим к идентификации измерений другого объекта. При выполнении неравенства (11) заключаем, что измеренные с трех ИП пеленги относятся к одному и тому же объекту.

Исключением из этого является случайное событие (рис.1), порождающее ложную засечку. В общем случае может быть найдено по нескольку пеленгов в плоскости Р12 и плоскости Р13. Назовем эти пеленги выделенными.

Размещено на http://www.allbest.ru/

Для утверждения того факта, что два пеленга из числа выделенных от ИП-2 и ИП-3 определяют направления на один и тот же объект, необходимо, чтобы они принадлежали (с точностью до погрешностей измерений) одной плоскости Р23 (или Р32). Но если в этой же плоскости окажется еще хотя бы один пеленг (с ИП-2 или ИП-3), то однозначной идентификации не произойдет.

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

Заключение

Очевидно, что результаты идентификации по предложенным алгоритмам ИГП и ИПП наряду с другими условиями в значительной степени будут зависеть от условий их применения, т. е. от расположения ИП и характеристик используемых СИ, высоты и дальности размещения объекта регистрации, его формы и концентрации в нем элементов. Эти зависимости представляются достаточно сложными, имеющими стохастический характер, поэтому для их изучения необходимо разработать экспериментальную статистическую модель с алгоритмическим описанием причинно-следственных связей, что является дальнейшим развитием настоящей работы.

Литература

1. Molebny V., et al. Laser radar: historical prospective - from the East to the West // Opt. Eng. 2016. V. 56. № 3. P. 031220.

2. Малашин М.С., Каминский Р.П., Борисов Ю.Б. Основы проекти-рования лазерных локационных систем. М.: Высшая школа, 1983. 207 с.

3. Васильев С.В. и др. Перспективные методы и средства траекторных измерений. Казань: Новое знание, 2005. 128 с.

4. Natanson O.G., et al. Reflectometry in open and fiber mediums: technology transfer // Proc. SPIE. 2005. V. 5854. pp. 205-214.

5. Ефимов В.О. и др. Оптическая рефракция и модельные методы учета ее влияния на характеристики дальномерного тракта лазерного локатора слежения // Инженерный вестник Дона. 2017. №3. URL: ivdon.ru/magazine/archive/ n3y2017/4314.

6. Ефимов В.О. и др. Дисперсионный метод учета статической и динамической рефракции для дальномерного тракта лазерного локатора слежения // Инженерный вестник Дона. 2017. №3. URL: ivdon.ru/magazine/ archive/n3y2017/4315.

7. Ефимов В.О. и др. Методы компенсации влияния атмосферной турбулентности на точность измерения угловых координат в лазерных локаторах слежения // Научно-технический вестник Поволжья. 2017. №4. C. 69-71.

8. Morozov O.G., et al. Two-frequency scanning LFM lidars: theory and applications // Proc. of SPIE. 2002. V. 4539. pp. 158-168.

9. Ильин Г.И., Морозов О.Г., Польский Ю.Е. ЛЧМ-лидар с преобразованием частоты // Оптика атмосферы и океана. 1995. Т. 8. № 12. С. 1871.

10. Ильин Г.И., Морозов О.Г., Польский Ю.Е. Блок импульсной накачки лидарных комплексов // Оптика атмосферы и океана. 1995. Т. 8. № 5. С. 762.

11. Ильин Г.И. и др. Исследования ЛЧМ-лидара с преобразованием частоты // Оптика атмосферы и океана. 1997. Т. 10. № 2. С. 435.

12. Ильин Г.И., Морозов О.Г., Польский Ю.Е. Особенности построения электрооптических амплитудно-фазовых формирователей двухчастотного лазерного излучения для дифференциальных ЛЧМ-лидаров // Оптика атмосферы и океана. 1998. Т. 11. № 5. С. 513.

13. Ильин Г.И., Морозов О.Г., Польский Ю.Е. Применение амплитудно-фазового преобразования частоты лазерного излучения для создания специальных схем ЧМ-лидаров // Оптика атмосферы и океана. 1999. Т. 12. № 4. С. 360.

Аннотация

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

Ключевые слова: беспилотный летательный аппарат, группа объектов, идентификация пеленгов, горизонтальный пеленг, пространственный пеленг, критерии идентификации, ложная засечка.

Размещено на Allbest.ru

...