Уравнение времени

Можно встретить определение уравнения времени как разницы «местного истинного солнечного времени» и «местного среднего солнечного времени» (в англоязычной литературе — local apparent solar time и local mean solar time). Данное определение формально более точно, но не влияет на результат, так как для любой конкретной точки на Земле эта разница одинакова.

Кроме того, не следует путать ни местное истинное, ни местное среднее солнечное время с официальным местным временем (standard time).

В отличие от звёзд, чьё видимое суточное движение практически равномерно и обусловлено только вращением Земли вокруг своей оси, суточное движение Солнца не равномерно, так как обусловлено и вращением Земли вокруг своей оси, и обращением Земли вокруг Солнца, и наклоном земной оси к плоскости орбиты Земли.

Уравнение времени как сумма составляющих

Вверху слева: влияние эксцентриситета; вверху справа: влияние наклона эклиптики; внизу слева: сумма — уравнение времени; внизу справа: позиция истинного Солнца относительно среднего Солнца.
(Графики приведены в соответствии с «инвертированным» определением уравнения времени, принятым в англоязычной литературе.)

Кривая уравнения времени является суммой двух периодических кривых — с периодами 1 год и 6 месяцев. Практически синусоидальная кривая с годичным периодом обусловлена неравномерным движением Солнца по эклиптике. Эта часть уравнения времени называется уравнением центра или уравнением от эксцентриситета. Синусоида с периодом 6 месяцев представляет разность времён, вызванную наклоном эклиптики к небесному экватору и называется уравнением от наклона эклиптики[1].

Уравнение времени обращается в ноль четыре раза в году — приблизительно 15 апреля, 13 июня, 1 сентября и 25 декабря[3].

Соответственно, в каждое время года существует свой максимум уравнения времени: около 12 февраля — +14,3 мин, 15 мая — −3,8 мин, 27 июля — +6,4 мин и 4 ноября — −16,4 мин. Точные величины уравнения времени даются в астрономических ежегодниках.

Может применяться как дополнительная функция в некоторых моделях часов.

Уравнение можно аппроксимировать отрезком ряда Фурье как сумму двух синусоидальных кривых с периодами, соответственно, один год и шесть месяцев:

${\displaystyle E=7.53\cos(B)+1.5\sin(B)-9.87\sin(2B)}$,

где

${\displaystyle B=360(N-81)/365}$, если углы выражаются в градусах,

или

${\displaystyle B=2\pi (N-81)/365}$, если углы выражаются в радианах.

Там, где ${\displaystyle N}$ — номер дня в году, например:

${\displaystyle N=1}$ на 1 января,

${\displaystyle N=2}$ на 2 января,

и так далее.

#!/usr/bin/env ruby

puts puts "\n\t\== You're using Ruby ver." + RUBY_VERSION + "\t\=="

=begin
Equation of Time calculation
    *** No guarantees are implied. Use at your own risk ***

    Written by E. Sevastyanov, 2017-05-14; clarified by E. Sevastyanov, 2022-01-20

    Based on "Equation of time" WikiPedia article as of 2016-11-28
    (which describes angles in a bewildering mixture of degrees and radians)
    and
    Del Smith, 2016-11-29

    It appears to give a good result, but I make no claims for accuracy.
=end

pi = (Math::PI) # pi
t = Time.now.getutc
yd = t.yday
yy = t.year
delta = yd - 1 # (Текущий день года - 1)
a = t.to_a; 
delta = delta + a[2] / 24.0 + a[1] / 1440.0 + a[0] / 86400.0 # Поправка на дробную часть дня

np = case yy #The number np is the number of days from 1 January to the date of the Earth's perihelion. (http://www.astropixels.com/ephemeris/perap2001.html)
when 2017 ; 3.59583333333333
when 2018 ; 2.23263888888889
when 2019 ; 2.22222222222222
when 2020 ; 4.325
when 2021 ; 1.57708333333333
when 2022 ; 2.28819444444444
when 2023 ; 3.67847222222222
when 2024 ; 2.02708333333333
when 2025 ; 3.56111111111111
when 2026 ; 2.71944444444444
when 2027 ; 2.10625
when 2028 ; 4.51944444444444
when 2029 ; 1.75902777777778
when 2030 ; 2.425 
else; 2
end 

lambda = 23.4372 * pi / 180; # Earth's inclination in radians
omega = 2 * pi / 365.25636 # angular velocity of annual revolution (radians/day)
alpha = omega * ((delta + 10) % 365) # angle in (mean) circular orbit, solar year starts 21. Dec
beta = alpha + 0.03340560188317 * Math.sin(omega * ((delta - np) % 365)) # angle in elliptical orbit, from perigee  (radians)
gamma = (alpha - Math.atan(Math.tan(beta) / Math.cos(lambda))) / pi # angular correction
eot = (43200 * (gamma - gamma.round)) # equation of time in seconds (as The United States Naval Observatory defines it by placing itself at the center of the universe). 
puts "EOT = #{ eot * -1 } секунд" # (as "how many seconds the clock is ahead (+) or behind (−) the apparent sun")

• Аналемма

• Величина колебаний уравнения времени в течение года на портале Гринвичской королевской обсерватории.
• Уравнение времени на сегодняшний день — визуализация.
• Образец построения графика уравнения времени, где прорисованы:

1 — составляющая уравнения времени, определяемая неравномерностью движения Земли по орбите,

2 — составляющая уравнения времени, определяемая наклоном эклиптики к экватору,

3 — уравнение времени.