儒略日(JDN)是从公元前4713年1月1日中午UTC起算的整数天数,C++标准库未提供因其专注公历且不视其为通用需求,需用Meeus公式等整数算法手动计算。

什么是儒略日,为什么 C++ 标准库不直接提供
儒略日(Julian Day Number, JDN)是从公元前 4713 年 1 月 1 日中午(UTC)开始连续计数的天数,常用于天文计算和跨历法时间比对。C++ 标准库(包括 <chrono> 和 <ctime>)只处理公历(Gregorian calendar),且不暴露儒略日换算逻辑——它不认为这是“通用时间需求”,所以得自己算或借助可靠公式。
别用 std::time() + 手动天数累加,误差会随年份放大;也别依赖系统时区转换再折算,UTC 正午这个基准点很容易搞错。
用已验证的整数算法直接算 JDN(推荐)
最稳妥的方式是用由 Jean Meeus 在《Astronomical Algorithms》中给出的整数公式(适用于公历 1582 年 10 月 15 日之后,即格里高利历启用后;若需支持更早日期,得切换儒略历公式,但绝大多数工程场景不需要)。
- 输入必须是 UTC 时间下的年、月、日(注意:儒略日以 UTC 中午为日界,不是午夜)
- 若你只有本地时间,先用
std::gmtime()或 C++20 的std::chrono::utc_clock转成 UTC,再取日期部分 - 月份按 1–12,年份可正可负(BC 年份用 0 表示 1 BC,-1 表示 2 BC,依此类推)
示例(C++11 起可用):
立即学习“C++免费学习笔记(深入)”;
int julian_day_number(int year, int month, int day) {
if (month <= 2) {
year--;
month += 12;
}
int a = year / 100;
int b = 2 - a + a / 4;
return static_cast<int>(365.25 * (year + 4716)) +
static_cast<int>(30.6001 * (month + 1)) +
day + b - 1524;
}注意:365.25 和 30.6001 是浮点近似,但对公元 1 年至今所有日期都给出精确整数结果。若追求纯整数运算(避免浮点舍入风险),可用等价整数形式:
int jdn = (1461 * (year + 4800 + (month - 14) / 12)) / 4 +
(367 * (month - 2 - 12 * ((month - 14) / 12))) / 12 -
(3 * ((year + 4900 + (month - 14) / 12) / 100)) / 4 +
day - 32075;获取“当前”JDN 的陷阱:UTC 中午 vs 系统时间
儒略日从 UTC 每日 12:00 开始,不是 00:00。如果你直接拿 std::time(nullptr) 得到秒数,再除以 86400 取整,会得到“儒略日序数(Julian Date, JD)”的小数部分错误——它把日界设在了午夜。
组合式C++代码评审方案,融合静态分析、AI推理、多轮迭代评审和C++专项检查,适用于PR审查、增量代码审查、全项目评审和代码质量评分,触发词包括review cpp、cpp代码评审、C++review、代码审查。
正确做法是:先拿到当前 UTC 时间的 std::tm,然后构造当天 12:00 UTC 对应的秒数,再转成 JDN。
- 用
std::gmtime()获取 UTC 的tm结构体 - 把
tm_hour设为12,tm_min和tm_sec设为0,再用std::mktime()(注意:它假定输入是本地时区!所以得先转成 UTC 秒数) - 更安全的做法(尤其 C++20):用
std::chrono::utc_clock::now()→ 转sys_days→ 加 12h → 再调用上面的整数函数
简单稳健方案(C++11):
auto now = std::time(nullptr); auto tm_utc = *std::gmtime(&now); tm_utc.tm_hour = 12; tm_utc.tm_min = tm_utc.tm_sec = 0; // 注意:mktime 期望本地时区,所以得手动补偿时区偏移 long offset_sec = std::mktime(&tm_utc) - now; auto jdn = julian_day_number(tm_utc.tm_year + 1900, tm_utc.tm_mon + 1, tm_utc.tm_mday);
其实只要确保传入的年月日是 UTC 当日即可,julian_day_number 函数本身不关心时间,只认日期。所以最简方式就是:
auto t = std::time(nullptr); auto tm = *std::gmtime(&t); int jdn = julian_day_number(tm.tm_year + 1900, tm.tm_mon + 1, tm.tm_mday);
——因为儒略日定义的是“自公元前 4713-01-01T12:00 UTC 起的整日数”,只要日期对,JDN 就对;小数部分(JD)才需要精确到秒,而 JDN 是整数。
Windows 和嵌入式平台要注意的细节
某些 Windows CRT 实现中,std::gmtime() 在 year < 1970 时可能返回 nullptr;嵌入式 libc(如 newlib)可能不支持负年份或 BC 年份。这时不要硬传 tm_year = -4713,而是限定使用范围(比如只处理 1900–2100 年)。
- 检查
std::gmtime返回值是否为 null,否则崩溃 - 避免用
std::put_time或strftime解析 BC 年份——它们通常不支持 - 如果项目允许第三方依赖,
date.h(Howard Hinnant’s date library)提供floor<days>(utc_clock::now()) + days{12}后转 JDN,更健壮
JDN 计算本身很轻量,但“当前”二字背后藏着时区、日界、历法三重校验。别图省事跳过 UTC 转换,也别信网上随手搜的“一行代码 JDN”。

















