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

Автор: Зубов Н.Е., Микрин Е.А., Негодяев С.С., Рябченко В.Н., Богачев А.В., Воробьева Е.А.

Журнал: Труды Московского физико-технического института @trudy-mipt

Статья в выпуске: 4 (20) т.5, 2013 года.

Бесплатный доступ

Рассматривается задача гравитационной разгрузки кинетического момента инерционных исполнительных органов космического аппарата во всех трех каналах управления для круговых орбит с использованием метода точного размещения полюсов. Синтезированы в аналитическом виде законы управления гравитационной разгрузки, однозначно определяемые параметрами объекта и задаваемыми корнями характеристического уравнения.

Инерционные исполнительные органы, разгрузка кинетического момента, обратная связь по состоянию, замкнутая система, метод точного размещения полюсов, ортогональный делитель нуля

Короткий адрес: https://sciup.org/142185944

IDR: 142185944

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

В работе [1] рассмотрена, проблема, разгрузки кинетического момента, инерционных исполнительных органов КА в канале тангажа. Следует заметить, что из-за. автономности данного канала, управления относительно двух других существенно упрощается поиск решения задачи, которая значительно усложняется для гироскопически связанных каналов крен-рысканье. Главное усложнение заключается в том, что для канала, тангажа, система, управления относится к классу SIMO-систем, тогда, как каналы крен-рысканье к MIMO-системам. Предложенные в работах [1, 2, 3] методы синтеза модального управления КА, основанные на ленточной теории [1], обобщенной формулы Аккермана [2] и метода точного размещения полюсов [3] принципиально способны получить аналитическое решение задачи разгрузки кинетического момента для круговых орбит в каналах крен-рысканье с использованием гравитационного момента. Исследования, проведенные авторами в этом направлении, показали, что наиболее предпочтительным в смысле компактности аналитического решения является метод точного размещения полюсов.

В предлагаемой работе с применением указанного метода, модифицированного с учетом результатов работ [4, 5], для круговых орбит получено решение задачи гравитационной разгрузки кинетического момента, инерционных исполнительных органов КА во всех трех каналах управления.

1. Математическая модель КА

Будем рассматривать разгрузку инерционных исполнительных органов (ИИО) без расхода. рабочего тела. Для этой цели используется гравитационный момент [6].

Линеаризованные уравнения для круговых орбит около углового положения (0,0, 0) вращательного движения КА при наличии в каждом канале управления ИИО с учетом гравитационного момента, в соответствии с [6, 7] имеют вид

Jx7 + 4 (Л - Jy) ... + (Jx + Jy - J.) '' = —hy ш о - h_x ,

^J y ^ ^— ⁽ ^J x ⁺ ^J y ^— ^J . ) ^Ш 0 7 ^{+ ( J}. ^— ^J x) Ш^ = ҺхШ о ^— ^ y , C¹-¹)

J.Ө — 3 (Jy — Jx) Ш2Ө = —h, где Jx, Jy, Jz — главные центральные моменты инерции КА; hx, hy, hz — проекции кинетического момента ИИО на оси связанного базиса; шо — орбитальная угловая скорость движения КА (ш = [ 0 0 —шо ]Т); у, ф, Ө — углы отклонения связанного базиса по крену вокруг местной вертикали и от местной вертикали по рысканью и тангажу.

Применяя теорему об изменении кинетического момента отдельно к корпусу КА и отдельно к ИИО на основании системы (1), имеем

^J x 7 ^{+ 4} ( ^J z - ^J y ) ^ш о 7 ⁺ ⁽ ^J x + ^J y - ^J z ) ^ш О ^ф = - ^u x, ^J y ^ф - ⁽ ^J x ⁺ ^J y - ^J z ) ^ш 0 7 ⁺ ⁽ ^J z - ^J x) ^ш о ^ф = -^u y , ^h x ^{+ h} y ^ш О = ^u x, ^h y - ^һ х ^ш 0 = ^u y ,

(1 . 2)

^J z ^Ө - ³ ⁽ ^J y - ^J x) ^ш 0 ^Ө = - ^u z , hz = Uz .

Здесь u_x,Uy,u_z — моменты реакции в подшипниках маховиков, через которые осуществляется как управляющее воздействие на корпус КА с целью поддержания его ориентации, так и одновременное регулирование (в нашем случае) разгрузка накопленного кинетического момента ИИО). Определим вектор состояния КА следующей матрицей-столбцом:

X = [ Ж1

= [ 7 ⁷

Ж2 | Жз | Ж4 | Ж5 | Ж6 | Ж7 | Ж8 | Ж9 | Жщ | Ж_И | Ж 12 ] ^Т =

•

ф ф hx

^h y

J hx dt

t / 0

•

hydt Ө Ө h.

' z

t T

J hzdt .

Тогда система (1.2) с учетом (1.3) в векторно-матричном виде запишется так:

(1 . 3)

(1 . 4)

где

x = Ax + Bu ,

A =

0	1	0	00	00	0	0	0	0	0
« 21	0	0	« 24 0	00	0	0	0	0	0
0	0	0	10	00	0	0	0	0	0
0	« 42	« 43	00	00	0	0	0	0	0
0	0	0	00	ш о 0	0	0	0	0	0
0	0	0	0 — ш о	00	0	0	0	0	0
0	0	0	01	00	0	0	0	0	0
0	0	0	00	10	0	0	0	0	0
0	0	0	00	00	0	0	1	0	0
0	0	0	00	00	0	« 99	0	0	0
0	0	0	00	00	0	0	0	0	0
0	0	0	00	00	0	0	0	1	0
			0	00
			1 J .	00
			0	00
			0 -	- T 0 ^Jy
			1	00
		B =	0	10
			0	00
			0	00
			0	00
			0	0 — T ^J2
			0	01
			0	00

^4ш 0 ,

« 21

024 = шо,

^J x + ^J y

^J y

(1 . 5)

^О 42

—

^J X + ^J y

J y

^{— J} z ------^0 ,

0 43 —

^{— J} z

^J y

W 2 ,

^O 99 ^—

—

J .

Ш 0 .

2. Алгоритм точного размещения полюсов

Задана линейная многомерная динамическая модель КА (1.4), в которой x G R¹², u G R³, R - множество действительных чисел. Матрица B G R^12x3 имеет полный ранг по столбцам (входные сигналы линейно независимы), а матрица A G R^12x12 заведомо неустойчива, т.е. множество ее собственных значений (спектр)

eig (A) — {Ai G C : det (AI12 — A) — 0} , являющееся одновременно множеством корней (полюсов) характеристического полинома, eig (A) — roots(det (AI12 — A), где I12 — единичная матрица размера 12 х 12, C — множество комплексных чисел (комплексная плоскость), обязательно включает такие Аг G C, что Re (Ai) > 0.

Система. (1.4) полностью управляемая, т.е. выполняются эквивалентные условия [3]:

rank ( B I AB I • • • I A9B ) — 12,(2.1)

VAi G eig (A) : rank ( A — АД12 | B ) — 12.(2.2)

Требуется найти закон управления с обратной связью u — — Kx,(2.3)

где K G R3x12 — матрица регулятора по состоянию, чтобы все элементы множества eig (A — BK) — roots (det (AI12 — A + BK))

лежали в открытой левой полуплоскости C, т.е. Re (Ai) < 0. В отличие от [3] введем модифицированную многоуровневую декомпозицию системы (1.4) следующего вида, (для простоты продолжая считать, что все матрицы Bi имеют полный ранг по столбцам):

нулевой уровень A0 — A, Bo — B, (2.4) первый уровень A1 — Bq AoBO--, B1 — Bo- Ao Bo , (2.5) второй уровень A2 — B-A1B--, B2 — B-A1B1, (2.6) третий (конечный) уровень A3 — B-A2B--, B3 — B-A2B2, (2.7) где B- — 2-полуобратная матрица для B- [4, 5], т.е. матрица, удовлетворяющая условиям регулярности

B - B -- B - — B -, B -- B - B -^- — B -^-. (2 . 8)

l ill

Таким образом, справедлива.

Теорема 1. Пусть система (1.4) полностью управляемая и выполнена многоуровневая декомпозиция (2.4) - (2.7), где все матрицы B ^-- удовлетворяют условиям регулярности (2.8), а матрица K G R^rxn удовлетворяет формулам

K — K o — B - A — Ф o B o , B ^- — KX + B + , (2 . 9)

K1 = B-A1 - Ф1B-, B- = K2B1" + B+, (2.10) К = B- A2 - Ф2B2, B- = K3B2" + B2+ , (2.11) Кз = B-Аз - ФзB-, (2.12) тогда выполняется тождество eig (A - BK) = ^ еід(Фі-і).

(2 . 13)

C =1

Поскольку в нашем случае выполняется условие симметричности

(B^-B^ T = B^-B^ , ( B^B^- ^ T = b^^- O C C C

тогда имеем B«±- = B«±+, (2.14) где B^+ — пссвдообратііая матрица. для B^.

Согласно многоуровневой декомпозиции для МІМО-системы (1.4) с матрицами (1.5) вычислим делитель нуля и псевдообратную матрицу для нулевого уровня. Имеем

Для рассмотрения первого уровня декомпозиции найдем псевдообратную матрицу для делителя нуля B^. которая

^ +

B 0 =

в данном случае запишется						так:
1	0	0	0	0	0	0	0	0 ^⎤
0	0	J x	0	0	0	0	0	0
		( J X +1)
0	1	0	0	0	0	0	0	0 ^⎥⎥
0	0	0	^J	0	0	0	0	₀ ⎥_⎥
			( J _y ² +1)
0	0	1	0	0	0	0	0	0 ^⎥⎥
		( J X +1)
0	0	0	1	0	0	0	0	0 ^⎥⎥
			( J _y ² +1)
0	0	0	0	1	0	0	0	0 ^⎥
0	0	0	0	0	1	0	0	₀ ⎥_⎥
0	0	0	0	0	0	1	0	0 ^⎥
0	0	0	0	0	0	0	J x	0 ^⎥
							( J _z ² +1)
0	0	0	0	0	0	0	1	0 ^⎥
							( J 2 +1)
0	0	0	0	0	0	0	0	1

	1	0	0	0	0	0	0	0	0	0	0	0
	⎢ 0	0	1	0	0	0	0	0	0	0	0	0
	⎢_⎢0	Jx	0	0	1	0	0	0	0	0	0	0
	⎢_⎢0	0	0	^J y	0	1	0	0	0	0	0	0
B^ =	⎢_⎢0	0	0	0	0	0	1	0	0	0	0	0
	⎢₀	0	0	0	0	0	0	1	0	0	0	0
	⎢_⎢0	0	0	0	0	0	0	0	1	0	0	0
	⎢₀	0	0	0	0	0	0	0	0	J.	1	0
	0	0	0	0	0	0	0	0	0	0	0	1

⎡₀

(j/+ 1) 0 0 1 ( j X +1) 0 ⁰⁰⁰⁰ 0 0

B o⁺ = ⁰

^{0 0} (j/+ 1) ^{0 1} (/Дд ^{0 0 0 0 0 0}

000 0 0 0 0 0 ( // +1) 1 y^, °J

Следовательно, матрицы		первого	уровня декомпозиции будут выглядеть следующим
образом:	0	0	J x	0	00	0	0	0 "
			(JX+1)
	0	0	0	^J y	00	0	0	0
				(JX+1)
	^Jx^o21	0	0	^ o + J _x J _y a 24 (J2+1)	00	0	0	0
	0	JxO 43	Ш 0 -J_x Jy« 42	0	00	0	0	0
A1 = B^AoB^"⁺=			(JX+1)
	0	0		0	00	0	0	0
			(JX+1)
	0	0	0	1	00	0	0	0
				(JX+1)
	0	0	0	0	00	0	J _z	0
							(J2+1)
	0	0	0	0	0 0 J.099		0	0
	0	0	0	0	00	0	1	0
							(J2+1)
			1 ^- JX	0	0
			0 0	1 ^J y шо - ^Jx^⁴	0 0
			—шо - N⁴²	0	0
B1	= B q A0B0 =		1	0	0	.
			0	1	0
			0	0	1 ■ . . J _z
			0	0	0
			0	0	1

Далее, действуя аналогичным образом, выполним многоуровневую декомпозицию системы (1.4) по формулам (2.5) - (2.7).

Зададим матрицы Ф = Фо, Ф1, Ф2, Ф3, фигурирующие в методе задания полюсов, в следующем диагональном виде:

^/ 01

Ф0 = 0

^/ 21

Ф2 = 0

0 0 1 Г /11 0 0"

/о2 0 , Ф1 = 0 /12 0 ,

0 /03 J L 0 0 /13 _

0 0 1 Г /31 0 0

/22 0 , Ф3 = 0 /32 0 .

0 /23 J L 0 0 /33

Здесь (/01, /11; /21, /31; /02, /12; /22, /32; /03, /13; /23, /33) - в общем случае комплексные числа, подчиненные определенным правилам их формирования [3].

Опуская промежуточные выкладки, основанные на. использовании соотношений (2.10) - (2.12), в соответствии с (2.9) матрица, коэффициентов решения задачи разгрузки кинетического момента. ИИО имеет вид

K = Ko = Bo A0 - Φ0B0- = K11 ^12 ^13 K14 ^15 ^16 K17 0 0 0 0 01 K21 ^22 ^23 ^24 ^25 ^26 0 ^28 0 0 0 0 (2.15) , 0 0 0 0 0 0 0 0 ^39 ^310 ^311 ^312 где

° 21 ⁽ ^J x⁽ ^o 21 + ^/11^/21 + ^/11^/31 + ^/21^/31⁾ + ^J y 042^0)

^U ,.,2

(Щх^шо ⁺ Jy a 42 W 0 ⁺J_xa 21 ⁾

/o1(^J_x°21 °43^(/11 ⁺ ^/ 21 ⁺ ^/ 31 ) ⁺ ⁽ ^J _x ^o 43 ⁺ ^J _x° 24 ° 42

О43( ^ш0 + Jy °42ш0+^х021)

^J y °42^0^)/11^/21^/31^{) ]},

(2.16)

^(J x ⁰ 21 ^{+ J} y ⁰ 42 ^ 0 ^)(/ 01 ^{+ /} 11 ^{+ /} 21 ^{+ /} 31⁾ .

^K12 = ^{— J}x 2,7...., 7 „ X ------⁺

^(J x ^ o ^{+ J} y ⁰ 42 ^ 0 ^{+ J} x ⁰ 21 )

Л^ЛнУ'пУ^ ^{+ /} 01 ^/ 11 ^/ 31 ^{+ /} 01 ^/ 21 ^/ 31 ^{+ /} 11 ^/ 21 ^/ 31 ⁾ 1

+ / T 9 . T . T X ,

(2 . 17)

⁽ ^J x Ш 2 + ^J y 0-42 ^ 0 ⁺ ^J x ⁰ 21 ⁾

_r „ ( J _x 0 24

^ 13 = ^J x ^[/ 01------

^J y ^ш 0 ^)(/ 11 ^/ 21 ^{+ /} 11 ^/ 31 ⁺ / 21 / 31 )

Jy ^43

^(J x ^ o + ^J y ° 42 W o + ^J _x ° 21 )

^J y ^ 0 ⁽ ° 43 ^(/ 11 ^{+ /} 21 ⁺ ^/ 31 ⁾ ⁺ АхЛхЛх )

^J x ⁰ 24 ^/ 11 ^/ 21 ^/ 31

(2 . 17)

^O 43 ^(J x ^ 2 ^{+ J} y ° 42 ^ 0 ^+J x ° 21 ⁾

] ,

^J x ⁰ 24 ⁽⁰ 21 ⁰ 43 ^{— /} 01 ^/ 11 ^/ 21 ^/ 31⁾ .

^A14 = ^Jx^[ „ 2 , 7 „ , , . I \ ⁺

⁰ 43 ^(J x ^ 0 ^{+ J} y ⁰ 42 ^ 0 ^+J x ⁰ 21 )

+ ^J y Ш д ⁽⁰ 43 ⁽⁰ 43 ^{+ 0} 24 ⁰ 42 ^{+ /} 01 ^/ 11 ^{+ /} 01 ^/ 21 ^{+ /} 01 ^/ 31 ^{+ /} 11 ^/ 21 ^{+ /} 11 ^/ 31 ⁺ ^ 21 ^/ 31 ^{) + /} 01 ^/ 11 ^/ 21 ^/ 31 ⁾

^O 43 ^(J x ^ 2 ^{+ J} y ° 42 ^ 0 ^+J x ° 21 ⁾

] ,

K 15 = Д^[

Ш д'(^/ 01 +/ 11 +У 21 + У31) — (^/О1 / 11 ^/ 21+У01 / 11 ^/ 31+У01 ^/ 21 ^/ 31+/п ^/ 21 ^/ 31)

⁽А^ш 0 ^{+ J} y ⁰ 42 ^Ш 0 ^{+ J} x ^o 21 )

(2 . 18)

] , (2 . 19)

K 16 = JxJy ^[

° 43 ⁽ ^/ 11 ^/ 21 ^{+ /} 11 ^/ 31 ^{+ /} 21 ^/ 31

W 2 ) ⁺ ^/ 01 ⁽ ^o 43 ⁽ ^/ 11 ^{+ /} 21 ^{+ /} 31 ^{) +} ЛхАхАх )

^J x ° 24 ^/ 01 ^/ 11 ^/ 21 ^/ 31

^J y O 43 ^(J x ^ 0 ^{+ J} y ° 42 ^ 0 ^+J x ° 21 ⁾

^J y O 43 ^(J x ^ 0 ⁺ ^J y ° 42 ^ 0 ^+J x ° 21 ⁾ ^

K 17 = ^J x ^[

^/ 01 ^/ 11 ^/ 21 ^/ 31 ^(J y ^ p ^{— J} x ⁰ 24 ^ 0 ^{+ J} y ⁰ 43 )

(2 . 20)

^J y ° 43 ^(J x ^ 0 ⁺ ^J y ° 42 ^ 0 ^+J x ° 21 ⁾

] ,

(2 . 22)

^J y ⁰ 42 ^/ 12 ^/ 22 ^/ 32 ^{+ J} x ^ 0 ⁽⁰ 21 ^(/ 12 ^{+ /} 22 ^{+ /} 32 ^{) + /} 12 ^/ 22 ^/ 32 ⁾

²¹ = ^Jy[ T 2 ⁺

^J y ^ш 0 ^J x ⁰ 24 ^ 0 ^+J y ⁰ 43

^(J y O^MM ⁺ ^/12 ^/ 32 ^{+ /} 22 ^/ 32 ^{) + J} x ^ ₀ C ⁰ 21 ⁺ ^/12 ^/ 22 ^{+ /} 12 ^/ 32 ^{+ /} 22 ^/ 32 ⁾

^+/ 02------------------------------------------------------------------

(2 . 23)

^ 22 = ^J y ^[

^J y ^ 0

^J _y 0 42 ( 0 21 0 43 ^{— /} 02 ^/ 12 ^/ 22 ^/ 32 ⁾

⁰ 21 ^(J y ^Ш 0 ^J x ⁰ 24 ^ 0 ^+J y ⁰ 43 )

^J x ⁰ 24 ^ o ^{+ J} y ⁰ 4з

] ,

^J x ^ 0 ⁽ ° 21 ⁽ ° 21 ^{+ 0} 24 ⁰ 42 ^{+ /} 02 ^/ 12 ^{+ /} 02 ^/ 22 ^{+ /} 02 ^/ 32 ^{+ /} 12 ^/ 22 ^{+ /} 12 ^/ 32 ^{+ /} 22 ^/ 32 ^{) + /} 02 ^/ 12 ^/ 22 ^/ 32 ⁾

⁰ 21 ^(J y ^Ш 0 ^J x ⁰ 24 ^ 0 ^+J y ⁰ 43 )

] ,

⁰ 43 ^(J y ( ° 43 ^{+ /} 12 ^/ 22 ^{+ /} 12 ^/ 32 ^{+ /} 22 ^/ 32 ⁾

J_x0 24 W ₀ )

(2 . 24)

^² 3 ^J y ^[ 2

^(J y ^Ш 0 ^J x ⁰ 24 ^ 0 ^+J y ⁰ 43 )

+ ^/ 02 ^(J y ⁰ 21 ⁰ 43 ^(/ 12 ⁺ f 22 ⁺ ^/ 32 ) ^{+ (J} y ⁰ 21 ^{+ J} y ⁰ 24 ⁰ 42 ⁺ ^J x ⁰ 24 W ₀ ⁾ ^/ 12 ^/ 22 ^/ 32^

(2 . 25)

0 21 ( J y Ш 0

^{— J} x ⁰ 24 ^ 0 ^{+ J} y ⁰ 43 )

] ,

^ 24 = ^- ^J y ^[

^J y ( ⁰ 43 ^(/ 02 ⁺ K2 ^{+ /} 22 ^{+ /} 32 ^{) + /} 02 ^/ 12 ^/ 22 ^{+ /} 02 ^/ 12 ^/ 32 ⁺ Л^^^ ⁺ A^^^ )

J x ® 24 ^ o ( / 02 + /12 + /22 + /32)

^J y ^Ш 0

^J x ⁰ 24 ^ 0 ^+J y ⁰ 43

^J y ^Ш 0 ] ,

^J x ⁰ 24 ^ 0 ^+J y ⁰ 43

^ 26 = ^J y ^[

^J x ⁰ 21 ^ 0 ^(/ 12 ^/ 22 ^{+ /} 12 ^/ 32 ^{+ /} 22 ^/ 32 ^— W 2 )

^25 = ^Jy^[ _T / _T 2 7 _n \ ⁺

^J x ⁰ 21 ^(J y ^Ш 0 ^J x ⁰ 24 ^ 0 ^+J y ⁰ 43 )

+ ^/ 02 ^(J y ⁰ 24 ^/ 12 ^/ 22 ^/ 32 ^{+ J} x ^ 0 ( ⁰ 21 ^(/ 12 ^{+ /} 22 ^{+ /} 32 ^{) + /} 12 ^/ 22 ^/ 32 ⁾⁾⁾

^J x ^o 21 ^(J y ^ш 0

(2 . 26)

(2 . 27)

^{— J} x ⁰ 24 ^ 0 ^+J y ⁰ 43 )

] ,

Ш 2 ^(/ 02 ⁺ Л ⁺ ^22 ^{+ /} 32 ^{) —} (A^A²^ ^{+ /} 02 ^/ 12 ^/ 32 ^{+ /} 02 ^/ 22 ^/ 32 ⁺ A2A2A2)

^(J y ^W 0 ^{— J} x ⁰ 24 ^ 0 ^{+ J} y ⁰ 43 )

^K 28 = J_y ^[

/ о2 / 12 / 22 / з2 ( ^ _ж Ш 0 + ^J y 0 42 ^ 0 +^ 0 21 )

^J x ^O 21 ⁽ ^J y ^ш 0 ^J x ^O 24 ^U 0 ⁺ ^J y ^O 43 )

] ,

(2 . 29)

K 3g = J, 0 99 (1 + ^M23 ⁺ ^Ы33 ⁺ 5 23 5 33 ) + J, /_ю ( / 13 + / 23 + / 13 + ^ 13 ^ 23^ ) , (2 . 30)

O gg O gg

/ 03 / 13 / 23 ⁺ / 03 / 13 / 33 ⁺ / 03 / 23 / 33 ⁺ / 13 / 23 / 33

^K 310 = ^- ^J , ^(/ 03 ⁺ / 13 ⁺ / 23 ⁺ / 13 +-- ) , ⁽ 2 . ³¹⁾

K 311

^/ 03 ^(/ 13 ^/ 23 ^{+ /} 13 ^/ 33 ^{+ /} 23 ^/ 33 ^{) + /} 13 ^/ 23 ^/ 33 O gg

^K 312 =

^/ 03 ^/ 13 ^/ 23 ^/ 33

^O gg

(2 . 32)

(2 . 33)

Итак, выражения (2.15) - (2.33) определяют аналитическое решение задачи разгрузки кинетического момента. ИИО во всех трех каналах управления. Они имеют относительно простой и «физичный» смысл и однозначно определяются параметрами орбиты, массоинерционными характеристиками КА, а. также значениями корней характеристического полинома:

det (AI12 - A + BK) = (а - Ах) ••• (а - А12) = Д (а - А^ , г=1

в качестве которых, аналогично, как это было предложено в [3], могут быть использованы корни нормированного полинома. Баттерворта. 12-го порядка.

Соответственно регулятор с матрицей обратной связи (2.15) естественным образом распадается на. два. независимых регулятора: регулятор в каналах крен-рыскание

7 - ^ = [

K 11

K 21

^K 12 ^K 13 ^K 14 ^K 15 ^K 16 ^K 17 ⁰

^K 22 ^K 23 ^K 24 ^K 25 ^K 26 ^{0 K} 28

(2 . 34)

и регулятор в канале тангажа.

^К = [ ^K 3g ^K 310 ^K 311 ^K 312 ] . ⁽² . ³⁵⁾

Следует заметить, что матрица, коэффициентов усиления канала, тангажа. (2.35) с выражениями (2.30) - (2.33), если коэффициенты характеристического уравнения выразить через корни, совпадает с решением, приведенным в [1], что и естественно, поскольку иного быть не может, т.к. канал тангажа, является автономным и принадлежит к классу SIMO-систем.

Преимущества, данного аналитического решения заключаются в следующем. Во-первых, оно пригодно для всех типов КА, дает возможность проводить исследования, варьируя значения полюсов из области их устойчивости. Во-вторых, позволяет ликвидировать проблему плохой обусловленности матрицы А, поскольку для КА значения ее компонент 021, 024 , 042 , 043, Ogg могут отличаться на два и более порядков друг от друга, что, как показали исследования авторов, затрудняет или даже делает практически невозможным получение решения с использованием не символьной, а. числовой матрицы А . Кроме этого, наличие аналитических выражений коэффициентов обратной связи существенно уменьшает затраты вычислительных ресурсов бортовой вычислительной системы.

3. Заключение

В работе для круговых орбит в линеаризованной постановке с использованием метода. точного размещения полюсов получено аналитическое решение задачи гравитационной разгрузки кинетического момента, инерционных исполнительных органов КА во всех трех каналах управления и соответственно показаны преимущества, этого решения, которое не требует настройки под определенный тип КА и может быть легко реализовано в бортовой ЭВМ.

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

Богачев А.В., Воробьева Е.А., Зубов Н.Е., Микрин Е.А., Мисриханов М.Ш., Рябченко В.Н., Тимаков С.Н. Разгрузка кинетического момента инерционных исполнительных органов космического аппарата в канале тангажа//Изв. РАН. Теория и системы управления. -2011. -Т. 3. -С. 125-132
Воробьева Е.А., Зубов Н.Е., Микрин Е.А., Мисриханов М.Ш., Рябченко В.Н., Тимаков С.Н. Синтез стабилизирующего управления космическим аппаратом на основе обобщенной формулы Аккермана//Изв. РАН. Теория и системы управления. -2011. -Т. 1 -С. 96-106.
Зубов Н.Е., Микрин Е.А., Мисриханов М.Ш., Рябченко В.Н. Синтез развязывающих законов стабилизации орбитальной ориентации космического аппарата//Изв. РАН. Теория и системы управления. -2012. -Т. 1
Rao C.R., Mitra S.K. Generalized inverse of matrices and its applications. -N.Y.: John Wiley and Sons, 1971
Ben-Israel A., Greville T.N.E. Generalized Inverses: Theory and Applications. -Springer-Verlag. New York, Inc. 2003
Зубов Н.Е., Микрин Е.А., Мисриханов М.Ш., Рябченко В.Н., Тимаков С.Н., Черемных Е.А. Идентификация положения равновесной ориентации международной космической станции как задача матричного пополнения с устойчивостью//Изв. РАН. Теория и системы управления. -2012. -Т. 2
Раушенбах Б.В., Токарь Е.Н. Управление ориентацией космических аппаратов. -М.: Наука, 1974

Еще

Статья научная