cpdesignpl.m

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

function [amatpl] = cpdesignpl (xsat,ysat,zsat,station,npl);
%CPDESIGN: Create the design-matrix for pseudolite-case
%
% The function creates the design-matrix, which is necessary for the
% computation of DOP-values and reliability measures. Note that this
% function computes the design matrix for a single epoch, therefore
% the input-arguments (xsat,ysat,zsat,station,elev) should be given
% for that single epoch as well.
%
% Syntax:
%    [amatpl] = cpdesignpl (xsat,ysat,zsat,station,npl);
%
% Input arguments:
%    xsat    - X-positions of satellites (WGS'84, Earth-fixed)
%    ysat    - Y-positions of satellites (WGS'84, Earth-fixed)
%    zsat    - Z-positions of satellites (WGS'84, Earth-fixed)
%    station - Station coordinates (WGS'84, Earth-fixed)
%    npl     - number of pseudolites
%
% Output arguments:
%    amatpl    - Design-matrix

% ----------------------------------------------------------------------
% File.....: cpdesign.m
% Date.....: 12-OCT-2000
% Version..: 2.0sv
% Author...: Peter Joosten
%            Mathematical Geodesy and Positioning
%            Delft University of Technology
% Edited by: Sandra Verhagen
%            11-DEC-2000
% ----------------------------------------------------------------------

xsat = xsat - station(1);
ysat = ysat - station(2);
zsat = zsat - station(3);
dist = sqrt(xsat.*xsat + ysat.*ysat + zsat.*zsat);

amatpl = zeros (npl,4);
  
for j = 1:npl;
    
  amatpl(j,1) = xsat(j)/dist(j);
  amatpl(j,2) = ysat(j)/dist(j);
  amatpl(j,3) = zsat(j)/dist(j);
  amatpl(j,4) = -1;
 
end;