JavaScript实现公历日期转儒略日(JDN)算法代码

在天文学计算、卫星轨道模拟以及大跨度历史日期推演中，传统公历的闰年交替和大小月差异会显著增加跨周期天数差的计算复杂度。儒略日（Julian Day Number, JDN）通过将时间定义为一个连续累加的整数天数计数体系，从根本上消除了历法断代与月份规则的干扰。本文将完整介绍基于经典 Fliegel-Van Flandern 与 Jean Meeus 改进模型的算法原理，并给出可在浏览器与 Node.js 环境中稳定运行的纯 JavaScript 实现代码。

一、什么是儒略日与公历转换的本质问题

儒略日数（Julian Day Number，通常缩写为 JDN）是以格林威治标准时间公元前 4713 年 1 月 1 日中午 12:00（平午）为原点（第 0 天）开始连续计算的整数天数系统。而包含当天日内小数部分的时间表示则被称为儒略日（Julian Date，缩写为 JD）。

实现公历到 JDN 转换的核心难点在于两点：其一是公历每个月的天数并不均匀（如 31、28/29、31、30 等），导致月份累计天数难以用简单的连续乘除法表示；其二是 2 月份作为闰月位于年初附近，如果发生闰年，后续所有月份的累计天数都会产生位移偏移。

图1：月份前移逻辑与 30.6001 系数拟合原理

二、JDN 核心算法推导与计算管道

历史上著名的天文学家 Jean Meeus 在其专著《Astronomical Algorithms》中总结了公历转换为儒略日的精确闭式公式。整个计算管道可以划分为清晰的四个阶段：

1. 阶段一：年月重映射
   若输入的月份 $M \le 2$，令年份 $Y' = Y - 1$，月份 $M' = M + 12$；若 $M > 2$，则 $Y' = Y$，$M' = M$。
2. 阶段二：格里高利历世纪闰年修正项（$B$）
   公历规定：公元年份能被 4 整除为闰年，但逢百年不闰（不能被 100 整除），逢四百年又闰（能被 400 整除）。因此，世纪修正项计算为：

   const A = Math.floor(Y' / 100);
   const B = 2 - A + Math.floor(A / 4);

   注：如果是 1582 年 10 月 4 日及以前的纯儒略历，此项修正值 $B = 0$。现代系统处理外推公历（Proleptic Gregorian Calendar）时均保留该修正项。
3. 阶段三：基准年日积分累加
   以公元前 4716 年作为基准起点，每四年增加一个闰日，平均每年 365.25 天：

   const term1 = Math.floor(365.25 * (Y' + 4716));

4. 阶段四：月份累计与偏移常数合并
   加入月份累计天数、当月日（$D$）、世纪修正项 $B$ 并减去基准原点偏移量 1524：

   const JDN = term1 + Math.floor(30.6001 * (M' + 1)) + D + B - 1524;

图2：JDN 算法四步执行管道与数据流向

三、JavaScript 完整代码实现

以下是符合现代 JavaScript（ES6+）规范的完整实现代码。代码模块包含三个层次的函数：

• gregorianToJDN(year, month, day)：针对常规年月日的纯整型 JDN 计算；
• gregorianToJD(year, month, day, hour, minute, second)：支持时、分、秒日内时间分数的浮点数 JD 计算；
• dateToJD(date)：支持接收原生 JavaScript Date 对象，自动转换时区并提取 UTC 分量。
• jdnToGregorian(jdn)：逆运算函数，将儒略日数还原为公历年月日，用于闭环自测。

/**
 * 公历日期转儒略日数 (Julian Day Number, JDN)
 * 采用 Jean Meeus 经典天文学算法
 * 
 * @param {number} year  公历年份 (整数，如 2024)
 * @param {number} month 公历月份 (1 - 12)
 * @param {number} day   公历日 (1 - 31)
 * @returns {number} 整数儒略日数 (以当日正午 12:00 UT 为基准)
 */
function gregorianToJDN(year, month, day) {
  // 必须确保参数为整数类型
  let y = Math.floor(year);
  let m = Math.floor(month);
  const d = Math.floor(day);

  // 边界校验
  if (isNaN(y) || isNaN(m) || isNaN(d) || m < 1 || m > 12 || d < 1 || d > 31) {
    throw new RangeError("无效的公历日期参数: 请输入合法的年月日数值");
  }

  // 阶段一：1月与2月前移至前一年作为13月与14月
  if (m <= 2) {
    y -= 1;
    m += 12;
  }

  // 阶段二：格里高利历世纪闰年修正项 B
  const a = Math.floor(y / 100);
  const b = 2 - a + Math.floor(a / 4);

  // 阶段三 & 四：年分天数累计 + 月份天数累计 + 当月日 + 修正值 - 1524
  const jdn = Math.floor(365.25 * (y + 4716)) +
              Math.floor(30.6001 * (m + 1)) +
              d + b - 1524;

  return jdn;
}

/**
 * 公历日期带时分秒转儒略日 (Julian Date, JD)
 * 包含日内连续时间分数，正午 12:00:00 UT 对应 .000000
 * 
 * @param {number} year   年份
 * @param {number} month  月份 (1 - 12)
 * @param {number} day    日期 (1 - 31)
 * @param {number} hour   小时 (0 - 23，世界时 UT)
 * @param {number} minute 分钟 (0 - 59)
 * @param {number} second 秒数 (0 - 59.999...)
 * @returns {number} 浮点数儒略日 (JD)
 */
function gregorianToJD(year, month, day, hour = 12, minute = 0, second = 0) {
  const jdn = gregorianToJDN(year, month, day);

  // 计算相对于正午 12:00:00 UT 的日内时间偏移比例
  const dayFraction = (hour - 12) / 24 + minute / 1440 + second / 86400;

  return jdn + dayFraction;
}

/**
 * 直接将 JavaScript 原生 Date 对象转换为儒略日 (JD)
 * 严格基于 UTC 国际标准时间，规避本地时区 offset 的干扰
 * 
 * @param {Date} date 原生 JavaScript 日期对象
 * @returns {number} 浮点数儒略日 (JD)
 */
function dateToJD(date) {
  if (!(date instanceof Date) || isNaN(date.getTime())) {
    throw new TypeError("必须传入有效的 Date 对象实例");
  }

  const year = date.getUTCFullYear();
  const month = date.getUTCMonth() + 1; // getUTCMonth 返回 0-11
  const day = date.getUTCDate();
  const hour = date.getUTCHours();
  const minute = date.getUTCMinutes();
  const second = date.getUTCSeconds() + date.getUTCMilliseconds() / 1000;

  return gregorianToJD(year, month, day, hour, minute, second);
}

/**
 * 儒略日数 (JDN) 逆运算还原为公历年月日
 * 用于双向校验
 * 
 * @param {number} jdn 整数儒略日数
 * @returns {{year: number, month: number, day: number}} 公历日期对象
 */
function jdnToGregorian(jdn) {
  const l = Math.floor(jdn) + 68569;
  const n = Math.floor((4 * l) / 146097);
  const l1 = l - Math.floor((146097 * n + 3) / 4);
  const i = Math.floor((4000 * (l1 + 1)) / 1461001);
  const l2 = l1 - Math.floor((1461 * i) / 4) + 31;
  const j = Math.floor((80 * l2) / 2447);
  const day = l2 - Math.floor((2447 * j) / 80);
  const l3 = Math.floor(j / 11);
  const month = j + 2 - 12 * l3;
  const year = 100 * (n - 49) + i + l3;

  return { year, month, day };
}

四、基准测试与历史关键节点校验

为了确保算法的准确性，必须针对天文计算中公认的标准基准点（Epoch）以及极易出现边界错误的闰年、世纪年份进行全面核验。

图3：多时间体系与关键天文历元对照

五、工程实践中的注意事项与避坑指南

在将该算法整合到实际业务工程中时，开发者往往会遇到时区偏差、JavaScript 浮点数截断以及数值溢出等常见问题。理解以下要点可帮助避免隐蔽的生产环境缺陷：

图4：工程实践细节与算法可靠性要点

六、总结

儒略日（JDN）通过将不规则的年、月、日历法结构映射为一条平直的时间数轴，极大地简化了日期跨度计算与天体运动周期的推导。通过本文的 JavaScript 实现，开发者可以在仅依赖 Math.floor 的前提下，以极高的性能与数学严密性完成双向转换。在实际项目中，请务必注意时区隔离（统一使用 UTC）以及避免使用 32 位位运算符替代数学舍入函数，从而确保算法在各种极端日期与负年份场景下的精确与稳健。