stringtranslate.com

Ecuación de Kepler

Soluciones de la ecuación de Kepler para cinco excentricidades diferentes entre 0 y 1

En mecánica orbital , la ecuación de Kepler relaciona varias propiedades geométricas de la órbita de un cuerpo sujeto a una fuerza central .

Fue derivada por Johannes Kepler en 1609 en el Capítulo 60 de su Astronomia nova , [1] [2] y en el libro V de su Epítome de la astronomía copernicana (1621) Kepler propuso una solución iterativa a la ecuación. [3] [4] Esta ecuación y su solución, sin embargo, aparecieron por primera vez en una obra del siglo IX de Habash al-Hasib al-Marwazi , que trataba problemas de paralaje. [5] [6] [7] [8] La ecuación ha jugado un papel importante en la historia de la física y las matemáticas, particularmente en la mecánica celeste clásica .

Ecuación

La ecuación de Kepler es

donde es la anomalía media , es la anomalía excéntrica y es la excentricidad .

La "anomalía excéntrica" ​​es útil para calcular la posición de un punto que se mueve en una órbita kepleriana. Por ejemplo, si el cuerpo pasa por el periastrón en las coordenadas , , en el tiempo , entonces para averiguar la posición del cuerpo en cualquier tiempo, primero se calcula la anomalía media a partir del tiempo y el movimiento medio mediante la fórmula , luego se resuelve la ecuación de Kepler anterior para obtener , luego se obtienen las coordenadas de:

donde está el semieje mayor , el semieje menor .

La ecuación de Kepler es una ecuación trascendental porque el seno es una función trascendental y no se puede resolver algebraicamente . Por lo general, se requieren análisis numéricos y desarrollos en serie para evaluar .

Formas alternativas

Existen varias formas de la ecuación de Kepler. Cada una de ellas está asociada a un tipo específico de órbita. La ecuación estándar de Kepler se utiliza para órbitas elípticas ( ). La ecuación hiperbólica de Kepler se utiliza para trayectorias hiperbólicas ( ). La ecuación radial de Kepler se utiliza para trayectorias lineales (radiales) ( ). La ecuación de Barker se utiliza para trayectorias parabólicas ( ).

Cuando , la órbita es circular. Al aumentar, el círculo se vuelve elíptico. Cuando , hay cuatro posibilidades:

Un valor ligeramente superior a 1 da como resultado una órbita hiperbólica con un ángulo de giro de poco menos de 180 grados. Si se aumenta más, el ángulo de giro se reduce y, a medida que se tiende al infinito, la órbita se convierte en una línea recta de longitud infinita.

Ecuación hiperbólica de Kepler

La ecuación hiperbólica de Kepler es:

donde es la anomalía excéntrica hiperbólica. Esta ecuación se obtiene redefiniendo M como la raíz cuadrada de −1 multiplicado por el lado derecho de la ecuación elíptica:

(en el que ahora es imaginario) y luego reemplazando por .

Ecuaciones radiales de Kepler

La ecuación radial de Kepler para el caso en que el objeto no tiene suficiente energía para escapar es:

donde es proporcional al tiempo y es proporcional a la distancia desde el centro de atracción a lo largo del rayo y alcanza el valor 1 en la distancia máxima. Esta ecuación se obtiene multiplicando la ecuación de Kepler por 1/2 y estableciendo en 1:

y luego hacer la sustitución

La ecuación radial para cuando el objeto tiene suficiente energía para escapar es:

Cuando la energía es exactamente la mínima necesaria para escapar, entonces el tiempo es simplemente proporcional a la distancia elevada a la potencia 3/2.

Problema inverso

Calcular para un valor dado de es sencillo. Sin embargo, resolver para cuando se da puede ser considerablemente más complicado. No existe una solución en forma cerrada . Resolver para es más o menos equivalente a resolver para la anomalía verdadera, o la diferencia entre la anomalía verdadera y la anomalía media, que se denomina " ecuación del centro ".

Se puede escribir una expresión de serie infinita para la solución de la ecuación de Kepler usando la inversión de Lagrange , pero la serie no converge para todas las combinaciones de y (ver más abajo).

La confusión sobre la solubilidad de la ecuación de Kepler ha persistido en la literatura durante cuatro siglos. [9] El propio Kepler expresó dudas sobre la posibilidad de encontrar una solución general:

Estoy suficientemente convencido de que no se puede resolver a priori (la ecuación de Kepler), debido a la diferente naturaleza del arco y del seno. Pero si me equivoco y alguien me indica el camino, será para mí el gran Apolonio .

—  Johannes Kepler [10]

La expansión de la serie de Fourier (con respecto a ) utilizando funciones de Bessel es [11] [12] [13]

Con respecto a , es una serie de Kapteyn .

Ecuación inversa de Kepler

La ecuación inversa de Kepler es la solución de la ecuación de Kepler para todos los valores reales de :

Evaluando esto obtenemos:

Estas series se pueden reproducir en Mathematica con la operación InverseSeries.

InverseSeries[Series[M - Sin[M], {M, 0, 10}]]
InverseSeries[Series[M - e Sin[M], {M, 0, 10}]]

Estas funciones son series simples de Maclaurin . Estas representaciones de funciones trascendentales en series de Taylor se consideran definiciones de esas funciones. Por lo tanto, esta solución es una definición formal de la ecuación inversa de Kepler. Sin embargo, no es una función completa de en un valor dado distinto de cero . De hecho, la derivada

tiende a cero en un conjunto infinito de números complejos cuando el más cercano a cero está en y en estos dos puntos

(donde se toma que el cosh inverso es positivo), y tiende al infinito en estos valores de . Esto significa que el radio de convergencia de la serie de Maclaurin es y la serie no convergerá para valores de mayores que este. La serie también se puede utilizar para el caso hiperbólico, en cuyo caso el radio de convergencia es La serie para cuando converge cuando .

Si bien esta solución es la más sencilla en cierto sentido matemático, [¿ cuál? ] , existen otras soluciones que son preferibles para la mayoría de las aplicaciones. Alternativamente, la ecuación de Kepler se puede resolver numéricamente.

La solución fue encontrada por Karl Stumpff en 1968, [14] pero su importancia no fue reconocida. [15] [ aclaración necesaria ]

También se puede escribir una serie de Maclaurin en . Esta serie no converge cuando es mayor que el límite de Laplace (aproximadamente 0,66), independientemente del valor de (a menos que sea un múltiplo de ), pero converge para todos si es menor que el límite de Laplace. Los coeficientes de la serie, distintos del primero (que es simplemente ), dependen de de forma periódica con período .

Ecuación radial inversa de Kepler

La ecuación radial inversa de Kepler ( ) para el caso en que el objeto no tiene suficiente energía para escapar se puede escribir de manera similar como:

Evaluando esto obtenemos:

Para obtener este resultado usando Mathematica :

InverseSeries[Series[ArcSin[Sqrt[t]] - Sqrt[(1 - t) t], {t, 0, 15}]]

Aproximación numérica del problema inverso

El método de Newton

Para la mayoría de las aplicaciones, el problema inverso se puede calcular numéricamente encontrando la raíz de la función:

Esto se puede hacer iterativamente a través del método de Newton :

Nótese que y están en unidades de radianes en este cálculo. Esta iteración se repite hasta que se obtiene la precisión deseada (por ejemplo, cuando < precisión deseada). Para la mayoría de las órbitas elípticas, un valor inicial de es suficiente. Para órbitas con , se puede utilizar un valor inicial de. Numerosos trabajos desarrollaron conjeturas iniciales precisas (pero también más complejas). [16] Si es idénticamente 1, entonces la derivada de , que está en el denominador del método de Newton, puede acercarse a cero, lo que hace que los métodos basados ​​en derivadas como Newton-Raphson, secante o regula falsi sean numéricamente inestables. En ese caso, el método de bisección proporcionará convergencia garantizada, particularmente porque la solución puede limitarse en un pequeño intervalo inicial. En las computadoras modernas, es posible lograr 4 o 5 dígitos de precisión en 17 a 18 iteraciones. [17] Se puede utilizar un enfoque similar para la forma hiperbólica de la ecuación de Kepler. [18] : 66–67  En el caso de una trayectoria parabólica, se utiliza la ecuación de Barker .

Iteración de punto fijo

Un método relacionado comienza por señalar que . Sustituir repetidamente la expresión de la derecha por la de la derecha produce un algoritmo de iteración de punto fijo simple para evaluar . Este método es idéntico a la solución de Kepler de 1621. [4] En pseudocódigo :

función E ( e , M , n ) E = M para k = 1 a n E = M + e * sen E siguiente k devuelve E                      

El número de iteraciones, , depende del valor de . La forma hiperbólica tiene, de manera similar , .

Este método está relacionado con la solución del método de Newton anterior en que

Para realizar el primer pedido en pequeñas cantidades y ,

.

Véase también

Referencias

  1. ^ Kepler, Johannes (1609). "LX. Methodus, ex hac Physica, hoc est genuina & verissima hipothesi, extruendi utramque partem æquationis, & distantias genuinas: quorum utrumque simul per vicariam fieri hactenus non potuit. argumentum falsæ hipotheseos". Astronomia Nova Aitiologētos, Seu Physica Coelestis, tradita commentariis De Motibus Stellæ Martis, Ex observeibus GV Tychonis Brahe (en latín). págs. 299–300.
  2. ^ Aaboe, Asger (2001). Episodios de la historia temprana de la astronomía. Springer. pp. 146-147. ISBN 978-0-387-95136-2.
  3. ^ Kepler, Johannes (1621). "Libri V. Pars altera". Epitome astronomiæ Copernicanæ usitatâ formâ Quaestionum & Responsionum conscripta, inq; VII. Libros digesta, quorum tres hi priores sunt de Doctrina Sphæricâ (en latín). págs. 695–696.
  4. ^ ab Swerdlow, Noel M. (2000). "Solución iterativa de Kepler a la ecuación de Kepler". Revista de Historia de la Astronomía . 31 (4): 339–341. Código Bibliográfico :2000JHA....31..339S. doi :10.1177/002182860003100404. S2CID  116599258.
  5. ^ Colwell, Peter (1993). Resolviendo la ecuación de Kepler a lo largo de tres siglos. Willmann-Bell. pág. 4. ISBN 978-0-943396-40-8.
  6. ^ Dutka, J. (1 de julio de 1997). "Una nota sobre la "ecuación de Kepler"". Archivo de Historia de las Ciencias Exactas . 51 (1): 59–65. Bibcode :1997AHES...51...59D. doi :10.1007/BF00376451. S2CID  122568981.
  7. ^ North, John (15 de julio de 2008). Cosmos: una historia ilustrada de la astronomía y la cosmología. University of Chicago Press. ISBN 978-0-226-59441-5.
  8. ^ Livingston, John W. (14 de diciembre de 2017). El auge de la ciencia en el Islam y Occidente: de la herencia compartida a la separación de caminos, siglos VIII al XIX. Routledge. ISBN 978-1-351-58926-0.
  9. ^ Se afirma a menudo que la ecuación de Kepler "no puede resolverse analíticamente "; véase, por ejemplo, aquí. Otros autores afirman que no puede resolverse en absoluto; véase, por ejemplo, Madabushi VK Chari; Sheppard Joel Salon; Numerical Methods in Electromagnetism , Academic Press, San Diego, CA, EE. UU., 2000, ISBN 0-12-615760-X , pág. 659. 
  10. ^ "Mihi ſufficit credere, ſolvi a priori non poſſe, propter arcus & ſinus ετερογενειαν. Erranti mihi, quicumque viam monſtraverit, es erit mihi magnus Apollonius". Hall, Asaph (mayo de 1883). "El problema de Kepler". Anales de Matemáticas . 10 (3): 65–66. doi :10.2307/2635832. JSTOR  2635832.
  11. ^ Fitzpatrick, Philip Matthew (1970). Principios de la mecánica celeste . Academic Press. ISBN 0-12-257950-X.
  12. ^ Colwell, Peter (enero de 1992). "Funciones de Bessel y ecuación de Kepler". The American Mathematical Monthly . 99 (1): 45–48. doi :10.2307/2324547. ISSN  0002-9890. JSTOR  2324547.
  13. ^ Boyd, John P. (2007). "Búsqueda de raíces para una ecuación trascendental sin una primera aproximación: polinomio de la ecuación de Kepler a través de la ecuación polinómica del seno de Chebyshev". Matemáticas numéricas aplicadas . 57 (1): 12–18. doi :10.1016/j.apnum.2005.11.010.
  14. ^ Stumpff, Karl (1 de junio de 1968). "Sobre la aplicación de las series de Lie a los problemas de la mecánica celeste". Nota técnica D-4460 de la NASA. {{cite journal}}: Requiere citar revista |journal=( ayuda )
  15. ^ Colwell, Peter (1993). Resolviendo la ecuación de Kepler a lo largo de tres siglos . Willmann–Bell. pág. 43. ISBN 0-943396-40-9.
  16. ^ Odell, AW; Gooding, RH (1986). "Procedimientos para resolver la ecuación de Kepler". Mecánica celeste . 38 (4). Springer Science and Business Media LLC: 307–334. Bibcode :1986CeMec..38..307O. doi :10.1007/bf01238923. ISSN  1572-9478. S2CID  120179781.
  17. ^ Keister, Adrian. "El análisis numérico para hallar la altura de un segmento circular". Wineman Technology . Wineman Technology, Inc . Consultado el 28 de diciembre de 2019 .
  18. ^ Pfleger, Thomas; Montenbrück, Oliver (1998). Astronomía en la computadora personal (Tercera ed.). Berlín, Heidelberg: Springer. ISBN 978-3-662-03349-4.

Enlaces externos