Драконовы узлы Луны
Чтобы построить приближенную математическую модель лунных узлов на основе периодов обращения и начальных положений, нам нужно учесть важный астрономический факт: в идеально изолированной системе (если бы Солнце не влияло на орбиту Луны гравитацией) линия узлов была бы неподвижна в пространстве. Однако из-за солнечного возмущения плоскость орбиты Луны вращается (прецессирует). Линия узлов движется назад (регрессирует) с периодом примерно 
лет (или ~ 798).
лет (или ~ 798).Период прецессии узлов (TΩ) аналитически можно грубо оценить через периоды обращения Луны (Tm) и Солнца (Ts) по формуле Ньютона-Лапласа:
(где i ~ 5.14° — наклон орбиты Луны к эклиптике).
(где i ~ 5.14° — наклон орбиты Луны к эклиптике).
Ниже приведен алгоритм и тригонометрические формулы для расчета долготы Восходящего узла (Ω) и Истинного узла с учетом первой (самой главной) тригонометрической поправки.
Шаг 1. Вводные данные и константы
Определим переменные, исходя из вашего условия:
- Ts — сидерический период обращения Солнца (Земли) ~ 365.256 суток.
- Tm — сидерический период обращения Луны ~ 27.322 суток.
- TΩ — период вращения самих узлов ~ 679$ суток (движутся в обратную сторону).
- t — время (в сутках), прошедшее с опорной (конкретной) даты (t = 0).
Каждому периоду соответствует своя средняя угловая скорость (в градусах в сутки):
- Солнца:
- Луны:
- Узла:
На конкретную дату (t=0) заданы начальные средние долготы: Ls0 (Солнца), Lm0 (Луны) и Ω0 (Восходящего узла).
Шаг 2. Расчет текущих средних положений (Линейная часть)
Для любого момента времени t средние долготы планет и узла рассчитываются линейно:
- Средняя долгота Солнца:
- Средняя долгота Луны:
- Средняя долгота Восходящего узла (Средний узел):
Поскольку узел движется в противоположную сторону (навстречу движению Луны), знак перед скоростью будет минус:
Шаг 3. Тригонометрические формулы поправок (Истинный узел)
Средний узел смещается равномерно. Но гравитация Солнца то ускоряет, то замедляет плоскость лунной орбиты. Чтобы превратить Средний узел в Истинный, нужно рассчитать периодические колебания (эвекцию плоскости орбиты).
Для этого вводятся два фундаментальных тригонометрических аргумента:
- Удлинение Луны (D) — угловое расстояние между Луной и Солнцем:
D = Lm(t) – Ls(t) - Аргумент широты Луны (F) — расстояние Луны от её среднего узла:
F = Lm(t) – Ωmean(t)
Главное тригонометрическое колебание узла зависит от синуса угла между Солнцем и линией узлов. Математически эта поправка выводится через комбинацию аргументов (2D – 2F).
Итоговая приближенная тригонометрическая формула Истинного узла (Ωtrue):
Если подставить значения D и F, то тригонометрический аргумент упрощается до связи Солнца и Узла:
Получаем классическую формулу первого приближения:
Как это работает на практике?
- Когда Солнце находится точно на линии узлов (Ω = Ls), синус равен нулю, и Истинный узел совпадает со Средним. Именно в эти периоды происходят коридоры затмений.
- Когда Солнце находится под углом 45° к узлам, синус равен 1 (или -1), и Истинный узел отклоняется от Среднего на максимальную величину — около 1.483° (или 1°29′).
Эта базовая тригонометрическая формула покрывает более 90% всех отклонений лунного узла и дает точность, достаточную для большинства прикладных задач.
Если вы хотите перевести эти формулы в программный код, укажите, на каком языке программирования его написать, и я подготовлю рабочий скрипт.
