¿Cuáles son los problemas que deben tenerse en cuenta al tratar de encontrar la posición planetaria mediante el uso de efemérides jpl?

Actualmente estoy tratando de determinar la posición del sol y la luna a partir de jpl ephemeris DE200, que se refiere al ecuador dinámico y al equinoccio de 2000 y usa el tiempo tdb. Estoy usando el sistema j2000 ECI. Las cosas que me preocupan son las siguientes:

1.¿Se da la posición de la Luna con respecto a la Tierra o con respecto al baricentro? Además, cuando encontramos la posición del sol usando la posición del baricentro tierra-luna y la coordenada lunar, ¿será la posición del sol con respecto a la Tierra?

2. Si estoy usando la hora UTC, ¿cómo obtener exactamente el TDB o TAi, ya que está oscilando? (Después de eso, quiero obtener el JD para TDB)

3.Si estoy usando el sistema de coordenadas J2000 ECI, ¿debo tener en cuenta el efecto de la nutación y la precesión? ¿Y como hacer eso?

DE200 es bastante antiguo, es posible que desee utilizar DE431 en su lugar. Podría decirse que la "mejor" manera de usar estos archivos es a través de SPICE: naif.jpl.nasa.gov/naif
Si no recuerdo mal, el libro Astronomical Algorithms de Jean Meeus tiene capítulos sobre esto con descripciones de cómo hacerlo con diferentes niveles de precisión, según sus necesidades.

Respuestas (1)

No menciona con qué está tratando de acceder a esto, pero recomendaría encarecidamente AstroPy para manejar las efemérides y manejar las coordenadas y los sistemas de tiempo. Para sus preguntas:

  1. La Luna y la Tierra se dan en relación con el baricentro Tierra-Luna (centro de masa). Por separado, puede obtener los vectores del Baricentro del Sistema Solar al Baricentro de la Tierra-Luna. Entonces, si desea SSB->Moon, deberá realizar dos llamadas a las efemérides y realizar una suma de vectores.

  2. Depende de la precisión con la que desee las cosas y de la rapidez con la que se mueva lo que desea. Probablemente puedas usar Terrestrial Time (TT) con las efemérides y quieras estar horriblemente fuera de lugar para cualquier otra cosa que no sea la Luna. El desplazamiento TT-UTC (actualmente 67,184 segundos) consta de dos partes; una parte variable de UTC a TAI, actualmente 37 segundos, que consiste (en su mayoría) en el número de segundos bisiestos y una parte fija (32.184 segundos) de TAI a TT. TDB varía alrededor de TT en aproximadamente 2 milisegundos y es algo periódico (dominado por la elipticidad de la órbita de la Tierra y Júpiter y Saturno), pero hacerlo correctamente requiere una expansión de 791 términos (ver, por ejemplo, http://star-www.rl.ac.uk/docs /sun67.htx/sun67ss154.html#Q1-158-946 )

  3. Sí, si desea hacer algo relacionado con la Tierra, como predecir posiciones de algo en el cielo, deberá tener en cuenta el sesgo del marco, la precesión y la nutación junto con la rotación de la Tierra/tiempo sideral para alinear el marco de referencia celeste ( ICRS/J2000) con el de la Tierra (ITRS).

Recomiendo la Circular 179 de USNO ( http://aa.usno.navy.mil/publications/docs/Circular_179.php ), que tiene muchos detalles sobre cómo hacer esto en el nuevo marco de la IAU. El libro de Meeus es muy bueno (también lo tengo) pero se basa en los métodos antiguos de ecuador y equinoccio de fecha IAU1984.

Estoy tratando de implementarlo usando matlab. La idea es que estoy tratando de construir un propagador de órbita para determinar la posición de 2 satélites en el marco ECI, usando como valores iniciales los vectores de posición y velocidad de los satélites. El propagador tiene en cuenta SRP, arrastre, j2 y el tercer cuerpo (sol, luna). Entonces, lo que entiendo de lo que dices es que si quiero tener en cuenta la precesión y la nutación, primero tendré que transformar el posición y velocidades de los satélites en ECEF y después de eso para tener en cuenta la precesión, la nutación y luego volver a ECI.
Eres valiente... ;-) Creo que necesitarás algún tipo de código de precesión y nutación para pasar de un cuadro a otro. Por lo que he visto del código satelital como el de Bill Gray , todos los cálculos son el ecuador medio y el equinoccio de la fecha, por lo que debe ir de una forma u otra a través de la matriz de precesión/nutación o su transposición para comparar/agregar el vector de posición del satélite a Vectores Sol/Luna que en J2000 significarán ecuador/equinoccio (ICRS). (O convierta la pos. Sol/Luna para que signifique el ecuador de la fecha y ECI con un cambio u origen)
¿Cuáles son los problemas que deben tenerse en cuenta al tratar de encontrar la posición planetaria mediante el uso de efemérides jpl?
RespuestasAquíPolítica de privacidad1y, 0v