1

我正在一个地面站上跟踪我们未来的立方体卫星，需要做的一件事是计算它经过头顶时的频率多普勒频移。为了测试这一点，我们正在使用 ISS TLE。我正在使用 Skyfield，它非常有帮助，但我遇到了一个我似乎无法弄清楚的简单问题。我需要物体的速度，这很容易获得，但我还需要相对于我的位置的方向。我假设它是速度矢量的一部分，因为速度既是大小又是方向。也许我在代码中遗漏了一些显而易见的东西，目前修复它的方式是在两个时间点获取距离，并确定它是关闭还是疏远。然后我简单地将向量乘以 -1（如果它接近），否则乘以 1。

diff = satObsDiff.at(ts_now)
diff1 = satObsDiff.at(ts_next)

velocity = diff.speed().km_per_s * 1000 #converts km to m
print("Velocity: ")
print(velocity)
adjusted_velocity = velocity

range1 = diff.distance().km
range2 = diff1.distance().km
change = (range1 - range2)*1000

direction = 1

if change >= 0:
    direction = 1
else:
    direction = -1

observed_freq = ((C/(C + (adjusted_velocity * direction))) * emitted_freq)
4

1 回答 1

1

如果两个向量指向相同的方向，则它们的点积为正，如果它们成直角，则为零，如果它们彼此背离，则它们的点积为负。所以我怀疑你可以通过计算相对位置与速度的点积并检查它是否大于或小于零来判断从一颗后退的卫星接近：

d = np.dot(velocity.km, diff.position.km)
if d > 0:
   print('receding')
else:
   print('approaching')

如果您想知道范围变化的速度，您也可以折腾一些分区：

https://stackoverflow.com/a/55226228/85360

于 2019-11-15T02:08:57.743 回答