pldopsmatrix.m

Sky lab calculate position of satellite in GNSS systems and calculate elevation and azimuth them and

function pldopsmatrix (tsat,dops,maxdop);
%PLDOPS: Plot DOP's (Dilution of Precision)
%
% The function creates a plot showing the different DOP'values, including
% GDOP, PDOP, TDOP, HDOP, FDOP, LDOP and VDOP as a function of time
%
% Syntax:
%    pldops (tsat,gdop,pdop,tdop,hdop,fdop,ldop,vdop,maxdop);
%
% Input arguments:
%    tsat   - Time (seconds)
%    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.....: pldops.m
% Date.....: 28-MAY-1999
% Version..: 1.0
% Author...: Peter Joosten
%            Mathematical Geodesy and Positioning
%            Delft University of Technology
% ----------------------------------------------------------------------



% --------------------------------------
% --- Create the figure, set options ---
% --------------------------------------

h = figure('Units','normalized', ...
	   'Position',[0.2 0.5 0.6 0.3], ...
	   'Color',[0 0 0], ...
	   'Name','DOP-values vs. time', ...
	   'NumberTitle','off');

h = axes ('Units','normalized', ...
	  'Position',[0.05 0.10 0.90 0.80], ...
	  'Color',[0 0 0], ...
	  'XColor',[1 1 1], ...
	  'YColor',[1 1 1], ...
	  'Xgrid','on', ...
	  'Ygrid','on', ...
	  'Box','on');

hold on;

% -----------------------------------------
% --- Limits, grid, labels and such ... ---
% -----------------------------------------

xlabel ('Time [sec]');
ylabel ('DOP-values');

ymin = 0.9*min(min(dops));
ymax = 1.1*max(max(dops));
if nargin >= 3; ymax = maxdop; end;

set (gca,'Xlim',[min(tsat) max(tsat)]);
set (gca,'Ylim',[ymin ymax]);

% -------------------
% --- Actual plot ---
% -------------------
gdop=dops(1,:);
h = plot (tsat,gdop,'y-');
set (h,'LineWidth',2);

idx = round((1/8)*length(tsat));
h   = text(tsat(idx),gdop(idx)+0.1,'GDOP');
set (h,'Color','yellow');
      
pdop=dops(2,:);
h = plot (tsat,pdop,'g-');
set (h,'LineWidth',2);

idx = round((2/8)*length(tsat));
h   = text(tsat(idx),pdop(idx)+0.1,'PDOP');
set (h,'Color','green');

tdop=dops(3,:);
h = plot (tsat,tdop,'c-');
set (h,'LineWidth',2);

idx = round((3/8)*length(tsat));
h   = text(tsat(idx),tdop(idx)+0.1,'TDOP');
set (h,'Color','cyan');

hdop=dops(4,:);
h = plot (tsat,hdop,'r-');
set (h,'LineWidth',2);

idx = round((4/8)*length(tsat));
h   = text(tsat(idx),hdop(idx)+0.1,'HDOP');
set (h,'Color','red');

fdop=dops(5,:);
h = plot (tsat,fdop,'b-');
set (h,'LineWidth',2);

idx = round((5/8)*length(tsat));
h   = text(tsat(idx),fdop(idx)+0.1,'FDOP');
set (h,'Color','blue');

ldop=dops(6,:);
h = plot (tsat,ldop,'w-');
set (h,'LineWidth',2);

idx = round((6/8)*length(tsat));
h   = text(tsat(idx),ldop(idx)+0.1,'LDOP');
set (h,'Color','white');

vdop=dops(7,:);
h = plot (tsat,vdop,'m-');
set (h,'LineWidth',2);
idx = round((7/8)*length(tsat));
h   = text(tsat(idx),vdop(idx)+0.1,'VDOP');
set (h,'Color','magenta');

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