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


Уравнение времени — разница между средним солнечным временем (ССВ) и истинным солнечным временем (ИСВ), то есть УВ = ССВ − ИСВ. Эта разница в каждый конкретный момент времени одинакова для наблюдателя в любой точке Земли. Уравнение времени можно узнать из специализированных астрономических изданий, астрономических программ или вычислить по формуле, приведённой ниже.

В таких изданиях, как «Астрономический календарь», уравнение времени определяется как разность часовых углов среднего экваториального солнца и истинного солнца, то есть при таком определении УВ = ССВ − ИСВ.

В англоязычных изданиях часто применяется «инвертированное» определение уравнения времени: УВ = ИСВ − ССВ, то есть разница между истинным солнечным временем и средним солнечным временем.

Пояснения к определению

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

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

Неравномерность движения истинного Солнца

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

Влияние эллиптичности орбиты

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

Влияние наклона земной оси

Поскольку плоскость эклиптики наклонена к плоскости небесного экватора, имеет место следующее явление:

  • Солнце вблизи солнцестояний (зимнего и летнего) движется почти параллельно небесному экватору, и скорость его перемещения практически полностью вычитается из суточного движения небесной сферы — результирующая скорость изменения часового угла Солнца минимальна;
  • Солнце вблизи равноденствий (осеннего и весеннего) движется под максимальным углом к небесному экватору, и скорость его перемещения лишь частично вычитается из суточного движения небесной сферы — результирующая скорость изменения часового угла Солнца максимальна.

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

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

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

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

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

Расчёт

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

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

где

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

или

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

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

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

и так далее.

Программа расчета на Ruby для текущей даты

#!/usr/bin/env ruby puts puts " == You're using Ruby ver." + RUBY_VERSION + " ==" =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")

Имя:*
E-Mail:
Комментарий: