Ecuación del centro

Vista simulada de un objeto en una órbita elíptica, como se ve desde el foco de la órbita. La vista gira de acuerdo con la anomalía media, por lo que el objeto parece oscilar hacia adelante y hacia atrás en esta posición media con la ecuación del centro. El objeto también parece hacerse cada vez más pequeño a medida que se aleja más y más debido a la excentricidad de la órbita. Un marcador (rojo) muestra la posición relativa del ápside.

En astrodinámica, es posible abordar el problema de los dos cuerpos Kepleriano mediante la ecuación del centro, considerando la diferencia angular entre la posición real de un cuerpo en su órbita elíptica y la posición que ocuparía si su movimiento fuera uniforme, en una órbita circular del mismo período. Se define como la diferencia entre la anomalía verdadera, Plantilla:Math, y la anomalía media, Plantilla:Math, y se expresa típicamente como una función de la anomalía media, Plantilla:Math y de la excentricidad orbital, Plantilla:Math.

Discusión

Desde la antigüedad, el problema de predecir los movimientos de los cuerpos celestes se ha simplificado al reducirlo a uno de un solo cuerpo en órbita alrededor de otro. Al calcular la posición del cuerpo alrededor de su órbita, a menudo es conveniente comenzar asumiendo un movimiento circular. Esta primera aproximación es simplemente una velocidad angular constante multiplicada por una cantidad de tiempo. Hay varios métodos de proceder para corregir la posición circular aproximada a la producida por el movimiento elíptico, muchos de ellos complejos, y muchos que implican la solución de la Ecuación de Kepler. En contraste, la ecuación del centro es uno de los métodos más fáciles de aplicar.

En casos de excentricidad orbital pequeña, la posición dada por la ecuación del centro puede ser casi tan precisa como cualquier otro método para resolver el problema. Muchas órbitas de interés, como las de los cuerpos en el Sistema solar o las de los satélites artificiales que giran alrededor de la Tierra, tienen casi estas órbitas circulares. A medida que la excentricidad se hace mayor y las órbitas son más elípticas, la exactitud de la ecuación disminuye, fallando por completo en los valores más altos, y por lo tanto, no se usa para tales órbitas.

La ecuación en su forma moderna se puede truncar a cualquier nivel arbitrario de precisión, y cuando se limita a solo los términos más importantes, puede producir una aproximación calculada fácilmente de la posición verdadera cuando la precisión total no es importante. Tales aproximaciones se pueden usar, por ejemplo, como valores iniciales para soluciones iterativas de la ecuación de Kepler, o en el cálculo de tiempos de ajuste, que debido a los efectos atmosféricos no pueden predecirse con mucha precisión.

En la antigua Grecia, en particular Hiparco de Nicea, se conocía la ecuación del centro como prostaféresis, aunque su comprensión de la geometría del movimiento de los planetas no era la misma. La palabra ecuación (del latín, aequatio, -onis) en el sentido presente proviene de la astronomía. Johannes Kepler utilizó el término y lo definió como aquella cantidad variable determinada por el cálculo que debe sumarse o restarse del movimiento promedio para obtener el movimiento verdadero. En astronomía, el término ecuación de tiempo tiene un significado similar. La ecuación del centro en su forma moderna se desarrolló como parte del análisis de perturbaciones, es decir, el estudio de los efectos de un tercer cuerpo en el movimiento de dos cuerpos.

Expansión en serie

Error máximo de la expansión en serie de la ecuación del centro, expresada en radianes, en función de la excentricidad orbital (eje inferior) y de la potencia de la excentricidad orbital e, en la que la serie se trunca (eje derecho). Téngase en cuenta que con pequeñas excentricidades (lado izquierdo del gráfico), la serie no necesita ser desarrollada a potencias elevadas para producir resultados precisos.

Ecuación del centro expandida en serie en función de la anomalía media para varias excentricidades, con la expresión truncada en e⁷ para todas las curvas. Téngase en cuenta que la ecuación truncada falla con excentricidades altas, y produce una curva oscilante.

En el movimiento kepleriano, las coordenadas del cuerpo vuelven a los mismos valores con cada órbita, que es la definición de una función periódica. Dichas funciones se pueden expresar como series de cualquier variable angular con crecimiento continuo, y la variable de mayor interés es la anomalía media, Plantilla:Math. Debido a que aumenta de manera uniforme con el tiempo, expresar cualquier otra variable como una serie en la anomalía media es esencialmente lo mismo que expresarlo en términos del tiempo. Cuando la excentricidad Plantilla:Math de la órbita tiene un valor pequeño, los coeficientes de la serie se pueden desarrollar en términos de potencias de Plantilla:Math. Téngase en cuenta que, si bien estas series se pueden presentar en forma truncada, en realidad representan una suma de un número infinito de términos.

La serie para Plantilla:Math, la anomalía verdadera, se puede expresar de manera más conveniente en términos de Plantilla:Math, Plantilla:Math y una función de Bessel de primera especie,

$\nu =M+2\sum _{s=1}^{\mathrm{\infty }}\frac{1}{s}\{J_s(se)+\sum _{p=1}^{\mathrm{\infty }}{\beta }^p[J_{s-p}(se)+J_{s+p}(se)]\}\sin sM,$

donde

$J_n(se)$ es la función de Bessel y

$\beta =\frac{1}{e}(1-\sqrt{1-e^2}).$

El resultado está en radianes.

Las funciones de Bessel se pueden ampliar en potencias de Plantilla:Math en la forma

$J_n(x)=\frac{1}{n!}{\left(\frac{x}{2}\right)}^n\sum _{m=0}^{\mathrm{\infty }}(-1{)}^m\frac{{\left(\frac{x}{2}\right)}^{2m}}{m!\prod _{k=1}^m(n+k)}$

y Plantilla:Math por,

${\beta }^m={\left(\frac{e}{2}\right)}^m[1+m\sum _{n=1}^{\mathrm{\infty }}\frac{(2n+m-1)!}{n!(n+m)!}{\left(\frac{e}{2}\right)}^{2n}].$

Sustituyendo y reduciendo, la ecuación para Plantilla:Math se vuelve (truncada en el orden Plantilla:Math),

$\begin{array}{cc}\nu =M& +(2e-\frac{1}{4}{e}^3+\frac{5}{96}{e}^5+\frac{107}{4608}{e}^7)\sin M\\ & +(\frac{5}{4}{e}^2-\frac{11}{24}{e}^4+\frac{17}{192}{e}^6)\sin 2M\\ & +(\frac{13}{12}{e}^3-\frac{43}{64}{e}^5+\frac{95}{512}{e}^7)\sin 3M\\ & +(\frac{103}{96}{e}^4-\frac{451}{480}{e}^6)\sin 4M\\ & +(\frac{1097}{960}{e}^5-\frac{5957}{4608}{e}^7)\sin 5M\\ & +\frac{1223}{960}{e}^6\sin 6M+\frac{47273}{32256}{e}^7\sin 7M+\cdots \end{array}$

y según la definición, llevando Plantilla:Math hacia el lado izquierdo de la expresión, Plantilla:Equation box 1 da la ecuación del centro.

Esta ecuación se deduce a veces de una manera alternativa y se presenta en términos de potencias de Plantilla:Math con coeficientes en funciones de $sinM$ (truncadas en el orden Plantilla:Math),

$\begin{array}{cc}\nu =M& +2e\sin M+\frac{5}{4}{e}^2\sin 2M\\ & +\frac{e^3}{12}(13\sin 3M-3\sin M)\\ & +\frac{e^4}{96}(103\sin 4M-44\sin 2M)\\ & +\frac{e^5}{960}(1097\sin 5M-645\sin 3M+50\sin M)\\ & +\frac{e^6}{960}(1223\sin 6M-902\sin 4M+85\sin 2M)+\cdots \end{array}$

que es idéntica a la fórmula anterior.

Para Plantilla:Math pequeña, la serie converge rápidamente. Si Plantilla:Math excede 0.6627..., diverge para algunos valores de Plantilla:Math, efecto descubierto por Pierre-Simon Laplace.

Ejemplos

	Excentricidad orbital	Ecuación del centro máxima (serie truncada como se muestra)
		e7	e3	e2
Venus	0.006777	0.7766°	0.7766°	0.7766°
Tierra	0.01671	1.915°	1.915°	1.915°
Saturno	0.05386	6.174°	6.174°	6.186°
Marte	0.09339	10.71°	10.71°	10.77°
Mercurio	0.2056	23.68°	23.77°	23.28°

Véase también

• Mecánica celeste
• Problema de los dos cuerpos gravitacional
• Órbita de Kepler
• Problema de Kepler
• Problema de los dos cuerpos

Lecturas adicionales

• Marth, A. (1890). On the computation of the equation of the centre in elliptical orbits of moderate eccentricities. Monthly Notices of the Royal Astronomical Society, Vol. 50, p. 502. Gives the equation of the center to order e¹⁰.
• Morrison, J. (1883). On the computation of the eccentric anomaly, equation of the centre and radius vector of a planet, in terms of the mean anomaly and eccentricity. Monthly Notices of the Royal Astronomical Society, Vol. 43, p. 345. Gives the equation of the center to order e¹².
• Morrison, J. (1883). Errata. Monthly Notices of the Royal Astronomical Society, Vol. 43, p. 494.