mjd2date.m

来自「Modified Julian date conversion.」· M 代码 · 共 25 行

function [day,month,year]=mjd2date(mjd)

a=fix(mjd)+2400001.0;

part_of_day=mjd-fix(mjd);

if a<2299161.0
   b=0.0;                            % for Julijan calendar
   c=a+1524.0;
else
   b=fix((a-1867216.25)/36524.25);   % for Gregorian calendar
   c=a+b-fix(b/4)+1525.0;
end

d=fix((c-121.1)/365.25);
e=fix(365.25*d);
f=fix((c-e)/30.6001);

day=c-e-fix(30.6001*f)+part_of_day;
month=f-1-12*fix(f/14);
year=d-4715.0-fix((7+month)/10);