Файл: Шаталов, В. А. Применение ЭВМ в системе управления космическим аппаратом.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 21.10.2024

Просмотров: 135

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

где матрица (а) определяется

элементами системы

уравне­

ний (3. 101).

 

 

 

Согласно уравнениям (3. 106),

(3.107) вектор измерения сол­

нечного датчика может быть представлен как

 

5 = 5 В

 

 

(3. 108)

Так как кососимметричная матрица

направляющих косинусов

в уравнении (3. 108) является особой

(ее определитель

равен

нулю), то уравнение (3. 108) не может быть решено относительно углов у, Ф, -ф. Физически это означает, что солнечный датчик определяет лишь направление на'Солнце, но не позволяет опре­ делить угол поворота КА относительно этого направления. Из уравнения (3. 108) может быть получена с высокой точностью линейная комбинация углов у, О, O', что позволит уточнить угло­ вое положение КА.

В цифровой системе отсчета показания солнечного датчика используются наравне с показаниями датчиков горизонта и угло­ вой скорости.

Вектор состояния цифровой системы ориентации с солнечным

датчиком имеет вид

 

 

S0 = SQ{у, &, <!>, ув, Ъв , «>т, <о», «ч, Уг, ft,., УА, Уа,

&д),

109)

 

(3.

где ув, “&в — составляющие вектора S0, наблюдаются при помо­ щи солнечного датчика.

Остальные составляющие соотношения (3. 109) описаны при пояснении зависимости (3. 105).

Вектор измерений связан с переменными состояния соотно­ шением (3. 105). Например, при учете измерений только от дат­ чика горизонта и солнечного датчика вектор измерений углового положения КА

 

Е = Е { ут, #г, Sx, S y, Sz),

(3.110)

где

уг, 'О'г — показания датчика горизонта;

 

 

Sx, Sy, Sz — разности

между направляющими

косинусами

 

осей х',

у', г ’, измеренными солнечным датчи­

 

ком в связанной системе координат и предвари­

 

тельно рассчитанными в орбитальной системе

 

отсчета

на основе информации

об эфемери­

 

дах КА.

 

 

Матрица Ни в данном случае имеет вид

104


1

 

0

0

1

0

0

0

0

1

0

1

1

0

0

0

0

 

1

0

0

1

0

0

0

0

1

0

 

0

1

0

0

*^0г

$ 0 У

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

~ $ 0 х

*-*3X2

*■ *3X2

^ 3 X 2

!

^ 3 X2

*-*3X2

$ 0 У

 

$ 0 х

0

 

 

 

 

 

 

 

 

 

 

 

 

(3.111)

а ковариационная матрица вектора ha помех считается заданной. Для вычисления элементов матрицы Н8 для текущего момен­ та времени ^ необходимо определить косинусы вектора направ­

ления на Солнце в орбитальной системе координат. Из рис. 3. 32 следует

50^.=cos ту, cos i cos — sin ii3 cos (Зэ;

Soy=

— sin i siri рэ;

(3. 112)

S0z=

sin Ti3 cos / sin рэ— cos т)эcos Рэ.

 

Углы т|э и рэ отсчитываются от плоскости эклиптики, а не от экваториальной плоскости.

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

Рис. 3. 32. К

определению ориентации Солнца

Рис. 3. 33. Структурная схема цифровой

8 орбитальной

системе координат

системы определения углов

105


и датчиков угловых скоростей. Структурная схема цифровой системы определения углов приведена на рис. 3. 33.

Анализ приведенных методов ориентации при помощи систем, включающих БЦВМ, показывает следующее:

математическое обеспечение БЦВМ должно предусматри­ вать возможность действий с матрицами (их умножение, транс­ понирование, обращение) и тригонометрическими функциями;

— хранение элементов матриц предъявляет к ЗУ повышенные

требования.

текущей ориентации

КА

решаются,

Задачи

определения

в реальном масштабе времени.

 

 

 

 

Т а б л и ц а

3.1

 

 

 

 

 

 

 

 

С р е д н я я к в а д р а т и ч е с к а я о ш и б к а ,

 

 

 

 

р а д ( г р а д )

 

 

С и с т е м а

 

 

 

 

 

 

 

 

о с ь к р е н а

о сь т а н г а ж а

о с ь р ы с к а н и я

Обычная

аналоговая гирокомпас-

0,00174

0,00437

0,00384

ная

 

 

( 0 , 1 0 )

(0,25)

( 0 , 2 2 )

Цифровая

без солнечного

датчика

0,00157

0,00419

0 , 0 2 1

 

 

 

(0,09)

(0,24)

( 0 , 1 2 )

Цифровая

с солнечным датчиком

0,000523

0,00035

0,000523

 

 

 

(0,03)

( 0 , 0 2 )

(0,03)

Т а б л и ц а

3. 2

 

 

 

 

 

Ошибки датчика горизонта

 

 

 

 

 

 

И с то ч н и к

 

 

 

В ел и ч и н а ,

 

 

 

 

р а д ( г р а д )

 

 

 

 

 

 

Смещение датчика относительно осей тангажа и крена

 

0,000872

 

 

 

 

 

 

(0,05)

Случайные помехи по осям тангажа и крена

 

 

0,00174

 

 

 

 

 

 

(0,1)

Атмосферные помехи орбитальной частоты по осям тангажа

0,000523

и крена (максимальная амплитуда)

 

 

 

(0,03)

Упругие

тепловые деформации конструкции по

осям

(мак­

 

симальная амплитуда):

 

 

 

 

 

тангажа

 

 

 

 

0,00437

 

 

 

 

 

 

(0,25)

крена

 

 

 

 

 

0,0

 

 

 

 

 

 

(0,0)

рыскания

 

 

 

 

0,000698

 

 

 

 

 

 

(0,04)

106


В табл. 3. 1 приведены сравнительные данные по точности определения углов у, •&, тф при исходных ошибках системы ори­ ентации (табл. 3.2—3.4). Расчетные данные получены для слу­ чая движения КА по околоземной орбите с высотой перигея и апогея соответственно 140 и 370 км, наклонения — я/2 рад (—90°), аргумента перигея орбиты — 0 рад (0°), начальной истинной аномалии —1,393 рад (—80°) и угла наклонения на Солнце в геоцентрической системе координат —0,196 рад

( - 10°) [8].

Т а б л и ц а

3. 3

 

 

Основные ошибки при ориентации на Солнце

 

 

 

И с т о ч н и к

В е л и ч и н а

Смещение солнечного датчика (по всем осям)

0,000174 рад

 

 

 

(0,01°)

Случайные помехи

 

0,000174 рад

 

 

 

(0,01°)

Ошибки эфемериды:

 

 

в направлении

полета

1500 м

перпендикулярно

направлению полета

150 м

ошибка высоты

 

150 м

Т а б л и ц а

3. 4

 

 

Ошибки датчика угловой скорости

 

 

 

И с т о ч н и к

В е л и ч и н а

Смещение датчика

(по всем осям)

0,000872 рад

 

 

 

(0,05°)

Случайные помехи

 

0,00872 рад/ч

 

 

 

(0,5 град/ч)

Уходы нуля (по всем осям)

0,00872 рад/ч

 

 

 

(0,5 град/ч)

Максимальная амплитуда упругих тепловых

дефор­

маций конструкции по осям

 

тангажа

 

0,00437 рад

 

 

 

(0,25°)

крена

 

 

0,0 рад

 

 

 

(0,0°)

рыскания

 

0,000698 рад

 

 

 

(0,04°)

107


Астроориентация КА

В последние годы широкое распространение получает ориен­ тация КА, основанная на измерениях опорных светил (звезд). При этом осуществляется автоматическая обработка данных при помощи БЦВМ. Источником информации (измерений) может служить, например, автоматическое оптико-электронное устрой­ ство ориентации, состоящее из объектива, оптического модуля­ тора, фотоэлектронного умножителя и электронного преобразую­ щего блока [11]. При решении задачи ориентации при помощи БЦВМ по данным трехосного оптико-электронного устройства наибольшую сложность представляет распознавание звездного узора при выборе опорных звезд. В общем случае звезды можно распознать по трем характеристикам: по их относительному положению на небесной сфере, по их яркости и наблюдаемому: цвету свечения. ■

Наиболее простой из них — первый (по распознаванию звезд­ ного узора). При этом возникает противоречие при построении системы ориентации: чем меньше опорных звезд, тем больше; требования к оптической системе, одновременно меньше требо­ вания к памяти БЦВМ, хранящей в ЗУ угловые расстояния между выбранными для ориентации опорными светилами.

На рис. 3. 34 показана зависимость между возможными еди­ ничными ошибками в измерении углового расстояния и числом

Допустимая ошибка измерения углового расстояния,

р а д град

Рис. 3.34. График зависимости долу* стимой ошибки измерения углового расстояния звезд от числа звезд, уча­ ствующих в процессе распознавания

108

7

l l l l l l l i l l l

lllll

Угол поворот а сканирую ­ щ е го уст рой ст ва ; в

а)

I I I

N 11

0 i Qj

Он

5)

Рис. 3. 35. Операции распознавания опорных звезд:

б— угол поворота сканирую щ его устройства

звезд, угловое расстояние между которыми измеряется [11]. Вели­ чина допустимой единичной ошибки возрастает пропорционально числу звезд, участвующих в процессе распознавания звездного узора. Выделенная область соответствует испытаниям оптической системы с углом зрения 0,802 рад (46°) и вероятности неодно­ значного распознавания рисунка звезд, равной 0,99.

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

1.Накопление данных угловых измерений. В запоминающее устройство БЦВМ вводится информация, состоящая из импуль­ сов, соответствующих моментам появления звезд, а также фоно­ вым шумам (рис. 3. 35, а).

2.Селекция при помощи сравнивающего алгоритма. При ска­ нировании небосвода оптико-электронным устройством в момент попадания в приемное устройство излучения звезд генери­

руются импульсы,

разделенные

двумя

равными

интервалами

А ц и Ajk (см.

рис. 3. 35, б)

Затем

БЦВМ

ведет поиск

по каталогу угловых измерений, отбирая такие результаты, кото­ рые обладают свойствами Д;у- = Д;-*-[-Ё^, где Sg <с Аи , Дjk.

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

3. Определение положения цели. Для каждого набора из трех углов БЦВМ определяет номерные координаты звезды — цели.

4. Вычисление угловых расстояний между звездами. Машина рассчитывает и хранит в запоминающем устройстве расстояния

109