example2.m

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

%EXAMPLE2: Demonstrate capabilities of SkyLab (Skyplot for Matlab)
%
% This script demonstrates some of the capabilities of Skylab, the Matlab
% version of SkyPlot. This demo reads an ephemeris-file, computes
% satellite-positions and converts them to azimuths/elevations for two choosen
% stations (Delft/Greenwich in this case). Finally the satellites in common
% view of both stations are determined

% ----------------------------------------------------------------------
% File.....: example2.m
% Date.....: 31-MAY-1999
% 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('../data/yuma59.txt'));

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

station1 = [3924697.206 301131.833  5001903.289]; % Delft, Geodesy building
station2 = [4275717.804      0.000  4716876.330]; % Greenwich

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

tsat = mktsat ('11-OCT-2000 12:00:00','11-OCT-2000 14:00:00',300);

% ------------------------
% --- Cutoff elevation ---
% ------------------------

cutoff = 10;

% -----------------------------------------------------------------
% --- Find satellites in common view and write results to table ---
% -----------------------------------------------------------------

intervis ('',tsat,station1,station2,eph,cutoff);

% ------------------------------
% --- End of script EXAMPLE2 ---
% ------------------------------