Разработка алгоритма калибровки инерциальной навигационной системы на двухосном испытательном стенде

advertisement
Электронный журнал «Труды МАИ». Выпуск № 63
www.mai.ru/science/trudy/
УДК 681.2.089:681.2.087
Разработка алгоритма калибровки инерциальной навигационной
системы на двухосном испытательном стенде
К.К. Веремеенко, И.А. Галай
Аннотация
В статье рассматриваются методика и алгоритм калибровки бесплатформенной инерциальной навигационной системы (БИНС) авиационного класса точности по координатной, скоростной и угловой информации. В качестве входных данных для алгоритма калибровки используются сигналы с выхода системы. Приведенные методика и алгоритм калибровки
БИНС реализуются на двухосном поворотном столе. Приведены результаты имитационного
моделирования процесса калибровки с использованием разработанного алгоритма.
Ключевые слова: калибровка; поворотный стол; гироскоп; акселерометр; математическая
модель; алгоритм; бесплатформенная инерциальная навигационная система.
Одна из ключевых задач при создании бесплатформенных инерциальных навигационных систем (БИНС) - обеспечение их требуемой точности. Решить ее можно только обеспечив необходимый уровень инструментальных ошибок гироскопов и акселерометров. Один из
способов решения этой задачи – это калибровка датчиков, которая заключается в нахождении коэффициентов модели ошибок ее датчиков поступательного и углового движения или
калибровочных коэффициентов. Существует два подхода к проведению испытаний. Различаются они по типу информации, используемой в алгоритмах идентификации калибровочных коэффициентов. В первом случае информация напрямую снимается с блока чувствительных элементов БИНС, во втором подходе используется выходная информация БИНС о
широте, долготе и высоте местоположения, составляющих путевой скорости и скорости изменения высоты, а также информация об углах ориентации подвижного объекта. Реализация
первого метода может быть затруднена, если блок уже установлен в инерциальную навигационную систему и его выходы недоступны. Такие модули калибруются отдельно от инер1
циальной навигационной системы. Учитывая это, первый способ пригоден только для систем
собственной разработки, где имеются выходы инерциальных датчиков, или отдельно закупаемых трехмерных инерциальных модулей. Второй подход заключается в снятии сигналов
непосредственно с выходов системы. Располагая данными об ошибке определения навигационных параметров и параметров ориентации, задавшись моделью ошибок БИНС и моделью
ошибок гироскопов и акселерометров, рассчитываются калибровочные коэффициенты для
датчиков первичной информации, входящих в состав системы. Идея использования математической модели ошибок БИНС при калибровке заключается в нахождении соотношений, с
помощью которых по измеренным выходным сигналам системы можно определить калибровочные коэффициенты инерциальных измерителей. Реализация такого способа стала возможной с появлением прецизионных поворотных столов, которые определяют ориентацию
БИНС относительно точно выставленной в горизонте и в азимуте планшайбы динамического
симулятора. В статье предлагается вариант алгоритма калибровки инерциальной навигационной системы, в котором производится обработка информации с выхода системы и для
идентификации параметров математических моделей ошибок инерциальных измерителей
применяется оптимальный фильтр Калмана. Структурная схема процесса калибровки отражена на рисунке 1.
  X 
X

X , KA, KГ

 A, K
Г
X,K
Модель ошибок
Ф, R, Г, Q
KГ , KA
ОФК
X
Модель ошибок
гироскопов и
акселерометров
, K Г
 n, K A
Z
 ,  , h, Vn , VE
Уравнение ошибок
БИНС
 ,  , 
V
Параметры сценария
моделирования
Формирование шумов
корректора
НУ
интегрирования
Рисунок 1. Схема калибровки БИНС.
Предлагаемый алгоритм использует выходные данные о местоположении, скорости и
2
углах, которые поступают с выхода БИНС. В качестве алгоритма оптимальной обработки
информации применяется оптимальный фильтр Калмана (ОФК). Чтобы иметь возможность
применить ОФК, уравнения системы должны быть записаны в пространстве состояний. Полная модель содержит уравнение состояния системы и уравнение измерений в следующей
форме:
𝑥 = 𝐹𝑥 + 𝐺𝑤 ,
𝑧 = 𝐻𝑥 + 𝑣
(1)
где 𝑥- вектор состояния;𝑤 - вектор шумов системы;𝑣 - вектор шумов измерений;𝐹 - матрица
динамики системы;𝐺 - матрица шумов системы; 𝐻 - матрица измерений; 𝑧 - разность показаний ИНС и эталонного оборудования.
Модель ошибок БИНС в скалярно матричной форме записи имеет вид (2) - (3) [1].
Модель позволяет исследовать связь между инструментальными погрешностями первичных
измерителей, неточностями задания начальных условий и ошибками БИНС в определении
местоположения, скорости и углов ориентации объекта.
0
𝑥1
2
2
2
Ω
+
Ω
𝑥2
𝑦
𝑧 − 𝜔0
𝑑 𝑥3
0
=
− Ω𝑧 − Ω𝑥 Ω𝑦
𝑑𝑡 𝑥4
𝑥5
0
𝑥6
Ω𝑦 − Ω𝑥 Ω𝑧
𝑥1
0
𝑥2
0
𝑥3
0
∙ 𝑥 +
−𝛾
4
𝑥5
0
𝑥6
𝛽
1
0
0
−2Ω𝑧
0
2Ω𝑦
0
𝛾
0
0
0
𝛼
0
−𝛽
0
𝛼
0
0
0
𝛼
𝑑
𝛽 = −Ω𝑧
𝑑𝑡 𝛾
Ω𝑦
0
Ω 𝑧 − Ω𝑥 Ω𝑦
0
2
Ω𝑥 + Ω2𝑧 − 𝜔02
0
− Ω𝑥 − Ω𝑦 Ω𝑧
0
0
0
0
0
0
0
0
0
0
0
0
Ω𝑧
0
−Ω𝑥
0
2Ω𝑧
1
0
0
2Ω𝑥
0
− Ω𝑦 − Ω𝑥 Ω𝑧
0
Ω𝑥 − Ω𝑦 Ω𝑧
0
2
2𝜔0 + Ω2𝑥 + Ω2𝑦
0
0 𝑛𝑥
∆𝑛
0 𝑛𝑦
𝑥
0
0 𝑛𝑧
∙
+ ∆𝑛 ;
0
0
𝑦
0
0
0
0
0
∆𝑛𝑧
−Ω𝑦 𝛼
∆Ω𝑥
Ω𝑥 ∙ 𝛽 + ∆Ω𝑦 ,
𝛾
0
∆Ω𝑧
где𝑥1 , 𝑥3 ,𝑥5 – ошибки определения координат и высоты ИНС, м;
𝑥2 , 𝑥4 , 𝑥6 - производные компонент 𝑥1 , 𝑥3 , 𝑥5 , м/с;
3
0
−2Ω𝑦
0
2Ω𝑥
1
0
(2)
(3)
 ,  ,  – угловые погрешности ориентации измерительного трѐхгранника относительного
  [  ,  , ]
вычисленного, компоненты вектора угла малого поворота 
, рад;
 x ,  y ,  z – погрешности гироскопов, приведѐнные к осям вращающейся системы координат, в которой ведѐтся решение уравнений ошибок, рад/с;
∆𝑛𝑥 , ∆𝑛𝑦 , ∆𝑛𝑧 – погрешности акселерометров, приведѐнные к осям вращающейся системы ко-
ординат, в которой ведѐтся решение уравнений ошибок, м/с2;
n x , n y , n z – проекции вектора кажущегося ускорения центра масс объекта, м/с2;

 X ,Y , Z , 
X ,Y ,Z - проекции вектора абсолютной угловой скорости вращения выбранного трѐх-
гранника и его производные, рад/с2;
 0 – частота Шулера (ωo=1,24 10-3c-1 ).
Для обеспечения процесса калибровки в приведенных соотношениях проекции погрешностей гироскопов и акселерометров  x ,  y ,  z и ∆𝑛𝑥 , ∆𝑛𝑦 , ∆𝑛𝑧 необходимо выразить непосредственно через их инструментальные ошибки, которые приводятся к необходимой вращающейся системе координат согласно формулам:
x  C111  C12 2  C133 ;
 y  C211  C22 2  C233 ;
z  C311  C32 2  C333 ;
nx  C11n1  C12 n2  C13n3 ;
(4)
ny  C21n1  C22 n2  C23n3 ;
nz  C31n1  C32 n2  C33n3 ;
здесь C11..C33 - элементы матрицы направляющих косинусов между связанной с объектом
системы координат и выбранной вращающейся системой координат.
Инструментальные ошибки датчиков представляются своими математическими моделями. В данной статье для иллюстрации взят один из наиболее часто встречающихся вариантов представления таких моделей [1]:
4
1  m0 x  m1x   ny  nz   m3 x  x  1;
(5)
2  m0 y  m1 y   nx  nz   m3 y  y  2 ;
3  m0 z  m1z   nx  ny   m3 z  z  3 ;
n1  a1x  a2 x nx   n1;
(6)
n2  a1 y  a2 y ny   n2 ;
n3  a1z  a2 z nz   n3 ;
где m0 x , m0 y , m0 z - постоянный дрейф соответствующего гироскопа; m1x , m1y , m1z - дрейф
гироскопа, зависящий от первой степени перегрузки (в разных типах гироскопа причины такой зависимости могут быть разными); m3x , m3 y , m3z - дрейф гироскопа, вызванный нестабильностью масштабного коэффициента; a1x , a1y , a1z - постоянная составляющая ошибки
акселерометра; a2 x , a2 y , a2 z - параметры, характеризующие ошибку акселерометра из-за неточности масштабного коэффициента; nx , ny , nz - проекции вектора кажущегося ускорения
на оси соответствующих гироскопов;  x ,  y ,  z - проекции абсолютной угловой скорости
вращения на оси чувствительности гироскопов.
Вектор состояния системы, объединяющий ошибки БИНС и коэффициенты модели
ошибок датчиков, имеет размерность 24:
X   x1 , x2 , x3 , x4 , x5 , x6 ,  ,  ,  , a1x , a1 y , a1z , a2 x , a2 y , a2 z , m0 x , m0 y , m0 z , m1x , m1 y , m1z , m3 x , m3 y , m3 z 
Т
(7)
Вектор шумов системы W с учетом принятой модели (5), (6) имеет вид:
W   n1  n2  n3 1 2 3 
T
Матрица интенсивности шумов системы 𝑄 имеет вид:
5
(8)
2
𝜎11
0
0
𝑄=
0
0
0
0
2
𝜎22
0
0
0
0
0
0
2
𝜎33
0
0
0
0
0
0
2
𝜎44
0
0
0
0
0
0
2
𝜎55
0
0
0
0
∙ 𝑇,
0
0
2
𝜎66
(9)
2
2
2
2
2
где  112 ,  22
,  33
- СКО случайной ошибки акселерометров;  44
,  55
,  66
- СКО случайной
ошибки гироскопов; Т–шаг дискретизации алгоритмов оценивания.
Матрица динамики системы Fв рассмотренном случае будет иметь вид:
F
F   11
0159
F13 
.
0159 
F12
0156
(10)
Блочные элементы этой матрицы определяются соотношениями:
0
1

 2
2
2
0
  y   z  0 

0
0

  
     2
z
x
y
z

F11  
0
0

   
2 y
 
y
x
z

0
0


0
0

0
0

0
 
z
  x y 
0
0
0
0
0
2 z
   
 
y
x
z
2 y
0
 nz
1
0
0
0
0
2 x
nz
0
1
0
0
0
ny
nx
0
0
0
0
 z
y
z
0
 x
0
(     )
0
0
0
2
x
2
z
2
0
 
0
 0
C
 11 C12
 0
0

C21 C22
F12   0
0

C31 C32
 0
0

0
 0
 0
0

 2
2
0
0
0
0
0
C13
0
C23
0
C33
0
0
0
  y z 
0
     2
 
x
y
z
x
0
0
0
x
  2x   2y 
0
0
0
0
C11nx
0
C21nx
0
C31nx
0
0
0
6
0
C12 n y
0
C22 n y
0
C32 n y
0
0
0
0 
C13nz 
0 

C23nz 
0 ;

C33nz 
0 

0 
0 
0 

ny 
0 
 nx 

0 ;

0 

 y 
x 

0 
(12)
(11)
0
 0
 0
0

 0
0

0
 0
 0
0
F13  
0
 0
C
 11 C12

C21 C22

C31 C32
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
C13
C11   n y  nz  C12   nx  nz  C13   nx  n y  C11 x
C12 y
C23
C21   n y  nz  C22   nx  nz  C23   nx  n y  C21 x
C22 y
C33
C31   n y  nz  C32   nx  nz  C33   nx  n y  C31 x
C32  y







;

C13 z 

C23 z 

C33 z 
0
0
0
0
0
0
(13)
Матрица шумов системы G может быть записана через блочные матрицы:
G 
G   11  ,
0156 
(14)
где матрица G11 представима в виде:
0
 0
C
 11 C12
 0
0

C21 C22
G11   0
0

C31 C32
 0
0

0
 0
 0
0

0
C13
0
C23
0
C33
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
C11 C12
C21 C22
C31 C32
0 
0 
0 

0 
0 ;

0 
C13 

C23 
C33 
(15)
Матрица связи измерений и вектора состояния 𝐻имеет вид:
𝐻 = 𝐻11
09∗15 ,
где матрица 𝐻11 равна:
7
(16)
1

0
0

( R  h)

1

0
0

( R  h) cos 

0
0
0

 
 H   tan  1
 z
x
 ( R  h)

H
H11  
0
0


( R  h)

0
0
0


1

0
0

( R  h)


1

0
0
( R  h)


1

 tan 
0
0
 ( R  h)
0
0
0
0
0
1
0
y
1  x
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 ;


1 0 0 0

0 1 0 0


0 0 1 0


0 0 0 1

(17)
здесь𝜑 - широта местоположения поворотного стола с установленной на нем БИНС.
Матрицу шумов измерений запишем в следующей форме:
𝑉 = 𝛿𝜑эт
𝛿𝜆эт
𝛿ℎэт
𝛿𝑉𝑛_эт
𝛿𝑉𝑒_эт
𝛿ℎэт
𝛿𝛼эт
𝛿𝛽эт
𝑇
𝛿𝛾эт ,
(18)
где 𝛿𝜑эт, 𝛿𝜆эт, 𝛿ℎэт - ошибки определения эталонного значения соответственно широты, долготы и высоты места проведения эксперимента, получаемые в результате измерения местоположения стенда с предустановленным прибором высокопрецизионным устройством;𝛿𝑉𝑛_эт
- ошибка определения Северной составляющей путевой скорости измерителями поворотного
стенда; 𝛿𝑉𝑒_эт - ошибка определения Восточной составляющей путевой скорости измерителями поворотного стенда; 𝛿ℎэт - ошибка определения изменения высоты измерителями стенда; 𝛿𝛼эт,𝛿𝛽эт,𝛿𝛾эт- ошибки поворотного стенда в определении ориентации.
Матрицу интенсивности шумов измерений представим в виде:
8
𝑄=
2
𝜎𝛿𝜑
0
0
0
0
0
0
0
0
0
0
0
2
𝜎𝛿𝜆
0
0
0
2
𝜎𝛿ℎ
0
0
0
2
𝜎𝛿𝑉
𝑛
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
2
𝜎𝛿𝑉
𝑒
0
0
0
0 ∙ 𝑇,
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
2
𝜎𝛿ℎ
0
0
0
𝜎𝛼2
0
0
0
𝜎𝛽2
0
0
0
0
0
0
0
0
0
0
0
𝜎𝛾2
(19)
2
2
2
2
2
2
где 𝜎𝛿𝜑
, 𝜎𝛿𝜆
, 𝜎𝛿ℎ
, 𝜎𝛿𝑉
, 𝜎𝛿𝑉
, 𝜎𝛿ℎ
, 𝜎𝛼2 , 𝜎𝛽2 , 𝜎𝛾2 - значения СКО ошибок эталонных значений ши𝑛
𝑒
роты, долготы, высоты, скоростей и углов соответственно.
На основе модели (7)-(19) можно организовать процедуру оптимального оценивания
параметров вектора состояния, содержащего и коэффициенты модели ошибок датчиков.
ОФК реализуется в дискретной форме. При этом уравнения дискретной стохастической системы записываются в виде:
𝑋𝑘 = Ф𝑋𝑘−1 + Г𝑊𝑘−1
𝑍𝑘 = 𝐻𝑋𝑘 + 𝑉𝑘
(20)
Воспользуемся условием статической эквивалентности [2], чтобы перейти к интенсивностям дискретных шумов системы Wи шумов измеренийV:
𝑄1 =
𝑄 𝑡
;
𝑇
𝑅1 =
𝑅 𝑡
,
𝑇
(21)
(22)
где R1, Q1– матрицы интенсивностей шумов дискретной системы, T– шаг дискретизации. Для
нестационарной системы интенсивности шумов системы и шумов измеренийQ1 и R1 вычисляются на каждом шаге дискретизации.
Для записи дискретной формы оптимального фильтра Калмана необходимо использовать переходные матрицы динамики системы и шумов системы -Ф и Г, которые представляются в виде матричных рядов от матрицы динамики F и матрицы шумов системы G. В рас9
сматриваемом случае используется разложение в матричный ряд до второго члена (линейного относительно шага интегрирования):
Ф = 𝐸 + 𝐹𝑇;
(23)
Г = 𝐸 + 0.5 ∙ 𝐹𝑇 ∙ 𝐺𝑇,
(24)
где Е – единичная матрица.
Начальная ковариационная матрица 𝑃0 определяется начальным состоянием системы:
𝑃0 = 𝑐𝑜𝑣 𝑥𝑖 𝑡0 ,
(25)
где 𝑥𝑖 𝑡0 элементы вектора состояния 𝑋 при t=0.
Алгоритм дискретного фильтра Калмана взят в форме Джозефа, что обеспечивает симметричность матрицы ковариаций:
𝑥𝑘 = Ф𝑥𝑘−1 + 𝐾𝑘 𝑍𝑘 − 𝐻Ф𝑥𝑘−1
𝑆𝑘 = Ф𝑃𝑘−1 ФТ + Г𝑄1 ГТ
𝐾𝑘 = 𝑆𝑘 𝐻 Т 𝐻𝑆𝑘 𝐻 Т + 𝑅1 −1
𝑃𝑘 = 𝐸 − 𝐾𝑘 𝐻 𝑆𝑘 𝐸 − 𝐾𝑘 𝐻 Т + 𝐾𝑘 𝑅1 𝐾𝑘Т
(26)
Решением этой системы уравнений является оценка вектора состояния системы𝑥, в состав которого входят оценки значений калибровочных коэффициентов.
В ходе испытаний, с целью обеспечения наблюдаемости большинства составляющих
вектора состояния, инерциальную навигационную систему необходимо устанавливать в различные положения. После каждой перестановки системы производится начальная выставка.
Время испытаний выбирается произвольно, но обычно каждая итерация составляет не менее
5 минут. Результатом испытаний являются калибровочные коэффициенты модели ошибок
гироскопов и акселерометров.
Для исследования свойств предложенного алгоритма был разработан программный
комплекс, предназначенный для проведения имитационных исследований процесса калибровки БИНС и проверки работоспособности разработанных алгоритмов. Для реализации
программного комплекса выбрана среда разработки MatLab.
Были проведены серии имитационного моделирования на разработанном программном
10
комплексе. Полученные результаты проанализированы и представлены виде графиков ниже
на рисунках 2-4. В качестве примера, представлены графики ошибок оценок постоянных составляющих погрешности акселерометра и предсказанные их СКО. Предлагаемый в статье
пример рассматривает БИНС на лазерных гироскопах и маятниковых акселерометрах высокой точности.
Рис. 2. Ошибки оценки постоянных погрешностей акселерометров и их СКО.
Первое положение БИНС.
Первым этапом калибровки является установка БИНС на предварительно отгоризонтированную и ориентированную в азимуте планшайбу поворотного стенда и ее начальная
выставка, после чего происходит вращение поворотного стола по определенной траектории.
В качестве примера результатов проведенного имитационного моделирования рассмотрены
графики ошибки оценки постоянных составляющих погрешностей акселерометров и их СКО
на каждом этапе калибровки. В первом положении системы на поворотном основании проекция вектора ускорения свободного падения полностью проецируется на измерительную
ось одного из акселерометров блока. Вследствие этого, исходя из анализа графика на рисунке 2, для одного из акселерометров блока на первом этапе калибровки не вычисляется оценка
постоянной составляющей погрешности. На рисунке 2 видно, что ошибка оценки для двух
акселерометров блока за 100 секунд устанавливается на уровне 0.0001 м/c2.
11
Рис. 3. Ошибки оценки постоянных погрешностей акселерометров и их СКО.
Второе положение БИНС.
На втором этапе калибровки (рисунок 3), после проведения переориентация системы
и ее повторной начальной выставки, происходит оценка ошибки не наблюдавшегося в первом случае акселерометра с точностью порядка 0.001 м/c2. Отметим, что СКО ошибки оценки постоянной составляющей ошибки одного из акселерометров блока несколько выше, чем
у аналогичных датчиков на конец второго этапа калибровки. Графики для третьего положения БИНС на планшайбе динамического симулятора приведены на рисунке 4. Аналогично,
для нового этапа калибровки произведена очередная переориентация БИНС и ее последующая начальная выставка.
Рис. 4. Ошибки оценки постоянных погрешностей акселерометров и их СКО.
Третье положение БИНС.
12
Анализ графика на рисунке 4 показывает, что СКО оценки постоянной составляющей
погрешности акселерометра, оценившейся на предшествующем этапе недостаточно хорошо,
в процессе третьего этапа калибровки существенно улучшилась по сравнению с предшествующим шагом калибровки, достигнув значения аналогичного параметра у других акселерометров блока, оцененных на первом и втором этапах. Таким образом, к концу третьего этапа
калибровки БИНС найдены оценки постоянных составляющих погрешностей акселерометра
с точностью 0.00001 м/с2, удовлетворяющей требованиям калибровки маятниковых акселерометров высокой степени точности. Аналогичный анализбыл проведен для оценок всех
составляющих вектора состояния (7). Имитационное моделирование подтвердило, что
разработанный алгоритм позволяет провести калибровку БИНС с достаточной точностью в
рамках заявленных моделей для перспективных разработок систем ориентации и навигации.
Однако, следует отметить, что данные результаты получены методом имитационного
моделирования, когда многое определяется принятыми математическими моделями, а
настройки оптимального фильтра, как правило, адекватны задаваемым возмущениям. В
рамках будущих работ, с целью проверки работоспособности разработанного алгоритма в
реальных условиях, будет проведено полунатурное моделирование на специализированном
прецизионном оборудовании кафедры 305 Московского авиационного института (НИУ).
Список использованной литературы.
1. Веремеенко К.К., Желтов С.Ю. и др. Современные информационные технологии в задачах
навигации и наведения беспилотных маневренных летательных аппаратов. // М.: ФИЗМАТЛИТ, 2009 г., 556 с.
2. Ривкин С.С., Ивановский Р.И., Костров А.В. Статическая оптимизация навигационных
систем. Л., Судостроение, 1976 г., 280 с.
3. Галай И.А., Зимин Р.Ю. Анализ характеристик двухосного поворотного стенда ком-
пании ActiDyn. // сборник тезисов докладов “Инновации в авиации и космонавтике –
2011”, 2011 г., 151 с.
4. Галай И.А. Программно-математическое обеспечение лаборатории испытаний интег-
рированных навигационных систем. // Тезисы докладов Международного научнотехнического семинара “Современные технологии в задачах управления, автоматики
и обработки информации”, г. Алушта, 2011 г.
5. А.Е. Федоров, Д.А. Рекунов, С.Е. Переляев, Ю.Н. Челноков. Калибровка блока инерциальных чувствительных элементов и моделирование автономного режима функционирования
13
инерциальной системы на базе монолитного трехкомпонентного лазерного гироскопа. // Новости навигации, №3, 2010 г.
6. Тарановский Д.О. Метод калибровки блока маятниковых поплавковых акселеромет-
ров корабельной инерциальной навигационной системы. // С.-Пб.: Автореферат, 2009.
7. Николаев С.Г. Калибровка бесплатформенных инерциальных навигационных сис-
тем. // Приборостроение, 2009 г.
Сведения об авторах.
Веремеенко Константин Константинович, доцент, заместитель декана Московского авиационного института (национального исследовательского университета), к.т.н.
МАИ, Волоколамское ш., 4, Москва, А-80, ГСП-3, 125993, +7-499-158-12-02, e-mail:
nio3@mai.ru
Галай Игорь Анатольевич, аспирант Московского авиационного института (национального
исследовательского университета).
МАИ, Волоколамское ш., 4, Москва, А-80, ГСП-3, 125993, e-mail: igorgalay@me.com
14
Download