Тема: Основы теории горизонтальных астроориентаторов

  • Вид работы:
    Контрольная работа
  • Предмет:
    Авиация и космонавтика
  • Язык:
    Русский
  • Формат файла:
    MS Word
  • Размер файла:
    91,48 Кб
Основы теории горизонтальных астроориентаторов
Основы теории горизонтальных астроориентаторов
Вы можете узнать стоимость помощи в написании студенческой работы.
Помощь в написании работы, которую точно примут!














Контрольная работа

ОСНОВЫ ТЕОРИИ ГОРИЗОНТАЛЬНЫХ АСТРООРИЕНТАТОРОВ

Содержание

1. Принципы получения информации, необходимой для вычисления координат и курса

. Алгоритмы определения координат и курса по информации о высотах двух звезд

. Анализ погрешностей астроориентатора

Литература

Приложение

1. Принципы получения информации, необходимой для вычисления координат и курса

Горизонтальными астроориентаторами (АО) или астроориентаторами горизонтальной системы координат называют устройства, в которых вычисление координат точки места и курса ЛА осуществляется путем использования информации об измеренных горизонтальных координатах двух светил или величин, функционально связанных с этими координатами.

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

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

В этом случае с датчиков углов, установленных на осях подвеса астросекстантов, можно получить информацию об углах a1; b1 и a2; b2, характеризующих положение оптических осей телескопов, относительно трехгранника x1; y1; z1 с ортами i1, j1, k1, связанного с платформой. Если платформа идеально строит плоскость горизонта и оптические оси телескопов совпадают с линиями визирования светил, т.е.

j1 = r ; q = S11; q22 = S22 , (1)

где r - орт вертикали, то при этом можно полагать, что

b1 = h1; b2 + h2 , (2)

где h1 и h2 - высоты звезд.

В действительности платформа реализует плоскость горизонта, с некоторыми погрешностями, также с погрешностями происходит и процесс слежения за светилами. В силу этого b1 = H1; b2 = H2, где H1 и H2 - измеренные значения углов высоты светил.

В том случае, когда астросекстанты установлены на платформе, стабилизированной в азимуте, например так, что ось i1 ориентирована на Север, то при выполнении условий (1) можно считать a1 = AC1; a2 = АС2, где АС1 и АС2 - азимуты первой и второй звезд.

В общем же случае углы a1 и a2 есть курсовые углы светил, характеризующие поворот астросекстантов относительно одной из осей платформенного трехгранника. Используя сигналы с датчиков углов b1,2, можно вычислить функции sin H1 и sin H2, определяющие величины направляющих косинусов между ортами оптических осей телескопов и ортом j1 = r1 платформы,

j1 · q11 = r1· q11 = sin H1; j1 · q22 = r1· q22 = sin H2. (3)

Второй подход предусматривает установку двух астросекстантов (или одного, поочередно пеленгующего две звезды) на некоторой базовой плате, смонтированной на корпусе летательного аппарата. В этом случае могут быть измерены углы b1, a1 и b2, которые теперь характеризуют положение оптических осей телескопов относительно трехгранника, жестко связанного с установочной базой. Если принять, что трехгранник осей базы совпадает с осями трехгранника, связанного с ЛА, то углы b1,2 и a1,2 определяют положение оптических осей телескопов в осях связанного трехгранника. Информацию об углах, характеризующих положение единичного орта вертикали места в осях трехгранника, связанного с JlA, можно получить с датчиков углов, установленных на осях подвеса построителя вертикали (гировертикали, инерциальной вертикали и т.п.).

Таким образом, в этом случае в осях связанного с ЛА базиса iC , jC , kC могут быть определены три единичных вектора

q11 = iC cos b1 cos a1 + jC sin b1 + kC cos b1 sin a1;22 = iC cos m 2 cos n2 + jC sin m2 + kC cos m2 sin n2a (4) 1 = iC r11 + jC r12 + kC r13 ,

где r 1 - орт построенного направления вертикали.

Проекции вектора r1 на оси связанного с ЛА трехгранника можно представить так:

r11 = - sin J2 cos g2; r12 = cos J2 cos g2;

r13 = sin g2, (5)

где J2 и g2 - углы поворота связанного базиса относительно трехгранника построителя вертикали.

Информация об этих углах поступает в виде сигналов с датчиков, установленных на осях подвеса построителя.

Очевидно, что после того, как вычислены направляющие косинусы ортов q11, q 22 и r1 можно вычислить следующие скалярные произведения:

r × q11 = r11 cos b1 cos a1 + r12 sin b1 + r13 cos b1 sin a1 ;× q12 = r11 cos b2 cos a2 + r12 sin b2 + r13 cos b2 sin a2 .. (6)

Обозначим r × q11 = sin H1; r × q12 = sin H2 и запишем (6) в виде

sin H1 = r11 cos b1 cos a1 + r 12 sin b1 + r 13 cos b1 cos a1 ;H2 = r11 cos b2 cos a2 + r12 sin b2 + r13 cos b2 sin a2 . (7)

B (7) углы H1 и Я2 можно рассматривать как вычисленные значения углов высоты первой и второй звезд.

Из сравнения выражений (3) и (7) следует вывод, что, в конечном счете, вне зависимости от того, как установлены астросекстанты: на платформе или «на корпусе» объекта, можно получить информацию о величинах направляющих косинусов между ортами оптических осей телескопов и ортом построенного направления вертикали.

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

При размещении секстантов «на корпусе» не возникает дополнительных трудностей, связанных с конструированием построителя вертикали или курсовертикали. Увеличение числа вычислительных операций, имеющее место в этом случае, при наличии БЦВМ не является принципиальным. Однако к следящим системам астросекстантов потребуется предъявить более жесткие требования по быстродействию. Кроме того, в этом случае появляются дополни- тельные источники погрешностей, обусловленные тем, что блок астросекстантов и построитель вертикали (курсовертикаль) устанавливаются в различных точках летательного аппарата. Если объект является абсолютно жестким телом, то дополнительные погрешности могут быть вызваны несовпадением базисов, относительно которых измеряются углы y2, J2, g2 и b1,2, a1,2 .

В случае «нежесткого» объекта появляются погрешности, обусловленные тем фактом, что измерение углов y2, J2, g2 и b1,2, a1,2 осуществляется относительно базисов, изменяющих взаимную ориентацию, например, за счет изгибных колебаний летательного аппарата.

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

. Алгоритмы определения координат и курса по информации о высотах двух звезд

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

Поскольку на подвижных объектах определение азимута требует наличия точного измерителя курса, то в литературе [3] рассмотрены алгоритмы, предусматривающие вычисление углов азимута светил. Измерение разности азимутов двух звезд осуществить на подвижном объекте проще, чем измерение самих углов, поскольку для определения этой величины следует измерить только курсовые углы двух звезд.

Однако метод высот двух светил позволяет вычислить координаты точки места вообще без привлечения информации о каких-либо азимутальных координатах светил. Для использования этого метода достаточно располагать значениями углов высот двух звезд, полученных путем измерений, либо значениями sin h1 и sin h2.

Обратимся к рассмотрению алгоритмов определения координат точки места ЛА в этом случае. Воспользуемся первым выражением (20) и применительно к двум звездам запишем следующую совокупность равенств:

sin h1 = sin d1 sin j + cos d1 cos j cos (S - a1);

sin h2 = sin d2 sin j + cos d2 cos j cos (S - a2). (8)

Если полагать, что путем измерений определены значения функций sin h1 и sin h2 для пеленгуемой определенной пары звезд и в силу этого известны их экваториальные координаты d1,2 и a1,2, то совокупность равенств (8) можно рассматривать как систему трансцендентных алгебраических уравнений, относительно переменных j и S.

Решая при помощи выбранного итерационного метода в бортовом вычислителе систему (8), можно определить координаты точки, в которой измерены высоты звезд.

Для того чтобы система уравнений вида (8) имела решения, необходимо, чтобы входящие в нее уравнения были линейно независимы. Условие независимости уравнений (8) сводится к выполнению требования


Находя соответствующие частные производные и раскрывая определитель, получаем следующее выражение для якобиана:

W = [cos d1 sin d2 sin (S - a1) - sin d1 cos d2 sin (S - a2) ] cos2 j +

+ (1/2) cos d1 cos d2 sin (a2 - a1) sin 2j

Анализ полученного выражения позволяет выявить некоторые особенности системы уравнений (8). Так, очевидно, что при j = p/2 якобиан равен нулю и уравнения системы несовместны.

Этот математический результат отражает тот факт, что при j = p/2 точка места совпадает с Северным полюсом и понятие звездного времени, а также и долготы места теряет смысл.

Рассмотрим, какой вид принимает якобиан в некоторых случаях полета.

Якобиан системы уравнений (8)равен нулю, если полет происходит в плоскости, определяемой векторами S11; S22, либо в некоторый момент времени траектория полета ЛА пересекает эту плоскость. При этом имеет место равенство азимутов звезд: A1 = A2.

Таким образом, помимо особой точки j = p/2 в пространстве имеется и особая плоскость (плоскость, проходящая через векторы S11; S22), в которой нарушается совместность уравнений (8). Этот результат объясняется тем, что при нахождении ЛА в особой плоскости имеет место равенство h2 = h1 + r, где r есть постоянный угол между векторами S11; S22 .

При полете в этой плоскости, в принципе, достаточно измерять высоту одной из звезд, так как измерение высоты второй звезды не дает новой информации. По измерению, например, h1 однозначно определяется положение ЛА в плоскости S11; S22, а положение последней в пространстве (относительно базиса r1, r2, r3) всегда известно.

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

Полагаем, что углы склонений d1,2 и прямых восхождений a1,2 выбранной для пеленгации пары звезд известны с приемлемой точностью. Примем также, что в результате измерений определены

sin H1 = sin h1; sin H2 = sin h2

Кроме того, будем считать, что известны приближенные значения координат точки места jп и Sп = T + lп.

Знание приближенных значений координат необходимо для того, чтобы процесс итерационного решения можно было бы осуществить и, кроме того, для исключения неоднозначности решения системы (8), поскольку, как нетрудно видеть, эта система может иметь две пары решений: j; S и p - j; p - S.

Информацию о величинах jп и lп можно получить от астрономического вычислителя координат или от инерциальной системы. Располагая величинами д1)2, d1,2, a1,2, jп и Sп, в бортовой ЦВМ можно вычислить

sin hB = sin d1 sin jп + cos d1 cos jп cos (Sп - a1);

sin h2B = sin d2 sin jп + cos d2 cos jп cos (Sп - a2)

Представим теперь функции sin h1 и sin h2 в виде разложения в ряд в окрестностях точки с координатами jп, Sп. Ограничиваясь тремя членами ряда, имеем


д sin h1/ дj = sin d1 cos j - cos d1 sin j cos (S - a1);

д sin h2/ дj = sin d2 cos j - cos d2 sin j cos (S - a2);

д sin h1/ дS = - cos d1 cos j sin (S - a1); (15)

д sin h2/ дS = - cos d2 cos j sin (S - a2);

и их значения в точке с координатами jп, Sп всегда могут быть вычислены. Положим, что в вычислителе сформированы следующие разности:

D1 = sin h - sin h1, D2 = sin h - sin h2. (16)


Выражения (17) представляют собой систему линейных алгебраических уравнений относительно неизвестных Dj и DS, которые имеют смысл итерационных поправок к приближенным значениям координат.

Уравнения (17) разрешимы, если их якобиан:


не равен нулю. Нетрудно видеть, что свойства якобиана системы (18) аналогичны свойствам якобиана системы (8) и уравнения (17) разрешимы во всех, за исключением отмеченных ранее, случаях.

Решения системы (17) имеют вид

;


Используя Dj и DS, определенные по (19), можно найти уточненные значения координат

j1п = jп + Dj; S1п = S п + DS

Получив по (20) значения координат, итерационную процедуру можно повторить и в результате найти новые значения поправок и еще раз уточнить координаты, полученные после первого итерационного цикла.

Итерационный процесс продолжается до тех пор, пока не будут выполнены условия

|D1| £ a ; |D2| £ a ,

где a - число, задаваемое из соображений требуемой точности вычисления координат.

Таким образом, на каждом итерационном такте в ЦВМ выполняются следующие операции: вычисляются значения sin h1B и sin h; находятся D1 и D2; определяются значения соответствующих частных производных; вычисляется по формуле значение якобиана системы; определяются по формулам (19) значения Dj и DS; находятся уточненные значения координат места по формулам (20).

Рассмотренная итерационная процедура соответствует использованию при решении системы уравнений (8) так называемого метода Ньютона. Эта процедура несколько упростится, если значения соответствующих частных производных и якобиана системы вычислить лишь один раз на первом итерационном цикле, используя для этой цели значения jп и Sп. Тогда в последующих итерационных циклах поправки Dj и DS будут вычисляться проще, в частности, в выражениях (19) не потребуется выполнять операции деления. Изложенное изменение процедуры соответствует построению решений системы (8) по так называемому модифицированному методу Ньютона.

. Анализ погрешностей астроориентатора

Перейдем к анализу погрешностей рассматриваемого метода определения координат. Если полагать, что значения функций sin h1 и sin h2 определены точно, то путем решения системы уравнений (8), в принципе, можно определить координаты j и S с любой, наперед заданной, точностью. Если же путем измерений получены значения функций sin H1 и sin H2, то в БЦВМ будут найдены координаты Ф и L, являющиеся решениями системы

sin H1 = sin d1 sin Ф + cos d1 cos Ф cos (L - a1);H2 = sin d2 sin Ф + cos d2 cos Ф cos (L - a2)

при условии, что значения углов склонений и прямых восхождений пеленгуемых звезд известны с приемлемой точностью.

Вычитая из системы уравнений (21) уравнения (8), находим:

D1 = sin H1 - sin h1 = (sin Ф - sin j) sin d1 + + [cos Ф cos (L - a1) - cos j cos (S - a1) ] cos d1;

D2 = sin H2 - sin h2 = (sin Ф - sin j) sin d2 + (22)

+ [cos Ф cos (L - a2) - cos j cos (S - a2) ] cos d2 .

Представим найденные из (21) координаты в виде Ф = j + dj; L = S + dS, где dj и dS - погрешности в вычислении координат.

Из (22) с точностью до величин второго порядка малости имеем

D1 = [sin d1 cos j - cos d1 sin j cos (S - a1) ] dj -

- cos d1 cos j sin (S - a1) dS;

D2 = [sin d2 cos j - cos d2 sin j cos (S - a2) ] dj -

- cos d2 cos j sin (S - a2) dS . (23)

Выражения (23) можно рассматривать как систему уравнений относительно переменных dj и dS. Разрешая эту систему, получаем

, (24)

где W определяется выражением (10).

Совокупность формул (24) позволяет на уровне первого приближения найти погрешности в определении углов j и S, вызванные погрешностями в определении величин sin h1 и sin h2.

Получим формулы для определения величин D1 и D2. При этом рассмотрим влияние на погрешности в определении значений sin h1 и sin h2 таких факторов, как неточное построение вертикали, наличие погрешностей фотоследящих систем и датчиков углов, установленных на осях подвесов телескопов. Положим первоначально, что фотоследящие системы обеспечивают точное слежение и что датчики углов a1;2 и b1;2 идеальны, но вертикаль построена с погрешностями.

Будем считать, что отклонение платформы построителя от плоскости горизонта определяется углами J1 и g1, а для ортов горизонтного платформенного базиса сохраним обозначения i10, j10, k10. Матрица направляющих косинусов между базисом i1, j1, k1 и построенным платформой базисом i10, j10, k10. доопределяется выражением

(25)


Поскольку углы J1 и g1 характеризуют отклонение плоскости платформы от горизонта, за счет поворота по осям внутренней и наружной рамок подвеса, их можно принять малыми и считать

; . (26)

Используя (25) и (26), находим:

J = j10 - i10J1 + k10 g1, cosJ . (27)

В осях трехгранника i10, j10, k10 при условии, что светила запеленгованы, векторы q11 и q22 представим так:

q 11 = S11 = i10 cos h1 cos a + j10 sin h1 + k10 cos h1 sin a; 22 = S22 = i10 cos h2 cos a + j10 sin h2 + k10 cos h2 sin a,

где a, a - курсовые углы звезд в осях платформы построителя при его горизонтальном положении.

Теперь, используя (26) и (27), найдем следующие скалярные произведения:

j1 · q 11 = sin h1 - J1 cos h1 cos a+ g1 cos J cos h1 sin a;1 · q 22 = sin h2 - J1 cos h2 cos a+ g1 cos J cos h2 sin a.

В осях трехгранника i1, j1, k1 имеем

j1 · q 11 = sin b1; j1 · q 22 = sin b2 . (30)

D11 = sin b1 - sin h1 = - J1 cos h1 cos a + g1 cos J cos h1 sin a(31)21 = sin b2 - sin h2 = - J1 cos h2 cos a + g1 cos J cos h2 sin a.

Выражения (31) позволяют найти величины D1 и D2, обусловленные погрешностями построителя вертикали.

Получим теперь выражения для величин D12 и D22, характеризующих погрешности в определении функций sin h1 и sin h2, вызванные погрешностями фотоследящих систем (ФСС) и датчиков углов. Качественно эта группа погрешностей приводит к тому, что либо от оптической оси соответствующего телескопа не совпадает с линией направления на светило, либо при совпадении этих ортов, сдатчиков углов снимаются сигналы об углах a и b, не соответствующих действительному положению оптической оси.

Представим орты оптических осей телескопов так:

q11 = S11 + Dq11; q22 = S22 + Dq22 , (32)

где векторы q11 и q22 характеризуют погрешности, вызванные погрешностями фотоследящих систем и датчиков углов. Умножая (32) скалярно на /10, находим

j1 · q11 = sin h1 + j10 · Dq11 ;

j1 · q22 = sin h2 + j10 · Dq22 .

В свою очередь, проекции векторов Dg11 и Dg11 на ось j10 можно представить в виде

j10 · D q11 = Db1 cos h1, j10 · D q22 = Db2 cos h1 ,

где Db1,2 - суммарные углы, обусловленные как погрешностями ФСС, так и погрешностями датчиков углов b1,2;

Db1 = Db11 + Db12; Db2 = Db21 + Db22,

где Db11; Db21 - погрешности соответствующих фотоследящих систем; Db12; Db22 - погрешности датчиков углов.

Запишем теперь выражения (33) в виде

j10 · q22 = sin h2 + (Db21 + Db22) cos h2.

Учитывая, что в алгоритме вычисления координат используется информация о функциях j10 · q11 = sin H1; j10 · q22 = sin H2, представим (34) так:

D1,2 = sin H1 - sin h1 = (Db11 + Db12) cos h1; (35)

D2,2 = sin H2 - sin h2 = (Db21 + Db22) cos h2 .

Выражения (35) позволяют найти величины D12 и D22. Используя полученные формулы, по (24) можно оценить влияние погрешностей построителя вертикали, ФСС и датчиков углов на погрешности определения координат или оценить совокупное влияние этих факторов.

Отметим одно важное обстоятельство, которое следует принимать во внимание при оценке погрешностей по выражениям (24). Поскольку входящие в (24) частные производные и якобиан W зависят от района полета JlA, то и численные значения погрешностей определения координат при постоянных D1,2 будут зависеть от координат точки места. В частности, приняв, что D1 = D2 = D, из (24) получаем

. (36)

Подставляя в последние равенства выражения для частных производных и якобиана из (15) и (11) находим, что

(37)

(38)

Из (38), в частности, следует, что в особой точке системы уравнений (8) при j = p/2 выражение для погрешности dS также имеет точку разрыва. Этот результат отражает тот факт, что в окрестностях особой точки уравнений (8) определение координаты S будет выполняться с существенными погрешностями, а в самой особой·точке определение этой координаты вообще теряет смысл.

Таким образом, можно заключить, что рассматриваемый метод определения координат имеет ограничения на применение, связанные с районами полета JlA. При полетах в окрестностях географических полюсов следует либо отказываться от определения координаты S, либо переходить на определение координат точки места в другой координатной сетке, например ортодромической.

Помимо рассмотренных погрешностей определения координат, вызванных неточным определением величин sin h1 и sin h2, на точность их определения повлияют погрешности выбранного численного метода решения системы уравнений (8) и погрешности, вызванные представлением чисел в БЦВМ в разрядной сетке конечной длины. Аналитическое рассмотрение влияния этой группы погрешностей на погрешность определения координат, как правило, выполнить не удается. Оценить влияние указанных факторов можно путем моделирования на универсальной ЭЦВМ. При этом составляется программа решения системы уравнений (8) тем методом, который предполагается использовать в бортовой ЦВМ и решение задачи выполняется при представлении в ЭЦВМ чисел в разрядной сетке различной длины. Влияние такого фактора, как представление в БЦВМ чисел с фиксированной запятой, также может быть учтено программой.

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

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

Рассмотрим возможные алгоритмы определения курса летательного аппарата. Задача в этом случае сводится к вычислению азимута одной из пеленгуемых звезд и к последующему определению курса с использованием сигнала датчика курсового угла этой звезды. Причем, если в случае установки астросекстантов на построителе вертикали сигнал о курсовом угле звезды, определенном в осях платформы построителя, может быть получен непосредственно с соответствующего датчика, то при установке секстантов «на корпусе» JlA сигналы датчиков углов a1,2 должны быть подвергнуты обработке по определенному алгоритму с целью получения информации о курсовом угле. Таким образом, алгоритм определения курса в общем случае содержит две основные части: алгоритм вычисления азимута и алгоритм вычисления курсового угла.

При вычислении азимута может быть использован алгоритм, рассмотренный в разд. 3 применительно к горизонтальному астрокомпасу

(39)

Поскольку углы j и S вычислены при решении задачи определения координат, а значение cos h можно вычислить по формуле , то по (39) достаточно определить знаки функций cos АC и sin AC и по сочетанию этих знаков отнести угол AC к определенной четверти. Величина же угла АC в этой четверти определяется путем вычисления функции

, (40)

при значениях аргумента от 0 до p/4 либо функции

, (41)

при аргументе, изменяющемся от p/4 до p/2.

Анализ погрешностей определения угла азимута по алгоритмам (40) и (41) может быть выполнен с использованием формул (81), (82). При этом величины dj dS определяются из выражений (24).

Алгоритм вычисления курсового угла светила при установке астросекстантов на корпусе летательного аппарата строится следующим образом. В осях базиса iс; jс; kc единичный вектор оптической оси секстанта q11 представляется так:

q11 = ic cos a1 sin b1 + ic sin b1 + kc sin a1 cos b1.

В осях платформенного базиса i1; j1; k1 тот же вектор можно представить в виде

q11 = i1 cos a cos H1 + j1 sin H1 + k1 sin a cos Н1 ,

где a - курсовой угол светила в осях платформы построителя вертикали.

Матрицу направляющих косинусов между связанным и платформенным базисом представим так:

,

где матрицы [J2] и [g2] имеют вид

; ,

J2, g2 - углы, информация о которых снимается с датчиков углов платформы.

Имеет место следующее матричное равенство:

(42)

Используя (42), получаем

. (43)


Совокупность равенств (43) либо выражение (44) позволяют вычислить курсовой угол светила в осях платформы и могут рассматриваться как алгоритм его вычисления.

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

Для случая, когда построитель вертикали идеален, (42) можно представить в виде

, (45)

где h1 и aк - значения угла высоты и курсового угла светила, определяемые при отсутствии погрешностей построителя вертикали.

Равенство (42) запишем так:

. (46)

Вычитая из (46) выражение (45), получаем

. (47)

В свою очередь,

D q11 = cos H1 cos a1 - cos h1cos aк;

D q12 = sin H1 - sin h1 ; (48)

D q13 = cos H1 sin a1 - cos sin aк .

Теперь из (48) нетрудно получить следующие, справедливые с точностью до величин второго порядка малости относительно D q13, выражения для оценки погрешностей.

sin a1 - sin aк = (D q13 - D q12 sin h1 sin aк) sec h1;

cos a1 - cos aк = (D q11 + D q12 sin h1 cos aк) sec h1. (49)

В формулах (49) величины Dq определяются по (47) и являются функциями погрешностей построителя вертикали.

Курс летательного аппарата по вычисленным значениям углов азимута и курсового угла светила определяется так:

y = АC + aк ,

если используется построитель вертикали, платформа которого не имеет азимутальной оси. Если используется курсовертикаль, например, со свободной в азимуте ориентацией платформы.

y = АC + aк + y2 , (50)

где y2 - угол поворота платформы относительно корпуса в азимутальной плоскости, информация о котором поступает с датчика угла, установленного на вертикальной оси подвеса платформы.

В том случае, когда в курсовертикали обеспечивается ориентация осей платформы по странам света, должно выполняться равенство

звезда астроориентатор алгоритм координата

АC + aк = y2 . (51)

Если (51) не выполняется, то этот факт можно рассматривать как свидетельство наличия азимутального дрейфа платформы курсовертикали и, сформировав сигнал D y2 = y2 - (AC + aк), скорректировать ее азимутальную ориентацию.

Литература

1. Коваленко В.В., Лысов А.Н. Малогабаритная инерциальная система. Учебное пособие. 2010.

. Журнал «Новости электроники» КОМПЭЛ., выпуск № 4., 2011 // #"justify">Приложение


Рис. 2. Положение особой плоскости

Похожие работы

 

Не нашел материала для курсовой или диплома?
Пишем качественные работы
Без плагиата!