Труды МАИ. Выпуск №84
www.mai.ru/science/trudy/
УДК 527:519.8
О стендовой калибровке авиационных бескарданных инерциальных
навигационных систем
Вавилова Н.Б.*, Васинёва И.А.**, Парусников Н.А.
Московский государственный университет имени М.В. Ломоносова, Ленинские
горы д. 1, стр. 52, Москва, 119899, Россия *e-mail: nb-vavilova@yandex. ru **e-mail: [email protected]
Аннотация
Рассматривается задача калибровки бескарданной инерциальной навигационной системы (БИНС). Под калибровкой понимается определение параметров инструментальных погрешностей чувствительных элементов при помощи проведения экспериментов на специализированных стендах. Методика калибровки, позволяющая оценивать параметры инструментальных погрешностей инерциальных датчиков только при помощи показаний этих датчиков, переносится на калибровку на точных стендах. Показано, что дополнительная информация, доставляемая датчиками точного стенда, заметно улучшает точность калибровки.
Ключевые слова: бескарданная инерциальная навигационная система, калибровка, испытательный стенд, инструментальные погрешности, фильтр Калмана.
Введение
Основой приборного комплекса, решающего задачу навигации самолетов, служат в настоящее время инерциальные навигационные системы (ИНС). Получили распространение два типа таких систем: системы с горизонтируемой платформой и бескарданные инерциальные навигационные системы (БИНС). В настоящее время используются, как правило, системы второго типа.
Калибровка чувствительных элементов - необходимый этап технологического производства навигационных систем, предшествующий основным режимам эксплуатации ИНС - режимам начальной выставки и навигации. Калибровка состоит в определении параметров математической модели инструментальных погрешностей инерциальных датчиков с целью последующей компенсации этих погрешностей в режиме навигации. Математическая модель погрешностей задается априорно.
Традиционные способы калибровки, осуществляемые на специализированных стендах, состоят в сложной последовательности поворотов БИНС. При этом обычно навигационная система последовательно устанавливается в различных положениях. Такой эксперимент сложно организуется, потому что требуется обратная связь от стенда, обработка данных проходит в процессе эксперимента и необходимо изменение плана в зависимости от результатов промежуточных вычислений. Кроме того, в некоторых известных методиках калибровочными сигналами для датчиков
угловых скоростей (ДУС) служит угловая скорость Земли, что не обеспечивает хорошую обусловленность задачи оценки. Иногда также калибровка осуществляется по отдельности для ньютонометров (акселерометров) и ДУС, и тогда возникает задача согласования координатных трехгранников, связанных с разными типами датчиков.
В лаборатории управления и навигации МГУ имени М.В.Ломоносова разработаны методы калибровки, лишенные указанных недостатков. Применительно к грубым стендам они описаны в [1], [2]. Под условным термином "грубые" принимаются стенды, в которых не используется угловая информация платформы стенда относительно корпуса стенда. Ниже эти методы распространяются на методы калибровки на точных стендах, в которых такая информация используется. Но в этом случае приходится учитывать погрешности измерения углов и погрешности синхронизации информационных потоков во времени. При математическом описании задачи используются обозначения, системы координат и представления, изложенные в [3].
Математическое описание задачи
Прежде чем рассматривать конкретные варианты решения задачи калибровки, выведем из [3] общие соотношения решения такой задачи. Первичную инерциальную информацию доставляют три ньютонометра и три датчика угловой скорости (ДУС). Предполагается, что их оси чувствительности близки к осям трехгранника №\22 23, называемого приборным. М - местоположение
чувствительной массы приведенного трехосного ньютонометра. В проекциях на оси приборного трехгранника измеряется внешняя сила / 2 и его угловая скорость со2:
/2 = /2 + /
2>
= ^ ~Уг,
т
где А/2 = (А/2 , А/22, А/2 у - вектор погрешности измерений ньютонометров,
т
у 2 = (У21,у22У3) - вектор погрешности измерений ДУС.
Оси приборного трехгранника определяются следующим образом. Ось М21 выберем так, чтобы она совпадала с направлением оси чувствительности ньютонометра, который назван первым. Ось М22 выберем в плоскости, образованной осями чувствительности первого и второго ньютонометров, так чтобы ось М22 была ортогональна М2^. Ось М23 составляет с осями М2\, М22 правый ортогональный трехгранник. Угол между осью чувствительности второго ньютонометра и осью М22 и угол между осью чувствительности третьего ньютонометра и осью М23 предполагаются малыми. Собственные инструментальные погрешности каждого из ньютонометров включают в себя ошибку нулевого сигнала (ошибку нуля), погрешности масштабов, перекосы и высокочастотную составляющую, которая считается белым шумом. С учетом сказанного, вектор инструментальных
т
погрешностей А/2 = /'2 - /2 = (А/21, А/22, А/2з) имеет следующий вид
А/2 = А/2 +Г/2 +А/2, где Г =
'гп О о
Г21 г22 0 Г31 г32 г33
/ = (/ , А/2 , / )Т
21 22 23
высокочастотные погрешности типа белого шума. Для погрешностей ДУС
0 V
принимается аналогичная модель: = + ©щ + у2, где 0 =
©11 012 013 021 ©22 ©23 ч031 ©32 ©33 у
Введем исходную инерциальную систему отсчета О^^^ • О - геометрический
центр Земли. О^ - ось вращения Земли, направленная на северный полюс. О^^ -
плоскость Земного экватора. Ось О^ направлена на точку весеннего равноденствия.
Л2 - матрица ориентации приборного трехгранника относительно инерциального. Информацией об ориентации приборного трехгранника служит матрица Лу.
Ориентация трехгранника Му, называемого модельным, определяется
кинематическим соотношением: А у = со'х А у, где со =
г 0 (Щ -Щ 0 щ
КЩ2 -Щ1 0 у
Ориентацию приборного трехгранника относительно модельного определим
Т
вектором малого поворота (Зу = (/(,(2,(3) . Уравнение имеет вид:
0у = са'2/Зу + . Проще и удобнее уравнение, описывающее поведение вектора
( представлять в проекциях на оси опорного трехгранника, жестко связанного с географической вертикалью. Ось Мх1 касается проходящей через точку М параллели и направлена на восток, ось Мх2 лежит в меридиальной плоскости и направлена на север, ось Мх3 противоположна по направлению вектору силы тяжести. Вектор угловой скорости трехгранника Мх обозначим и х. Обозначим
через Ly матрицу ориентации модельного трехгранника, являющегося числовым образом трехгранника Mz, относительно Mx:
Ly = a>'zLy - Lyux , Ly (t0 ) = L0 , (*) тогда px = LyTpy,
¡3x=ux/3x + vx, vx = I<yVz ■ Информация, доставляемая ньютонометрами,
позволяет нам образовать вектор
Wx = fx - (0,0, g)T = Px (0,0, g)T + LyTAfz.
БИНС устанавливается на платформу стенда и может совершать некоторое программное угловое движение. Задача определения инструментальных погрешностей сводится к решению задачи оценивания, при этом используется алгоритм калмановской фильтрации. Возможность определения инструментальных параметров зависит от программных значений со z. Ив такой постановке задача калибровки решена в [2]. Для максимально возможной обусловленности решения задачи предложена процедура калибровки, включающая в себя три цикла. БИНС последовательно устанавливается на платформе стенда в трех различных положениях. Каждая из приборных осей последовательно совмещается с осью вращения стенда. Рассматривается только вариант, когда ось вращения стенда лежит в горизонтальной плоскости. Заранее ясно, что результаты калибровки при вращении вокруг вертикальной оси будут иметь существенно меньшую точность. Причина в том, что измеряемая в осях приборного трехгранника сила f с точностью до инструментальных погрешностей неподвижна относительно приборного трехгранника, в то время как в случае горизонтальной оси вращения эта сила
совершает относительно этого трехгранника значительные эволюции. Предполагается для определенности, что ось вращения близка к направлению оси Мх. Тогда на каждом из циклов с осью Мх совмещаются М^,Мг2,Мг3. Также предполагается, что установка осуществлена в соответствии с видом начальной матрицы Ьу, приведенным ниже.
Первый цикл.
Платформа таким образом вращается вокруг оси М7 приборного трехгранника, что матрица Ьу имеет вид:
г1 0 01 г1 0 0 1
Ьу (¿о) = 0 1 0 , Ьу - 0 соб ¥ Бт ¥
у0 0 1, v 0 - Бт ¥ соб ¥у
Здесь ¥(1;) - угол, на который поворачивается платформа стенда против часовой стрелки.
Второй цикл.
С осью Мх\ совпадает с точностью до погрешности установки ось М2.
'0 0 11 Г 0 - Бт ¥ соб ¥1
Ьу (¿0) - 1 0 0 , Ьу - 1 0 0
V0 1 0, V 0 соб ¥ Бт ¥у
Третий цикл.
С осью Мх\ совпадает с точностью до погрешности установки ось М3.
Г 0 1 01 Г 0 соб ¥ Бт ¥1
Ьу (¿0) - 0 0 1 , Ьу - 0 - Бт ¥ соб ¥
V1 0 0, V 1 0 0 ,
Угловая скорость выбирается в виде кусочно-постоянной функции, так что на каждом интервале постоянства можно написать ¥(t) = Qt. В соответствии с принятой программой вращения имеют место соотношения: на первом цикле сС = (Q,0,0)r, на втором цикле с/ = (0, Q,0)r, на третьем цикле с/ = (0,0, Q)r, где Q = const.
Заметим во избежание недоразумений, что приведенные выражения для I
используются здесь только для анализа возможностей калибровки. В реальности при калибровке эти матрицы доставляются решением уравнений Пуассона (*).
Для точных стендов существует возможность расширить вектор коррекции за счет дополнительной информации об углах ориентации платформы стенда относительно корпуса стенда. Выведем соответствующие соотношения, используя некоторые дополнительные трехгранники. Базовый трехгранник Mp, жестко связанный с корпусом стенда. Обычно стенд устанавливается на фундамент и ориентируется в географической координатной сетке, то есть в идеале совпадает с
трехгранником Mx0. Ориентацию трехгранника Mp относительно Mx0 определим
Т
постоянным вектором малого поворота Yp = (Yp, Yp , Yp ) . Для некоторых высокоточных стендов величина Yp оказывается настолько малой, что ей можно
пренебречь. Трехгранник Mz , жестко связанный с платформой стенда. При установке корпуса БИНС на платформе приборный трехгранник Mz, с точностью до
погрешности установки, совмещается с трехгранником Mz* . Ориентация
*
приборного трехгранника Mz относительно трехгранника Mz* определяется
8
Т
постоянным вектором малого поворота 32 = (£2 ,322,32) . Матрица ориентации
•г-
платформенного трехгранника стенда М2 относительного базового трехгранника
стенда Мр - Ь *. В свою очередь матрица Ь * может быть определена величинами
* * *
к*, к*, а~3 - углами последовательных трех поворотов, эквивалентных углам курса у, тангажа в и крена У в авиации. Имеют место соотношения Ь * = ¿3^^, где
Ь =
{ * * л
еовк* Бтк* О * *
- Бтк еоБЛ! О
О
О1
, Ь2
О О
О
о
л
О
* *
еоБк* Бтк2
* *
- б1ик2 еоБк* у
{ * п ■ *\
еоБк3 О - Бтк3
, Ь
О 1
О
О ообк
3 У
Выходной информацией стендовых угловых датчиков служат величины К*' = К* + р1, к* ' = к* + р2, к* ' = к* + Р3, где р1, р2, р3 - инструментальные
^^ ^^ ^^
погрешности в измерениях этих углов. Величины к* ,к* , позволяют построить
^^ ^^ ^^
матрицу Ь' * равную Ь(к* ,к* ,к* ) путем замены в матрице Ь * величин
1»
к* ( = 1,2,3) на их измеренные значения. Последняя матрица определяет ориентацию трехгранника М2* . Ориентацию трехгранника М2* определим
относительно трехгранника М2 вектором малого поворота ¡л2 = 2,^г3)т.
Положим АЬ* = Ь * - Ь *. Тогда ]й2 = АЬ*. Ь . или далее будут использоваться
проекции вектора Ц на оси трехгранника Мр или МхО. Поэтому ¡лх = Ь\ /и2 или
2
т
Т Л
¡их = Ь и2Ь * . Получим выражение для вектора¡лх через ошибки измерения углов
2 2
р1, р2, р3, его порождающие. Используя разложение Тейлора с оставлением только
линеиных слагаемых, получим
<тт <тт rj-i rj-i
¡ux = L* al** = L LL ( l3L2 aL + L ALL + AL3L2L ):
rp ^ T T / * T T T / *
L sl / 5a* p1 + l l2 sl2 / 5a2Lp2 + L l2 l3 sl3 / 5a3l2 L1p3 , = GÍP + G2P2 + ^зРз
где
Q =
^ 0 1 0Л -10 0 v 0 0 0,
0 0
0 0
Sin A - COS A
- Sin Aj
*
COSA* 0
, G3
0
*
sin A**
*
sin A** 0
* * * COS A* COSA2 sin A* COSA**
* *Л
COSA* COSA**
* *
-SinA* COSA** 0
V
После преобразования получим
* * *
Ju1 = - Sin A1 COS A2P3 + COS A1 P2, J2 = COS A* COS A*P3 + Sin A*p2, (**) j3 = Sin A*p3 + p1.
Далее при учете p, p, p может быть использовано обратное преобразование:
Pl = J3 - SinA2P3,
Р2 = Sin A* J2 + COS A* J1, (***)
* * *
P3 COS A* = COS A1J2 - Sin A1 J.
Кроме того, очень важно учитывать то обстоятельство, что данные,
доставляемые стендом, и данные БИНС не являются синхронизированными, каждые
имеют свою шкалу времени и разную дискретность. Обозначим запаздывание данных стенда относительно данных БИНС через т.
Дополнительный вектор измерения, доставляемый при помощи информации стенда, образуем в двух вариантах: при помощи матрицы ориентации или при помощи углов.
Первый способ формирования дополнительных измерений - при помощи матриц ориентации.
Сначала рассмотрим выражение - разность матриц, вычисленных при помощи информации БИНС и при помощи углов от стенда, без учета ошибок синхронизации.
ЬУ - Ь2* = Ь2 + 02Ь2 - Ь * - &2Ь * • ' 2 2
Распишем входящую в это выражение разность:
Ь2 - Ь * = Ь - Ь™ + Ь™ - Ь * = -Ь у п + 8Ь * •
2 Ь2р + Ь2р Ь * = Ь2Ур +82Ь *•
Здесь Ь - матрица ориентации приборного трехгранника БИНС М2 относительно базового трехгранника стенда Мр . Таким образом,
ЬУ - Ь2* = @2Ь2 - &2Ь * - Ь2Ур + 82Ь * •
г г
С точностью до членов второго порядка малости можно представить:
Ьу - Ь2* = @2Ьу - ^2Ьу - ЬуУр + 82Ьу •
Образуем векторы измерений в проекции на оси трехгранника М и в проекции на
*
оси трехгранника Мр, эти векторы обозначим через Ж; и Жр. Положим с учетом запаздывания
/v * т1
Wz - (ь у (I) - ь; *« -т))\1у ц).
Оставляя первый порядок малости по г и используя для выражения производной матрицы Ь уравнение Пуассона (*), получим:
* - К - Ь; * ур + 4 - +
т т
Далее имеем Жр - Ь; * * или Жр - Ь;* * Ь;*. Из последних соотношений получим: Жр - -К -ур + ЬТ; * 4 Ь; * - МХ +
Итак, векторы измерений имеют вид:
Ж; -К; - Ь ;*Ур + £; - + 0;Г, Жр -Кх -Ур + Ь;* -Лх + 0 хг. Здесь - (0,0,0)г на всех трех циклах.
Эквивалентная форма построения вектора измерений. Матрицу Ь; можно
* * *
рассматривать как функцию трех углов к*,к*,к*. Модельные значения этой
111 Т Т Т Т I
матрицы Ьу (к ,к* ,кз ). Вектор измерения Жр - (Ж[ ,Жз ) образуем следующим образом:
*г
Ж1 = к[{О - А (Г - г) = А^ - А +
= 4(0 - *2 - т) = ~ />2 +
Щ* = СОБ к'2 (0(^*3 (0 ~ КЪ (/ ~ Г)) = СОБ ЛГ2 (7)(ДЛГ3 - /?3 ) + СОБ АГ^ЛГ3 (7)"Г,
12
где ÁK¿ = к ' -к*, i -1,2,3..
Здесь множитель cosk2 вводится, как это следует из дальнейшего для того, чтобы избежать вырождения при возможных эволюциях платформы стенда, когда угол к'2 оказывается близким к 90°.
Для того, чтобы выразить ák1, ák2, ák3 через компоненты малых углов Рх УР А можно воспользоваться выводом уравнений (***) Аналогично тому, как
было определено преобразование между ошибками измерений углов и
компонентами вектора малого поворота, получим:
*
w1 =-Рз-Yp3-$x3 -Р1 -qx3г,
W2* = - sin к'фг - cos к'фх + sin к[ур2 + cos к[ур1 - sin к[8х2 - cos k[Sx1
-Р2 - sin K[Q.x2Г - cosk1Q
W3* = - cosk{^2 + sinK^ + cosk"1^p2 - sinKVp1 - COS 2 + sinK1Ax1 - cosk2p3 - cosk1Qx2r + sinK"1Qx1r, sx = LTySz.
Далее предполагаем, что используются высокоточные стенды, типа [4], поэтому параметрами у, р можно пренебречь.
Итак, задача калибровки сводится к оцениванию постоянных величин
V0 ,&ij,Áf0 ,Tij а ,г, i, j = 1,2,3, при помощи измерений Wz и дополнительных zi J zi J> zi
измерений, доставляемых точным стендом, например, в форме Wp.
Анализ наблюдаемости
При выборе алгоритма калибровки и анализе точности основной инструмент исследования - ковариационный анализ. Тем более важно, что он позволяет учесть
уровень шумов. Строго говоря, ковариационного анализа оказывается для этого и достаточно. Но, следуя традициям, и, может быть, для лучшего понимания задачи предварительно проведем анализ наблюдаемости. При этом мы, конечно, понимаем ограниченность такого анализа для решения задачи.
В начале проведем такоИ анализ применительно к измерениям Wx, в которых не используются измерения, связанные со стендом. Поскольку калибровочным сигналом является угловая скорость вращения стенда, которая на несколько порядков больше угловоИ скорости вращения Земли, то величиной ux можно пренебречь. Поэтому имеем:
Px=LyT{vQz+®co'z).
При анализе наблюдаемости будем использовать то обстоятельство, что если наблюдается какая-то переменная, то наблюдается ее производная. Это позволяет исключить величины р из вектора коррекции. Анализ наблюдаемости сводится таким образом к анализу соотношений:
Wx = -gx /3x-LyTQ!z(Afz° + TLygx) + LyTmzLygx, где = (0,0,gf.
Первый цикл:
zj = - g cos Qt (v® +0*1^) + g Sin Qt (V3 +©3jQ),
z2 = g(v? + 011Q) - Q sin QtAf0 - Q COS QtAf0 - gQ cos2Qtr33 -
2
- gQ sin Qt (Г** +Г32) + gQ sin Qt cos Qt (Г ** -Г32) Z3 = -Q sin Qtf0 +Q cos QtAf2 -Qg cos2Qtr32 + gQ sin2Qt (Г 22 -Г33).
Второй цикл:
z1 = -g COS at(v° + 032a) + g Sin at(V +012Q-Qr21),
z2 = g (v0 + ©22a) - a cos at Áf10 - a sin atÁfO+ag cos2at (r33 - rn) -
-Qg sin2atT3!,
z3 = a cos atÁf30 - a sin atÁf +ag sin2at (Г33 - r^)+ag cos2atr31
Третий цикл:
zx = -g cos at(v0 + ©13a - ar31)+g sin atV+©23a - ar32),
z2 = g (v30 + ©33a) - a cos atÁf20 - a sin atÁf0 + +ag cos2at (rn - r 22) --ag sin2at r21,
z3 = a cos atÁf0 - a sin at f0+ag sin2at (rn -r22)+ag cos2atr 21.
Набор функций 1,sin at,cos Qt,sin2Qt,cos2Qt линейно независимый, поэтому наблюдаемыми являются постоянные множители перед этими функциями. Кроме того, в процессе калибровки величина a может менять знак. Легко убедиться в том, что наблюдаемы все параметры принятых моделей инструментальных погрешностей. При более тонком анализе следует учитывать порядок коэффициентов перед этими функциями, но как раз такой учет дает ковариационный анализ.
Анализ наблюдаемости с использованием дополнительной информации от стенда. На интервалах, где угловая скорость a постоянна, продифференцируем измерения. С учетом вышесказанного анализ наблюдаемости с использованием дополнительной информации проводится для соотношений:
Wp=Px-\Tz*co'Sz
На первом цикле:
Wp1 =vi0 +011Q,
Wpi = cos Qt(v2 +021Q) - sin Qt(v3 +031Q) -Qsin QtS2 -Q cos Qt43, Wp3 = sin Qt(v2 + 021Q) + cos Qt(v3 +031Q) + Q cos QtS2 -Q sin QtS3.
На втором цикле:
Wp1 =V20 + 022Q,
Wp2 = cos Qt (v3° + 032Q) - sin Qt (у0 + 012Q) - Q cos Qts1 - Q sin Qts3,
Wp = sin Qt (v30 + 032Q) + cos Qt (v0 + 012Q) - Q sin Qt^ + Q cos Qt43. На третьем цикле:
Wp1 =V30 + 033Q,
Wp2 = cos Qt(v? + 013Q) - sin Qt(vj + 023Q) - Q sin Qts1 - Q cos Qts2,
Wp3 = sin Qt (v0 + 013Q) + cos Qt (v2 + 023Q) + Q cos Qt^ - Q sin Qts2.
Аналогично рассуждениям, проведенным для анализа наблюдаемости без использования дополнительной информации от стенда, заметим, что параметры погрешностей стенда 8г= (4,42,4)T являются коэффициентами при независимых функциях. Следовательно, они наблюдаемы.
Здесь еще раз надо сказать, что вывод о приемлемости алгоритма доставляет ковариационный анализ.
Ковариационный анализ
При ковариационном анализе были выбраны следующие близкие к реальным априорные характеристики для параметров инструментальных погрешностей в виде их среднеквадратических значений. а 0 = 0.5 °/час, а^ = 0,02 м/с2 ,
аг/7 =а@ц = 10 , агу = ащ = 3'Предполагается, что между собой они не
коррелированы. При моделировании стандартные отклонения белых шумов в ДУС а 8 полагаются равными 0.1 °/час на частоте 1 Гц, а стандартные отклонения шумов
ньютонометров а 8 = 0,001 м/с2 на частоте 1 Гц. Погрешности информации от
стенда: ст§. =1°,стт = 0.05 с, а^ 5 =15" на частоте 1 Гц. Задаются программные
движения стенда (демонстрационный вариант) такие, что продолжительность каждого цикла 15 минут. Выбирается кусочно-постоянная угловая скорость. Период на одном цикле выбирается так, чтобы за половину длительности цикла происходило изменение угла с угловой скоростью 10°/с, далее изменение этого угла со скоростью -10°/с.
Результаты оценивания без использования информации от стенда
,°/час А/0, м/с2 гй Г ' г У,
а 0.01 5.1 -10_7 0.1 9-10_5 1.36-10_5 0.1
Результаты оценивания с использованием информации от стенда
у<0 ,°/час и У, А/0, м/с2 Гц Г ' г у,
а 0.003 9.7 -10_8 0.05 7.7-10_5 1.36-10_5 0.05
Для того чтобы провести сравнение между собой результатов калибровки с использованием измерений от точного стенда и результатов без использования таковых, проведено моделирование дисперсионного уравнения ошибок автономной навигации на траектории "змейка". В начальные условия для ковариационной матрицы ошибок БИНС в блоке инструментальных погрешностей подставлялись значения, полученные в результате калибровки, а ошибки выставки также определялись результатами калибровки. Моделировался полет самолета по данной траектории в течение 1 часа. Критерием качества калибровки, также как и в [5],
навигационного эллипсоида, Л^, АЛ - ошибки в определении широты и долготы. После калибровки без использования измерений точного стенда ошибка автономной навигации составила 1890 м, при использовании новых измерений - 1100 м. Таким образом, использование новых измерений позволяет без усложнения плана калибровки и увеличения времени повысить качество калибровки БИНС и улучшить точность автономной навигации в 1.5 раза.
1. Результаты ковариационного анализа предложенных методов калибровки без учета информации от стенда и с ее использованием показывают их высокую эффективность, поэтому эти методы могут быть рекомендованы к реализации.
была выбрана величина
, где а - длина большой полуоси
Выводы:
2. Дополнительная информация, доставляемая стендом, несколько улучшает результат калибровки, и, по возможности, ее необходимо использовать.
Библиографический список
1. Парусников Н.А. Задача калибровки бескарданной инерциальной навигационной системы на стенде // Известия РАН. Механика твердого тела. 2009. № 4. С.3-6.
2. Голован А.А., Парусников Н.А. Математические основы навигационных систем. Часть II. - М.: МАКС Пресс, 2012. - 170 с.
3. Голован А.А., Парусников Н.А. Математические основы навигационных систем. Часть I. - М.: МАКС Пресс, 2011. - 132 с.
4. URL: http: //www.acutronic.com/
5. Васинёва И.А., Кальченко А.О. Анализ точности калибровки бескарданной инерциальной навигационной системы в полете в зависимости от некоторых типов эволюций самолета // Вестник Московского университета. Математика. Механика. 2014. № 1. С. 65-68.