轨道测定

轨道测定(OD, Orbit Determination, 亦称轨道确定或轨道决定)是估算行星、小行星、彗星、月球、行星卫星、人造卫星、和太空船等天体绕行其引力源的轨道的技术。透过确定天体的轨道元素，不仅可以推测天体未来的位置，并透过观测来验证；也可以知道还未被发现前的位置。太空船进行星际旅行时，需要不断便换轨道，也需要确定变换后的轨道，以便准确航向目的地。

是取得一系列资料以送入轨道测定的算法。通常一位地基的观测者的观测资料包括时间标记、方位角、高度角、斜距和/或范围率值。因为肉眼的观测不能满足精密定轨的需求，所以都要使用望远镜或雷达装置。

轨道测定之后，数学的推演技术可以用于预测物体未来的轨道位置。随着时间的推移，物体的实际轨道路径往往会偏离预期的路径（尤其是天体的摄动是很难预测的，像是大气阻力等）；新的轨道测定使用新的观测，并有助于重新角准轨道的知识。

美国和做伙的国家，范围广泛的光学、和雷达的资源，允许联合太空作战中心观测与搜集地球轨道上所有的物体。这些观测用于新的轨道计算和测定，以及维护卫星目录的总体精度。防撞计算可以使用这些资料来计算一个轨道上的物体与另一个轨道上的物体碰撞的概率。如果在目前轨道上的碰撞风险是不能接受的，卫星的营运单位可能会决定调整轨道（如果碰撞的概率很低，它是不可能调整轨道的。因为这样做将会导致卫星的推进剂迅速耗尽）。当观测的数量和品质提高，轨道测定技术的准确性也会提高，就会减少提醒卫星营运单位注意的假警报。其它国家，包括俄罗斯和中国，都有类似的追踪资源。

轨道确定中有两个著名的问题，要分别在不同状况下求得绕轨天体的轨道：

高斯问题现在可以用吉伯斯 Gibbs 在 1890 年左右发明的向量进行处理。

轨道测定最基本的原理是由轨道物体的状态向量，即其位置向量及速度向量求出它的几个轨道要素。吉伯斯法主要是透过几个辅助向量，求得物体的速度向量。如此就可以透过状态向量得出轨道要素。以下将介绍关于这几个辅助向量的重要的定理，并与一些轨道参数连结。最终，将求得其速度向量，以便与进行轨道测定。

由轨道焦点连结到3个观察点，可形成 3 个位置向量${\displaystyle \stackrel{\to <!--\; \to \; -->}{r_1}}$, ${\displaystyle \stackrel{\to <!--\; \to \; -->}{r_2}}$ 及 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{r_3}}$ 并定义出轨道平面。这些向量包含充足的轨道讯息，因此可以透过最小平方法来修正轨道的误差。在轨道面上我们可以定义出垂直于这个平面的单位向量 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{k}}$，指向近地点方向的单位向量 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{i}}$ 及与其正交的方向单位向量 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{j}=\stackrel{\to <!--\; \to \; -->}{k}\times <!--\; \times \; -->\stackrel{\to <!--\; \to \; -->}{i}}$。

由三个位置向量所定义出来的 Gauss-Gibbs 向量 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{G}}$ 指向 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{j}}$ 的方向 (半短轴方向)，其中 :

透过牛顿力学解刻普勒问题的过程中，会得出一个离心率向量 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{e}}$，其大小与轨道离心率相同，方向则指向近心点方向，${\displaystyle \stackrel{\to <!--\; \to \; -->}{e}=e\stackrel{\to <!--\; \to \; -->}{i}}$。

在开普勒轨道中，离心率向量是一个不变量，形式为 :

由牛顿力学所推导出来的的轨道极座标式为 :

其中，${\displaystyle re\cos <!--\; \; -->(\mathrm{\theta <!--\; \theta \; -->})=\stackrel{\to <!--\; \to \; -->}{e}\cdot <!--\; \cdot \; -->\stackrel{\to <!--\; \to \; -->}{r}}$因此，有以下关系式 :${\displaystyle p-<!--\; -\; -->r=\stackrel{\to <!--\; \to \; -->}{e}\cdot <!--\; \cdot \; -->\stackrel{\to <!--\; \to \; -->}{r}}$.

计算 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{G}}$ 与 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{e}}$ 的内积，并利用以上关系式，即可证明两者互相垂直。

故 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{G}}$ 与 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{j}}$ (半短轴方向) 同向，可表示为 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{G}=\Vert <!--\; \Vert \; -->\stackrel{\to <!--\; \to \; -->}{G}\Vert <!--\; \Vert \; -->\stackrel{\to <!--\; \to \; -->}{j}}$。

令三个位置向量所定义出来的面积向量为:

则:

且轨道离心率 (eccntricity) 为:

分别计算向量 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{A}}$ 各分量与 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{e}}$ 的外积, 并考虑 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{e}\cdot <!--\; \cdot \; -->\stackrel{\to <!--\; \to \; -->}{r}=p-<!--\; -\; -->r}$, 可得:

将 3 个乘积相加，消去含 ${\displaystyle p}$ 的同值异号分项，并合并相同向量的系数，即得,

由于三个向量 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{r_1}\times <!--\; \times \; -->\stackrel{\to <!--\; \to \; -->}{r_2},\stackrel{\to <!--\; \to \; -->}{r_2}\times <!--\; \times \; -->\stackrel{\to <!--\; \to \; -->}{r_3},\stackrel{\to <!--\; \to \; -->}{r_3}\times <!--\; \times \; -->\stackrel{\to <!--\; \to \; -->}{r_1}}$ 和他们的总和 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{A}}$ 均垂直于轨道平面。所以

移项即可证

最后，令三个位置向量所定义出来的加权体积向量为:

则, 轨道的半正焦弦 (semi-latus rectum) 的长度 ${\displaystyle p}$ 可由体积向量与面积向量求得:

且绕轨物体的相对角动量 (specific angular momentum), ${\displaystyle h}$, 为:

由于三个位置向量共平面，因此它们可以写成：

在此，${\displaystyle \stackrel{\to <!--\; \to \; -->}{k}}$ 是垂直于轨道平面的单位向量，并假设它具有与角动量向量相同的方向：

由于三个向量互相独立，故存在系数 ${\displaystyle ({\mathrm{\lambda <!--\; \lambda \; -->}}_1,{\mathrm{\lambda <!--\; \lambda \; -->}}_2,{\mathrm{\lambda <!--\; \lambda \; -->}}_3)}$，使得它们的线性组合为零向量：

将此方程与 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{e}}$ 求内积, 并考虑 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{e}\cdot <!--\; \cdot \; -->\stackrel{\to <!--\; \to \; -->}{r}=p-<!--\; -\; -->r}$，则有：

可知

如果上式与 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{r_1},\stackrel{\to <!--\; \to \; -->}{r_2},\stackrel{\to <!--\; \to \; -->}{r_3}}$ 分别求外积。则有：

展开上式，且消去同值异号项目，合并系数后, 得:

求解此联立方程式，（并假定 k 为任选的比例常数）可得:

将这些系数代入 ${\displaystyle p}$ 的参数式，结果就是,

又, 由牛顿力学所推导出来的运动轨迹方程式可知,

因此, 透过 ${\displaystyle p}$ 的桥接, 可以得到 ${\displaystyle h}$ 与 的关系：

而 这两个由观测位置所定义出来的辅助向量, 也与轨道的几何性质 ${\displaystyle p}$ 及运动力学的参数 ${\displaystyle h}$ 巧妙地结合在一起。

位置向量所对应的速度向量可以透过离心率向量计算出来。方法是透过 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{k}}$ 与 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{e}}$ 的外积，取得速度向量的表示式。

计算速度向量的步骤如下：

因此，以重力参数和位置向量来表示，我们有以下的速度向量方程式：

(${\displaystyle \mathrm{\mu <!--\; \mu \; -->}}$ 为标准重力参数 standard gravitational parameter).

由前面的定理可知,

且

故

总结以上结果，速度向量 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{v}}$ 与 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{G},\stackrel{\to <!--\; \to \; -->}{A},\stackrel{\to <!--\; \to \; -->}{V}}$ 的关系可表达为以下方程式：

这也可以有另一种证明方式。方法是利用 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{G}=\stackrel{\to <!--\; \to \; -->}{A}\times <!--\; \times \; -->\stackrel{\to <!--\; \to \; -->}{e}}$ 的关系及 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{e}}$ 与 ${\displaystyle }$ 的关系，找出 ${\displaystyle }$ 与 ${\displaystyle }$ 的可能关系。证明如下：

因此，速度向量也可以表示为：

由先前的定理已知 ${\displaystyle h}$ 与 ${\displaystyle }$ 有关。故可代入:

最终，透过三个位置向量所定义的辅助向量 ${\displaystyle \stackrel{\to <!--\; \to \; -->}{G},\stackrel{\to <!--\; \to \; -->}{A},\stackrel{\to <!--\; \to \; -->}{V}}$, 可以将速度向量表示为三个位置向量的函数：

轨道确定的基本任务是由轨道状态向量 ，确定一个轨道物体相对于其中心物体参考框架的古典轨道元素或开普勒元素,${\displaystyle a,e,i,\mathrm{\Omega <!--\; \Omega \; -->},\mathrm{\omega <!--\; \omega \; -->},\mathrm{\nu <!--\; \nu \; -->}}$ 。中心天体是万有引力的来源，如太阳、地球、月球和其他行星。而轨道天体则包括围绕太阳的行星、围绕地球的人造卫星和围绕行星的太空船等。牛顿运动定律对轨道物体的轨迹, 即开普勒轨道, 有很好的解释。

由一个状态向量确定轨道的步骤摘要如下：