init_satdata.m

来自「gps“可用性”matlab代码 Matlab Algorithm Avail」· M 代码 · 共 40 行

function satdata = init_satdata(geodata, alm_param, satdata, t)

%*************************************************************************
%*     Copyright c 2001 The board of trustees of the Leland Stanford     *
%*                      Junior University. All rights reserved.          *
%*     This script file may be distributed and used freely, provided     *
%*     this copyright notice is always kept with it.                     *
%*                                                                       *
%*     Questions and comments should be directed to Todd Walter at:      *
%*     [e-mail]                                              *
%*************************************************************************
%
global COL_SAT_PRN COL_SAT_XYZ COL_SAT_XYZDOT COL_SAT_UDREI ...
        COL_SAT_COV COL_SAT_SCALEF COL_SAT_MAX
global CONST_H_GEO

% get sat positions
if isempty(satdata), % get from almanac
	ngps=size(alm_param,1);
    if(~isempty(geodata))
      ngeo=size(geodata,1);
      geo = ngps + (1:ngeo);
      satdata(geo,:)=repmat(NaN,ngeo,COL_SAT_MAX);
      satdata(geo,COL_SAT_PRN) = geodata(:,1);
      satdata(geo,COL_SAT_XYZ) = llh2xyz([zeros(ngeo,1) geodata(:,2) ...
	                                        CONST_H_GEO*ones(ngeo,1)]);
      satdata(geo,COL_SAT_XYZDOT) = zeros(ngeo,3);
  end
end
[prn,satxyz,satvel] = alm2satposvel(t,alm_param);

ngps = size(prn,1);
gps=1:ngps;
satdata(gps,:) = repmat(NaN,ngps,COL_SAT_MAX);
satdata(gps,COL_SAT_PRN) = prn;
satdata(gps,COL_SAT_XYZ) = satxyz;
satdata(gps,COL_SAT_XYZDOT) = satvel;