《宋代定朔计算精度分析_滕艳辉.pdf》由会员分享，可在线阅读，更多相关《宋代定朔计算精度分析_滕艳辉.pdf（9页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、2014年7月咸阳师范学院学报Jul.2014第29卷 第4期Journal of Xianyang Normal UniversityVol.29 No.4科技史与科技传播研究宋代定朔计算精度分析滕艳辉(咸阳师范学院 数学与信息科学学院，陕西 咸阳 712000)摘要：定朔是中国古代天文学中一个十分重要的概念。梳理了宋代历法中定朔推步的过程，严格按照历法术文将 纪元历 等11部历法的定朔算法编写成计算机程序，并计算出各部历法行用期间的所有定朔时刻，与定朔的真实值相比较，得出各部历法推算定朔的客观精度。使用古今两种观点对历法的精度进行分析，得到在古人眼中各部历法定朔的计算精度。利用周期函数法分
2、析定朔误差的来源，得到各历法的经朔误差。关键词：定朔；精度分析；纪元历；统天历中图分类号：N09；P1-092文献标识码：A文章编号：1672-2914（2014）04-0059-09AnAccuracyAnalysis of the Real New MoonAlgorithm in Song DynastyTENG Yan-hui(School of Mathematics and Information Science,Xianyang Normal University,Xianyang712000,Shaanxi,China)Abstract:The real new moon is
3、 an extremely important concept in Chinese ancient astronomy.Inthis paper,the calculation process of real new moon in calendars of the Song Dynasty was introduced atfirst.Then the real new moon algorithms in 11 calendars including Jiyuan Calendar were compiled intocomputer programs strictly accordin
4、g to their prediction texts.The paper calculated all the real newmoon moments by using the methods of each calendar.The objective accuracy of these calendars is esti-mated by comparing these calculations with true values,and the calculating accuracy is obtained fromthe perspectives of the ancients t
5、hrough the use of both ancient and modern viewpoints.Finally,the er-ror of mean new noon of each calendar can be calculated by analyzing the sources of the real new moonerror with a method of periodic function.Key words:the real new moon;accuracy analysis;Jiyuan Calendar;Tongtian Calendar收稿日期：2014-0
6、3-21基金项目：国家自然科学科研基金项目（11073016）；咸阳师范学院科研基金项目（12XSYK009）。作者简介：滕艳辉（1983-），男，内蒙古赤峰市人，咸阳师范学院数学与信息科学学院讲师，研究方向为中国数理天文学。定朔是古人对平朔进行日月不均匀性修正，得到的真朔。定朔的使用提高了历谱编排和日月五星计算的精度。1950年代以来，一些研究者在定朔算法及相关的研究方面作出了重要的工作。文献1-5对定朔算法进行了详细的解读；文献6-7对定朔算法的精度进行了比较；滕艳辉等8-10着重研究了宋代历法中的定朔计算方法，并对一些历法的精度进行了分析。本文将在这些研究基础上对宋代的定朔算法和计算精度
7、进行全面总结和更深入的探索。本文完全遵循现存宋代历法定朔推步的全过程，借助计算机，计算出每一部历法所有月份的定朔时刻，并与其真实的定朔值比较，研究历法中定朔算法的计算精度，最后结合推算误差，对影响各历法推算定朔的因素进行探讨。1宋代历法定朔算法综述两宋期间，共行用了18部官方历法，其中有14部完整地保存在 宋史 中。这些历法普遍使用定朔算法。历法推算定朔的基本方法是：先根据基本天文常数等数据计算出要求得年月的经朔时刻，再根据此时太阳、月亮的相关数据求得日月不均匀性的修正数值，最后在经朔上加入太阳和月亮的不均匀性的改正则可求得定朔时刻。其算法的基本模型是8：Trm=Tpm+(s+m)/A（日）（
8、1）式中Trm表示定朔时刻，Tpm表示经朔时刻，s和m则分别表示日、月改正数。具体每一部历法，在求经朔和日月改正数的方法上会略有不同。崇天历 使用传统的“积年法”求经朔，使用表格法求日月改正数。北宋的 纪元历乾元历 和南宋的 会元历 统元历 淳熙历 乾道历 成天历 开禧历 共8部现存历法的定朔算法与 崇天历 的相同。明天历 和 观天历 在经朔计算上与 崇天历 同，而日月改正数则使用“相减相乘”法直接用公式计算。南宋 统天历 求日月改正数方法与 崇天历 相同，而求经朔上则使用其首创的“岁实消长法”。它不使用传统的上元积年，其近地点与冬至点、合朔点也不重合。宋初的 仪天历和 应天历 由于包含了伊斯
9、兰历法的内容，其定朔计算方法与中国传统的算法有些区别，而对其定朔算法，令人满意地详细解读仍未出现。由于使用上元积年，利用古代历法进行定朔推算，计算量是很大的。但定朔算法只涉及了加减乘除运算，其模型十分简洁，现代电子计算机却很胜任这样的工作。我们将各历法的定朔算法模型编写成计算机程序，就可以方便地计算出每一个定朔时刻。为了提高结果的计算精度，程序化的过程严格按照历法中给出的算法进行。这样，使用程序推算就相当于在复原古代历法的当时推算情况。本文选择了Visual Basic语言作为实现宋代历法定朔算法的语言，编写应用程序。用此程序可以方便地求出任意月份的定朔时刻。92 定朔计算的精度分析目前的程序
10、，已经能够对除宋初的 应天历 仪天历 二历和 明天历 外的其余11部现存宋代历法进行定朔的推算。现代天文学利用迭代法能对定朔时刻进行精确计算，其计算值可以作为定朔的真实值9。使用 三千五百年历日天象 中“合朔满月表”所给出的定朔计算结果作为定朔的真实值11。为了方便比较，将真实值的结果转化成以小时为单位，并将历法计算值转化成依北京时间所表示的时刻。设时差为，根据公式（1）得Trmbj=Trm+=Tpm+(s+m)/A+。（2）经计算，对于 崇天历 和 乾元历=27.6 min（注）；对于 统天历 开禧历 和 成天历 三历=-0.8 min；对于 纪元历 等其它5历=22.6 min（见文献8）
11、。利用公式(2)，就可以计算出时差修正后的定朔值。2.1 定朔推算的整体水平推算值与真实值的差距（误差）可以用来衡量历法精度的高低。使用误差的标准差、平均值和变化趋势进行描述和评价。历法值偏离真实值的平均大小用误差的平均值来反映；偏离的范围和程度用误差绝对值的标准差描述；给出一条拟合直线，其斜率可以反映推算精度随时间变化的趋势和程度8。在图1至11中，列出了各历定朔计算误差的拟合直线。这里必须指出，在文献8-9中，作者已经分析了宋代 纪元历 等历法的定朔计算精度。但是由于作者没有将历法中的某些术文和天文常数进行校正，直接依据 宋史 律历志 所记载的原文进行推算，因此得到的历法定朔计算值是错误的
12、。本文中，我们详细校正了这些历法中每一段术文和每个常数的正确性，并严格检验程序的正确性，确保历法计算结果的真实可信。在文献10中，作者改正了历法常数和程序，得到的 开禧历 的定朔计算精度是正确的，本文在讨论 开禧历 精度时，将直接引用此结果。11部历法总共2 770个定朔的推算中，最大误差为2.498 9 h；若按“四刻以下为近”12作为划分依据，则误差在1 h以下的可认为满足要求，占95.60%；而误差大于1 h的有122次，占4.40%。定朔计算较为精确的，按“二刻以下为亲”判断12，误差小于0.5 h的有1 934次，所占比例为69.82%。在所有推算的误差中，平均为0.064 4 h，
13、约合4 min；若取其绝对值均值为0.391 5 h，约为23.5 min；标准差是0.488 7。各部历法的详细统计数据如表1所示。把这11部历法分为行用时间较长的（一般超过10 a）和行用时间较短的两组，并具体分析每一部历法的定朔计算精度。乾元历 定朔计算误差绝对值均值为0.476 6 h，合28.6 min，其余指标都在所有11部历法平均水平之下。最大误差为2.498 9 h，这是11部历法唯一一 应天历 和 仪天历 的定朔包含伊斯兰天文学的内容，与中国传统历法定朔计算方法不同，我们拟用专门的文字进行讨论，本文暂不进行计算；由于使用 明天历 的天文常数计算得到其月亮近地点位置与其余宋代历
14、法所计算的位置相差近30度，我们怀疑 明天历 的常数或是术文有误，还需要进一步斟酌，故暂时还未编写 明天历 定朔计算程序。古代历法计算出的定朔值是以当时帝都或观测台所在的位置为标准的，即采用当地的地方时计时，而北京时间是120E的地方时。各造历地点的地理经度分别是：阳城为113.1E，俊仪为114.35E，临安为120.2E。60咸阳师范学院学报第29卷0601201802403210-1-2误差/h=0.0007x+0.176,R2=0.0069（a）982985988991994997100010.50-0.5平均值绝对值均值误差/h（b）图1乾元历 定朔误差逐年变化图公元纪年定朔推算的次
15、数表1纪元历 等11部历法的推朔精度次大于2 h的误差。乾元历 误差大于1 h的有21次，几乎都大于0，说明历法值晚于真实值，即“推算后天”。乾元历 推算误差情况如图1（a）所示。取每一年计算的平均值并分类统计，结果在图1（b），其余历法的描述与 乾元历 类似。崇天历 定朔误差平均值接近于0，而绝对值均值略大于所有历法平均值，为0.399 7 h，合24 min。历法在40年间最大误差只有1.4 h，并且精度每年降低仅为0.13 s，说明单从定朔上看，崇天历 还是很优秀的。其定朔误差图如图2所示。历法乾元历崇天历观天历纪元历统元历乾道历淳熙历会元历统天历开禧历成天历总计定朔个数23559312
16、430939611117399111545742770误差最值/h2.498 91.403 7-1.658 71.132 81.167 01.247 11.542 41.526 71.310 11.102 01.034 22.498 9标准差0.556 20.499 30.464 60.464 40.459 70.440 30.426 40.431 10.455 00.409 90.417 30.488 7平均值/h0.256 10.010 9-0.134 5-0.141 2-0.000 60.275 40.299 00.558 40.082 40.141 90.188 10.064 4绝对值
17、均值/h0.476 60.399 70.374 50.386 50.375 70.418 60.405 30.611 40.368 00.337 90.374 40.391 51 h个数21346147712153101122拟合线斜率0.000 70.00000 3-0.001 3-0.000 30.0000 70.002 40.0000 7-0.003 6-0.002 3-0.000 1-0.005 70100200300400500600210-1-2误差/h=0.000003x+0.176,R2=0.000001（a）102410321040104810561064107410.50-
18、0.5平均值绝对值均值误差/h（b）图2崇天历 定朔推算误差图公元纪年定朔推算的次数平均值绝对值均值第4期滕艳辉：宋代定朔计算精度分析61纪元历 的各推算指标接近11部历法的平均值。其行用前期的推算好于后期，纪元历 行用前期（11051127年），有11次误差大于1 h，后期（11331135 年）有 3 次。其定朔推算误差绝对值均值为0.385 5 h，约23 min。纪元历 定朔精度如图3所示。统元历 的定朔推算中，误差绝对值在1 h以上的仅有7次，是所有历法中最少的；有279次误差在0.5 h以下，超过70%。其误差平均值仅为-0.000 6 h，合2 s多；取绝对值均值是0.376 7
19、 h，约为22.5 min。统元历 推朔精度明显高于北宋时期的任何一部历法。在 统天历 重测天文常数前，推算定朔方面，统元历 是最优的一部历法。其定朔推算精度如图4。淳熙历 的推算误差中，有12次在1 h以上；116次在0.5 h以下。淳熙历 为负值的误差明显少于为正值的，尤其是在0.5 h以上的那些，历法推算定朔后天。其误差最大值为1.542 4 h，这是比较大的。淳熙历 的推朔精度倒退到了北宋初期的水平。其误差分布如图5所示。开禧历 各项精度指标都在平均水平之上。历法行用了44 a的545次定朔计算中，只有10次推算误差大于1 h；而在0.5 h以下的精密推算超过七成。0501001502
20、00250300210-1-2误差/h=0.0003x-0.095,R2=0.0032（a）11061111111611211126113110.50-0.5-1误差/h（b）图3纪元历 定朔推算误差图公元纪年定朔推算的次数平均值绝对值均值050100150200250300350400210-1-2误差/h=0.00007x-0.0148,R2=0.0003（a）1136 11401144 1148 1152 11561160 1164 116810.50-0.5误差/h（b）图4统元历 定朔推算误差图公元纪年定朔推算的次数平均值绝对值均值0306090120150180210-1误差/h=
21、0.00007x-0.2933,R2=0.00006（a）1177117911811183118511871189119110.50-0.5误差/h（b）图5淳熙历 定朔推算误差图公元纪年定朔推算的次数平均值绝对值均值62咸阳师范学院学报第29卷在推算定朔方面，开禧历 无疑是宋代最优秀的历法。其定朔计算精度如图6。2.2观天历 等5部历法定朔推算精度观天历 共有124次定朔计算，误差大于1 h的有6次，小于0.5 h的有95次，其定朔推算精度处于中等水平，略好于 崇天历 和 纪元历，与 统元历 的相当。然而其最大误差却达到了-1.658 7 h，仅次于 乾元历 的最大值。由此可知，推算定朔平均
22、精度方面，北宋历法要低于南宋历法。因这些历法行用时间较短，仅给出每次定朔误差的描述图，不再求每一年的平均值。图7描述了 观天历 的定朔误差推算精度。乾道历 仅有的111次定朔推算中，误差在1 h以上的就有7次，而在0.5 h以下的却仅有76次。乾道历 推朔精度降低到了 崇天历 时的水平，比其前面的 统元历 要低很多。会元历 的精度更低，其误差绝对值平均达到了0.611 4 h，是所有历法中最低的。会元历 共99个推算误差中，只有41次小于 0.5 h，却有15次大于1 h。两历的定朔推算精度如图8和图9所示。统天历 精确测量了天文常数，也提高了其定朔计算的精度。统天历 误差绝对值均值为0.36
23、8 0 h，只有3次大于1 h的误差。其后面的 成天历，各项指标基本与 统天历 相当。授时历 以前，在推算定朔的精度上 统天历 之后的宋代历法达到了顶点。图10和11分别描述了 统天历 成天历 2历的定朔误差。所有的11部历法中，误差平均值最小值为 统元历 的0.000 6 h，接近0；崇天历 淳熙历 统天历3历也在0.1 h以下；最大平均值是 会元历，超过了30 min。而取绝对值的误差平均值 开禧历 最小，仍是 会元历 最大。误差标准差最小的是 开禧历，而除 乾元历 外均小于0.5。而若考查历法的080160240320400480560210-1-2误差/h=0.0001x+0.174,
24、R2=0.002（a）120812151222122912361243125010.50-0.5误差/h（b）图6开禧历 定朔推算误差图公元纪年定朔推算的次数平均值绝对值均值12141618110112110-1-2误差/h图7观天历 定朔推算误差图定朔推算的次数=-0.0013u+0.0536,R2=0.011163146617691106210-1误差/h图8乾道历 定朔推算误差图定朔推算的次数=-0.0024x+0.1393,R2=0.03161163146617691210-1误差/h定朔推算的次数=-0.0036x+0.7409,R2=0.0591图9会元历 定朔推算误差图第4期滕艳
25、辉：宋代定朔计算精度分析63平均值绝对值均值最大值精度发展趋势，可以考察其误差拟合直线的斜率，斜率为正值表明其历法推算值会随时间推移早于真实值，斜率为负，则表明会晚于真实值。斜率绝对值以成天历 会元历 乾道历 3历较大，而以 崇天历的最小。计算得到，成天历 的精度变化为每年降低4.104 min，而 崇天历 则为每年0.129 6 s。各部历法精度指标比较如图12所示。3古人眼中的历法定朔推算精度制造历法所需的基本天文常数是在观象台（天文台）得到的，而历法是在帝都使用的，如果两地的位置距离较远，造成两者经度或纬度上的差异过大，那么同样的推算结果，依现代和古人两种视角去评价，得到的结论就可能不同
26、。我们先前得到精度可能不是古人认为的这样9。南宋的 统元历 等4历行用地帝都临安与观测点俊仪上相差了近6的经度，折合24 min。这样，使用式（2）计算修正后的定朔时，如果取为北京时间与造历观测点地方时的时间差，计算结果就是依现代人的眼光来衡量的客观精度；如果取为北京时间与行用帝都地方时的时间差，则计算结果就是依当时人视角的历法推算精度9。纪元历统天历 开禧历 成天历 4历，帝都与观测点一致，两种情形下的值相同。对于 乾元历 崇天历 观天历 3历取值为22.6 min；其余南宋历法取为-0.8 min9。古人眼中的推算误差就是在已经得到的误差上加上或是减去观测地与帝都的时差后的所得。例如，统元
27、历 等南宋前期历法，古人认为的误差是在已得误差上减去俊仪临安两地时差0.40 h。由于 乾元历 崇天历 观天历 3历观测点与帝都经度差仅为1.2，古今两种方式考查精度差别不显著。纪元历 精度是0.375 8 h，达到北宋历法的最高水平。但是若以帝都临安为依据计算，它在南宋行用时，误差竟达到了0.732 3 h，约合44 min，定朔计算严重后天。37次定朔计算中，就有10次误差超过1 h，这在当时人们眼中是很低的精度。当对比表1和表2中有关数据时，发现虽然观测纪元历 行用前期，观测点俊仪就位于其都汴梁，而行用后期，其都城在监安，仍要分别讨论其精度。北宋帝都汴梁地理经度为114.35E，南宋帝都
28、临安为120.2E。210-1误差/h=0.0023x+0.2102,R2=0.02611163146617691106图10 统天历 定朔推算误差图3210-1乾元崇天观天纪元统元乾道淳熙会元统天开禧成天图12宋代11部历法的定朔推算精度比较图图13古人眼中的定朔推算精度比较图定朔推算的次数111213141516171210-1误差/h定朔推算的次数=-0.0036x+0.7409,R2=0.0591图11成天历 定朔推算误差图误差/h平均值绝对值均值最大值3210-1乾元崇天观天纪元前期纪元后期统元乾道淳熙会元统天开禧成天误差/h64咸阳师范学院学报第29卷历法乾元历1崇天历观天历纪元历
29、/前2纪元历/后统元历乾道历淳熙历会元历统天历开禧历成天历总计定朔个数2355931242723739611117399111545741 768平均值0.176 1-0.069 1-0.214 5-0.1277-0.634 2-0.393 9-0.117 9-0.081 00.165 10.082 40.141 90.188 10.064 4绝对平均值0.457 30.407 70.392 00.375 80.732 30.500 60.362 80.353 00.373 90.368 00.337 90.374 40.391 5最大值2.418 91.346 31.738 71.129 6
30、1.526 11.423 01.231 11.162 41.133 31.310 11.102 01.034 22.498 91h个数1935911103624231011170.5h个数1444099019012111751327080407511 043表2纪元历 等11部历法的推朔精度附注1：为了便于比较，将 乾元历 等3部历法的数据自表1中摘录至此。附注2：纪元历 行用前期选择1106-1127年的数据，依行用地在开封推算；其行用后期选择1133-1135年数据，依行用地在临安推算。点与帝都不一致的历法，其误差差距很大，但它们也表现出了规律性。表1中误差较小的 统元历 和 纪元历 二历
31、在表2中的精度却很差。相反，表1中精度最低的 会元历 在表2中却要好很多。这种差异和比较，我们能够对 统天历 以前的南宋历法改革及当时天文学家的认识进行如下推测。南宋初年，由于继续使用的 纪元历 是北宋历法，所依据观测数据是俊仪的，不加修正，直接在临安使用产生较大误差是显然的，这样就不得不改历。统元历 制造者陈得一使用的观测数据仍是崇宁年间的，尽管他提高了历法的基本天文常数的精度，使用了比 纪元历 更优的回归年和近点月等常数13，使得 统元历 定朔推算的实际精度好于 纪元历。然而他可能并没有认识到地点的东西差距与天象发生时间的关系，因此在实际使用时，统元历 并不理想，严重后天，因为 统元历 精
32、确推算出俊仪的天象，而临安的实际天象在时间上已经超前了24 min。统元历 在当时人们眼中的精度很差，甚至还不如 乾元历。而刘孝荣等在改历的过程中，就意识到了东西两地存在时间上的差异，并且应该知道时间上西方的要晚于东方的。于是，在 乾道历 淳熙历 会元历 3历的编制中，他们对部分天文常数做了调整，向前推移了天象发生时间，更接近于实际天象在临安发生的时间。而这样做使得历法推算肯定会后于在俊仪的真实天象。表1中，这3历误差的平均值均小于0就反映了这一点。由于地球自转，经度相异的两地，地方时不同，经度每差15时间大约相差1 h。这样，经度不同的两地某种天象就会有不同的发生时刻（地方时）。中国古代对地
33、理经纬度概念的认识要比对经度的认识早一些。因为冬至时正午地面影长在不同纬度的差异较为明显，早在唐代 大衍历 时就有关于不同纬度正午太阳高度及日长差异的算法。而在古代中原地区由经度差仅导致几十分钟的地方时差异，到元代才出现对经度概念的认识。随元太祖西征时，耶律楚材用金大明历 推算某个时刻在撒马尔罕应该有月食出现，但却要推迟发生。他意识到这种天象的时间差是由东西距离差所造成的，于是提出“里差”的概念，来解决时差问题，并进献了 庚午元历14。而证明在此之前中国人认识到经度或是时差的证据还没有出现。从前面的分析，我们可以猜想，耶律楚材以前的南宋历法家可能就已经清楚了东西时间差异的问题。不管人们当时是否
34、知道该现象的原因，但他们已经能够解决历法上东西时间差异的问题。在北宋第4期滕艳辉：宋代定朔计算精度分析65此处 开禧历 的经朔误差与文献9得到的有微小的差异，这是由于本文的程序是在文献9基础上又略作了优化，但这种优化是十分微小的。开封很优秀的 纪元历，在临安使用时就会严重后天，根据 纪元历 而制造的 统元历 仍是后天，并且更为严重。然而，只需向前推移推算结果，就能提高在临安的推算精度，也只能向前移，而不能向后推。以上分析，我们甚至可以做出结论，在南宋 乾道历时，中国天文学家已经能够知道位于西方的地点，要比位于东方地点的时间更晚，这是人类认识时间的一大进步。然而，这些也仅是我们的推测，还需要找到
35、史料上证据来证明。4 定朔计算的误差分析历法精度可以使我们了解当时历法的优劣及推算天象的准确程度，而我们现在还需要关注的是为什么会有这样的精度，哪些算法或参数是影响精度的主要因素。由于定朔是在经朔的基础上加入日月改正数，那么定朔的误差就可以分为经朔推算的误差和日月改正数的误差。文献9以 开禧历 为例，讨论了影响其定朔精度的各种因素，并给出一种周期函数法，可以求得其经朔误差对定朔的影响。这种周期函数法仅适用于行用时间超过9年或112个朔望月的历法，并且在历法行用初期所得的“经朔误差”更接近于实际情况。这样，能够使用周期函数法求经朔的历法只有 崇天历 等6部历法。对于这6部历法，我们先求得各部历法
36、行用初年的9年（112个朔望月）的定朔误差的平均值，然后再此基础上减去期间56个朔望月误差的累积值，得到各部历法的经朔误差。例如，开禧历 的行用初年起9年的定朔误差平均值为8.080 9 min，其朔望月常数为29.530 591 716 d，可知其朔望月常数误差是0.398 2 s，那么 开禧历 的经朔推算误差以约4 s每年的速度增大。这样，除去其朔望月误差的累积值，其行用初年经朔误差的平均值应该是7.709 2 min。实际上，这与误差拟合直线的截距很接近。纪元历 等6部历法的经朔误差列于表3中。历法乾元历崇天历纪元历统元历淳熙历开禧历朔望月历法值/d29.530 592 0029.530
37、 594 9029.530 589 8529.530 591 6329.530 595 7429.530 591 72朔望月真实值/d29.530 587 0229.530 586 7029.530 586 8829.530 586 9529.530 587 0429.530 587 11常数误差累积值/h0.034 1270.011 0170.003 9870.006 2910.011 6940.006 194定朔9年均值/h0.253 968-0.002 786-0.109 838-0.020 6970.321 0650.134 681经朔误差/h0.219 842-0.013 800-0
38、.113 830-0.026 9900.309 3720.128 488拟合线截距0.176 0-0.011 8-0.095 5-0.014 80.293 30.174 0表3纪元历 等6部历法的朔望月常数与经朔误差各部历法经朔误差的差异很大，淳熙历 最大，达到了18 min，而 崇天历 的还不到1 min。与定朔误差对比，可以得到经朔误差到定朔的影响有多大；由定朔误差减去经朔误差就是由其他因素所引起的误差。6部历法中，淳熙历 的其他误差仅有5 min，而纪元历 崇天历 和 统元历 的可以达到30 min，这三部历法的经朔误差为负，其他历法的经朔误差为正。经朔误差是影响定朔的一个主要因素，但是
39、在经朔误差较低的情况下，定朔仍然很高，如 淳熙历，而经朔误差很高的 崇天历 的定朔误差却要大一些。这就说明，其他因素的作用仍是十分重要的。5 总结与展望定朔算法是一个构造简单但很重要的算法。中国古代推朔的基本思路是先计算出一个平均值，再设计算法修正平均值，得到最终计算值。宋代历法推算定朔的精度是：平均误差为4 min，误差绝对值的平均值为 0.391 5 h，约合 23.5 min，最大误差为2.498 9 h，满足当时推朔精度要求的应该在九成以上。就每部历法而言，会元历 为最劣，开禧历 为最佳。进一步指出，当历法观测地与使用地不同时，古人做不到客观地判断定朔精度。于是，我们分别使66咸阳师范
40、学院学报第29卷用现代人和古人两种视角去评价历法的精度。那么，在古人看来，纪元历 与 统元历 的定朔计算精度都比较低，而 淳熙历 和 会元历 的精度则明显提高。这些历法中，仍然以 开禧历 为最优。虽然已经掌握了宋代11部历法的定朔计算精度，但我们目前还没有完全弄清楚 明天历 应天历 仪天历 这3部历法推算定朔的术文和常数。至于其它历史时期的历法定朔精度问题更是没有全面展开，在程序化和精度分析方面还有很多工作要做。参考文献：1刘金沂，赵澄秋麟德历 定朔计算法M/薄树人，刘金沂.中国天文学史文集：第3集.北京：科学出版社，1984：36-882薄树人纪元历 解M/薄树人文集.合肥：中国科学技术大学
41、出版社，2003：369-4473张培瑜授时历 定朔日躔及历书推步M/薄树人，刘金沂.中国天文学史文集：第6集.北京：科学出版社，1994：30-584薮内清.隋唐历法史研究M.东京：三省堂，1944：61-825李勇，张培瑜.中国古历定朔推步综述J.天文学进展，1996，14（1）：66-766尚晓清，曲安京大衍历 定朔算法及程序说明J.西北大学学报：自然科学版，1999，29（3）：193-1967张培瑜，王桂芬，陈月英，等宣明历 定朔计算和历书研究J.紫金山天文台台刊，1992，11（2）：121-1558滕艳辉,王鹏云纪元历 定朔算法模型及分析J.西北大学学报：自然科学版，2008，40（5）：855-8589滕艳辉，王鹏云纪元历 等8部宋代历法的定朔推步及精度分析J.中国科技史杂志，2009，30（1）：55-6410滕艳辉，唐 泉开禧历 定朔计算精度及分析J.自然科学史研究，2011，30(1)：19-2711张培瑜.三千五百年历日天象M.郑州：河南教育出版社，199012中华书局编辑部.历代天文等志汇编M.北京：中华书局，1976：268813陈美东.古历新探M.沈阳：辽宁教育出版社，1995：221-22714孙小淳.从“里差”看地球、地理经度概念之传入中国J.自然科学史研究，1998，17（4）：304-311第4期滕艳辉：宋代定朔计算精度分析67