力学时的误差:
星历表的计算一般是通过力学时来计算
首先是力学时的转换,moira基于swisseph来计算,swisseph对力学时有一套拟合公式,随着时间的增大力学时误差逐渐变大
力学时(TT)与协调时(UTC)的差值,由观测获取,其值可以到www.iers.org查询
路径为IERS /Data / Products / Tools Earth orientation data
Long term earth orientation data->EOP C01 IAU2000 (1900-now)
国际原子时(TAI)协调世界时(UTC)与世界时(UT1)
可以查到UT1-TAI的值
2005.30 -32.5872399
2005.35 -32.6041043
2005.40 -32.6149983
2023.25 -37.0239469
2023.30 -37.0237672
2023.35 -37.0225477
2023.40 -37.0160727
= -37 s
通过iers查询UT1-UTC的值为在2005年的差值为0.6秒左右,在2023年为0.001秒以下,忽略不计,如果该差值大于0.9秒,授时机构会增加一秒闰秒来减小误差。
因此可以近似认为UTC=UT1,上述的UT1-TAI值可以近似为UTC-TAI
UTC时间时我们通过互联网/电视/手机获取到的时间,可以认为就是看到的0时区的时钟时间
再根据TT = TAI + 32.184 seconds,其差值为固定的32.184秒
以2023.4年的力学时为例:UTC-TAI = -37.016,TAI = UTC + 37.016
TT = TAI + 32.184 = 37.016 + 32.184 = UTC + 69.20秒,这个69.20秒一般称为
可以得到结论:2005年清明前后差值delta T(ΔT)为64.788
2023年清明前后差值delta T(ΔT)为69.204
然后来看moira计算得到的delta T(ΔT):
2005年清明 delta T(ΔT)为65.225, 误差为0.5秒左右
2023年清明 delta T(ΔT)为86.593, 误差为17秒左右
然后来看未修正前的二十四节气
moira计算2005年清明交节时间为2005年4月5日0时34分15秒
2023年清明交节时间为2023年4月5日9时12分45秒
同时在互联网上查询得到清明节交节的信息:
2005年0时34分17秒交节,误差2秒
中国天文年历显示清明2023年9时13分交节,误差15秒
于是修改moira源码,校正delta T的误差得到新计算的清明交节时间:
2005年清明交节时间为2005年4月5日0时34分15秒,误差不变
2023年清明交节时间为2023年4月5日9时13分2秒,误差减小15秒
结论:由于swisseph的delta T拟合公式的精度与实际上的delta差了10秒左右,导致计算出来的数据可能会有十几秒的误差,但也在可接受范围内,因此moira的力学时误差在大部分情况下可以忽略不计