Уравнение Кеплера¶

Уравнение Кеплера — это трансцендентное уравнение, связывающее положение небесного тела на эллиптической орбите со временем его движения. Уравнение является фундаментальным в небесной механике и позволяет определить истинную аномалию (угол между радиус-вектором тела и направлением на перицентр) в любой момент времени, зная среднюю аномалию и эксцентриситет орбиты. Впервые выведено Иоганном Кеплером в 1609 году в работе «Новая астрономия» при описании движения Марса.

¶История открытия

Иоганн Кеплер, анализируя многолетние наблюдения Тихо Браге за движением Марса, сформулировал три закона движения планет. Второй закон Кеплера (закон площадей) гласит, что радиус-вектор планеты заметает равные площади за равные промежутки времени. Однако для практического вычисления положения планеты на орбите требовалось связать время движения с угловыми параметрами. Кеплер вывел уравнение, которое в современной записи имеет вид:

\[ M = E - e \sin E \]

где \(M\) — средняя аномалия, \(E\) — эксцентрическая аномалия, \(e\) — эксцентриситет орбиты. Кеплер решил это уравнение численно для конкретных значений, используя метод последовательных приближений. В 1650-х годах уравнение было формализовано и введено в математический обиход.

¶Математическая формулировка

¶Основные параметры орбиты

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

• Средняя аномалия \(M\) — угол, который прошло бы тело, двигаясь равномерно по окружности с периодом, равным орбитальному. Выражается как \(M = n(t - t_0)\), где \(n\) — среднее движение (2π/T), \(t\) — текущее время, \(t_0\) — время прохождения перицентра.
• Эксцентрическая аномалия \(E\) — вспомогательный угол, определяемый через геометрическое построение: из центра эллипса проводится окружность радиусом, равным большой полуоси, и проекция точки на орбите на эту окружность даёт угол \(E\).
• Истинная аномалия \(\nu\) — угол между радиус-вектором тела и направлением на перицентр.

¶Уравнение Кеплера

Уравнение Кеплера для эллиптической орбиты имеет вид:

\[ M = E - e \sin E \]

Для параболических орбит (\(e = 1\)) используется уравнение Баркера, а для гиперболических (\(e > 1\)) — гиперболическое уравнение Кеплера:

\[ M = e \sinh H - H \]

где \(H\) — гиперболическая эксцентрическая аномалия.

¶Связь с истинной аномалией

После нахождения \(E\) истинная аномалия \(\nu\) вычисляется по формулам:

\[ \tan\frac{\nu}{2} = \sqrt{\frac{1+e}{1-e}} \tan\frac{E}{2} \]

или через тригонометрические функции:

\[ \nu = 2 \arctan\left( \sqrt{\frac{1+e}{1-e}} \tan\frac{E}{2} \right) \]

¶Решение уравнения Кеплера

Уравнение Кеплера является трансцендентным и не имеет аналитического решения в элементарных функциях. Для его решения применяются численные и приближённые методы.

¶Метод последовательных приближений (итераций)

Наиболее распространённый метод — итерация по формуле:

\[ E_{n+1} = M + e \sin E_n \]

Начальное приближение обычно берётся \(E_0 = M\) (для малых эксцентриситетов) или \(E_0 = M + e \sin M\) (более точное). Метод сходится для всех \(e < 1\), причём скорость сходимости линейна.

¶Метод Ньютона

Для ускорения сходимости используется метод Ньютона (касательных):

\[ E_{n+1} = E_n - \frac{E_n - e \sin E_n - M}{1 - e \cos E_n} \]

Этот метод имеет квадратичную сходимость и требует 3–5 итераций для достижения точности \(10^{-12}\) при типичных эксцентриситетах.

¶Разложение в ряд

Для малых эксцентриситетов (\(e \ll 1\)) уравнение можно решить разложением в ряд по степеням \(e\):

\[ E = M + e \sin M + \frac{e^2}{2} \sin 2M + \frac{e^3}{8} (3 \sin 3M - \sin M) + \dots \]

Этот ряд сходится для \(e < 0,6627\) (предел Лапласа).

¶Специальные функции

В XIX веке Фридрих Бессель предложил решение через функции Бесселя:

\[ E = M + \sum_{k=1}^{\infty} \frac{2}{k} J_k(ke) \sin(kM) \]

где \(J_k\) — функция Бесселя первого рода. Этот ряд сходится для всех \(e < 1\).

¶Применение в астрономии и космонавтике

¶Расчёт эфемерид

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

¶Орбитальные манёвры

В космонавтике уравнение применяется для расчёта времени перелёта между двумя точками орбиты (задача Ламберта), а также для планирования манёвров перехода между орбитами (например, гомановский переход).

¶Двойные звёзды

Для визуально-двойных звёзд уравнение Кеплера используется для определения орбитальных параметров по наблюдениям взаимного положения компонентов.

¶Обобщения и модификации

¶Релятивистское уравнение Кеплера

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

\[ M = E - e \sin E + \frac{3GM}{c^2 a(1-e^2)} (E - e \sin E) \]

где \(G\) — гравитационная постоянная, \(c\) — скорость света, \(a\) — большая полуось.

¶Уравнение для возмущённого движения

При наличии возмущений (гравитация других тел, негравитационные силы) уравнение Кеплера становится приближённым, и для точных расчётов используются методы численного интегрирования (например, метод Рунге — Кутты).

¶Интересные факты

• Кеплер первоначально вывел уравнение в геометрической форме, без использования тригонометрических функций. Современная запись появилась позже, благодаря работам Эйлера и Лагранжа.
• Уравнение Кеплера — один из первых примеров трансцендентного уравнения, решённого численно. Кеплер сам выполнил 70 итераций для одного из значений.
• Для круговых орбит (\(e = 0\)) уравнение вырождается в тривиальное \(M = E\).
• Предел сходимости ряда Лапласа (\(e < 0,6627\)) связан с радиусом сходимости разложения и был строго доказан только в XX веке.
• В 2020-х годах уравнение Кеплера продолжает использоваться в программном обеспечении для навигации космических аппаратов, включая миссии к Луне и Марсу.

¶Критика и ограничения

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

¶Источники

1. Кеплер И. Новая астрономия (Astronomia Nova). — 1609.
2. Брумберг В. А. Релятивистская небесная механика. — М.: Наука, 1972.
3. Мюррей К., Дермотт С. Динамика Солнечной системы. — М.: Физматлит, 2010.
4. Danby J. M. A. Fundamentals of Celestial Mechanics. — Willmann-Bell, 1988.
5. Vallado D. A. Fundamentals of Astrodynamics and Applications. — Microcosm Press, 2007.