以下是引用秦汉昌在2009-1-8 23:42:00的发言:
他们的算法不够理想,是指他们的算法不够精确,还是不够简单呢?
算法是简单的,至少那篇文章中给的是粗算.文中也只给了太阳黄经的粗算公式.所以张先生说误差(近二三百年内)在10分钟之内.
十分钟的误差对于日常生活来说无关紧要的,但对于计算天文日历,似乎不够精确了.比如算出来的节气时刻为23:59,若实际精确时刻为0:01,仅仅二分钟误差,日期就不在同一天了.
也许当年的计算机不够普及(文章是94年1月刊出的,那么最晚也是93年写的了),大多数人只好用计算器来辅助计算数据了.
不过,张先生是紫金山天文台的,应当有较高配置的计算机了.(个人电脑93年时586已经流行了。大型计算机在天文台也应当有的,不大可能用手工(计算器)来排历谱的.)
扯远了。回到秦汉昌(这个名字倒有点像笔名)朋友以前提的用黄帝纪日取代儒略日来计算天体物理这个话题来。
我看到的计算太阳月亮位置的公式都与儒略日有关,比如(引用张培瑜先生的计算公式)
Lo=(1006908.04″+129602768.13″T+1.089″T2)/3600
DL1=((6910.057″-17.24″T)×sinM+72.338″×sin(2M))/3600
DL2=-(20.49″+17.2″sinQ)/3600
L=
Lo+
DL1+ DL2
其中:T=(JD-2415020)/36525
M=357.528°+35999.05°E
Q=125.045°-1934.136E
E=(JD-2451545)/36525.
JD——所要求太阳黄经的日期时刻对应的儒略日及小数。
L——太阳黄经
Lo——太阳平黄经
DL1——中心差修正
DL2——光行差、章动修正
如果秦兄直接把黄帝纪日转化为儒略日,那么公式中的系数可以直接套用。(但这样做已经失去了黄帝纪日本身的意义)。公式中2415020,2451545都是比较特殊的日子(儒略日),一个是1899年12月31日12点。一个是2000年1月1日12点。
如果要利用黄帝纪日,那么秦兄要自己推导出一个公式系数,才能避开儒略日这个关键字。
但愿秦兄有耐心去做这件事,等着你的黄帝纪日版的万年历早日推出。