В задачах пеленгового сопровождения объектов построение траектории наблюдаемого объекта ведется, как правило, по неполным данным, с помощью методов математического моделирования [3, 4, 8]. Задача построения траектории - это задача восстановления функции по ограниченному числу данных наблюдения. При пассивном сопровождении маневрирующих объектов необходимо одновременно рассматривать аппроксимирующие функции различной сложности, моделирующие разные типы маневра.
На сегодняшний день существует несколько математических методов, пригодных для реализации указанных алгоритмов, при этом каждый из них обладает рядом существенных ограничений, что сильно усложняет их практическую реализацию [6, 1, 2].
Применение методов системного анализа позволило нам сформулировать новую постановку задачи определения КПДО для нелинейно движущихся объектов, позволяющую снять большую часть приведенных ограничений.
В основу постановки задачи положено нахождение функции координат объекта от времени. Это позволяет перейти от сложившейся практики непосредственного определения параметров движения к нахождению функции траектории произвольного вида, оптимально аппроксимирующей результаты дискретных наблюдений.
Очевидно, что любые кинематические параметры движения объекта, являющиеся по сути производными функциями координат от времени, могут быть легко получены из данной функции.
Постановка задачи
Условия задачи подразумевают, что объект наблюдения (ОН) движется в двумерном пространстве по гладкой траектории. Траектория движения объекта представляет собой функцию вектора координат от времени и заключает в себе всю полноту информации о положении, параметрах и характере движения наблюдаемого объекта.
Наблюдатель, осуществляющий слежение за объектом, движется по гладкой траектории, также представляющей собой функцию вектора координат от времени. Данную функцию будем полагать известной и адекватной реальным положению и параметрам движения наблюдателя.
В дискретные моменты времени t1, выбранные на равномерной сетке с началом координат t0 и шагом Δt, наблюдатель осуществляет измерение угла пеленга на объект наблюдения P(t). Под углом пеленга понимается угол между направлением на север и направлением на объект наблюдения.
Наблюдение угла пеленга производится с некоторой ошибкой, называемой ошибкой измерений и считающейся распределенной по нормальному закону распределения [4].
Необходимо по данным наблюдения восстановить траекторию движения цели с заданной точностью.
Анализ постановки задачи
Будем рассматривать координаты объекта наблюдения и поступающие в систему данные о результатах измерений пеленгов и координатах наблюдателя во времени.
Объект наблюдения движется непрерывно, в то время как наблюдения осуществляются дискретно. Тогда задача построения траектории объекта наблюдения представляет собой задачу построения функции, аппроксимирующей некие дискретные значения, полученные в результате обработки наблюдений.
Следует заметить, что мы не имеем возможности наблюдать координаты движущегося объекта непосредственно, а регистрируем лишь углы пеленга, являющиеся функцией одновременно координат наблюдателя и наблюдаемого объекта.
В математическом смысле данная задача представляет собой задачу восстановления функции по ограниченному числу данных наблюдения.
Для решения этой задачи требуется ввести функцию, реализующую связь данных наблюдения с оцениваемыми параметрами.
В большинстве существующих методов в качестве математической модели, определяющей данную связь, используются уравнения кинематики, а рассматриваемые в них кинематические параметры существенно ограничивают множество моделей движения наблюдаемого объекта [5, 7].
При рассмотрении процесса движения как функции траектории кинематические параметры движения, выступающие в качестве производных, исключаются из процесса решения и не оказывают влияние на решение задачи.
В этой ситуации связь между данными наблюдения и состоянием наблюдаемой системы осуществляется посредством задания уравнения прямой, проходящей через позиции наблюдателя и ОН и являющейся линией пеленга (рисунок 1).
Рисунок 1. Определение линии пеленга на ОН,
- координаты наблюдателя, - координаты ОН, P - угол пеленга на ОН
Уравнение прямой, в данном случае, имеет вид:
(1)
где - линейный угловой коэффициент.
Уравнение (1) является уравнением прямой, в практическом смысле представляющее собой уравнение линии пеленга на наблюдаемый объект.
Данное уравнение лежит в основе математической модели предложенного далее метода нахождения траектории объекта наблюдения.
Следующим шагом в построении математической модели будет выбор вида функции, максимально приближенно описывающей оцениваемые параметры.
Следует отметить, что функция траектории, в силу физической природы описываемого объекта, является гладкой и непрерывной. Для описания подобных функций в научной практике широко используются методы аппроксимации искомой функции полиномами.
Среди широкого семейства аппроксимирующих полиномов мы выбрали семейство ортогональных многочленов Чебышева, обладающие рядом важных для рассматриваемой задачи свойств.
Рассмотрение альтернативных систем ортогональных полиномов показало отсутствие существенного, с точки зрения оценки точности, различия в реализации предложенного метода на их основе.
Метод N-полиномов
Будем полагать, что наблюдатель движется по траектории, описываемой функциями .
Для нахождения функции траектории движения объекта построим аппроксимирующие функции координат от времени наблюдения , выражая их посредством линейную комбинацию ортогональных многочленов, например полиномов:
(2)
где:
- полиномы Чебышева 1-го или 2-го рода [8], либо полиномы Лежандра;
- приведенное время наблюдения.
Сетку дискретизации для каждого из измерений рассчитаем по формуле:
(3)
где:
ti- время i- того наблюдения;
t1- время первого наблюдения;
tm- время последнего наблюдения.
В ходе решения задачи мы получаем данные о положении наблюдаемого объекта в виде углов пеленга Pi, которые затем преобразуются в угловые коэффициенты ki для уравнения (1) следующим образом:
(4)
В силу вычислительных особенностей тригонометрических функций присутствующий в данных наблюдения углов пеленга шум оказывает неравномерное воздействие на точность вычислений при различных значениях аргумента функции.
При этом максимальное искажение углового коэффициента происходит при углах наблюдения, близких к 90°.
Для минимизации влияния шумовых возмущений для каждого уравнения наблюдения будем осуществлять тождественное преобразование координат, поворачивающее базовую систему координат таким образом, чтобы углы наблюдения в новой системе координат были близки к 0°.
Для этого осуществим поворот координат на угол α, который будем выбирать таким образом, чтобы . При этом угловой коэффициент ki в уравнении (1) примет вид:
(5)
Запишем значения базовых координат объекта наблюдения в виде полярных координат:
(6)
где Ri- расстояние между наблюдателем и объектом наблюдения.
Запишем повернутые координаты объекта как:
(7)
Осуществим тождественное преобразование для системы (7):
(8)
Подставив (6) в (8), получим выражение, описывающее координаты объекта наблюдения в системе повернутых координат:
(9)
Аналогичные преобразования выполним и для координат наблюдателя.
Переписав уравнение (1) для повернутой системы координат и сгруппировав известные члены в правую часть, получим:
(10)
Учитывая (5) и подставив (9) в (10), получим:
(11)
Обозначим
Для удобства запишем (11) как:
(12)
Подставив аппроксимирующие функции (2) в (12) для каждого из произведенных наблюдений, а затем, преобразовав полученную систему к матричному виду, получим матричное уравнение:
(13)
где
(14)
(15)
(16)
Здесь (14) представляет собой матрицу наблюдений, а (16) вектор столбец, содержащий координаты наблюдателя.
Решив систему (13) относительно X, мы найдем коэффициенты .
Подставив найденные коэффициенты (2), мы получим искомые функции, определяющие траекторию движения объекта наблюдения.
Применение системы ортогональных полиномов Чебышева позволяет вести наблюдения за объектом на равномерной временной сетке, при этом количество измерений может превышать 2*n, где n - степень полинома, аппроксимирующего траекторию объекта.
В данном случае необходимо привести матрицу A к квадратному виду. Для решения данной задачи воспользуемся широко применяемым методом наименьших квадратов, который позволит осуществить дополнительную фильтрацию измерений. В результате система (13) примет вид:
(17)
Дальнейшее решение системы осуществляется аналогично решению системы (13).
Сравнительный анализ точности определение параметров движения для исследуемых методов
Далее приведем основные результаты сравнительного моделирования ошибки определения дистанции до наблюдаемого объекта по методу N-пеленгов и Методу N- полиномов.
На рисунках 2-5 представлены результаты для различных типов движения объекта наблюдения.
Рисунок 2. Ошибка определения дистанции до ОН при равномерном прямолинейном движении.
Заштрихованная поверхность - метод N-пеленгов; Каркасная поверхность - метод N-полиномов
Рисунок 3. Ошибка определения дистанции до ОН при равноускоренном прямолинейном движении.
Заштрихованная поверхность - метод N-пеленгов; каркасная поверхность - метод N-полиномов
Рисунок 4. Ошибка определения дистанции до ОН при нелинейном движении.
Заштрихованная поверхность - метод N-пеленгов; каркасная поверхность - метод N-полиномов
Рисунок 5. Ошибка определения дистанции до ОН при равномерном движении с изменением курса.
Заштрихованная поверхность - метод N-пеленгов; Каркасная поверхность - метод N-полиномов
Распределение осей на графиках 2-5:
Ось N - порядковый номер дискретного наблюдения;
Ось S - максимальный уровень шума при наблюдении угла пеленга (в угловых минутах);
Ось E - ошибка определения дистанции (в % от дальности до ОН).
Для каждого, представленного на графике типа движения, расчет произведен по данным серии из 1000 вычислительных экспериментов.
Параметры экспериментов: моделирование шума определения угла наблюдения по нормальному закону с макс. уровнем шума 60', время наблюдения - 600 с, период дискретных наблюдений 15 с, начальная дистанция до наблюдаемого объекта 3000 м.
Рис. 2 показывает практическую эквивалентность исследуемых методов для случая равномерного прямолинейного движения ОН. Из приведенных на рис. 3-5 графиков видно, существенное улучшение точности определения дистанции при использовании метода N-полиномов для нелинейно движущегося объекта.
Обсуждение результатов
В данной работе автором предложен метод, позволяющий повысить точность определения параметров нелинейно движущихся объектов с использованием только угломерной информации.
Предложенный метод относится к классу геометрических методов и обладает значительной вычислительной простотой по сравнению с методами другой природы. Метод разработан с учетом основных факторов, препятствующих решению задачи, и ориентирован на учет нелинейного движения наблюдаемого объекта.
Метод позволяет:
- осуществлять адекватную оценку траекторий ОН различного уровня сложности, в том числе и нелинейных траекторий;
- осуществлять последующий анализ таких параметров как: вектор скорости, ускорения, скорости изменения ускорения ОН, путем анализа продуцируемой ими функции координат ОН от времени;
- осуществлять непрерывное решение, в не зависимости от параметров движения ОН;
- формировать более точные оценки КПДО по сравнению с применяемым на практике методом N-пеленгов в случаях нелинейного движения наблюдаемого объекта.
Рецензенты:
- Потапенко А. А., доктор физико-математических наук, профессор, заведующий кафедрой математики, ГОУ ВПО «Северо-Западный государственный заочный технический университет», г. Санкт-Петербург.
- Федорцов А. Б., доктор физико-математических наук, профессор, заведующий кафедрой физики, ГОУ ВПО «Северо-Западный государственный заочный технический университет», г. Санкт-Петербург.