Орбітальні елементи Кеплера — пояснення a, e, i, Ω, ω, ν

3D Simulations

🛸 Орбітальна механіка · Астродинаміка

📅 Квіт. 2026 ⏱ ~11 хв читання 🟡 Середній рівень · Останнє оновлення: 3 липня 2026 р.

Орбітальні елементи Кеплера — шість чисел, що описують будь-яку орбіту

Шість чисел — a, e, i, Ω, ω, ν — повністю задають будь-яку кеплерівську орбіту задачі двох тіл, чи то колова траєкторія супутника на висоті 400 км, комета, що пірнає з зовнішньої частини Сонячної системи, чи мільярдний космічний апарат на шляху до Юпітера. Дізнайтеся, що геометрично означає кожен елемент, як їх вимірюють, як перетворити їх на положення та швидкість і як вони повільно змінюються з часом.

Автор: Команда MySimulator · Редакційна перевірка: Редакція MySimulator

Чому саме шість елементів?

Закон всесвітнього тяжіння Ньютона — це звичайне диференціальне рівняння другого порядку у трьох вимірах. Його загальний розв'язок містить 2 × 3 = 6 сталих інтегрування. Ці шість сталих можна вибрати багатьма способами; кеплерівська параметризація обирає їх такими, щоб вони були геометрично змістовними: три описують форму та розмір орбіти, а три описують її орієнтацію в просторі, плюс положення тіла вздовж неї.

Рівносильно, миттєвий стан тіла на орбіті вимагає трьох координат положення (r = [x, y, z]) та трьох компонент швидкості (v = [vx, vy, vz]) у конкретну епоху. Ці шість чисел називають вектором стану. Кеплерівські елементи — це просто інший, більш інтуїтивний базис для того самого шестивимірного простору.

6 Сталих інтегрування для гравітаційного ЗДР задачі двох тіл

3 Елементи, що описують форму та розмір орбіти (a, e, період)

3 Елементи, що орієнтують орбіту в інерціальному просторі (i, Ω, ω)

1 Елемент, що задає положення вздовж орбіти в задану епоху (ν або M чи t = t₀)

Припущення двох тіл: Кеплерівські елементи точні лише в чистій системі двох тіл (одна точкова маса обертається навколо іншої). Реальні орбіти відхиляються через сплюснутість Землі, гравітацію третього тіла, атмосферний опір та тиск випромінювання. Тоді елементи стають оскулюючою орбітою — найкраще наближеним кеплерівським еліпсом у кожну мить, що повільно змінюється з часом.

Шість кеплерівських елементів

Стандартні оскулюючі орбітальні елементи означено нижче. Опорна площина та напрямок змінюються залежно від застосування: для супутників Землі опорною площиною є екваторіальна площина, а опорним напрямком — точка весняного рівнодення (γ); для геліоцентричних орбіт обидві відносять до системи відліку екліптики/J2000.

a

Велика піввісь

Половина довжини найдовшого діаметра еліпса. Безпосередньо пов'язана з орбітальним періодом T та енергією E через T² ∝ a³ (третій закон Кеплера) та E = −GM/(2a).

Діапазон: 0 < a < ∞ (км або а.о.)

e

Ексцентриситет

Форма конічного перерізу: e = 0 коло → 0 < e < 1 еліпс → e = 1 парабола (траєкторія втечі) → e > 1 гіпербола. Відношення c/a, де c — відстань до фокуса.

Діапазон: 0 ≤ e < 1 (замкнені орбіти)

i

Нахил

Кут між орбітальною площиною та екваторіальною (чи екліптичною) площиною. i = 0 — пряма екваторіальна; i = 90° — полярна; i > 90° — ретроградна (обертання протилежне до обертання Землі).

Діапазон: 0° ≤ i ≤ 180°

Ω

RAAN — пряме сходження висхідного вузла

Кут від точки весняного рівнодення γ до висхідного вузла (де орбіта перетинає екваторіальну площину, прямуючи на північ), виміряний в екваторіальній площині. Визначає поворот орбітальної площини навколо полярної осі.

Діапазон: 0° ≤ Ω < 360°

ω

Аргумент перицентру

Кут від висхідного вузла до перицентру (точки найближчого підходу), виміряний в орбітальній площині. Орієнтує еліпс усередині орбітальної площини. Для екваторіальних орбіт (i = 0), Ω + ω → довгота перицентру ˜ω.

Діапазон: 0° ≤ ω < 360°

ν

Істинна аномалія

Кут від перицентру до поточного положення супутника, виміряний у фокусі. Це єдиний залежний від часу елемент. У перицентрі ν = 0; в апоцентрі ν = 180°. Пов'язаний з часом через рівняння Кеплера.

Діапазон: 0° ≤ ν < 360°

Відстані перицентру та апоцентру

Щойно a та e відомі, найближча та найдальша точки випливають негайно:

rₚ = a (1 − e) ← перицентр (перигей / перигелій) rₛ = a (1 + e) ← апоцентр (апогей / афелій) період T = 2π √(a³ / μ) ← третій закон Кеплера (μ = GM)

Угоди про назви: Перицентр та апоцентр отримують різні префікси залежно від центрального тіла — перигей/апогей (Земля), перигелій/афелій (Сонце), перийовій/апойовій (Юпітер), перицентр/апоцентр (загальний випадок).

Рівняння Кеплера та поширення орбіти

Щоб знайти, де тіло перебуває в конкретний момент часу t, ми використовуємо ланцюжок з трьох аномалій:

Середня аномалія M: зростає лінійно з часом. M = n(t − t₀), де n = 2π/T — середній рух, а t₀ — епоха проходження перицентру.
Ексцентрична аномалія E: геометричний кут, означений на допоміжному колі, описаному навколо еліпса.
Істинна аномалія ν: фактичний кут у фокусі від перицентру до тіла.

Розв'язування рівняння Кеплера

Зв'язок між середньою та ексцентричною аномаліями — це трансцендентне рівняння Кеплера:

M = E − e sin(E) ← рівняння Кеплера

Його неможливо обернути аналітично. Стандартний числовий розв'язок використовує ітерацію Ньютона–Рафсона, починаючи з E₀ = M:

Eₙ₊₁ = Eₙ − (Eₙ − e sin Eₙ − M) / (1 − e cos Eₙ)

Збіжність зазвичай досягається за 5–10 ітерацій для e < 0.9, а для орбіт з високим ексцентриситетом (комети, e → 1) — методом Лаґерра.

Від ексцентричної до істинної аномалії

tan(ν/2) = √((1+e)/(1−e)) ⋅ tan(E/2) або рівносильно: cos ν = (cos E − e) / (1 − e cos E) sin ν = √(1−e²) sin E / (1 − e cos E)

Орбітальний радіус у будь-якій точці випливає з рівняння конічного перерізу (vis-viva в кутовій формі):

r = a (1 − e²) / (1 + e cos ν) ← рівняння орбіти (фокус у початку координат)

Рівняння vis-viva дає швидкість на будь-якому радіусі: v² = μ(2/r − 1/a), незалежно від ексцентриситету, щойно відомі a та r.

Перетворення на положення та швидкість (вектор стану)

Повне перетворення орбітальних елементів на вектор стану ‘r, v’ в інерціальній системі відліку — це триетапний процес:

Розв'яжіть рівняння Кеплера для ексцентричної аномалії E, потім знайдіть істинну аномалію ν.
Знайдіть положення та швидкість у перифокальній системі (початок у фокусі, вісь x — до перицентру, вісь z — нормаль до орбітальної площини):

r⃗ₘₑₙᴿ = [r cos ν, r sin ν, 0] ← r = a(1−e²)/(1+e cos ν) v⃗ₘₑₙᴿ = (μ/h) [−sin ν, e + cos ν, 0] ← h = √(μ a (1−e²))

Поверніть до інерціальної системи (IJK) за допомогою трьох послідовних поворотів:

R = R₃(−Ω) ⋅ R₁(−i) ⋅ R₃(−ω) де R₃(θ) повертає навколо z, R₁(θ) — навколо x r⃗ = R ⋅ r⃗ₘₑₙᴿ , v⃗ = R ⋅ v⃗ₘₑₙᴿ

Обернена задача — перетворення вектора стану на орбітальні елементи — вимагає обчислення моменту імпульсу h = r × v, вектора ексцентриситету e = v × h/μ − r̂ і подальшого виведення кожного елемента геометрично. Це процедура, яку виконують мережі стеження за супутниками за радарними вимірюваннями.

TLE — двохрядковий набір елементів: Найпоширеніший формат для розповсюдження орбітальних даних супутників. TLE кодує видозмінений набір середніх елементів (за моделлю поширення SGP4/SDP4, яка враховує J2 та опір) як два рядки по 69 символів. Понад 27 000 об'єктів відстежують у такий спосіб у Мережі космічного спостереження США.

Приклади реальних орбіт

МКС a = 6 731 км · e ≈ 0.0001 · i = 51.6° · T ≈ 92 хв

GPS a = 26 560 км · e ≈ 0 · i = 55° · T = 12 год (півзоряний)

GEO a = 42 164 км · e ≈ 0 · i ≈ 0° · T = рівно 24 год

Місяць a = 384 400 км · e = 0.055 · i = 5.1° · T = 27.3 доби

Особливі типи орбіт

Сонячно-синхронна орбіта (SSO): Нахил ~98° робить керовану J2 прецесію вузлів рівно одним обертом на рік, тож площина орбіти завжди обернена під тим самим кутом до Сонця — ідеально для дистанційного зондування.
Орбіта «Молнія»: e ≈ 0.74, i ≈ 63.4°, T = 12 год. Нахил 63.4° «заморожує» апсидальну прецесію (ω ̇ = 0), тож апогей залишається над північною півкулею ~8 годин за прохід.
Орбіта «Тундра»: як «Молнія», але T = 24 год, що дає одне вікно покриття високих широт на добу.
Заморожена орбіта: комбінація e та ω, підібрана так, що збурення J2 та J3 взаємно знищуються, зберігаючи висоту перицентру сталою. Mars Reconnaissance Orbiter використовує заморожену орбіту.

Комета Галлея — екстремальний еліпс

Комета Галлея (1P/Halley) має a = 17.9 а.о. та e = 0.967, що дає перигелій 0.59 а.о. (усередині орбіти Венери) та афелій 35.1 а.о. (за Нептуном). Її період ~75 років, проте її швидкість у перицентрі ~54 км/с набагато вища за орбітальну швидкість Землі ~30 км/с — усе тому, що vis-viva вимагає ½mv² + Eₘₙₙ = стала.

Збурення та J2

Реальні орбіти відхиляються від чистого кеплерівського руху через кілька ефектів. Найважливіший для супутників Землі — J2, член сплюснутості в геопотенціалі Землі (екваторіальний радіус Землі перевищує її полярний радіус на ~21 км).

Вікові ефекти J2

Ω̇ = − (3/2) n J₂ (R⊕/p)² cos i ← регресія вузлів ω̇ = (3/4) n J₂ (R⊕/p)² (5 cos²i − 1) ← апсидальна прецесія де p = a(1−e²), J₂ = 1.08263×10⁻³, R⊕ = 6371 км

Прирівнявши ω̇ = 0, отримуємо критичний нахил 5 cos²i − 1 = 0, який дає i = 63.4° або 116.6°. Саме тому радянські супутники «Молнія» усі використовують 63.4°.

Інші збурення

Атмосферний опір: Вікове зменшення a та e, що зрештою призводить до входу в атмосферу. МКС втрачає ~2 км висоти на місяць без піднімання орбіти.
Гравітація третього тіла: Місяць і Сонце спричиняють припливи, що збурюють високоорбітальні траєкторії (GEO, HEO) на масштабах років.
Тиск сонячного випромінювання (SRP): Плоска панель площею 1 м² на 1 а.о. зазнає ~4.6×10−6 Н/м². Критично для великих космічних апаратів із сонячним вітрилом.
Загальнорелятивістська прецесія: Перигелій Меркурія зміщується на 43″/століття понад ньютонівські передбачення — одне з перших підтверджень ЗТВ.

Оскулюючі елементи: За наявності збурень ми означуємо орбітальні елементи в кожну мить як елементи найкраще наближеного кеплерівського еліпса (дотичного до фактичної траєкторії за положенням і швидкістю). Ці «оскулюючі» елементи повільно змінюються з часом і є величинами, що табулюються в ефемеридах.

Спробуйте симуляції: Орбітальна механіка дозволяє візуалізувати прецесію та опір у реальному часі; Відхилення астероїда показує, як мала Δv змінює всі шість елементів одночасно.