张若昀和杜鹃吻戏视频:儒略日计算

参考：http://korea.3dapple.com/user1/216/archives/2004/840.shtml
转贴原文——

一个计算儒略日的c语言函数

根据Fundamental Astronomy上介绍的算法写的。

float jde(int Y,int M,int D,int hour,int min,int sec)/*计算儒略日历书时；*/
{
long int f,g,mid1,mid2;
float J,JDE,A;
if(M>=3)
{
f=Y;
g=M;
}
if(M==1||M==2)
{
f=Y-1;
g=M=12;
};
mid1=floor(365.25*f);
mid2=floor(30.6001*(g+1));
A=2-floor(f/100)+floor(f/400);
J=mid1+mid2+D+A+1720994.5;
JDE=J+hour/24+min/1440+sec/86400;
return JDE;

};

求出给定年(I)，月(J)，日(K)的儒略日:
儒略日= K - 32075 + 1461 * (I + 4800 + (J-14)/12)/4+367*(J-2-(J-14)/12*12)/12-3*((I+4900+(J-14)/12)/100)/4

例如1986年12月25日的儒略日(I=1986，J=12，K=25)=2446791

儒略日（Julian day）是指由公元前4713年1月1日，协调世界时中午12时开始所经过的天数，多为天文学家采用，用以作为天文学的单一历法，把不同历法的年表统一起来。

儒略日是一种不用年月的长期纪日法，简写为JD。是由法国学者Joseph Justus Scliger（1540–1609）在1583年所创，这名称是为了纪念他的父亲——意大利学者Julius Caesar Scaliger（1484–1558）。

儒略日的起点订在公元前4713年（天文学上记为 -4712 年）1月1日格林威治时间平午（世界时12:00），即JD 0 指定为 4713 B.C. 1月1日12:00 UT到4713 B.C. 1月2日12:00 UT的24小时。每一天赋予了一个唯一的数字，顺数而下，如：1996年1月1日12:00:00的儒略日是2450084。这个日期是考虑了太阳、月亮的运行周期，以及当时收税的间隔而订出来的。Joseph Scliger定义儒略周期为7980年，是因28、19、15的最小公倍数为28×19×15=7980。其中：

28年为一太阳周期(solar cycle)，经过一太阳周期，则星期的日序与月的日序会重复。

19年为一太阴周期，或称默冬章(Metonic cycle)，因235朔望月=19回归年，经过一太阴周期则阴历月年的日序重复。

15年为一小纪(indiction cycle)，此为罗马皇帝君士坦丁（Constantine）所颁，每15年评定财产价值以供课税，成为古罗马用的一个纪元单位，

故以7980年为一儒略周期，而所选的起点公元前4713年，则是这三个循环周期同时开始的最近年份。

以儒略日计日是为方便计算年代相隔久远或不同历法的两事件所间隔的日数。

由于儒略日数字位数太多，国际天文学联合会于1973年采用简化儒略日（MJD），其定义为 MJD = JD– 2400000.5。MJD相应的起点是1858年11月17日世界时0时。