如何在Pyephem中跟踪特斯拉跑车？

Question

如你所知，SpaceX向太空发射了一辆特斯拉跑车，它现在正在绕太阳运行。

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

我如何在Pyephem中跟踪它？你可以在喷气推进实验室的地平线网站https://ssd.jpl.nasa.gov/horizons.cgi?CGISESSID=c3cbd47fbf603007d1b627107c28962f&s_body=1#top上找到一些轨道元素，搜索特斯拉

看起来相关的数据是：

 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？我对地心轨道器做了很多跟踪，但没有做日心跟踪。

谢谢!

Answer 1

下面是它的工作原理。创建新的EllipticalBody：

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

然后设置参数(来自here和here的参数定义)：

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()

但即便如此，我也不知道M或a应该是什么。HORIZONS应该将这些参数分别输出为MA和A，但从数据中看它们并不明显。理论上，这些可以从TP (近日点朱利安日期)和QR (近日点距离)推导出来，这是给定的(因为HORIZONS根据用户手册从MA和A内部计算它们)，但我不完全确定这是如何做到的。

一些可能性:从XEphem文档中，我学到了一些基本公式：

由于来自地平线的a=1.3299017090067253)

• M=n(E-T) (T是来自HORIZONS)

• n=0.9856076686/P

• P=a**(3/2)

• M=6.991725093611749

的TP )，

• a=q/(1-e) (由于来自地平线的q=QR

问题是，虽然这看起来很好，但当你把它放进去时，数字与地平线所说的相去甚远。