Movatterモバイル変換


[0]ホーム

URL:


Jump to content
WikipediaThe Free Encyclopedia
Search

Orbital elements

From Wikipedia, the free encyclopedia
Parameters that define a specific orbit
This articleneeds additional citations forverification. Please helpimprove this article byadding citations to reliable sources. Unsourced material may be challenged and removed.
Find sources: "Orbital elements" – news ·newspapers ·books ·scholar ·JSTOR
(April 2025) (Learn how and when to remove this message)

Orbital elements are theparameters required to uniquely identify a specificorbit. Incelestial mechanics these elements are considered intwo-body systems using aKepler orbit. There are many different ways to mathematically describe the same orbit, but certain schemes are commonly used inastronomy andorbital mechanics.

A real orbit and its elements change over time due to gravitationalperturbations by other objects and the effects ofgeneral relativity. A Kepler orbit is an idealized, mathematical approximation of the orbit at a particular time.

When viewed from aninertial frame, two orbiting bodies trace out distinct trajectories. Each of these trajectories has itsfocus at the commoncenter of mass. When viewed from a non-inertial frame centered on one of the bodies, only the trajectory of the opposite body is apparent; Keplerian elements describe these non-inertial trajectories. An orbit has two sets of Keplerian elements depending on which body is used as the point of reference. The reference body (usually the most massive) is called theprimary, the other body is called thesecondary. The primary does not necessarily possess more mass than the secondary, and even when the bodies are of equal mass, the orbital elements depend on the choice of the primary.

Orbital elements can be obtained fromorbital state vectors (position and velocity vectors along with time and magnitude of acceleration) by manual transformations or with computer software through a process known asorbit determination.[1]

It non-closed orbits exist, although these are typically referred to as trajectories and not orbits, as they are not periodic. The same elements used to describe closed orbits can also typically be used to represent open trajectories.

Common orbital elements by type

[edit]

Required parameters

[edit]

In general, eight parameters are necessary to unambiguously define an arbitrary and unperturbed orbit. This is because the problem contains eightdegrees of freedom. These correspond to the three spatialdimensions which define position (x,y,z in aCartesian coordinate system), the velocity in each of these dimensions, the magnitude ofacceleration (only magnitude is needed, as the direction is always opposite the position vector), and the current time (epoch). The mass orstandard gravitational parameter of the central body can specified instead of the acceleration, as one can be used to find the other given the position vector through the relationa=μ/r2{\displaystyle a=\mu /r^{2}}. These parameters can be described asorbital state vectors, but this is often an inconvenient and opaque way to represent an orbit, which is why orbital elements are commonly used instead.

When describing an orbit with orbital elements, typically two are needed to describe the size and shape of the trajectory, three are needed describe the rotation of the orbit, one is needed to describe the speed of motion, and two elements are needed to describe the position of the body around its orbit along with the epoch time at which this occurs. However, if the epoch time is chosen to be the time at which the position-describing element of choice (e.g. mean anomaly) is equal to some constant (usually zero), then said element can be omitted, meaning only seven elements are required in total.

Commonly only 6 variables are specified for a given orbit, as the motion-describing variable can be the mass or standard gravitational parameter of the central body, which is often already known and does not need specifying, and the epoch time can be considered part of the reference frame and not as a distinct element. However, in any case, 8 values will need to be known, regardless of how they are categorized.

Additionally, certain elements can be omitted if they are not required for the desired application (e.g. both epoch elements and the motion element are not needed if only the shape and orientation need to be known).

Size and shape describing parameters

[edit]

Two parameters are required to describe the size and the shape of an orbit. Generally any two of these values can be used to calculate any other (as described below), so the choice of which to use is one of preference and the particular use case.

  • Eccentricity (e) — shape of the ellipse, describing how much it deviates from a perfect a circle. An eccentricity of zero describes a perfect circle, values less than 1 describe an ellipse, values greater than 1 describe a hyperbolic trajectory, and a value of exactly 1 describes a parabola.[2]
  • Semi-major axis (a) — half the distance between theapoapsis and periapsis (long axis of the ellipse). This value is positive for elliptical orbits, infinity for parabolic trajectories, and negative for hyperbolic trajectories, which can hinder its usability when working with different types of trajectories.[3]
  • Semi-minor axis (b) — half the short axis of the ellipse. This value shares the same limitations as with the semi-major axis.
  • Semi-parameter (p) — the width of the orbit at the primary focus (at atrue anomaly ofπ/2 or 90°). This value is useful for its use in theorbit equation, which can return the distance from the central body givenp and the true anomaly for any type of orbit or trajectory. This value is also commonly referred to as the semi-latus rectum, and given the symbol. Additionally, this value will always be defined and positive unlike the semi-major and semi-minor axes.[3]
  • Apoapsis (ra) — The farthest point in the orbit from the central body (at a true anomaly ofπ or 180°). This quantity is undefined (or infinity) for parabolic and hyperbolic trajectories, as they continue moving away from the central body forever. This value is sometimes given the symbolQ.[2]
  • Periapsis (rp) — The closest point in the orbit from the central body (at a true anomaly of 0). Unlike with apoapsis, this quantity is defined for all orbit types. This value is sometimes given the symbolq.[2]

For perfectly circular orbits, there are no points on the orbit that can be described as either the apoapsis or periapsis, as they all have the same distance from the central body. Additionally it is common to see the affix for apoapsis and periapsis changed depending on the central body (e.g. apogee and perigee for orbits of theEarth, and aphelion and perihelion for orbits of theSun).

Other parameters can also be used to describe the size and shape of an orbit such as thelinear eccentricity,flattening, andfocal parameter, but the use of these is limited.

Relations between elements

[edit]
Further information:Conic section,Apsis, andSemi-major and semi-minor axes

This section contains the common relations between these orbital elements, but more relations can be derived through manipulations of one or more of these equations. The variable names used here are consistent with the ones described above.

Eccentricity can be found using the semi-minor and semi-major axes like so:
     e=1b2a2{\displaystyle e={\sqrt {1-{\frac {b^{2}}{a^{2}}}}}} whena>0{\displaystyle a>0},e=1+b2a2{\displaystyle e={\sqrt {1+{\frac {b^{2}}{a^{2}}}}}} whena<0{\displaystyle a<0}

Eccentricity can also be found using the apoapsis and periapsis through this relation:
     e=rarpra+rp{\displaystyle e={\frac {r_{a}-r_{p}}{r_{a}+r_{p}}}}

The semi-major axis can be found using the fact that the line that connects the apoapsis to the center of the conic, and from the center to the periapsis both combined span the length of the conic, and thus the major axis. This is then divided by 2 to get the semi-major axis.
     a=rp+ra2{\displaystyle a={\frac {r_{p}+r_{a}}{2}}}

The semi-minor axis can be found using the semi-major axis and eccentricity through the following relations. Two formula are needed to avoid taking thesquare root of a negative number.
     b=a1e2{\displaystyle b=a{\sqrt {1-e^{2}}}} whene<1{\displaystyle e<1},b=ae21{\displaystyle b=a{\sqrt {e^{2}-1}}} whene>1{\displaystyle e>1}

The semi-parameter can be found using the semi-major axis and eccentricity like so:
     p=a(1e2){\displaystyle p=a\left(1-e^{2}\right)}

Apoapsis can be found using the following equation, which is a form of theorbit equation solved forν=π{\displaystyle \nu =\pi }.
     ra=p1e{\displaystyle r_{a}={\frac {p}{1-e}}} , whene<1{\displaystyle e<1}

Periapsis can be found using the following equation, which, as with the equation for apoapsis, is a form of theorbit equation instead solved forν=0{\displaystyle \nu =0}.
     rp=p1+e{\displaystyle r_{p}={\frac {p}{1+e}}}

Rotation describing elements

[edit]
In this diagram, theorbital plane (yellow) intersects a reference plane (gray). For Earth-orbiting satellites, the reference plane is usually the Earth's equatorial plane, and for satellites in solar orbits it is theecliptic plane. The intersection is called theline of nodes, as it connects the reference body (the primary) with the ascending and descending nodes. The reference body and thevernal point (♈︎) establish a reference direction and, together with the reference plane, they establish a reference frame.

Three parameters are required to describe the orientation of the plane of the orbit, and the orientation of the orbit within that plane.

  • Inclination (i) — vertical tilt of the ellipse with respect to the reference plane, typically theequator of the central body, measured at theascending node (where the orbit passes crosses the reference plane, represented by the green anglei in the diagram). Inclinations near zero indicateequatorial orbits, and inclinations near 90° indicatepolar orbits. Inclinations from 90 to 180° are typically used to denoteretrograde orbits.
  • Longitude of the ascending node (Ω) — describes the angle from theascending node of the orbit ( in the diagram) to the reference frame's reference direction (♈︎ in the diagram). This is measured in the reference plane, and is shown as the green angleΩ in the diagram. This quantity is undefined for perfectly equatorial (coplanar) orbits, but is often set to zero instead by convention.[3] This quantity is also sometimes referred to as the right ascension of the ascending node (or RAAN).
  • Argument of periapsis (ω) — defines the orientation of the ellipse in the orbital plane, as an angle measured from the ascending node to the periapsis (the closest point the satellite body comes to the primary body around which it orbits), the purple angleω in the diagram. This quantity is undefined for circular orbits, but is often set to zero instead by convention.[3]

These three elements together can be described asEuler angles defining the orientation of the orbit relative to the reference coordinate system. Although these three are the most common, other elements do exist, and are useful to describe other properties of the orbit.

  • Longitude of periapsis (ϖ) — describes the angle between the vernal point and the periapsis, measured in the reference plane. This can be described as the sum of the longitude of the ascending node and the argument of periapsis:ϖ=Ω+ω{\displaystyle \varpi =\Omega +\omega }. Unlike the longitude of the ascending node, this value is defined for orbits where the inclination is zero.

Motion over time describing elements

[edit]

One parameter is required to describe the speed of motion of the orbiting object around the central body. However, this can be omitted if only a description of the shape of the orbit is required. Various quantities that do not directly describe a speed can be used to satisfy this condition, and it is possible to convert from one to any other (formula below).

  • Mean motion (n) — quantity that describes the averageangular speed of the orbiting body, measured as an angle per unit time. For non-circular orbits, the actual angular speed is not constant, so the mean motion will not describe a physical angle. Instead this corresponds to a change in themean anomaly, which indeed increases linearly with time.
  • Orbital period (P) — the time it takes for the orbiting body to complete one full revolution around the central body. This quantity is undefined for parabolic and hyperbolic trajectories, as they are non-periodic.
  • Standard gravitational parameter (μ) — quantity equal to the mass of the central body times thegravitational constantG. This quantity is often used instead of mass, as it can be easier to measure with precision than either mass orG, and will need to be calculated in any case in order to find the acceleration due to gravity. This is also often not included as part of orbital element lists, as it can assumed to be known based on the central body.
  • Mass of the central body (M) — the mass of only the central body can be used, as in most cases the mass of the orbiting body is insignificant and does not meaningfully influence the trajectory. However, when this is not the case (e.g.binary stars), the mass of the2-body system can be used instead.

Relations between elements

[edit]

This section contains the common relations between the set of orbital elements described above, but more relations can be derived through manipulations of one or more of these equations. The variable names used here are consistent with the ones described above.

Mean motion can be calculated using the standard gravitational parameter and the semi-major axis of the orbit (μ can be substituted forGM). This equation returns the mean motion in radians, and will need to be converted ifn is desired to be in a different unit.
     n=μa3{\displaystyle n={\sqrt {\frac {\mu }{a^{3}}}}} whena>0{\displaystyle a>0},n=μa3{\displaystyle n={\sqrt {\frac {\mu }{-a^{3}}}}} whena<0{\displaystyle a<0}

Because the semi-major axis is related to the mean motion and standard gravitational parameter, it can be calculated without being specified. This is especially useful ifμ is assumed to be known, as thenn can be used to calculatea, and likewise for specifyinga. This can allow one less element to specified.

Orbital period can be found fromn given the fact that the mean motion can be described as a frequency (number of orbits per unit time), which is the inverse of period.
     P=2πn{\displaystyle P={\frac {2\pi }{n}}}ifn is in radians, orP=360n{\displaystyle P={\frac {360^{\circ }}{n}}} ifn is in degrees.

The standard gravitational parameter can be found given the mean motion and the semi-major axis through the following relation (assuming thatn is in radians):
     μ=n2a3{\displaystyle \mu =n^{2}a^{3}}

The mass of the central body can be found given the standard gravitational parameter using a rearrangement of its definition as the product of the mass and the gravitational constant.
     M=μG{\displaystyle M={\frac {\mu }{G}}}

Epoch describing elements

[edit]

Two elements are needed to describe the position of the body around its orbit, and the time at which this occurs. If this time is defined to be at a point where the specific position variable is a designated constant (usually zero), then it does not need to be specified.

  • Epoch (t0) — time at which one of the below elements is defined. Alternatively this is the point in time where the orbital elements were measured. Sometimes the epoch time is considered as part of the reference frame and is not listed as a distinct element.
  • Time of periapsis passage (T0) — time at which the orbiting body is at periapsis. This is also when the mean anomaly and true anomaly (and others) are zero, so they do not need to be defined. This value is not defined for circular orbits, as they do not have a uniquely defined point of periapsis.
  • Mean anomaly at epoch (M0) — mean anomaly at the epoch time. Mean anomaly is a mathematically convenient angle that increases linearly with time as if the orbit were perfectly circular. Zero is defined as being at periapsis, and one period spans 2π radians. The rate at which the mean anomaly increases is equal to the mean motionn. Because this angle is relative to periapsis, it is not defined for circular orbits.
  • Mean longitude at epoch (L0) — mean longitude at the epoch time. Mean longitude is similar to mean anomaly, in that it increases linearly with time and does not represent the real angular displacement. Unlike with mean anomaly, mean longitude is defined relative to the vernal point, which means it is defined for circular orbits.
  • Eccentric anomaly at epoch (E0) — the eccentric anomaly at the epoch time. Eccentric anomaly is defined at the angular displacement along the auxiliary circle of the ellipse (circle tangent to the ellipse both at apses). This value takes into account the varying speed of objects in elliptical orbits, but does not account for the elliptical shape of the orbit. As such, it still does not correspond to the real angular displacement of the orbiting body. Like with mean anomaly and true anomaly, the eccentric anomaly is measured relative to periapsis, and is not defined for circular orbits. The eccentric anomaly is also not defined for parabolic and hyperbolic trajectories, and instead the parabolic anomaly or hyperbolic anomaly are used.[3]
  • True anomaly at epoch (ν0{\displaystyle \nu _{0}}) — angle that represents the real angular displacement of the orbiting body at the epoch time, taking into account the varying speed and elliptical shape of an orbit. Like with mean anomaly, true anomaly is measured relative to periapsis, and thus it has the same limitations with circular orbits.
  • True longitude at epoch (l0) — the angular displacement of the orbiting body at the epoch time. Unlike with the true anomaly, the true longitude is measured relative to the vernal point, so it can be defined for circular orbits.
  • Mean argument of latitude (uM0) at epoch — the angular displacement of the orbiting body at the epoch time. Mean argument of latitude is similar to the mean anomaly and mean longitude, but instead it is measured relative to the ascending node. This means while it is well defined for circular orbits, it is not for equatorial orbits.[3]
  • Argument of latitude at epoch (u0) — the angular displacement of the orbiting body at the epoch time. This angle is measured relative to the ascending node, so while it is defined for circular orbits, it is not defined for equatorial orbits.

These elements are also used to describe the position of an object in its orbit in a more general context, and are not limited to describing the state at an epoch time.

Relations between elements

[edit]

This section contains the common relations between the set of orbital elements described above, but more relations can be derived through manipulations of one or more of these equations. The variable names used here are consistent with the ones described above. These formulae also hold true for conversions between these elements in general.

Epoch can be found given the time of periapsis passage, the mean anomaly at epoch, and mean motion like so:
     t0=T0+M0n{\displaystyle t_{0}=T_{0}+{\frac {M_{0}}{n}}}

Time of periapsis passage can be found from the epoch, mean anomaly at epoch, and mean motion by re-arranging the previous equation like so:
     T0=t0M0n{\displaystyle T_{0}=t_{0}-{\frac {M_{0}}{n}}}

Mean anomaly can be found from the eccentric anomaly and eccentricity using Kepler's equation like so:
     M=EesinE{\displaystyle M=E-e\sin E}

Mean longitude can be found using the mean anomaly at epoch and the longitude of periapsis.
     L=M+ϖ{\displaystyle L=M+\varpi } orL=M+ω+Ω{\displaystyle L=M+\omega +\Omega }

Eccentric anomaly can be found with the mean anomaly and eccentricity usingKepler's equation through various means, such as iterative calculations or numerical solutions (for some values ofe). Kepler's equation is given as
     E=M+esinE{\displaystyle E=M+e\sin E},

and can be solved through aroot-finding algorithm (usuallyNewton's Method) like so:
     En+1=En+MEn+esin(En)1ecos(En){\displaystyle E_{n+1}=E_{n}+{\frac {M-E_{n}+e\sin(E_{n})}{1-e\cos(E_{n})}}}

Typically a starting guess of eitherM{\displaystyle M},Me{\displaystyle M-e},M+e{\displaystyle M+e}, orM+esinM{\displaystyle M+e\sin M} are used.[3][4] This iteration can be repeated until a desired level of tolerance is reached.

True anomaly can be found from the eccentric anomaly and through the following relations. The quadrant of the solution can be resolved using anatan2(y,x) function.[3]
     sinν=1e2sinE1ecos(E),cosν=cosEe1ecosE{\displaystyle \sin \nu ={\frac {{\sqrt {1-e^{2}}}\sin E}{1-e\cos \left(E\right)}},\cos \nu ={\frac {\cos E-e}{1-e\cos E}}}

True longitude can be found using the true anomaly and longitude of periapsis through the following relation:
     l=ν+ϖ{\displaystyle l=\nu +\varpi } orl=ν+ω+Ω{\displaystyle l=\nu +\omega +\Omega }

Mean argument of latitude can be calculated using the mean anomaly and argument of periapsis like so:
     uM=Ω+M{\displaystyle u_{M}=\Omega +M}

Argument of latitude can be found using the true anomaly and argument of periapsis like so:
     u=ν+Ω{\displaystyle u=\nu +\Omega }

Common sets of elements

[edit]

Classical Keplerian elements

[edit]

While in theory, any set of elements that meets the requirements above can be used to describe an orbit, in practice, certain sets are much more common than others.

The most common elements used to describe the size and shape of the orbit are the semi-major axis (a), and the eccentricity (e). Sometimes the semi-parameter (p) is used instead ofa, as the semi-major axis is infinite for parabolic trajectories, and thus cannot be used.[3][2]

It is common to specify the period (P) or mean motion (n) instead of the semi-major axis in Keplerian element sets, as each can be computed from the other provided thestandard gravitational parameter (μ{\displaystyle \mu }) is known for the central body though the relations above.

For the epoch, the epoch time (t) along with themean anomaly (M0),mean longitude (L0),true anomaly (ν0{\displaystyle \nu _{0}}) or (rarely) theeccentric anomaly (E0) are often used. The time of periapsis passage (T0) is also sometimes used for this purpose.[2]

It is also quite common to see either the mean anomaly or the mean longitude expressed directly, without eitherM0 orL0 as intermediary steps, as alinear function of time. This method of expression will consolidate the mean motion as the slope of this linear equation. An example of this is provided below:
     M(t)=M0+n(tt0){\displaystyle M(t)=M_{0}+n(t-t_{0})}

Elements by body type

[edit]

The choice of elements can differ depending on the type of astronomical body. The eccentricity (e) and either the semi-major axis (a) or the distance of periapsis (q) are used to specify the shape and size of an orbit. The longitude of the ascending node (Ω) the inclination (i) and the argument of periapsis (ω) or thelongitude of periapsis (ϖ) specify the orientation of the orbit in its plane. Either the Mean longitude at epoch (L0) the mean anomaly at epoch (M0) or the time of periapsis passage (T0) are used to specify a known point in the orbit. The choices made depend whether the vernal equinox or the node are used as the primary reference.[5][6]

Sets of orbital elements
ObjectElements used
Major planete,a,i,Ω,ϖ,L0
Comete,q,i, Ω,ω,T0
Asteroide,a,i, Ω,ω,M0

Two-line elements

[edit]
Main article:Two-line element set

Orbital elements can be encoded as text in a number of formats. The most common of them is theNASA /NORAD"two-line elements" (TLE) format,[7] originally designed for use with 80 column punched cards, but still in use because it is the most common format, and 80-character ASCII records can be handled efficiently by modern databases.

The two-line element format lists the eccentricity (e), inclination (i), longitude of the ascending node (Ω), argument of periapsis (ω), mean motion (n), epoch (t0), and mean anomaly at epoch (M0).[7][3] Since the format is primarily meant for orbits of the Earth, the standard gravitational parameter (μ), can be assumed and used to calculate the semi-major axis with the mean motion via the relations above.

Depending on the application and object orbit, the data derived from TLEs older than 30 days can become unreliable. Orbital positions can be calculated from TLEs throughsimplified perturbation models (SGP4 /SDP4 / SGP8 / SDP8).[8]

Example of a two-line element:[9]

1 27651U 03004A   07083.49636287  .00000119  00000-0  30706-4 0  26922 27651 039.9951 132.2059 0025931 073.4582 286.9047 14.81909376225249

Delaunay variables

[edit]

The Delaunay orbital elements were introduced byCharles-Eugène Delaunay during his study of the motion of theMoon.[10] Commonly calledDelaunay variables, they are a set ofcanonical variables, which areaction-angle coordinates. The angles are simple sums of some of the Keplerian angles, and are often referred to with different symbols than other in applications like so:

along with their respectiveconjugate momenta,L,G, andH.[11] The momentaL,G, andH are theaction variables and are more elaborate combinations of the Keplerian elementsa,e, andi.

Delaunay variables are used to simplify perturbative calculations in celestial mechanics, for example while investigating theKozai–Lidov oscillations in hierarchical triple systems.[11] The advantage of the Delaunay variables is that they remain well defined and non-singular (except forh, which can be tolerated) even for circular and equatorial orbits.

Euler angle transformations

[edit]

The anglesΩ,i,ω are theEuler angles (corresponding toα,β,γ in the notation used in that article) characterizing the orientation of the coordinate system

,ŷ, from the inertial coordinate frameÎ,Ĵ,

where:

  • Î,Ĵ is in the equatorial plane of the central body.Î is in the direction of the vernal equinox.Ĵ is perpendicular toÎ and withÎ defines the reference plane. is perpendicular to the reference plane. Orbital elements of bodies (planets, comets, asteroids, ...) in the Solar System usually theecliptic as that plane.
  • ,ŷ are in the orbital plane and with in the direction to thepericenter (periapsis). is perpendicular to the plane of the orbit.ŷ is mutually perpendicular to and.

Then, the transformation from theÎ,Ĵ, coordinate frame to the,ŷ, frame with the Euler anglesΩ,i,ω is:x1=cosΩcosωsinΩcosisinω ;x2=sinΩcosω+cosΩcosisinω ;x3=sinisinω;y1=cosΩsinωsinΩcosicosω ;y2=sinΩsinω+cosΩcosicosω ;y3=sinicosω ;z1=sinisinΩ ;z2=sinicosΩ ;z3=cosi ;{\displaystyle {\begin{aligned}x_{1}&=\cos \Omega \cdot \cos \omega -\sin \Omega \cdot \cos i\cdot \sin \omega \ ;\\x_{2}&=\sin \Omega \cdot \cos \omega +\cos \Omega \cdot \cos i\cdot \sin \omega \ ;\\x_{3}&=\sin i\cdot \sin \omega ;\\\,\\y_{1}&=-\cos \Omega \cdot \sin \omega -\sin \Omega \cdot \cos i\cdot \cos \omega \ ;\\y_{2}&=-\sin \Omega \cdot \sin \omega +\cos \Omega \cdot \cos i\cdot \cos \omega \ ;\\y_{3}&=\sin i\cdot \cos \omega \ ;\\\,\\z_{1}&=\sin i\cdot \sin \Omega \ ;\\z_{2}&=-\sin i\cdot \cos \Omega \ ;\\z_{3}&=\cos i\ ;\\\end{aligned}}}[x1x2x3y1y2y3z1z2z3]=[cosωsinω0sinωcosω0001][1000cosisini0sinicosi][cosΩsinΩ0sinΩcosΩ0001];{\displaystyle {\begin{bmatrix}x_{1}&x_{2}&x_{3}\\y_{1}&y_{2}&y_{3}\\z_{1}&z_{2}&z_{3}\end{bmatrix}}={\begin{bmatrix}\cos \omega &\sin \omega &0\\-\sin \omega &\cos \omega &0\\0&0&1\end{bmatrix}}\,{\begin{bmatrix}1&0&0\\0&\cos i&\sin i\\0&-\sin i&\cos i\end{bmatrix}}\,{\begin{bmatrix}\cos \Omega &\sin \Omega &0\\-\sin \Omega &\cos \Omega &0\\0&0&1\end{bmatrix}}\,;} wherex^=x1I^+x2J^+x3K^ ;y^=y1I^+y2J^+y3K^ ;z^=z1I^+z2J^+z3K^ .{\displaystyle {\begin{aligned}\mathbf {\hat {x}} &=x_{1}\mathbf {\hat {I}} +x_{2}\mathbf {\hat {J}} +x_{3}\mathbf {\hat {K}} ~;\\\mathbf {\hat {y}} &=y_{1}\mathbf {\hat {I}} +y_{2}\mathbf {\hat {J}} +y_{3}\mathbf {\hat {K}} ~;\\\mathbf {\hat {z}} &=z_{1}\mathbf {\hat {I}} +z_{2}\mathbf {\hat {J}} +z_{3}\mathbf {\hat {K}} ~.\\\end{aligned}}}

The inverse transformation, which computes the 3 coordinates in the I-J-K system given the 3 (or 2) coordinates in the x-y-z system, is represented by the inverse matrix. According to the rules ofmatrix algebra, the inverse matrix of the product of the 3 rotation matrices is obtained by inverting the order of the three matrices and switching the signs of the three Euler angles.

That is,

[i1i2i3j1j2j3k1k2k3]=[cosΩsinΩ0sinΩcosΩ0001][1000cosisini0sinicosi][cosωsinω0sinωcosω0001];{\displaystyle {\begin{bmatrix}i_{1}&i_{2}&i_{3}\\j_{1}&j_{2}&j_{3}\\k_{1}&k_{2}&k_{3}\end{bmatrix}}={\begin{bmatrix}\cos \Omega &-\sin \Omega &0\\\sin \Omega &\cos \Omega &0\\0&0&1\end{bmatrix}}\,{\begin{bmatrix}1&0&0\\0&\cos i&-\sin i\\0&\sin i&\cos i\end{bmatrix}}\,{\begin{bmatrix}\cos \omega &-\sin \omega &0\\\sin \omega &\cos \omega &0\\0&0&1\end{bmatrix}}\,;} whereI^=i1x^+i2y^+i3z^ ;J^=j1x^+j2y^+j3z^ ;K^=k1x^+k2y^+k3z^ .{\displaystyle {\begin{aligned}\mathbf {\hat {I}} &=i_{1}\mathbf {\hat {x}} +i_{2}\mathbf {\hat {y}} +i_{3}\mathbf {\hat {z}} ~;\\\mathbf {\hat {J}} &=j_{1}\mathbf {\hat {x}} +j_{2}\mathbf {\hat {y}} +j_{3}\mathbf {\hat {z}} ~;\\\mathbf {\hat {K}} &=k_{1}\mathbf {\hat {x}} +k_{2}\mathbf {\hat {y}} +k_{3}\mathbf {\hat {z}} ~.\\\end{aligned}}}

The transformation from,ŷ, to Euler anglesΩ,i,ω is:Ω=arg(z2,z1)i=arg(z3,z12+z22)ω=arg(y3,x3){\displaystyle {\begin{aligned}\Omega &=\operatorname {arg} \left(-z_{2},z_{1}\right)\\i&=\operatorname {arg} \left(z_{3},{\sqrt {{z_{1}}^{2}+{z_{2}}^{2}}}\right)\\\omega &=\operatorname {arg} \left(y_{3},x_{3}\right)\\\end{aligned}}} wherearg(x,y) signifies the polar argument that can be computed with the standard functionatan2(y,x) available in many programming languages.

Perturbations and elemental variance

[edit]
Main article:Perturbation (astronomy)

Unperturbed,two-body,Newtonian orbits are alwaysconic sections, so the Keplerian elements define an unchangingellipse,parabola, orhyperbola. Real orbits have perturbations, so a given set of Keplerian elements accurately describes an orbit only at the epoch. Evolution of the orbital elements takes place due to thegravitational pull of bodies other than the primary, thenon-sphericity of the primary,atmospheric drag,relativistic effects,radiation pressure,electromagnetic forces, and so on.

Keplerian elements can often be used to produce useful predictions at times near the epoch. Alternatively, real trajectories can be modeled as a sequence of Keplerian orbits thatosculate ("kiss" or touch) the real trajectory. They can also be described by the so-calledplanetary equations, differential equations which come in different forms developed byLagrange,Gauss,Delaunay,Poincaré, orHill.

See also

[edit]

References

[edit]
  1. ^For example, with"VEC2TLE".amsat.org. Archived fromthe original on 20 May 2016. Retrieved19 June 2013.
  2. ^abcdeWeber, Bryan."Orbital Mechanics".Orbital Mechanics. "Orbital Nomenclature", and "Classical Orbital Elements". Retrieved21 February 2025.{{cite web}}: CS1 maint: url-status (link)
  3. ^abcdefghijVallado, David A. (2022).Fundamentals of astrodynamics and applications. Space technology library (4th ed.). Torrance, CA: Microcosm Press. pp. 41–112.ISBN 978-1-881883-18-0.
  4. ^Standish, E. Myles; Williams, James G."Approximate Positions of the Planets".NASA Solar System Dynamics. Retrieved20 February 2025.{{cite web}}: CS1 maint: url-status (link)
  5. ^Green, Robin M. (1985).Spherical Astronomy. Cambridge University Press.ISBN 978-0-521-23988-2.
  6. ^Danby, J. M. A. (1962).Fundamentals of Celestial Mechanics. Willmann-Bell.ISBN 978-0-943396-20-0.{{cite book}}:ISBN / Date incompatibility (help)
  7. ^abKelso, T.S."FAQs: Two-line element set format".celestrak.com. CelesTrak.Archived from the original on 26 March 2016. Retrieved15 June 2016.
  8. ^Seidelmann, K.P., ed. (1992).Explanatory Supplement to the Astronomical Almanac (1st ed.). Mill Valley, CA: University Science Books.
  9. ^"SORCE".Heavens-Above.com. orbit data. Archived fromthe original on 27 September 2007.
  10. ^Aubin, David (2014). "Delaunay, Charles-Eugène".Biographical Encyclopedia of Astronomers. New York City: Springer New York. pp. 548–549.doi:10.1007/978-1-4419-9917-7_347.ISBN 978-1-4419-9916-0.
  11. ^abShevchenko, Ivan (2017).The Lidov–Kozai effect: applications in exoplanet research and dynamical astronomy. Cham: Springer.ISBN 978-3-319-43522-0.

External links

[edit]
Wikibooks has a book on the topic of:Astrodynamics/Classical Orbit Elements
Gravitationalorbits
Types
General
Geocentric
About
other points
Parameters
  • Shape
  • Size
Orientation
Position
Variation
Maneuvers
Orbital
mechanics
Portals:
Retrieved from "https://en.wikipedia.org/w/index.php?title=Orbital_elements&oldid=1287259106"
Category:
Hidden categories:

[8]ページ先頭

©2009-2025 Movatter.jp