Прикладная математика и механика, 2023, T. 87, № 6, стр. 915-953

Кватернионная регуляризация особенностей моделей астродинамики, порождаемых гравитационными силами (обзор)

Ю. Н. Челноков 1*

1 Институт проблем точной механики и управления РАН
Саратов, Россия

* E-mail: ChelnokovYuN@gmail.com

Поступила в редакцию 04.09.2023
После доработки 11.10.2023
Принята к публикации 15.10.2023

Полный текст (PDF)

Аннотация

Излагается регуляризация особенностей дифференциальных уравнений возмущенной пространственной задачи двух тел, порождаемых гравитационными силами, с использованием четырехмерных переменных и матриц Кустаанхеймо–Штифеля, а также кватернионная регуляризация уравнений этой задачи, предложенная автором и имеющая ряд преимуществ перед матричной регуляризацией Кустаанхеймо–Штифеля. Дается аналитический обзор работ, посвященных кватернионной регуляризации указанных особенностей с использованием переменных Кустаанхеймо–Штифеля, которая уникальна в совместной регуляризации, линеаризации и увеличении размерности для трехмерных кеплеровских систем. Рассмотрен предложенный автором новый метод регуляризации уравнений возмущенной пространственной задачи двух тел, основанный на использовании идеальных прямоугольных координат Ганзена, переменных Леви-Чивита и параметров Эйлера (Родрига–Гамильтона), а также на использовании в качестве дополнительных переменных кеплеровской энергии и реального времени и новой независимой переменной Зундмана. Приведены регулярные кватернионные уравнения в переменных Леви-Чивита и параметрах Эйлера этой задачи, которые имеют не только хорошо известные достоинства матричных уравнений в переменных Кустаанхеймо–Штифеля, но и обладают своими дополнительными достоинствами.

Ключевые слова: механика космического полета (астродинамика), возмущенная пространственная задача двух тел, регуляризация особенностей, порождаемых гравитационными силами, идеальная система координат, уравнения орбитального движения, переменные Кустаанхеймо–Штифеля, параметры Эйлера (Родрига–Гамильтона), координаты Ганзена, переменные Леви-Чивита, кватернион

1. Введение. Работа носит обзорный аналитический характер и дополняет собой нашу обзорную работу [1]. Обсуждается проблема регуляризации дифференциальных уравнений возмущенной пространственной задачи двух тел, лежащих в основе небесной механики и механики космического полета (астродинамики): проблема устранения особенностей типа сингулярностей (деления на ноль), которые порождаются действующими на небесное или космическое тело ньютоновскими гравитационными силами и которые осложняют аналитическое и численное исследование движения тела вблизи гравитирующих тел или его движения по сильно вытянутым орбитам. Для получения различных регулярных уравнений используются четырехмерные параметры Эйлера (Родрига–Гамильтона), четырехмерные переменные Кустаанхеймо–Штифеля (KS-переменные), двухмерные идеальные прямоугольные координаты Ганзена, регулярные двухмерные переменные Леви-Чивита и кватернионы Гамильтона.

В нашей работе [1] излагаются кватернионные и бикватернионные методы описания движения, модели теории конечных перемещений и регулярной кинематики твердого тела, основанные на использовании четырехмерных вещественных и дуальных параметров Эйлера (Родрига–Гамильтона). Эти модели кинематики не имеют особенностей, порождаемых использованием для описания движения классических углов Эйлера–Крылова и их дуальных аналогов. Обсуждается проблема регуляризации дифференциальных уравнений возмущенной пространственной задачи двух тел с помощью использования четырехмерных параметров Эйлера, четырехмерных переменных Кустаанхеймо–Штифеля и кватернионов Гамильтона. Излагается история проблемы регуляризации и регулярные уравнения Кустаанхеймо–Штифеля, нашедшие широкое применение в небесной механике и астродинамике. Излагаются предложенные нами кватернионные методы регуляризации, имеющие ряд преимуществ перед матричной регуляризацией Кустаанхеймо–Штифеля, и различные, предложенные нами, регулярные кватернионные уравнения возмущенной пространственной задачи двух тел (как для абсолютного, так и для относительного движения). Приведены результаты сравнительного исследования точности численного интегрирования различных форм регуляризованных уравнений небесной механики и астродинамики в переменных Кустаанхеймо–Штифеля и ньютоновских уравнений в декартовых координатах, показывающие, что точность численного интегрирования регуляризованных уравнений в переменных Кустаанхеймо–Штифеля значительно выше (на несколько порядков) точности численного интегрирования ньютоновских уравнений.

В настоящей работе кратко описывается матричная регуляризация Кустаанхеймо–Штифеля (KS-регуляризация) особенностей дифференциальных уравнений возмущенной пространственной задачи двух тел, порождаемых гравитационными силами, с использованием четырехмерных переменных Кустаанхеймо–Штифеля. Показывается, что предложенный нами кватернионный метод регуляризации уравнений этой задачи имеет ряд преимуществ перед матричной регуляризацией Кустаанхеймо–Штифеля. Дается аналитический анализ работ, посвященных кватернионной регуляризации указанных особенностей дифференциальных уравнений возмущенной пространственной задачи двух тел с использованием четырехмерных переменных Кустаанхеймо–Штифеля. Многие из этих работ опубликованы в ведущих зарубежных изданиях.

В работе также излагается регуляризация Леви-Чивита уравнений возмущенной плоской задачи двух тел (уравнений плоского движения). Изложен предложенный нами новый метод регуляризации уравнений возмущенной пространственной задачи двух тел, основанный на использовании идеальной системы координат, двухмерных идеальных прямоугольных координат Ганзена и двухмерных переменных Леви-Чивита, описывающих движение второго (рассматриваемого) тела в идеальной системе координат, в которой уравнения пространственного движения принимают вид уравнений плоского движения, а также основанный на использовании четырехмерных параметров Эйлера (Родрига–Гамильтона) и кватерниона Гамильтона, характеризующих ориентацию идеальной системы координат в инерциальной системе координат. Эти параметры Эйлера и кватернион ориентации идеальной системы координат являются скалярными и кватернионным оскулирующими элементами орбиты изучаемого (второго) тела (медленно изменяющимися переменными).

Леви-Чивита в отношении своих попыток обобщить предложенную им знаменитую регуляризацию уравнений плоской задачи двух тел на пространственную задачу позже признал [2]: “Проблема в пространстве долго сопротивлялась моим усилиям, т.к. я пытался подойти к ней с помощью аналогичных изменений координат…”.

Штифель и Шейфеле в своей широко цитируемой книге [3] отмечали, что Леви-Чивита приложил много усилий, чтобы найти обобщение своего метода регуляризации дифференциальных уравнений плоского движения в задаче двух тел на общую пространственную задачу двух тел, но безуспешно. В работе [4] (см. также книгу [5]) говорится, что из-за фундаментальных трудностей, первоначально разъясненных Хопфом [6] и Гурвицем [7], невозможно обобщить преобразование Леви-Чивита к эквивалентному набору трехмерных переменных (на случай трехмерного пространства).

Тем не менее, автором статьи [8] было показано, что регуляризация Леви-Чивита может быть с успехом использована для построения регулярных уравнений возмущенной пространственной задачи двух тел. Эта регуляризация регулярных уравнений возмущенной пространственной задачи двух тел, предложенная в [8], излагается в нашей обзорной работе.

Для получения регулярных уравнений исходные дифференциальные уравнения возмущенной пространственной задачи двух тел в декартовых координатах записываются нами в двух системах координат: орбитальной и идеальной системах координат (название системы координат “идеальная”, по-видимому, было введено Deprit [9]) с использованием кватернионного оскулирующего элемента орбиты, компоненты которого – параметры Эйлера – являются медленно изменяющимися переменными и характеризуют инерциальную ориентацию в пространстве идеальной системы координат. В этих уравнениях переменными являются расстояние r от центра масс второго тела до центра масс первого тела, производная dr/dt (проекция ${{{v}}_{1}}$ вектора скорости v центра масс второго тела на направление радиус-вектора r второго (изучаемого) тела), модуль $c = \left| {{\mathbf{r}} \times {\mathbf{v}}} \right|$ вектора c момента орбитальной скорости второго тела, обобщенная истинная аномалия φ и параметры Эйлера (Родрига–Гамильтона) ${{\Lambda }_{j}}$ ( j = 0, 1, 2, 3) (скалярные оскулирующие элементы орбиты).

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

Затем в полученных уравнениях нами осуществляется переход от координат Ганзена к переменным Леви-Чивита и переход к новой независимой переменной τ в соответствии с дифференциальным соотношением dt = rdτ (преобразованием времени Зундмана [10]). Кроме этого, также в качестве дополнительной переменной вводится кеплеровская энергия h, через которую выражаются коэффициенты дифференциальных уравнений в переменных Леви-Чивита, и дополнение полученных уравнений дифференциальными уравнениями для кеплеровской энергии и времени.

В итоге получаются скалярные и кватернионные регулярные уравнения движения второго тела (регулярные уравнения возмущенной пространственной задачи двух тел), полученные в работе [8], в которых совместно используются переменных Леви-Чивита и параметры Эйлера (Родрига–Гамильтона).

Эти регулярные уравнения имеют не только хорошо известные достоинства уравнений Кустаанхеймо–Штифеля (регулярность, линейность в новом времени для кеплеровских движений, близость к линейным уравнениям для возмущенных движений), но и обладают своими дополнительными достоинствами:

1) для невозмущенного эллиптического кеплеровского движения изучаемого тела они эквивалентны уравнениям движения не четырехмерного одночастотного гармонического осциллятора, как в случае Кустаанхеймо–Штифеля, а уравнениям движения двухмерного одночастотного гармонического осциллятора; кватернион ориентации идеальной системы координат, в которой записаны эти уравнения движения, в этом случае остается постоянным;

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

Отметим, однако, что эти уравнения не пригодны для исследования редко встречающихся прямолинейных орбит, когда модуль вектора момента орбитальной скорости второго тела $c = \left| {{\mathbf{r}} \times {\mathbf{v}}} \right|$ обращается в ноль, поскольку кватернионное дифференциальное уравнение ориентации идеальной системы координат в этом случае вырождается (в знаменателях коэффициентов этого уравнения присутствует величина c). В статье указывается, как можно избавиться от этого недостатка излагаемых регулярных уравнений возмущенной пространственной задачи двух тел, в которых совместно используются переменные Леви-Чивита и параметры Эйлера, и приводятся регулярные уравнения в этих переменных, не имеющие указанного недостатка, но являющиеся более сложными.

2. Проблема регуляризации уравнений небесной механики и механики космического полета (астродинамики). В основе небесной механики и астродинамики лежат ньютоновские дифференциальные уравнения возмущенной пространственной задачи двух тел. Эти уравнения вырождаются при соударении второго (изучаемого) тела с первым (центральным) телом (при равенстве нулю расстояния между телами), что делает использование этих уравнений неудобным при изучении движения второго тела в малой окрестности центрального тела или его движения по сильно вытянутым орбитам. Сингулярность в начале координат создает в задаче двух тел не только теоретические, но и практические (вычислительные) трудности. Устранение особенностей типа сингулярности (деления на ноль) классических (ньютоновских) уравнений небесной механики и астродинамики, порождаемых силами гравитации, получило название “регуляризация” (Леви-Чивита 1920), а уравнения, не имеющие этих особенностей, называются регулярными. Среди методов регуляризации и регулярных моделей прикладной небесной механики и астродинамики в последнее время широкое распространение получили кватернионные методы и модели, основанные на использовании гиперкомплексных переменных – кватернионов Гамильтона, компонентами (элементами) которых являются четырехмерные переменные Кустаанхеймо–Штифеля или четырехмерные параметры Эйлера (Родрига–Гамильтона). Эти методы и модели имеют ряд преимуществ аналитического и вычислительного характера перед другими методами и моделями.

Изучению различных аспектов кватернионной регуляризации дифференциальных уравнений возмущенной пространственной задачи двух тел с использованием четырехмерных переменных Кустаанхеймо–Штифеля (KS-переменных) посвящены работы [1128], а также работы автора статьи [8, 2942].

Приводятся [4250] результаты сравнения численного решения уравнений орбитального движения небесных и космических тел в переменных Кустаанхеймо–Штифеля, параметрах Эйлера и в других переменных, которые свидетельствуют об эффективности использования KS-переменных и параметров Эйлера в задачах небесной механики и астродинамики.

Проведено [51, 52] сравнительное исследование точности численного интегрирования классических ньютоновских дифференциальных уравнений пространственной ограниченной задачи трех тел (Земля, Луна и космический аппарат) в декартовых координатах и построенных автором статьи [53] регулярных кватернионных дифференциальных уравнений этой задачи в KS-переменных, принимающих вид регулярных кватернионных уравнений возмущенной пространственной задачи двух тел в случае отсутствия поля тяготения Луны. Регулярные кватернионные уравнения в KS-переменных показали значительно более высокую точность, чем уравнения в декартовых координатах: для круговой орбиты точность оказалась выше на 2 порядка, для возмущенных эллиптических орбит со средним эксцентриситетом – на 4 порядка, для возмущенной эллиптической орбиты с высоким эксцентриситетом – на 7 порядков.

Отметим, что в книге Бордовицыной [43] приведены результаты численных исследований решений уравнений невозмущенной и возмущенной пространственной задачи двух тел (решений уравнений невозмущенного и возмущенного движения ИСЗ) ряда авторов с использованием уравнений в KS-переменных и уравнений в декартовых координатах, демонстрирующие преимущество уравнений в KS-переменных перед уравнениями в декартовых координатах (в смысле точности их численного интегрирования). Сравнение этих результатов с нашими результатами [51, 52] показало, что они в целом согласуются между собой. Отметим также, что Бордовицыной и Шарковским использованы регулярные уравнения в переменных Кустаанхеймо–Штифеля в канонической, а не в осцилляторной форме, выведенные ими с использованием соответствующего гамильтониана.

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

3. Регуляризация Кустаанхеймо–Штифеля уравнений возмущенной пространственной задачи двух тел и регуляризация Леви-Чивита уравнений возмущенной плоской задачи двух тел. Векторное и скалярные дифференциальные уравнения возмущенной пространственной задачи двух тел имеют вид

(3.1)
$\begin{gathered} \frac{{{{d}^{2}}{\mathbf{r}}}}{{d{{t}^{2}}}} + \frac{{f\left( {m + M} \right)}}{{{{r}^{3}}}}{\mathbf{r}} = {\mathbf{p}}\left( {t,{\mathbf{r}},\frac{{d{\mathbf{r}}}}{{dt}}} \right);\quad r = \left| {\mathbf{r}} \right| \\ \frac{{{{d}^{2}}{{\xi }_{1}}}}{{d{{t}^{2}}}} + \frac{{f\left( {m + M} \right)}}{{{{r}^{3}}}}{{\xi }_{1}} = {{p}_{1}},\quad \frac{{{{d}^{2}}{{\xi }_{2}}}}{{d{{t}^{2}}}} + \frac{{f\left( {m + M} \right)}}{{{{r}^{3}}}}{{\xi }_{2}} = {{p}_{2}} \\ \frac{{{{d}^{2}}{{\xi }_{3}}}}{{d{{t}^{2}}}} + \frac{{f\left( {m + M} \right)}}{{{{r}^{3}}}}{{\xi }_{3}} = {{p}_{3}};\quad r = \sqrt {\xi _{1}^{2} + \xi _{2}^{2} + \xi _{3}^{2}} , \\ \end{gathered} $
где r – радиус-вектор центра масс второго (изучаемого) тела, проводимый из центра масс первого (центрального) тела, $r = \left| {\mathbf{r}} \right|$ – расстояние между телами, m и M – массы второго и первого тел; f – гравитационная постоянная; p – вектор возмущающего ускорения центра масс второго тела (вектор-функция времени t, радиус-вектора r и вектора скорости ${\mathbf{v}} = d{\mathbf{r}}{\text{/}}dt$ второго тела в системе координат ξ, имеющей начало в центре масс первого тела и координатные оси, параллельные осям инерциальной системы координат), t – время; ξk – декартовые координаты центра масс изучаемого тела в системе координат ξ, ${{p}_{k}}$ = ${{p}_{k}}\left( {t,{{\xi }_{1}},{{\xi }_{2}},{{\xi }_{3}},{{{\dot {\xi }}}_{1}},{{{\dot {\xi }}}_{2}},{{{\dot {\xi }}}_{3}}} \right)$ $\left( {k = 1,2,3} \right)$ – проекции возмущающего ускорения p на оси инерциальной системы координат, совпадающие с его проекциями на оси системы координат ξ, верхняя точка – символ дифференцирования по времени t.

Уравнения (3.1) вырождаются при соударении второго тела с центральным телом (при равенстве нулю расстояния r между телами), что делает использование этих уравнений неудобным при изучении движения второго тела в малой окрестности центрального тела или его движения по сильно вытянутым орбитам. Сингулярность в начале координат, как уже отмечалось, создает не только теоретические, но и практические (вычислительные) трудности.

Проблема устранения указанной особенности, известная в небесной механике и астродинамике как проблема регуляризации дифференциальных уравнений возмущенной задачи двух тел, восходит к Эйлеру [54] и Леви-Чивита [2, 55, 56]. Наиболее эффективная регуляризация уравнений возмущенной пространственной задачи двух тел (спинорная или KS-регуляризация) была предложена Кустаанхеймо и Штифелем [57, 58]. Она наиболее полно изложена в монографии [3].

В регуляризации Кустаанхеймо использованы достоинства методов теории спиноров: вместо одной комплексной переменной теории Леви-Чивита была взята пара комплексных чисел. На языке вещественного анализа это эквивалентно введению четырех параметров ${{u}_{1}}$, ${{u}_{2}}$, ${{u}_{3}}$, ${{u}_{4}} = {{u}_{0}}$.

В регуляризации Кустаанхеймо–Штифеля, изложенной в книге [3], использована обобщенная матрица Леви-Чивита, названная KS-матрицей. Эта матрица (обозначается нами как $L\left( {{{{\mathbf{u}}}_{{KS}}}} \right)$) – четырехмерная квадратная матрица, содержащая в левом верхнем углу двухмерную квадратную матрицу Леви-Чивита и имеющая вид

(3.2)
$L\left( {{{{\mathbf{u}}}_{{KS}}}} \right) = \left( {\begin{array}{*{20}{c}} {{{u}_{1}}}&{ - {{u}_{2}}}&{ - {{u}_{3}}}&{{{u}_{0}}} \\ {{{u}_{2}}}&{{{u}_{1}}}&{ - {{u}_{0}}}&{ - {{u}_{3}}} \\ {{{u}_{3}}}&{{{u}_{0}}}&{{{u}_{1}}}&{{{u}_{2}}} \\ {{{u}_{0}}}&{ - {{u}_{3}}}&{{{u}_{2}}}&{ - {{u}_{1}}} \end{array}} \right)$

В этой матрице ${{u}_{j}}$j = 0, 1, 2, 3) – переменные Кустаанхеймо–Штифеля, называемые KS-переменными (вместо обозначения ${{u}_{4}}$ одной из переменных Кустаанхеймо–Штифеля здесь и далее нами используется обозначение ${{u}_{0}}$).

В матричной записи преобразование Кустаанхеймо–Штифеля имеет вид (3.3) [3]:

(3.3)
$\left( {\begin{array}{*{20}{c}} {{{\xi }_{1}}} \\ {{{\xi }_{2}}} \\ {{{\xi }_{3}}} \\ 0 \end{array}} \right) = \left( {\begin{array}{*{20}{c}} {{{u}_{1}}}&{ - {{u}_{2}}}&{ - {{u}_{3}}}&{{{u}_{0}}} \\ {{{u}_{2}}}&{{{u}_{1}}}&{ - {{u}_{0}}}&{ - {{u}_{3}}} \\ {{{u}_{3}}}&{{{u}_{0}}}&{{{u}_{1}}}&{{{u}_{2}}} \\ {{{u}_{0}}}&{ - {{u}_{3}}}&{{{u}_{2}}}&{ - {{u}_{1}}} \end{array}} \right)\left( {\begin{array}{*{20}{c}} {{{u}_{1}}} \\ {{{u}_{2}}} \\ {{{u}_{3}}} \\ {{{u}_{0}}} \end{array}} \right) = L\left( {{{{\mathbf{u}}}_{{KS}}}} \right)\left( {\begin{array}{*{20}{c}} {{{u}_{1}}} \\ {{{u}_{2}}} \\ {{{u}_{3}}} \\ {{{u}_{0}}} \end{array}} \right)$

В скалярной записи связи декартовых координат ξk с переменными Кустаанхеймо–Штифеля ${{u}_{j}}$ имеют вид

(3.4)
${{\xi }_{1}} = u_{0}^{2} + u_{1}^{2} - u_{2}^{2} - u_{3}^{2},\quad {{\xi }_{2}} = 2\left( {{{u}_{1}}{{u}_{2}} - {{u}_{0}}{{u}_{3}}} \right),\quad {{\xi }_{3}} = 2\left( {{{u}_{1}}{{u}_{3}} + {{u}_{0}}{{u}_{2}}} \right),$
что с точностью до перестановки индексов совпадает с отображением Хопфа [59].

Таким образом, в основе регуляризации Кустаанхеймо–Штифеля лежит нелинейное неоднозначное преобразование декартовых координат изучаемого тела, которое основывается на переходе от трехмерных декартовых координат ξk к новым четырехмерным переменным Кустаанхеймо–Штифеля uj (т.е., на переходе от трехмерного пространства декартовых координат к новому четырехмерному пространству).

Уравнения Кустаанхеймо–Штифеля в скалярной записи имеют вид [3]

(3.5)
$\frac{{{{d}^{2}}{{u}_{j}}}}{{d{{\tau }^{2}}}} - \frac{1}{2}h{{u}_{j}} = \frac{1}{2}r{{q}_{j}};\quad j = 0,1,2,3$
(3.6)
$\frac{{dh}}{{d\tau }} = 2\left( {{{q}_{0}}\frac{{d{{u}_{0}}}}{{d\tau }} + {{q}_{1}}\frac{{d{{u}_{1}}}}{{d\tau }} + {{q}_{2}}\frac{{d{{u}_{2}}}}{{d\tau }} + {{q}_{3}}\frac{{d{{u}_{3}}}}{{d\tau }}} \right)$
(3.7)
$\frac{{dt}}{{d\tau }} = r,\quad r = \left| {\mathbf{r}} \right| = u_{0}^{2} + u_{1}^{2} + u_{2}^{2} + u_{3}^{2}$
(3.8)
$\begin{gathered} {{q}_{0}} = {{u}_{0}}{{p}_{1}} - {{u}_{3}}{{p}_{2}} + {{u}_{2}}{{p}_{3}},\quad {{q}_{1}} = {{u}_{1}}{{p}_{1}} + {{u}_{2}}{{p}_{2}} + {{u}_{3}}{{p}_{3}} \\ {{q}_{2}} = - {{u}_{2}}{{p}_{1}} + {{u}_{1}}{{p}_{2}} + {{u}_{0}}{{p}_{3}},\quad {{q}_{3}} = - {{u}_{3}}{{p}_{1}} - {{u}_{0}}{{p}_{2}} + {{u}_{1}}{{p}_{3}} \\ \end{gathered} $

Здесь τ – новая независимая переменная, называемая фиктивным временем, связанная с временем t дифференциальным соотношением dt = rdτ (преобразованием времени Зундмана [10]), h – кеплеровская энергия единицы массы изучаемого тела, рассматриваемая как дополнительная переменная и определяемая соотношением

(3.9)
$h = \frac{1}{2}{{{v}}^{2}} - f\left( {m + M} \right)\frac{1}{r};\quad {v} = \left| {\mathbf{v}} \right|,\quad {\mathbf{v}} = \frac{{d{\mathbf{r}}}}{{dt}}$

Время t также рассматривается как дополнительная (зависимая) переменная.

Уравнения Кустаанхеймо–Штифеля (3.5)–(3.8) образуют систему десяти обыкновенных нелинейных, в общем случае нестационарных, дифференциальных уравнений относительно четырех KS-переменных uj, их первых производных duj/dτ по новой независимой переменной τ, кеплеровской энергии h и времени t.

Эти скалярные уравнения эквивалентны матричным уравнениям [3]

(3.10)
$\begin{gathered} \frac{{{{d}^{2}}{{{\mathbf{u}}}_{{ks}}}}}{{d{{\tau }^{2}}}} - \frac{1}{2}h{{{\mathbf{u}}}_{{ks}}} = \frac{1}{2}rL\left( {{{{\mathbf{u}}}_{{ks}}}} \right){{{\mathbf{P}}}_{{ks}}},\quad \frac{{dh}}{{d\tau }} = - 2\left( {\frac{{d{{{\mathbf{u}}}_{{ks}}}}}{{d\tau }},{{{\left( {L\left( {{{{\mathbf{u}}}_{{ks}}}} \right)} \right)}}^{{\text{T}}}}{{{\mathbf{P}}}_{{ks}}}} \right) \\ \frac{{dt}}{{d\tau }} = ({{{\mathbf{u}}}_{{ks}}},{{{\mathbf{u}}}_{{ks}}})\quad {{{\mathbf{u}}}_{{ks}}} = \left( {{{u}_{1}},{{u}_{2}},{{u}_{3}},{{u}_{0}}} \right),\quad {{{\mathbf{P}}}_{{ks}}} = \left( {{{p}_{1}},{{p}_{2}},{{p}_{3}},0} \right), \\ \end{gathered} $
где ${{{\mathbf{u}}}_{{ks}}}$ – четырехмерный вектор-столбец KS-переменных, L(uks) – KS-матрица, определяемая соотношением (3.2), ${{{\mathbf{P}}}_{{ks}}}$ – четырехмерный вектор-столбец, сопоставляемый трехмерному вектору возмущающего ускорения p, $\left( {{\mathbf{a}},{\mathbf{b}}} \right)$ – скалярное произведение четырехмерных вектор-столбцов a и b.

Уравнения Кустаанхеймо–Штифеля обобщают уравнения Леви-Чивита уравнений возмущенной плоской задачи двух тел, имеющие в скалярной записи вид

(3.11)
$\frac{{{{d}^{2}}{{u}_{j}}}}{{d{{\tau }^{2}}}} - \frac{1}{2}h{{u}_{j}} = \frac{1}{2}r{{q}_{j}};\quad j = 1,2$
(3.12)
$\frac{{dh}}{{d\tau }} = 2\left( {{{q}_{1}}\frac{{d{{u}_{1}}}}{{d\tau }} + {{q}_{2}}\frac{{d{{u}_{2}}}}{{d\tau }}} \right);\quad \frac{{dt}}{{d\tau }} = r,\quad r = \left| {\mathbf{r}} \right| = u_{1}^{2} + u_{2}^{2}$
(3.13)
${{q}_{1}} = {{u}_{1}}{{p}_{1}} + {{u}_{2}}{{p}_{2}},\quad {{q}_{2}} = - {{u}_{2}}{{p}_{1}} + {{u}_{1}}{{p}_{2}}$

Здесь ${{u}_{1}}$ и ${{u}_{2}}$ – переменные Леви-Чивита, связанные с декартовыми координатами ${{\xi }_{1}}$ и ${{\xi }_{2}}$ соотношениями

(3.14)
${{\xi }_{1}} = u_{1}^{2} - u_{2}^{2},\quad {{\xi }_{2}} = 2{{u}_{1}}{{u}_{2}},\quad {{\left( {\xi _{1}^{2} + \xi _{2}^{2}} \right)}^{{1{\text{/}}2}}} = r = \left| {\mathbf{r}} \right| = u_{1}^{2} + u_{2}^{2}$

В случае плоского движения четырехмерная квадратная матрица $L\left( {{{{\mathbf{u}}}_{{KS}}}} \right)$ Кустаанхеймо–Штифеля, определяемая соотношением (3.2), переходит в двухмерную квадратную матрицу Леви-Чивита, имеющую вид

(3.15)
$L\left( {{{{\mathbf{u}}}_{{lc}}}} \right) = \left( {\begin{array}{*{20}{c}} {{{u}_{1}}}&{ - {{u}_{2}}} \\ {{{u}_{2}}}&{{{u}_{1}}} \end{array}} \right)$

В матричной записи преобразование Леви-Чивита (3.14) имеет вид (3.16) [3]:

(3.16)
$\left( {\begin{array}{*{20}{c}} {{{\xi }_{1}}} \\ {{{\xi }_{2}}} \end{array}} \right) = \left( {\begin{array}{*{20}{c}} {{{u}_{1}}}&{ - {{u}_{2}}} \\ {{{u}_{2}}}&{{{u}_{1}}} \end{array}} \right)\;\left( {\begin{array}{*{20}{c}} {{{u}_{1}}} \\ {{{u}_{2}}} \end{array}} \right) = L\left( {{{{\mathbf{u}}}_{{lc}}}} \right)\;\left( {\begin{array}{*{20}{c}} {{{u}_{1}}} \\ {{{u}_{2}}} \end{array}} \right)$

В комплексной записи преобразование Леви-Чивита принимает вид (3.17):

(3.17)
${{\xi }_{1}} + i{{\xi }_{2}} = {{\left( {{{u}_{1}} + i{{u}_{2}}} \right)}^{2}}$
и устанавливает отображение параметрической плоскости переменных Леви-Чивита ${{u}_{1}}$ и ${{u}_{2}}$ на физическую плоскость переменных ${{\xi }_{1}}$ и ${{\xi }_{2}}$ (декартовых координат).

Скалярные регулярные уравнения возмущенной плоской задачи двух тел (3.11)–(3.13) Леви-Чивита эквивалентны его матричным уравнениям [3]

(3.18)
$\begin{gathered} \frac{{{{d}^{2}}{{{\mathbf{u}}}_{{lc}}}}}{{d{{\tau }^{2}}}} - \frac{1}{2}h{{{\mathbf{u}}}_{{lc}}} = \frac{1}{2}rL\left( {{{{\mathbf{u}}}_{{lc}}}} \right){{{\mathbf{P}}}_{{lc}}},\quad \frac{{dh}}{{d\tau }} = - 2\left( {\frac{{d{{{\mathbf{u}}}_{{lc}}}}}{{d\tau }},{{{\left( {L\left( {{{{\mathbf{u}}}_{{lc}}}} \right)} \right)}}^{{\text{T}}}}{{{\mathbf{P}}}_{{lc}}}} \right);\quad \frac{{dt}}{{d\tau }} = ({{{\mathbf{u}}}_{{lc}}},{{{\mathbf{u}}}_{{lc}}}) \\ {{{\mathbf{u}}}_{{lc}}} = \left( {{{u}_{1}},{{u}_{2}}} \right),\quad {{{\mathbf{P}}}_{{lc}}} = \left( {{{p}_{1}},{{p}_{2}}} \right), \\ \end{gathered} $
где ${{{\mathbf{u}}}_{{lc}}}$ – двухмерный вектор-столбец переменных Леви-Чивита, L(ulc) – матрица Леви-Чивита, определяемая соотношением (3.15), ${{{\mathbf{P}}}_{{lc}}}$ – двухмерный вектор-столбец, сопоставляемый двухмерному вектору возмущающего ускорения p, $\left( {{\mathbf{a}},{\mathbf{b}}} \right)$ – скалярное произведение двухмерных вектор-столбцов a и b.

Отметим основные хорошо известные достоинства уравнений Кустаанхеймо–Штифеля (3.5)–(3.7):

− они, в отличие от ньютоновских уравнений (3.1), регулярны в центре притяжения, когда расстояние между телами $r = 0$;

− линейны во времени τ для невозмущенных кеплеровских движений (в отличие от существенно нелинейных ньютоновских уравнений) и имеют в этом случае вид системы четырех независимых линейных дифференциальных уравнений второго порядка с одинаковыми постоянными коэффициентами, равными половинной кеплеровской энергии h, взятой со знаком “–”:

$\frac{{{{d}^{2}}{{u}_{j}}}}{{d{{\tau }^{2}}}} - \frac{1}{2}h{{u}_{j}} = 0;\quad h = {\text{const,}}\quad j = 0,1,2,3$

− для эллиптического кеплеровского движения, когда кеплеровская энергия h = = const < 0, эти уравнения эквивалентны уравнениям движения четырехмерного одночастотного гармонического осциллятора во времени τ, квадрат частоты которого равен половине кеплеровской энергии, взятой со знаком минус;

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

− близки к линейным уравнениям для возмущенных кеплеровских движений;

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

Эти свойства регулярных уравнений позволили разработать эффективные методы нахождения решений в аналитической или численной форме таких трудных для классических методов задач как исследование движения вблизи притягивающих масс или движения по орбитам с большими эксцентриситетами. Штифелем и Шейфеле [3], Бордовицыной и Шарковским [43] показано, что точность численного интегрирования регуляризованных уравнений в переменных Кустаанхеймо–Штифеля значительно выше точности интегрирования ньютоновских уравнений. Так, Бордовицыной и Шарковским показано, что использование регулярных уравнений в переменных Кустаанхеймо–Штифеля позволяет повысить точность численного решения ряда задач небесной механики и астродинамики, например, задачи о движении искусственного спутника Земли (ИСЗ) по орбитам с большими эксцентриситетами, от трех до пяти порядков по сравнению с решениями, полученными при использовании классических ньютоновских уравнений. Сравнение этих результатов с нашими, приведенными выше, показало, как уже отмечалось выше, что они в целом согласуются между собой.

Fukushima было также показано [45, 46], что KS-регуляризация приводит к очень эффективной схеме интегрирования уравнений орбитального движения, повышающей точность и скорость численного интегрирования. Это связано не только со структурой уравнений, но также с использованием нескольких методов, которые приносят важные преимущества численной схеме. Проводится [46] численное сравнение четырех схем регуляризации трехмерной задачи двух тел в условиях возмущения: регуляризации KS, Шперлинга–Бюрде (SB), Бюрде–Феррандиса (BF) и трехмерное расширение регуляризации Леви-Чивита (LC). В аннотации [46] сказано: КS и расширенная -регуляризация с масштабированием энергии Кеплера обеспечивают наилучшую экономическую эффективность при интеграции почти всех возмущенных задач двух тел. Подчеркнем, что и Fukushima [46] сравнивает семь схем, описанных в § 1, а также “нерегуляризованную обработку, а именно прямое интегрирование в декартовых координатах”. Для каждой из этих формулировок им проводится тестовая интеграция Икара, охватывающая около 1 миллиона лет, и измеряется время его выполнения.

4. Кватернионная регуляризация и регулярные кватернионные уравнения возмущенной пространственной задачи двух тел, предложенные автором статьи. В основе знаменитой регуляризации Кустаанхеймо–Штифеля лежит, как уже отмечалось, нелинейное неоднозначное преобразование декартовых координат изучаемого тела, названное KS-преобразованием. Причем это преобразование состоит в переходе от трехмерного пространства декартовых координат к четырехмерному пространству новых координат (к четырехмерным переменным Кустаанхеймо–Штифеля). Поэтому, по мнению Штифеля и Шейфеля прямой вывод регулярных уравнений в трехмерном (т.е. пространственном) случае невозможен ([3], стр. 29).

В книге [3] постулируется матричное регулярное уравнение пространственной задачи двух тел (первое из уравнений (3.10)), записанное ими по аналогии с матричным регулярным уравнением Леви-Чивита плоского движения (3.18), и с помощью нескольких теорем доказывается, что при этом удовлетворяется старое векторное ньютоновское уравнение (первое из уравнений (3.1)).

Вскоре после открытия KS-преобразования было рассмотрено использование кватернионов Гамильтона (четырехмерных гиперкомплексных чисел) и четырехмерных кватернионных матриц для регуляризации уравнений пространственной задачи двух тел, поскольку четырехмерность пространства новых координат делало естественным использование для такой регуляризации кватернионов Гамильтона и четырехмерных кватернионных матриц. Однако, Штифель и Шейфеле полностью отвергли эту идею, написав в одиннадцатой главе своей книги [3], посвященной геометрии KS-преобразования, что “Любая попытка заменить теорию KS-матриц более популярной теорией кватернионных матриц приводит поэтому к неудаче или, во всяком случае, к очень громоздкому формализму”. Поэтому вместо хорошо известных в математике четырехмерных кватернионных матриц Штифель и Шейфеле в [3] для построения теории регуляризации используют новые четырехмерные квадратные матрицы, названные ими KS-матрицами, имеющие вид (3.2).

Приведенное утверждение Штифеля и Шейфеле было опровергнуто, по-видимому впервые, автором статьи [2932]. Известный западный ученый Вальдвогель по поводу цитированного выше высказывания Штифеля и Шейфеле о бесперспективности использования в теории регуляризации кватернионных матриц в работе [20] говорит: “Это утверждение было впервые опровергнуто Челноковым (1981), который представил теорию регуляризации пространственной задачи Кеплера, используя геометрические представления во вращающейся системе координат и кватернионные матрицы. В серии статей (напр., 1992 и 1999) тем же автором была расширена теория кватернионной регуляризации и приведены практические применения.”

Автором статьи было показано [2932] ([29] – с использованием классических кватернионных матриц, [30] – с использованием кватернионов Гамильтона), что в действительности кватернионный подход к регуляризации

– позволяет дать прямой и наглядный вывод регулярных уравнений в KS-переменных (что, как уже отмечалось, ставилось Штифелем и Шейфеле под сомнение ([3], cтр. 29) из-за неоднозначности KS-преобразования,

− позволяет дать наглядные геометрическую и кинематическую интерпретации регуляризующему KS преобразованию;

− раскрывает геометрический смысл его неоднозначности;

− позволяет получить более общие регулярные уравнения возмущенной пространственной задачи двух тел, частным случаем которых являются регулярные уравнения Кустаанхеймо–Штифеля.

Для построения кватернионных регулярных уравнений возмущенной пространственной задачи двух тел в [29] была использована четырехмерная кватернионная матрица

(4.1)
$n\left\{ {\mathbf{u}} \right\} = \left( {\begin{array}{*{20}{c}} {{{u}_{0}}}&{ - {{u}_{1}}}&{ - {{u}_{2}}}&{ - {{u}_{3}}} \\ {{{u}_{1}}}&{{{u}_{0}}}&{{{u}_{3}}}&{ - {{u}_{2}}} \\ {{{u}_{2}}}&{ - {{u}_{3}}}&{{{u}_{0}}}&{{{u}_{1}}} \\ {{{u}_{3}}}&{{{u}_{2}}}&{ - {{u}_{1}}}&{{{u}_{0}}} \end{array}} \right),$
элементами которой являются переменные Кустаанхеймо–Штифеля uj.

Эта матрица отличается от KS-матрицы (3.2). Если под элементами uj  этой матрицы понимать четырехмерные параметры Эйлера (Родрига–Гамильтона) ${{\lambda }_{j}}$, то она будет являться хорошо известной кватернионной матрицей поворота, описывающей вращение в трехмерном пространстве.

В работе [30] для этих целей использован кватернион Гамильтона u:

(4.2)
${\mathbf{u}} = {{u}_{0}} + {{u}_{1}}{\mathbf{i}} + {{u}_{2}}{\mathbf{j}} + {{u}_{3}}{\mathbf{k}},$
элементами которого являются переменные Кустаанхеймо–Штифеля uj.

В кватернионной записи связи декартовых координат ξk с переменными Кустаанхеймо–Штифеля uj имеют вид

(4.3)
${{{\mathbf{r}}}_{\xi }} = {{\xi }_{1}}{\mathbf{i}} + {{\xi }_{2}}{\mathbf{j}} + {{\xi }_{3}}{\mathbf{k}} = {\mathbf{\bar {u}}} \circ {\mathbf{i}} \circ {\mathbf{u}};\quad {\mathbf{\bar {u}}} = {{u}_{0}} - {{u}_{1}}{\mathbf{i}} - {{u}_{2}}{\mathbf{j}} - {{u}_{3}}{\mathbf{k}},$
приведенный Штифелем и Шейфеле в их книге [3].

Показано [29, 30], что регуляризующее KS-преобразование координат заключается в переходе от трехмерных декартовых координат центра масс второго тела в инерциальной системе координат к новым четырехмерным переменным, которые являются нормированными определенным образом компонентами сопряженного кватерниона поворота ${\mathbf{\bar {\lambda }}}$ = ${{\lambda }_{0}} - {{\lambda }_{1}}{\mathbf{i}}$${{\lambda }_{2}}{\mathbf{j}} - {{\lambda }_{3}}{\mathbf{k}}$, характеризующего ориентацию вращающейся системы координат η в инерциальной системе координат. Ось ${{\eta }_{1}}$ этой системы координат направлена вдоль радиус-вектора r центра масс второго тела, а ее начало находится в центре масс этого тела. Нормирующий множитель равен квадратному корню из расстояния r от центра масс второго тела до центра притяжения, поэтому кватернионные переменные u и λ связаны кватернионным соотношением ${\mathbf{u}} = \sqrt r {\kern 1pt} {\mathbf{\bar {\lambda }}}$, которые эквивалентны скалярным соотношениям: ${{u}_{0}} = {{r}^{{1/2}}}{{\lambda }_{0}}$, ${{u}_{i}} = - {{r}^{{1/2}}}{{\lambda }_{i}}$, $i = 1,2,3$.

Также было показано, что билинейное соотношение Кустаанхеймо–Штифеля

(4.4)
${{u}_{1}}\frac{{d{{u}_{0}}}}{{d\tau }} - {{u}_{0}}\frac{{d{{u}_{1}}}}{{d\tau }} + {{u}_{3}}\frac{{d{{u}_{2}}}}{{d\tau }} - {{u}_{2}}\frac{{d{{u}_{3}}}}{{d\tau }} = 0$
связывающее между собой KS-переменные uj и их первые производные по новой независимой переменной τ, накладывает на движение трехгранника η дополнительное (неголономное) условие, заключающееся в равенстве нулю проекции ${{\omega }_{1}}$ вектора ω абсолютной угловой скорости трехгранника η на направление радиус-вектора r (ось ${{\eta }_{1}}$):

$\begin{gathered} {{\omega }_{1}} = 2\left( {{{\lambda }_{0}}\frac{{d{{\lambda }_{1}}}}{{dt}} - {{\lambda }_{1}}\frac{{d{{\lambda }_{0}}}}{{dt}} - {{\lambda }_{2}}\frac{{d{{\lambda }_{3}}}}{{dt}} + {{\lambda }_{3}}\frac{{d{{\lambda }_{2}}}}{{dt}}} \right) = \\ = \frac{2}{r}\left( {{{u}_{1}}\frac{{d{{u}_{0}}}}{{d\tau }} - {{u}_{0}}\frac{{d{{u}_{1}}}}{{d\tau }} + {{u}_{3}}\frac{{d{{u}_{2}}}}{{d\tau }} - {{u}_{2}}\frac{{d{{u}_{3}}}}{{d\tau }}} \right) = 0 \\ \end{gathered} $

Отметим, что по словам Штифеля и Шейфеле ([3], стр. 29) соотношение вида (4.4):

${{u}_{4}}{{{v}}_{1}} - {{u}_{3}}{{{v}}_{2}} + {{u}_{2}}{{{v}}_{3}} - {{u}_{1}}{{{v}}_{4}} = 0,$
где uj и ${{{v}}_{j}}$ – компоненты четырехмерных векторов u и $v$, “называется билинейным соотношением и играет основную роль в нашем построении регулярной небесной механики”.

Таким образом, переход в уравнениях пространственной задачи двух тел от декартовых координат центра масс второго тела к KS-переменным фактически означает запись этих уравнений во вращающейся системе координат η с использованием в качестве параметров ориентации этой системы координат четырехмерных параметров Эйлера (Родрига–Гамильтона) ${{\lambda }_{j}}$, являющихся компонентами кватерниона поворота λ этой системы координат.

Дальнейшие преобразования этих уравнений связаны с нормировкой параметров Эйлера ${{\lambda }_{j}}$ (кватерниона поворота λ) с помощью множителя $\sqrt r $ и с переходом в них к переменным Кустаанхеймо–Штифеля uj, а также с введением в качестве дополнительных зависимых переменных кеплеровской энергии h и времени t и с переходом к новой независимой переменной Зундмана τ.

Отметим, что автором были получены [29, 30] более общие (в сравнении с уравнениями Кустаанхеймо–Штифеля) матричные (с использованием кватернионных матриц) и кватернионные регулярные уравнения возмущенной пространственной задачи двух тел в KS-переменных в предположении, что вышеуказанное билинейное соотношение (4.4) не выполняется. Эти уравнения содержат дополнительные слагаемые, в которых присутствуют проекции ${{\omega }_{1}}$ и ${{\varepsilon }_{1}}$ векторов угловой скорости и углового ускорения сопровождающего трехгранника η на направление радиус-вектора r центра масс второго тела (одна из этих проекций (${{\omega }_{1}}$ или ${{\varepsilon }_{1}}$) является произвольно задаваемым параметром) и являются более сложными.

Матричные регулярные уравнения возмущенной пространственной задачи двух тел в переменных Кустаанхеймо–Штифеля в общем случае, когда билинейное соотношение (4.4) не выполнется, имеют вид [29]

$\left( {\begin{array}{*{20}{c}} {u_{1}^{{''}}} \\ {u_{0}^{{''}}} \\ { - u_{3}^{{''}}} \\ {u_{2}^{{''}}} \end{array}} \right) - \frac{1}{2}{{\omega }_{1}}\left( {\begin{array}{*{20}{c}} 0&{ - {{\xi }_{1}} - 3r}&{ - {{\xi }_{2}}}&{ - {{\xi }_{3}}} \\ {{{\xi }_{1}} + 3r}&0&{ - {{\xi }_{3}}}&{{{\xi }_{2}}} \\ {{{\xi }_{2}}}&{{{\xi }_{3}}}&0&{ - {{\xi }_{1}} + 3r} \\ {{{\xi }_{3}}}&{ - {{\xi }_{2}}}&{{{\xi }_{1}} - 3r}&0 \end{array}} \right)\left( {\begin{array}{*{20}{c}} {u_{1}^{'}} \\ {u_{0}^{'}} \\ { - u_{3}^{'}} \\ {u_{2}^{'}} \end{array}} \right) - $
$ - \;\frac{1}{2}\left( {\begin{array}{*{20}{c}} {h + \omega _{1}^{2}{{r}^{2}}}&{ - {{r}^{2}}{{\varepsilon }_{1}}}&0&0 \\ {{{r}^{2}}{{\varepsilon }_{1}}}&{h + \omega _{1}^{2}{{r}^{2}}}&0&0 \\ 0&0&{h + \omega _{1}^{2}{{r}^{2}}}&{{{r}^{2}}{{\varepsilon }_{1}}} \\ 0&0&{ - {{r}^{2}}{{\varepsilon }_{1}}}&{h + \omega _{1}^{2}{{r}^{2}}} \end{array}} \right)\left( {\begin{array}{*{20}{c}} {{{u}_{1}}} \\ {{{u}_{0}}} \\ { - {{u}_{3}}} \\ {{{u}_{2}}} \end{array}} \right) = \frac{1}{2}r\left( {\begin{array}{*{20}{c}} {{{q}_{1}}} \\ {{{q}_{0}}} \\ { - {{q}_{3}}} \\ {{{q}_{2}}} \end{array}} \right)$
(4.5)
$\begin{gathered} h{\kern 1pt} ' = 2\left( {{{q}_{0}}u_{0}^{'} + {{q}_{1}}u_{1}^{'} + {{q}_{2}}u_{2}^{'} + {{q}_{3}}u_{3}^{'}} \right),\quad t{\kern 1pt} ' = r,\quad r = \left| {\mathbf{r}} \right| = u_{0}^{2} + u_{1}^{2} + u_{2}^{2} + u_{3}^{2} \\ {{q}_{0}} = {{u}_{0}}{{p}_{1}} - {{u}_{3}}{{p}_{2}} + {{u}_{2}}{{p}_{3}},\quad {{q}_{1}} = {{u}_{1}}{{p}_{1}} + {{u}_{2}}{{p}_{2}} + {{u}_{3}}{{p}_{3}}, \\ \end{gathered} $
${{q}_{2}} = - {{u}_{2}}{{p}_{1}} + {{u}_{1}}{{p}_{2}} + {{u}_{0}}{{p}_{3}},\quad {{q}_{3}} = - {{u}_{3}}{{p}_{1}} - {{u}_{0}}{{p}_{2}} + {{u}_{1}}{{p}_{3}}$
${{\xi }_{1}} = u_{0}^{2} + u_{1}^{2} - u_{2}^{2} - u_{3}^{2},\quad {{\xi }_{2}} = 2\left( {{{u}_{1}}{{u}_{2}} - {{u}_{0}}{{u}_{3}}} \right),\quad {{\xi }_{3}} = 2\left( {{{u}_{1}}{{u}_{3}} + {{u}_{0}}{{u}_{2}}} \right),$
где ${{\omega }_{1}}$ и ${{\varepsilon }_{1}} = d{{\omega }_{1}}{\text{/}}dt$ – проекции векторов угловой скорости и углового ускорения трехгранника η на направление радиус-вектора r (одна из этих проекций (${{\omega }_{1}}$ или ${{\varepsilon }_{1}}$) является произвольно задаваемым параметром), верхний штрих, по-прежнему, – символ дифференцирования по независимой переменной τ.

Кватернионные регулярные уравнения этой задачи в общем случае в переменных Кустаанхеймо–Штифеля имеют вид [30] (в этой работе эти уравнения записаны в другой кватернионной форме, соответствующей приведенной выше матричной записи этих уравнений (4.5), когда кватернионная переменная ${\mathbf{\alpha }}$, введенная в этой работе и соответствующая вектор-столбцу $\left( {{{u}_{1}},{{u}_{0}}, - {{u}_{3}},{{u}_{2}}} \right)$, равна ${\mathbf{\bar {u}}} \circ {\mathbf{i}}$):

(4.6)
$\begin{gathered} \frac{{{{d}^{2}}{\mathbf{u}}}}{{d{{\tau }^{2}}}} + \frac{3}{2}r{{\omega }_{1}}{\mathbf{i}} \circ \frac{{d{\mathbf{u}}}}{{d\tau }} + \frac{1}{2}\left( {{{r}^{2}}{{\varepsilon }_{1}} + {{\omega }_{1}}\frac{{d{\mathbf{u}}}}{{d\tau }} \circ {\mathbf{\bar {u}}}} \right) \circ {\mathbf{i}} \circ {\mathbf{u}} - \frac{1}{2}\left( {h + {{r}^{2}}\omega _{1}^{2}} \right){\mathbf{u}} = - \frac{1}{2}r{\mathbf{i}} \circ {\mathbf{u}} \circ {{{\mathbf{p}}}_{\xi }} \\ \frac{{dh}}{{d\tau }} = {{{\mathbf{p}}}_{\xi }} \cdot \left[ {\frac{d}{{d\tau }}\left( {{\mathbf{\bar {u}}} \circ {\mathbf{i}} \circ {\mathbf{u}}} \right)} \right];\quad \frac{{dt}}{{d\tau }} = r,\quad \frac{{d{{\omega }_{1}}}}{{d\tau }} = r{{\varepsilon }_{1}} \\ {\mathbf{u}} = {{u}_{0}} + {{u}_{1}}{\mathbf{i}} + {{u}_{2}}{\mathbf{j}} + {{u}_{3}}{\mathbf{k}},\quad r = \left| {\mathbf{r}} \right| = {\mathbf{\bar {u}}} \circ {\mathbf{u}} = {\mathbf{u}} \circ {\mathbf{\bar {u}}} = u_{0}^{2} + u_{1}^{2} + u_{2}^{2} + u_{3}^{2} \\ {{{\mathbf{p}}}_{\xi }} = {{p}_{1}}{\mathbf{i}} + {{p}_{2}}{\mathbf{j}} + {{p}_{3}}{\mathbf{k}},\quad {{p}_{k}} = {{p}_{k}}\left( {t,{{\xi }_{1}},{{\xi }_{2}},{{\xi }_{3}},{{{\dot {\xi }}}_{1}},{{{\dot {\xi }}}_{2}},{{{\dot {\xi }}}_{3}}} \right) \\ \end{gathered} $

В этих уравнениях u, h, t – неизвестные функции независимой переменной τ, ${\mathbf{\bar {u}}}$ – кватернион, сопряженный кватерниону u, проекции ${{p}_{k}}$ вектора p возмущающего ускорения на оси инерциальной системы координат являются заданными функциями времени t, декартовых координат и проекций вектора скорости центра масс второго тела на оси опорной системы координат ξ, которые могут быть представлены как функции времени t и переменных ${{u}_{j}}$, $u_{j}^{'}$; ${{\omega }_{1}}$ и ${{\varepsilon }_{1}} = d{{\omega }_{1}}{\text{/}}dt$ – проекции векторов угловой скорости и углового ускорения трехгранника η на ось ${{\eta }_{1}}$, заданные как функции переменных t и ${{u}_{j}}$, $u_{j}^{'}$; центральная точка означает скалярное произведение.

Для нахождения проекций радиус-вектора r и вектора скорости ${\mathbf{v}} = {\mathbf{\dot {r}}}$ центра масс изучаемого тела на оси инерциальной системы координат (координат ξk и их производных ${{\dot {\xi }}_{k}}$ по времени t) через переменные uj и $u_{j}^{'}$ необходимо воспользоваться кватернионными соотношениями

${{{\mathbf{r}}}_{\xi }} = {{\xi }_{1}}{\mathbf{i}} + {{\xi }_{2}}{\mathbf{j}} + {{\xi }_{3}}{\mathbf{k}} = {\mathbf{\bar {u}}} \circ {\mathbf{i}} \circ {\mathbf{u}},\quad {{{\mathbf{v}}}_{\xi }} = {{\dot {\xi }}_{1}}{\mathbf{i}} + {{\dot {\xi }}_{2}}{\mathbf{j}} + {{\dot {\xi }}_{3}}{\mathbf{k}} = 2{{r}^{{ - 1}}}{\mathbf{\bar {u}}} \circ {\mathbf{i}} \circ {\mathbf{u}}{\kern 1pt} '$

Таким образом, нами были получены более общие (в сравнении с уравнениями Кустаанхеймо–Штифеля) регулярные матричные и кватернионные уравнения возмущенной пространственной задачи двух тел. Эти уравнения сложнее уравнений Кустаанхеймо–Штифеля. Они, в первую очередь, имеют теоретический интерес, демонстрируя тот факт, что регуляризация достигается и в том случае, когда ${{\omega }_{1}} \ne 0$, т.е. когда не выполняется билинейное соотношение, являющееся одним из основных в теории регуляризации Кустаанхеймо–Штифеля. Возможно, эти уравнения будет целесообразно использовать для высокоточных численных расчетов, поскольку при их интегрировании не требуется выполнения билинейного соотношения (4.4), которое неизбежно будет нарушаться в процессе численного интегрирования из-за методических и вычислительных погрешностей.

Полагая ${{\omega }_{1}} = 0$ и ${{\varepsilon }_{1}} = 0$, из уравнений (4.6) получаем нашу кватернионную форму регулярных уравнений Кустаанхеймо–Штифеля:

(4.7)
$\begin{gathered} \frac{{{{d}^{2}}{\mathbf{u}}}}{{d{{\tau }^{2}}}} - \frac{1}{2}h{\mathbf{u}} = - \frac{1}{2}r{\mathbf{i}} \circ {\mathbf{u}} \circ {{{\mathbf{p}}}_{\xi }},\quad \frac{{dh}}{{d\tau }} = {\text{2}}\operatorname{scal} \left( {\frac{{d{\mathbf{\bar {u}}}}}{{d\tau }} \circ {\mathbf{q}}} \right);\quad \frac{{dt}}{{d\tau }} = r = {\mathbf{u}} \circ {\mathbf{\bar {u}}} \\ q = {{q}_{1}}{\mathbf{i}} + {{q}_{2}}{\mathbf{j}} + {{q}_{3}}{\mathbf{k}} = - {\mathbf{i}} \circ {\mathbf{u}} \circ {{{\mathbf{p}}}_{\xi }},\quad {{{\mathbf{p}}}_{\xi }} = {{p}_{1}}{\mathbf{i}} + {{p}_{2}}{\mathbf{j}} + {{p}_{3}}{\mathbf{k}} \\ \end{gathered} $

Здесь scal(•) – скалярная часть кватернионного произведения ${\mathbf{\bar {u}}}{\kern 1pt} '\; \circ {\mathbf{q}}$, ${{{\mathbf{p}}}_{\xi }}$ – кватернион действующих возмущений.

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

Автором также были получены другие регулярные кватернионные уравнения возмущенной пространственной задачи двух тел и искусственного спутника Земли (ИСЗ) в новых четырехмерных переменных: модифицированных переменных Кустаанхеймо–Штифеля, введенных в работе [34] (см. также [38]). Эти уравнения обладают всеми достоинствами выше приведенных матричных и кватернионных уравнений в переменных Кустаанхеймо–Штифеля, но имеют более простую и симметричную структуру для движения второго тела (например, космического аппарата (спутника)) в гравитационном поле Земли, в описании которого учитываются не только центральная, но и зональные, тессеральные и секториальные гармоники [34, 41].

Более простые и симметричные структуры уравнений приводят к более эффективным вычислительным алгоритмам при численном интегрировании дифференциальных уравнений движения спутника. Удобство и эффективность использования полученных уравнений движения спутника в модифицированных переменных Кустаанхеймо–Штифеля для аналитического исследования движения спутника показано в [41] на примере движения спутника в гравитационном поле Земли, в описании которого учитываются его центральная (ньютоновская) и зональные гармоники. В этой работе найдены первые интегралы уравнений движения спутника в модифицированных переменных в указанном случае, предложены замены переменных и преобразования этих уравнений, позволившие получить для изучения движения спутника замкнутые системы дифференциальных уравнений меньшей размерности, в частности, системы уравнений четвертого и третьего порядков.

Предложен [32, 34] способ однозначного нахождения значений параметров Эйлера и переменных Кустаанхеймо–Штифеля по заданным в рассматриваемый момент времени значениям декартовых координат и их первых производных (проекциям вектора скорости). Этот способ имеет в сравнении со способом неоднозначного нахождения переменных Кустаанхеймо–Штифеля [3] существенные преимущества.

5. Работы по регуляризации Кустаанхеймо–Штифеля (KS-регуляризации) и кватернионной регуляризации уравнений задачи двух тел других авторов. В работе Velte [11], посвященной профессору Отто Фольку, дается с помощью кватернионов в векторном изложении новый вывод KS-преобразования, действующего из четырехмерного пространства параметров в трехмерное физическое пространство. Использование кватернионов в векторном изложении (с использованием векторной записи для кватернионов со скалярной составляющей ξ ∈ R и векторной компонентой x ∈ R3) позволило автору дать каждому шагу в выводе KS-преобразования непосредственную геометрическую интерпретацию. В частности, KS-преобразование представлено как преобразование Леви-Чивита, сформулированное в повернутой системе координат. Отмечается, что здесь вполне естественным образом вступают (появляются) углы Эйлера. В качестве простого приложения приводятся явные формулы для KS-преобразования в случае эллиптического движения Кеплера в пространстве.

Во введеннии отмечается, что Volk и Waldvogel указали [60] на геометрическую интерпретацию KS-преобразования, включающую эйлеровы углы повернутой системы координат. Velte также отмечает, что “В этой статье Volk (1973) указал, что KS-преобразование уже встречается в письме Эйлера к Кристиану Гольдбаху в виде особого случая так называемой идентичности Эйлера”.

В завершение вывода KS-преобразования говорится: “Таким образом, мы нашли, что преобразование Кустаанхеймо и Штифеля можно вывести следующим образом. Во-первых, в плоскости движения вводится и формулируется посредством кватернионов преобразование Леви-Чивита. Введя теперь повернутую систему координат, получим преобразование Кустаанхеймо и Штифеля”.

Далее в разделе “Связь с углами Эйлера” отмечается, что “До сих пор мы описывали поворот системы координат с помощью кватерниона D, заданного в виде (9), содержащим половинный угол поворота и единичный вектор оси вращения. Но можно также выразить угол поворота и единичный вектор оси вращения с помощью трех эйлеровых углов (долготы узла, наклона плоскости движения и угла между линией узла и направлением единичного вектора el в плоскости движения)”. Поэтому компоненты кватерниона D являются функциями трех указанных углов Эйлера.

Отметим, что в статье [11], на наш взгляд, Velte фактически реализована прозрачная идея, заключающаяся в том, что если в системе координат X1X2X3 точка совершает плоское движение в плоскости X1X2 с координатами x1 и x2, то в системе координат $X_{1}^{'}X_{2}^{'}X_{3}^{'}$, повернутой на произвольный угол δ относительно исходной системы коооординат X1X2X3 вокруг оси, направление которой не совпадает с осью X3 и задается единичным вектором d, точка будет совершать пространственное (трехмерное) движение в системе координат $X_{1}^{'}X_{2}^{'}X_{3}^{'}$ с тремя координатами $x_{1}^{'}$, $x_{2}^{'}$ и $x_{3}^{'}$, которые могут быть выражены через четырехмерные переменные Кустаанхеймо–Штифеля uj. Именно такие формулы, связывающие новые координаты $x_{1}^{'}$, $x_{2}^{'}$ и $x_{3}^{'}$ точки с переменными Кустаанхеймо–Штифеля uj, в итоге получены из формул, связывающих декартовы координаты x1 и x2 в системе координат X1X2X3 с переменными Леви-Чивита ξ и η, в результате проделанных выкладок с использованием кватернионов (в том числе с использованием введенного постоянного произвольного кватерниона $D = \cos \frac{\delta }{2}$ + $\sin \frac{\delta }{2}{\mathbf{d}}$, описывающего поворот исходной системы координат (переход от системы координат X1X2X3 к системе координат $X_{1}^{'}X_{2}^{'}X_{3}^{'}$). Постоянство кватерниона $D$ в статье не оговорено, но это следует из использованного Velte термина “поворот” и проводимых выкладок. При этом полученные переменные Кустаанхеймо–Штифеля uj выражены им в явном виде через переменные Леви-Чивита ξ и η и компоненты dj кватерниона D. В четвертом разделе им также представлены формулы, с помощью которых компоненты dj могут быть найдены с помощью углов Эйлера. Также для эллиптического движения Кеплера в повернутой системе координат $X_{1}^{'}X_{2}^{'}X_{3}^{'}$ получены формулы для переменных Кустаанхеймо–Штифеля uj, выраженные через параметры орбиты, эксцентрическую аномалию E и компоненты dj кватерниона D.

Подчеркнем, что Velte-преобразование Кустаанхеймо и Штифеля (KS-преобразование) представлено как преобразование Леви-Чивита, сформулированное в повернутой на постоянный угол δ системе координат. Нами [29, 30] показано, как уже отмечалось в разд. 4, что регуляризующее KS-преобразование координат заключается в переходе от трехмерных декартовых координат центра масс второго (изучаемого) тела в инерциальной системе координат к новым четырехмерным переменным Кустаанхеймо–Штифеля ${{u}_{j}}$, которые являются нормированными определенным образом компонентами сопряженного кватерниона вращения ${\mathbf{\bar {\lambda }}} = {{\lambda }_{0}} - {{\lambda }_{1}}{\mathbf{i}}$${{\lambda }_{2}}{\mathbf{j}} - {{\lambda }_{3}}{\mathbf{k}}$, где ${{\lambda }_{j}}$ – параметры Эйлера. Эти параметры Эйлера ${{\lambda }_{j}}$ и кватернион вращения ${\mathbf{\lambda }}$ характеризуют ориентацию вращающейся (подвижной) системы координат η в инерциальной системе координат. Ось ${{\eta }_{1}}$ этой системы координат направлена вдоль радиус-вектора r центра масс второго тела, а ее начало находится в центре масс этого тела. Нормирующий множитель равен квадратному корню из расстояния r от центра масс второго тела до центра притяжения, поэтому кватернионные переменные u и λ связаны соотношением ${\mathbf{u}} = \sqrt r {\kern 1pt} {\mathbf{\bar {\lambda }}}$.

Подробный анализ работы Velte проводится также в связи с тем, что в следующем разделе излагаются полученные нами регулярные уравнения возмущенной пространственной задачи двух тел [8], в состав которых входят регулярные дифференциальные уравнения движения в модифицированных переменных Леви-Чивита второго порядка, описывающие движение второго (изучаемого) тела в идеальной системе координат и имеющие вид регулярных уравнений Леви-Чивита плоского движения, и регулярное кватернионное дифференциальное уравнение в четырехмерных параметрах Эйлера первого порядка, описывающих изменение ориентациии идеальной системы координат в инерциальной системе координат в процессе движения тела. Совместное использование модифицированных переменных Леви-Чивита и параметров Эйлера позволило нам расширить знаменитую регуляризацию Леви-Чивита плоского движения на пространственное движение с использованием идеальной системы координат.

В работе Vivarelli [12] KS-преобразование, введенное Кустаанхеймо и Штифелем в небесную механику, сформулировано в терминах гиперкомплексных чисел как произведение кватерниона и его антиинволюции. Следовательно, онo представляет особый морфизм вещественной алгебры кватернионов, имеющий трехмерное вещественное линейное подпространство, а также представляет собой естественное обобщение преобразования Леви-Чивита. Показано, что кватернионная матрица произведения приводит к KS-матрице; билинейное соотношение Кустаанхеймо–Штифеля и два тождества, которые играют центральную роль в KS-теории, легко выводятся с использованием кватернионов. Дается подходящее кватернионное калибровочное преобразование, которое приводит к хорошо известному расслоению четырехмерного пространства. В дополнение приводится несколько геометрических интерпретаций.

Отметим, что кватернион q, используемый автором [12], имеет вид

$q = {{u}_{1}} + {{u}_{2}}i + {{u}_{3}}j + {{u}_{4}}k,$
где i, j, k – векторные мнимые единицы Гамильтона.

Преобразование Кустаанхеймо–Штифеля представляется в следующей кватернионной форме:

$x = {{x}_{1}} + {{x}_{2}}i + {{x}_{3}}j = q{{q}_{*}},\quad {{q}_{*}} = k\bar {q}{{k}^{{ - 1}}} = {{u}_{1}} + {{u}_{2}}i + {{u}_{3}}j - {{u}_{4}}k,$
где ${{q}_{*}}$ – кватернион, антиинволюцивный кватерниону q.

Это представление отличается от кватернионного представления преобразования Кустаанхеймо–Штифеля

${{{\mathbf{r}}}_{\xi }} = {{\xi }_{1}}{\mathbf{i}} + {{\xi }_{2}}{\mathbf{j}} + {{\xi }_{3}}{\mathbf{k}} = {\mathbf{\bar {u}}} \circ {\mathbf{i}} \circ {\mathbf{u}};\quad {\mathbf{\bar {u}}} = {{u}_{0}} - {{u}_{1}}{\mathbf{i}} - {{u}_{2}}{\mathbf{j}} - {{u}_{3}}{\mathbf{k}},$
приведенного в [3] и используемого в наших работах.

В работах Vivarelli [13, 14] развивается геометрический и физический взгляд на векторное произведение двух кватернионов и рассматривается связь трех классических задач механики через гиперкомплексное KS-преобразование.

Исследование возмущенного кеплеровского движения проводится Штифелем и Шейфеле в [3] не только с использованием регулярных уравнений в осцилляторной форме и методов теории колебаний, но также с использованием регулярных уравнений в канонической форме, для чего ими разработана теория канонического KS-преобразования. Такой канонический подход к проблеме регуляризации, использующий KS-преобразование, развит в работах Лидова [61, 62] и Лидова и Ляховой [63] и широко используется в настоящее время. В более поздней работе [64] рассмотрено применение в теории регуляризации канонических уравнений задачи двух тел обобщенной KS-матрицы и связанных с ней преобразований.

Отметим также работу Шагова [15]. В ней получены дифференциальные уравнения движения космического аппарата в рамках возмущенной пространственной задачи двух тел, записанные в орбитальной системе координат и использующие для описания его движения в инерциальном пространстве кватернион поворота, нормированный посредством множителя, равного квадратному корню из модуля c вектора момента скорости спутника. В этих уравнениях в качестве независимой переменной использована переменная τ, связанная с временем t дифференциальным соотношением dτ = = (1/2)cr–2dt (здесь r – расстояние до центра Земли). Уравнения во “времени” τ линейны для невозмущенного кеплеровского движения спутника.

В работе Deprit, Elipe и Ferrer “Линеаризация: Лаплас против Штифеля” [16] матричная KS-теория Кустаанхеймо–Штифеля изложена в терминах кватернионов. Авторы рассматривают KS-преобразование независимо от его возможного применения к кеплеровским системам, уточняют несколько теорем, сформулированных Штифелем. В аннотации сказано: “с одной стороны, мы отказываемся от формализма матричной теории, чтобы продолжить исключительно в контексте алгебры кватернионов; с другой стороны, мы объясняем, как в иерархии гиперкомплексных систем как KS-преобразование, так и классическое проективное разложение возникают путем удвоения преобразования Леви-Чивита”.

Во введении этой работы отмечается: “Будучи свободными от вычислительных сервитутов (вычислительного рабства, вычислительной зависимости), мы даже поставили перед собой задачу переустановить всю KS-теорию в терминах кватернионов. К тому времени, как Штифель и Шейфеле завершили свою монографию, они осознали тесную связь между своим матричным формализмом и теорией кватернионов. Было предложено, чтобы они воспользовались этим; они отреагировали на предложение в чрезмерных терминах (Штифель и др., 1971, p. 286). Действительно ли они считают, что переход от матриц к кватернионам приведет к “провалу или, по крайней мере, к очень громоздкому формализму”?. Несмотря на суровые предсказания Штифеля, мы приняли вызов. Мы потерпели неудачу? Читатель – это наше жюри. Построение KS-преобразования как эманации альтернативной билинейной формы над алгеброй затрат кватернионов не сложнее, чем матричный формализм Штифеля и Шейфеле. Кроме того, мы находим награды в упражнениях: теоремы обострены, некоторые в значительной степени; доказательства сокращены; общий дизайн Штифеля значительно улучшился в отношении его глобального и внутреннего значения, не говоря уже о том, чтобы обеспечить стиль программирования для манипулирования кватернионами с помощью процессоров символов общего назначения.”

В работе Deprit, Elipe и Ferrer также рассматривается линеаризация уравнений движения, записанных в цилиндрической, сферической и орбитальной системах координат, посредством введения новых переменных и введения новой независимой переменной вместо времени t. В этой работе в качестве альтернативы KS-преобразованию было предложено собственное преобразование: DEF-преобразование. В отношении этого преобразования и известного BF-преобразования, также рассмотренного ими, отмечается следующее:

“Об альтернативах KS-преобразованию Штифель и Шейфеле (1971, p. 288) выпустили предупреждение, чуть ли не предписание: “авторы убеждены, что поиск других преобразований […] не очень многообещающий”. Многим читателям их предзнаменование создало впечатление, что KS-методика уникальна в совместной регуляризации, линеаризации и увеличении размерности для трехмерных кеплеровских систем. Факты отвергают иск. Кустанхаймо, Штифель и Шейфеле никогда не упоминали о решающем шаге, который Фок (1935, 1936) сделал в этом направлении тридцать пять лет назад, даже в их кратком упоминании того же шага, предпринятого, но независимо, Мозером (Moser 1970). Здесь у нас есть место для последней записи в конкурсе: BF-преобразование, предложенное Burdet (1969) для координатной части и завершенное Ferrfandiz (1986a,b, 1987, 1988) для моментальной части. Мы дополняем его нашим собственным преобразованием, DEF-преобразованием, которое, как мы утверждаем, одинаково хорошо выполняет все задачи KS-преобразования – линеаризацию, регуляризацию и каноничность – хотя, как мы склонны верить, более простым и более интуитивным способом (раздел 4.1). По общему признанию, конструкция включает в себя тяжелые алгебраические манипуляции, однако не более, чем в случае с KS- или BF-преобразованием. Кроме того, мы передаем эту работу процессору символов.”

Восьмимерное DEF-преобразование декартовых координат и проекций вектора скорости материальной точки имеет (в обозначениях авторов этого преобразования) следующий вид [16]:

$\begin{gathered} {\mathbf{x}} = {{u}_{0}}{\mathbf{u}} \\ {\mathbf{X}} = {{U}_{0}}{\mathbf{u}} + \frac{1}{{{{u}_{0}}}}\left( {{\mathbf{u}} \times {\mathbf{U}}} \right) \times {\mathbf{u}} = {{U}_{0}}{\mathbf{u}} + \frac{1}{{{{u}_{0}}}}\left( {\left( {{\mathbf{u}} \cdot {\mathbf{u}}} \right){\mathbf{U}} - \left( {{\mathbf{u}} \cdot {\mathbf{U}}} \right){\mathbf{u}}} \right), \\ \end{gathered} $
где x = r – радиус-вектор точки, X = v = ${\mathbf{\dot {r}}}$ – ее вектор скорости; u0, U0 и u, U – новые скалярные и трехмерные векторные переменные.

Независимая переменная (время t) заменяется на обобщенную истинную анома-лию  f, такую, что

$u_{0}^{2}df = {{\beta }^{2}}Qdt;\quad Q = \left\| {\mathbf{Q}} \right\|,\quad {\mathbf{Q}} = {\mathbf{u}} \times {\mathbf{U}} = \frac{1}{{{{\beta }^{2}}}}{\mathbf{x}} \times {\mathbf{X}},$
где β – параметр.

DEF-преобразование имеет много общего с BF-преобразованием (Burdet–Ferrandiz transformation), имеющим вид

$\begin{gathered} {\mathbf{x}} = {{u}_{0}}{\mathbf{u}} \\ {\mathbf{X}} = \frac{1}{{{{{\left\| {\mathbf{u}} \right\|}}^{2}}}}\left( {{{U}_{0}}{\mathbf{u}} + \frac{1}{{{{u}_{0}}}}\left( {{\mathbf{u}} \times {\mathbf{U}}} \right) \times {\mathbf{u}}} \right) = \left( {{{U}_{0}}{\mathbf{u}} + \frac{1}{{{{u}_{0}}}}\left( {\left( {{\mathbf{u}} \cdot {\mathbf{u}}} \right){\mathbf{U}} - \left( {{\mathbf{u}} \cdot {\mathbf{U}}} \right){\mathbf{u}}} \right)} \right) \\ \end{gathered} $

Оно отличается от него отсутствием множителя ${\text{1/}}{{\left\| {\mathbf{u}} \right\|}^{2}}$ и является более простым.

Уравнения движения материальной точки в ньютоновском гравитационном поле в DEF-переменных (при отсутствии возмущений) после перехода к обобщенной истинной аномалии  f в качестве независимой переменной принимают следующий вид [16]:

$\begin{gathered} \frac{{d{{u}_{0}}}}{{df}} = \frac{{u_{0}^{2}}}{Q}{{U}_{0}},\quad \frac{{d{{U}_{0}}}}{{df}} = \frac{Q}{{{{u}_{0}}}} - \frac{\mu }{{{{\beta }^{3}}Q}} \\ \frac{{d{\mathbf{u}}}}{{df}} = \frac{1}{Q}{\mathbf{Q}} \times {\mathbf{u}},\quad \frac{{d{\mathbf{U}}}}{{df}} = \frac{1}{Q}{\mathbf{Q}} \times {\mathbf{U}} - \frac{{u_{0}^{2}}}{{{{\beta }^{4}}}}\left( {2h + 3\frac{\mu }{{\beta {{u}_{0}}}}} \right){\mathbf{u}};\quad \frac{{dt}}{{df}} = \frac{{u_{0}^{2}}}{{{{\beta }^{2}}Q}}, \\ \end{gathered} $
где μ – гравитационная постоянная, Q = const – “угловой момент”, h – кеплеровская энергия (постоянная величина).

Дифференцируя третье уравнение приведенной системы уравнений по переменной  f, получаем [16] линейное дифференциальное уравнение второго порядка относительно векторной переменной u:

$\frac{{d{{{\mathbf{u}}}^{2}}}}{{d{{f}^{2}}}} + {\mathbf{u}} = 0$

Для новой скалярной переменной σ = Q2/(μu0), введенной вместо переменной u0, получаем, учитывая первое и второе уравнения вышеприведенной системы, следующее линейное дифференциальное уравнение относительно скалярной переменной σ [16]:

$\frac{{d{{\sigma }^{2}}}}{{d{{f}^{2}}}} + \sigma = \frac{1}{{{{\beta }^{3}}}}$

После получения двух последних уравнений ([16], p. 191) говорится “Это завершает линеаризацию кеплеровских систем в трех измерениях посредством DEF-преобразования”.

Никаких комментариев в отношении второго уравнения и основного четвертого динамического уравнения движения материальной точки (уравнения первого порядка для векторной переменной U) полученной в [16] системы уравнений не дается, хотя видно, что четвертое уравнение содержит в правой части векторное слагаемое

$\frac{{u_{0}^{2}}}{{{{\beta }^{4}}}}\left( {2h + 3\frac{\mu }{{\beta {{u}_{0}}}}} \right){\mathbf{u}},$
которое после подстановки в него аналитических решений u0 = u0f ) и u = uf ) (функций обобщенной истинной аномали f ), полученных в результате интегрирования приведенных последних двух линейных уравнений, принимает вид нестационарного выражения, явно зависящего от “времени” f.

Уравнение для переменной U принимает при этом вид векторного дифференциального неоднородного линейного уравнения, однородная часть которого имеет постоянные коффициенты, а неоднородная часть явно зависит от “времени” f  (обобщенной истинной аномалии f ), что не позволяет построить аналитическое решение этого уравнения в простой форме.

Кватернионные уравнения движения материальной точки в ньютоновском гравитационном поле в переменных Кустаанхеймо–Штифеля имеют вид [30] (в них четырехмерная кватернионная переменная u отлична от векторной переменной u DEF)

$\begin{gathered} \frac{{{{d}^{2}}{\mathbf{u}}}}{{d{{\tau }^{2}}}} - \frac{1}{2}h{\mathbf{u}} = 0;\quad \frac{{dt}}{{d\tau }} = r;\quad h = {\text{const}} \\ {\mathbf{u}} = {{u}_{0}} + {{u}_{1}}{\mathbf{i}} + {{u}_{2}}{\mathbf{j}} + {{u}_{3}}{\mathbf{k}},\quad r = \left| {\mathbf{r}} \right| = {\mathbf{\bar {u}}} \circ {\mathbf{u}} = {\mathbf{u}} \circ {\mathbf{\bar {u}}} = u_{0}^{2} + u_{1}^{2} + u_{2}^{2} + u_{3}^{2} \\ \end{gathered} $

Здесь первое кватернионное уравнение эквивалентно системе четырех независимых линейных однородных дифференциальных уравнений второго порядка относительно компонент uj  кватернионной переменной u (переменных Кустаанхеймо–Штифеля [3]):

$\frac{{{{d}^{2}}{{u}_{j}}}}{{d{{\tau }^{2}}}} - \frac{1}{2}h{{u}_{j}} = 0;\quad h = {\text{const,}}\quad j = 0,1,2,3$

Эти уравнения имеют постоянные коэффициенты, равные половинной постоянной кеплеровской энергии h материальной точки, взятой со знаком “минус”, и легко интегрируются в элементарных функциях в случае эллиптического кеплеровского движения, когда h < 0, или, для любого знака энергии h, в функциях Штумпфа. В случае эллиптического кеплеровского движения основное кватернионное уравнение эквивалентно уравнению движения четырехмерного одночастотного осциллятора. Уравнение для времени t может быть проинтегрировано отдельно от этих уравнений.

Радиус-вектора r и вектор скорости ${\mathbf{v}} = {\mathbf{\dot {r}}}$ материальной точки в инерциальной системе координат связаны с переменными u и $d{\mathbf{u}}{\text{/}}d\tau $ кватернионными соотношениями

${\mathbf{r}} = {\mathbf{\bar {u}}} \circ {\mathbf{i}} \circ {\mathbf{u}},\quad {\mathbf{v}} = \frac{2}{r}{\mathbf{\bar {u}}} \circ {\mathbf{i}} \circ \frac{{d{\mathbf{u}}}}{{d\tau }}$

На наш взгляд, из сравнения приведенных уравнений видно преимущество уравнений в KS-переменных перед уравнениями в DEF-переменных. Отметим, что из кватернионных уравнений возмущенного движения материальной точки в KS-переменных получаются [37] уравнения в кватернионных оскулирующих элементах (в кватернионных медленных переменных), удобные для исследования возмущенного эллиптического движения материальной точки.

Среди других ранних зарубежных работ отметим статьи Vrbik [17, 18], в которых продемонстрирована эффективность применения кватернионов к решению возмущенной задачи Кеплера. Так приведено доказательство формул [18] для построения пертурбативного решения возмущенной задачи Кеплера с использованием алгебры кватернионов в формулировке уравнений Кустанхеймо–Штифеля. Главным преимуществом такого подхода является удаление из соответствующего решения быстрых колебаний (в случае консервативных сил) и малых делителей (в случае зависящих от времени сил).

Среди более поздних работ отметим статьи Вальдфогеля [19, 20], имеющего совместные работы с Штифелем [65, 66]. Утверждается [20], что “кватернионы для регуляризации небесной механики – верный (правильный) путь” и что кватернионы “являются идеальным инструментом для описания и разработки теории пространственной регуляризации в небесной механике”.

Вальдвогель [20] пишет: “Данная статья подкрепляет это утверждение. Мы начнем с краткого введения в алгебру кватернионов, затем мы опишем процедуру регуляризации и ее следствия в элегантной форме. Кроме того, будет приведен альтернативный вывод теории движения Кеплера, основанный на регуляризации. Также мы рассмотрим регуляризацию ограниченной пространственной задачи трех тел, т.е. пространственное обобщение преобразования Бирхоффа. В завершение будет приведено описание возмущенного движения Кеплера в регуляризованных переменных”.

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

Укажем основные особенности кватернионного метода регуляризации Вальдфогеля [20], а именно, для регуляризации предлагается использовать “звездно сопряженный” кватернион (“star conjugate of the quaternion”)

${\mathbf{u}}{\kern 1pt} * = - k{\mathbf{\bar {u}}}k = {{u}_{0}} + i{{u}_{1}} + j{{u}_{2}} - k{{u}_{3}}$
$x = {{x}_{1}} + {{x}_{2}}i + {{x}_{3}}j = q{{q}_{*}},\quad {{q}_{*}} = k\bar {q}{{k}^{{ - 1}}} = {{u}_{1}} + {{u}_{2}}i + {{u}_{3}}j - {{u}_{4}}k,$
где ${\mathbf{u}} = {{u}_{0}} + i{{u}_{1}} + j{{u}_{2}} + k{{u}_{3}}$, а также отображение

${\mathbf{u}} \in {\text{U}} \to {\mathbf{x}} = {\mathbf{uu}}{\kern 1pt} *$

В этом отображении используется нетрадиционное представление трехмерного вектора x кватернионом x = x0 + ix1 + jx2 с нулевой k-компонентой (отметим, что специальный символ “$ \circ $” кватернионного произведения Вальдфогелем не используется). Такой кватернион x является формальным обобщением (наращиванием) комплексной переменной x = x0 + ix1, использованной Леви-Чивита в теории регуляризации уравнений плоского движения.

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

${\mathbf{x}} = {\mathbf{uu}}{\kern 1pt} * = - {\mathbf{u}}k{\mathbf{\bar {u}}}k$

В скалярной записи из последней формулы имеем

${{x}_{0}} = u_{0}^{2} - u_{1}^{2} - u_{2}^{2} + u_{3}^{2},\quad {{x}_{1}} = 2\left( {{{u}_{0}}{{u}_{1}} - {{u}_{2}}{{u}_{3}}} \right),\quad {{x}_{2}} = 2\left( {{{u}_{0}}{{u}_{2}} + {{u}_{1}}{{u}_{3}}} \right),$
“что в точности является KS-преобразованием в его классической форме или – до перестановки индексов – отображением Хопфа” [20].

Отметим, что кватернионы Вальдфогеля x, u и ${\mathbf{u}}{\kern 1pt} *$ совпадают с кватернионами Vivarelli [12] x, q и $q{\kern 1pt} *$ с точностью до обозначений индексов их компонент.

В классической теории кватернионов трехмерному вектору x ставится в соответствие кватернион ${\mathbf{x}} = {{x}_{1}}{\mathbf{i}} + {{x}_{2}}{\mathbf{j}} + {{x}_{3}}{\mathbf{k}}$ с нулевой скалярной частью. В работах автора статьи для регуляризации используются кватернионные переменные ${\mathbf{u}} = {{u}_{0}} + {{u}_{1}}{\mathbf{i}} + {{u}_{2}}{\mathbf{j}}$ + + u3k и ${\mathbf{u}}{\kern 1pt} * = u_{0}^{*} + u_{1}^{*}{\mathbf{i}} + u_{2}^{*}{\mathbf{j}} + u_{3}^{*}{\mathbf{k}}$, не совпадающия (по смыслу) с кватернионной переменной Вальдфогеля, и кватернион ${\mathbf{x}} = {{x}_{1}}{\mathbf{i}} + {{x}_{2}}{\mathbf{j}} + {{x}_{3}}{\mathbf{k}}$ с нулевой скалярной частью. В этих работах используется отображение

${\mathbf{r}} = {{x}_{1}}{\mathbf{i}} + {{x}_{2}}{\mathbf{j}} + {{x}_{3}}{\mathbf{k}} = {\mathbf{\bar {u}}} \circ {\mathbf{i}} \circ {\mathbf{u}},$
а также отображение

${\mathbf{r}} = {{x}_{1}}{\mathbf{i}} + {{x}_{2}}{\mathbf{j}} + {{x}_{3}}{\mathbf{k}} = {\mathbf{\bar {u}}}{\kern 1pt} * \circ \;{\mathbf{k}} \circ {\mathbf{u}}{\kern 1pt} *$

В скалярной записи первое из этих отображений дает в точности преобразование Кустаанхеймо–Штифеля

${{x}_{1}} = u_{0}^{2} + u_{1}^{2} - u_{2}^{2} - u_{3}^{2},\quad {{x}_{2}} = 2\left( {{{u}_{1}}{{u}_{2}} - {{u}_{0}}{{u}_{3}}} \right),\quad {{x}_{3}} = 2\left( {{{u}_{1}}{{u}_{3}} + {{u}_{0}}{{u}_{2}}} \right),$
которое по своей форме отличается от выше приведенного преобразования Вальдфогеля.

В скалярной записи второе из этих отображений дает преобразование координат

${{x}_{1}} = 2\left( {u_{1}^{*}u_{3}^{*} - u_{0}^{*}u_{2}^{*}} \right),\quad {{x}_{2}} = 2\left( {u_{2}^{*}u_{3}^{*} + u_{0}^{*}u_{1}^{*}} \right),\quad {{x}_{3}} = u_{0}^{{*2}} - u_{1}^{{*2}} - u_{2}^{{*2}} + u_{3}^{{*2}},$
отличное от преобразования Кустаанхеймо–Штифеля, и позволяет, как показано нами [41], получить регуляризованные кватернионные уравнения возмущенного движения искусственного спутника Земли в ее гравитационном поле в четырехмерных переменных $u_{j}^{*}$, которые обладают всеми достоинствами уравнений движения спутника в переменных Кустаанхеймо–Штифеля ${{u}_{j}}$ (также полученных нами), но имеют более простую и симметричную структуру в случае учета в потенциале гравитационного поля Земли не только его центральной (ньютоновской) составляющей, но и других составляющих: зональных, тессеральных и секториальных гармоник поля тяготения Земли.

Отметим полученное элегантное кватернионное представление [19, 20] пространственного отображения Биркгофа [67], используемого в теории регуляризации уравнений ограниченной задачи трех тел. Такое же преобразование, известное как отображение Жуковского, используется в аэродинамике для отображения поперечных разрезов аэродинамической поверхности до приблизительной круговой формы. В работе [20] это преобразование названо преобразованием Жуковского–Биркгофа (Joukowsky–Birkhoff mapping). Это представление приведено там в качестве дополнения к более ранним работам по теории регуляризации [65, 68].

Saha показано [21], как преобразование Кустанхеймо–Штифеля, записанное в кватернионной форме, можно интерпретировать как вращение в трех измерениях с использованием оси и угла вращения. Отметим, что такая интерпретация преобразования Кустанхеймо–Штифеля была дана значительно ранее в работах автора статьи [29, 30].

Отмечается [21], что “преобразование Кустанхеймо–Штифеля превращает гравитационную проблему двух тел в гармонический осциллятор при переходе к четырем измерениям. В дополнение к интересам математической физики, KS-преобразование оказалось очень полезным при моделировании N-тел, где оно помогает справляться с близкими встречами. И все же формализм остается несколько загадочным, а роль дополнительного измерения особенно загадочна. В этой статье показано, как базовое преобразование можно интерпретировать как вращение в трех измерениях. Например, если повернуть телескоп из зенита к выбранной звезде за один оборот, то можно представить ось вращения и угол как KS-преобразование звезды. Неединственность оси вращения кодирует дополнительное измерение. Эта геометрическая интерпретация становится очевидной при написании KS-преобразований в кватернионной форме, что также помогает получить краткие выражения для регуляризованных уравнений движения”.

В статье [22] регуляризация Кустанхеймо–Штифеля пространственной задачи Кеплера представлена в симплектическом и кватернионном подходах.

Рассматриваются [23] устойчивость и хаос в пространстве Кустаанхеймо–Штифеля, индуцированные расслоением Хопфа. Даны [24] полностью регулярные и универсальные решения проблемы относительного движения космического корабля в двух формах, полученных из регуляризаций KS и Сперлинга–Бурде (SB). В полученных решениях нет особенностей, и на их форму не влияет тип эталонной орбиты (круговая, эллиптическая, параболическая или гиперболическая). Кроме того, решения задачи даны в компактных тензорных выражениях и непосредственно относятся к вектору начального состояния корабля-лидера. Формулировки SB и KS вводят фиктивное время посредством преобразования Сундмана (Sundman transformation). Из-за использования альтернативной независимой переменной решения строятся на основе теории асинхронного относительного движения. Этот метод упрощает необходимые производные. Замкнутые выражения частных производных орбитального движения по начальному состоянию приведены явно.

Среди последних работ в области преобразований Леви-Чивита и Кустанхеймо–Штифеля отметим работы Breiter и Langner [2527]. Отмечается [27], что “KS-преобразование приобрело популярность в матрично-векторной формулировке Кустанхеймо и Штифеля (1965), но его гораздо проще интерпретировать и обобщать на языке кватернионной алгебры, очень тесно связанной с оригинальной спинорной формулировкой Кустанхеймо (1964)”. Излагается KS-преобразование в кватернионной форме, предлагается его обобщенное определение, развивается геометрическая интерпретация KS-переменных, предложенная в [21], рассматривается билинейная форма Кустанхеймо–Штифеля и ее обобщение. Однако отметим, что основная цель [27] – “вывести альтернативный набор переменных действие–угол, который не основан на понятии плоскости орбиты (таким образом, избегая сингулярностей, когда орбита вырождается в прямой отрезок), и проверить его на некоторой хорошо известной астрономической проблеме (на проблеме Лидова–Козаи)”.

В связи с этим говорится “Но тем, кто хочет извлечь выгоду из богатства канонического формализма, требуется набор переменных “действие–угол” регуляризованной задачи Кеплера. Первый шаг в этом направлении можно найти в монографии Штифеля и Шейфеле (1971), где симплектические полярные координаты вводятся для каждой отдельной степени свободы. Однако этот подход не учитывает вырождение проблемы и, следовательно, не подходит для методов возмущения на основе усреднения. Более того, не было предпринято никаких попыток связать этот набор с ограничением, известным как “билинейный инвариант”, эффективно сокращающим систему до трех степеней свободы. Обе проблемы были решены Чжао (2015), который предложил переменные “LCF” (предположительно названные в честь Леви-Чивита (1906) и Фейоса (2001)). В его подходе движение в KS-переменных рассматривается в колеблющейся “плоскости Леви-Чивита” (Деприт и др., 1994) как проблема с двумя степенями свободы. Третья степень свободы добавляется парой переменных действие–угол, ориентирующих плоскость. Избыточная четвертая степень скрыта в определении плоскости Леви-Чивита. Преобразованный кеплеровский гамильтониан зависит от одной переменной действия, два других действия тесно связаны с угловым моментом и его проекцией на полярную ось. Интересно, что результат идентичен “изоэнергетическим переменным”, найденным Леви-Чивита (1913) без регуляризации.”

Задача Лидова–Козаи анализируется [27] в терминах LKS-переменных, которые позволяют непосредственно исследовать устойчивость для всех равновесий, кроме круговой экваториальной и полярной радиальной орбит.

Предложен [28] альтернативный угловой подход к KS-преобразованию, предложена его интерпретация через симметрию.

Во многих работах проводилось сравнение точности численного решения регулярных уравнений в переменных Кустаанхеймо–Штифеля с решениями других регулярных уравнений, предложенных в этих работах. В частности, в [47] такое сравнение проводилось с уравнениями, предложенными в этой работе. В их состав входят уравнения в параметрах Эйлера, характеризующие ориентацию орбитальной системы координат. В качестве переменных используются расстояние, его производная по новой независимой переменной (а затем величины, являющиеся постоянными интегрирования в решениях уравнений невозмущенного движения для этих переменных), параметры Эйлера, а также момент количества орбитального движения.

В работах Fukushima (2005, 2007) [45, 46], ученого из японской Национальной астрономической обсерватории, показано, что KS-регуляризация приводит к очень эффективной схеме интегрирования уравнений орбитального движения, повышающей точность и скорость численного интегрирования. Это связано не только со структурой уравнений, но также с использованием нескольких методов, которые приносят важные преимущества численной схеме.

Приводится [46] численное сравнение четырех схем регуляризации трехмерной задачи двух тел в условиях возмущения: регуляризации SB, KS, BF и трехмерное расширение регуляризации LC. В аннотации статьи [46] сказано: КS и расширенная -регуляризация с масштабированием энергии Кеплера обеспечивают наилучшую экономическую эффективность при интеграции почти всех возмущенных задач двух тел. Подчеркнем, что Fukushima сравнивает семь схем, описанных в § 1 его статьи, а также “нерегуляризованную обработку, а именно прямое интегрирование в декартовых координатах”. Для каждой из этих формулировок им проводится тестовая интеграция Икара, охватывающая около 1 миллиона лет, и измеряется время его выполнения.

В качестве переменных состояния нового специального метода возмущений для задачи двух тел предложены [48] cемь пространственных элементов (в том числе элемент, обратный удвоенной полной энергии, и два первых интеграла (first integrals of the unperturbed motion)) и элемент времени. Новые элементы сохраняют нулевой эксцентриситет и наклон, а также отрицательные значения полной энергии. Предложены уравнения движения, записанные в орбитальной системе координат, для описания ориентации которой используются параметры Эйлера. Уравнения для параметров Эйлера записываются в скалярной и матричной формах. Сравниваются результаты численных решений предлагаемых уравнений и других уравнений, в том числе уравнений в KS-переменных. Из приводимых результатов численных решений можно сделать вывод, что для решаемых задач лучшую точность дают уравнения в KS-переменных и уравнения, предложенные авторами, причем решения новых уравнений имеют меньшие погрешности, чем решения уравнений в KS-переменных.

Показано [49], что специальные методы возмущения, основанные на регуляризованных формулировках, могут конкурировать и даже работать лучше, чем полуаналитические методы для изучения долгосрочного движения (порядка десятилетий) объектов, вращающихся вокруг Земли. Обращается внимание на то, что для такого рода применений формулировка Коуэлла (Cowell) никогда не используется из-за малых требуемых размеров шагов интегрирования, что вызывает сильное накопление ошибки округления и длительное время вычислений. Авторы этой работы разработали код Fortran, названный THALASSA, который включает метод Коуэлла, EDromo, регуляризацию KS, и набор регулярных элементов, которые были получены Штифелем и Шейфеле (1971, раздел 19) из уравнений в KS-переменных. Сложный числовой решатель, названный LSODAR (решатель для обыкновенных дифференциальных уравнений с автоматическим поиском корней), был включен для интеграции дифференциальных уравнений движения.

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

Представлен [50] новый метод вычисления орбит в возмущенной задаче двух тел: векторы положения и скорости движущегося объекта в декартовых координатах заменяются восемью орбитальными элементами, то есть константами невозмущенного движения. Два из них связаны с радиальным движением, следующие четыре – параметры Эйлера, задающие ориентацию промежуточной системы координат, эволюция которой отслеживает ориентацию плоскости орбиты и направление отсчета на ней. Полная энергия и элемент времени дополняют вектор состояния. Численные тесты включены для оценки эффективности предлагаемого специального метода возмущений. На примере орбитального движения двух комет показывается, что вычисления с использованием этого метода намного более точные и более быстрые, чем классические вычисления орбит с декартовыми координатами.

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

Рассмотрим предложенную нами кватернионную регуляризацию [8] уравнений возмущенной пространственной задачи двух тел с использованием переменных Леви-Чивита. Она основана на использовании двухмерных идеальных прямоугольных координат Ганзена, двухмерных переменных Леви-Чивита, четырехмерных параметров Эйлера (Родрига–Гамильтона) и кватерниона вращения Гамильтона, описывающих ориентацию идеальной системы координат [9], а также основана на использовании в качестве дополнительных переменных кеплеровской энергии и реального времени и на использовании новой независимой переменной Зундмана. Эти уравнения имеют не только хорошо известные достоинства уравнений Кустаанхеймо–Штифеля, но и имеют свои дополнительные достоинства.

6.1. Уравнения возмущенной пространственной задачи двух тел, записанные в орбитальной и идеальной системах координат с использованием кватернионного оскулирующего элемента орбитального движения. Эти уравнения имеют вид [8, 33, 34]

(6.1)
$\frac{{d{{{v}}_{1}}}}{{dt}} = \frac{{{{c}^{2}}}}{{{{r}^{3}}}} - \frac{\mu }{{{{r}^{2}}}} + {{p}_{{{\text{1orb}}}}},\quad \frac{{dr}}{{dt}} = {{{v}}_{1}},\quad \frac{{dc}}{{dt}} = r{{p}_{{{\text{2orb}}}}},\quad \frac{{d\varphi }}{{dt}} = \frac{c}{{{{r}^{2}}}}$
(6.2)
$\left. \begin{gathered} 2\frac{{d{{\Lambda }_{0}}}}{{dt}} = - {{\Omega }_{1}}{{\Lambda }_{1}} - {{\Omega }_{2}}{{\Lambda }_{2}},\quad 2\frac{{d{{\Lambda }_{1}}}}{{dt}} = {{\Omega }_{1}}{{\Lambda }_{0}} - {{\Omega }_{2}}{{\Lambda }_{3}}, \hfill \\ 2\frac{{d{{\Lambda }_{2}}}}{{dt}} = {{\Omega }_{2}}{{\Lambda }_{0}} + {{\Omega }_{1}}{{\Lambda }_{3}},\quad 2\frac{{d{{\Lambda }_{3}}}}{{dt}} = {{\Omega }_{2}}{{\Lambda }_{1}} - {{\Omega }_{1}}{{\Lambda }_{2}} \hfill \\ \end{gathered} \right\}$
(6.3)
${{\Omega }_{1}} = \frac{r}{c}{{p}_{{{\text{3orb}}}}}\cos \varphi ,\quad {{\Omega }_{2}} = \frac{r}{c}{{p}_{{{\text{3orb}}}}}\sin \varphi ,\quad {{\Omega }_{3}} = 0$
(6.4)
$2\frac{{d{\mathbf{\Lambda }}}}{{dt}} = {\mathbf{\Lambda }} \circ {{{\mathbf{\Omega }}}_{{{\text{id}}}}};\quad {{{\mathbf{\Omega }}}_{{{\text{id}}}}} = {{\Omega }_{1}}{\mathbf{i}} + {{\Omega }_{2}}{\mathbf{j}} = \frac{r}{c}{{p}_{{{\text{3orb}}}}}\left( {\cos \varphi {\mathbf{i}} + \sin \varphi {\mathbf{j}}} \right),\quad {{\Omega }_{3}} = 0$

Уравнение (6.4) является кватернионной записью скалярной системы дифференциальных уравнений (6.2) и (6.3) в параметрах Эйлера ${{\Lambda }_{j}}$, характеризующих ориентацию идеальной системы координат в инерциальной системе координат. Кватернионная переменная ${\mathbf{\Lambda }}$ в этом уравнении имеет смысл кватернионного оскулирующего элемента орбиты второго (изучаемого) тела: при равенстве нулю составляющей ${{p}_{{{\text{3orb}}}}}$ вектора p возмущающего ускорения центра масс второго тела, перпендикулярной плоскости, проходящей через радиус-вектор r и вектор v скорости этого тела, то есть, плоскости мгновенной орбиты второго тела, кватернион ${\mathbf{\Lambda }} = {\mathbf{const}}$.

Уравнения (6.1)(6.3) или (6.1) и (6.4), (6.3) – уравнения движения второго тела, записанные в двух вращающихся системах координат: в орбитальной системе координат ${{\eta }_{{{\text{orb}}}}}$ и в идеальной системе координат ${{\eta }_{{{\text{id}}}}}$. Уравнения (6.1) записаны в орбитальной системе координат, а уравнения (6.2) и (6.3) (или (6.4), (6.3)) – в идеальной системе координат.

Ось ${{\eta }_{{{\text{1orb}}}}}$ орбитальной системы координат ${{\eta }_{{{\text{orb}}}}}$, имеющей начало в центре масс второго тела, направлена по радиус-вектору r, а ее ось ${{\eta }_{{{\text{3orb}}}}}$ направлена вдоль вектора ${\mathbf{c}} = {\mathbf{r}} \times {\mathbf{v}}$ момента орбитальной скорости центра масс второго тела. Ось ${{\eta }_{{{\text{3id}}}}}$ идеальной системы координат ${{\eta }_{{{\text{id}}}}}$ параллельна оси ${{\eta }_{{{\text{3orb}}}}}$ орбитальной системы координат, а координатные оси ${{\eta }_{{{\text{1id}}}}}$ и ${{\eta }_{{{\text{2id}}}}}$ идеальной системы координат лежат в плоскости координатных осей ${{\eta }_{{{\text{1orb}}}}}$ и ${{\eta }_{{{\text{2orb}}}}}$ орбитальной системы координат и получаются из них поворотом вокруг оси ${{\eta }_{{{\text{3orb}}}}}$ на угол $\varphi $ по ходу часовой стрелки.

В уравнения (6.1)(6.3) или (6.1) и (6.4), (6.3) переменными являются расстояние r от центра масс второго тела до центра масс первого тела, производная dr/dt (проекция ${{{v}}_{1}}$ вектора скорости v центра масс второго тела на направление радиус-вектора r (на ось ${{\eta }_{{{\text{1orb}}}}}$ орбитальной системы координат)), модуль вектора момента орбитальной скорости второго тела $c = \left| {{\mathbf{r}} \times {\mathbf{v}}} \right|$, обобщенная истинная аномалия φ и параметры Эйлера (Родрига–Гамильтона) ${{\Lambda }_{j}}$j = 0, 1, 2, 3), характеризующие ориентацию идеальной системы координат ${{\eta }_{{{\text{id}}}}}$ в инерциальной системе координат ξ.

Фигурирующие в этих уравнениях величины ${{p}_{{k\operatorname{orb} }}}$ являются проекциями вектора p возмущающего ускорения центра масс второго тела на оси орбитальной системы координат ${{\eta }_{{{\text{orb}}}}}$.

Вектор Ω абсолютной угловой скорости идеальной системы координат ${{\eta }_{{{\text{id}}}}}$ параллелен радиус-вектору r центра масс второго тела и определяется формулой

${\mathbf{\Omega }} = \frac{{{{p}_{{{\text{3orb}}}}}}}{c}{\mathbf{r}};\quad c = \left| {{\mathbf{r}} \times {\mathbf{v}}} \right|$

Проекции ${{\Omega }_{i}}$ вектора Ω на оси идеальной системы координат ${{\eta }_{{{\text{id}}}}}$ запишем в виде

${{\Omega }_{1}} = \frac{r}{c}{{p}_{{{\text{3orb}}}}}\cos \varphi = \frac{1}{c}{{{\rm H}}_{1}}{{p}_{{{\text{3orb}}}}},\quad {{\Omega }_{2}} = \frac{r}{c}{{p}_{{{\text{3orb}}}}}\sin \varphi = \frac{1}{c}{{{\rm H}}_{2}}{{p}_{{{\text{3orb}}}}},\quad {{\Omega }_{3}} = 0,$
где ${{H}_{1}} = r\cos \varphi $, ${{H}_{2}} = r\sin \varphi $ – проекции радиус-вектора r на оси системы координат ${{\eta }_{{{\text{id}}}}}$.

Декартовые координаты ξk в инерциальной системе координат и проекции ${{{v}}_{k}} = {{{v}}_{{k\operatorname{orb} }}}$ вектора скорости центра масс второго тела на оси орбитальной системы координат находятся через указанные переменные по формулам

${{\xi }_{1}} = r\left( {\lambda _{0}^{2} + \lambda _{1}^{2} - \lambda _{2}^{2} - \lambda _{3}^{2}} \right),\quad {{\xi }_{2}} = 2r\left( {{{\lambda }_{1}}{{\lambda }_{2}} + {{\lambda }_{0}}{{\lambda }_{3}}} \right),\quad {{\xi }_{3}} = 2r\left( {{{\lambda }_{1}}{{\lambda }_{3}} - {{\lambda }_{0}}{{\lambda }_{2}}} \right)$
${{{v}}_{1}} = {{{v}}_{{{\text{1orb}}}}} = \frac{{dr}}{{dt}},\quad {{{v}}_{2}} = {{{v}}_{{{\text{2orb}}}}} = \frac{c}{r},\quad {{{v}}_{3}} = {{{v}}_{{{\text{3orb}}}}} = 0,$
где ${{\lambda }_{j}}$ – параметры Эйлера, характеризующих ориентацию орбитальной системы координат в инерциальной системе координат, которые предварительно находятся через переменные ${{\Lambda }_{j}}$ с использованием формул
${{\Lambda }_{0}} = {{\lambda }_{0}}{{\varphi }_{0}} + {{\lambda }_{3}}{{\varphi }_{3}},\quad {{\Lambda }_{1}} = {{\lambda }_{1}}{{\varphi }_{0}} - {{\lambda }_{2}}{{\varphi }_{3}},\quad {{\Lambda }_{2}} = {{\lambda }_{2}}{{\varphi }_{0}} + {{\lambda }_{1}}{{\varphi }_{3}},\quad {{\Lambda }_{3}} = - {{\lambda }_{0}}{{\varphi }_{3}} + {{\lambda }_{3}}{{\varphi }_{0}},$
в которых ${{\varphi }_{0}} = \cos \left( {\varphi {\text{/}}2} \right)$, ${{\varphi }_{3}} = \sin \left( {\varphi {\text{/}}2} \right)$, или с использованием кватернионных формул

${\mathbf{\Lambda }} = {\mathbf{\lambda }} \circ \left( {\cos \frac{\varphi }{2} - \sin \frac{\varphi }{2}{\mathbf{k}}} \right);\quad {\mathbf{\lambda }} = {\mathbf{\Lambda }} \circ \left( {\cos \frac{\varphi }{2} + \sin \frac{\varphi }{2}{\mathbf{k}}} \right)$

Проекции ${{{v}}_{{k\xi }}}$ вектора скорости центра масс второго тела на оси инерциальной системы координат ξ определяются кватернионными соотношениями

${{{\mathbf{v}}}_{\xi }} = {{{v}}_{{1\xi }}}{\mathbf{i}} + {{{v}}_{{2\xi }}}{\mathbf{j}} + {{{v}}_{{3\xi }}}{\mathbf{k}} = {{\dot {\xi }}_{1}}{\mathbf{i}} + {{\dot {\xi }}_{2}}{\mathbf{j}} + {{\dot {\xi }}_{3}}{\mathbf{k}} = {\mathbf{\lambda }} \circ {{{\mathbf{v}}}_{{\eta \operatorname{orb} }}} \circ {\mathbf{\bar {\lambda }}} = {\mathbf{\lambda }} \circ \left( {\dot {r}{\mathbf{i}} + \left( {c{\text{/}}r} \right){\mathbf{j}}} \right) \circ {\mathbf{\bar {\lambda }}},$
а связи проекций ${{p}_{{k\xi }}}$ вектора p возмущающего ускорения на оси инерциальной системы координат с их проекциями ${{p}_{{k\operatorname{orb} }}}$ на оси орбитальной системы координат определяются кватернионными соотношениями перепроектирования

$\begin{gathered} {{{\mathbf{p}}}_{\xi }} = {{p}_{{1\xi }}}{\mathbf{i}} + {{p}_{{2\xi }}}{\mathbf{j}} + {{p}_{{3\xi }}}{\mathbf{k}} = {\mathbf{\lambda }} \circ {{{\mathbf{p}}}_{{{\text{orb}}}}} \circ {\mathbf{\bar {\lambda }}} \\ {{{\mathbf{p}}}_{{{\text{orb}}}}} = {{p}_{{{\text{1orb}}}}}{\mathbf{i}} + {{p}_{{{\text{2orb}}}}}{\mathbf{j}} + {{p}_{{{\text{3orb}}}}}{\mathbf{k}} = {\mathbf{\bar {\lambda }}} \circ {{{\mathbf{p}}}_{\xi }} \circ {\mathbf{\lambda }} \\ \end{gathered} $

Отметим, что уравнения возмущенного кеплеровского движения, записанные в идеальной системе координат и состоящие из уравнений для полярных координат и уравнений для угловых переменных, описывающих ориентацию идеальной системы координат, были получены Andoyer (1923) [69] и Musen (1959) [70].

6.2. Уравнения возмущенной пространственной задачи двух тел, записанные в идеальной системе координат с использованием идеальных прямоугольных координат Ганзена и кватернионного оскулирующего элемента орбитального движения. Введем двухмерные идеальные прямоугольные координаты Ганзена ${{H}_{1}}$, ${{H}_{2}}$, являющиеся проекциями радиус-вектора r центра масс второго тела на оси идеальной системы координат ${{\eta }_{{{\text{id}}}}}$, связанные с переменными r и φ (полярными координатами) соотношениями

(6.5)
${{H}_{1}} = r\cos \varphi ,\quad {{H}_{2}} = r\sin \varphi ,\quad {{H}_{3}} = 0$

Дифференцируя эти соотношения дважды по времени и используя уравнения (6.1), получим вместо уравнений (6.1)–(6.3) следующие уравнения:

(6.6)
$\frac{{{{d}^{2}}{{H}_{1}}}}{{d{{t}^{2}}}} + \frac{\mu }{{{{r}^{3}}}}{{H}_{1}} = {{p}_{{{\text{1id}}}}},\quad \frac{{{{d}^{2}}{{H}_{2}}}}{{d{{t}^{2}}}} + \frac{\mu }{{{{r}^{3}}}}{{H}_{2}} = {{p}_{{{\text{2id}}}}},\quad {{H}_{3}} = 0,\quad {{r}^{2}} = H_{1}^{2} + H_{2}^{2}$
(6.7)
$\left. \begin{gathered} 2\frac{{d{{\Lambda }_{0}}}}{{dt}} = - {{\Omega }_{1}}{{\Lambda }_{1}} - {{\Omega }_{2}}{{\Lambda }_{2}},\quad 2\frac{{d{{\Lambda }_{1}}}}{{dt}} = {{\Omega }_{1}}{{\Lambda }_{0}} - {{\Omega }_{2}}{{\Lambda }_{3}} \hfill \\ 2\frac{{d{{\Lambda }_{2}}}}{{dt}} = {{\Omega }_{2}}{{\Lambda }_{0}} + {{\Omega }_{1}}{{\Lambda }_{3}},\quad 2\frac{{d{{\Lambda }_{3}}}}{{dt}} = {{\Omega }_{2}}{{\Lambda }_{1}} - {{\Omega }_{1}}{{\Lambda }_{2}} \hfill \\ \end{gathered} \right\}$
(6.8)
${{\Omega }_{1}} = \left( {{{H}_{1}}{\text{/}}c} \right){{p}_{{{\text{3id}}}}},\quad {{\Omega }_{2}} = \left( {{{H}_{2}}{\text{/}}c} \right){{p}_{{{\text{3id}}}}},\quad c = {{H}_{1}}{{\dot {H}}_{2}} - {{H}_{2}}{{\dot {H}}_{1}},\quad {{p}_{{{\text{3id}}}}} = {{p}_{{{\text{3orb}}}}}$

Уравнения (6.7) и (6.8) в кватернионной записи имеют вид

(6.9)
$2\frac{{d{\mathbf{\Lambda }}}}{{dt}} = {\mathbf{\Lambda }} \circ {{{\mathbf{\Omega }}}_{{{\text{id}}}}},$
где

$\begin{gathered} {\mathbf{\Lambda }} = {{\Lambda }_{0}} + {{\Lambda }_{1}}{\mathbf{i}} + {{\Lambda }_{2}}{\mathbf{j}} + {{\Lambda }_{3}}{\mathbf{k}} \\ {{{\mathbf{\Omega }}}_{{{\text{id}}}}} = {{\Omega }_{1}}{\mathbf{i}} + {{\Omega }_{2}}{\mathbf{j}} = \left( {{{p}_{{{\text{3id}}}}}{\text{/}}c} \right)\left( {{{H}_{1}}{\mathbf{i}} + {{H}_{2}}{\mathbf{j}}} \right) = {{p}_{{{\text{3id}}}}}{{\left( {{{H}_{1}}{{{\dot {H}}}_{2}} - {{H}_{2}}{{{\dot {H}}}_{1}}} \right)}^{{ - 1}}}\left( {{{H}_{1}}{\mathbf{i}} + {{H}_{2}}{\mathbf{j}}} \right) \\ \end{gathered} $

Декартовые координаты ξk центра масс второго тела в инерциальной системе координат и проекции ${{{v}}_{{k\xi }}}$ вектора абсолютной скорости v центра масс второго тела на оси инерциальной системы координат находятся через идеальные координаты Ганзена ${{H}_{1}}$, ${{H}_{2}}$ и их производные ${{\dot {H}}_{1}}$, ${{\dot {H}}_{2}}$ по кватернионным формулам

(6.10)
${{{\mathbf{r}}}_{{in}}} = {\mathbf{\Lambda }} \circ {{{\mathbf{r}}}_{{{\text{id}}}}} \circ {\mathbf{\bar {\Lambda }}},\quad {{{\mathbf{r}}}_{{in}}} = {{{\mathbf{r}}}_{\xi }} = {{\xi }_{1}}{\mathbf{i}} + {{\xi }_{2}}{\mathbf{j}} + {{\xi }_{3}}{\mathbf{k}},\quad {{{\mathbf{r}}}_{{{\text{id}}}}} = {{H}_{1}}{\mathbf{i}} + {{H}_{2}}{\mathbf{j}}$
(6.11)
${{{\mathbf{v}}}_{{in}}} = {\mathbf{\Lambda }} \circ {{{\mathbf{v}}}_{{{\text{id}}}}} \circ {\mathbf{\bar {\Lambda }}},\quad {{{\mathbf{v}}}_{{in}}} = {{{\mathbf{v}}}_{\xi }} = {{\dot {\xi }}_{1}}{\mathbf{i}} + {{\dot {\xi }}_{2}}{\mathbf{j}} + {{\dot {\xi }}_{3}}{\mathbf{k}},\quad {{{\mathbf{v}}}_{{{\text{id}}}}} = {{\dot {H}}_{1}}{\mathbf{i}} + {{\dot {H}}_{2}}{\mathbf{j}}$

Проекции ${{p}_{{k\operatorname{id} }}}$ вектора p возмущающего ускорения на оси идеальной системы координат ${{\eta }_{{{\text{id}}}}}$ связаны с его проекциями ${{p}_{{k\xi }}}$ на оси инерциальной системы координат ξ соотношениями перепроектирования

(6.12)
${{{\mathbf{p}}}_{{{\text{id}}}}} = {\mathbf{\bar {\Lambda }}} \circ {{{\mathbf{p}}}_{{in}}} \circ {\mathbf{\Lambda }},\quad {{{\mathbf{p}}}_{{{\text{id}}}}} = {{p}_{{{\text{1id}}}}}{\mathbf{i}} + {{p}_{{{\text{2id}}}}}{\mathbf{j}} + {{p}_{{{\text{3id}}}}}{\mathbf{k}},\quad {{{\mathbf{p}}}_{{in}}} = {{{\mathbf{p}}}_{\xi }} = {{p}_{{1\xi }}}{\mathbf{i}} + {{p}_{{2\xi }}}{\mathbf{j}} + {{p}_{{3\xi }}}{\mathbf{k}}$

Уравнения (6.6)(6.8) или (6.6), (6.9) являются уравнениями движения второго тела, записанными в идеальной системе координат ${{\eta }_{{{\text{id}}}}}$. В этих уравнениях переменными являются двухмерные идеальные прямоугольные координаты Ганзена ${{H}_{1}}$, ${{H}_{{\text{2}}}}$, их первые производные по времени ${{\dot {H}}_{1}}$, ${{\dot {H}}_{2}}$ и параметры Эйлера ${{\Lambda }_{j}}$j = 0, 1,2,3), характеризующие ориентацию идеальной системы координат ${{\eta }_{{{\text{id}}}}}$ в инерциальной системе координат ξ. Фигурирующие в этих уравнениях величины ${{p}_{{i\operatorname{id} }}}$ являются проекциями вектора p возмущающего ускорения центра масс второго тела на оси идеальной системы координат ${{\eta }_{{{\text{id}}}}}$ (${{p}_{{{\text{3id}}}}} = {{p}_{{{\text{3orb}}}}}$).

Отметим, что скалярные уравнения возмущенного кеплеровского движения в ганзеновских координатах и параметрах Эйлера, записанные в идеальной системе координат, ранее были получены в других формах и другими способами Deprit (1976) [9] и Брумбергом (1980) [71]. Уравнения (6.6)–(8.8) получены автором статьи [8] независимо и согласуются с уравнениями [9] и [71].

6.3. Регулярные кватернионные уравнения возмущенной пространственной задачи двух тел в переменных Леви-Чивита и параметрах Эйлера. Введем систему координат $O{{\eta }_{{{\text{id}}}}}$, оси которой параллельны осям идеальной системы координат ${{\eta }_{{{\text{id}}}}}$, а начало находится в центре масс O первого тела. Ориентация системы координат $O{{\eta }_{{{\text{orb}}}}}$, оси которой параллельны осям орбитальной системы координат ${{\eta }_{{{\text{orb}}}}}$, а начало находится в центре O, в идеальной системе координат ${{\eta }_{{{\text{id}}}}}$ (а также и в системе координат $O{{\eta }_{{{\text{id}}}}}$) характеризуется кватернионом поворота ${\mathbf{\Phi }}$, имеющим вид

${\mathbf{\Phi }} = \cos \frac{\varphi }{2} + \sin \frac{\varphi }{2}{\mathbf{k}}$

Компоненты ${{\Phi }_{j}}$j = 0, 1, 2, 3) этого кватерниона определяются соотношениями

${{\Phi }_{0}} = \cos \frac{\varphi }{2},\quad {{\Phi }_{1}} = {{\Phi }_{2}} = 0,\quad {{\Phi }_{3}} = \sin \frac{\varphi }{2}$
и являются параметрами Эйлера, характеризующими ориентацию орбитальной системы координат ${{\eta }_{{{\text{orb}}}}}$ в идеальной системе координат ${{\eta }_{{{\text{id}}}}}$.

Ранее введенные двухмерные идеальные прямоугольные координаты Ганзена ${{H}_{1}}$, ${{H}_{2}}$, определяемые соотношениями (6.5), являются декартовыми координатами центра масс второго тела в системе координат $O{{\eta }_{{{\text{id}}}}}$ и связаны с переменными r и ${{\Phi }_{j}}$ соотношениями

${{H}_{1}} = r\cos \varphi = r\left( {\Phi _{0}^{2} - \Phi _{3}^{2}} \right),\quad {{H}_{2}} = r\sin \varphi = 2r{{\Phi }_{0}}{{\Phi }_{3}}\quad ({{H}_{3}} = 0)$

Введем переменные Леви-Чивита

${{U}_{0}} = \sqrt r {{\Phi }_{0}},\quad {{U}_{3}} = - \sqrt r {{\Phi }_{3}},$
связанные с ганзеновскими координатами соотношениями

(6.13)
${{H}_{1}} = U_{0}^{2} - U_{3}^{2},\quad {{H}_{2}} = - 2{{U}_{0}}{{U}_{3}}$

Введение знака “–” в соотношении для переменной Леви-Чивита ${{U}_{3}} = - \sqrt r {{\Phi }_{3}}$ и, как следствие, для координаты Ганзена ${{H}_{2}} = - 2{{U}_{0}}{{U}_{3}}$ объясняется следующими соображениями. Для пространственного движения связи трехмерных декартовых координат ξk с четырехмерными переменными Кустаанхеймо–Штифеля ${{u}_{j}}$ имеют вид соотношений (3.4):

${{\xi }_{1}} = u_{0}^{2} + u_{1}^{2} - u_{2}^{2} - u_{3}^{2},\quad {{\xi }_{2}} = 2\left( {{{u}_{1}}{{u}_{2}} - {{u}_{0}}{{u}_{3}}} \right),\quad {{\xi }_{3}} = 2\left( {{{u}_{1}}{{u}_{3}} + {{u}_{0}}{{u}_{2}}} \right)$

Для плоского движения координата ${{\xi }_{3}} = 0$. Двухмерные декартовые координаты ξ1 и ξ2 в этом случае совпадают с координатами Ганзена ${{H}_{1}}$ и ${{H}_{2}}$ (напомним, что координата Ганзена ${{H}_{3}} = 0$): ${{\xi }_{1}} = {{H}_{1}}$ и ${{\xi }_{2}} = {{H}_{2}}$ (идеальная система координат в случае плоского движения может трактоваться как опорная система кординат ξ).

Нулевому значению координаты ${{\xi }_{3}}$ соответствуют, как видно из (3.4), нулевые значения либо переменных Кустаанхеймо–Штифеля ${{u}_{1}}$ и ${{u}_{2}}$, либо переменных ${{u}_{0}}$ и ${{u}_{3}}$. При этом переменными Леви-Чивита будут либо переменные ${{u}_{0}} = {{U}_{0}}$ и ${{u}_{3}} = {{U}_{3}}$, либо переменные ${{u}_{1}} = {{U}_{1}}$ и ${{u}_{2}} = {{U}_{2}}$.

Выбирая первый вариант выбора переменных ${{u}_{0}} = {{U}_{0}}$ и ${{u}_{3}} = {{U}_{3}}$ в качестве переменных Леви-Чивита, получим, учитывая, что ${{u}_{1}} = {{u}_{2}} = 0$

${{\xi }_{1}} = {{H}_{1}} = u_{0}^{2} - u_{3}^{2} = U_{0}^{2} - U_{3}^{2},\quad {{\xi }_{2}} = {{H}_{1}} = - 2{{u}_{0}}{{u}_{3}} = - 2{{U}_{0}}{{U}_{3}}$

Для второго варианта выбора переменных ${{u}_{1}} = {{U}_{1}}$ и ${{u}_{2}} = {{U}_{2}}$ в качестве переменных Леви-Чивита имеем соответственно

${{\xi }_{1}} = {{H}_{1}} = u_{1}^{2} - u_{2}^{2} = U_{1}^{2} - U_{2}^{2},\quad {{\xi }_{2}} = {{H}_{1}} = 2{{u}_{1}}{{u}_{2}} = 2{{U}_{1}}{{U}_{2}}$

Нами, таким образом, выше был выбран первый вариант определения переменных Леви-Чивита: ${{U}_{0}} = \sqrt r {{\Phi }_{0}}$, ${{U}_{3}} = - \sqrt r {{\Phi }_{3}}$, соответствующий описанию движения орбитальной системы координат относительно идеальной системы координат в параметрах Эйлера ${{\Phi }_{0}}$ и ${{\Phi }_{3}}$.

Отметим, что уравнениям Леви-Чивита (3.11)–(3.14) (они приведены в книге Штифеля и Шейфеле [3]) соответствует второй вариант выбора определения переменных Леви-Чивита.

Дополнительно дадим другое объяснение введения знака “–” в соотношении для переменной Леви-Чивита ${{U}_{3}} = - \sqrt r {{\Phi }_{3}}$ и, как следствие, для координаты Ганзена ${{H}_{2}} = - 2{{U}_{0}}{{U}_{3}}$. По аналогии с соотношениями

${{u}_{0}} = {{r}^{{1{\text{/}}2}}}{{\lambda }_{0}},\quad {{u}_{i}} = - {{r}^{{1{\text{/}}2}}}{{\lambda }_{i}};\quad i = 1,2,3,$
определяющими переменные Кустаанхеймо–Штифеля ${{u}_{j}}$ для пространственного движения через параметры Эйлера ${{\lambda }_{j}}$ и расстояние r до центра притяжения, для движения в идеальной системе координат (эквивалентного в этом случаю плоскому движению) будем иметь соотношения для переменных Леви-Чивита ${{U}_{0}}$ и ${{U}_{3}}$:
${{U}_{0}} = {{r}^{{1{\text{/}}2}}}{{\Phi }_{0}},\quad {{U}_{3}} = - {{r}^{{1{\text{/}}2}}}{{\Phi }_{3}}\quad \left( {{{U}_{1}} = {{U}_{2}} = 0} \right),$
которые определяют движение центра масс второго тела в идеальной системе координат. В этих соотношениях роль параметров Эйлера ${{\lambda }_{j}}$ играют выше введенные параметры Эйлера Φj:
${{\Phi }_{0}} = \cos \left( {\varphi {\text{/}}2} \right),\quad {{\Phi }_{1}} = {{\Phi }_{2}} = 0,\quad {{\Phi }_{3}} = \sin \left( {\varphi {\text{/}}2} \right),$
характеризующие ориентацию орбитальной системы координат ${{\eta }_{{{\text{orb}}}}}$ в идеальной системе координат. В соотношениях для переменных Леви-Чивита ${{U}_{0}}$ и ${{U}_{3}}$, как видно, присутствует знак “–” в соотношении для переменной ${{U}_{3}}$. Выражения для координат Ганзена принимают при этом вид
${{H}_{1}} = r\left( {\Phi _{0}^{2} - \Phi _{3}^{2}} \right) = U_{0}^{2} - U_{3}^{2},\quad {{H}_{2}} = 2r{{\Phi }_{0}}{{\Phi }_{3}} = - 2{{U}_{0}}{{U}_{3}},$
согласующийся с формулами (3.4) в случае плоского движения, когда ${{u}_{1}} = 0$, ${{u}_{2}} = 0$. В этих выражениях присутствует знак “–” в выражении для координаты ${{H}_{2}}$.

Проекции вектора скорости центра масс второго тела на оси системы координат O${{\eta }_{{{\text{id}}}}}$ (а также на оси идеальной системы координат ${{\eta }_{{{\text{id}}}}}$) связаны с производными по времени от переменных Леви-Чивита соотношениями

(6.14)
$\begin{gathered} {{{v}}_{{{\text{1id}}}}} = {{{v}}_{1}}\cos \varphi - {{{v}}_{2}}\sin \varphi = {{{\dot {H}}}_{1}} = 2\left( {{{U}_{0}}{{{\dot {U}}}_{0}} - {{U}_{3}}{{{\dot {U}}}_{3}}} \right) \\ {{{v}}_{{{\text{2id}}}}} = {{{v}}_{1}}\sin \varphi + {{{v}}_{2}}\cos \varphi = {{{\dot {H}}}_{2}} = - 2\left( {{{U}_{3}}{{{\dot {U}}}_{0}} + {{U}_{0}}{{{\dot {U}}}_{3}}} \right) \\ {{{v}}_{{{\text{3id}}}}} = {{{v}}_{3}} = {{{\dot {H}}}_{3}} = 0,\quad {{{v}}_{k}} = {{{v}}_{{k\operatorname{orb} }}} \\ \end{gathered} $

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

(6.15)
$\begin{gathered} h = (1{\text{/}}2){{{v}}^{{\text{2}}}} - \mu {\text{/}}r = (1{\text{/}}2)\left( {\dot {H}_{1}^{2} + \dot {H}_{2}^{2}} \right) - \mu {{\left( {H_{1}^{2} + H_{2}^{2}} \right)}^{{ - 1/2}}} \\ c = \left| {{\mathbf{r}} \times {\mathbf{v}}} \right| = {{H}_{1}}{{{\dot {H}}}_{2}} - {{H}_{2}}{{{\dot {H}}}_{1}}, \\ \end{gathered} $
в новых переменных Леви-Чивиты ${{U}_{0}}$ и ${{U}_{3}}$ принимают вид

(6.16)
$h = 2r\left( {\dot {U}_{0}^{2} + \dot {U}_{3}^{2}} \right) - \mu {\text{/}}r,\quad r = U_{0}^{2} + U_{3}^{2},\quad c = 2\left( {U_{0}^{2} + U_{3}^{2}} \right)\left( {{{U}_{3}}{{{\dot {U}}}_{0}} - {{U}_{0}}{{{\dot {U}}}_{3}}} \right)$

Перейдем в уравнениях движения второго тела (6.6)–(6.8), записанным в идеальной системе координат с использованием идеальных прямоугольных координат Ганзена ${{H}_{1}}$, ${{H}_{2}}$, к переменным Леви-Чивита ${{U}_{0}}$, ${{U}_{3}}$ по формулам (6.13), (6.14) и к новой независимой переменной τ в соответствии с дифференциальным соотношением Зундмана dt = rdτ. Введем также в качестве дополнительной переменной кеплеровскую энергию h, определяемую соотношениями (6.15), (6.16) и удовлетворяющую дифференциальному уравнению $dh{\text{/}}dt = {\mathbf{p}} \cdot {\mathbf{v}}$. В итоге получим следующие скалярные регулярные дифференциальные уравнения движения второго тела (регулярные уравнения возмущенной пространственной задачи двух тел) [8]:

(6.17)
$\frac{{{{d}^{2}}{{U}_{0}}}}{{d{{\tau }^{2}}}} - \frac{1}{2}h{{U}_{0}} = \frac{1}{2}r{{Q}_{0}},\quad \frac{{{{d}^{2}}{{U}_{3}}}}{{d{{\tau }^{2}}}} - \frac{1}{2}h{{U}_{3}} = \frac{1}{2}r{{Q}_{3}}$
(6.18)
$\frac{{dh}}{{d\tau }} = 2\left( {{{Q}_{0}}\frac{{d{{U}_{0}}}}{{d\tau }} + {{Q}_{3}}\frac{{d{{U}_{3}}}}{{d\tau }}} \right)$
(6.19)
$\left. \begin{gathered} 2\frac{{d{{\Lambda }_{0}}}}{{d\tau }} = - r\left( {{{\Omega }_{1}}{{\Lambda }_{1}} + {{\Omega }_{2}}{{\Lambda }_{2}}} \right),\quad 2\frac{{d{{\Lambda }_{1}}}}{{d\tau }} = r\left( {{{\Omega }_{1}}{{\Lambda }_{0}} - {{\Omega }_{2}}{{\Lambda }_{3}}} \right) \\ 2\frac{{d{{\Lambda }_{2}}}}{{d\tau }} = r\left( {{{\Omega }_{2}}{{\Lambda }_{0}} + {{\Omega }_{1}}{{\Lambda }_{3}}} \right),\quad 2\frac{{d{{\Lambda }_{3}}}}{{d\tau }} = r\left( {{{\Omega }_{2}}{{\Lambda }_{1}} - {{\Omega }_{1}}{{\Lambda }_{2}}} \right) \\ \end{gathered} \right\}$
(6.20)
$\frac{{dt}}{{d\tau }} = r$

Здесь

$r = \left| {\mathbf{r}} \right| = U_{0}^{2} + U_{3}^{2},\quad {{Q}_{0}} = {{U}_{0}}{{p}_{{{\text{1id}}}}} - {{U}_{3}}{{p}_{{{\text{2id}}}}},\quad {{Q}_{3}} = - {{U}_{3}}{{p}_{{{\text{1id}}}}} - {{U}_{0}}{{p}_{{{\text{2id}}}}}$
${{\Omega }_{1}} = \frac{1}{c}\left( {U_{0}^{2} - U_{3}^{2}} \right){{p}_{{{\text{3id}}}}},\quad {{\Omega }_{2}} = - \frac{2}{c}{{U}_{0}}{{U}_{3}}{{p}_{{{\text{3id}}}}};\quad c = 2\left( {{{U}_{3}}\frac{{d{{U}_{0}}}}{{d\tau }} - {{U}_{0}}\frac{{d{{U}_{3}}}}{{d\tau }}} \right),\quad {{p}_{{{\text{3id}}}}} = {{p}_{{{\text{3orb}}}}}$

Эти уравнения необходимо дополнить кватернионными соотношениями перепроектирования (6.10)–(6.12).

Отметим, что регулярные уравнения (6.17)(6.20) содержат, как подсистему, уравнения (6.17), (6.18), (6.20), имеющие вид регулярных уравнений Леви-Чивита плоской задачи двух тел. Эта подсистема уравнений совпадают с уравнениями Леви-Чивита плоской задачи двух тел (3.11)–(3.14), если в них формально ввести переобозначения ${{U}_{0}} = {{u}_{1}}$, ${{U}_{3}} = - {{u}_{2}}$ и ${{Q}_{0}} = q$, ${{Q}_{3}} = - {{q}_{2}}$.

Скалярные регулярные уравнения движения второго тела (6.17)–(6.20) (регулярные уравнения возмущенной пространственной задачи двух тел) в кватернионной записи имеют следующий вид [8]:

(6.21)
$\frac{{{{d}^{2}}{\mathbf{U}}}}{{d{{\tau }^{2}}}} - \frac{1}{2}h{\mathbf{U}} = \frac{1}{2}r{\mathbf{Q}},\quad \frac{{dh}}{{d\tau }} = 2\operatorname{scal} \left( {\frac{{d{\mathbf{\bar {U}}}}}{{d\tau }} \circ {\mathbf{Q}}} \right)$
(6.22)
$2\frac{{d{\mathbf{\Lambda }}}}{{d\tau }} = r{\mathbf{\Lambda }} \circ {{{\mathbf{\Omega }}}_{{{\text{id}}}}} = \frac{r}{c}{{p}_{{{\text{3id}}}}}{\mathbf{\Lambda }} \circ \left( {\left( {U_{0}^{2} - U_{3}^{2}} \right){\mathbf{i}} - 2{{U}_{0}}{{U}_{3}}{\mathbf{j}}} \right);\quad \frac{{dt}}{{d\tau }} = r$

Здесь

${\mathbf{U}} = {{U}_{0}} + {{U}_{3}}{\mathbf{k}},\quad {\mathbf{\bar {U}}} = {{U}_{0}} - {{U}_{3}}{\mathbf{k}},\quad {\mathbf{\Lambda }} = {{\Lambda }_{0}} + {{\Lambda }_{1}}{\mathbf{i}} + {{\Lambda }_{2}}{\mathbf{j}} + {{\Lambda }_{3}}{\mathbf{k}}$
$r = \left| {\mathbf{r}} \right| = \left\| {\mathbf{U}} \right\| = {\mathbf{U}} \circ {\mathbf{\bar {U}}} = {\mathbf{\bar {U}}} \circ {\mathbf{U}} = U_{0}^{2} + U_{3}^{2},\quad c = 2\left( {{{U}_{3}}\frac{{d{{U}_{0}}}}{{d\tau }} - {{U}_{0}}\frac{{d{{U}_{3}}}}{{d\tau }}} \right)$
${\mathbf{Q}} = - {\mathbf{i}} \circ {\mathbf{U}} \circ {{{\mathbf{P}}}_{{{\text{id}}}}},\quad {{{\mathbf{P}}}_{{{\text{id}}}}} = {{p}_{{{\text{1id}}}}}{\mathbf{i}} + {{p}_{{{\text{2id}}}}}{\mathbf{j}}$

Отображения радиус-вектора r и вектора скорости v центра масс второго тела на оси идеальной и инерциальной систем координат определяются соотношениями

${{{\mathbf{r}}}_{{{\text{id}}}}} = {{{\rm H}}_{1}}{\mathbf{i}} + {{{\rm H}}_{2}}{\mathbf{j}} = {\mathbf{\bar {U}}} \circ {\mathbf{i}} \circ {\mathbf{U}},\quad {{{\mathbf{r}}}_{{in}}} = {\mathbf{\Lambda }} \circ {{{\mathbf{r}}}_{{{\text{id}}}}} \circ {\mathbf{\bar {\Lambda }}},\quad {{{\mathbf{r}}}_{{in}}} = {{{\mathbf{r}}}_{\xi }} = {{\xi }_{1}}{\mathbf{i}} + {{\xi }_{2}}{\mathbf{j}} + {{\xi }_{3}}{\mathbf{k}}$
${{{\mathbf{v}}}_{{{\text{id}}}}} = \frac{{d{{{\mathbf{r}}}_{{{\text{id}}}}}}}{{dt}} = \frac{{d{{{\rm H}}_{1}}}}{{dt}}{\mathbf{i}} + \frac{{d{{{\rm H}}_{2}}}}{{dt}}{\mathbf{k}} = 2{\mathbf{\bar {U}}} \circ {\mathbf{i}} \circ \frac{{d{\mathbf{U}}}}{{dt}} = \frac{2}{r}{\mathbf{\bar {U}}} \circ {\mathbf{i}} \circ \frac{{d{\mathbf{U}}}}{{d\tau }}$
${{{\mathbf{v}}}_{{in}}} = {\mathbf{\Lambda }} \circ {{{\mathbf{v}}}_{{{\text{id}}}}} \circ {\mathbf{\bar {\Lambda }}},\quad {{{\mathbf{v}}}_{{in}}} = {{{\mathbf{v}}}_{\xi }} = {{\dot {\xi }}_{1}}{\mathbf{i}} + {{\dot {\xi }}_{2}}{\mathbf{j}} + {{\dot {\xi }}_{3}}{\mathbf{k}}$

Уравнения (6.17)(6.20) или (6.21) и (6.22) являются регулярными уравнениями возмущенной пространственной задачи двух тел, построенными с использованием идеальных прямоугольных координат Ганзена. В скалярных уравнениях (6.17)(6.20) регулярными переменными являются переменные Леви-Чивита U0 и U3, описывающие движение центра масс второго тела в идеальной системе координат O${{\eta }_{{{\text{id}}}}}$, кеплеровская энергия h, время t и параметры Эйлера ${{\Lambda }_{j}}$j = 0, 1, 2, 3), характеризующие ориентацию идеальной системы координат в инерциальной системе координат. В регулярных кватернионных уравнениях (6.21) и (6.22) переменными являются двухмерный кватернион U, описывающий движение центра масс второго тела в идеальной системе координат O${{\eta }_{{{\text{id}}}}}$, кеплеровская энергия h, время t и четырехмерный кватернион (кватернионный оскулирующий элемент) ${\mathbf{\Lambda }}$, характеризующий ориентацию идеальной системы координат в инерциальной системе координат. Фигурирующие в этих уравнениях величины ${{p}_{{{\text{1id}}}}}$, ${{p}_{{{\text{2id}}}}}$ являются проекциями вектора p возмущающего ускорения центра масс второго тела на оси идеальной системы координат. Эти величины находятся через проекции вектора p на оси инерциальной системы координат с помощью выше приведенных соотношений перепроектирования.

Регулярные уравнения (6.17)(6.20) или (6.21) и (6.22) возмущенной пространственной задачи двух тел образуют систему нелинейных нестационарных дифференциальных уравнений десятого порядка (такую же размерность имеют регулярные уравнения Кустаанхеймо–Штифеля) и обладают всеми достоинствами уравнений Кустаанхеймо–Штифеля:

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

$\frac{{{{d}^{2}}{{U}_{i}}}}{{d{{\tau }^{2}}}} - \frac{1}{2}h{{U}_{i}} = 0,\quad i = 0,3;\quad h = {\text{const,}}\quad {{\Lambda }_{j}} = {\text{const,}}\quad j = 0,1,2,3$
или кватернионный вид
$\frac{{{{d}^{2}}{\mathbf{U}}}}{{d{{\tau }^{2}}}} - \frac{1}{2}h{\mathbf{U}} = 0;\quad h = {\text{const,}}\quad \Lambda = {\text{ }}{\mathbf{const}}$
(для эллиптического кеплеровского движения, когда кеплеровская энергия h < 0, эти уравнения эквивалентны уравнениям движения двухмерного одночастотного гармонического осциллятора, квадрат частоты которого равен половине кеплеровской энергии, взятой со знаком “минус”);

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

− близки к линейным уравнениям для возмущенных кеплеровских движений;

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

Вместе с тем эти регулярные уравнения имеют существенные отличия:

1) для невозмущенного эллиптического кеплеровского движения изучаемого тела они эквивалентны уравнениям движения не четырехмерного одночастотного гармонического осциллятора, как в случае Кустаанхеймо–Штифеля, а уравнениям движения двухмерного одночастотного гармонического осциллятора, поскольку кватернион ориентации идеальной системы координат, в которой записаны эти уравнения орбитального движения, в этом случае остается постоянным;

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

Отметим, однако, что эти уравнения не пригодны для исследования прямолинейных орбит, когда модуль вектора момента орбитальной скорости второго тела $c = \left| {{\mathbf{r}} \times {\mathbf{v}}} \right|$ обращается в ноль, поскольку кватернионное дифференциальное уравнение ориентации идеальной системы координат в этом случае вырождается (из-за присутствия в этом уравнении знаменателя c).

От этого недостатка уравнений (6.17)–(6.20) или (6.21) и (6.22) можно избавиться, переходя в них от независимой переменной τ к новой независимой переменной $\tau {\kern 1pt} *$ в соответствии с дифференциальным соотношением $d\tau = cd\tau {\kern 1pt} *$ и дополняя полученные уравнения дифференциальным уравнением для переменной c. В итоге получаем уравнения, не имеющие указанной особенности:

$\frac{{d{\mathbf{U}}}}{{d{{\tau }}{\kern 1pt} *}} = с{\mathbf{S}},\quad \frac{{d{\mathbf{S}}}}{{d{{\tau }}{\kern 1pt} *}} - \frac{1}{2}сh{\mathbf{U}} = \frac{1}{2}сr{\mathbf{Q}}$
$\frac{{dc}}{{d{{\tau }}{\kern 1pt} *}} = сr\left[ {2{{p}_{{{\text{1id}}}}}{{U}_{0}}{{U}_{3}} + {{p}_{{{\text{2id}}}}}\left( {U_{0}^{2} - U_{3}^{2}} \right)} \right]$
$\frac{{dh}}{{d{{\tau }}{\kern 1pt} *}} = 2\operatorname{scal} \left( {\frac{{d{\mathbf{\bar {U}}}}}{{d{{\tau }}{\kern 1pt} *}} \circ {\mathbf{Q}}} \right)$
$2\frac{{d{\mathbf{\Lambda }}}}{{d{{\tau }}{\kern 1pt} *}} = сr{\mathbf{\Lambda }} \circ {{{\mathbf{\Omega }}}_{{{\text{id}}}}} = r{{p}_{{{\text{3id}}}}}{\mathbf{\Lambda }} \circ \left( {\left( {U_{0}^{2} - U_{3}^{2}} \right){\mathbf{i}} - 2{{U}_{0}}{{U}_{3}}{\mathbf{j}}} \right)$
$\frac{{dt}}{{d{{\tau }}{\kern 1pt} *}} = сr$
${\mathbf{Q}} = - {\mathbf{i}} \circ {\mathbf{U}} \circ {{{\mathbf{P}}}_{{{\text{id}}}}},\quad {{{\mathbf{P}}}_{{{\text{id}}}}} = {{p}_{{{\text{1id}}}}}{\mathbf{i}} + {{p}_{{{\text{2id}}}}}{\mathbf{j}}$
${{{\mathbf{p}}}_{{id}}} = {\mathbf{\bar {\Lambda }}} \circ {{{\mathbf{p}}}_{{in}}} \circ {\mathbf{\Lambda }},\quad {{{\mathbf{p}}}_{{{\text{id}}}}} = {{p}_{{{\text{1id}}}}}{\mathbf{i}} + {{p}_{{{\text{2id}}}}}{\mathbf{j}} + {{p}_{{{\text{3id}}}}}{\mathbf{k}},\quad {{{\mathbf{p}}}_{{in}}} = {{{\mathbf{p}}}_{\xi }} = {{p}_{{1\xi }}}{\mathbf{i}} + {{p}_{{2\xi }}}{\mathbf{j}} + {{p}_{{3\xi }}}{\mathbf{k}}$

В этих регулярных кватернионных уравнениях переменными, по-прежнему, являются двухмерный кватернион U, кеплеровская энергия h, время t, четырехмерный кватернион ${\mathbf{\Lambda }}$, а также двухмерный кватернион S и модуль c вектора момента орбитальной скорости. Эти уравнения сложнее регулярных кватернионных уравнений (6.21) и (6.22).

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

В заключение отметим, что Брумбергом [71] описано применение параметров Эйлера к выводу уравнений возмущенного движения пространственной задачи двух тел в ганзеновских координатах и указано на возможность дальнейшего преобразования полученных им уравнений возмущенного движения с использованием параболических координат Леви-Чивита.

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

Изложен предложенный автором новый метод регуляризации дифференциальных уравнений возмущенной пространственной задачи двух тел, основанный на использовании двухмерных идеальных прямоугольных координат Ганзена и двухмерных переменных Леви-Чивита, описывающих движение второго (рассматриваемого) тела в идеальной системе координат, в которой уравнения пространственного движения принимают вид уравнений плоского движения, а также основанный на использовании четырехмерных параметров Эйлера (Родрига–Гамильтона) и кватерниона Гамильтона, характеризующих ориентацию идеальной системы координат в инерциальной системе координат.

Изложены полученные с помощью этого метода регулярные скалярные и кватернионные уравнения возмущенной пространственной задачи двух тел. В скалярных уравнениях переменными являются переменные Леви-Чивита, кеплеровская энергия, время и указанные параметры Эйлера. В кватернионных уравнениях переменными являются двухмерный кватернион, описывающий движение центра масс второго тела в идеальной системе координат (его компоненты – переменные Леви-Чивита), кеплеровская энергия, время и четырехмернный кватернион ориентации идеальной системы координат (его компоненты – параметры Эйлера). Используемые в качестве переменных параметры Эйлера и кватернион ориентации идеальной системы координат являются скалярными и кватернионным оскулирующими элементами орбиты изучаемого (второго) тела (медленно изменяющимися переменными).

Изложенные регулярные уравнения возмущенной пространственной задачи двух тел, в которых совместно используются переменные Леви-Чивита и параметры Эйлера, образуют в общем случае систему нелинейных нестационарных дифференциальных уравнений десятого порядка (такую же размерность имеют и регулярные уравнения Кустаанхеймо–Штифеля). Они, также как и уравнения Кустаанхеймо–Штифеля, регулярны в центре притяжения (в отличие от нерегулярных ньютоновских уравнений), линейны для невозмущенных кеплеровских движений (в отличие от существенно нелинейных ньютоновских уравнений для этих движений); позволяют выработать единый подход к изучению всех трех типов кеплеровского движения (эллиптического, гиперболического, параболического) с использованием функций Штумпфа, близки к линейным уравнениям для возмущенных кеплеровских движений, позволяют представить правые части дифференциальных уравнений движения небесных и космических тел в полиномиальной форме, удобной для их численного решения.

Однако эти регулярные уравнения имеют следующие существенные отличия от регулярных уравнений в переменных Кустаанхеймо–Штифеля:

− для невозмущенного эллиптического кеплеровского движения изучаемого тела они эквивалентны уравнениям движения не четырехмерного одночастотного гармонического осциллятора, как в случае Кустаанхеймо–Штифеля, а уравнениям движения двухмерного одночастотного гармонического осциллятора, так как для этого случая движения тела кватернион ориентации идеальной системы координат, в которой записаны эти уравнения движения, и его компоненты (параметры Эйлера) остаются постоянными. Поэтому скалярные и кватернионное дифференциальные уравнения в параметрах Эйлера, описывающие изменение ориентации идеальной системы координат, выпадают в этом случае из рассмотрения;

− для возмущенного движения изучаемого тела кватернион ориентации идеальной системы координат является кватернионным оскулирующим элементом (т.е. медленно изменяющейся кватернионной переменной), а параметры Эйлера (компоненты этого кватерниона) – скалярными оскулирующими элементами, что также является полезным свойством этих уравнений, позволяющим эффективно использовать методы нелинейной механики.

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

Во введении было отмечено, что Леви-Чивита (1920) в отношении своих попыток обобщить предложенную им знаменитую регуляризацию уравнений плоской задачи двух тел на пространственную задачу позже признал их неудачу. Штифель и Шейфеле в своей книге (1971) также отмечали, что Леви-Чивита приложил много усилий, чтобы найти обобщение своего метода регуляризации дифференциальных уравнений плоского движения в задаче двух тел на общую пространственную задачу двух тел, но безуспешно. Aarseth и Zare (1974), а также Aarseth (2003) отмечали, что из-за фундаментальных трудностей, первоначально разъясненных Хопфом и Гурвицем, невозможно обобщить преобразование Леви-Чивита к эквивалентному набору трехмерных переменных (на случай трехмерного пространства).

Изложенные в разд. 3 результаты, полученные ранее автором [8], показывают, что регуляризация дифференциальных уравнений возмущенного плоского движения в задаче двух тел, предложенная Леви-Чивита, может быть с успехом использована для регуляризации дифференциальных уравнений возмущенного пространственного движения в задаче двух тел, если исходные уравнения возмущенного пространственного движения записать в идеальной системе координат, в которой они принимают вид уравнений возмущенного плоского движения, и затем получаемые уравнения дополнить скалярными или кватернионным уравнениями, описывающими изменение ориентации идеальной системы координат в параметрах Эйлера в процессе возмущенного движения тела, а также дополнить уравнениями для кеплеровской энергии и времени, рассматриваемыми в качестве дополнительных переменных (при этом в качестве независимой переменной используется независимая переменная Зундмана).

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

Работа выполнена в рамках темы FFNM-2022-0007.

Список литературы

  1. Челноков Ю.Н. Кватернионные и бикватернионные методы и регулярные модели аналитической механики (обзор) // ПММ. 2023. Т. 87. № 4. С. 519–556.

  2. Levi-Civita T. Sur la regularization du probleme des trois corps // Acta Math. 1920. V. 42. P. 99–144. https://doi.org/10.1007/BF02418577

  3. Штифель Е., Шейфеле Г. Линейная и регулярная небесная механика. M.: Наука, 1975. 303 с.

  4. Aarseth S.J., Zare K.A. Regularization of the three-body problem // Celest. Mech. 1974. V. 10. P. 185–205.

  5. Aarseth S.J. Gravitational N-Body Simulations. Cambridge: Univ. Press, 2003. 408 p.

  6. Hopf H. Uber die Abbildung der dreidimensionalen Sphare auf die Kugelflache // Math. Ann. 1931. V. 104. P. 637–665.

  7. Hurwitz A. Mathematische Werke. Vol. 2. Basel: Birkhauser, 1933.

  8. Челноков Ю.Н. Кватернионная регуляризация в небесной механике и астродинамике и управление траекторным движением. II // Космич. исслед. 2014. Т. 52. № 4. С. 322–336. https://doi.org/10.7868/S0023420614030029

  9. Deprit A. Ideal frames for perturbed keplerian motions // Celest. Mech. 1976. V. 13. № 2. P. 253–263.

  10. Sundman K.F. Memoire sur le probleme des trois crops // Acta Math. 1912. V. 36. P. 105–179.

  11. Velte W. Concerning the regularizing KS-transformation // Celest. Mech. 1978. V. 17. P. 395–403.

  12. Vivarelli M.D. The KS-transformation in hypercomplex form // Celest. Mech. 1983. V. 29. P. 45–50.

  13. Vivarelli M.D. Geometrical and physical outlook on the cross product of two quaternions // Celest. Mech. 1988. V. 41. P. 359–370.

  14. Vivarelli M.D. On the connection among three classical mechanical problems via the hypercomplex KS-transformation // Celest. Mech.&Dyn. Astron. 1991. V. 50. P. 109–124.

  15. Shagov O.B. On two types of equations of motion of an artificial Earth satellite in oscillatory form // Mech. Solids. 1990. № 2. P. 3–8.

  16. Deprit A., Elipe A., Ferrer S. Linearization: Laplace vs. Stiefel // Celest. Mech.&Dyn. Astron., 1994. V. 58. P. 151–201.

  17. Vrbik J. Celestial mechanics via quaternions // Canad. J. Phys. 1994. V. 72. P. 141–146.

  18. Vrbik J. Perturbed Kepler problem in quaternionic form // J. Phys. A: Math.&General., 1995. V. 28. P. 193–198.

  19. Waldvogel J. Quaternions and the perturbed Kepler problem // Celest. Mech.&Dyn. Astr. 2006. V. 95. P. 201–212.

  20. Waldvogel J. Quaternions for regularizing Celestial Mechanics: the right way // Celest. Mech.&Dyn. Astr. 2008. V. 102. № 1. P. 149–162.

  21. Saha P. Interpreting the Kustaanheimo-Stiefel transform in gravitational dynamics // Monthly Notices Roy. Astron. Soc. 2009. V. 400. P. 228–231. https://doi.org/10.1111/j.1365-2966.2009.15437.x.arXiv:0803.4441

  22. Zhao L. Kustaanheimo–Stiefel regularization and the quadrupolar conjugacy // R.&C. Dyn., 2015. V. 20. № 1. P. 19–36. https://doi.org/10.1134/S1560354715010025

  23. Roa J., Urrutxua H., Pelaez J. Stability and chaos in Kustaanheimo-Stiefel space induced by the Hopf fibration // Monthly Notices Roy. Astron. Soc. 2016. V. 459. № 3. P. 2444–2454. https://doi.org/10.1093/mnras/stw780.arXiv:1604.06673

  24. Roa J., Pelaez J. The theory of asynchronous relative motion II: universal and regular solutions // Celest. Mech.&Dyn. Astron. 2017. V. 127. pp. 343–368.

  25. Breiter S., Langner K. Kustaanheimo–Stiefel transformation with an arbitrary defining vector // Celest. Mech.&Dyn. Astron. 2017. V. 128. P. 323–342.

  26. Breiter S., Langner K. The extended Lissajous–Levi-Civita transformation // Celest. Mech.&Dyn. Astron. 2018. V. 130. Art. № 68. https://doi.org/10.1007/s10569-018-9862-4

  27. Breiter S., Langner K. The Lissajous–Kustaanheimo–Stiefel transformation // Celest. Mech.&Dyn. Astron. 2019. V. 131. Art. № 9. https://doi.org/10.1007/s10569-018-9862-4

  28. Ferrer S., Crespo F. Alternative angle-based approach to the KS-Map. An interpretation through symmetry // J. Geom. Mech. 2018. V. 10. № 3. P. 359–372.

  29. Челноков Ю.Н. К регуляризации уравнений пространственной задачи двух тел // Изв. АН СССР. МТТ. 1981. № 6. С. 12–21.

  30. Челноков Ю.Н. О регулярных уравнениях пространственной задачи двух тел // Изв. АН СССР. МТТ. 1984. № 1. С. 151–158.

  31. Челноков Ю.Н. Кватернионные методы в задачах возмущенного центрального движения материальной точки. Ч. 1: Общая теория. Приложения к задаче регуляризации и к задаче о движении ИСЗ. М.: 1985. 36 с. Деп. в ВИНИТИ 13.12.85. № 218628-В.

  32. Челноков Ю.Н. Кватернионные методы в задачах возмущенного центрального движения материальной точки. Ч. 2: Пространственная задача невозмущенного центрального движения. Задача с начальными условиями. М.: 1985. 18 с. Деп. в ВИНИТИ 13.22.85. № 8629-В.

  33. Челноков Ю.Н. Применение кватернионов в теории орбитального движения искусственного спутника. I // Космич. исслед. 1992. Т. 30. Вып. 6. С. 759–770.

  34. Челноков Ю.Н. Применение кватернионов в теории орбитального движения искусственного спутника. II // Космич. исслед. 1993. Т. 31. Вып. 3. С. 3–15.

  35. Челноков Ю.Н. Кватернионная регуляризация и стабилизация возмущенного центрального движения. Ч. 1 // Изв. РАН. МТТ. 1993. № 1. С. 20–30.

  36. Челноков Ю.Н. Кватернионная регуляризация и стабилизация возмущенного центрального движения. Ч. 2 // Изв. РАН. МТТ. 1993. № 2. С. 3–11.

  37. Челноков Ю.Н. Анализ оптимального управления движением точки в гравитационном поле с использованием кватернионов // Изв. РАН. ТиСУ. 2007. № 5. С. 18–44.

  38. Челноков Ю.Н. Кватернионная регуляризация в небесной механике и астродинамике и управление траекторным движением. I // Космич. исслед. 2013. Т. 51. № 5. С. 389–401. https://doi.org/10.7868/S0023420613050026

  39. Челноков Ю.Н. Кватернионная регуляризация в небесной механике и астродинамике и управление траекторным движением. III // Космич. исслед. 2015. Т. 53. № 5. С. 430–446. https://doi.org/10.7868/S0023420615050040

  40. Челноков Ю.Н. Возмущенная пространственная задача двух тел: регулярные кватернионные уравнения относительного движения // ПММ. 2018. Т. 82. № 6. С. 721–733. https://doi.org/10.31857/S003282350002736-9

  41. Челноков Ю.Н. Кватернионные уравнения возмущенного движения искусственного спутника Земли // Космич. исслед. 2019. Т. 57. № 2. С. 117–131. https://doi.org/10.1134/S002342061902002X

  42. Chelnokov Yu.N. Quaternion methods and models of regular celestial mechanics and astrodynamics // Appl. Math.&Mech. 2022. V. 43. № 1. P. 21–80. https://doi.org/10.1007/s10483-021-2797-9

  43. Бордовицына Т.В. Современные численные методы в задачах небесной механики. М.: Наука, 1984. 136 с.

  44. Бордовицына Т.В., Авдюшев В.А. Теория движения искусственных спутников Земли. Аналитические и численные методы. Томск: Изд-во Том. ун-та, 2007. 178 с.

  45. Fukushima T. Efficient orbit integration by linear transformation for Kustaanheimo–Stiefel regularization // Astron. J. 2005. V. 129. № 5. Art. № 2496. https://doi.org/10.1086/429546

  46. Fukushima T. Numerical comparison of two-body regularizations // Astron. J. 2007. V. 133. № 6. Art. № 2815.

  47. Pelaez J., Hedo J.M., Rodriguez P.A. A special perturbation method in orbital dynamics // Celest. Mech.&Dyn. Astron. 2007. V. 97. P. 131–150. https://doi.org/10.1007/s10569-006-9056-3

  48. Bau G., Bombardelli C., Pelaez J., Lorenzini E. Non-singular orbital elements for special perturbations in the two-body problem // Monthly Notices Roy. Astron. Soc. 2015. V. 454. P. 2890–2908.

  49. Amato D., Bombardelli C., Bau G., Morand V., Rosengren A.J. Non-averaged regularized formulations as an alternative to semi-analytical orbit propagation methods // Celest. Mech.&Dyn. Astron. 2019. V. 131. Art. № 21. https://doi.org/10.1007/s10569-019-9897-1

  50. Bau G., Roa J. Uniform formulation for orbit computation: the intermediate elements // Celest. Mech.&Dyn. Astron. 2020. V. 132. Art. № 10. https://doi.org/10.1007/s10569-020-9952-y

  51. Челноков Ю.Н., Логинов М.Ю. Новые кватернионные модели регулярной механики космического полета и их приложения в задачах прогноза движения космических тел и инерциальной навигации в космосе // Сб. матер.: XXVIII С.-Петербургская межд. конф. по интегрированным навигационным системам. С.-Петербург, 2021. С. 292–295.

  52. Челноков Ю.Н., Сапунков Я.Г., Логинов М.Ю., Щекутьев А.Ф. Прогноз и коррекция орбитального движения космического аппарата с использованием регулярных кватернионных уравнений и их решений в переменных Кустаанхеймо–Штифеля и изохронных производных // ПММ. 2023. Т. 87. Вып. 2. С. 124–156.

  53. Челноков Ю.Н. Кватернионная регуляризация уравнений возмущенной пространственной ограниченной задачи трех тел. I // Изв. РАН. МТТ. 2017. № 6. С. 24–54.

  54. Euler L. De motu rectilineo trium corporum se mutuo attrahentium // Nov. Comm. Petrop. 1765. V. 11. P. 144–151.

  55. Levi-Civita T. Traettorie singolari ed urbi nel problema ristretto dei tre corpi // Ann. Mat. Pura Appl. 1904. V. 9. P. 1–32.

  56. Levi-Civita T. Sur la resolution qualitative du probleme restreint des trois corps // Opere Math. 1956. № 2. P. 411–417.

  57. Kustaanheimo P. Spinor regularization of the Kepler motion // Ann. Univ. Turku. 1964. V. 73. P. 3–7. https://doi.org/10.1086/518165

  58. Kustaanheimo P., Stiefel E. Perturbation theory of Kepler motion based on spinor regularization // J. Reine Anqew. Math. 1965. V. 218. P. 204–219.

  59. Hopf H. Uber die Abbildung der dreidimensionalen Sphare auf die Kugelflache // Math. Ann. 1931. V. 104. P. 637–665.

  60. Volk O. Concerning the derivation of the KS-transformation // Celest. Mech. 1973. V. 8. P. 297–305.

  61. Лидов М.Л. Увеличение размерности гамильтоновых систем. KS-преобразование, использование частных интегралов // Космич. исслед. 1982. Т. 20. № 2. С. 163–176.

  62. Лидов М.Л. Метод построения семейств пространственных периодических орбит в задаче Хилла // Космич. исслед. 1982. Т. 20. № 6. С. 787–807.

  63. Лидов М.Л., Ляхова В.А. Семейства пространственных периодических орбит задачи Хилла и их устойчивость // Космич. исслед. 1983. Т. 21. № 1. С. 3–11.

  64. Полещиков С.М. Регуляризация канонических уравнений задачи двух тел с помощью обобщенной KS-матрицы // Космич. исслед. 1999. Т. 37. № 3. С. 322–328.

  65. Stiefel E.L., Waldvogel J. Generalisation de la regularisation de Birkhoff pour le mouvement du mobile dans l’espace a trois dimensions // Comptes Rendus Hebdomadaires des Seances de Lacademie des Sciences. 1965, Paris.

  66. Stiefel E., Rossler M., Waldvogel J., Burdet C.A. Methods of regularization for computing orbits in celestial mechanics // NASA Contractor Rep. NASA CR-769. 1967. P. 88–115.

  67. Birkhoff G.D. The restricted problem of three bodies // Rendiconti del Circolo Matematico di Palermo (1884–1940). 1915. V. 39 (1). P. 265–334.

  68. Waldvogel J. Die Verallgemeinerung der Birkhoff-Regularisierung fur das raumliche Dreikorperproblem // Bull. Astron. Ser. 3. 1967. V. 2. № 2. P. 295–341.

  69. Andoyer H. Cours de Mecanigue Celeste. Paris: Gauthier-Vilars, 1923.

  70. Musen P. Application of Hansen’s theory to the motion of an artificial satellite in the gravitational field of the Earth // J. Geoph. Res. 1959. V. 64. P. 2271–2279.

  71. Брумберг В.А. Аналитические алгоритмы небесной механики. М.: Наука, 1980. 208 с.

Дополнительные материалы отсутствуют.