УДК 519.688
А. В. Швецов1, А. А. Ложкин1, И. А. Юхновский1, Р. Р. Хайрулин1, В. В. Бенсон1,
1 2 А. Н. Бебенин , И. Е. Волкова
ПРОГНОЗИРОВАНИЕ ДВИЖЕНИЯ ЦЕЛИ ПО ДАННЫМ ТЕЛЕАВТОМАТА
СОПРОВОЖДЕНИЯ
ОАО «ЦНИИ Буревестник»1,
Нижегородский государственный технический университет им. Р. Е. Алексеева2
Рассмотрен вопрос о прогнозировании движения цели по данным, полученным с телеавтомата сопровождения. Алгоритм на основе метода наименьших квадратов позволяет предсказывать поведение целей, время полета снаряда до которых не превышает нескольких секунд, а также удаленных, но слабо маневрирующих целей.
Ключевые слова: зенитная артиллерия, прогнозирование, метод наименьших квадратов.
Достоверное прогнозирование движения цели представляет значительный интерес при разработке систем управления огнем зенитной артиллерии. Ошибки прогнозирования являются одними из основных и непосредственно влияют на точность стрельбы и, следовательно, эффективность систем.
Задача прогнозирования движения цели сводится к определению закона движения цели по данным наблюдения за ней в течение некоторого интервала времени и предсказания ее положения в пространстве спустя определенный промежуток времени. Последний определяется временем полета снаряда до точки встречи с целью.
Стоит отметить, что большинство целей зенитной артиллерии являются управляемыми и, следовательно, способны изменять закон своего движения в процессе полета и в том числе за время полета снаряда до цели. Поскольку, тем не менее, точка встречи рассчитывается до выстрела, то учесть эти изменения для неуправляемых снарядов не представляется возможным. Поэтому при постановке задачи мы будем исходить из того, что за прогнозируемое время закон движения цели не изменяется, а движение цели полностью определяется ее поведением в предшествующие моменты времени (интервал наблюдения).
Такое предположение, по крайней мере, должно быть справедливо для прогноза на небольшие промежутки времени вперед (порядка нескольких секунд) и для достаточно слабо маневрирующих целей. Последнее обстоятельство связано еще и с тем, что воздушные цели, движущиеся на больших скоростях, достаточно инерционны - их маневрирование ограничено максимально допустимыми перегрузками.
Для задач зенитной артиллерии удобно работать в системе координат местного горизонта, поскольку наблюдения за целью происходят именно в ней, а не в системе криволинейных (географических) координат. Задача об определении установок стрельбы (0, у, ¿) по движущимся целям (рис.1) может быть сведена к решению системы уравнений:
х(0, у, *) - хс (г) = 0
<у(0, У, *) - Ус (г) = 0, (1)
2(0, у, г) - 2с (г) = 0
где х, у, 1 - координаты центра масс снаряда в полете; хс, ус, 1С - координаты цели. Здесь используется стартовая система координат Оху1, ось Оу которой направлена по нормали к поверхности земного эллипсоида в точке расположения орудия и, соответственно, плоскость Ох1 является плоскостью местного горизонта. Система уравнений (1) должна быть решена относительно неизвестных установок стрельбы (0, у, ¿) с заданной точностью, где 0 - угол
© Швецов А. В., Ложкин А. А., Юхновский И. А., Хайрулин Р. Р., Бенсон В. В., Бебенин А. Н., Волкова И. Е., 2014.
между вектором начальной скорости снаряда и и плоскостью горизонта, у - угол между проекцией скорости на плоскость углов стрельбы, t - время полета снаряда до цели. Зависимость координат центра масс снаряда r(t) = (x(t), y(t), z(t)) при заданных значениях углов наведения орудия (0, у) является решением системы уравнений движения твердого тела в пространстве [1, 2], в которых учтены аэродинамические силы и моменты, действующие на снаряд в полете, сила тяжести и сила Кориолиса, связанная с суточным вращением Земли.
Данные о положении цели в момент времени ti представлены двумя углами 0ci и yci, которые приходят с телеавтомата сопровождения цели, и расстоянием до цели rci, измеряемом лазерным дальномером. В реальности эти приборы всегда обладают некоторой погрешностью измерений, поэтому если положить, что rci, 0ci и yci - истинные значения в момент времени ti, то
Д r . = r . — ~ . Cl 'ci 'Cl '
Д0 ci = 0 ci — ~ci ' ДШ ci = Ш ci — V Ci '
(2)
- соответствующие ошибки измерений. А радиус-вектор цели в системе координат местного горизонта, измеряемый в момент времени имеет следующие компоненты:
х . = r . cos 0 . cos ш . , ci ci ci v ci '
y ■ = r .sin0 . , Jci ci ci '
z . = — r . cos 0 .sin ш . . ci ci ci v ci
(3)
На рис. 1 начало координат О совпадает с орудием. Радиус-вектор г = (х, у, £) характеризует положение центра масс снаряда в пространстве, тогда как радиус-вектор гс = (хс, ус, гс) - положение цели. Вектор и - вектор скорости центра масс снаряда в начальный момент времени.
В качестве исходных данных для прогнозирования положения цели используются координаты цели (хсг-, уСг, zc¡) в последовательные N моменты времени (/) от до tf. Закон движения цели аппроксимируется полиномиальной зависимостью
гс ^) = ТРп (/) + Ют (/) + кОк (/), (4)
где I, _/, к - орты системы координат местного горизонта (рис.1). Как уже отмечалось, предполагается, что цель в течение достаточного для ее поражения промежутка времени после момента f движется по тому же закону гс(/), что и в интервале времени от до tf. Отметим, что порядки аппроксимирующих полиномов Рп(/), Qm(t), 0к(/) в общем случае не совпадают. Каждый из полиномов Рп(/), Qm(t), 0$) определяется соответственно (п+1), (т+1) и (к+1) коэффициентами
Рис. 1. Система координат местного горизонта
к
n m к
Pn (0= 2 a/ , Qm (t) = Z ь/ , Gk (0= 2
c t i
(5)
* = 0 * * = 0 * * = 0 * Для прогнозирования положения цели в настоящей работе используется метод наименьших квадратов [3, 4, 6]. Последний основан на отыскании минимума функционала
£(а„,...а ,Ь„,...Ь ,с„,...с, )= ^г ((■)-г ■)2 = \ 0 п 0 т 0 к> ^сч' а'
I = 1
= Х((Т (,)-.с,)2 (',)-У*)2 4к(',)-Ъ)2)
I = 1
(6)
Данный функционал минимален для набора коэффициентов а0,...ап, Ъ0,...Ът, с0,...ск, являющихся решением системы (п+т+к+3) уравнений
-= 0, I = 0,...п
да.
I
-= 0, I = 0,...т ,
дЬ.
I
-= о, I = 0,...к,
дс.
I
(7)
которая приводится к системе линейных уравнений относительно неизвестных коэффициентов полиномов:
п
I а I = 0
т
I Ь I = 0
N
Л
I £+1
^ = 1
' N I х\ I = 1
N
= I х Л-!, 1 = 0,...п
. .. а I I = 1
I +1
N
= I У л-1, 1 = 0,...т,
. /а I I = 1
к ( N
Л
I с I I = 0 1 {I = 1
I +1 _
N
= I г Л], 1 = 0,...к. . .. а I I = 1
(8а) (8б) (8в)
Рис. 2. Траектория цели в период наблюдения [¿0, $ (сплошная линия) и прогнозируемая траектория цели (пунктирная линия)
Каждая из систем уравнений (8а), (8б), (8в) может быть разрешена независимо. Поскольку число определяемых коэффициентов полиномов и, соответственно, число уравнений (8а)-(8в) обычно невелико, то для решения системы уравнений используется метод Гаусса [5]. Так как координаты цели определены неточно (2), то рассчитанные в (8а)-(8в) коэффициенты а/, Ъ/, С1 характеризуются определенными ошибками относительно истинных значений.
По этой причине очевидно, что необходимо ограничиваться полиномами небольшой степени, так как члены высокого порядка ait1 (l >> 1) быстро возрастают с ростом t и соответственно увеличивают ошибку прогнозирования положения цели в пространстве. Разумеется жестко задавать порядки аппроксимирующих полиномов Pn(t), Qm(t), Gk(t) не имеет смысла, поскольку заранее неизвестно по какому закону будет двигаться цель. С этой целью оптимальные для аппроксимации порядки полиномов должны быть определены расчетным путем по накопленным в результате наблюдения за целью данным - (xci, yci, zci), i = 1,.. .N.
В настоящей работе используются полиномы наинизших степеней, для которых выполняется условие:
\Pn fe)-xc\ <8, \Qm {t,)-yc\ <8, \Gk {t,)-zc\ <8, Vt, e|to, if J, (9)
Порядки аппроксимирующих полиномов определяются путем последовательного решения систем уравнений (8а), (8в) и проверки условия (9).
На рис. 2 представлен результат прогнозирования движения цели. Наблюдение за целью происходит в течении 4 с. (сплошная линия на рис.2). Истинные координаты цели изме-
~ о ~
няются по закону хс = -7000 + 200 t, yc = 5000 - 10t, zc = 0. Для расчетов использовались за-
шумленные данные xci = xci + A xci, yci = yci + A yci, zci = zci + A zci , где A xci, A yci, A zci - случайные ошибки измерения координат с равномерным распределением в интервале ( -s/2, 8/2 ). Прогноз осуществляется на промежуток времени 8 с. вперед (пунктирная линия на рис.2). Значения 8 выбирались в интервале от 1 до 10 м. Ошибки определения координат цели относительно истинных составили значения меньшие 8 для координат xc и zc и величину порядка 48 для координаты yc. Это можно считать хорошим результатом, с учетом того, что время прогноза превышало время наблюдения в два раза.
Следует отметить, что описанный в работе алгоритм может успешно применяться для прогнозирования движения слабо маневрирующих целей, находящихся на значительных дистанциях, а также целей, испытывающих значительные перегрузки, но время полета снаряда до которых не превышает 3-4 с. В дальнейшем для увеличения точности прогноза планируется применить фильтрацию Калмана [6].
Библиографический список
1. ГОСТ В 24288-80. - М., Госстандарт, 1980.
2. McCoy, R. L. Modem Exterior Ballistics / R. L. McCoy - Atglen, Schiffer Military History, 2012.
3. Линник, Ю. В. Метод наименьших квадратов и основы математико-статистической теории обработки наблюдений / Ю. В. Линник. - М., 1962.
4. Зайдель, А. Н. Ошибки измерений физических величин / А. Н. Зайдель. - Л.: Наука, 1974.
5. Бахвалов, Н.С. Численные методы / Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков. - М.: Бином. Лаборатория знаний, 2011.
6. Фарина, А. Цифровая обработка радиолокационной информации / А. Фарина, Ф. Студер. -М.: Радио и связь, 1993.
Дата nocmyrnenun в редакцию 11.12.2014
A. V. Shvetsov, A. A. Lozhkin, I. A .Juhnowski, R.R. Khairulin,V. V. Benson,
A. N. Bebenin, I. E. Volkova
PREDICTION OF TARGET MOTION BY DATA OF AN OPTICAL TRACKING SYSTEM
Nizhny Novgorod state technical university n.a. R. E. Alexeev, Joint Stock Company Central Research Institute "BUREVESTNIK"
The forecasting algorithm of target motion is discussed for the case of data obtained by an optical tracking system. The algorithm based on the method of least squares allows to predict a target behavior if the time flight of a projectile to the target does not exceed a few seconds as well as for the case of slowly maneuvering remote targets.
Key words: flak, forecasting algorithm, least-squares method.