Code covered by the BSD License
-
CalcEA(M,ecc,tol)
Orbit eccentric anomaly, Kepler's equation keplers equation
-
Groundtrack(Kepler,GMSTo,Tf,f...
Orbit groundtrack plot Latitude longitude lat long
-
Hohmann(R_init,R_fin,U)
Orbit Hohmann transfer
-
JD(yr,day)
Julian Date
-
KeplerCOE(Ro,Vo,dT,U,tol)
Orbit Kepler position velocity
-
NodeChange(dO,inc,Vinit)
Node change right ascension of the ascending node RAAN raan orbit
-
R1(x)
Rotation matrix direction cosine matrix
-
R2(x)
Rotation matrix direction cosine matrix
-
R3(x)
Rotation matrix direction cosine matrix
-
RVtoLatLong(ECEF)
orbit radius velocity latitude longitude ECEF
-
TwoBody(t,X,U)
Two body Orbit gravity
-
dInc(V,dI,fpa)
Inclination change orbit gravity
-
dVdI(R_init,R_fin,Inc,U,Tol)
Inclination change velocity change orbit hohmann transfer
-
ecef2eci(ECEF, GST, V_ECEF)
Orbit ECEF ECI Coordinate conversion
-
eci2ecef(ECI, GST, V_ECI)
Orbit ECEF ECI Coordinate conversion
-
elorb(R,V,U,tol)
Kepler orbital elements ECI Position orbit conversion
-
nuFromM(M,ecc,tol)
Kepler Orbit Anomaly true mean
-
nuFromTp(Tp,ecc,n,tol)
Kepler Orbit Anomaly true time periapse perigee
-
plotorb(ECEF, V_ECEF, mu, Rbo...
Orbit gravity plot orbit spherical
-
randv(a,ecc,inc,Omega,w,nu,U)
Kepler orbital elements ECI Position orbit conversion
-
topo(ECEF, lat, long, h, Rp)
Orbit range elevation azimuth position ground station site latitude longitude
-
zeroTo360(x,unit)
Angle reduce reduction degrees radians
-
Constants.m
-
View all files
from
Orbital Mechanics Library
by Richard Rieber
A compilation of all of the functions I wrote for my orbital mechanics class
|
| zeroTo360(x,unit)
|
%Angle reduce reduction degrees radians
% Richard Rieber
% October 1, 2009
% rrieber@gmail.com
%
% function y = zeroTo360(x,unit)
%
% Purpose: This function reduces an angle to the range of 0 - 360 degrees
% or 0 - 2*pi radians.
%
% Inputs: x - Angle to be reduced, may be an array of angles
% unit - Boolean, 1 for radians, 0 for degrees, defaults to
% degrees [OPTIONAL]
%
% Output: y - Reduced angle
%
function y = zeroTo360(x,unit)
if nargin == 1
unit = 0;
elseif nargin > 2
error('Too many inputs')
end
if unit
deg = 2*pi;
else
deg = 360;
end
y = zeros(1,length(x));
for j = 1:length(x)
if (x(j) >= deg)
x(j) = x(j) - fix(x(j)/deg)*deg;
elseif (x(j) < 0)
x(j) = x(j) - (fix(x(j)/deg) - 1)*deg;
end
y(j) = x(j);
end
|
|
Contact us at files@mathworks.com