способ определения линейной скорости летательного аппарата на основе использования нерадиолокационной информации

Формула изобретения

Способ определения линейной скорости летательного аппарата на основе использования нерадиолокационной информации, заключающийся в том, что с помощью оптико-электронной системы по двум последовательно полученным изображениям летательного аппарата определяют значения углов рыскания летательного аппарата _k, _k-1 и углов тангажа летательного аппарата _k, _k-1 в k-й и (k-1)-й моменты времени, соответственно, используя выражения:





где и - углы пространственной ориентации летательного аппарата в лучевой системе координат радиоэлектронной системы сопровождения в k-й момент времени, информация о которых поступает от телевизионной системы, входящей в состав оптико-электронной системы;

_гk и _вk - углы пеленга летательного аппарата в горизонтальной и вертикальной плоскостях, соответственно, в k-й момент времени, информация о которых поступает от теплопеленгатора, входящего в состав оптико-электронной системы,

на основе вычисленных значений углов рыскания летательного аппарата _k, _k-1 и углов тангажа летательного аппарата _k, _k-1, в k-й и (k-1)-й моменты времени, соответственно, определяют угловую скорость изменения угла рыскания летательного аппарата _k и угловую скорость изменения угла тангажа летательного аппарата _k в k-й момент времени, используя выражения:





где - интервал обработки изображений летательного аппарата оптико-электронной системой;

на основе вычисленных значений угловой скорости изменения угла рыскания летательного аппарата _k и угловой скорости изменения угла тангажа летательного аппарата _k в k-й момент времени с помощью функциональных статистических зависимостей (функций регрессии) определяют проекции вектора линейной скорости летательного аппарата на горизонтальную плоскость V_гk и вертикальную плоскость V_вk в k-й момент времени:





с использованием которых определяют значение линейной скорости летательного аппарата в k-й момент времени:

Описание изобретения к патенту

Изобретение относится к области радиотехники, в частности, к радиоэлектронным системам измерения координат и параметров движения и может быть использовано в бортовых и наземных радиоэлектронных системах сопровождения (РЭСС).

Известен способ измерения радиальной составляющей линейной скорости летательного аппарата (ЛА) с использованием доплеровского смещения частоты (см. Финкельштейн, М.И. Основы радиолокации. - М.: Радио и связь, 1983. - 386 с).

Сущность данного способа состоит в следующем. Определение радиальной составляющей линейной скорости ЛА V_p основано на измерении радиолокационной станцией (РЛС) доплеровского частотного сдвига f_д принимаемых сигналов относительно колебаний, генерируемых передатчиком РЛС:

где

- рабочая длина волны РЛС; f₀ - частота отраженных от ЛА колебаний; f - частота излучаемых РЛС колебаний.

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

Известен способ определения путевой скорости полета воздушной цели в наземной РЛС (Грачев В.В., Кейн В.М. Радиотехнические средства управления воздушным движением. - М.: Транспорт, 1975. - С.281-285).

Сущность способа состоит в следующем. Наземная РЛС с известными координатами, работающая в импульсном режиме, осуществляет круговое сканирование узким лучом диаграммы направленности антенны в горизонтальной плоскости. Измеренные РЛС координаты воздушной цели в полярной системе координат «азимут-дальность» переводятся в координаты (X, Y) прямоугольной системы координат. Для измерения путевой скорости воздушной цели вычисляют координаты (x_k-1 , y_k-1) в (k-1)-й момент времени и координаты (x _k, y_k) воздушной цели в k-й момент времени, а модуль путевой скорости определяется из выражения:

где t_k и t_k-1 - k-й и (k-1)-й моменты времени соответственно; (t_k -t_k-1) - период измерения координат (x_k-1 , y_k-1) и (x_k, y_k).

Недостатком данного способа является зависимость точности определения путевой скорости полета воздушной цели от точности вычисления координат (x_k-1, y_k-1) и (x_k , y_k), а также величины периода измерения координат (t_k-t_k-1).

Также известен способ определения путевой скорости воздушной цели в наземной радиолокационной станции (РЛС) (Князев И.Н., Князев Р.И. Способ определения путевой скорости воздушной цели в наземной радиолокационной станции. - Патент № 2273032 от 27.03.06. Бюллетень № 9. - 7 с.).

Сущность предлагаемого способа определения путевой скорости воздушной цели в наземной РЛС заключается в том, что измеряют частоту Доплера отраженных от движущейся цели сигналов в наземной РЛС f_дРЛС, а также частоту Доплера f_дR отраженных сигналов в дополнительной приемной позиции R, разнесенной в пространстве относительно наземной РЛС на базовое расстояние R_B, измеряют угол между направлениями «дополнительная приемная позиция R - цель Ц» и «дополнительная приемная позиция R - РЛС», угол между направлениями «РЛС - цель Ц» и «РЛС - дополнительная приемная позиция R», вычисляют бистатический угол =180°-( + ), а путевая скорость полета воздушной цели определяется по формуле:

где - рабочая длина волны, используемая в наземной РЛС; f _дРЛС - частота Доплера, измеряемая в наземной РЛС; f _дR - частота Доплера, измеряемая в дополнительной приемной позиции R; - бистатический угол между направлениями «РЛС - цель Ц» и «дополнительная приемная позиция R - цель Ц».

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

Ни один из рассмотренных способов не может быть принят за прототип предлагаемого способа определения линейной скорости ЛА.

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

Сущность предлагаемого способа измерения линейной скорости ЛА V_k в k-й момент времени на основе использования нерадиолокационной информации заключается в том, что с помощью оптико-электронной системы (ОЭС), входящей в состав РЭСС, по двум последовательно полученным изображениям ЛА определяют значения углов рыскания ЛА _k, _k-1, и углов тангажа ЛА _k, _k-1, в k-й и (k-1)-й моменты времени, соответственно, используя выражения:

где и - углы пространственной ориентации ЛА в лучевой системе координат РЭСС в k-й момент времени, информация о которых поступает от телевизионной системы (ТС), входящей в состав ОЭС; _гk и _вk - углы пеленга ЛА в горизонтальной и вертикальной плоскостях, соответственно, в k-й момент времени, информация о которых поступает от теплопеленгатора (ТП), входящего в состав ОЭС; на основе вычисленных значений углов рыскания ЛА _k, _k-1 и углов тангажа ЛА _k, _k-1, в k-й и (k-1)-й моменты времени, соответственно, определяют угловую скорость изменения угла рыскания ЛА _k и угловую скорость изменения угла тангажа ЛА _k в k-й момент времени, используя выражения:

где - интервал обработки изображений ЛА ОЭС;

на основе вычисленных значений угловой скорости изменения угла рыскания ЛА _k и угловой скорости изменения угла тангажа ЛА _k в k-й момент времени с помощью функциональных статистических зависимостей (функций регрессии) определяют проекции вектора линейной скорости ЛА на горизонтальную плоскость V_гk и вертикальную плоскость V_вk в k-й момент времени:

с использованием которых определяют значение линейной скорости ЛА в k-й момент времени:

Сущность способа поясняется следующим. Линейную скорость ЛА, движущегося по криволинейной траектории, можно описать выражением (Яворский Б.М., Детлаф А.А. Справочник по физике. - М.: Наука, 1985. - 512 с.):

где _k - абсолютное значение угловой скорости перемещения ЛА; r_k - радиус окружности, по дуге которой движется ЛА в k-й момент времени.

Так как ЛА может маневрировать и в горизонтальной, и в вертикальной плоскостях, то следует определять значения проекций вектора линейной скорости ЛА на горизонтальную V_гk и вертикальную V_вk плоскости:

где _гk, _вk - абсолютные значения угловых скоростей перемещения ЛА в горизонтальной и вертикальной плоскостях, соответственно, в k-й момент времени; r_гk, r_вk - проекции радиуса окружности, по дуге которой движется ЛА в k-й момент времени, на горизонтальную и вертикальную плоскости, соответственно.

Абсолютные значения угловых скоростей перемещения ЛА в горизонтальной и вертикальной плоскостях являются ни чем иным, как абсолютными значениями угловых скоростей изменения угла поворота _k и угла наклона _k траектории ЛА _k и _k в k-й момент времени, соответственно:

В РЭСС нет возможности измерять углы _k и _k, а, следовательно, и абсолютные значения угловых скоростей их изменения _k и _k. Однако есть возможность по получаемым от ОЭС в реальном масштабе времени изображениям сопровождаемого ЛА определять углы его рыскания _k и тангажа _k, а значит и абсолютные значения угловых скоростей их изменения _k и _k.

В силу малости углов атаки _k и скольжения _k при выполнении большинства маневров современными ЛА правомерно допущение о совпадении направления вектора линейной скорости ЛА и его продольной оси. Поэтому принимается, что _k= _k, _k= _k и _k= _k, _k= _k.

Тогда значения проекций вектора линейной скорости ЛА на горизонтальную плоскость V_гk и вертикальную плоскость V_вk в k-й момент времени определяются как:

Из выражений (17) и (18) видно, что при известных проекциях радиуса окружности, по дуге которой движется ЛА в k-й момент времени, на горизонтальную плоскость r_гk и вертикальную плоскость r_вk, проекции вектора линейной скорости ЛА на горизонтальную плоскость V_гk и вертикальную плоскость V_вk связаны детерминированной функциональной зависимостью с абсолютными значениями угловых скоростей изменения его угла рыскания _k и тангажа _k, соответственно.

Однако в РЭСС информация о параметрах вращательного движения сопровождаемого ЛА r _гk и r_вk отсутствует. этом случае предполагается, что значения V_гk и V_вk являются возможными значениями соответствующих случайных величин (СВ) V_г и V_в, а значения _k и _k возможными значениями соответствующих СВ и . Тогда связь СВ V_г и V_в с соответствующими СВ и можно описать статистическими функциональными зависимостями - функциями регрессии (Елисеева И.И., Юбзашев М.М. Общая теория статистики / под ред. чл.-корр. РАН И.И. Елисеевой. - М.: Финансы и статистика, 1995. - 368 с.):

где и - условные средние СВ V_г и V_в (средние арифметические значения СВ V_г и V_в при фиксированных значениях СВ и , соответственно).

Исследование вида функций регрессии (19) и (20) проводилось на основе сформированных по реальной полетной информации выборок СВ V_г, V _в, и . Кроме того, для исследования связи (19) осуществлялось разбиение выборок СВ V_г и в зависимости от условия: =0 (маневр ЛА осуществляется в горизонтальной, либо наклонной плоскостях) или 0 (маневр ЛА осуществляется по пространственной траектории). По аналогии, для исследования связи (20) проводилось разбиение выборок СВ V_в и в зависимости от условия: =0 (маневр ЛА осуществляется в вертикальной плоскости) или 0 (маневр ЛА осуществляется по пространственной траектории).

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

Один из вариантов программной реализации предлагаемого способа представлен на фигуре 1 при помощи логической схемы алгоритма функционирования программы определения линейной скорости ЛА на основе использования нерадиолокационной информации.

Логическая схема состоит из блока «Данные» I; блока «Подготовка» II; блоков «Процесс» III, VI-X; блоков «Решение» IV, V; блока «Запоминаемые данные» XI; блоков «Оперативное запоминающее устройство» XII, XIII.

В блоке I устанавливаются исходные значения параметров , ₀ и ₀. На вход блока III поступают значения , , _гk и _вk в k-й момент времени от внешних источников информации. Вычисление выражений, определяющих значения _k, _k, _k и тангажа _k в k-й момент времени производится в блоке III. В блоке IV в зависимости от выполнения условия _k=0 принимается решение о вычислении значения V _гk в k-й момент времени в блоке VI или блоке VII. В блоке V в зависимости от выполнения условия _k=0 принимается решение о вычислении значения V _вk в k-й момент времени в блоке VIII или блоке IX. Вычисления в блоке VI (блоке VII) и в блоке VIII (блоке IX) осуществляются параллельно. Результаты вычислений блока VI (блока VII) и блока VIII (блока IX) используются в блоке Х для вычисления значения V_k.

Результаты вычислений блока III и блока Х запоминаются в блоке XI для использования на следующем шаге вычислений в (k+1)-й момент времени (блок II). Вычисленное в блоке Х значение V_k в k-й момент времени запоминается в оперативном запоминающем устройстве (блок XII). Функционирование блоков III-XII осуществляется до момента сброса сопровождения (блок XIII).