根据2021年天文年历计算非表列的其他时刻的太阳赤经纬

1、利用世界时计算力学时

目的：由于2021年中国天文年历中给出的太阳视赤纬、视赤经是以力学时0h为准

前提知识点：

（1）世界时：我们日常使用的时间系统。以地球对着太阳自转一周作为单位，叫做1昼夜或1日，再细分为小时、分、秒。1秒是1日的86400分之一长。

英国格林尼治天文台所在的时区成为“世界时”。在同一时刻，与北京时间的关系为：

世界时 = 北京时-8h。

（2）力学时：力学时是依据天体动力学理论推算的时间。

（3）力学时与世界时之差表示为：

世界时换算力学时的改正值事先无法获得精确的长期推测值，从IERS网站（http://www.iers.org/）可以获得他们的精确值和长达一年的预推值，本年历（2021）的采用值为。

2、利用插值法（贝塞尔内插公式）计算非表列时刻

算法的详细描述在2021年中国天文年历中，此处只设计所用的部分。

（1）使用至二次差，此时

这包括天文年历各表的大部分，如太阳，月亮每两小时一值的视赤经和视赤纬等。

其中：为表中时刻的数值，t代表所要计算的时刻，为内插引数， , ,  为二次差，利用附表5载的二次差订正

更加详细的计算过程请查看，2021年中国天文年历 page：555-557。

matlab的计算代码后续补上。

参考文献：

[1]2021年中国天文年历

[2]天文普及年历的使用方法_张敏、张元东

[3]用 天 文 测 量 简 历 精 确 计 算 太 阳 位 置 的 方 法