utc2gpstnum.m

来自「Sky lab calculate position of satellite 」· M 代码 · 共 63 行

function [gpsweek,secs] = utc2gpstnum (T);
%UTC2GPST: Convert time/date string to GPSWEEK / Seconds
%
% This routine converts a time and date, given as a string in the form of
% 'dd-mon-yyyy hh:mm:ss.sss' into the gpsweek and seconds into this
% week. If the global variable GPSWEEK (uppercase) is defined, the seconds
% will be relative to this week.
%
% Syntax:
%    [gpsweek,secs] = utc2gpst (T);
%
% Input arguments:
%    T      : vector like [year month day hours minutes seconds] (all numerical)
%
% Output arguments:
%    gpsweek: gpsweek number
%    secs   : Seconds relative to this week
%

% ----------------------------------------------------------------------
% File.....: utc2gpst.m
% Date.....: 27-MAY-1999
% Version..: 1.0
% Author...: Peter Joosten
%            Mathematical Geodesy and Positioning
%            Delft University of Technology
% ----------------------------------------------------------------------

% ---------------------------------------
% --- Declare global variable GPSWEEK ---
% ---------------------------------------

global GPSWEEK;

% ----------------------
% --- Convert vector ---
% ----------------------


gpst    = mjd2gps(ymd2mjd([T(1) T(2) T(3)]));
gpst(2) = gpst(2) + 3600*T(4) + 60*T(5) + T(6);

gpsweek = gpst(1);
secs    = gpst(2);

% -------------------------------------------------------------
% --- if GPSWEEK was defined, adapt the results accordingly ---
% -------------------------------------------------------------

if length(GPSWEEK) ~= 0;

  secs    = secs + (gpsweek - GPSWEEK) * 7*24*3600;
  gpsweek = GPSWEEK;

end ;

% --------------------------------
% --- End of function utc2gpst ---
% --------------------------------