中国古代的日食时差算法.doc

石河子大学学报 (自然科学版) Journal of Shihezi University (Natural Science) 第 23 卷 第 4 期 2005 年 8 月 Vol . 23 No . 4 Aug. 2005 文章编号 :100727383 (2005) 0420416206 中国古代的日食时差算法 曲安京 , 唐 泉 (西北大学数学系 , 陕西 西安 710069) 摘要 : 因月亮视差引起的对日食食甚时刻的修正值 ,即定朔时刻与日食食甚时刻之差 ,中国古代称为时差 。本文讨 论了中国古代时差算法的沿革及其算法模型 ,并与理论算法作了比较 ,分析了其精度误差 。从理论上来说 ,时差算 法是关于太阳的时角和黄经的二元函数 ,但是中国古代的时差算法只考虑了时角因素 ,而完全没有考虑到太阳的 黄经变化对时差的影响 ,因此中国古代的时差算法是一个不完全的算法 。本文的讨论结果表明 ,中国古代的时差 算法虽然是一个有缺陷的算法 ,但其结果和理论结果还是比较接近 ,特别是和二至时的理论结果相比 ,误差只达到 2min 左右 。 关键词 : 时差 ; 食甚时刻 ; 定朔时刻 ; 精度分析 中图分类号 : P12091. 2 文献标识码 : A 地球赤道半径 R0 = 6378. 16km ,则月亮的地平视差 H0 = R0 ΠR = 0. 0165925 = 0 °. 9507 。 时差的天文意义与数值模型 由于月亮视差的影响 ,日食的视食甚时刻 ,并不 是定朔时刻 。从《宣明历》开始设计的日食三差中的 时差算法 ,就是为了修正这个差值 ,具体地说 ,中国 古代日食计算中的时差 ,就是日食的视食甚时刻与 定朔时刻之差 。 如图 1 所示 , 假设 E 为地心 , O 为观测者 , EO = R0 表示地球的赤道半径 , EM = R 表示地球到天 体 M 的距离 , Z 为天顶 。则 ∠M EZ = z 表示天体的 天顶距 , ∠OM E = H 表示其视差 。由于 ∠MO Z = H + z ,而 R sin z = tan ( H + z) ( Rcos z - R0 ) , 且 R 远大于地球半径 R0 ,令 H0 = R0 ΠR ,则有 , 1 图 1 天体的视差 日食时差算法的天文模型 ,可以借助图 2 进行 说明 :假设点 Q 、Q’分别表示黄道与白道 、视白道的 交点 , MM′A 与 MM″S 分别表示过月亮 M 的地平经 圈与黄经圈 。 点 M 、S 分别表示定朔时刻的月亮 、太 阳 。由 于月亮的视差 ,此时 ,我们在地球上看到的月亮的实 际位置为点 M′, 令 MM′= H = H0 sinz , 此数即为月 亮的视差 。点 M″、S 分别表示日食的食甚时刻的月 亮 、太阳 ,月亮从点 M′运动到点 M′的时?? ,就是日 R0 sin z R0 sin z H≈tan H = R - R cos Z≈ R ≈ H0 sin z 。 0 由此可见 ,天体的视差是其天顶距的函数 。当 天体位于天顶时 ,视差 H = 0 ; 当天体出现在地平线 上时 ,视差 H = H0 ,达到最大值 。 由于赤道半径相对于日地平均距离 ( 1 个天文 单位) 非常小 ,因此 ,太阳视差也是很小的 ,其地平视 差仅为 8″18 ,可以忽略不计 。对日食产生影响的 ,主 要是月亮的视差 。按月地平均距离 R = 384401km , 收稿日期 : 2005204225 基金项目 :国家自然科学基金 作者简介 : 曲安京 (1962) ,男 ,教授 ,从事精密科学史研究 。  食的时差 。 令 I = ∠MQS 表示黄白道夹角 ,这个数值比较 小 ,因此 ,三角形 MM′M″可以近似地看作直角三角 形 。令 F = ∠MAS 表示过月亮 M 的地平经圈与黄 道的交角 ,则 F = ∠MM′M″+ I≈ ∠MM′M″, 于是 , M′M″≈ Hcos F = H0 sin zcos F , 关数据的计算来推算Δt 值 。 如图 3 所示 ,点 P 表示北赤极 ,点 N 表示黄赤 道交点 (如春分点) ,点 S 表示太阳 , h = ∠S PZ 表示 其时角 ,点 Z 表示天顶 。令 C = ∠BSN , = ∠PS Z , 则角 F = C + C′。又设 δ= BS 、λ= NS 、z = S Z 分别 表示太阳 S 的赤纬 、黄经 、天顶距 ,ε = ∠BNS 为黄 赤大距 ,φ = 90°- PZ 表示观测点的地理纬度 ,则根 据球面三角的关系 ,可以得到如下的方程组 : sinδ= sinλsinε; cos z = sin