JavaScript计算公元前日期及格里高利历切换的儒略日

一、儒略日的本质与两大历法核心边界

儒略日是由法国学者约瑟夫·斯卡利杰（Joseph Scaliger）于1583年提出的一种连续计时系统。它定义为自格林尼治平午（世界时 12:00:00）公元前4713年1月1日（儒略历）起算所经过的连续天数。儒略日的优势在于没有任何年份、月份、闰秒或历法变迁的阶跃，任何两点之间的时间间隔都可以直接通过加减天数得出。

常规的时间库和 JavaScript 内置的 Date 对象无法胜任历史天文的儒略日计算，原因就在于无法自动处理以下两个历法边界：

• 天文学纪年与历史纪年的偏差：历史学中不存在“公元0年”，公元1年（1 AD）的前一年直接被称为公元前1年（1 BC）。而数学和天文学为了保证加减连续性，将公元前1年定义为第0年（Year 0），公元前2年为第 -1 年，公元前 N 年在天文坐标下应表示为 -(N - 1) 年。
• 1582年格里高利历切换的10天跳变：儒略历因每4年固定置闰一次，年均长度为365.25天，与回归年（约365.2422天）相比每年累积约11分钟误差。到1582年时春分点已提前10天。教皇格里高利十三世下令：1582年10月4日（星期四）的次日直接调整为1582年10月15日（星期五）。这意味着1582年10月5日至14日这10个日历日期在历史上根本不存在，但儒略日必须保持数学上的严格单调递增。

图1：天文学年代坐标系统与1582年教皇格里高利历改历分界模型

二、儒略日计算的标准推导公式

在天文学计算中，广泛采用比利时天文算法专家让·梅乌斯（Jean Meeus）在《天体测量计算方法》（Astronomical Algorithms）中推荐的标准算法。为确保1月和2月能平滑承接前一年的闰年周期，算法首先将1月和2月视作上一年的第13月和第14月。

1. 月份平移与基础修正

设输入公历年月日为 Y, M, D，时间（时分秒）折算为当天的日小数 UT：

UT = hour / 24 + minute / 1440 + second / 86400
DayFraction = D + UT

如果 M ≤ 2，则执行：

Y = Y - 1
M = M + 12

这种平移使得闰日（2月29日）恰好成为整年的最后一天，极大简化了年内积日的推算。

2. 历法分支判断与世纪修正项

我们需要根据原始日期确定其处于哪一种历法：

• 若日期在 1582年10月15日 00:00:00 及之后（采用格里高利历），需要计算世纪闰年修正项 B：

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

• 若日期在 1582年10月4日 23:59:59 及之前（采用儒略历），则不需要世纪闰年扣除：

  B = 0

3. 核心积日多项式

儒略日的标准连续计算公式为：

JD = Math.floor(365.25 * (Y + 4716)) + Math.floor(30.6001 * (M + 1)) + DayFraction + B - 1524.5

公式中的各常数含义如下：

• 365.25：儒略历平均年长，配合 Math.floor 每4年自动吸收一个闰日。
• 4716：将计算原点平移至公元前4713年，使得公元前日期的年数在加和后依然保持正数。
• 30.6001：用于拟合3月至次年2月各月份天数变化的多项式系数（采用 30.6001 而不是 30.6 是为了抵消计算机浮点数表示 30.6 时可能产生的下溢舍入误差）。
• -1524.5：基准对齐常数。因为儒略日以正午12:00（0.5日）为日界起点，减去 1524.5 可将午夜 00:00 的日小数准确对应至 .5。

图2：儒略日核心算法与历法分支判断流程图

三、闰年规则与历法突变的数学本质

历法切换的根本原因在于回归年天数与历法平均年长的微小差额。在儒略历中，只要年份能被4整除即为闰年，导致每400年中出现了100个闰年。而格里高利历引入了“整百年必须能被400整除才闰”的规定，使400年中的闰年数减少至97个，将历年平均长度修正为 365.2425 天。

下表清晰对比了两个历法体系在世纪年上的置闰差异：

图3：儒略历与格里高利历置闰规则及误差补偿机理

四、JavaScript 完整算法实现

在编写 JavaScript 算法时，有三个致命的语法与实现陷阱必须严格规避：

1. 严禁使用 Math.trunc() 或按位或 | 0 代替 Math.floor()：在处理公元前（负数年）时，Math.floor(-1.5) 返回 -2，而 Math.trunc(-1.5) 返回 -1。梅乌斯公式推导严格依赖向负无穷方向取整，若使用向零取整将导致公元前积日出现巨大跳变。
2. 非法日期的防御拦截：1582年10月5日到10月14日之间为历史空窗期，程序应当在此处主动检测并报错或返回警示，避免输入无意义日期。
3. 公元前标志的统一重载：支持传入负数年份（天文学纪年）或通过显式选项参数指定为历史 BC。

以下是具备工业级健壮性的 JavaScript 完整实现代码：

/**
 * 计算指定日期的儒略日 (Julian Day)
 * 
 * @param {Object} options 输入参数对象
 * @param {number} options.year 年份（天文纪年可传入0或负数；若isHistoricalBC为true，则传入正整数如4713表示前4713年）
 * @param {number} options.month 月份 (1 - 12)
 * @param {number} options.day 日期 (1 - 31)
 * @param {number} [options.hour=0] 时 (0 - 23)
 * @param {number} [options.minute=0] 分 (0 - 59)
 * @param {number} [options.second=0] 秒 (0 - 59.999)
 * @param {boolean} [options.isHistoricalBC=false] 是否为历史年代学纪年 (无0年，1 BC -> 0年，2 BC -> -1年)
 * @param {string} [options.calendar='auto'] 'auto' (1582年自动分界), 'julian' (强制纯儒略历), 'gregorian' (强制纯外推格里高利历)
 * @returns {number} 计算得到的儒略日数 (精确到浮点数天)
 */
function calculateJulianDay(options) {
  let {
    year,
    month,
    day,
    hour = 0,
    minute = 0,
    second = 0,
    isHistoricalBC = false,
    calendar = 'auto'
  } = options;

  // 1. 历史公元前纪年转换为天文学数学纪年
  let y = year;
  if (isHistoricalBC) {
    if (y <= 0) {
      throw new RangeError("历史公元前纪年中的年份参数必须为正整数 (例如 1 代表 1 BC)");
    }
    y = -(y - 1);
  }

  let m = month;
  const d = day;

  // 2. 检查 1582 年历法跳跃的非法日期区间
  if (calendar === 'auto' && y === 1582 && m === 10) {
    if (d > 4 && d < 15) {
      throw new RangeError("1582年10月5日至10月14日在格里高利历切换中被抹除，属于历史非法日期！");
    }
  }

  // 3. 计算日内小数部分 (世界时 UT 分数)
  const dayFraction = d + (hour + minute / 60 + second / 3600) / 24;

  // 4. 若为 1 月或 2 月，视作上一年的 13 月或 14 月
  if (m <= 2) {
    y = y - 1;
    m = m + 12;
  }

  // 5. 判断历法分支并计算世纪修正项 B
  let b = 0;
  let isGregorian = false;

  if (calendar === 'gregorian') {
    isGregorian = true;
  } else if (calendar === 'julian') {
    isGregorian = false;
  } else {
    // 自动判定：1582年10月15日及以后为格里高利历
    // 注意：这里的原始输入是 options.year, options.month, options.day
    const origY = isHistoricalBC ? -(year - 1) : year;
    if (
      origY > 1582 ||
      (origY === 1582 && month > 10) ||
      (origY === 1582 && month === 10 && day >= 15)
    ) {
      isGregorian = true;
    }
  }

  if (isGregorian) {
    // 必须使用 Math.floor 确保负数年份正确向下取整
    const a = Math.floor(y / 100);
    b = 2 - a + Math.floor(a / 4);
  }

  // 6. 套用梅乌斯标准多项式计算连续积日
  // 注意：所有向下取整必须严格使用 Math.floor
  const jd = Math.floor(365.25 * (y + 4716)) +
             Math.floor(30.6001 * (m + 1)) +
             dayFraction +
             b -
             1524.5;

  return jd;
}

五、测试用例验证与结果解析

为验证上述算法在关键天文分界和公元前边界的正确性，设计以下覆盖典型历史节点的基准测试用例：

// 测试 1：天文标准元元 J2000.0 (2000年1月1日 12:00:00 UT)
const testJ2000 = calculateJulianDay({
  year: 2000, month: 1, day: 1, hour: 12, minute: 0, second: 0
});
console.log("J2000.0 期望值: 2451545.0, 实际值:", testJ2000);

// 测试 2：儒略历最后一天 (1582年10月4日 00:00:00 UT)
const testJulianEnd = calculateJulianDay({
  year: 1582, month: 10, day: 4, hour: 0, minute: 0, second: 0
});
console.log("1582-10-04 期望值: 2299159.5, 实际值:", testJulianEnd);

// 测试 3：格里高利历第一天 (1582年10月15日 00:00:00 UT)
const testGregorianStart = calculateJulianDay({
  year: 1582, month: 10, day: 15, hour: 0, minute: 0, second: 0
});
console.log("1582-10-15 期望值: 2299160.5, 实际值:", testGregorianStart);
console.log("历法断层后相差天数 (应恰好为 1.0 天):", testGregorianStart - testJulianEnd);

// 测试 4：公元前1年 (天文0年) 1月1日 00:00:00 UT
const test1BC = calculateJulianDay({
  year: 1, month: 1, day: 1, hour: 0, minute: 0, second: 0, isHistoricalBC: true
});
console.log("1 BC 01-01 期望值: 1721423.5, 实际值:", test1BC);

// 测试 5：儒略日计时起点 (历史 4713 BC 1月1日 12:00:00 UT)
const testEpoch = calculateJulianDay({
  year: 4713, month: 1, day: 1, hour: 12, minute: 0, second: 0, isHistoricalBC: true
});
console.log("儒略历原点 4713 BC 期望值: 0.0, 实际值:", testEpoch);

测试运行后输出表明，1582年10月4日到15日虽然日期数字骤增了11，但算出的儒略日数从 2299159.5 递增到 2299160.5，其差额严格等于 1.0 天，保证了时间轴的完全连续。同时，4713 BC 算得 0.0，完美契合斯卡利杰的历法原点定义。

图4：测试矩阵与工程避坑核心原则架构总结

六、实践中的精度与时区考量

总结

利用 JavaScript 计算跨越两千多年的儒略日，关键在于理顺纪年体系的对应关系和1582年改历时的分界机制。通过将历史公元前年份映射为天文学负数坐标、在1-2月执行跨年平移，并严格采用 Math.floor 进行世纪分界修正，即可构建出一套高精度、全天候适用的连续计时函数，为天文学模拟、历史天象反推与历法转换提供可靠的底层支持。