vp90ecef.m

GPS software toolbox for GPS receiver development

%                              vp90ecef.m
%  Scope:   This MATLAB macro performs the transformation from PZ-90 ("Parametri
%           Zemli 1990", i.e. earth's parameters 1990) frame to ECEF WGS-84 
%           frame coordinates for a given position vector with components 
%           expressed in meters.
%  Usage:   pecef = vp90ecef(pecef)
%  Description of parameters:
%           ppz90  -  input, PZ-90 frame position vector, with components in 
%                     meters
%           pecef  -  output, ECEF (WGS-84) position vector with components in 
%                     meters
%  Reference: 
%           [1] Misra, P. N., Abbott, R. I., Gaposchkin, E. M., Integrated use
%               of GPS and GLONASS: Transformation between WGS84 and  PZ-90. 
%               Proceedings of the 9th International Technical Meeting of the 
%               Satellite Division of the Institute of Navigation, ION GPS-96, 
%               Kansas City, September 16-20, 1996, pp. 307-314. 
%  Last update:  04/04/00
%  Copyright (C) 1997-00 by LL Consulting. All Rights Reserved.

function   pecef = vp90ecef(ppz90)

rotangle = 1.9e-6;       %  rotation angle in radians
d = 2.5;                 %  translation distance in meters

pecef = zeros(3,1);
pecef(1) = ppz90(1) - ppz90(2) * rotangle;
pecef(2) = ppz90(1) * rotangle + ppz90(2) + 2.5;
pecef(3) = ppz90(3);