%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 buildingstation2 = [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 ---% ------------------------------