Code covered by the BSD License
- Delta_R_Trop=Error_Tropos...This Function approximate Troposspheric Group Delay Base on
- ECEF2GPS(Pos)
ECEF2LLA - convert earth-centered earth-fixed (ECEF)
- XYZ2ENU(A,Phi,Lambda)
- [Delta_I]=Error_Ionospher...This Function approximate Ionospheric Group Delay
- [E,A]=Calc_Azimuth_Elevat...This Function Compute Azimuth and Elevation of satellite from reciever
- [G,Delta_X,Pos_Rcv_n,B]=G...This Function use Ephemeris Data and Calculate satellite Position
- [Pos_xyz_Mat,Orbit_parame...This Function use Ephemeris Data and Calculate satellite Position
- d=Distance(Pos_SV,Pos_Rcv);
- dTclk_Offset=Error_Satell...This Function calculate Satellite Offest Clock Error Base on GPS Theory *
- dTclk_Rel=Error_Satellite...This Function calculate Satellite Clock Error Base on
- f=Kepler_Eq(x,M,e)
- plot_Orbit(Orbit_paramete...This Function show the satellite and its orbit+Earth geometry
- main.m
- main_request.m
-
View all files
from
GPS Navigation Toolbox
by moein mehrtash
This program include
Principle of Radio Navigation Calculation and GPS position error corrections.
|
| Delta_R_Trop=Error_Tropospheric_Hopfield(T_amb,P_amb,P_vap,Pos_Rcv,Pos_SV)
|
%This Function approximate Troposspheric Group Delay Base on
%application . edited by B. Parkinson,J. Spilker, P.Enge, AIAA,1996
%CopyRight By Moein Mehrtash
%**************************************************************************
% Written by Moein Mehrtash, Concordia University, 3/21/2008 *
% Email: moeinmehrtash@yahoo.com *
%**************************************************************************
% Reference:"GPS Theory and application",edited by B.Parkinson,J.Spilker, *
%**************************************************************************
%Input
% T_amb:'C =>At reciever antenna location
% P_amb:hPa =>At reciever antenna location
% P_vap:hPa =>Water vapore pressure at reciever antenna location
% Pos_Rcv : XYZ position of reciever (Meter)
% Pos_SV : XYZ matrix position of GPS satellites (Meter)
%Output:
% Delta_R_Trop: m =>Tropospheric Error Correction
%**************************************************************************
function Delta_R_Trop=Error_Tropospheric_Hopfield(T_amb,P_amb,P_vap,Pos_Rcv,Pos_SV)
S=size(Pos_SV);
m=S(1);n=S(2);
for i=1:m
[E,A0]=Calc_Azimuth_Elevation(Pos_Rcv,Pos_SV(i,:));
El(i)=E; %Elevation Rad
A(i)=A0; %Azimoth Rad
end
%Zenith Hydrostatic Delay
Kd=1.55208*10^(-4)*P_amb*(40136+148.72*T_amb)/(T_amb+273.16);
%Zenith Wet Delay
Kw=-.282*P_vap/(T_amb+273.16)+8307.2*P_vap/(T_amb+273.16)^2;
for i=1:m
Denom1(i)=sin(sqrt(El(i)^2+1.904*10^-3));
Denom2(i)=sin(sqrt(El(i)^2+.6854*10^-3));
%Troposhpheric Delay Correctoion
Delta_R_Trop(i)=Kd/Denom1(i)+Kw/Denom2(i); % Meter
end
|
|
Contact us