...

среда, 15 января 2014 г.

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

Статьи по небесной механике пользуются на Хабре некоторой популярностью, поэтому я решил рассказать

об одном фундаментальном уравнении движения, а именно, уравнении Кеплера.

Как известно, финитное движение небесных тел в Солнечной системе происходит по эллипсу. Однако, если необходимо

установить, в какой точке небесное тело находится в заданный момент времени,

этой информации недостаточно и надо воспользоваться уравнением Кеплера.

Выведем это уравнение.

Напомню, что эллипс — это сплюснутая на величину


image


окружность. Здесь e — эксцентриситет.

Площадь эллипса равна


image


a — большая полуось.

Согласно второму закону Кеплера, который гласит, что площадь,

заметаемая радиус-вектором небесного тела пропорциональна времени,

можно написать


image


где T — период обращения, t — текущее время,

t0 — момент времени прохождения перигелия (ближайшей к Солнцу точки орбиты).


Найдем теперь зависимость этой площади от положения небесного тела на орбите,

то есть от величины радиус-вектора r и истинной аномалии — угла между перигелием и небесным телом,

если смотреть с Солнца.


Для этого введем дополнительную переменную — эксцентрическую аномалию E.

Направим ось x от Солнца (начала координат), которое находится в одном из фокусов эллипса, в сторону перигелия

по линии апсид (прямой, соединяющей перицентр с апоцентром).

Создадим также вспомогательную окружность с радиусом, равным большой полуоси эллипса.


image


Абсцисса небесного тела в точке P

будет равна


image


ордината будет, соответственно, равна


image


Таким образом, получаем


image


image


Для вычисления площади S рассмотрим две вспомогательные фигуры:

сектор окружности P'-O-Per:


image


Сектор эллипса P-O-Per:


image


Треугольник O-O'-P:


image


Сектор P-O'-Per:


image


Отсюда, прощадь S равна:


image


Сравнивая это выражение с выведенным ранее из второго закона Кеплера,

получаем:


image


Это и есть уравнение Кеплера.


Заметим, что это трансцендентное уравнение относительно E и получить явное решение

в общем случае не удается.


Для вычисления координат небесного тела необходимо найти эксцентрическую аномалию по известному M и e, a затем определить радиус-вектор и угол и, если необходимо x и y по приведенным выше формулам.


Интересным вопросом остается нахождение методов решения уравнения Кеплера.

Этим занимались лучшие умы человечества на протяжении последних четырехсот лет.

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

требует отдельной статьи.


This entry passed through the Full-Text RSS service — if this is your content and you're reading it on someone else's site, please read the FAQ at fivefilters.org/content-only/faq.php#publishers.


Комментариев нет:

Отправить комментарий