Определение оптимальных конфигураций звезд, визируемых системой астроориентации космических аппаратов Текст научной статьи по специальности «Математика»

УДК 629.78.018:621.397.13

ОПРЕДЕЛЕНИЕ ОПТИМАЛЬНЫХ КОНФИГУРАЦИЙ ЗВЕЗД, ВИЗИРУЕМЫХ СИСТЕМОЙ АСТРООРИЕНТАЦИИ КОСМИЧЕСКИХ АППАРАТОВ

Б.М. Суховилов

DEFINITION OF THE OPTIMAL STAR CONFIGURATIONS SIGHTED BY THE ASTROORIENTATION SYSTEM OF SPACECRAFT

B.M. Syhovilov

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

Ключевые слова: астроориентация, астровизир, конфигурация звезд.

The methods of increasing the pointing accuracy of the target facility of spacecrafts are offered. They take into account the type of the problems solved with the help of the target facilities. The offered methods are based on selection of the pointing accuracy of the target facility criterion and definition of optimal for these criterion star configurations sighted by the astroorientation system of spacecrafts.

Keywords: astroorientation, astrovizir, star configuration.

В связи с повышением требований к точности ориентации целевой аппаратуры (ЦА) космических аппаратов (КА) находят широкое применение системы астроориентации КА, обеспечивающие в настоящее время наивысшую точность. При этом дальнейшее улучшение точности ориентации ЦА КА связано не только с повышением приборной точности астродатчиков системы астроориентации КА, но и с выбором оптимальных конфигураций звезд, визируемых АД, для заданных критериев наведения целевой аппаратуры КА.

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

В работе [1] рассмотрен класс статистических критериев точности наведения ЦА КА, среди которых наибольший интерес представляют: обобщенный критерий К0, учитывающий погрешности наведения ЦА относительно всех трех осей связанной системы координат (ССК) КА, и частный критерий , применимый для ЦА, имеющей осевую симметрию (антенны, телескопы и т.д.), и учитывающий только модуль возможного откло-

нения ЦА от заданного направления. Формулы, определяющие данные критерии, записываются в следующем виде:

Ko=Sp(P£); (1)

К, = m[|AV|2] = VT [Sp(PE)I3x3 - Р,] V, (2)

где Ps - ковариационная матрица вектора е, определяющего отклонения ССК КА от своего истинного положения; Sp(-) - след матрицы; М[-] -

математическое ожидание; |AV| - модуль вектора

отклонения оси ЦА от истинного положения; V -вектор направляющего косинуса оси ЦА в ССК КА; 13х3 - единичная матрица размером 3><3.

В работе [2] для случая использования равноточных АД получено значение ковариационной матрицы РЕ в следующем виде:

pË=-y(M3x3-BBTr\ (3)

где о - среднеквадратичное отклонение (СКО) погрешности измерения АД направления на звезду; N - количество визируемых звезд;

- матрица, составленная из векторов направляющих косинусов визируемых звезд;

Суховилов Борис Максович - д.т.н., заведующий ка- Syhovilov Boris Maksovich - PhD, head of informatics

федрой информатики ЮУрГУ; suhov@mf.susu.ac.ru. department of SUSU; suhov@inf.susu.ac.ru.

Р? - направляющий косинус i-й звезды, визируемой АД, в ССК КА.

Как следует из рассмотрения формулы (3), значение ковариационной матрицы РЕ и соответственно критериев К0 и Kj зависят не только от приборной точности АД, определяемой ст, но и от конфигурации визируемых системой астроориентации КА звезд, определяемой значением матрицы (М3х3 -ВВТ). Осуществляя подстановку ковариационной матрицы Р8 из (3) в формулы (1) и (2), решим задачу определения оптимальной конфигурации визируемых системой астроориентации КА звезд последовательно по отношению к критериям К0 и К,. Решение указанных задач позволит оптимальным образом расположить комплект АД на борту КА, а также назначить эффективный план астроизмерений в процессе определения ориентации целевой аппаратуры КА.

Сначала рассмотрим два важных частных случая компоновки системы астроориентации КА двумя и тремя АД, каждый из которых визирует одну звезду. Случай двух АД, визирующих по одной звезде, - это вариант минимального состава измерителей, при котором можно определить трехосную ориентацию КА. Случай трех АД соответствует варианту оценки трехосной ориентации КА на основе широко применяемых простых эвристических алгоритмов, например, таких, как TRIAD Algorithm [3].

Пусть система астроориентации аппарата оснащена двумя АД, каждый из которых визирует одну звезду. Определим значение критерия К0.

Для этого рассмотрим матрицы ВВТ и ВТВ . Ранг матриц ВВТ и ВТВ равен 2. Матрица ВТВ имеет размерность 2x2 и, следовательно, два отличных от нуля собственных значения. Матрица ВВТ имеет размерность 3x3, и поскольку её ранг равен двум, то одно из собственных значений матрицы ВВТ равно нулю. Известно, что не нулевые собственные значения матриц ВВТ и ВТВ всегда совпадают [4].

Проведем разложение матриц ВВТ и ВТВ по собственным векторам и собственным значениям:

BTB = U2jUT; BBt=WE2Wt, (4)

где U, W - матрицы, составленные из собственных векторов соответственно матриц ВВТ и ВТВ, где

1 0 0'

^1 0 0'

II Л- 2 .s2 = 0 Х2 0

L* 0 0 0

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

Не нулевые собственные значения матриц ВВТ и ВТВ равны. Проще найти собственные

ВТВ =

значения матрицы ВТВ, чем матрицы ВВТ. Кроме того, элементами матрицы ВТВ являются физически просто трактуемые величины - косинусы углов между визируемыми звездами:

1 с12

_с12 1 .

где с12 - косинус угла между звездами, визируемыми первым и вторым АД.

Найдем собственные значения ВТВ из характеристического уравнения с1е<ВтВ-М2х2) = 0, где 12х2 _ единичная матрица размером 2x2.

det

1-Х с12

с12 1-Х

= 0.

Решая (5), получаем: (^1 = 1 ~ с12;

%2 ~ 1 ^\2 •

(5)

(6)

При N = 2 с учётом того, что 21зхз = ^(213х3)\Ут, запишем значение критерия К0 в следующем виде:

2

=-у8р[(213х3 “£2)-1]- (7)

В (7) учтено, что преобразование подобия (213х3 ~ ^2 )_1 \^Т не изменяет след матрицы. Подставляя в формулу (7) значения Хи Х2 из (6), получим окончательное выражение для критерия К0 в случае оснащения системы астроориентации двумя АД в виде:

2

V — g *0“ 2

2 1 -cf2

(8)

Из (8) следует, что значение К0 минимально при с12 = 0, что соответствует случаю, когда АД визируют звёзды во взаимно перпендикулярных направлениях.

Определим значение критерия К0 при оснащении системы астроориентации КА тремя АД. В этом случае матрица ВТВ запишется в следующем виде: сп с13 1 с23

ВТВ =

1

С,2 _с13 с23

1

а след матрицы (313х3-ВВт) 1 составит:

-1

1

1

1

Sp (313х3-ВВт =-------+-

LV ’ J 3-Л.1 з-а,2 з-хъ

— ^27 — + Х2 + А,3)+(А,]А,2 + A,jA,3 + А,2А,3))х

х^27—9( А.| +Х-2 + А*3)+3(А^Х2 + A.jA,3 + А,2А,3) — —A.jA*2A*3 ) ,

(9)

где X.J, Х2, А-з - собственные значения матриц

ВВти ВТВ.

Запишем характеристическое уравнение для определения собственных значений Д2 Д3:

1-Х с12 с13

det с12 1-Х с23 = х3-зх2 +

с13 с23 1-Я.

+А,(3 — с\з — С%2 — С\ъ ) + С\2 + с23 + с13 ~

2^2^13^23 1 — 0.

Используя известные зависимости между корнями полинома и его коэффициентами [4], имеем:

^*1 ^2 ^3 =

Х^Х2 + Л-2^3 "Ь = 3 — С[2 — с2з — %; (Ю)

Х,]^2^3 = 1 + 2с12с2зс13 — С]22 — с|3 - С]23.

Используя соотношения (9) и (10), запишем значение критерия К0 в виде:

^о=—•

9 4- A,jA,2 ^“1 ^-3 "Ь ^-2^”3 2 3(А,|А.2 + A-jA-з + А-2А<з) — XjA/2A.3

1 о Л2 _2 _2

12-С12 -С13 -с23

2 2 2 “• (П) 2 8—2(с12 + с23 + С]3 )—2с12с13с2з

Взяв частные производные по переменным с12< с13> с23 0Т функции (11) и приравняв их к нулю, получаем систему нелинейных уравнений, решение которой составляет условие минимума критерия К0 в виде:

с12 = 0;

'13

= 0;

с2 з = 0.

Таким образом, в случае оснащения системы астроориентации тремя АД, каждый из которых визирует одну звезду, для обеспечения минимума критерия К0, звезды должны визироваться АД во взаимно перпендикулярных направлениях.

Получим общее выражение для критерия К0 при оснащении системы астроориентации КА N астродатчиками.

Для этого случая характеристический полином матрицы ВТВ представляет:

det(Ы№<лr-В1 В) = X -С1ХХ"~' +а2*. -

~а3Хы~3 +...+(-1)ыаы =0.

Поскольку коэффшщент аг характеристического полинома матрицы ВТВ равен сумме всех главных миноров порядка г этой матрицы, значения коэффициентов ах, а2, аг могут быть представлены в следующем виде:

(12)

а, = 8р(В В) = Ы;

ЛГ-1 N

«2 = 1 Ю-ф;

/=1 ;=/+1

ЛГ-2 ЛГ-1 N

аз ~ ]Е Е Е (1"*'2с)ус,^сд — Су — с(£ — Сд).

1=1 J=l+\ к=7+1

Используя известные зависимости между корнями Хх,Х2...ХИ характеристического полинома и его коэффициентами, с учетом того, что Х) = 0,1 = 4...И, имеем:

+ Х2 + ^з = \

^1^2 + Х2Х3 + ^1^з = а2\ (13)

^^2^3 = аз'

Сравнивая (12) и (13), получим соотношения между Хх,Х2, Х3 и косинусами углов с,у между/-й

иу-й звездами в виде:

А-1 + Х2 + Х3 =

ЛМ N

Х{Х2 +А,2А,3 ч-^А-з = ^ ^ (1 -Су);

,=1 7=/+1

Х^Х2Х3 —

N-2 N-1 N

= Е Е Е +2счс>кс1к - 4 - 4 - <7* )•

/=1 у=|+1 к=]+\

Используя формулы (1), (3) и (4), запишем значение критерия К0 в случае оснащения системы астроориентации N АД в виде:

2 (

(14)

*п=-

1

1

1

\N~X\

ЗУ

(15)

2 + А.^А»з + Х2Х3 ) — Л-2^3

Осуществляя подстановку из (14) в (15), получаем общее выражение (16) для критерия К0 при оснащении системы астроориентации КА N астродатчиками.

Подстановкой N = 2, N = 3 можно убедиться, что значение К0, вычисленное по формуле (16), и соответственно формулам (8) и (11) идентичны.

ЛГ-1 N

*2 + Е Ео-Ф ,=i ,=,+1

Nx

к«=-

N-1 ЛГ

хЁ Е о-Ф-

;=1 7=,+1

N-2 ЛГ-1 N N

“ZEE (1+2су%сд-^-4-сд)

(=1 j=i+IK=j+\ J

Для исследования оптимальных свойств критерия К0 была разработана программа в среде MATLAB. В качестве алгоритма оптимизации использовался метод деформируемого многогранника Нелдера-Мида. В качестве переменных опти-

(16)

мизации были приняты угловые сферические координаты звезд, расположенных на сфере единичного радиуса. Для исключения неопределенности решения система звезд фиксировалась на сфере посредством приравнивания нулю двух угловых сферических координат 1-й звезды и одной угловой сферической координаты 2-й. В результате численной оптимизации при последовательно увеличиваемом количестве звезд было установлено, что минимальное значение критерия К0, равное

9 ст2

К0 =-----, достигается при условии А., = Х2 = .

4#

Аналогичный результат получен в работе [2].

Отметим, что разработанная программа позволяет рассчитать оптимальные конфигурации произвольного количества визируемых звезд или бортовую конфигурацию произвольного количества АД, каждый из которых визирует одну звезду. В таблице приведен пример рассчитанных оптимальных конфигураций визируемых звезд в количестве от 2 до 10. Конфигурация определяется сферическими координатами широты ф и долготы 0 направления визирования АД звезды в ССК КА.

Данные в таблице могут использоваться для оптимального размещения на борту КА нескольких АД, участвующих в определении ориентации КА по критерию К0, учитывающему погрешности наведения ЦА относительно всех трех осей ССК КА.

Перейдем к рассмотрению критерия К] и решим задачу определения оптимальной по отношению К критерию К! конфигурации звезд, визируемых системой астроориентации КА. Критерий К] равен дисперсии модуля вектора отклонения оси ЦА от истинного направления. Он полезен, как отмечалось ранее, для анализа задач, выполняемых широко распространенным классом целевой аппаратуры КА, имеющей осевую симметрию.

Запишем значение критерия К, в терминах собственных значений йг (с/, >й1>с1ъ'>Щ

матрицы ковариации Ре в следующем виде:

й?2 + 0 0

К^У7-»!/ 0 ¿,+¿3 0 ч»т V, (17)

0 0 <1\ + й?2

где \|» - матрица, составленная из собственных векторов матрицы Ре.

Поскольку критерий (17) представляет собой квадратичную форму, минимальное значение К,:

(18)

и достигается, когда вектор направляющего косинуса V оси ЦА совпадает с направлением собственного вектора 1|/], соответствующего наибольшему собственному значению <1Х матрицы ковариации РЕ. С другой стороны, минимальное зна-

чение критерия Kj с учетом того, что РЕ = ^-f \ где f = (М3х3 -ВВТ), составляет: minjK.J =Sp(P,)-y4.frV, -

-w>-v£.

1

(19)

где £3 =-----наименьшее собственное значение f.

d\

Следовательно, критерий Kj минимален, когда v|/, = V и при этом \|/i также является собственным вектором, соответствующим наименьшему собственному значению матрицы f.

Поскольку аналогией матрице f в механике является матрица инерции системы материальных точек, расположенных на сфере единичного радиуса, и имеющих направляющие косинусы звезд, составляющие столбцы матрицы В, то вектор vj/, совпадает с главной осью инерции, относительно которой момент инерции указанной системы материальных точек минимален.

Аналогия с механической системой материальных точек позволяет достаточно просто конструировать оптимальные по отношению к критерию К] конфигурации визируемых АД звезд. Введем ограничение на максимальное значение косинуса угла между визирными осями /-го и j-го

АД в виде max {с,-, } < с.

При этом в случае оснащения системы астроориентации, например, двумя АД, каждый из которых визирует одну звезду, оптимальной конфигурацией являются две звезды, косинус угла между которыми равен максимально возможному значению -с, а ось ЦА представляет ось симметрии для векторов направляющих косинусов этих звезд.

Определим min{Ki} при N = 2. Как следует из (6), минимальное собственное значение матрицы f равно 4з = 1—И • Поэтому минимальное значение критерия Kj с учетом (8), (19) в случае оснащения системы астроориентации двумя АД, каждый из которых визирует одну звезду, составит:

2 (-min{Ki}=Y

1

2 1-

1

1-с

о

~2

ui _L'

2 + 1 + |с|

Пусть система астроориентации оснащена N АД, каждый из которых визирует одну звезду. Снимем ограничение на максимальное значение косинуса угла между визирными осями АД и определим, в этом случае, минимальное значение критерия К].

Обозначим ц2> Из (щ ^ Цг - ^з - ~ с°б-

ственные значения матрицы ВВТ.

Конфигурации звезд, оптимальные для критерия К0

Количество АД Угол, ° АД 1 АД 2 АДЗ АД 4 АД 5 АД 6 АД 7 АД 8 АД 9 АД Ю

2 Ф 0 0 — - - - - - - -

0 0 90 — - — - - - — -

3 Ф 0 0 90 — — — — - — -

0 0 90 90 - - - - - -

4 Ф 0 0 120 60 - - - - - -

0 0 109,47 109,47 250,53 — - - — — -

5 Ф 0 0 125,02 28,558 74,875 - — — - —

0 0 66,549 90,484 131,53 254,8 - - - -

6 ф 0 0 95,822 44,059 88,033 0 - - - -

0 0 90,299 77,271 159,42 269,14 74,134 - - — -

7 Ф 0 0 120,83 39,022 90,837 -0,002 69,135 - - -

0 0 94,767 75,248 141,47 290,72 61,399 57,085 - — -

8 Ф 0 0 129,15 25,322 94,901 0 52,939 77,674 -

0 0 54,722 91,954 106,81 280,78 117,26 10,54 78,911 -

9 Ф 0 0 121,49 122,2 77,072 0 10,429 56,201 53,911 —

0 0 98,15 105,21 112,76 282,17 108,32 -26,298 55,796 45,085 -

10 Ф 0 0 133,62 157,45 43,587 -0,018 -53,68 60,277 71,009 95,382

0 0 135,69 116,13 117,58 294,92 115,1 54,842 65,609 73,114 46,774

Поскольку Ц] + ц2 + Цз = N и Р-19 М-2’ М^з “ п0* ложительные числа, собственные значения матрицы {~ %\> %г< $з (§1 ^ \г ^ 4з ^ °). равные ^ = ЛГ-ц3, \2 = N -р.2, £з = > изменяются в диапазоне

от 0 до N в зависимости от конфигурации визируемых системой астроориентации КА звезд. При этом будет минимально, когда щ -»■ N, а ц2 0» ц3 -> 0. Это имеет место, когда визирные оси всех АД близки к совмещению с осью ЦА. Собственные значения матрицы Г в этом случае составляют:

Хх=и+^,

■ \2 =И + Ш2\ (20)

Дз=ДЛГ3,

где Ш, ,¡ = 1,2,3 - бесконечно малые числа.

Поскольку между собственными значениями матриц РЕ и (существуют следующие соотношения:

. °2 1 •

«1 =-------

' 2

л о1 1 •

а? = — •—;

2 2 $2

(I =— —

3 2 4/

то минимальное значение К) с учетом (18), (20) при Ш, —> 0 :

1

уЫ + Щ

N + AN-

2 у

ст

тш{К,1= lim —

№¡->0

Таким образом, при оснащении системы астроориентации КА N астродатчиками, минимум дисперсии модуля вектора отклонения оси ЦА от истинного направления достигается, когда все АД визируют звезду в направлении, противоположном или совпадающим с направлением оси ЦА КА.

Литература

1. Коффи, Т. С. Критерий ошибки ориентации полезной осесимметричной нагрузки КА / Т. С. Коффри //Астронавтика иракетодинамша. -1984. -№ 20. - С. 25-32.

2. Катаргин, М. Ю. Алгоритм среднеквадратичной оценки ориентации космических аппаратов / М. Ю. Катаргин // Космические исследования. - 1986. — Т. 24, вып. 6. — С. 826-830.

3. Shuster, М. D. Tree-Axis Attitude Determination from Vector Observations / M D. Shuster, S. D. Oh // Journal of Guidance and Control. - 1981. - Vol. 4, Jbl.-P. 70-77.

4. Воеводин, В. В. Матрицы и вычисления / В. В. Воеводин, Ю. А. Кузнецов. — М.: Гл. ред. физ.-мат. лит-ры. -1984. - 320 с.

5. Смирнов, В. И. Курс высшей математики: т. 1 /В. И. Смирнов. - М.: Наука. -1974. - 497 с.

Поступила в редакцию 27 сентября 2007г.