matlab年积日转GPS周，年积日转换为年月日，年月日转化为儒略日

年月日转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