Como se menciona en esta respuesta , ha habido varias preguntas sobre los datos representados en el sitio whereisroadster dentro del sitio Space Exploration Stack Exchange.
Se eliminó uno sobre distancias incorrectas (pero todavía es visible para usuarios de mayor reputación aquí ), otro sobre roadster parado durante once días en el chat , y un tercero ¿Por qué este gráfico muestra que la velocidad de Starman en relación con el Sol fluctúa tan a menudo? ha sido muy bien recibido y altamente votado. Hay muchas otras preguntas y respuestas que hacen referencia a este sitio aquí también.
Pensé en hacer una revisión al azar para ver cómo está.
Agarré la salida de Horizons para el objeto Roadster como un objetivo utilizando posiciones heliocéntricas y geocéntricas como origen para obtener los vectores de posición y velocidad de esos orígenes. 8 p.m. ( 20:00:00
) aquí en UTC+8 es 12:00:00
UTC, que es cuando la fecha juliana pasa a cero (al menos en decenas de segundos).
Me perdí la captura de pantalla exactamente a la hora, pero clavé la de 20:01:00
. He aplastado las ventanas para que ambas encajen en un ancho de 640 píxeles para que se muestren bien en este sitio SE sin manipulaciones. time.gov es de NIST e intenta sincronizarse a través de la conexión a través del tiempo de ping.
Actualicé las ventanas unos minutos antes solo para asegurarme de que las cosas se actualizaran internamente.
Sin embargo, el sitio whereisroadster obtiene una distancia significativamente diferente entre la Tierra y el Roadster que la que da Horizions, y se supone que este es el motor de datos del que el sitio deriva sus datos.
¿Por qué está tan fuera de lugar?
pos_rel_Earth = np.array([-1.064812011604053E+07, -8.608116821810877E+06, -2.355098230429459E+06])
r_rel_Earth = np.sqrt((pos_rel_Earth**2).sum())
print r_rel_Earth
13893474.1839
Esa es una diferencia del 1% en la distancia, ~100,000 kilómetros
Con respecto al Sol:
pos_rel_Sun = np.array([-1.591862416706663E+08, -2.299106450849950E+07, -2.353784956702188E+06])
r_rel_Sun = np.sqrt((pos_rel_Sun**2).sum())
print r_rel_Sun, v_rel_Sun
160855179.861 31.3143577616
Eso también es una diferencia de ~100,000 kilómetros.
Ephemeris / WWW_USER Mon Mar 26 11:34:36 2018 Pasadena, USA / Horizons
*******************************************************************************
Target body name: SpaceX Roadster (spacecraft) (-143205) {source: tesla_s9}
Center body name: Earth (399) {source: DE431mx}
Center-site name: BODY CENTER
*******************************************************************************
Start time : A.D. 2018-Mar-26 12:00:00.0000 TDB
Stop time : A.D. 2018-Mar-26 12:01:00.0000 TDB
Step-size : 1 minutes
*******************************************************************************
Center geodetic : 0.00000000,0.00000000,0.0000000 {E-lon(deg),Lat(deg),Alt(km)}
Center cylindric: 0.00000000,0.00000000,0.0000000 {E-lon(deg),Dxy(km),Dz(km)}
Center radii : 6378.1 x 6378.1 x 6356.8 km {Equator, meridian, pole}
Output units : KM-S
Output type : GEOMETRIC cartesian states
Output format : 2 (position and velocity)
Reference frame : ICRF/J2000.0
Coordinate systm: Ecliptic and Mean Equinox of Reference Epoch
$$SOE
2458204.000000000, A.D. 2018-Mar-26 12:00:00.0000, -1.064791354593495E+07, -8.608024506708045E+06, -2.355071522438325E+06, -3.442810011828843E+00, -1.538588913186599E+00, -4.451354722214971E-01,
2458204.000694444, A.D. 2018-Mar-26 12:01:00.0000, -1.064812011604053E+07, -8.608116821810877E+06, -2.355098230429459E+06, -3.442860173393305E+00, -1.538581181359120E+00, -4.451308989758778E-01,
$$EOE
Ephemeris / WWW_USER Mon Mar 26 11:55:30 2018 Pasadena, USA / Horizons
*******************************************************************************
Target body name: SpaceX Roadster (spacecraft) (-143205) {source: tesla_s9}
Center body name: Sun (10) {source: DE431mx}
Center-site name: BODY CENTER
*******************************************************************************
Start time : A.D. 2018-Mar-26 12:00:00.0000 TDB
Stop time : A.D. 2018-Mar-26 12:01:00.0000 TDB
Step-size : 1 minutes
*******************************************************************************
Center geodetic : 0.00000000,0.00000000,0.0000000 {E-lon(deg),Lat(deg),Alt(km)}
Center cylindric: 0.00000000,0.00000000,0.0000000 {E-lon(deg),Dxy(km),Dz(km)}
Center radii : 696000.0 x 696000.0 x 696000.0 k{Equator, meridian, pole}
Output units : KM-S
Output type : GEOMETRIC cartesian states
Output format : 2 (position and velocity)
Reference frame : ICRF/J2000.0
Coordinate systm: Ecliptic and Mean Equinox of Reference Epoch
$$SOE
2458204.000000000, A.D. 2018-Mar-26 12:00:00.0000, -1.591861789254073E+08, -2.298918688322970E+07, -2.353758256452293E+06, -1.045906639864063E+00, -3.129377652430419E+01, -4.450064271997949E-01,
2458204.000694444, A.D. 2018-Mar-26 12:01:00.0000, -1.591862416706663E+08, -2.299106450849950E+07, -2.353784956702188E+06, -1.045601993918644E+00, -3.129373246832911E+01, -4.450019025748677E-01,
$$EOE
El sitio es una aproximación. Para reducir la cantidad de datos enviados, cada objeto tiene sus datos almacenados en relación con el Sol. Esto es lo que permite que las parcelas se incluyan con la misma fuente de datos que se utilizó para encontrar la distancia. Los datos se generan utilizando el modo "Vector". Además, utiliza la interpolación lineal entre puntos día a día.
Dicho esto, las cantidades de error son demasiado grandes para que se deban a eso, que en el peor de los casos solo deberían ser unos pocos cientos de km, no 100,000 km o más. Resulta que hay una corrección para las zonas horarias en el código. A la velocidad actual de 3,6 km/s, eso significa que si las zonas horarias son incorrectas, el error podría ser de 12960 km/hora. Elimine el código que corrige incorrectamente las zonas horarias y tendrá el valor correcto.
interp1d
devuelve un útil objeto de interpolación de spline que funciona muy bien!
Edlotiade
UH oh
SF.
UH oh
SF.
UH oh