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 .
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:
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:
una trayectoria parabólica,
una trayectoria que va y viene a lo largo de un segmento de línea desde el centro de atracción hasta un punto situado a cierta distancia,
una trayectoria que entra o sale a lo largo de un rayo infinito que emana del centro de atracción, con su velocidad yendo a cero con la distancia
o una trayectoria a lo largo de un rayo, pero con una velocidad que no llega a cero con la distancia.
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 .
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-eSin[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.
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 2π ), 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 2π .
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:
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 ,
^ Aaboe, Asger (2001). Episodios de la historia temprana de la astronomía. Springer. pp. 146-147. ISBN978-0-387-95136-2.
^ 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.
^ Colwell, Peter (1993). Resolviendo la ecuación de Kepler a lo largo de tres siglos. Willmann-Bell. pág. 4. ISBN978-0-943396-40-8.
^ 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.
^ North, John (15 de julio de 2008). Cosmos: una historia ilustrada de la astronomía y la cosmología. University of Chicago Press. ISBN978-0-226-59441-5.
^ 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. ISBN978-1-351-58926-0.
^ 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.
^ "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.
^ Fitzpatrick, Philip Matthew (1970). Principios de la mecánica celeste . Academic Press. ISBN0-12-257950-X.
^ 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.
^ 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.
^ 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 )
^ Colwell, Peter (1993). Resolviendo la ecuación de Kepler a lo largo de tres siglos . Willmann–Bell. pág. 43. ISBN0-943396-40-9.
^ 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.
^ 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 .
^ Pfleger, Thomas; Montenbrück, Oliver (1998). Astronomía en la computadora personal (Tercera ed.). Berlín, Heidelberg: Springer. ISBN978-3-662-03349-4.
Enlaces externos
Wikimedia Commons alberga una categoría multimedia sobre Ecuación de Kepler .
Danby, John M.; Burkardt, Thomas M. (1983). "La solución de la ecuación de Kepler. I". Mecánica celeste . 31 (2): 95–107. Código Bibliográfico :1983CeMec..31...95D. doi :10.1007/BF01686811. S2CID 189832421.
Conway, Bruce A. (1986). "Un algoritmo mejorado debido a Laguerre para la solución de la ecuación de Kepler". 24th Aerospace Sciences Meeting . doi :10.2514/6.1986-84.
Mikkola, Seppo (1987). "Una aproximación cúbica para la ecuación de Kepler" (PDF) . Mecánica celeste . 40 (3): 329–334. Código Bibliográfico :1987CeMec..40..329M. doi :10.1007/BF01235850. S2CID 122237945.
Nijenhuis, Albert (1991). "Resolución de la ecuación de Kepler con alta eficiencia y precisión". Mecánica celeste y astronomía dinámica . 51 (4): 319–330. Código Bibliográfico :1991CeMDA..51..319N. doi :10.1007/BF00052925. S2CID 121845017.
Markley, F. Landis (1995). "Solucionador de ecuaciones de Kepler". Mecánica celeste y astronomía dinámica . 63 (1): 101–111. Código Bibliográfico :1995CeMDA..63..101M. doi :10.1007/BF00691917. S2CID 120405765.
Fukushima, Toshio (1996). "Un método para resolver la ecuación de Kepler sin evaluaciones de funciones trascendentales". Mecánica celeste y astronomía dinámica . 66 (3): 309–319. Código Bibliográfico :1996CeMDA..66..309F. doi :10.1007/BF00049384. S2CID 120352687.
Charles, Edgar D.; Tatum, Jeremy B. (1997). "La convergencia de la iteración de Newton-Raphson con la ecuación de Kepler". Mecánica celeste y astronomía dinámica . 69 (4): 357–372. Código Bibliográfico :1997CeMDA..69..357C. doi :10.1023/A:1008200607490. S2CID 118637706.
Stumpf, Laura (1999). "Comportamiento caótico en la función iterativa de Newton asociada con la ecuación de Kepler". Mecánica celeste y astronomía dinámica . 74 (2): 95–109. Bibcode :1999CeMDA..74...95S. doi :10.1023/A:1008339416143. S2CID 122491746.
Palacios, Manuel (2002). "Ecuación de Kepler y método acelerado de Newton". Revista de Matemática Computacional y Aplicada . 138 (2): 335–346. Bibcode :2002JCoAM.138..335P. doi : 10.1016/S0377-0427(01)00369-7 .
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.
Pál, András (2009). "Una solución analítica para el problema de Kepler". Monthly Notices of the Royal Astronomical Society . 396 (3): 1737–1742. arXiv : 0904.0324 . Código Bibliográfico :2009MNRAS.396.1737P. doi : 10.1111/j.1365-2966.2009.14853.x .
Esmaelzadeh, Reza; Ghadiri, Hossein (2014). "Arrancador apropiado para resolver la ecuación de Kepler". Revista Internacional de Aplicaciones Informáticas . 89 (7): 31–38. Bibcode :2014IJCA...89g..31E. doi : 10.5120/15517-4394 .
Zechmeister, Mathias (2018). "Método similar a CORDIC para resolver la ecuación de Kepler". Astronomía y Astrofísica . 619 : A128. arXiv : 1808.07062 . Bibcode :2018A&A...619A.128Z. doi : 10.1051/0004-6361/201833162 .