- •2. Относительное движение
- •2.1. Уравнения движения двух летательных аппаратов.
- •2.2. Уравнения движения в орбитальной системе координат.
- •2.3. Решение и приведение дифференциальных уравнений относительного движения двух ла.
- •2.4. Относительное движение двух летательных аппаратов в других системах координат.
- •2.5. Траектория относительно опорного движения. Смещение движения.
- •2.6. Вектор эксцентриситета.
- •2.7. Уравнения прогноза относительного движения двух ла, каждый из которых движется по эллиптической орбите.
2.6. Вектор эксцентриситета.
Рассмотрим более подробно величину е, которая по компонентам вектора состояния в моментопределяется формулой
Если мы спрогнозируем вектор состояния на момент и по компонентам нового вектора состояния вновь определиме, то получим
Существенным является то, что величина е для относительного движения, заданного начальным вектором состояния, постоянна, хотя ее компоненты меняются во времени:
т.е. можно трактовать как компоненты векторае, модуль которого постоянен, а направление относительно осей Ох и Оу характеризуется переменным углом отсчитываемым в соответствии с рис. 2.2.
Вектор е вращается в ОСК с угловой скоростью изменения угла , т.е. со скоростью. Заметим, что относительно инерциального пространства орбитальная СК (или цилиндрическая СК) вращается также с угловой скоростью, но в направлении, обратном направлению вращения векторае. Следовательно, можно утверждать, что направление вектора е неизменно в инерциальном пространстве и изменение его проекций в ОСК происходит только за счет вращения самой ОСК.
Для того чтобы понять, как расположен вектор е относительно кеплеровой возмущенной орбиты в инерциальном пространстве, заметим, что при достигается перигей,а привозмущенной орбиты. Но приимеем,т.е. – векторе в момент прохождения перигея направлен противоположно оси Оу (рис. 2.5). Как следует из рисунка, угол фактически равен фазовому углу, который пройден ЛА по возмущенной орбите от момента последнего прохождения перигея, т.е. он равен истинной аномалии ЛА
Рис. 2.4. Вектор эксцентриситета в декартовой объектоцентрической системе координат.
Рис. 2.5. Положение вектора эксцентриситета относительно кеплеровой возмущенной орбиты.
Вектор е имеет размерность длины и в соответствии с ранее полученными результатами его компонента характеризует высотный сдвиг возмущенного движения относительно круговой орбиты, радиус которой равен большой полуоси возмущенной орбиты
Выражения для перигейного и апогейного расстояний возмущенной орбиты могут быть теперь записаны в виде
Из кеплеровой теории известны соотношения
Сравнение выражений для , полученных при рассмотрении относительного движения и в кеплеровой теории, сразу определяет зависимость между эксцентриситетом относительного движенияе и кеплеровым эксцентриситетом :
2.7. Уравнения прогноза относительного движения двух ла, каждый из которых движется по эллиптической орбите.
До сих пор, говоря об относительном движении двух ЛА, мы полагали, что орбита одного из них является круговой. Можно ли использовать полученные результаты для описания относительного движения двух ЛА, каждый из которых движется по эллиптической орбите? Оказывается можно, и эта возможность предоставляется нам благодаря тому, что мы рассмотрели относительное движение в линейном приближении. Действительно, пусть два ЛА движутся по эллиптическим орбитам и пусть для обоих ЛА можно указать такое фиктивное опорное круговое движение и такую связанную с ним орбитальную систему координат, что будут выполняться все необходимые условия для решения задачи в линейном приближении, т.е. соответствующие отклонения координат обоих ЛА от начала системы отсчета будут достаточно малы (это одновременно означает и малость рассогласований координат самих ЛА). Для каждого ЛА мы вправе записать уравнение прогноза:
где - матрица прогноза, элементы которой являются функциями времени прогноза и угловой скорости опорного движения;- векторы состояния относительных движений летательных аппаратов в начальный момент времени;— векторы состояния относительных движений летательных аппаратов в конечный момент времени. Из разности этих уравнений нетрудно получить
где — векторы состояния движения первого ЛА относительно второго соответственно в начальный и конечный моменты времени.
Естественно, что при указанном относительном движении невозможно установить простую связь между параметрами относительного движения и кеплеровыми элементами, поскольку ие в этом случае характеризуют возмущение не круговой, а эллиптической орбиты.
Завершая описание линейной теории относительного движения, остановимся на свойствах матрицы прогноза , уравнения (2.18). Запишем это уравнение в виде
Где векторы состояния относительного движения соответственно в начальный моменти конечный момент— матрица, осуществляющая прогноз вектора состояния на временной интервал
Получим вектор состояния в момент иным образом, осуществив сначала прогноз на промежуточный моменти продолжив его затем до момента:
Таким образом, результирующая матрица прогноза на временной интервалесть результат произведения матриц прогнозапри осуществлении прогноза последовательно на интервали затем на интервалЭтот вывод можно обобщить следующим образом:
До сих пор рассматривался прогноз вектора состояния относительного движения на некоторый момент в будущем, но, очевидно, прогноз может быть осуществлен и на заданный момент в прошлом — для этого достаточно определить элементы матрицы прогноза для отрицательного временного интервала. При этом
При решении ряда задач может оказаться более удобным рассматривать относительное движение в какой-либо системе координат, отличающейся от использованных нами орбитальных вращающихся систем (ОСК или ЦСК). Если преобразование вектора состояния R при переходе к новой системе координат является линейным и задается постоянной матрицей А, то легко получить матрицу прогноза относительного движения в новой системе координат:
В качестве примера подобного преобразования укажем на преобразование "замораживания" вектора состояния, заданного в ОСК. В результате такого преобразования мы приходим к ОСКЗ (ОСК "замороженная"), оси которой в любой момент времени совпадают с осями ОСК и вместе с тем в каждый текущий момент неподвижны относительно инерциального пространства. Матрицу А этого преобразования можно представить в блочном виде
а результат преобразования вектора состояния R, заданного в ОСК, — в виде
Исходя из характера преобразования, можно сказать, что вектор состояния r есть результат проектирования разности абсолютных векторов состояния возмущенного и опорного движений на оси ОСК.
Матрица прогноза полученная по формуледает следующий аналог решения (2.6) для ОСКЗ: