MJD与日期的相互转化（matlab）

1.日期转换为mjd

先将日期转换成jd，jd再减去2400000.5就变成了mjd。在转换成jd的时候，要注意在1582.10.5之前使用的是儒略历，自此之后使用的是格里高利历，这两者在转换成jd的时候计算公式有些不同。

clear all;
clc;
y=input('请输入年份：');
m=input('请输入月份：');
d=input('请输入日：');
while m<=2
    y=y-1;
    m=m+12;
end
a=floor(y/100);
if (y<1582)||((y==1582)&&(m<10))||((y==1582)&&(m==10)&&(d<5))
    b=0;%儒略历时，b=0
else 
    b=2-a+floor(a/4);%格里高利历时b的取值
end
jd=floor(365.25*(y+4716))+floor(30.6001*(m+1))+d+b-1524.5;
mjd=jd-2400000.5;
fprintf('mjd:%f',mjd)

 

2.mjd转换成日期

首先将mjd转换成jd再转换成