7.5.4. Прямолинейные или почти прямолинейные орбиты
В принципе описанным выше методом можно пользоваться при любых значениях эксцентриситета и наклонения. Однако, когда орбита прямолинейная или почти прямолинейная, эффективность метода падает. Объясняется это тем, что в ряде используемых соотношений, и в частности в уравнении (7.18), величина h стоит
знаменателе. Поэтому в случае, когда все три компоненты
становятся меньше некоторого предела, следует, чтобы не терять точность, перейти к другим переменным. В качестве таких новых переменных можно взять полярные координаты
, для которых справедливы следующие соотношения:
Величины
— это эклиптические долгота и широта (или прямое восхождение и склонение) соответственно.
Кроме того, имеем
(7.55)
Таким образом, если даны
, то при помощи (7.54) и (7.55) можно вычислить
С другой стороны, имеем
Заметим, что знак
совпадает со знаком
, а соотношения
позволяют выбрать верное значение
Кроме того,
и, следовательно, по
можно найти
.
Продифференцировав первое уравнение (7.57) и воспользовавшись (7.3), получаем
где
Надо сказать, что
равно нулю только в случае прямолинейной орбиты.
Снова обозначая индексами k величины, относящиеся к невозмущенной кеплеровской орбите, можно написать
Возмущение величины радиуса-вектора равно
так что
Дважды дифференцируя это соотношение, получаем
и, следовательно,
Дифференцируя второе и третье соотношения (7.57) и обозначая возмущения величин а и Р индексом
, т. е.
получаем
Здесь индекс k указывает на то, что величины, заключенные в в скобки, соответствуют невозмущенному кеплеровскому движению.
Предположим, что опорная орбита все время является прямолинейной. Тогда в отличие от предыдущего случая (
) ее нельзя получить из оскулирующей орбиты, если только истинная орбита сама не является в точности прямолинейной. Если орбита прямолинейная, то
и
так что уравнения (7.58)-(7.60) принимают вид соответственно
где
Сама процедура вычислений очень похожа на соответствующую процедуру в случае орбиты, отличающейся от прямолинейной. При
по величинам
получаем
и вычисляем правые части уравнений (7.62)-(7.64). Если спрямление опорной орбиты производится в конце каждого шага, то
так что
в начале шага. Внутри шага
, вообще говоря, отклоняется от
. Теперь для того, чтобы получить значения
при
проинтегрируем нашу систему на один шаг.
Чтобы определить
при
надо вспомнить, что опорная орбита, все время является в точности прямолинейной. Поэтому
имеют те же значения, что и при
. Оставшиеся величины
в течение шага меняются и могут быть получены из соответствующих соотношений для прямолинейных эллипса, параболы или гиперболы (см. разд. 4.8).
Выбор в качестве опорной орбиты на данном шаге прямолинейных эллипса, параболы или гиперболы зависит от того, будет ли в начале шага энергия
отрицательной, нулевой или положительной величиной.
Значения
при
вычисляются по формулам
при
находятся из уравнений (7.54) и (7.55).
Теперь можно улучшить опорную орбиту. Если оскулирующая орбита все еще близка к прямолинейной, то в качестве новой опорной траектории опять берется прямолинейная орбита. Ее параметры в начале нового шага задаются следующим образом:
Новое значение энергии
определяет тип новой прямолинейной опорной орбиты (эллипс, парабола или гипербола). Когда момент количества движения h станет достаточно большим, следует использовать метод, разработанный для случая ненулевого к.