七政四余排盘软件moira误差分析之力学时误差

力学时的误差:

星历表的计算一般是通过力学时来计算

首先是力学时的转换,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的力学时误差在大部分情况下可以忽略不计

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容