example1.m

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

%EXAMPLE1: Demonstrate capabilities of SkyPlot (Matlab-version)
%
% This script demonstrates some of the capabilities of the Matlab version
% of SkyPlot. This demo reads an ephemeris-file, computes
% satellite-positions and converts them to azimuths/elevations for a
% choosen station (Delft in this case). The results are presented as tables
% and figures. 

% ----------------------------------------------------------------------
% File.....: example1.m
% Date.....: 03-JAN-2000
% Version..: 1.0
% Author...: Peter Joosten
%            Mathematical Geodesy and Positioning
%            Delft University of Technology
% ----------------------------------------------------------------------

% ----------------------------------------------
% --- Initialize, close all previous figures ---
% ----------------------------------------------

close all;
clear all;

% -------------------------------------------------------------
% --- Read ephemeris (RINEX nav-file or YUMA almanac)       ---
% --- and delete double occurences of SVN's                 ---
% -------------------------------------------------------------

eph = sortrows(rdyuma('yuma79.txt'));

% -----------------------------------------
% --- Station coordinates (WGS84 / XYZ) ---
% -----------------------------------------

station(1) = 3923625; % Delft (Geodesy-building)
station(2) =  298462;
station(3) = 5002803;
cutoff     = 10;      % Cutoff elevation in degrees

% -----------------
% --- Time/span ---
% -----------------

tsat = mktsat ('28-FEB-2001 12:00:00','28-FEB-2001 14:00:00',300);

% -------------------------------------------------------
% --- Compute satellite positions (XYZ WGS84 & AZ/EL) ---
% -------------------------------------------------------

[xsat,ysat,zsat,azim,elev,rotmat]    = cpaziele (tsat,eph,station);
[gdop,pdop,tdop,hdop,fdop,ldop,vdop] = cpdops (tsat,xsat,ysat,zsat,station,elev,cutoff,rotmat);
[mdbs, snrx, nabx, naby, nabz, nabt] = cprels (tsat,xsat,ysat,zsat,station,elev,cutoff);

% -----------------------------------
% --- Create plot(s) and table(s) ---
% -----------------------------------

%mktable1  ('../data/table01.txt',eph,tsat,azim,elev,xsat,ysat,zsat,cutoff);
%mktable2  ('../data/table02.txt',tsat,elev,cutoff,gdop,pdop,tdop,hdop,fdop,ldop,vdop);
%mktable3  ('../data/table03.txt',tsat,eph,elev,pdop,cutoff);
%mktable4  ('',eph,tsat,azim,elev,mdbs,snrx,nabx,naby,nabz,nabt,cutoff);

pleletim (eph,tsat,elev,cutoff);
plazitim (eph,tsat,azim,elev,cutoff);
plskyplt (eph,azim,elev,cutoff);
plnumsat (eph,tsat,elev,cutoff);
plvissat (eph,tsat,elev,cutoff);
pleleazi (eph,azim,elev,cutoff);
pldops   (tsat,gdop,pdop,tdop,hdop,fdop,ldop,vdop);
plworld  (eph,tsat,xsat,ysat,zsat);

% ------------------------------
% --- End of script EXAMPLE1 ---
% ------------------------------