www.mai.ru/science/trudy/
Труды МАИ. Выпуск № 91
УДК 629.3.05
Определение навигационных параметров беспилотного летательного аппарата на базе фотоизображения и инерциальных измерений
1 Ä 1 ÄÄ Iftftft
Антонов Д.А. , Жарков М.В. , Кузнецов И.М. ,
Лйййй 1 Ä Ä Ä Ä Ä
Лунев Е.М. , Пронькин А.Н.
1 Московский авиационный институт (национальный исследовательский университет), МАИ, Волоколамское шоссе, 4, Москва, A-80, ГСП-3, 125993, Россия Объединенная авиастроительная корпорация «Центр комплексирования», Авиационный переулок, 5, Москва, 125139, Россия *e-mail: [email protected] **e-mail: [email protected]. ***e-mail: im [email protected] ****e-mail: [email protected] *****e-mail: [email protected]
Аннотация
В статье предложен вариант комплексирования инерциальной навигационной системы (ИНС) и системы определения параметров ориентации и навигации на базе фотограмметрической обработки изображения искусственных наземных ориентиров (ФС) для использования на этапе посадки беспилотного летательного аппарата (БПЛА). Представлены математические модели погрешностей ИНС и фотограмметрической системы, необходимые комплексной обработки информации. Приводятся результаты имитационного моделирования для анализа работоспособности и основных характеристик предложенных технических решений.
Ключевые слова: беспилотный летательный аппарат, инерциальная навигационная система, фотограмметрия, комплексная обработка информации, система посадки.
Постановка задачи
Одним из наиболее активно растущих секторов авиационной отрасли является сектор беспилотных авиационных систем. При этом БПЛА самолетного типа средней дальности пользуются повышенным спросом у потенциальных заказчиков в связи с привлекательным соотношением стоимости и функциональных возможностей.
Пилотажно-навигационный комплекс (ПНК) БПЛА указанного класса должен отвечать не только жестким требованиям к точности определения пилотажно-навигационных параметров, но и к обеспечению характеристик доступности, непрерывности и целостности навигационного обеспечения полета. При этом на массогабаритные характеристики, энергопотребление и, как правило, стоимость бортового оборудования ПНК накладываются ограничения, продиктованные спецификой БПЛА этого типа.
Для преодоления обозначенных ограничений в качестве первичных
измерителей в составе ПНК БПЛА используются бюджетные и малогабаритные
датчики и системы, а для удовлетворения требований к характеристикам ПНК его
программное обеспечение функционирует на основе развитого математического
аппарата комплексной обработки первичной информации [1]. При этом одним из
наиболее сложных этапов полета с точки зрения требований к его навигационному
2
обеспечению является этап посадки. Возможность полной автоматизации управления БПЛА на этом этапе полета во-многом зависит от характеристик ПНК.
Подавляющее большинство существующих систем обеспечения автоматической посадки БПЛА предполагает использование спутниковой навигационной системы (СНС) в комплексе с инерциальной навигационной системой [2]. Однако известно, что сигналы СНС подвержены помехам. Поэтому целесообразно исследовать возможности комплексирования инерциальной навигационной системы (ИНС) с системами, свободными от недостатков СНС, для анализа потенциальной точности определения навигационных параметров. В качестве такой системы в работе выбрана система определения параметров ориентации и навигации БПЛА на базе фотограмметрической обработки изображения искусственных наземных ориентиров [3, 4].
Система определения параметров ориентации и навигации БПЛА на базе фотограмметрической обработки изображения искусственных наземных ориентиров (фотограмметрическая система - ФС) в силу принципа своего действия обычно используется именно на этапе посадки. Для функционирования ФС используется изображение от камеры, которая, как правило, является частью полезной нагрузки БПЛА. Иллюстрация принципа действия ФС приведена на рисунке 1 (а - наклон линии глиссады, Г - плоскость глиссады, перпендикулярная плоскости точки касания Т).
навигации БПЛА на базе фотограмметрической обработки изображения
искусственных наземных ориентиров Камера БПЛА позволяет непрерывно получать изображение подстилающей поверхности, на которой расположены ориентиры. Координаты ориентиров известны в горизонтальной системе координат Оху2. На основе изображения определяется относительное пространственное положение камеры 5 (5х, 5у, 52) и ее угловое положение. Так как взаимное расположение системы координат камеры 5х'у'2' и связанной системы координат Схсус2с известно, определяются углы курса, тангажа и крена, а также координаты БПЛА.
Для функционирования ФС на этапе посадки камера должна иметь возможность фиксировать изображение трех наземных ориентиров Я, О и В (рисунок 1), координаты которых известны в системе координат, связанной с точкой
касания. Относительные координаты и угловое положение камеры могут быть вычислены с использованием метода, описанного в [3, 5, 6].
Для осуществления комплексной обработки навигационной информации, во-первых, необходимо описать математические модели погрешностей каждой из систем, подлежащих комплексированию, во-вторых, разработать математическое обеспечение комплексной обработки информации.
Модель погрешностей инерциальной навигационной системы
Для описания математической модели погрешностей инерциальной навигационной системы была использована модель, предложенная в [7]. Эта модель имеет следующий вид:
х1
Х2
Ж хз
йг х4
х5
хб
0 0 0
0 0 0
0 0 0
^ + О1 -®02 ) (О 2 X )
+ О2 -®02
)
- (О7 + Ох • О2)
(О2 + Ох • о) (о2 + о
О-ОхО2) -О +0 0) О +¿2 + )
1 0 0
0 1 0
0 0 1
0 20 2 - 20^
- 20 2 0 2ОХ
2О7 - 2ОХ 0
Х^ 0 0 0 0
Х2 Х3 Х^ Х^ 0 0 0 0
0 0 0 а 0
+ 0 п -п У Пх • Р + ДПх
-п2 0 7 Апу
Х6 пу -ПХ 0 Апг
йг
а 0 О 2 -ог а АО
Р = -О 2 0 О X • Р + АО
7 ог -О X 0 8 АО
где х1гх2,х3 - погрешности БИНС в определении координат; х4,х5,х6 -производные компонент х1г х2, х3; а,р,у - угловые погрешности ориентации
х
9
г
измерительного трёхгранника БИНС относительного вычисленного, компоненты
5
вектора угла малого поворота 0 = [а,р,у]; bnx¡yy, - постоянные составляющие погрешностей акселерометров и гироскопов, приведённые к осям выбранной в качестве базовой системы координат; с0 - собственная частота колебаний ошибок
О 1
БИНС, частота Шулера (со0=1,25 10'c' ); QX 7 z, QX7 Z - проекции вектора
абсолютной угловой скорости вращения выбранного трёхгранника и его производные; ^ z - проекции вектора кажущегося ускорения центра масс объекта.
Полная погрешность БИНС в определении углов ориентации может быть представлена через углы а2, Р2, Y2 - ошибки ориентации вычисленного БИНС трехгранника относительно идеального географического:
«2 = а + а1; Р2 — Р - Р1' Y2 — Y - Yi,
«1 = Pi^^1' Yi — т1' tg^),
Pl Р2 К2
где а1 ,P1 ,Y1 погрешности ориентации вычисленного трехгранника относительно базового за счет наличия ошибок определения координат точки места; Ф - значение текущей широты; р1, р2 - радиусы кривизны, соответственно, меридионального и перпендикулярного ему сечений.
Связь приведенных ошибок ориентации БИНС и ошибок БИНС в определении углов крена, кура и тангажа может быть представлена следующим образом [7]:
^= ' «2 - О'' Р2,
бф = (sin (у) ■ cos(i|j) + cos (у) ■ sin(i|;) ■ tg(u)) ■ а2 + (sin (у) ■ sin(i|;) — cos (у) ■
С О S да ■ tg (и) ) ■ Р 2 - с О S (Y ■ Y 2,
6и = (с05(у) ■ С05(Ф) — 5т(у) 1 5ш(ф) 1 tg(u)) 1 (Х2 + (с05(у) 1 5ш(ф) + 5ш(у) 1 С08ф ■ tgu ■ (32—5ту -у2,
где - идеальные значения текущих углов истинного курса, тангажа и
крена; б ф, би , бу - погрешности определения углов ориентации алгоритма БИНС.
Приведение инструментальных ошибок гироскопов и акселерометров производится по следующим формулам:
Апх = С1 • Ащ + С12 • Ап2 + С13 • Ащ
Апу = С21 • Ащ + С22 • Ащ + С23 • Ап3 Ащ = С31 • Ащ + С32 • Ащ + С33 • Ап3
АО = Сп •АО + С12 •АО + С13 • АО3
АО„ = С21 •АО 1 + С22 •АО2 + С23 •АОз,
АО = С31 • АО + С32 • АО + С33 • АО3
где Су - элементы матрицы направляющих косинусов между связанной с БПЛА системой координат и выбранной вращающейся системой координат.
Модели погрешностей акселерометров и гироскопов при оценке потенциальной точности комплексной системы БИНС-ФС могут быть представлены в виде:
бП ; = ДП ; + бП 15 е, бп ; = Дщ + б П?° 1 5 6,
где , - шумовые составляющие погрешностей, описываемые
случайным гауссовским процессом типа "белого шума" с нулевым средним и известной интенсивностью; 1 = 1 , 2 , 3 - номер гироскопа или акселерометра.
Модель погрешностей фотограмметрической системы
Ниже приводится модель погрешностей системы определения параметров ориентации и навигации БПЛА на базе фотограмметрической обработки изображения искусственных наземных ориентиров для использования в алгоритмах комплексной обработки. Построение модели осуществлялось на основе результатов экспериментальных исследований погрешностей такой системы [5].
В ходе исследований была выявлена практически линейная зависимость погрешностей определения навигационных параметров от дальности до наземных ориентиров. Для иллюстрации этого факта на рисунке 2 для примера приведена погрешность определения одной из координат фотограмметрической системой в зависимости от времени приближения к точке касания взлетно-посадочной полосы. Коэффициенты линейной зависимости для разных навигационных параметров приведены в таблице 1. Диапазон дальностей, для которых применима приведенная математическая модель - от 0 до 1000 метров.
Ошибка определения у
200
150
3
сек
Рисунок 2 - Погрешность определения одной из координат фотограмметрической системой Погрешность определения навигационных параметров и параметров ориентации в ФС на постоянной дальности от наземных ориентиров были представлены в виде случайной составляющей с нормальным гауссовым распределением - белыми шумами с интенсивностями, приведенными в таблице 1.
Таблица 1 - Погрешности фотограмметрической системы
№ Погрешность Коэффициент Интенсивность Интенсивнос
п/п ухудшения шума на ть шума на
точности в максимальном минимально
зависимости удалении м удалении
от дальности
1 Горизонтальных координат 0,58 60 м 2 м
2 Высоты 0,435 45 м 1,5 м
3 Угла тангажа 0,014 1,5 град. 0,1 град.
4 Угла крена 0,025 3 град. 0,5 град.
5 Угла курса 0,0175 2 град. 0,25 град.
Комплексная обработка информации ИНС и ФС
Для комплексной обработки информации ИНС и ФС использована слабосвязанная схема комплексирования на основе оптимального фильтра Калмана [8]. В слабосвязанной схеме бесплатформенная ИНС (БИНС) и фотограмметрическая система (ФС) вырабатывают независимые решения. В связующем блоке фильтр Калмана на основании разности показаний БИНС и ФС формирует оценку вектора состояния, в который входят погрешности БИНС в определении навигационных параметров и инструментальные погрешности инерциального измерительного блока (ИИБ). Оценка вектора состояния
используется для коррекции БИНС и компенсации инструментальных погрешностей ИИБ (рисунок 3).
г
Рисунок 3 - Структурная схема слабосвязанной схемы комплексирования БИНС и ФС Оптимальный фильтр Калмана возможно применять при описании системы в пространстве состояний. Полная модель должна быть представлена уравнениями состояния и уравнениями измерений:
х = /<х + ви + а IV 2 = нх + V
где X - вектор состояния; и - вектор управления; ж - вектор шумов системы; V -
вектор шумов измерений; ^ - матрица динамики системы; в - матрица управления;
О - матрица шумов системы; Н - матрица измерений; 2 - вектор измерений.
11
Вектор состояния X при предложенной схеме комплексирования имеет вид:
Хт =[х1 х2 х3 х4 х5 х6 а Р у Д^ До2 Доз Д^ Дй2 Ди3].
Матрица динамики системы Р, соответствующая принятым вектору состояния и математической модели погрешностей БИНС, может быть записана следующим образом:
" 111 Р12 113
1 = 03x6 1 22 1 1 23
_06х6 0 6x3 06x6
где
Ри =
о о о
о2 + о2 -®о2
о о о
о 2 -о, о у
- (о 2 + о х о у) о 2 +о 2 -©
х у
о у -о х о 2
о о о
- (о у + о х о 2)
о х -о у о 2
- (о х +о у о 2) о 2 +о 2 + 2®
1 о о
о 1 о
о о 1
о 2о2 - 2о
- 2о 2 о 2ох
2оу - 2о х о
Рц =
0 0 0
0 0 0
0 0 0
0 - п п
2 у
П2 0 -п
- П, п 0
^13 =
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
с„ с с12 с с13 0 0 0
с с21 с с22 с с23 0 0 0
с с31 с с32 с с33 0 0 0
0 о2 -о у "0 0 0 с 11 с с12 с1
Р = 1 22 -о 2 0 о х 1 II 0 0 0 с с 21 с с22 с
у -о х 0 0 0 0 с с31 с с32 с3
'13 '23 '33
При этом вектор шумов системы ж :
Шт = {бП1П015е,бП2П015е,бПзП015е,бп1П015е,бп2П015е,бпзП015е}. Матрица О имеет вид:
2
о
2
о
х
х
о =
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
с 11 с с12 с с13 0 0 0
с с 21 с с22 с с23 0 0 0
с с13 с с32 с с33 0 0 0
0 0 0 с11 с с12 с с13
0 0 0 с с21 с с22 с с23
0 0 0 с с13 с с32 с с33
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
Вектор измерений 2 формируется как разность показаний БИНС и ФС:
2 ~ [х1 - дХФС Х2 - дуфС Ф ИНС - Ф ФС ЯИНС - дЯфС д§ИНС - д8фС ]
где дхФС, дуФС - ошибки в определении горизонтальных координат; ду ФС, финс - ошибки в определении угла курса ФС и ИНС, соответственно; диФС, МИНС -ошибка в определении угла тангажа ФС и ИНС, соответственно; д%ФС, д%шс -ошибка в определении угла крена ФС и ИНС, соответственно.
При выбранном векторе измерений вектор шумов измерений V примет вид:
vт = [^хфс дуфс ду фс д фс дёфс ].
Таким образом, приведенные выше математические модели погрешностей ИНС и ФС в совокупности со схемой и математическими соотношениями для осуществления комплексирования являются основой для построения комплексной навигационной системы ИНС-ФС.
Исследование характеристик комплексной системы
Исследование характеристик комплексной системы осуществлялось методом имитационного моделирования (структура программного комплекса имитационного моделирования приведена на рисунке 4). Точность комплексной системы, работающей на основе предложенных технических решений, определяется в первую очередь точностью оценивания вектора состояния в ОФК. Таким образом, сравнив оценку вектора состояния с точными значениями искомых ошибок, можно судить о точностных характеристиках комплексной системы. Кроме того, для того чтобы сделать вывод о работоспособности предложенной схемы комплексирования необходимо оценивать качество самого процесса фильтрации. Исходя из этого, при разработке программного комплекса, моделированию были подвергнуты следующие процессы:
- моделирование траектории движения БПЛА - для обеспечения информацией о проекциях вектора кажущегося ускорения п, вектора абсолютной угловой скорости и, линейной скорости и, высоты полета к, истинного курса ^, дальности до наземных ориентиров и другим параметрам;
- моделирование ошибок БИНС и фотограмметрической системы (в том числе ошибки вычисления параметров ориентации);
- построение измерений как разности ошибок БИНС и ФС;
- оптимальное оценивание вектора состояния комплекса;
- сравнение полученных оценок ОФК и реальных ошибок подсистем;
- оценивание качества фильтрации.
Рисунок 4 - Структура программного комплекса имитационного
моделирования
Моделирование проводилось с использованием траектории захода БПЛА на посадку по глиссаде. Максимальное удаление от ВПП - 1000 м. Ниже приведены результаты исследования, проведенного по двум сценариям. Сценарий 1 предусматривает использование гироскопов (датчиков угловой скорости - ДУС), точность которых соответствует современным микроэлектромеханическим системам (постоянный и случайный дрейфы - 1 град./ час), а сценарий 2 -волоконно-оптических гироскопов (постоянный и случайный дрейфы - 0,1 град./ час). Оба сценария предполагают наличие в составе БИНС акселерометров с постоянной и случайной погрешностью 0,01 м/с . Выбор такого состава датчиков продиктован спецификой задачи - навигация беспилотных ЛА малого и среднего класса.
Сценарий 1
На рисунке 5 показаны результаты моделирования ошибки оценки определения координат по широте и долготе и среднеквадратическое отклонение (СКО) ошибок оценивания. Из графиков видно, что основное оценивание заканчивается на 10 секунде. В установившемся режиме ошибка оценки широты находится на уровне 3 метров, СКО ошибки сходится к значению 7 метров. Ошибка оценки по долготе находится на уровне 10 метров, СКО ошибки сходится к 8 метрам на 100 секунде.
Рисунок 5 - Ошибки оценки определения координат по широте и долготе и СКО
ошибок оценивания На рисунке 6 приведена ошибка оценки определения северной и восточной составляющей путевой скорости и СКО ошибок оценивания. СКО ошибки оценивания для северной составляющей путевой скорости после 10 секунды значительно не меняется, для восточной составляющей путевой скорости - после 40 секунды. В установившемся режиме ошибка оценки северной составляющей путевой скорости находится на уровне 0,05 м/с, СКО ошибки сходится к значению
0,1 м/с. Ошибка оценки восточной составляющей путевой скорости находится на уровне 0,2 м/с, СКО ошибки сходится к 0,2 метрам на 100 секунде.
Рисунок 6 - Ошибки оценки определения северной и восточной составляющей путевой скорости и СКО ошибок оценивания На рисунке 7 показаны ошибки построения базового трехгранника по углам а и у и СКО ошибок оценивания этих углов. СКО ошибок оценивания углов с практической точки зрения не меняется существенно после 5 секунд. В установившемся режиме ошибка оценки углов а, у находится на уровне 0,1 градуса,
СКО ошибки также сходится к значению 0,1 градуса.
Время (сек.) Время (сек.)
Рисунок 7 - Ошибки построения базового трехгранника по углам а, у и СКО ошибок
оценивания этих углов На рисунке 8 приводятся результаты моделирования процесса оценки погрешностей гироскопов и СКО ошибок оценивания. Анализ результатов позволяет сделать вывод, что процесс оценивания затянут во времени. На заключительном этапе моделирования ошибка оценивания дрейфа уменьшается с 1 град/час до 0,2 град/час для гироскопа №1. СКО ошибки оценки дрейфа уменьшается с 1 град/час до 0,9 град/час. В условиях проведенного моделирования дрейф азимутального гироскопа №2 не оценивается.
Рисунок 8 - Ошибки оценки и СКО ошибок оценок постоянных составляющих
дрейфов ДУС
На рисунке 9 изображены результаты моделирования процесса оценки постоянных составляющих ошибок акселерометров. Процесс оценивания также, как
и в случае с гироскопами, затянут во времени. Оценивание постоянных
2 2
составляющих ошибок акселерометров уменьшается с 0,01 м/с до 0,004 м/с для
9 9
акселерометра №1. СКО ошибки оценки уменьшается с 0,01 м/с до 0,007 м/с . В условиях проведенного моделирования постоянная составляющая ошибки
акселерометра №2 не оценивается.
Рисунок 9 - Ошибки оценки и СКО ошибок оценок постоянных составляющих
ошибок акселерометров
Сценарий 2
На рисунке 10 показаны ошибки оценки определения координат по широте и долготе и СКО ошибок оценивания. После 25 секунды СКО ошибок оценивания значительно не меняется. В установившемся режиме ошибка оценки широты находится на уровне 5 метров, СКО ошибки сходится к значению 7 метров. Ошибка оценки по долготе находится на уровне 15 метров, СКО ошибки сходится к 8 метрам на 100 секунде.
Рисунок 10 - Ошибки оценки определения координат по широте и долготе и СКО
ошибок оценивания На рисунке 11 представлены ошибки оценки определения северной и восточной составляющей путевой скорости и СКО ошибок оценивания. С практической точки зрения СКО ошибок оценивания уже после первых трех секунд существенно не меняется. В установившемся режиме ошибка оценки находится на уровне 0,1 м/с, СКО ошибки также сходится к значению 0,1 м/с.
Рисунок 11 - Ошибки оценки определения северной и восточной составляющей путевой скорости и СКО ошибок оценивания На рисунке 12 показаны ошибки построения базового трехгранника по углам а, у и СКО ошибок оценивания этих углов. СКО ошибок оценивания углов сходится к установившемуся значению после 5 секунд. В установившемся режиме ошибка оценки углов а, у находится на уровне 0,1 градуса, СКО ошибки также сходится к значению 0,1 градуса.
Рисунок 12 - Ошибки построения базового трехгранника по углам а, у и СКО
ошибок оценивания этих углов На рисунке 13 показаны ошибки оценки показаний гироскопов и СКО ошибок оценивания. В условиях проведенного моделирования дрейф гироскопов №1 и №2 не оценивается.
Рисунок 13 - Ошибки оценки и СКО ошибок оценок постоянных составляющих дрейфов ДУС На рисунке 14 изображены ошибки оценки постоянных составляющих ошибок акселерометров. Процесс оценивания затянут во времени. Для акселерометра №1
ошибка оценки постоянных составляющих погрешностей акселерометров
2 2 2 уменьшается с 0,01 м/с до 0,004 м/с , а СКО ошибки оценки уменьшается с 0,01 м/с
до 0,007 м/с . В условиях проведенного моделирования постоянная составляющая
ошибки акселерометра №2 не оценивается.
Рисунок 14 - Ошибки оценки показаний акселерометров и СКО ошибок
оценивания
Заключение
В работе предложена схема и основные математические соотношения для комплексирования бесплатформенной инерциальной навигационной системы и системы определения параметров ориентации и навигации БПЛА на базе фотограмметрической обработки изображения искусственных наземных ориентиров. Использование в качестве корректора БИНС фотограмметрической системы повышает информационную избыточность пилотажно-навигационного комплекса БПЛА при возможных потерях сигналов спутниковых навигационных систем, что особенно важно на этапе посадки. Кроме того, фотограмметрическая система позволяет получать информацию о параметрах ориентации БПЛА, что значительно улучшает характеристики ПНК в определении углов курса, тангажа и крена.
Приведенные результаты имитационного моделирования демонстрируют работоспособность предложенного подхода к построению системы. Точность определения параметров ориентации находится на теоретически ожидаемом уровне. Важным результатом исследования можно назвать выявленный факт незначительного влияния точности гироскопов (в выбранном диапазоне их точности) на точность оценки вектора состояния при работе системы в комплексном режиме.
Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 14-08-01080 а.
Библиографический список
1. Алёшин Б.С., Афонин А.А., Веремеенко К.К. и др. Ориентация и навигация подвижных объектов: современные информационные технологии. - М.: ФИЗМАТЛИТ, 2006, - 421 с.
2. Алешин Б.С., Антонов Д.А., Веремеенко К.К., Зимин Р.Ю., Жарков М.В., Кузнецов И.М., Пронькин А.Н. Малогабаритный интегрированный навигационно-посадочный комплекс // Труды МАИ, 2012, № 54: http://www.mai.ru/science/trudy/published.php?ID=29692
3. Лунев Е.М. Повышение точности определения навигационных параметров беспилотного летательного аппарата на базе фотограмметрических измерений на
этапе посадки // Вестник Московского авиационного института. 2011. Т. 18. № 2. С. 150-159.
4. Кузнецов А.Г. Автоматизация процесса посадки малогабаритного беспилотного летательного аппарата в особых ситуациях // Труды МАИ, 2011, № 45: http://www.mai.ru/science/trudy/published.php?ID=25425
5. Лунев Е.М. Исследование сходимости нового алгоритма определения навигационных параметров беспилотного летательного аппарата на базе фотоизображения // Труды МАИ, 2010, № 45 http: //www.mai.ru/science/trudy/published.php?ID=25431 &PAGEN_2=2
6. Лунев Е.М., Павлова Н.В. Программно-алгоритмическое обеспечение для определения навигационных параметров беспилотного летательного аппарата на базе фотоизображения // Вестник Московского авиационного института. 2009.Т. 16. № 6. С. 111-119.
7. Веремеенко К.К., Желтов С.Ю., Ким Н.В., и др. Современные информационные технологии в задачах навигации и наведения беспилотных маневренных летательных аппаратов. - М.: ФИЗМАТЛИТ, 2009. - 556 с.
8. Степанов О.А. Интегрированные инерциально-спутниковые системы навигации // Гироскопия и навигация. 2002. № 1 (36). С. 23-45