儒略日转公历算法

1. % Gregorian Date from Julian Date 
2. % This calculation is valid for any Julian Day Number including negative JDN and 
3. % produces a Gregorian date (or possibly a proleptic Gregorian date)
5. Z = floor(JD0 - 1721118.5); % JD0 is the Julian Date in propagation 
6. R = (JD0 - 1721118.5 - Z);  % R is the fractional part of JD0
7. G = (Z - .25);  
8. A = (floor(G / 36524.25));  % Calculate the value of A which is the number of full centuries
9. B = (A - (A / 4)); % The value of B is this number of days minus a constant
10. year = (floor((B+G) / 365.25)); % Calculate the value of Y, the year in a calendar whose years start on March 1
11. C = (B + Z - floor(365.25 * year));  % Day count
12.          month = (fix((5 * C + 456) / 153));  % Month
13.          UT = (C - fix((153 * month - 457) / 5) + R); % Calculation for UTC
14.          day = floor(UT);  % Gregorian Day