Dia Juliano y día J2000
-
m3ntol
- Mensajes: 2058
- Registrado: 04 Ago 2006, 23:00
- Ubicación: Cúmulo de Virgo
- Contactar:
Dia Juliano y día J2000
Hola,
estoy intentando hacerme mis propios programas de astronomía y quiero hacer un algoritmo para pasar de ascensión recta y declinación a altura y azimut.
Tengo un algoritmo que me pide el J2000. Mi problema es que tengo una formula y unas tablas y según emplee unas u otras tengo un desfase de 1.5 días.
Esta es la formula:
d=367.0*Y-((7*(Y+(M+9)/12))/4)+(275*M/9)+D-730530+UT/24;
Y estas las tablas
Tabla A | Tabla B
Días al principio del | Días desde J2000 al
mes | principio de cada año
|
Mes Año Año | Año Días | Año Días
normal bisiesto| |
| |
Ene 0 0 | 1998 -731.5 | 2010 3651.5
Feb 31 31 | 1999 -366.5 | 2011 4016.5
Mar 59 60 | 2000 -1.5 | 2012 4381.5
Abr 90 91 | 2001 364.5 | 2013 4747.5
May 120 121 | 2002 729.5 | 2014 5112.5
Jun 151 152 | 2003 1094.5 | 2015 5477.5
Jul 181 182 | 2004 1459.5 | 2016 5842.5
Ago 212 213 | 2005 1825.5 | 2017 6208.5
Sep 243 244 | 2006 2190.5 | 2018 6573.5
Oct 273 274 | 2007 2555.5 | 2019 6938.5
Nov 304 305 | 2008 2920.5 | 2020 7303.5
Dic 334 335 | 2009 3286.5 | 2021 7669.5
No se cual es el bueno, por ejemplo para el 1 de noviembre de 2000 a las 18:27 con el primero obtengo un J2000 de 306.768750 mientras que con las tablas obtengo 305.26875, justo 1.5 dias menos ¿cual es el bueno?
estoy intentando hacerme mis propios programas de astronomía y quiero hacer un algoritmo para pasar de ascensión recta y declinación a altura y azimut.
Tengo un algoritmo que me pide el J2000. Mi problema es que tengo una formula y unas tablas y según emplee unas u otras tengo un desfase de 1.5 días.
Esta es la formula:
d=367.0*Y-((7*(Y+(M+9)/12))/4)+(275*M/9)+D-730530+UT/24;
Y estas las tablas
Tabla A | Tabla B
Días al principio del | Días desde J2000 al
mes | principio de cada año
|
Mes Año Año | Año Días | Año Días
normal bisiesto| |
| |
Ene 0 0 | 1998 -731.5 | 2010 3651.5
Feb 31 31 | 1999 -366.5 | 2011 4016.5
Mar 59 60 | 2000 -1.5 | 2012 4381.5
Abr 90 91 | 2001 364.5 | 2013 4747.5
May 120 121 | 2002 729.5 | 2014 5112.5
Jun 151 152 | 2003 1094.5 | 2015 5477.5
Jul 181 182 | 2004 1459.5 | 2016 5842.5
Ago 212 213 | 2005 1825.5 | 2017 6208.5
Sep 243 244 | 2006 2190.5 | 2018 6573.5
Oct 273 274 | 2007 2555.5 | 2019 6938.5
Nov 304 305 | 2008 2920.5 | 2020 7303.5
Dic 334 335 | 2009 3286.5 | 2021 7669.5
No se cual es el bueno, por ejemplo para el 1 de noviembre de 2000 a las 18:27 con el primero obtengo un J2000 de 306.768750 mientras que con las tablas obtengo 305.26875, justo 1.5 dias menos ¿cual es el bueno?
-
Jonatan
- Mensajes: 2150
- Registrado: 07 Ago 2005, 23:00
-
acafar
- Moderador
- Mensajes: 3715
- Registrado: 13 Sep 2005, 23:00
- Ubicación: OJ287
Hola,
No he repasado el valor que da la fórmula con tu ejemplo, pero si da el valor 306.768750 es que la fórmula es errónea, tiene alguna errata o te pide el día en calendario juliano y no gregoriano (que no me extrañaría).
En cambio el valor de las tablas sí que creo que es correcto.
Yo he utilizado a veces la fórmula de la wikipedia:
http://en.wikipedia.org/wiki/Julian_day
Para el ejemplo que citas da exactamente lo mismo que las tablas. Es decir,
JDN=2451850,
JD=2451850.26875, y por tanto
J2000 = JD - 2451545.0 = 305.26875
He mirado en alguna página que contiene un programa javascript para esto mismo y también coincide con la fórmula de la wikipedia y con tu resultado mediante las tablas.
Por cierto Jonatan que el valor que tú das es el día juliano correcto sin tener en cuenta la hora minuto y segundo. El J2000 se obtiene del día juliano con fracciones (teniendo en cuenta hora minuto y segundo), restándole el día juliano correspondiente al 1 de enero de 2000 al mediodía, que es 2451545.0.
saludos
No he repasado el valor que da la fórmula con tu ejemplo, pero si da el valor 306.768750 es que la fórmula es errónea, tiene alguna errata o te pide el día en calendario juliano y no gregoriano (que no me extrañaría).
En cambio el valor de las tablas sí que creo que es correcto.
Yo he utilizado a veces la fórmula de la wikipedia:
http://en.wikipedia.org/wiki/Julian_day
Para el ejemplo que citas da exactamente lo mismo que las tablas. Es decir,
JDN=2451850,
JD=2451850.26875, y por tanto
J2000 = JD - 2451545.0 = 305.26875
He mirado en alguna página que contiene un programa javascript para esto mismo y también coincide con la fórmula de la wikipedia y con tu resultado mediante las tablas.
Por cierto Jonatan que el valor que tú das es el día juliano correcto sin tener en cuenta la hora minuto y segundo. El J2000 se obtiene del día juliano con fracciones (teniendo en cuenta hora minuto y segundo), restándole el día juliano correspondiente al 1 de enero de 2000 al mediodía, que es 2451545.0.
saludos
