0

我正在尝试为我的轨道中的每个传播步骤获取笛卡尔位置和速度向量。我正在使用 OrbitalPy http://pythonhosted.org/OrbitalPy/使用经典的开普勒元素生成轨道。

根据文档，我应该能够从中获取状态向量（位置和速度）class orbital.utilities.StateVector，但我得到一个类型错误：new () 恰好需要 3 个参数（给定 2 个）

这是代码：

from scipy.constants import kilo

import orbital
from orbital import earth, KeplerianElements, Maneuver, plot, utilities
from orbital.utilities import Position, Velocity

import matplotlib.pyplot as plt
import numpy as np

#Orbit Setup
orbitPineapple = KeplerianElements.with_period(96 * 60, body=earth, i=(np.deg2rad(51.6)))
plot(orbitPineapple)
plt.show()
orbitPineapple

Out[23]: KeplerianElements(a=6945033.343911132,
                  e=0,
                  i=0.90058989402907408,
                  raan=0,
                  arg_pe=0,
                  M0=0.0,
                  body=orbital.bodies.earth,
                  ref_epoch=<Time object: scale='utc' format='jyear_str' value=J2000.000>)

prop1 = orbital.maneuver.PropagateAnomalyTo(M=1.00)
orbitX = orbitPineapple.apply_maneuver(prop1)
plot(orbitPineapple, title='Go Pineapple!')
plt.show()

orbital.utilities.StateVector(orbitPineapple)

TypeError                                 Traceback (most recent call last)
<ipython-input-53-91fb5303082b> in <module>()
      4 #print(orbital.utilities.StateVector.velocity(orbitPineapple))
      5 
----> 6 orbital.utilities.StateVector(orbitPineapple)
      7 #orbital.utilities.StateVector.position(orbitPineapple())
      8 
    TypeError: __new__() takes exactly 3 arguments (2 given)
4

2 回答 2

2

我不使用这个包，但错误很简单，可以诊断。从文档中您可以看到orbital.utilities.StateVector需要两个参数；一个用于“位置”，一个用于“速度”。当您这样做时，orbital.utilities.StateVector(orbitPineapple)您只提供一个参数 ( orbitPineapple)，其值将被视为代表“位置”。您还需要提供速度。

至于错误...takes exactly 3 arguments (2 given)，python 高估了类方法所需/传递的参数的数量，因为它在计算self参数时会考虑参数。例如：

class Testing(object):


    def __init__(self):
        self.a = 2

    def do_something(self, b):
        self.a += b

obj = Testing()
obj.do_something(2, 3) # Clearly passing only 2 arguments

给出：

TypeError: do_something() takes exactly 2 arguments (3 given)

因此，您可以将错误解读为“接受 2 个参数，但您只给出了 1 个”

于 2017-09-21T17:53:50.227 回答
0

事实证明，问题出在 OrbitalPy 上。只有在使用原始轨道名称时才能获得状态向量。

在这种情况下orbitPineapple.r，将返回位置 (x,y,z) 并orbitPineapple.v返回 (Vx,Vy,Vy)。

每次机动应用后，位置和速度矢量都会更新，只需使用与原始轨道名称完全相同的行print(orbitPineapple.r, orbitPineapple.v)。

此外，一个非常有用的功能可以为我节省几个小时，您只需键入一个变量或函数，name.然后按Tab 键，所有选项都会显示出来。

于 2017-10-03T19:03:32.407 回答