julday.m

此功能包用于各种GPS坐标和时间的转换

function jd = julday(y,m,d,h)
% JULDAY  Conversion of date as given by
% Input :
%         y ... year (four digits)
%         m ... month
%         d ... day
%         h ... hour and fraction hereof (for example, 13h 40m 23.0s -> 13.6731)
%         The conversion is only valid in the time span
%         from March 1900 to February 2100

% Output: Julian day(jd)
% See  Hofmann-Wellenhof et al., p. 41--42

%Kai Borre 05-18-96
%Copyright (c) by Kai Borre
%$Revision: 1.0 $  $Date: 1997/09/26 $
% Modified by Chon 25 Spet 2000

if m <= 2
   y = y-1; 
   m = m+12;
end
jd = floor(365.25*y)+floor(30.6001*(m+1))+d+h/24+1720981.5;
%%%%%%% end julday.m  %%%%%%%%%%%%