Уравнение времени
Благодаря нецентральности земной орбиты и отклонению оси по отношению к орбитальной плоскости, видимое движение Солнца меняется в течение года. Это приводит к отклонениям во временных значениях, полученных солнечными часами (истинное солнечное время) по сравнению с более точными часами, например, теми, которые выдают время, определённое видимым вращением неподвижных звёзд (звёздное время). Данная разность известна как уравнение времени.
SiderealTime выдаёт значения звёзного времени в смешанных величинах, используя MixedUnit.
![Click for copyable input](assets.ru/the-equation-of-time/In_47.png)
SiderealTime[]
![](assets.ru/the-equation-of-time/O_43.png)
![Click for copyable input](assets.ru/the-equation-of-time/In_48.png)
SiderealTime[];
QuantityUnit[%]
![](assets.ru/the-equation-of-time/O_44.png)
Звёздное время по сути является углом и, таким образом, может быть переведено в единицы измерения углов.
![Click for copyable input](assets.ru/the-equation-of-time/In_49.png)
UnitConvert[SiderealTime[], "AngularDegrees"]
![](assets.ru/the-equation-of-time/O_45.png)
Уравнение времени определено движением Солнца, заданным SunPosition. Рассчитайте недельные значения в полдень по горному времени.
![Click for copyable input](assets.ru/the-equation-of-time/In_50.png)
dates = DateRange[
DateObject[{2015, 1, 1}, TimeObject[{12, 0}, TimeZone -> 0],
TimeZone -> 0.`],
DateObject[{2015, 12, 31}, TimeObject[{12, 0}, TimeZone -> 0],
TimeZone -> 0.`], 7];
![Click for copyable input](assets.ru/the-equation-of-time/In_51.png)
sunpos = SunPosition[GeoPosition[{0, 0}], dates,
CelestialSystem -> "Equatorial"]
![](assets.ru/the-equation-of-time/O_46.png)
Используйте SiderealTime для расчёта прямого восхождения по меридиану Гринвича для этих же дат.
![Click for copyable input](assets.ru/the-equation-of-time/In_52.png)
stime = SiderealTime[GeoPosition[{0, 0}], dates]
![](assets.ru/the-equation-of-time/O_47.png)
И, наконец, найдите соответствующие значения уравнения времени как разницу между прямым восхождением Солнца и звёздным временем в полдень в течение года. Выразите результат в минутах.
![Click for copyable input](assets.ru/the-equation-of-time/In_53.png)
equationoftime =
UnitConvert[(TimeSeriesMap[First, sunpos] - stime),
"MinutesOfRightAscension"]
![](assets.ru/the-equation-of-time/O_48.png)
Абсолютные значения уравнения времени могут быть больше 15 минут.
![Click for copyable input](assets.ru/the-equation-of-time/In_54.png)
MinMax[equationoftime]
![](assets.ru/the-equation-of-time/O_49.png)
![](assets.ru/the-equation-of-time/O_50.png)