ut_jd.m

在gps定位中需要将时间转换为儒略日

function [jd,mjd]=ut_jd(Y,M,D)
%将民用年历转化为儒略日的函数
%time必须是一天之内以小时为单位的时间
if M<=2
    y=Y-1;m=M+12;
else
    y=Y;m=M;
end
jd=fix(365.25*y)+fix(30.6001*(m+1))+D+1720981.5;
mjd=jd-2400000.5;