Solucionador de Órbita de Kepler

Resuelve numéricamente la ecuación de Kepler E − e·sin(E) = M para la anomalía excéntrica E, la anomalía verdadera ν, el radio orbital y la velocidad relativa, con métodos de Newton/bisección/punto fijo, historial de convergencia y un diagrama SVG

Dada la anomalía media M (grados), la excentricidad e y el semieje mayor opcional a (UA), la herramienta resuelve la ecuación trascendente

M = E − e · sin(E)

para la anomalía excéntrica E. A partir de E obtiene:

  • Anomalía verdadera ν: tan(ν/2) = √((1+e)/(1−e)) · tan(E/2)
  • Radio orbital: r = a · (1 − e·cos E)
  • Velocidad relativa (vis-viva): v/v_circ = √(2 − r/a)

Tres métodos:

  • Newton-Raphson (recomendado): derivada analítica f′(E) = 1 − e·cos E, convergencia cuadrática. Para e alta (>0.8) usa E = π como valor inicial.
  • Bisección: robusto, encajonado en [M−e, M+e], convergencia lineal.
  • Punto fijo E ← M + e·sin E: el más simple, pero diverge cuando e → 1.

El informe dibuja la órbita en SVG (elipse centrada en el foco, con perihelio, afelio, foco y posición actual con su radio vector), muestra las magnitudes, clasifica la órbita (circular/elíptica/parabólica/hiperbólica) y el historial completo de iteraciones con residuos.

Resultados de ejemplo

1 Ejemplos

Resolver una órbita elíptica de excentricidad moderada

M=45°, e=0.3, a=1 UA; compara E, ν, r y el historial de Newton.

Kepler solution with orbit SVG and iteration history.
Ver parámetros de entrada
{ "M": 45, "e": 0.3, "a": 1, "method": "newton", "maxIter": 100, "tol": "1e-12" }

Datos clave

Categoría
Geografía y ciencia
Tipos de entrada
number, select
Tipo de salida
html
Cobertura de muestras
4
API disponible
Yes

Resumen

El Solucionador de Órbita de Kepler es una herramienta científica diseñada para resolver numéricamente la ecuación trascendente de Kepler E − e·sin(E) = M. A partir de la anomalía media, la excentricidad y el semieje mayor, calcula con precisión la anomalía excéntrica, la anomalía verdadera, el radio orbital y la velocidad relativa, proporcionando un diagrama orbital en formato SVG y el historial detallado de convergencia mediante métodos como Newton-Raphson, bisección o iteración de punto fijo.

Cuándo usarlo

  • Al calcular la posición y velocidad de un cuerpo celeste en su trayectoria orbital a partir de su anomalía media y excentricidad.
  • Al comparar la velocidad de convergencia y el comportamiento de los métodos numéricos de Newton-Raphson, bisección y punto fijo en órbitas de distinta excentricidad.
  • Al generar representaciones visuales SVG de órbitas elípticas o circulares para informes académicos o de mecánica celeste.

Cómo funciona

  • Introduzca la anomalía media (M) en grados, la excentricidad (e) de la órbita y, opcionalmente, el semieje mayor (a) en unidades astronómicas.
  • Seleccione el método numérico de resolución (Newton-Raphson, bisección o punto fijo), defina el límite de iteraciones y el nivel de tolerancia deseado.
  • La herramienta ejecuta el algoritmo seleccionado para resolver la ecuación trascendente y determinar la anomalía excéntrica (E).
  • Se calculan la anomalía verdadera, el radio orbital y la velocidad relativa, generando un gráfico SVG interactivo de la órbita y la tabla de convergencia.

Casos de uso

Estudiantes y docentes de astrofísica que necesitan visualizar la posición orbital y analizar el residuo de los métodos numéricos paso a paso.
Ingenieros aeroespaciales que requieren calcular de forma rápida la anomalía verdadera y el radio vector de un satélite en una órbita elíptica.
Aficionados a la astronomía que desean modelar órbitas planetarias y obtener la velocidad relativa de un cuerpo en un punto específico de su trayectoria.

Ejemplos

1. Cálculo de posición para un asteroide en órbita excéntrica

Astrofísico computacional
Contexto
Se requiere determinar la posición exacta de un asteroide con una excentricidad moderada de 0.4 y una anomalía media de 60 grados.
Problema
Resolver la ecuación de Kepler con alta precisión para obtener la anomalía verdadera y el radio orbital sin programar un resolvedor desde cero.
Cómo usarlo
Configure la anomalía media en 60, la excentricidad en 0.4 y el semieje mayor en 1.5 UA. Seleccione el método de Newton con una tolerancia de 1e-12.
Configuración de ejemplo
M = 60, e = 0.4, a = 1.5, method = 'newton', tol = '1e-12'
Resultado
Se obtiene la anomalía excéntrica y verdadera exactas, el radio orbital en UA, y un diagrama SVG que muestra la posición del asteroide respecto al foco.

2. Comparación de convergencia en órbita casi circular

Profesor de mecánica celeste
Contexto
Preparando material didáctico para mostrar cómo converge el método de punto fijo frente a Newton en órbitas de baja excentricidad como la de la Tierra.
Problema
Demostrar la diferencia en el número de iteraciones necesarias para alcanzar una tolerancia estricta.
Cómo usarlo
Introduzca una anomalía media de 45 grados y una excentricidad de 0.0167. Ejecute primero con el método de punto fijo y luego con Newton-Raphson.
Configuración de ejemplo
M = 45, e = 0.0167, a = 1.0, method = 'fixedpoint', tol = '1e-15'
Resultado
El sistema genera el historial completo de iteraciones mostrando que, para excentricidades bajas, el punto fijo converge rápidamente en pocos pasos.

Probar con muestras

image, svg, barcode

Hubs relacionados

Preguntas frecuentes

¿Qué método de resolución es el más recomendado?

El método de Newton-Raphson es el recomendado por su convergencia cuadrática rápida, utilizando una aproximación inicial adaptada para excentricidades altas.

¿Por qué diverge el método de punto fijo con excentricidades altas?

La iteración de punto fijo simple para la ecuación de Kepler diverge cuando la excentricidad (e) se aproxima a 1, debido a las limitaciones de convergencia del método en ese rango.

¿Qué ocurre si dejo el semieje mayor (a) en 0?

Si el semieje mayor se define como 0, el radio orbital se expresará y calculará en unidades relativas respecto al propio semieje mayor (a).

¿Cómo se calcula la anomalía verdadera a partir de la anomalía excéntrica?

Se obtiene mediante la relación trigonométrica tan(ν/2) = √((1+e)/(1-e)) * tan(E/2), que vincula directamente ambos ángulos.

¿Qué tipo de órbitas puede clasificar y graficar la herramienta?

Clasifica y representa gráficamente órbitas circulares y elípticas basándose en los valores de excentricidad introducidos.

Documentación de la API

Punto final de la solicitud

POST /es/api/tools/kepler-orbit-solver

Parámetros de la solicitud

Nombre del parámetro Tipo Requerido Descripción
M number -
e number -
a number No -
method select No -
maxIter number No -
tol select No -

Formato de respuesta

{
  "result": "
Processed HTML content
", "error": "Error message (optional)", "message": "Notification message (optional)", "metadata": { "key": "value" } }
HTML: HTML

Documentación de MCP

Agregue este herramienta a su configuración de servidor MCP:

{
  "mcpServers": {
    "elysiatools-kepler-orbit-solver": {
      "name": "kepler-orbit-solver",
      "description": "Resuelve numéricamente la ecuación de Kepler E − e·sin(E) = M para la anomalía excéntrica E, la anomalía verdadera ν, el radio orbital y la velocidad relativa, con métodos de Newton/bisección/punto fijo, historial de convergencia y un diagrama SVG",
      "baseUrl": "https://elysiatools.com/mcp/sse?toolId=kepler-orbit-solver",
      "command": "",
      "args": [],
      "env": {},
      "isActive": true,
      "type": "sse"
    }
  }
}

Puede encadenar múltiples herramientas, por ejemplo: `https://elysiatools.com/mcp/sse?toolId=png-to-webp,jpg-to-webp,gif-to-webp`, máximo 20 herramientas.

Si encuentra algún problema, por favor, póngase en contacto con nosotros en [email protected]