ephtest.m

这是国外关于卫星导航方面一书的源代码

function [x,y,z,tsat,rela]=ephtest(station,wie)
% -----------------------------------------------------------------------------
%  GPSLab - EPHTEST.M  orbit calculation and so on for all satellites represented 
%                      in the Rinex-NAV-File, referenced to  the WGS84 (ECEF or pseudo-ECI).
% -----------------------------------------------------------------------------
%
%  INPUT:   station    option for ephemeris file of either Ref.- or Roverstation
%           wie        option for calculation due to ECEF (=1) [default]
%                                    or due to pseudo-ECI (=0) system 
%
%  OUTPUT:  x,y,z,tsat,rela    look at MatLab-File eph2xyzn.m
%
%           about produced variables you get informed via file var_orb.txt 
%
% (c) iapg 1998 zeb, addenda: 19.2.1999 zeb
% -----------------------------------------------------------------------------

if nargin==1
   wie=1;
end

if wie~=1&wie~=0
   errordlg('Don''t do that! There are only the option values 0 or 1 for the 2nd variable (set to 1).', ...
      'GPSLab: Error during call of EPHTEST.M');
   wie=1;
end

if station == 1
   load s1eph;
elseif station == 2
	load s2eph;
	s1eph=s2eph;
   clear s2eph;
else
   warndlg(['unvalid variable "station": the value of "station" has to be ' ...
         '1 or 2 !'], 'GPSLab: Error during call of EPHTEST.M');
   return;
end

[lephs1,dummy]=size(s1eph);
if rem(lephs1,24)~=0
   errordlg(['Error in program EPHTEST.M during the reading of the ephemeris.',...
             ' The length of the files S1EPH,S2EPH does not coincide with the multiple of 24' ...
             ' and so seems defect. Check for example the RINEX-Import.'], ...   
   'GPSLab: Break - (c) iapg 1998 zeb');
   return;
end

% s1eph=s1eph(1:(6*24));   % Test mit nur 6 Satelliten

h=waitbar(0,'calculation of the 24h-arcs of the satellites');
for j=1:((length(s1eph))/24),
   for i=1:(24*6),
   [x(i,j),y(i,j),z(i,j),tsat(i,j),rela(i,j)]=eph2xyzn(-600+i*600,0,s1eph((j*24-23):(j*24)),wie);
	end
waitbar(j/((length(s1eph))/24));   
end
close(h);

%---------------------------- End of the universe --------------------------

% -----------------------------------------------------------------------------
% (c) iapg 1998 zeb
% -----------------------------------------------------------------------------