Способ определения кинематических параметров движения летательного аппарата
METHOD OF DETERMINING KINEMATIC PARAMETERS OF AIRCRAFT MOVEMENT

Область техники
Изобретение относится к области навигации и может быть использовано для определения угловых и пространственных координат, а также скоростей и ускорений летательного аппарата, например спутника.
Уровень техники
Известен способ определения высоты спутника по звездам, состоящий из каталога, содержащего последовательно соединенные звездный датчик, сигнальный процессор, который служит для обработки входного сигнала, полученного от звездного датчика, блок ранжирования звезд по яркости, блок идентификации звезд, к другому входу которого подключен блок памяти, в котором хранится каталог координат навигационных звезд (патент США №4680718 Method and apparatus of determining an attitude of a satellite (МПК B64G 1/36; G01C 21/24; G01S 3/78; G01S 3/782; G01S 3/785; G01S 5/16; G01V 8/10; G01S 3/786; (IPC1-7): G06F 7/56, опубл. 14.07.1987)).
Недостаток его состоит в том, что устройство по способу неспособно вычислять линейную скорость и ускорение спутника.
Известно устройство для определения ориентации аппарата по звездам, состоящее из блока памяти каталога навигационных звезд, звездного датчика, сигнального процессора, блока ранжирования сигнала по яркости, блока идентификации звезд, блока определения координат ориентации, блока памяти каталога перекрывающихся секторов, блока сравнения предварительных координат и блока выборки звезд и формирования рабочего каталога (полезная модель РФ №23979 УСТРОЙСТВО ДЛЯ ОПРЕДЕЛЕНИЯ ОРИЕНТАЦИИ КОСМИЧЕСКОГО АППАРАТА ПО ЗВЕЗДАМ (МПК G01J 1/20, G01C 21/24, опубл. 20.07.2002).
Недостаток его состоит в том, в нем также отсутствует возможность вычисления линейной скорости и ускорения летательного аппарата.
Наиболее близким техническим решением является метод обнаружения и определения пространственного местоположения перемещающихся воздушных объектов искусственного происхождения посредством измерения интенсивности изотропного реликтового излучения, включающий процесс последовательного измерения интенсивности изотропного реликтового излучения, приходящего в точку произведения метода, путем последовательного сканирования измерительными приборами небесной сферы, настроенным на частоту изотропного реликтового излучения (изобретение РФ №2563320 МЕТОД ОБНАРУЖЕНИЯ И ОПРЕДЕЛЕНИЯ ПРОСТРАНСТВЕННОГО МЕСТОПОЛОЖЕНИЯ ПЕРЕМЕЩАЮЩИХСЯ ВОЗДУШНЫХ ОБЪЕКТОВ ИСКУССТВЕННОГО ПРОИСХОЖДЕНИЯ ПОСРЕДСТВОМ ИЗМЕРЕНИЯ ИНТЕНСИВНОСТИ ИЗОТРОПНОГО КОСМИЧЕСКОГО РЕЛИКТОВОГО ИЗЛУЧЕНИЯ, МПК: G01C 21/02, опубл. 20.09.2015).
Недостатки его заключаются в том, что метод применим только в воздушном пространстве, для определения координат и вектора скорости воздушных объектов необходимо использовать несколько измерительных приборов, которые разнесены пространственно, метод может быть реализуем только в условиях прямой видимости воздушного объекта со стороны измерительных приборов, метод в значительной степени зависит от состояния атмосферы и не работает в условиях ливневого дождя, обильного снегопада, в области грозовых туч и при искусственной засветке атмосферы.
Раскрытие изобретения
Задача предлагаемого изобретения заключается в расширении области применимости способа (метода) с использованием реликтового излучения и уменьшении факторов, оказывающих возможное влияние на точность способа.
Задача решается и технический результат достигается тем, что способ определения кинематических параметров движения летательного аппарата включает в себя процедуру позиционирования реализующего способ устройства регистрации, процедуру измерения интенсивности реликтового излучения, вычисление координат движущегося объекта. При этом для определения кинематических параметров движения летательного аппарата, а именно: ориентации, положения, скорости и ускорения летательного аппарата в пространстве, вводят базу данных, хранящую в себе карту реликтового излучения и модель движения Земли, измерительные приборы устанавливают на летательном аппарате, устройства сканирования для определения температуры в заданном угловом диапазоне направляют в области пространства, где на карте реликтового излучения максимальны производные по углу и по температуре, измеряют спектральные плотности мощности излучения не менее чем в трех направлениях и не менее чем в трех спектральных диапазонах длин волн, вычисляют температуры реликтового излучения и при сопоставлении данных, получаемых с устройств вычисления температуры реликтового излучения с данными, содержащимися в базе данных, определяют кинематические параметры движения летательного аппарата: угловые координаты, скорости и ускорения, а также линейные скорости и ускорения.
При реализации способа определения кинематических параметров движения летательного аппарата используют установленные на борту ЛА измерители спектральной плотности мощности, установленные в устройстве сканирования в заданном угловом диапазоне, управляемые блоками управления, блоки вычисления температуры реликтового излучения, вычислитель и базу данных. База данных хранит в себе измеренную карту реликтового излучения. Сопоставление данных, получаемых с устройств вычисления температуры реликтового излучения с данными, содержащимися в базе, позволяют вычислять кинематические параметры движения летательного аппарата: угловые координаты, скорости и ускорения, а также линейные скорости и ускорения.
Перечень фигур
На фиг. 1 представлена структурная блок-схема устройства, реализующего способ определения кинематических параметров движения летательного аппарата.
На фиг. 2 представлен график точности в определении угловой координаты летательного аппарата, в зависимости от направления измерителя спектральной плотности мощности.
На фиг. 3 представлен график точности определения линейной скорости летательного аппарата в зависимости от направления измерителя спектральной плотности мощности.
1
1
2
2
3
1
На фиг. 4 представлена модель дипольной составляющей реликтового излучения, представленная в виде линий уровня. Отмечены три направления датчиков спектральной плотности мощности, для которых ведется расчет точности способа (b=41°, l=131°; b=49°, l=179°; b=-11°, l=119°).
Осуществление изобретения
Способ определения кинематических параметров движения летательного аппарата осуществляется в устройстве, состоящим из n каналов, в состав которых входят: n измерителей спектральной плотности мощности на различных длинах волн 1.1, 1.2…1.n, но не меньше трех, устройства сканирования в заданном угловом диапазоне 2.1, 2.2…2.n, блоки вычисления температуры реликтового излучения 3.1, 3.2…3.n, блоки управления сканирующими устройствами 6.1, 6.2…6.n, базы данных 4, вычислитель 5.
Предлагаемый способ реализуют следующим образом.
Устройство, реализующее способ определения кинематических параметров летательного аппарата на Земной поверхности, настраивают следующим образом: устройства сканирования направляют в области пространства, где на карте реликтового излучения максимальны производные по углу и по температуре. Таким образом обеспечивается минимальная ошибка в определении углового положения летательного аппарата, угловой скорости и линейной скорости. Конкретная настройка зависит от текущего положения Земли и погодных условий. Спектральная плотность мощности реликтового излучения, измеренная на разных длинах волн, с помощью 1.1, 1.2…1.n, оси которых располагаются не в одной плоскости, поступает на вход вычислителей 3.1, 3.2…3.n. Вычислители 3.1, 3.2…3.n преобразуют полученные на разных длинах волн значения к кривой изменения спектральной плотности мощности от длины волны, соответствующей излучению абсолютно черного тела с известной температурой. На основе полученной кривой на выходе 3.1, 3.2…3.n формируются значения температур реликтового излучения в заданном направлении. Вычисленные значения температур реликтового излучения поступают на вычислитель 5. На основе сравнения данных, пришедших с 3.1, 3.2…3.n, и данных, хранящихся в базе данных 4, формируются угловые координаты летательного аппарата и радиус-вектор в земной системе координат. Скорость летательного аппарата вычисляется на основе сравнения температуры реликтового излучения на карте и измеренной температурой на выходе вычислителей 3.1, 3.2…3.n. В случае покоя летательного аппарата измеренные температуры и температуры на карте будут равны. В случае движения летательного аппарата в каком-либо направлении изменение температуры в этом направлении будет пропорционально скорости движения. Ускорения определяют как изменение скоростей во времени. Вычислитель 5 имеет два выхода: на одном выходе формируется вектор кинематических параметров летательного аппарата, на другом выходе - управляющее воздействие, поступающее на вход блоков управления сканирующими устройствами 6.1, 6.2…6.n. Блоки управления сканирующими устройствами 6.1, 6.2…6.n осуществляют поворот устройств сканирования в заданном угловом диапазоне 2.1, 2.2…2.n согласно требуемой программе по минимизации сигнала рассогласования между текущей ориентацией измерителей спектральной плотности мощности и их ориентацией, при котором ошибка измерения углового положения летательного аппарата минимальна. Устройства сканирования в заданном угловом диапазоне 2.1, 2.2…2.n содержат в себе измерители спектральной плотности мощности 1.1, 1.2…1.n.
-6
-3
В основе способа определения кинематических параметров движения летательного аппарата лежит свойство реликтового излучения - дипольная анизотропия. Это свойство интерпретируется как следствие эффекта Доплера, возникающего при движении Солнечной системы относительно реликтового фона со скоростью примерно 370 км/с в сторону созвездия Льва [1]. Способ позволяет однозначно определить угловое положение, скорость, ускорение и радиус-вектор аппарата в пространстве на основе измерений спектральной плотности мощности по нескольким направлениям. Спектр реликтового излучения схож со спектром излучения абсолютно черного тела. В настоящее время определение температуры T+dT реликтового излучения в выбранном направлении возможно с точностью dT~10, а величина, при которой заметна дипольная анизотропия излучения, вносит свой вклад T'~10 [2], то есть существенное отличие на 3 порядка.
Температура излучения в выбранном направлении для тела, движущегося со скоростью υ, равняется:
,
0
где T - температура в выбранном направлении, T - температура реликтового излучения в ИСО, которая покоится относительно излучения, θ - угол между вектором υ и направлением наблюдения,
Возможную точность определения угла dθ можно рассчитать с учетом того, что точность измерения температуры реликтового излучения на данный момент составляет dT=0.000001 K.
Для определения dθ найдем производную Tabs по θ:
Выразим dθ:
0
Скорость движения Солнца в направлении созвездия Льва составляет υ=370 км/с, T=2.72548 К, а угол θ=(0; 360)°
Полученная зависимость точности определения угла dθ от θ имеет вид, представленный на фиг. 2.
min
min
Минимальная ошибка определения углового положения вдоль одного направления dθ=58,2ʺ достигается, если θ=90°, 270°. В интервале θ=(5,8…174,2)°∧(185,8…354,2)° ошибка не превышает 10dθ
Для нахождения ошибки определения скорости летательного аппарата dυ найдем
Отсюда dv:
Полученная зависимость dv от θ имеет вид, представленный на фиг. 3.
min
min
Минимальная ошибка определения скорости dv=0.11 км/с достигается, если θ=0°, 180°, 360°. В интервале θ=(0…84,38)°∧(95,62…264,38)°∧(274,38°…360)° ошибка не превышает 10dv.
Для оценки точности предлагаемого способа было произведено моделирование работы реализующего способ устройства с тремя датчиками спектральной плотности мощности. В качестве модели дипольной составляющей реликтового излучения была взята модель, предложенная в [3].
x
y
z
В галактических координатах распределение дипольной анизотропии характеризуется тремя амплитудами ΔT, ΔT и ΔT:
x
y
z
ΔT(l,b)=ΔTcoslcosb+ΔTsinlcosb+ΔTsinb,
l - галактическая широта, b - галактическая долгота (галактическая система координат)
Результирующие температуры реликтового излучения, воссозданные по этой модели, представлены в виде линий уровня на Фиг. 4. По осям отложена галактическая широта, долгота, для каждой линии уровня характерна температура в мкК, обозначенная на графике.
1
1
2
2
3
1
Для расчета точности устройства, реализующего предлагаемый способ, выберем три направления на карте реликтового излучения (b=41°, l=131°; b=49°, l=179°; b=-11°, l=119°).
Суммарная погрешность в позиционировании ЛА будет равна сумме погрешностей по каждому из направлений. Погрешности в определении позиционирования ЛА для каждого из этих направлений:
1
Δ=0,02°,
2
Δ=0,026°,
3
Δ=0,022°
Полученная точность Δ=2'22ʺ не является предельной для предлагаемого способа, ввиду выбора трех направлений неоптимальным образом. Также точность предлагаемого способа может быть увеличена использованием большего количества - n троек датчиков спектральной плотности мощности реликтового излучения.
Список использованной литературы
[1] Чернин А.Д. Звезды и физика, М.: Наука, 1984, с. 152-153.
[2] Насельский П.Д., Новиков Д.И., Новиков И.Д. Реликтовое излучение Вселенной. - М.: Наука, 2003. - 390 с.
[3] Kogut A., Banday A.J., Bennett C.L. et al, 1996а. ApJ. V. 470. P. 653.
[4] Kogut A., Banday A.J., Bennett C.L. et al. II 1996b. ApJ. V. 464. L29.