此圖顯示NASA在1962年如何處理軌道測量的數據(只有歷史上的興趣)。 軌道測定(Orbit Determination,亦稱軌道確定或軌道決定)是估算行星、小行星、彗星、月球、行星衛星、人造衛星、和太空船等天體繞行其引力源的軌道的技術。透過確定天體的軌道元素,不僅可以推測天體未來的位置,並透過觀測來驗證;也可以知道還未被發現前的位置。太空船進行星際旅行時,需要不斷變換軌道,也需要確定變換後的軌道,以便準確航向目的地。
觀測是取得一系列資料以送入軌道測定的演算法。通常一位地基的觀測者的觀測資料包括時間標記、方位角、高度角、斜距和/或範圍率值。因為肉眼的觀測不能滿足精密定軌的需求,所以都要使用望遠鏡或雷達裝置。
軌道測定之後,數學的推演技術可以用於預測物體未來的軌道位置。隨著時間的推移,物體的實際軌道路徑往往會偏離預期的路徑(尤其是天體的攝動是很難預測的,像是大氣阻力等);新的軌道測定使用新的觀測,並有助於重新角準軌道的知識。
美國和做伙的國家,範圍廣泛的光學、和雷達的資源,允許聯合太空作戰中心觀測與蒐集地球軌道上所有的物體。這些觀測用於新的軌道計算和測定,以及維護衛星目錄的總體精度。防撞計算可以使用這些資料來計算一個軌道上的物體與另一個軌道上的物體碰撞的概率。如果在目前軌道上的碰撞風險是不能接受的,衛星的營運單位可能會決定調整軌道(如果碰撞的概率很低,它是不可能調整軌道的。因為這樣做將會導致衛星的推進劑迅速耗盡)。當觀測的數量和品質提高,軌道測定技術的準確性也會提高,就會減少提醒衛星營運單位注意的假警報。其它國家,包括俄羅斯和中國,都有類似的追蹤資源。
軌道確定中有兩個著名的問題,要分別在不同狀況下求得繞軌天體的軌道:
- 高斯問題 (the problem of Gauss),是要從已知的 3 個連續位置,
,
和
,確定物體的運動軌道. 在1801年, 高斯透過收集到的三筆數據,於 1801 年 1 月重新找到了天文學家觀測後丟失的穀神星 (Ceres)。此舉讓高斯大大出名。因此,這個問題以他的名字命名。
- 藍伯特問題 (the problem of Lambert),在已知兩組連續位置和日期,{
,
} 及 {
,
},的情況下,求解物體的運動軌道。
高斯問題現在可以用吉伯斯 Gibbs 在 1890 年左右發明的向量進行處理。
軌道測定最基本的原理是由軌道物體的狀態向量,即其位置向量及速度向量求出它的幾個軌道要素。吉伯斯法主要是透過幾個輔助向量,求得物體的速度向量。如此就可以透過狀態向量得出軌道要素。以下將介紹關於這幾個輔助向量的重要的定理,並與一些軌道參數連結。最終,將求得其速度向量,以便與進行軌道測定。
由軌道焦點連結到3個觀察點,可形成 3 個位置向量
,
及
並定義出軌道平面。這些向量包含充足的軌道訊息,因此可以透過最小平方法來修正軌道的誤差。在軌道面上我們可以定義出垂直於這個平面的單位向量
,指向近地點方向的單位向量
及與其正交的方向單位向量
。
由三個位置向量所定義出來的 Gauss-Gibbs 向量
指向
的方向 (半短軸方向),其中 :
.
透過牛頓力學解刻普勒問題的過程中,會得出一個離心率向量
,其大小與軌道離心率相同,方向則指向近心點方向,
。
在克卜勒軌道中,離心率向量是一個不變量,形式為 :
,
為角動量 (angular momentum)
為相對角動量或比角動量 (specific angular momentum)
為標準重力參數 (standard gravitational parameter).
由牛頓力學所推導出來的軌道極座標式為 :
![{\displaystyle r={\frac {p}{1+e\cos(\theta )}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9f47c09eef337a01a12d4a0468f077002d620530)
![{\displaystyle p=r(1+e\cos(\theta ))}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b1e706b4d94f67523919742776af75407578b0d2)
其中,
因此,有以下關係式 :
.
計算
與
的內積,並利用以上關係式,即可證明兩者互相垂直。
.
故
與
(半短軸方向) 同向,可表示為
。
令三個位置向量所定義出來的面積向量為:
![{\displaystyle {\vec {A}}={\vec {r_{1}}}\times {\vec {r_{2}}}+{\vec {r_{2}}}\times {\vec {r_{3}}}+{\vec {r_{3}}}\times {\vec {r_{1}}}=\|{\vec {A}}\|\,{\vec {k}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ee7c1a9497589e6def98c79a0bb0d710e716d7af)
則:
![{\displaystyle {\begin{aligned}{\vec {A}}\times {\vec {e}}={\vec {G}}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6e34bc9fea376196ef81d48cb53ac1f3a891a272)
且軌道離心率 (eccntricity) 為:
.
分別計算向量
各分量與
的外積, 並考慮
, 可得:
.
將 3 個乘積相加,消去含
的同值異號分項,並合併相同向量的係數,即得,
![{\displaystyle {\begin{aligned}{\vec {A}}\times {\vec {e}}&=(r_{2}-r_{3}){\vec {r_{1}}}+(r_{3}-r_{1}){\vec {r_{2}}}+(r_{1}-r_{2}){\vec {r_{3}}}\\&={\vec {G}}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/57557a9c754ddc1475f50ced2c1964201d174cb9)
由於三個向量
和他們的總和
均垂直於軌道平面。所以
![{\displaystyle {\vec {A}}\times {\vec {e}}=\|{\vec {A}}\|\|{\vec {e}}\|\sin(90^{\circ })\,{\vec {j}}=\|{\vec {A}}\|\,\,e\,{\vec {j}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c6b17a9a897eb8075532c932ceae71b0e79a5f8d)
![{\displaystyle \|{\vec {A}}\times {\vec {e}}\|=\|{\vec {A}}\|\,\,e}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e12787045ed3ee6452b1f3cfb636ec904c2509f6)
移項即可証
.
最後,令三個位置向量所定義出來的加權體積向量為:
.
則, 軌道的半正焦弦 (semi-latus rectum) 的長度
可由體積向量與面積向量求得:
.
且繞軌物體的相對角動量 (specific angular momentum),
, 為:
.
由於三個位置向量共平面,因此它們可以寫成:
![{\displaystyle {\begin{aligned}{\vec {r_{1}}}\times {\vec {r_{2}}}&=a_{12}\,{\vec {k}}\Rightarrow a_{12}={\vec {r_{1}}}\times {\vec {r_{2}}}\cdot {\vec {k}}\\{\vec {r_{2}}}\times {\vec {r_{3}}}&=a_{23}\,{\vec {k}}\Rightarrow a_{23}={\vec {r_{2}}}\times {\vec {r_{3}}}\cdot {\vec {k}}\\{\vec {r_{3}}}\times {\vec {r_{1}}}&=a_{31}\,{\vec {k}}\Rightarrow a_{31}={\vec {r_{3}}}\times {\vec {r_{1}}}\cdot {\vec {k}}\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/42375bc0af320dd471a8c36e0abb6d73cce3d489)
在此,
是垂直於軌道平面的單位向量,並假設它具有與角動量向量相同的方向:
.
由於三個向量互相獨立,故存在係數
,使得它們的線性組合為零向量:
.
將此方程與
求內積, 並考慮
,則有:
,
可知
.
如果上式與
分別求外積。則有:
![{\displaystyle \lambda _{1}\,{\vec {r_{1}}}\times {\vec {r_{1}}}+\lambda _{2}\,{\vec {r_{2}}}\times {\vec {r_{1}}}+\lambda _{3}\,{\vec {r_{3}}}\times {\vec {r_{1}}}={\vec {0}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b9ff1b0e3c814b0685e0a6651b142dd0afbeb34b)
![{\displaystyle \lambda _{1}\,{\vec {r_{1}}}\times {\vec {r_{2}}}+\lambda _{2}\,{\vec {r_{2}}}\times {\vec {r_{2}}}+\lambda _{3}\,{\vec {r_{3}}}\times {\vec {r_{2}}}={\vec {0}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a96da26053d9c3ec70ae85d86b6b3101ee8e8835)
,
展開上式,且消去同值異號項目,合併係數後, 得:
.
求解此聯立方程式,(並假定 k 為任選的比例常數)可得:
.
將這些係數代入
的參數式,結果就是,
![{\displaystyle {\begin{aligned}p&={\frac {\lambda _{1}\,r_{1}+\lambda _{2}\,r_{2}+\lambda _{3}\,r_{3}}{\lambda _{1}+\lambda _{2}+\lambda _{3}}}\\&={\frac {({\vec {r_{2}}}\times {\vec {r_{3}}})\,r_{1}+({\vec {r_{3}}}\times {\vec {r_{1}}})\,r_{2}+({\vec {r_{1}}}\times {\vec {r_{2}}})\,r_{3}}{({\vec {r_{2}}}\times {\vec {r_{3}}})+({\vec {r_{3}}}\times {\vec {r_{1}}})+({\vec {r_{1}}}\times {\vec {r_{2}}})}}\\&={\frac {\vec {V}}{\vec {A}}}={\frac {\vec {\|V\|}}{\vec {\|A\|}}}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/78d7557118a073a6371a9b8623d74a750a902a54)
又, 由牛頓力學所推導出來的運動軌跡方程式可知,
![{\displaystyle p={\frac {h^{2}}{\mu }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ab8cbd30ddc42a6ea05d6b29b9f6306611b93bf5)
因此, 透過
的橋接, 可以得到
與 [
] 的關係:
![{\displaystyle {\begin{aligned}p&={\frac {\vec {V}}{\vec {A}}}={\frac {\|{\vec {V}}\|}{\|{\vec {A}}\|}}={\frac {h^{2}}{\mu }}\\\Rightarrow h&={\sqrt {\frac {\mu \|{\vec {V}}\|}{\|{\vec {A}}\|}}}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a49a1acebbd820ca8f5b6d09936dd62d3549558a)
而 [
] 這兩個由觀測位置所定義出來的輔助向量, 也與軌道的幾何性質
及運動力學的參數
巧妙地結合在一起。
位置向量所對應的速度向量可以透過離心率向量計算出來。 方法是透過
與
的外積,取得速度向量的表示式。
計算速度向量的步驟如下:
.
因此,以重力參數和位置向量來表示,我們有以下的速度向量方程式:
![{\displaystyle {\vec {v}}={\frac {\mu }{h}}(e{\vec {j}}+{\vec {k}}\times {\frac {\vec {r}}{r}})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/68c479e0dcb4ea122c6666f08fddafa83650d11a)
(
為標準重力參數 standard gravitational parameter).
由前面的定理可知,
,
且
![{\displaystyle h={\sqrt {\frac {\mu \|{\vec {V}}\|}{\|{\vec {A}}\|}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/522531c16b78a1a89b9a020d80d0a91b3af6c387)
故
.
總結以上結果,速度向量
與
的關係可表達為以下方程式:
.
這也可以有另一種證明方式。方法是利用
的關係及
與
的關係,找出
與
的可能關係。證明如下:
.
因此,速度向量也可以表示為:
.
由先前的定理已知
與
有關。故可代入:
.
最終,透過三個位置向量所定義的輔助向量
, 可以將速度向量表示為三個位置向量的函數:
![{\displaystyle {\begin{aligned}{\vec {v}}&={\frac {\mu }{\|{\vec {A}}\|\cdot h}}({\vec {G}}+{\vec {A}}\times {\frac {\vec {r}}{r}})\\&={\frac {\mu }{\|{\vec {A}}\|}}{\sqrt {\frac {\|{\vec {A}}\|}{\mu \|{\vec {V}}\|}}}({\vec {G}}+{\vec {A}}\times {\frac {\vec {r}}{r}})\\{\vec {v}}&={\sqrt {\frac {\mu }{\|{\vec {A}}\|\|V\|}}}({\vec {G}}+{\vec {A}}\times {\frac {\vec {r}}{r}}).\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a3269ea5399851b8b838035d1a4b92cb93f53b05)
軌道確定的基本任務是由軌道狀態向量 [
],確定一個軌道物體相對於其中心物體參考框架的古典軌道元素或克卜勒元素,
。中心天體是萬有引力的來源,如太陽、地球、月球和其他行星。而軌道天體則包括圍繞太陽的行星、圍繞地球的人造衛星和圍繞行星的太空船等。牛頓運動定律對軌道物體的軌跡, 即克卜勒軌道, 有很好的解釋。
由一個狀態向量確定軌道的步驟摘要如下:
- 由狀態向量計算軌道物體的相對角動量 (比角動量) (specific angular momentum)
:
![{\displaystyle {\vec {h}}={\vec {r}}\times {\vec {v}}=\left|{\vec {h}}\right|{\vec {k}}=h{\vec {k}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0fabfa954c38bb47b3842ec28bdb28f62e705da8)
- 在這裡
是軌道平面 z 軸的單位向量。比角動量是個軌道物體的常數向量。且它的方向垂直於軌道物體的軌道平面。
- 由
計算升交點向量 (ascending node vector)
。假設
代表参考平面Z轴的單位向量,
將垂直於中心體的参考平面,則:
![{\displaystyle {\vec {n}}={\vec {K}}\times {\vec {h}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/34a704cc1d094f7073dff710db61bd40fcb8691c)
- 升交點向量是從中心天體指向軌道平面升交點的向量。由於升交點線是軌道平面和參考平面的交線,故它同時垂直於參考平面向量 (
) 和軌道平面法向量 (
或者
)。因此,升交點向量可以由這兩個向量的外積來定義。
- 計算軌道的離心率向量 (eccentricity vector)
。離心率向量具有軌道離心率的大小,
,並指向軌道近心點的方向。這個方向通常被定義為軌道平面的 x 軸並且有一個單位向量
.根據運動定律,可表示為:
![{\displaystyle {\begin{aligned}{\vec {e}}&={{\vec {v}}\times {\vec {h}} \over {\mu }}-{{\vec {r}} \over {\left|{\vec {r}}\right|}}=e{\vec {i}}\\&=\left({{\left|{\vec {v}}\right|}^{2} \over {\mu }}-{1 \over {\left|{\vec {r}}\right|}}\right){\vec {r}}-{{\vec {r}}\cdot {\vec {v}} \over {\mu }}{\vec {v}}\\&={\frac {1}{\mu }}\left[\left({{\left|{\vec {v}}\right|}^{2}}-{\mu \over {\left|{\vec {r}}\right|}}\right){\vec {r}}-{({\vec {r}}\cdot {\vec {v}})}{\vec {v}}\right]\\e&=\left|{\vec {e}}\right|\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7eb685f04f99f78c67e2d073b04d0c195a695873)
- 在此
是質量為
之中心體的標準引力参數, 而
則是萬有引力常数。
![{\displaystyle p={\frac {h^{2}}{\mu }}=a(1-e^{2})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c9a3a9f04dba48b4598b584e28b9529275eda3d4)
, (如果
)。
![{\displaystyle {\begin{aligned}\cos(i)&={\frac {{\vec {K}}\cdot {\vec {h}}}{h}}={\frac {h_{K}}{h}}\\\Rightarrow i&=\arccos({\frac {{\vec {K}}\cdot {\vec {h}}}{h}}),i\in [0,180^{\circ }],\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0b8303be226f7e50c418498f63f958563bad317c)
- 在此
是
在參考框架的 Z 座標。
![{\displaystyle {\begin{aligned}\cos(\Omega )&={\frac {{\vec {I}}\cdot {\vec {n}}}{n}}={\frac {n_{I}}{n}}=\cos(360-\Omega )\\\Rightarrow \Omega &=\arccos({\frac {{\vec {I}}\cdot {\vec {n}}}{n}})=\Omega _{0},{\text{ or }}\\\Rightarrow \Omega &=360^{\circ }-\Omega _{0},{\text{ if }}n_{J}<0,\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dce38f972b3c34e40ddd9290bbba96685a6002ff)
- 在此
和
分别是
在参考框架中的 X 和 Y 座標。
- 請注意
,但
僅定義在 [0, 180] 度範圍。所以
代表的角度是模棱兩可的,因為在 [0,360] 度中有兩個角度,即
和
,都有相同的
值。所以,實際上它傳回的角度可能是
或者
. 因此,我們必需根據向量在被測量之平面上的 Y 坐標的正負號來進行象限的判斷。在本案例中,
的正負號可用於此例之判斷。
![{\displaystyle {\begin{aligned}\cos(\omega )&={\frac {{\vec {n}}\cdot {\vec {e}}}{ne}}=\cos(360-\omega )\\\Rightarrow \omega &=\arccos({\frac {{\vec {n}}\cdot {\vec {e}}}{ne}})=\omega _{0},{\text{ or }}\\\Rightarrow \omega &=360^{\circ }-\omega _{0},{\text{ if }}e_{K}<0,\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8c3bbd52a965c3fd9826025980bd7c9a928decbd)
- 在此
是
在参考框架中的 Z 座標。
![{\displaystyle {\begin{aligned}\cos(\nu )&={\frac {{\vec {e}}\cdot {\vec {r}}}{er}}=\cos(360-\nu )\\\Rightarrow \nu &=\arccos({\frac {{\vec {e}}\cdot {\vec {r}}}{er}})=\nu _{0},{\text{ or }}\\\Rightarrow \nu &=360^{\circ }-\nu _{0},{\text{ if }}{\vec {r}}\cdot {\vec {v}}<0.\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/423165c460c7112d218d70ee7bf17020e79c020d)
- 在此
的正負號可用於檢查
所在象限, 並修正
傳回的角度,因為它與飛行路徑角 (fly-path angle)
具有相同的正負號 .並且已知,
時,飛行路徑角的符号始终为正, 而當
時,始終為負 [1]。兩者關係在於
, 故
, 與
有相同正負號 .
![{\displaystyle {\begin{aligned}\cos(u)&={\frac {{\vec {n}}\cdot {\vec {r}}}{nr}}=\cos(360-u)\\\Rightarrow u&=\arccos({\frac {{\vec {n}}\cdot {\vec {r}}}{nr}})=u_{0},{\text{ or }}\\\Rightarrow u&=360^{\circ }-u_{0},{\text{ if }}r_{K}<0,\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fb1897366c73d7d956b4cd070317946a14d447b4)
- 在此
是
在参考框架中的 Z 座標。
- Curtis, H.; Orbital Mechanics for Engineering Students, Chapter 5; Elsevier (2005) ISBN 0-7506-6169-0.
- Taff, L.; Celestial Mechanics, Chapters 7, 8; Wiley-Interscience (1985) ISBN 0-471-89316-1.
- Bate, Mueller, White; Fundamentals of Astrodynamics, Chapters 2, 5; Dover (1971) ISBN 0-486-60061-0.
- Madonna, R.; Orbital Mechanics, Chapter 3; Krieger (1997) ISBN 0-89464-010-0.
- Schutz, Tapley, Born; Statistical Orbit Determination, Academic Press. ISBN 978-0126836301
- Orbit Determination and Satellite Navigation
- Satellite Orbit Determination