About Site Map Submit Contact Us Log in | Create an account
Create an account Log In
Average Rating
User Rating:
Visitors Rating:
My rating:

Write review
See full specifications

scripts default iconOrbital Elements from Position/Velocity Vectors (Scripts) Publisher's description

vec2orbElem(rs,vs,mus) converts positions (rs) and velocities (vs)

vec2orbElem(rs,vs,mus) converts positions (rs) and velocities (vs)
of bodies with gravitational parameters (mus) to Keplerian orbital elements.

rs: 3n x 1 stacked initial position vectors:
or 3 x n matrix of position vectors.
vs: 3n x 1 stacked initial velocity vectors or 3 x n matrix
mus: gravitational parameters (G*m_i) where G is the
gravitational constant and m_i is the mass of the ith body.
if all vectors represent the same body, mus may be a
a: semi-major axes
e: eccentricities
E: eccentric anomalies
I: inclinations
omega: arguments of pericenter
Omega: longitudes of ascending nodes
P: orbital periods
tau: time of periapsis crossing
A, B: orientation matrices (see Vinti, 1998)

All units must be complementary, i.e., if positions are in AUs, and
time is in days, dx0 must be in AU/day, mus must be in
AU^3/day^2 (these are the units in solarSystemData.mat).

The data in solarSystemData.mat was downloaded from JPL's System Web
Interface (http://ssd.jpl.nasa.gov/?horizons). It includes
positions for the planets, the sun and pluto (because I went to
grade school before 2006). Positions for planets with moons are for
the barycenters.

%solar system oribtal elements
ssdat = load('solarSystemData.mat');
rs = ssdat.p0(1:end-3) - repmat(ssdat.p0(end-2:end),9,1);
vs = ssdat.v0(1:end-3) - repmat(ssdat.v0(end-2:end),9,1);
mus = ssdat.mus(1:9) + ssdat.mus(10);
[a,e,E,I,omega,Omega,P,tau,A,B] = vec2orbElem(rs,vs,mus);
%convert back:
r = A*diag(cos(E) - e) + B*diag(sin(E));
rdot = (-A*diag(sin(E))+B*diag(cos(E)))*...
diag(sqrt(mus(:).'.*a.^-3)./(1 - e.*cos(E)))

System Requirements:

MATLAB 7.9 (2009b)
Program Release Status: New Release
Program Install Support: Install and Uninstall

Orbital Elements from Position/Velocity Vectors (Scripts) Tags:

Click on a tag to find related softwares

Is Orbital Elements from Position/Velocity Vectors (Scripts) your software?

Manage your software

Most Popular

scripts default icon ASK, OOK, FSK, BPSK, QPSK, 8PSK modulation 1.1
ASK, OOK, FSK, BPSK, QPSK, 8PSK modulation contain several functions for digital modulation simulation
scripts default icon Simulink Communication Labs 1.1
Simulink Communication Labs allows you to learn communication systems in greater depth.
scripts default icon M-QAM modulation and demodulation 1.1
M-QAM modulation and demodulation is the QAM modulation and demodulation tech.
scripts default icon LZW Compression/Decompression 1.1
LZW Compression/Decompression - Updated LZW compressor and decompressor with reasonable performance
scripts default icon InSPIRE utility to plot a 2D displacement field (Scripts) 1.0
This program plots the deformation field (displace vectors) contained in vector.txt.