年月日转GPS周,年积日转换为年月日,年月日转化为儒略日的函数(注意:函数一里面也有年月日转化为儒略日,如果想运行,注意函数套娃)
% 函数一:年积日转换为GPS周
function gps_week = y_m_d2gpsweek(year,doy)
% 年月日转化为儒略日
[year,mon,day]=doy2y_m_d(year, doy);
% 年月日转
mjd = y_m_d2julian_d(year, mon, day);
week = floor((mjd - 44244.0) / 7.0);
day = floor(mjd - 44244.0 - week * 7.0);
gps_week = [num2str(week),num2str(day)];
end
% 函数二:年积日转换为年月日
function [year,mon,day]=doy2y_m_d(year, doy)
day = doy;
mon = 0;
monthdays = [0, 31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31];
for i=1:13
monthday = monthdays(i);
if i==3 && leapyear(year) == 1
monthday = monthday + 1;
end
if day > monthday
day = day - monthday;
else
mon = i - 1;
break
end
end
end
% 函数三:年月日转化为儒略日
function mjd=y_m_d2julian_d(year, mon, day)
if mon <= 2
mon = mon + 12;
year = year - 1;
end
mjd = 365.25 * year - rem(365.25 * year , 1.0) - 679006.0;
mjd = mjd + floor(30.6001 * (mon + 1)) + 2.0 - ...
floor(year / 100.0) + floor(year / 400) + day;
end