如果你能耐心地从头到尾把这个讨论贴看一遍，你就会知道B4488表示公元前4488年。

而-4487是天文上使用的记年方法（可以说是国际通用的）。

公元1年，对应1年，

公元前1年，对应0年，

公元前2年，对应-1年，

余类推。

公元前4488年，对应-4487年。

 

没有错位的。

暑期考虑再出新版本，增强星历功能。

在牧夫天文论坛转了几圈，没有收获。不过突发奇想：用数值积分方法解决星历问题，应当也是很有趣的方法。

牧夫的朋友，大多不喜欢进行星历计算。

毕竟现成的天文方面的软件很多，且都是图文结合的。很容易查到相关的星历。

所以搞计算的朋友就少了。

 

许兄如果想用数值积分法进行星历的计算，如果是纯力学计算，估计岁差是很难用力学方法计算出来的。

不知许兄用什么方法来解决岁差的问题？

 

还有一个就是实际计算的精度问题。毕竟用现有的星体质量数据，其精度不见得有多高。计算得到的周期估计不会很精确。

（尽管数值积分程序本身的精度很高，但所用的星体的原始数据精度不高啊。）

 

请许兄用牧夫上的程序实算一下，看看你的数值积分程序计算得到的星体位置与你的寿星万年历计算得到的星体位置，二者相差多少？

我们可以用日食食甚时刻来验证。可以查《中国天文年历》关于日食时的一些参数。

 

岁差运动不一定要用数值方法来求解。岁差运动可以用动力学方法来解释，但我们不一定要使用纯数值方法来求解。一般的，岁差运动由观测值直接得到，并适当通过数理方法表达为多项式。

还有，数值积分不需要质量数据，需要的是G*M的值，这个精密数据在DE405星表中有提供。

数值积分方法（如我在牧夫的方法）可以准确反应轨道的变化，至于计算实际星历，确实还有很多问题要解决。

当然，月亮轨道就烦麻许多，二者距离太近，不能把地球看质点。

岁差运动不一定要用数值方法来求解。岁差运动可以用动力学方法来解释，但我们不一定要使用纯数值方法来求解。一般的，岁差运动由观测值直接得到，并适当通过数理方法表达为多项式。

还有，数值积分不需要质量数据，需要的是G*M的值，这个精密数据在DE405星表中有提供。

数值积分方法（如我在牧夫的方法）可以准确反应轨道的变化，至于计算实际星历，确实还有很多问题要解决。

当然，月亮轨道就烦麻许多，二者距离太近，不能把地球看质点。

 

新手机都用半年多了，今天才发现它原来还有GPS功能，挺高兴，在阳台测了一下，我家地标为

经度119.0957° = 119°05'45"，纬度+25.3181° = +25°19'05"，误差30米

 

zb9003发上来的地标是（应用在V4.41版）

经度119°05'，纬度 + 25°19'

 

区政府在我家东面500米左右，靠北约150米。地理坐标每1角秒约30米，如此看来zb9003的数据正好指秀屿区政府，十分准确。区政府办公位置是前几年者确定的，办公楼也刚盖六年，所以我估计zb9003的地标数据库应当是比较可靠的，数据也很新。

 

google卫星地图一直看不清我们这个地方，通过这次实验，终于知道这里的精确地标了。

我所说的质量就是指G*M,我在牧夫里,就是用这种方法计算行星位置的.再推算行星周期.

 

但就算用G*M来代表质量,有效数字并不太多.

 

另外,请问许兄,要不要考虑木星的卫星?土星的卫星?

木卫三、土卫六、木卫四、木卫一都是质量较大的卫星（大于月亮），它们对整个太阳系会不会产生足够大的影响（在几十年甚至几百年内）？

 

我用数值积分程序时，发现把木星及土星的卫星计算进去时，地球的轨道是稍稍有点变化的。

 

如果再考虑柯伊柏带的小天体（也有好几个矮行星，包括冥王星，阋神星（136199 Eris）（质量及体积都可能比冥王星要大），妊神星（136108 Haumea）（质量及体积比冥王星也小不了多少）。

那么，用数值积分来计算行星轨道，在几百年之后及上千年之后，轨道的误差或许也不小的。

 

因为这些海王星之外的天体，目前还没有完全观测清楚，不知道会出现什么样的干扰？