第一章 关于寿星万年历.doc

第一章 关于寿星万年历 一、《寿星万年历》概述 欢迎使用寿星万年历。寿星万年历是一款采用现代天文算法制作的农历历算程序，可以很方便地进行公历、农历、回历三历之间的转换。它提供公元-4712年到公元9999年的日期查询功能，其中-721年到1960年的农历数据已经与陈垣的《二十史朔闰表》、张培瑜的《三千五百年历日天象》、方诗铭的《中国史历日和中西历日对照表》核对；它还含有从公元前2000多年以前到今的基本年号。万年历中含有二千多个国内城市的经纬度，并且用户可根据自已的需要扩展经纬度数据。此外，本万年历还提供了历谱数据导出功能。 本软件是一款真正精准的年代跨度大的日历工具，可作为一般的实用日历工具，也可供历史学家放心地使用(最好读完此文，以免误解)；对历算工作者、天文爱好者均有较大的算法参考价值。 在天文学方面，《寿星万年历》以精确著称。过去几百年中，寿星万年历的误差是非常小的，节气时刻计算及日月合朔时刻的平均误差小于1秒。当使用力学时作为时间基础时，与现今最精密星历表比较，太阳坐标的最大可能误差为0.2角秒，月亮坐标的最大可能误差为3角秒，平均误差为最大可能误差的1/6，这使得精确的日月食计算成为可能。寿星万年历是目前世界上唯一一款能够精确计算日月食过程的万年历。软件中还可以查询日出日没、月出月没时间等，与国内外著名的天文软件或天文年历比较，结果完全一致。此外，万年历中还即时显示了“月亮被照亮部分比例”、太阳和月亮的地心黄道坐标、地心赤道坐标、站心地平坐标等。其中日月的地心黄道、赤道采用高精度算法；“月亮被照亮部分比例”的计算采用低精度的算法，但可以满足一般需求；地平坐标的计算已适当考虑了大气折射，由于大气的气压及气温变化具有许多不确定的因素，所以只考虑平均情况下的大气折射修正，在本软件中，地平坐标被描述为方位角与高度角，它是以观测点为中心的坐标，所以已经考虑了周日视差修正。屏幕显示的方位角从正南向西测量。方位角与高度角已转换到是站心地平坐标，并且做了视差修正，同时在地平真纬度大于0时进行大气折射修正。 [天文算法]：寿星万年历在进行月历推算时采用了先进的现代天文算法，它是基于2002年巴黎天文台西尔特处Jean CHAPRONT 和 Gerard FRANCOU发表的ELP/MPP02月球运动理论，在几百年范围内，该理论与当今世界上公认的最精密的DE405/406系列星历表仅相差3毫角秒，在几千年范围相差不超过4角秒。由于农历的计算不需要这么高的精度，我们做了一些截断处理及实质性的简化处理，但仍可以满足高精度历算的要求。程序中采用了1988年法国巴黎天文台Pierre Bretagnon与Gerard Francou发表的VSOP87行星运动理论计算太阳坐标。由于该理论与DE405/406星历表存在少量差异，我们对VSOP87进行少量修正，提高了计算精度。节气及合朔时刻问题在天文学层面看来，就是已知天体的坐标反求时间的问题。在VSOP87或ELP中，天体的坐标是时间的函数z = f(t)，所谓的求节气时刻或月相时刻就是已知z求t，显然这是在求解一个关于t的方程。伟大的英国天文学家物理学家牛顿给出了一种非常有效的迭代算法：牛顿求根法。用这种方法，求t所花费的时间仅是求f(t)花费时间的1.2——1.3倍，并且迭代引起的额外时间误差只零点零几秒。牛顿求根法在一般的《数值方法》书籍中均有介绍。 [日食月食]：万年历中提供的图表，可以精确的确定日月食情况。计算结果与《2008年中国天文年历》比对，日食误差在1秒左右，月食误差在3秒以内，但应注意到，《2008年中国天文年历》提供的月食精度只有0.1分，即+-3秒。事实上，受到地球大气的影响，月食的精度无法计算得很高，所以计算时把地球看成了圆球，而不是椭球。《寿星万年历》日月食预报及计算的本质就时把日月坐标以图表的形式展示出来，所以寿星万年历没有完整的计算日月食（日食中心线、地方食表、日食南北界及月球本影半影的计算除外），用户通过寿星万年历得到的日月食精确数据，仍是用户自已动手的结果。因为用户判断食甚比较困难一些，所以程序中给出了食甚、初亏等关键时间的计算。通过日食图表与用户互动观察，有利于学习一些比较专业的科普知识，体验一下DIY的乐趣。当然，如果用户觉得很有必要加上完整的日月食计算（如各种日食界线及相关图表），我可以考虑改进。另外：如果您利用寿星万年历得到的日食计算结果与媒体上公布的或《中国天文年历》上的数据相差2秒钟以上，多半是地理坐标数据的采用值不同或是媒体的错误造成的，请认真核对。比如《2008年中国天文年历》第526页，给出了2007年3月19日阿勒泰市(新疆)的日食计算及结果，其实这里有个低级错误，日期应改为2008年8月1日。 [日月出没]：受到各地的地形地貌、大气状态