УДК 629.78.051.062.2
управление ориентацией и разгрузка накопленного кинетического момента инерционных исполнительных органов космического аппарата на высокоэллиптической орбите
© 2017 г. Богачев А.в.1, воробьева Е.А.1, Зубов н.Е.12
1 Ракетно-космическая корпорация «Энергия» имени С.П. Королёва (РКК «Энергия») Ул. Ленина, 4А, г. Королёв, Московская обл., Российская Федерация, 141070, e-mail: post@rsce.ru
2Московский государственный технический университет им. Н.Э. Баумана
(МГТУ им. Н.Э. Баумана) 2-я Бауманская ул., 5, г. Москва, Российская Федерация, 105005, e-mail: mail@bmstu.ru
В статье рассматривается динамика углового движения космического аппарата на высокоэллиптической орбите. Разработан алгоритм управления угловым движением космического аппарата и разгрузки кинетического момента инерционных исполнительных органов космического аппарата на высокоэллиптической орбите, использующий для этого гравитационный момент. Для построения законов управления ориентацией и разгрузки накопленного кинетического момента космического аппарата применяется метод точного размещения полюсов. Получено численное решение задачи управления ориентацией и разгрузкой кинетического момента космического аппарата на высокоэллиптической орбите по всем трем каналам управления на участке орбиты, свободном от выполнения целевой задачи, в окрестности перигея. Также произведен анализ устойчивости системы линеаризованных дифференциальных уравнений, описывающих динамику поведения космического аппарата на высокоэллиптической орбите.
Ключевые слова: управление ориентацией, разгрузка кинетического момента, инерционные исполнительные органы, высокоэллиптическая орбита, космический аппарат, гравитационный момент, метод точного размещения полюсов.
attitude control and momentum management for the spacecraft with inertial actuators on high-elliptical orbit
Bogachev A.v.1, vorobyeva E.A.1, Zybov N.E.1, 2
1S.P. Korolev Rocket and Space Public Corporation Energia (RSC Energia) 4A Lenin str., Korolev, Moscow region, 141070, Russian Federation, e-mail:post@rsce.ru
2Bauman Moscow State Technical University (Bauman MSTU)
5 2nd Bauman str., Moscow, 105005, Russian Federation, e-mail: mail@bmstu.ru
The article deals with the dynamics of the angular motion of the spacecraft on the high-elliptical orbit. The attitude control and momentum management algorithm is designed for the spacecraft with inertial actuators on the high-elliptical orbit using gravity gradient torque. Accurate pole placement method is used to develop attitude control and momentum management control laws. A numerical solution is derived for attitude control low and momentum management of the spacecraft on the high-elliptical orbit on all three control channels in the vicinity of the perigee of the orbit. Also, the stability analysis is performed for the linearized differential equations describing the dynamics of the spacecraft on the high-elliptical orbit.
Key word: attitude control, momentum management, inertial actuator, high-elliptic orbit, spacecraft, gravity gradient forces, accurate pole placement method.
БО1А4ЕВ А.В. ВОРОБЬЕВА Е.А. ЗУБОВ Н.Е.
БОГАЧЕВ Алексей Викторович — кандидат технических наук, начальник сектора РКК «Энергия», e-mail: alexey.bogachev@rsce.ru
BOGACHEV Alexey Viktorovich — Candidate of Science (Engineering), Head of Subdepartment at RSC Energia, e-mail: alexey.bogachev@rsce.ru
ВОРОБЬЕВА Екатерина Андреевна — инженер-математик РКК «Энергия», e-mail: ekaterina.vorobeva@rsce.ru
VOROBYEVA Ekaterina Andreevna — Engineer-mathematician at RSC Energia, e-mail: ekaterina.vorobeva@rsce.ru
ЗУБОВ Николай Евгеньевич — доктор технических наук, профессор, главный научный сотрудник РКК «Энергия», декан факультета РКТ МГТУ им. Н.Э. Баумана, e-mail: nikolay.zubov@rsce.ru ZUBOV Nikolay Evgenyevich — Doctor of Science (Engineering), Professor, Chief research scientist at RSC Energia, Dean of the Faculty of Rocket and Space Technology of the Bauman MSTU, e-mail: nikolay.zubov@rsce.ru
Введение
В работе рассматривается движение космического аппарата (КА) с инерционными исполнительными органами (ИИО) на высокоэллиптической орбите (ВЭО). Для таких орбит характерно деление на два участка: рабочий (в окрестности апогея) и дежурный (в окрестности перигея). Далее в работе речь будет идти о дежурном участке, на котором осуществляются вспомогательные динамические операции. Одной из таких операций является разгрузка накопленного на ИИО кинетического момента. При насыщении системы ИИО КА применяют различные способы разгрузки накопленного кинетического момента. Обычно управление центром масс КА и разгрузка ИИО автоматических КА производится с использованием ракетных двигателей малой тяги, в т. ч. электроракетных [1]. Естественно, наиболее предпочтительными способами управления кинетическим моментом ИИО являются безрасходные способы, использующие моменты малых сил, действующих на КА [2-4].
Целью данной работы является решение задачи управления ориентацией и разгрузки кинетического момента КА с ИИО на ВЭО по всем трем каналам управления
с использованием гравитационного момента и с применением алгоритма точного размещения полюсов.
Математическая модель КА
Для вывода уравнений, описывающих динамику углового движения КА, введем две системы координат:
• орбитальная система координат (ОСК), начало которой совпадает с центром масс КА, ось OYОСК направлена вдоль радиуса-вектора, OXОСК лежит в плоскости кеплеровой орбиты и составляет с вектором линейной скорости острый угол, OZОСК дополняет систему до правой. Рабочая орбита и оси ОСК, а также два участка, на которые делится орбита: рабочий и дежурный, представлены на рис. 1.
Рис. 1. Рабочая орбита и оси орбитальной системы координат: — — дежурный участок; — — рабочий участок
• связанная с КА система координат (ССК) с началом в центр масс КА. Для упрощения математических выкладок оси ССК совпадают по направлению с главными центральными осями инерции КА.
Для перехода от ОСК к ССК необходимо выполнить три последовательных вращения (рис. 2). Первое вокруг начального положения оси OZССК (О^СК) — тангажный разворот на угол 0, второе вокруг промежуточного положения оси ОГССК (ОТ) — разворот по каналу рыскания на угол у и вокруг конечного
положения оси СХ^
разворот по крену на угол у.
Рис. 2. Взаимное положение осей орбитальной и связанной систем координат
Ниже представлены уравнения, описывающие движение КА около центра масс.
1. Кинематические соотношения в ква-тернионной форме [5]:
^ОСК Л ОСК° (Юабс)сск ^оХэСК0 ^ОСК, (1)
где Л£ск — кватернион разворота от ОСК к ориентации, требуемой для разгрузки; (®абс)сск — абсолютная угловая скорость в проекции на ССК; (®0)ОСк — орбитальная угловая скорость в проекции на ОСК.
2. Динамические уравнения углового движения КА, несущего вращающиеся массы, при воздействии гравитационного момента [6]:
/юб + юб х /ю,+ Н + юб х Я = М , (2)
^ абс абс ^ абс абс ^аъ ^ '
где ] — тензор инерции КА; И — вектор кинетического момента ИИО; М
гравитационный момент.
Далее:
• проецируем уравнение (2) на оси ССК;
• предполагаем малыми углы крена и рыскания, а также производные по времени от углов крена, рыскания и тангажа;
• применяем теорему об изменении кинетического момента отдельно к корпусу КА и отдельно к вращающимся элементам конструкции;
• производим замену переменных ш0, ю0, г на переменные и параметры кеплеро-вой орбиты (эксцентриситетом, истинной, эксцентрической и средней аномалией);
• линеаризуем полученную систему уравнений около углового положения (0, 0, 0О).
Таким образом, линеаризованные уравнения углового движения КА при наличии ИИО для эллиптических орбит при учете действия гравитационного момента в случае применения теоремы об изменении кинетического момента [7] в соответствии с работами [8, 9] имеют вид:
7уу 7
7хх
о о
(1+ есовV) -ь З-^е сое29С
7 + 7 -7
" хх ^ у у ^ 22
7х
0°-(1+ еС08 +
+
о
2е«ту - 3
А7 -7 VI
уу
7х
\ ^ хх /
-8т2 еп+о сов2 еп
7х
V =
7 7 ^
хх
V 7уу J
0 (1+ есое+ 300 81п20 IV +
7 + 7 - 7"
" хх ^ у у ^ 22
7
уу
Ш1 + есо8-
о
2 ^пу + 3
А7 -7 VI
" хх ^ 22 -I-
7у
V •-> уу /
-8т2 е0+Осо82 е0
2 о о
7у
уу
Н + .. 0(1 + есо%у)Н = и ,
х V Г у х
Н + есо= и ,
у V Г 3 х "
03
7уу 7хх
722
г
соз20^
7уу 7хх
ху \
722
О
51п20 2 —з е51 пу —, 0 г 3 Т
22
Ни
(3)
У
и
х
и
у
У
и
где I , I , I — главные центральные мо-
^ ^ XX ^ УУ ^ 22 г Г
менты инерции КА; к , к , к — проекции
кинетического момента
У 2
ИИО
на оси свя-
занного базиса; и , и , и — моменты реакции
7 X У 2 1 ^
в подшипниках маховиков, через которые осуществляются как управляющее воздействие на корпус КА с целью поддержания его ориентации, так и разгрузка накопленного кинетического момента ИИО; у, у, 0 — углы отклонения связанного базиса от орбитального по крену, рысканию и тангажу, соответственно; & — малый угол отклонения от углового положения (0, 0, 0О) в канале тангажа, & = 0 - 0О; г — радиус-вектор центра масс спутника относительно центра Земли; V — истинная аномалия; е — эксцентриситет орбиты; С = цМ; ц — гравитационная постоянная; М — масса Земли.
Из этих уравнений видно, что пространственное угловое движение КА разделяется на два: плоское движение в канале тангажа и неразделяющееся движение по каналам крена и рыскания.
При переходе в пространство состояний с вектором состояния х перепишем систему (3) таким образом, чтобы элементами матрицы М. в соответствующем канале управления являлись элементы, не относящиеся к этому каналу. Тогда система уравнений (3) распадется на три системы вида:
х. = А.х. + В и +М. , I = х, у, г (4)
г г г г г г 7 > \ /
' 0 ^ "1
А - ^Ч-04-0-Ч : в,- Л
1
0 1 0 0
аг21 0 0 0
0 0 0 0
0 0 1 0
0
V У
х. =
г
/ к1 &
Мг
0
т.п
т
,31
х21
у21
ЛУУ -^22
V Jxx У
Jxx -^22
0
V У
V
0
Ч У
V.
О О
— (1 + есоэу) + 3—с О520п
Г 3 ^ •* г 3 и
Л
V и уу У
ОО
V
«221 г 3
30 г'
ЛУУ хх
V Л22 У
соз(2е0);
тх2Г
Л + Л - Л
хх УУ 22 Jхх
/—+ бСО +
О
+ —
2 бэт V 3
Ч -/VI
УУ 22 1
УУ 22
V ^ хх у
тУ21
Л + Л - Л
хх УУ 22 УУ
— 51П290 4- Осо829П
2 0 0
у;
-00 (1 + есовV)у -
О
2^пV + 3
АI - /V1 ^
•/хх ^ -эт29„ + &С082еп Л о 0 0
V уу а
3О
т
221 2 г 3
ЛУУ - Лхх
УУ хх Л22
О
эт(2е^ 2 — ^ nv;
О
тх31= +бС08 v) кУ;
тУ31 = т!0 (1 + е^)^; т2з1 =
Построение законов управления
В пространстве состояний задача ставится следующим образом. Для линейной динамической системы со многими входами и выходами
х = Ах + Ви,
где х — «-мерный вектор, и — г-мерный вектор входа, матрица Б(п*г) имеет полный ранг по столбцам, а множество собственных значений матрицы А равно eig(А) = (А,, е С : det(АJn - А) = = О, I = 1 ,..., п}, требуется найти такой закон управления с обратной связью и = Кх,
где К — матрица регуляторов по состоянию, чтобы множество собственных значений матрицы замкнутой системы лежало в левой полуплоскости комплексной плоскости. В данном случае мы имеем:
х=А.х. + В и. + М., г = х, у, г,
г г г г г г' ' ^' '
и = К х.
г г г
(5)
В канале тангажа система (5) имеет вид:
х- (А + ВК)х + М.
2-2 2 2* 2 2
Рассмотрим случай, когда свободный член в уравнении для канала тангажа равен нулю:
тг21 = °,
3 с_
2 г 3
ЛУУ - Лх
УУ х
Л2
О
5т(2е0) 2— етпV- 0.
22
Тогда 0о= ^агсзт
4 5
3( 5 5)
УУ ^ XX'
\
-е51 nv
при
У
3 ( 5уу 5XX )
условии | sinv I <
4 1 е
^ 22
Далее будем рассматривать участок траектории в окрестности перигея, где выполняется условие
522
3 (] 5 )
уу
-eslnv
< 1.
(6)
Для каналов крена и рыскания законы управления ориентацией и разгрузкой в соответствии с работой [10] имеют вид:
г = х, у, (7)
где Ег — единичная матрица, а Ni = (1 1 1 1).
Таким образом, системы уравнений для каналов крена, рыскания и тангажа на участке
Г А \Б. \ -1 Г м \
и. -К.Х.+ [К!-Е] ¿ г г 1 ¿ | 1 1 N ~! "о г ! г ~ 0"
траектории, где выполняется условие (6), имеют вид:
Х = (А + ВК)х,, i = х, у, г, (8)
где КГ(К^Ка]Ка]К^.
Для обеспечения асимптотической устойчивости линейной динамической системы с постоянными коэффициентами необходимо и достаточно, чтобы спектр собственных значений замкнутой системы лежал в левой полуплоскости комплексной плоскости. В данной работе матрицы А. систем управлений имеют периодические коэффициенты, а матрицы замкнутых систем А. + В.К. являются матрицами с переменными коэффициентами. В этом случае для обеспечения устойчивости решений систем (8) в работе применялась теорема об асимптотической устойчивости линейной системы с переменными коэффициентами [11].
Для того чтобы применить эту теорему к системам (8), перепишем матрицу системы для каждого канала управления в виде суммы двух матриц: с постоянными Апост и переменными коэффициентами С(£):
х. =
0
К
«¿21л--—
5г
К
А
1 0 0
К 2 5 К.. К4 ¿а т
•• г г К 2 •• гг К3 ¿з • ¿г К4 ¿а
0 1 0
0
0 о
ООО
«¿21- «¿21 п ООО
ООО ООО
о
к- 5
¿1 ^ и
+ в
вгз'
к- 5..
¿2 ^ и
в.
в Л
Кз а
в, в,з
; а
вЮ (fi0 + fi1 + fi2 + fi3)' Pi1 fi0 fi1 + /«./¿2 + + filfi2 + Íi1 /¿3 + -/¿2 -/¿3'
Pi2 fi2 + /¿з + fiofi2fiз + /и/12-/гз); P¿з = /ю /и-/¿З.
Здесь К.. — элементы матрицы регуляторов по состоянию в соответствующем канале управления, вычисленные по формулам метода точного размещения полюсов [12], суть которого заключается в том, чтобы заменить спектр значений замкнутой системы на желаемый спектр; / — в соответствии с методом точного размещения полюсов, собственные числа матрицы А
' 1 ^ пост
(корни характеристического полинома); Рю, Рй, Ра, Рй определяются по формулам, которые выражают коэффициенты многочлена через его корни (формулы Виета); аШп -значение элемента ат в перигее (при V = 0°); i = х, у, г;] = 1, 2, 3, 4.
Из доказательства теоремы вытекает следующий критерий:
4с.. < пгнп
1г
Ее
-1 V
где еи — константа, ограничивающая норму матрицы С(£) и зависящая только от А ; i = х, у, г.
В качестве эталонного многочлена, расположение корней которого и является желаемым при замыкании системы с постоянной матрицей А , выбирался полином Баттер-
1 ^ пос^' 1 1
ворта четвертого порядка, т. е.
1Т11П
Ие /..
Ро,51П
Г \
п
Ч У
а
¡21л
¡21л
¿21 п
¡21п
Таким образом, для соответствующего канала управления рш — радиус полуокружности в левой полуплоскости с центром в начале координат, на которой лежат корни эталонного многочлена, определяется из условия:
р04>
4г
4С и
эш
г \
П
V8,
Определим значение константы е12 для канала танга^ка.
Под нормой матрицы С;(£) понимается следующее:
0 0 0 0
аг21 аг21П 0 0 0
0 0 0 0
0 0 0 0
йг21 аг21П|
\аг21 аг21п1
3 с
3 - 3 Л
^ уу ^ XX \ У
008(20 )
3С
3 - 3 Л
^ уу ^ XX
V У
3С
3 - 3 л
" у у " XX V 3 J
С1, -
6 с
'3 -3 '
" уу " XX
3,
V /
где гп — радиус-вектор центра масс КА относительно центра Земли в перигее. Тогда
24
Рсг > - Г
с
г3 - 3 л
^ уу ^ XX
V 3, ;
(9)
Действуя аналогичным образом, получим оценки радиусов р0х и р0у для каналов крена и рыскания, соответственно:
P0x >
8(4 + е)
81П
V8,
с
3уу - 3,,
уу ,,
Роу >
8(4 + е)
81П
ч8,
с
V JXX '
3XX
(10)
3у
(11)
Результаты моделирования
Математическое моделирование проводилось в среде ЫЛТЬЛБ при следующих условиях:
• начальные значения по кинетическому моменту по всем трем каналам управления принимались равными 10 Нмс;
• моделирование проводилось на участке орбиты, где значение истинной аномалии удовлетворяло условию (6);
• значения радиусов полуокружностей, на которых лежат собственные значения постоянных матриц ^пост для соответствующих каналов управления, выбирались, исходя из ограничений (9-11).
• значения углов отклонения ССК от ОСК по каналам крена и рыскания не превышали 20°, так как система уравнений (3) была получена в предположении малости этих углов. В то же время на угол отклонения в канале тангажа 0 это ограничение не накладывалось.
Численным методом определялось значение радиуса окружности по каждому каналу управления, исходя из ограничений (9-11), соответственно. Из всех возможных вариантов значения радиуса окружности в соответствующем канале управления выбиралось такое, при котором приращение кинетического момента было максимально при соблюдении ограничений в 20° на величину углов отклонения ССК от ОСК по каналам крена и рыскания.
Поиск управления по всем трем каналам осуществлялся с использованием выражения (7), в котором в качестве элементов матрицы М в каналах крена и рыскания брались элементы этих матриц на предыдущем вычислительном шаге. Моделирование объекта управления осуществлялось с использованием динамических уравнений Эйлера (2) и кинематических соотношений в кватернионной форме (1).
Ниже представлены графические результаты моделирования. На рис. 3-5 для каналов крена, рыскания и тангажа представлено по два графика, на верхних — изменение угла отклонения в соответствующем канале в градусах, на нижних — изменение кинетического момента в Нмс по соответствующей оси.
Из рис. 3-5 видно, что кинетический момент КА на участке орбиты в окрестности перигея разгружается, причем по каналам крена и тангажа кинетический момент разгружается полностью. При этом углы отклонения ССК от ОСК по каналам крена и рыскания не превышают 20°.
3
3
г
П
П
г
п
3
П
3
г
П
3
П
Рис. 3. Изменение кинетического момента инерционных исполнительных органов и угол отклонения в канале крена
Рис. 4. Изменение кинетического момента инерционных исполнительных органов и угол отклонения в канале рыскания
Рис. 5. Изменение кинетического момента инерционных исполнительных органов и угол отклонения в канале тангажа
Заключение
В работе произведен анализ устойчивости решения систем линейных дифференциальных уравнений, описывающих динамику поведения КА на высокоэллиптической
орбите. Получены численные решения задачи управления ориентацией и разгрузкой кинетического момента КА на высокоэллиптической орбите по всем трем каналам на участке орбиты в окрестности перигея.
Список литературы
1. Платонов В.Н. Одновременное управление движением центра масс и вокруг центра масс при маневрах космических аппаратов на геостационарной и высокоэллиптических орбитах с использованием электроракетных двигателей // Космическая техника и технологии. 2013. № 1. С. 56-65.
2. Воробьева Е.А., Зубов Н.Е., Микрин Е.А. Безрасходная разгрузка накопленного кинетического момента инерционных исполнительных органов автономного космического аппарата на высокоэллиптической орбите // Инженерный журнал: наука и инновации (электронное научно-техническое издание). Режим доступа: http:// engjournal.ru/catalog/it/nav/1072.html (дата обращения 01.06.2017 г.).
3. Богачев А.В. Управление движением КА «Ямал-100» с использованием светового давления // Ракетно-космическая техника. Труды. Серия XII. Королёв: РКК «Энергия». 2002. Выпуск 1. С. 114-120.
4. Богачев А.В. Управление кинетическим моментом космического аппарата на высокоэллиптической орбите с использованием гравитационного момента и момента сил светового давления // Ракетно-космическая техника. Труды. Серия XII. Королёв: РКК «Энергия». 2010. Выпуск 1. С. 57-60.
5. Браней В.Н. Лекции по теории бесплатформенных инерциальных навигационных систем управления. М.: МФТИ, 2009. 304 с.
6. Белецкий В.В. Движение искусственного спутника относительно центра масс. М.: Наука, 1965. 416 с.
7. Раушенбах Б.В., Токарь Е.Н. Управление ориентацией космических аппаратов. М.: Наука, 1974. 600 с.
8. Зубов Н.Е., Микрин Е.А., Рябченко В.Н. Матричные методы в теории и практике систем автоматического управления летательными аппаратами. М: МГТУ им. Баумана, 2016. 666 с.
9. Тимаков Н.С. Исследование управляемого углового движения космического аппарата на высокоэллиптической орбите // Сб. трудов IX конференции молодых ученых «Навигация и управление движением», 13-15 марта 2007 г., Санкт-Петербург. СПб.: ЦНИИ «Электроприбор», 2007. С. 330-336.
10. Справочник по теории автоматического управления / Под ред. А.А. Красовско-го. М.: Наука, 1987. 712 с.
11. Беллман Р. Теория устойчивости решений дифференциальных уравнений. М.: Иностранная литература, 1954. 216 с.
12. Зубов Н.Е., Микрин Е.А., Мисриха-нов М.Ш., Рябченко В.Н. Синтез развязывающих законов стабилизации орбитальной ориентации космического аппарата / Известия РАН. Теория и системы управления. 2012. № 1. С. 92-108.
Статья поступила в редакцию 06.04.2017 г.
Reference
1. Platonov V.N. Odnovremennoe upravlenie dvizheniem tsentra mass i vokrug tsentra mass pri manevrakh kosmicheskikh apparatov na geostatsionarnoi i vysokoellipticheskikh orbitakh s ispol'zovaniem elektroraketnykh dvigatelei [Simultaneous control of centre of mass and around centre of mass motion during SC maneuvers performance on geosynchronous and high elliptic orbits using electric rocket engines]. Kosmicheskaya tekhnika i tekhnologii, 2013, no. 1, pp. 56-65.
2. Vorob'eva E.A., Zubov N.E., Mikrin E.A. Bezraskhodnaya razgruzka nakoplennogo kineticheskogo momenta inertsionnykh ispolnitel'nykh organov avtonomnogo kosmicheskogo apparata na vysokoellipticheskoi orbite. Inzhenernyi zhurnal: nauka i innovatsii (elektronnoe nauchno-tekhnicheskoe izdanie) [Desaturation without propellant consumption of accumulated kinetic momentum of inertial control elements of a spacecraft free-lying in a high elliptical orbit // Engineering journal: science and innovations (electronic science and technology publication)]. Available at: http://engjournal.ru/ catalog/it/nav/1072.html (accessed 01.06.2017).
3. Bogachev A.V. Upravlenie dvizheniem KA «Yamal-100» s ispol'zovaniem svetovogo davleniya [Controlling the motion of Yamal-100 SC using light pressure]. Raketno-kosmicheskaya tekhnika. Trudy. Ser. XII. Korolev: RKK «Energiya»publ., 2002, issue 1,pp. 114-120.
4. Bogachev A.V. Upravlenie kineticheskim momentom kosmicheskogo apparata na vysokoellipticheskoi orbite s ispol'zovaniem gravitatsionnogo momenta i momenta sil svetovogo davleniya [Controlling angular momentum of a spacecraft in high elliptical orbit using gravitational torque and the moment of light pressure force]. Raketno-kosmicheskaya tekhnika. Trudy. Ser. XII. Korolev: RKK «Energiya» publ., 2010, issue 1, pp. 57-60.
5. Branei V.N. Lektsii po teorii besplatformennykh inertsial'nykh navigatsionnykh sistem upravleniya [Lectures on the theory of strap down inertial navigation control systems]. Moscow, MIPT publ, 2009.304 p.
6. Beletskii V.V. Dvizhenie iskusstvennogo sputnika otnositel'no tsentra mass [Artificial satellite motion about its center of mass]. Moscow, Nauka publ., 1965. 416 p.
7. Raushenbakh B.V., Tokar' E.N. Upravlenie orientatsiei kosmicheskikh apparatov [Spacecraft attitude control]. Moscow, Nauka publ., 1974. 600p.
8. Zubov N.E., Mikrin E.A., Ryabchenko V.N. Matrichnye metody v teorii i praktike sistem avtomaticheskogo upravleniya letatel'nymi apparatami [Matrix methods in the theory and practice of automatic control systems for lying vehicles]. Moscow, MGTUim. Baumanapubl., 2016. 666p.
9. Timakov N.S. Issledovanie upravlyaemogo uglovogo dvizheniya kosmicheskogo apparata na vysokoellipticheskoi orbite. Sb. trudov IX konferentsii molodykh uchenykh «Navigatsiya i upravlenie dvizheniem» [A study of controlled angular motion of a spacecraft in a high elliptical orbit. Collected papers of the IX conference of young scientists Navigation and Motion Control]. 13-15 March 2007, Saint-Petersburg. Saint-Petersburg: TsNII «Elektropribor» publ., 2007. Pp. 330-336.
10. Spravochnik po teorii avtomaticheskogo upravleniya [Handbook of automatic control theory]. Ed. A.A. Krasovskiy. Moscow, Nauka publ., 1987. 712 p.
11. Bellman R. Teoriya ustoichivosti reshenii differentsial'nykh uravnenii [Theory of stability of differential equations solutions]. Moscow, Inostrannaya literaturapubl., 1954.216p.
12. Zubov N.E., Mikrin E.A., Misrikhanov M.Sh, Ryabchenko V.N. Sintez razvyazyvayushchikh zakonov stabilizatsii orbital'noi orientatsii kosmicheskogo apparata [Synthesis of uncoupling laws of stabilization of orbital attitude of a spacecraft]. Izvestiya RAN. Teoriya i sistemy upravleniya, 2012, no. 1, pp. 92-108.