Question

正如您所知，SpaceX将太空特斯拉跑车推向太空，现在正在绕太阳运行。

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

我如何在Pyephem中追踪这个？您可以在JPL的Horizons网站上找到一些轨道元素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（Perihelion朱利安日期）和QR（Perihelion距离）中得出，它们是给出的（因为HORIZONS根据用户手册在内部从MA和A中计算出了这些值），不确定如何做到这一点。

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

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

情况是，尽管看起来不错，但将其放入时，数字与HORIZONS所说的相去甚远。

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

1 个答案: