Глава 3. УЛУЧШЕНИЕ ПЕРВОНАЧАЛЬНОЙ ОРБИТЫ

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

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

Наряду с учебными пособиями [1], [2] для более подробного ознакомления со способами улучшения орбит можно рекомендовать [7], [8].

§ 3.01. Дифференциальное исправление орбит. Постановка задачи

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

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

же моменты по элементам предварительной орбиты, отнесенной к моменту

Обозначим через ошибки вычисленных координат, так что

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

Зависимость от выражается с точностью до членов первого порядка формулами

где — вычисленные значения (индекс для сокращения записи опускаем).

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

где — элементы орбиты, а — значения частных производных по на момент Членами порядка выше первого в формулах (3.3.01), (3.3.02) относительно пренебрегают. В формулах (3.3.02) часто для

удобства вычислений в качестве принимают некоторые функции элементов.

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

После нахождения численных значений производных подставляют выражения (3.3.02) для в (3.3.01) и получают соотношения вида

где — численные коэффициенты. Эти соотношения представляют собой так называемые условные уравнения относительно решаемые по способу наименьших квадратов (с применением ЭВМ, если число наблюдений достаточно велико). См. ч. VII, гл. 4.