Как я могу отследить родстер Tesla в Pyephem? - PullRequest
Новая
       34

Как я могу отследить родстер Tesla в Pyephem?

0 голосов
/ 14 мая 2018

Как вы, наверное, знаете, SpaceX запустил родстер Tesla в космос, и теперь он вращается вокруг Солнца.

https://en.wikipedia.org/wiki/Elon_Musk%27s_Tesla_Roadster

Как я могу отследить это в Pyephem? Вы можете найти некоторые орбитальные элементы на сайте JPL's Horizons https://ssd.jpl.nasa.gov/horizons.cgi?CGISESSID=c3cbd47fbf603007d1b627107c28962f&s_body=1#top, search for Tesla

и похоже, что соответствующие данные: 

 EPOCH= 2458164.5
  EC= .2585469914787243  QR= .9860596231806226  TP= 2458153.620483722645
  OM= 317.3549094214575  W = 177.3203028023227  IN= 1.088451292866039
  A1= 2.960683526738534E-9  R0= 1.  ALN= 1.  NM= 2.  NK= 0.
  SRC= -2.057839421666802E-7  7.966781900129693E-9  -1.720426606925749E-9
       -4.722542923190676E-7  2.197679131968537E-9  -1.230413802372471E-6
       -2.500290306870021E-7 -3.361070889248183E-9  -1.765963020682463E-5
       -3.047907547965759E-7 -4.640202045440381E-7  -4.271481116360573E-9
        2.657789409005983E-5  1.726818074599357E-6  -1.359673746135991E-6
       -2.478836748687631E-5 -2.309863204867099E-8   -.0002351644867403515
       -1.875169281895894E-6 -2.063647245529267E-6  -1.670539551586607E-6
       -4.019207817588603E-6 -3.128134469402375E-9  -3.034540373576942E-5
        1.733661692209129E-7 -7.052327854535979E-7  -2.650181216776434E-7
       -1.310976135791957E-10
    H= 25.289 G= 0.15

Как мне подключить это к Pyephem? Я много слежу за геоцентрическими орбитальными аппаратами, но ничего гелиоцентрического.

Спасибо!

1 Ответ

0 голосов
/ 12 июля 2018

Вот как это должно работать. Создать новый EllipticalBody: 

tesla=ephem.EllipticalBody()
tesla.name='Tesla Roadster'

затем установите параметры (определения параметров из здесь и здесь ): 

tesla._epoch=2458164.5         #julian date of osculating elements
tesla._e=0.2585469914787243    #eccentricity
tesla._Om='317.3549094214575'  #longitude of ascending node wrt ecliptic
tesla._om='177.3203028023227'  #argument of perihelion wrt ecliptic
tesla._inc='1.088451292866039' #inclination of orbit plane wrt ecliptic
tesla._H=25.289                #absolute magnitude parameter
tesla._G=0.15                  #magnitude slope parameter
tesla._epoch_M=???             #probably the same as tesla._epoch
tesla._M=???                   #unknown quantity, see below
tesla._a=???                   #unknown quantity, see below

Возможно, эпоха неправильна, поскольку J2000 эквивалентен 2451545.0, а в PyEphem - 36525.0, и в этом случае, по моим расчетам, возможно 

tesla._epoch=43144.5

затем запустите compute: 

tesla.compute()

Но даже в этом случае, я не знаю, кем должны быть М или А. ГОРИЗОНТЫ должны выводить эти параметры как MA и A соответственно, но они не очевидны из данных. Теоретически, они могут быть получены из TP (Perihelion Julian Date) и QR (Perihelion Distance), которые даны (поскольку HORIZONS вычисляет их из MA и A для внутреннего использования в руководстве пользователя), но я не полностью уверен, как это сделать.

Некоторые возможности: из документации XEphem я выучил несколько основных формул:

  • a = q / (1-e) (поскольку q = QR из HORIZONS, a = 1.3299017090067253)
  • M = n (E-T) (T - TP от ГОРИЗОНТОВ)
  • п = 0,9856076686 / Р
  • P = а ** (3/2)
  • M = 6,991725093611749

Дело в том, что в то время как все выглядит хорошо и прекрасно, когда вы вставляете это, цифры сильно отличаются от того, что заявляет ГОРИЗОНТ.

...