01-07-2023
Уравне́ние Ке́плера описывает движение тела по эллиптической орбите в задаче двух тел и имеет вид:
где — эксцентрическая аномалия, — эксцентриситет орбиты, а — средняя аномалия.
Впервые это уравнение было получено астрономом Иоганном Кеплером в 1619 году. Играет значительную роль в небесной механике.
Содержание |
Уравнение Кеплера в классической форме описывает только движение только по эллиптическим орбитам, то есть при 0 ≤ ε < 1. Движение по гиперболическим орбитам (ε > 1) подчиняется гиперболическому уравнению Кеплера, сходному по форме с классическим. Движение по прямой линии (ε = 1) описывается радиальным уравнением Кеплера. Наконец, для описания движения по параболической орбите (ε = 1) используют уранение Баркера. При ε < 0 орбит не существует.
Рассмотрим движение тела по орбите в поле другого тела. Найдем зависимость положения тела на орбите от времени. Из II закона Кеплера следует, что
Здесь r — расстояние от до тела от гравитирующего центра, υ — истинная аномалия — угол между направлениями на перицентр орбиты и на тело, μ = GM0 — произведение постоянной тяготения на массу гравитирующего тела, a — большая полуось орбиты. Отсюда можно получить зависимость времени движения по орбите от истинной аномалии:
Здесь tp — время прохождение через перицентр.
Дальнейшее решение задачи зависит от типа орбиты, по которой движется тело.
Уравнение эллипса в полярных координатах имеет вид
Тогда уравнение для времени приобретает вид
Для того, чтобы взять интеграл вводят следующую подстановку:
Величина E называется эксцентрической аномалией. Благодаря такой подстановке интеграл легко берется. Получается следующее уравнение:
Величина (μ/a)1/2 является средней скоростью движения тела по орбите. В небесной механике для этой величины используется термин среднее движение. Произведение среднего движения на время называется средней аномалией M. Эта величина характеризует среднее смещение тела по орбите за период.
Таким образом получаем уравнение Кеплера для эллиптического движения:
Уравнение гиперболы в полярных координатах имеет тот же вид, что и уравнение эллипса. Значит, интеграл получается такой же по виду. Однако, использовать эксцентрическую аномалию в данном случае нельзя. Воспользуемся параметрическим представлением гиперболы: , . Тогда уравнение для гиперболы принимает вид
а связь между υ и H
Благодаря такой подстановке интеграл приобретает ту же форму, что и в случае с эллиптической орбитой. После произведения преобразований получаем гиперболическое уравнение Кеплера:
Величина H называется гиперболической эксцентрической аномалией. Поскольку , то последнее уравнение можно преобразовать следующим образом:
Отсюда видно, что E = iH.
Решение уравнения Кеплера в эллиптическом и гиперболическом случаях существует и единственно при любых вещественных M[1]. Для круговой орбиты (ε = 0) уравнение Кеплера принимает тривиальный вид М = E. В общем виде Уравнение Кеплера трансцендентное. Оно не решается в алгебраических функциях. Однако, его решение можно найти различными способами с помощью сходящихся рядов. Общее решение уравнения Кеплера можно записать с помощью рядов Фурье:
где
Этот ряд сходится, когда величина ε не превышает значения предела Лапласа.
Среди численных методов решения уравнения Кеплера часто используются метод неподвижной точки («метод простой итерации») и метод Ньютона[2]. Для эллиптического случая в методе неподвижной точки за начальное значение E0 можно взять M, а последовательные приближения имеют следующий вид[1]:
В гиперболическом случае метод неподвижной точки подобным образом использовать нельзя, однако этот метод даёт возможность вывести для такого случая другую формулу приближений (с гиперболическим арксинусом)[1]:
Уравнение Кеплера.