cpelev1.m

In the planning stage of a GNSS-measurement campaign it is useful to compute several design para-met

function [elev] = cpelev1 (xsat,ysat,zsat,plhstation,xyzstation);
%CPAZIELE: Compute elevation/azimuth of satellites/pseudolites
%
% This routine computes satellite positions in an earth-fixed WGS'84
% frame. If station coordinates are specified (in WGS'84), elevations and
% azimuths are computed as well.
%
% Syntax:
%    [elev] = cpaziele1 (xsat,ysat,zsat,plhstation,xyzstation)
%
% Input arguments:
%    xsat       - xsat(i) = X-coordinate (WGS'84) for satellite (i)
%    ysat       - ysat(i) = Y-coordinate (WGS'84) for satellite (i)
%    zsat       - zsat(i) = Z-coordinate (WGS'84) for satellite (i)
%    plhstation - WGS'84 coordinates for station (phi,lambda,height)
%    xyzstation - WGS'84 coordinates for station (x,y,z)
%
% Output arguments:
%    elev       - elev(i,j) = elevation for satellite (i)
% 

% ----------------------------------------------------------------------
% File.....: cpelev.m
% Date.....: 25-MAY-1999
% Version..: 1.0sv
% Author...: Peter Joosten
%            Mathematical Geodesy and Positioning
%            Delft University of Technology
% Edited by: Sandra Verhagen
%            28-jun-2000
% ----------------------------------------------------------------------

% ---------------------------
% --- Initialize matrices ---
% ---------------------------

elev = zeros(size(xsat,1),1);

% -------------------------------------------------
% --- Compute azimuths/elevations, if requested ---
% -------------------------------------------------
xsat = xsat - xyzstation(1);
ysat = ysat - xyzstation(2);
zsat = zsat - xyzstation(3);
 
rotmat(1,1) = - sin(plhstation(1)) * cos(plhstation(2));
rotmat(1,2) = - sin(plhstation(1)) * sin(plhstation(2));
rotmat(1,3) =   cos(plhstation(1))                     ;
rotmat(2,1) = -                      sin(plhstation(2));
rotmat(2,2) =                        cos(plhstation(2));
rotmat(2,3) =   0d0;
rotmat(3,1) = - cos(plhstation(1)) * cos(plhstation(2));
rotmat(3,2) = - cos(plhstation(1)) * sin(plhstation(2));
rotmat(3,3) = - sin(plhstation(1))                     ;

%for i = 1:size(xsat,1);
  
  neh  = rotmat * [xsat ysat zsat]';
  dist = sqrt(neh(1,:).^2 + neh(2,:).^2 + neh(3,:).^2);
  elev = (asin(-neh(3,:)./dist))*180/pi;
  
%end;

% --------------------------------
% --- End of function cpelev1  ---
% --------------------------------