mktable2.m

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

function mktable2 (file,tsat,elev,cutoff,gdop,pdop,tdop,hdop,fdop,ldop,vdop)
%MKTABLE2: Create a table with results (DOP's);
%
% The function create a table showing numerical values for azimuth 
% and elevation of the satellites, as well as the positions of the
% satellites in earth-fixed WGS'84 coordinates.
%
% Syntax:
%    mktable2 (file,tsat,elev,cutoff,gdop,pdop,tdop,hdop,fdop,ldop,vdop)
%
% Input arguments:
%    file   - Filename of file to be created
%             empty string ==> write to screen
%    tsat   - Times on which positions are computed
%    elev   - Elevations of satellites at times tsat
%    cutoff - Selected cutoff elevation
%    gdop   - Geometrical Dilution of Precision
%    pdop   - Positional  Dilution of Precision
%    tdop   - Time Dilution of Precision
%    hdop   - Horizontal Dilution of Precision
%    fdop   - Lattitude Dilution of Precision
%    ldop   - Longitude Dilution of Precision
%    vdop   - Vertical Dilution of Precision
%
% Output arguments:
%    none

% ----------------------------------------------------------------------
% File.....: mktable2.m
% Date.....: 12-OCT-2000
% Version..: 2.0
% Author...: Peter Joosten
%            Mathematical Geodesy and Positioning
%            Delft University of Technology
% ----------------------------------------------------------------------

% -----------------
% --- Open file ---
% -----------------

if length(file) ~= 0;
  fid = fopen (file,'wt'); 
else;
  fid = 1;
end;

% -------------------
% --- Write DOP's ---
% -------------------

fprintf (fid,'|--------------------------|--------|-------|-------|-------|-------|-------|-------|-------|\n');
fprintf (fid,'|          Date/time [UTC] | # sats |  GDOP |  PDOP |  TDOP |  HDOP |  FDOP |  LDOP |  VDOP |\n');
fprintf (fid,'|--------------------------|--------|-------|-------|-------|-------|-------|-------|-------|\n');

for i = 1:length(tsat);
  
    str = gpst2str (tsat(i));
    fprintf (fid,'| %24s | %6d | %5.2f | %5.2f | %5.2f | %5.2f | %5.2f | %5.2f | %5.2f |\n',...
    str,length(find(elev(:,i)>cutoff)),gdop(i),pdop(i),tdop(i),hdop(i),fdop(i),ldop(i),vdop(i));
    
end;

fprintf (fid,'|--------------------------|--------|-------|-------|-------|-------|-------|-------|-------|\n\n\n');

% ------------------
% --- Close file ---
% ------------------

if fid ~= 1;
  fclose (fid);
end;

% ----------------------
% --- End of routine ---
% ----------------------