pyephem - 如何追踪位于 Pyephem 的 Tesla Roadster？

Question

如您所知，SpaceX 将特斯拉跑车发射到太空，目前它正在绕太阳运行。

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

我如何在 Pyephem 中追踪这个？您可以在 JPL 的 Horizo​​ns 网站 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

然后运行计算:

tesla.compute()

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

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

• a=q/(1-e)(由于来自 HORIZONS 的 q=QR，a=1.3299017090067253)
• M=n(E-T)(T 是来自 HORIZONS 的 TP)
• n=0.9856076686/P
• P=a**(3/2)
• M=6.991725093611749

事实是，虽然这看起来一切都很好，但当你把它放进去时，数字与 HORIZONS 所说的相差甚远。