Чи можна передбачити рух планети за допомогою математики?

Ми можемо подивитися на планету сьогодні. Але як заздалегідь вказати, де вона буде через місяць, коли того неба ще ніхто не бачив? Самого малюнка орбіти замало: він показує шлях, проте не каже, коли планета дістанеться певної його ділянки.

Математика дозволяє передбачити положення планети, якщо задано модель руху та початковий стан. Надійність такого прогнозу перевіряють новими спостереженнями. Початковий стан включає положення і швидкість у визначений момент; швидкість має і величину, і напрям. JPL: як визначають орбіти й ефемериди.

До шляху треба додати годинник і причину руху

Якщо намалювати замкнену криву й поставити на ній крапку, лишиться багато можливостей. Тіло може рухатися в один чи інший бік, швидше чи повільніше. Навіть та сама початкова точка ще не визначає майбутнього. Треба знати, як тіло рухається зараз і за яким правилом цей рух змінюється.

У фізичній моделі таким правилом є закон руху під дією сил. Планета продовжує рух за інерцією, а тяжіння зорі змінює її швидкість. Для кругової орбіти тяжіння весь час повертає напрям швидкості до нового дотичного напрямку. За відповідної початкової швидкості відстань до зорі залишається сталою. Це спеціальний випадок, а не властивість будь-якого запуску тіла. OpenStax: механіка кругової орбіти.

Отже, причинний ланцюжок має вигляд: сили та початковий стан → зміна руху → положення в наступні моменти. У загальній задачі комп’ютер обчислює ці зміни послідовно. Для деяких спрощених випадків є готовий математичний розв’язок, який одразу пов’язує час із положенням. Його й використаємо, щоб навчитися читати прогноз.

Побудуймо одну уявну орбіту з явними умовами

Уявімо планету навколо зорі, значно масивнішої за неї. Нехтуємо рухом самої зорі та впливом інших тіл. Орбіта — коло радіусом 1 au, а один оберт триває 360 діб. Це спеціально підібрані навчальні числа для уявної системи, не рік Землі й не прогноз для Сонячної системи. Одна астрономічна одиниця, au, — приблизно 150 мільйонів кілометрів. JPL: визначення au.

Радіус, період і тяжіння зорі мають бути сумісними з круговим рухом: для заданої реальної зорі радіус і період не можна незалежно обрати будь-якими. Тут ми задаємо можливу ідеалізовану систему, щоб відділити читання графіка від складнішого розрахунку сил.

Виберімо нерухомі осі в площині орбіти. Зоря буде початком координат; додатний напрям $x$ — праворуч, $y$ — вгору. Дивимося перпендикулярно до цієї площини. У момент $t=0$ планета праворуч від зорі, а рух спрямований проти годинникової стрілки. Період — час повного оберту й повернення в те саме положення з тим самим напрямком руху.

За рівномірного кругового руху за чверть періоду планета проходить чверть кола. Через 90 діб вона буде вгорі, через 180 — ліворуч, через 270 — внизу, через 360 — знову праворуч. Ми отримали конкретні положення для конкретних моментів, а не лише форму шляху.

Графік координати — інший рисунок, ніж орбіта

На першому рисунку обидві осі описують простір: де планета праворуч чи ліворуч, вище чи нижче зорі. На другому горизонтальна вісь — час, а вертикальна — лише координата $x$ . Точка B позначає той самий момент на двох рисунках, але осі мають різний зміст.

Одна орбіта, два різні рисунки

Навчальна модель: R = 1 au, T = 360 діб. Це не орбіта Землі.

Вигляд на площину орбіти

Координата x із часом

Зоря — початок; осі нерухомі. Дивимося перпендикулярно до площини. Рух A → B → C → D → E.

Ті самі моменти в обох рисунках — обчислення моделі
Точка	t (доби)	Кут	x (au)	y (au)
A	0	0°	+1	0
B	90	90°	0	+1
C	180	180°	−1	0
D	270	270°	0	−1
E	360	360°	+1	0

x = 0 не означає, що планета опинилася біля зорі. У B та D відрізняється y. Крапки й зорю збільшено; точки не є спостереженнями.

На графіку $x$ спадає від +1 au до −1 au, а тоді зростає назад. Це не означає, що планета рухається вперед і назад по прямій. Так змінюється лише її проєкція на вісь $x$ , поки саме тіло обходить зорю по колу. Відстань до зорі весь час становить 1 au; від’ємною буває координата, а не відстань.

Особливо показові B та D. В обох випадках $x=0$ , проте в B планета вгорі, а в D — внизу. Однієї координати недостатньо, щоб відновити повне положення. Треба також знати $y$ або кут і радіус. Крім того, $x$ змінюється нерівномірно, хоча швидкість уздовж кола стала: біля правого й лівого країв рух спрямований переважно вздовж $y$ .

Відгалуження · З формуламиНеобов’язково · можна пропустити

Як час перетворюється на дві координати

Позначимо радіус $R$ , період $T$ , а кут від додатного напряму $x$ — $\theta$ . Для нашого початку відліку та напряму руху:

\theta(t)=360^\circ\frac{t}{T}.

Час $t$ і період $T$ мають бути в однакових одиницях. Їхнє відношення — частка оберту. Косинус і синус переводять кут у проєкції радіуса на дві осі:

\begin{aligned} x(t)&=R\cos\theta(t),\\ y(t)&=R\sin\theta(t). \end{aligned}

Це параметризація: спосіб задавати координати через спільний параметр, тут час. Косинусоїду на графіку отримано з цієї геометрії, а не підібрано на око до набору крапок. Сам запис синуса й косинуса ще не пояснює, яка сила підтримує орбіту.

Наприклад, через 45 діб:

\begin{aligned} \frac{t}{T}&=\frac{45}{360}=\frac18,\\ \theta&=45^\circ,\\ x&\approx0{,}707\,\mathrm{au},\\ y&\approx0{,}707\,\mathrm{au}. \end{aligned}

Планета перебуває праворуч угорі, а не на відстані 0,707 au від зорі: поєднання обох координат дає радіус 1 au. При $t=0$ формули повертають $(R,0)$ ; при $t=T$ — те саме положення після одного оберту. Якщо використовувати програмні функції sin/cos, їхній кут зазвичай задають у радіанах: замість $360^\circ$ у першій формулі беруть $2\pi$ .

Намальований прогноз ще треба зіставити з небом

У нашій схемі всі крапки обчислені за одним правилом. Їхня відповідність лінії перевіряє узгодженість рисунка й обчислень. Вона ще не показує, що якась реальна планета рухається саме так.

Для перевірки спершу фіксують модель і її параметри, потім порівнюють прогноз із новими вимірами. Якщо прогнозують положення на небі із Землі, треба обчислити відповідний напрям від спостерігача. Його не можна просто прирівняти до координати в системі із зорею в центрі. Для точного порівняння важливий і час, за який світло доходить до спостерігача. JPL Horizons: координати та таблиці для спостерігача.

Розбіжність оцінюють разом із невизначеністю вимірів. Вона може вказувати на неточно визначений початковий стан, неврахований вплив чи проблему спостереження. Нові дані допомагають уточнювати розв’язок; одного збігу недостатньо, щоб оголосити модель безпомилковою для всіх майбутніх моментів.

Відгалуження · ГлибшеНеобов’язково · можна пропустити

Чим підгонка відрізняється від перевірки

Уявімо, що період і початкове положення ми добирали так, щоб крива проходила близько до вже відомих вимірів. Це підгонка параметрів. Вона потрібна, адже початковий стан реального тіла невідомий із безмежною точністю. Проте ці самі виміри вже вплинули на вибір кривої; збіг із ними не є незалежною перевіркою.

Корисно запитати, чи узгодиться прогноз зі спостереженням, яке не використовували під час підгонки. І чи лишається модель придатною в інший час? Хороше відтворення минулої ділянки не гарантує вдалого продовження далеко за неї.

JPL окремо публікує наближені формули положень планет із підігнаними орбітальними елементами та швидкостями їхньої зміни. Для них задані часові інтервали й оцінки похибок; переносити їх за межі цих інтервалів не дозволяє сама форма формул. JPL: наближені положення планет.

Що додають справжні ефемериди

Реальні планети не обходять Сонце рівномірно по однакових колах. У двотільному наближенні орбіти описують еліпсами; ближче до Сонця планета рухається швидше. NASA: закони Кеплера. Взаємне тяжіння тіл додатково змінює рух, тому він не повторюється абсолютно точно. JPL: межі незмінної орбіти.

Ефемериди задають обчислені положення та інші характеристики руху для визначених моментів. Наприклад, JPL створила DE440 і DE441, чисельно розв’язуючи рівняння руху та узгоджуючи результат із наземними й космічними спостереженнями. Це датований приклад справжнього розрахунку, а не назва нашої кругової моделі. Park та інші, 2021: оригінальна праця про DE440/DE441.

Для запиту до системи на кшталт Horizons треба вказати тіло, момент, центр або спостерігача, систему координат та одиниці. «Де буде Марс?» без цих умов — неповне питання. Ефемерида є результатом моделі, обмеженої даними й точністю обчислень; вона не є фотографією майбутнього.

Повернімося до навчальної схеми: через 450 діб координата $x$ знову дорівнює нулю. Чи повернулася планета в початкову точку? Ні. 450 діб — один повний оберт і ще чверть. Планета в B: $x=0$ , $y=+1$ au. До A вона поверталася через 360 діб. Тут потрібні і час, і зміст координат, а не лише висота крапки на графіку.

Отже, передбачати рух можна: фізичні закони й початковий стан дають майбутні положення, а нові спостереження перевіряють точність. Графік робить цей зв’язок видимим лише тоді, коли ми розуміємо його осі та межі моделі. Погляньте ще раз на склад Сонячної системи: кожне додаткове тіло в цій картині може стати частиною точнішого розрахунку.