ut_jd.m

来自「在gps定位中需要将时间转换为儒略日」· M 代码 · 共 10 行

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;