太阳位置 (VSOP87)：VSOP87 定义了太阳系行星运动的解析解-matlab开发

VSOP87 的主要版本将行星的椭圆变量定义为时间的函数。 VSOP 数据提供以下椭圆变量： A 天文单位中的半长轴（AU，其中 1AU = 149597870.7km） L 平均经度k = e cos π 其中 e 是偏心率，π 是近日点的经度h = e sinπ q = sin i/2 cos Ω 其中 i 是倾角，Ω 是升交点的经度p = sin i/2 sin Ω 计算基于 J2000 纪元 2000-01-01 12:00:00 之后的儒略世纪时间。 T = (JDE - 2451545)/365250