我国古代历法中应用的内插法

《我国古代历法中应用的内插法》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、我国古代历法中应用的内插法梅荣照　　我国古代历法中应用的内插法　　已知函数f（x）在自变量是x1，x2，……xn时的对应值是f（x1）。f（x2），……f（xn），求xi和xi＋1之间的函数值的方法，称作内插法。如果xn是按等距离变化的，称自变数等间距内插法；如果xn是接不等距离变化的，称自变数不等间距内插法。例如f（x）＝x3，当x＝1，2，3，4，5，……时x3＝1，8，27，64，125，……求x＝4.26时x3＝（4.26）3的值，就可以应用等间距内插公式。等间距内插法的一般公式是：　　叫n级差分。

2、从n级差分的定义容易得到，当f（x）是一次函数时，二级差分是0；f（x）是二次函数时，三级差分是0；f（x）是n次函数时，n＋1级差分是0。　　我国古代历法工作者，为了制定一个好的历法，很早就应用内插法的公式。　　朔、望跟制定历法和计算日、月食有密切关系。怎样确定合朔的准确时刻，一直是历法中一个重要的项目。根据一个朔望月的平均日数来确定会朔时刻，叫做"平朔"。例如《古四分历》、西汉《太初历》和《三统历》就是用平朔来确定会朔时刻的。这种方法所取的分数不是过大，就是过小，年长日久，合朔发生的时刻不是提前就是推后

3、。最初历法工作者以为这种提前和推后是分数部分不精密造成的，他们曾经企图调整分数部分以求得准确的平均朔望月日数来解决这个问题。但是由于月球运动的速度是随时间不同而变化的，因此平朔的方法不可能得到真正的合朔时刻。东汉初年，公元一世纪初，天文学家发现了月球绕地球运行的变速运动。东汉建安十一年（公元206年），刘洪在《乾象历》中首次提出用一次内插的方法来确定合朔的时刻。确定合朔发生的真正时刻叫做"定朔"。要根据月球运行的速度来计算合朔，不仅需要知道每日月球的运行度数，还必须知道小于一日的月球运行度数。刘洪测出月球在

4、一个近点月（月球从近地点出发绕地球运行一周又回到近地点的时间间隔）里每日运行的度数。设日数是n，n日共行的度数是f（n），对n＋s（s＜1）日月球运行的度数，刘洪应用下列一次内插公式　　f（n＋s）＝f（n）＋s△　　进行计算，其中△是一级差分f（n＋1）－f（n）。刘洪以后，三国时期的杨伟，南北朝时期的何承天、祖冲之都是用这个公式计算月行度数的。因为月球运动速度一日之内就变化很大，f（n）不是一次函数，△2f（x）不等于0，因此上述公式只能得到不很精密的近似值。　　随着天文观测技术的发展，天文学家认识到太

5、阳运行速度的变化也影响到合朔时刻，因此提出测量日、月、五星视行度数的更加精密的年插公式。隋文帝开皇二十年（公元600年），天文学家刘焯在《皇极历》中提出一个推算月、月、五星视行度数的等间距二次内插公式：　　求太阳的视行度数时，ι是一节气的日数，求月行度数时，ι=1日。刘焯的公式虽比以前精密得多，但是由于节气日数l实际上不是按等间距变化，日、月、五星也不是作等加速运动（就是说三级差分不等于0），因此仍然存在缺点。这两个问题分别由唐代一行和元代郭守敬等解决了。唐玄宗开元十五年（公元727年），一行在他的《大衍历

6、》中提出了不等间距的二次内插公式：　　当ι1＝ι2时，和刘悼的等间距二次内插公式相同。　　元世祖至元十八年（公元1281年），郭守敬等的《授时历》在计算日、月、五星视行度数中考虑了日、月、五星运行的不等速运动情况，认为距离是时间的三次函数。不过《授时历》没有求出三次内插公式，而是用差分表来解决这个问题的。例如计算太阳每日的运行度数：《授时历》把冬至到春分（共88.91日）这一象限分成六段，测出每段太阳的实际运行度数，就可以算出以段为等间距的差分表。从表中知道，三级差分都相等而四级差分等于0，因此考虑　　f(

7、t)=d+at+bt2+ct2。　　实际上，f（0）＝0（第0段的运行度数是0），可见d＝0.所以　　f(t)=at+bt2+ct3。　　这样就可以变三次函数为二次函数　　应用二次内插公式便可以算出F（t）的具体表达式，从而得到　　f（t）＝tF（t）＝513.32t－2.46t2－0.0031t3。　　令t＝0，1，2，3，按差分的定义便可以求出f（0），△1，△2，△3（△4＝0），继续按差分定义，用加减法就可以得出以日为等间距的差分表①。　　从公元三世纪到公元十三世纪一千多年中，内插法在天文学家和数学

8、家中得到广泛的应用，除上面提到的以外，如唐代傅仁均的《戊寅历》（公元619年）、李淳风的《麟德历》（公元664年）、徐昂的《宣明历》（公元822年）、宋代秦九韶《数书九章》（公元1247年）的"缀术推星"、元代朱世杰《四元玉鉴》（公元1303年）的"如象招数"等等。特别需要提出的是朱世杰的工作。　　朱世杰"如象招数"②的问题是已知一级差分　　这里△、△2、△3、△4分别代表各级差分的第一个差分。朱世杰知道，公式的