earth_shape.txt

An Example Simulation of Differential GPS

function y = earth_shape()

% Contains earth parameters [major-axis, minor-axis, ellipticity,
% eccentricity, sidereal rate]

Re = 6378137.0; % semi-major axis (equatorial) radius in meters

Rp = 6356752.3142; % semi-minor axis (polar) radius in meters

ellipticity = (Re-Rp)/Re; % ellipticity

eccentricity = sqrt(ellipticity*(2-ellipticity));

Omega = 7.292115167e-5;

y = [Re, Rp, ellipticity, eccentricity, Omega];