所以我觉得，如何仅仅把十体（从太阳起到海王星加上月亮，正好十体）的G*M代入计算，估计最终的误差不一定小。可能精确度要小于DE405/DE406的星历。

恭喜许兄在地球上找到了自己的位置。

改天可以用导弹轰炸你了。

 

我有杭州的1：1万的地形图。地图上每1mm代表10米。已经相当精确了。

所以这个地理信息是不能外露的（属于国家机密啊。），万一不友好的国家发导弹过来，那么就完了。

 

有些地方，若我们公司去搞测量，那么可以拿到1：500的地形图。太恐怖，图上1mm代表50厘米。马路上一个窨井盖都可以从图上找到的。

刚才看了一下DE405中质点受力的表达式，很复杂，与我使用的质点受力表达式只有一项相同，即牛顿万有引力。DE405计算质点引力包含了广义相对论的改正，用级数展开式计算。还计算了3粒较大的小行星的作用力，此外，对与地球、月亮、火星，还考虑了其它小行星的作用力。DE405表达式复杂的产要原因是考虑了相对论。而且要求计算引力时要迭代计算完成，当然，我写的建个积分器计算加速度也是迭代计算完成的，思路与DE405的相同。

可见，计算不精确的主要原因有两个：没有考虑小行星的作用力（牛顿引力），其二没有考虑相对论。

看来，要想计算出与DE405相同的数据，还是用它们的积分器为好。

 

DE星历表我看到一句话

However,use of Newtonian accelerations for these terms is sufficiently accurate.

也许是说，用牛顿力学可以算得很准确，只是还不够精确。

我想，误差的主要原因来自于3个小行星Ceres,Pallas,and Vesta。有空的时候我认直计算一下。

 

经他细分析，主要误差也不是小行星，而是初始数据与模型不匹配。

支持，非常了不起

你过虑了。

现行历法已经考虑这个情况了。

所以规定冬至日一定安排在子月（建寅为正的历法，子月就是十一月）。

 

不存在冬至日漂移现象的。

至于雨水，历法上没有强制规定它是哪个月。但由于冬至日固定在子月，雨水日也不会漂移太远的。

 

所以现行农历是完全可以正常运行下去的。

 

我已经推算过，1000万年内，历法都可以不用修改。

5000万年内，历法稍作修改就可以继续使用。

 

我是按月亮每年远离地球3厘米，再考虑地球公转周期变短，自转周期变长的现行天文数据进行推算的。

当然，几万年几百万年后，天文上的观测数据或许会改变的。但到了那个时候，自然会有高人来修改历法。

也不用我们来担心历法能否使用的问题了。

就这水平还历改--冬至是原点--怎么漂移

综合考察农历，已经没有太多的实用价值。

建议，以现行格历为基础，以年为单位，标注24节气日期时刻、月朔月望日期时刻、年首日干支及其他天文事件就可以了。

无须对农历的月序及闰月进行编排。以最靠近立春的朔日为春节（或者由国家正式将春节定为立春日）。逐渐淡化农历的日序。至于民间历的编制不予引导和干涉。

真正对农事有影响习惯的，是节气。至于观潮，完全可以从朔日自行推算。

你的想法很好。

同样想法的人也很多，请看“历改委西电”的一些帖子，在“历改信息”专题里可以查到相关的建议。

 

最后，我建议，有志于改历的朋友，请号召一下，去找管理历法的部门。

目前官方管理历法的机构是“中国科学院紫金山天文台”。

若不去找官方的机构商谈改历之事，只是在民间讨论改历之事，估计现行历法不大可能更改的。

 

另外，国家已经把一些重要的农历日子作为“法定节假日”了。

比如春节，中秋节，比如端午节，是法定放假的日子。

若使用公历（格历）的月日序，那么这三个农历节日就没直观上的文化内涵了。

 

所以，想改历的话，还得找国务院去协商。请他们把放假的日子改一下。

 

偏题了。

今后若有讨论改历的话题，请另开贴吧。

 

这个帖子主要讨论许剑伟先生的寿星万年的程序及相关问题。