人造衛星與應用（進階）：擾動力學、軌道設計與鏈路預算

為什麼一顆衛星「放對地方」之後，還是會慢慢跑掉？

入門篇我們把人造衛星想成一顆在理想圓軌道上、繞著一個完美球形地球永遠轉下去的小石頭。這個圖像足以解釋為什麼地球同步衛星掛在赤道上方不動、為什麼低軌衛星每九十多分鐘繞地球一圈。但真實世界遠比這個乾淨的圖像複雜：地球不是完美的球、太陽和月球會拉扯衛星、稀薄的高層大氣會慢慢「吃掉」軌道、太陽光本身的壓力也會推著衛星偏移。

於是任務工程師面對的不是「把衛星放到軌道上」這一個動作，而是一個持續的博弈：軌道每天都在被各種微小的力擾動（perturbation），衛星必須週期性地點火修正，否則幾年內就會偏離崗位、甚至墜入大氣。這一篇我們就把入門篇略過的「擾動力學」攤開來看——它正是太陽同步軌道（sun-synchronous orbit）、衛星星系（constellation）設計、與離軌處置（deorbit）背後的核心。我們也會算一條真正的通訊鏈路預算（link budget），看看為什麼從三萬六千公里外傳回一張照片，訊號弱到難以置信卻仍能被解讀。

擾動力學：理想軌道之外的真實世界

兩體問題（two-body problem）給出的克卜勒橢圓，是把地球當成一個質點或完美球殼的結果。但真實衛星受到的加速度可以寫成

$$\ddot{\vec{r}} = -\frac{\mu}{r^3}\vec{r} + \vec{a}_{\text{pert}}$$

其中第一項是理想的平方反比引力（$\mu = GM_\oplus \approx 3.986\times10^{14}\ \mathrm{m^3/s^2}$），第二項 $\vec{a}_{\text{pert}}$ 是所有擾動的總和。對低軌（LEO）衛星而言，這些擾動依量級大致是：

• 地球扁率（$J_2$）：地球赤道半徑比極半徑大約 21 公里，這個「胖赤道」造成的非球形引力，其擾動加速度約為主引力的 $10^{-3}$，是 LEO 最大的擾動。
• 大氣阻力（atmospheric drag）：在 $400\ \mathrm{km}$ 高度仍有極稀薄的大氣，量級約 $10^{-5}$ 到 $10^{-4}$，且強烈受太陽活動影響。
• 日月引力（third-body）：太陽與月球的潮汐式拉扯，對 LEO 約 $10^{-6}$，但對高軌（GEO）反而變得相對重要。
• 太陽輻射壓（solar radiation pressure, SRP）：陽光光子打在衛星上的動量傳遞，約 $10^{-7}$，對面積大、質量輕的衛星（如太陽帆或大型太陽能板）不可忽略。

這些擾動不是隨機雜訊，而是有清楚物理結構的力。理解它們的關鍵，是不要把它們當成「干擾」要去消除，而是去問：它們會怎樣有系統地改變軌道的形狀與朝向？答案往往出人意料地有用——有些任務正是「故意」利用擾動，把它變成免費的工具。

$J_2$ 擾動：把地球的「胖赤道」變成設計工具

描述地球非球形引力最重要的修正項，是引力位能展開式中的 $J_2$ 項：

$$U(r,\phi) = -\frac{\mu}{r}\left[1 - J_2\left(\frac{R_\oplus}{r}\right)^2 \frac{3\sin^2\phi - 1}{2} + \cdots\right]$$

其中 $\phi$ 是地心緯度，$R_\oplus$ 是地球赤道半徑，係數 $J_2 \approx 1.0826\times10^{-3}$。這個看似很小的數字，卻會在每一圈軌道上累積出可觀的效果。經過軌道平均（averaging over one orbit），可以證明 $J_2$ 不改變半長軸 $a$、離心率 $e$、傾角 $i$ 的長期值，但會讓兩個角度元素持續漂移。

最重要的是升交點赤經（right ascension of the ascending node, $\Omega$）的歲差：

$$\dot{\Omega} = -\frac{3}{2}\,n\,J_2\left(\frac{R_\oplus}{p}\right)^2 \cos i$$

其中 $n = \sqrt{\mu/a^3}$ 是平均角速度，$p = a(1-e^2)$ 是半通徑。這條公式說的是：衛星的軌道平面會繞著地球自轉軸緩慢轉動，轉動方向與速率取決於傾角 $i$。

注意那個 $\cos i$。當 $i < 90^\circ$（順行軌道）時 $\dot\Omega < 0$，軌道面向西退行；當 $i > 90^\circ$（逆行軌道）時 $\dot\Omega > 0$，軌道面向東進動。這就給了任務設計者一個漂亮的把戲：如果我們選一個略大於 $90^\circ$ 的傾角，讓軌道面以恰好等於地球繞太陽公轉的角速率（每天約 $0.9856^\circ$）向東進動，那麼這個軌道面相對於太陽的方向就會永遠固定。

這就是太陽同步軌道（sun-synchronous orbit）。它讓衛星每次經過同一條緯度時，當地的太陽時間都相同——陰影長度一致、光照條件穩定，這正是地球觀測與氣象衛星夢寐以求的。它不是靠燃料維持的，而是把 $J_2$ 擾動「設計進去」的免費禮物。

動手算一下：太陽同步軌道的傾角

我們來算一顆高度 $h = 700\ \mathrm{km}$、近圓軌道（$e \approx 0$）衛星要做太陽同步，需要多大傾角。

先求半長軸：$a = R_\oplus + h = 6378 + 700 = 7078\ \mathrm{km}$。

平均角速度：

$$n = \sqrt{\frac{\mu}{a^3}} = \sqrt{\frac{3.986\times10^{14}}{(7.078\times10^6)^3}} \approx 1.060\times10^{-3}\ \mathrm{rad/s}$$

換算成每天的角度：$n \times 86400 \approx 91.6\ \mathrm{rad/day} \approx 5247^\circ/\mathrm{day}$（即約 14.6 圈/天，合理）。

我們要的歲差率是 $\dot\Omega = +0.9856^\circ/\mathrm{day} = +1.991\times10^{-7}\ \mathrm{rad/s}$。代入歲差公式（$e=0$ 時 $p=a$）：

$$\cos i = -\frac{2\,\dot\Omega}{3\,n\,J_2}\left(\frac{a}{R_\oplus}\right)^2$$

$$\cos i = -\frac{2\times(1.991\times10^{-7})}{3\times(1.060\times10^{-3})\times(1.0826\times10^{-3})}\times\left(\frac{7078}{6378}\right)^2$$

分母為 $3.443\times10^{-6}$，第一個分數約 $-0.1156$，再乘以 $(1.110)^2 = 1.232$，得

$$\cos i \approx -0.1424 \quad\Rightarrow\quad i \approx 98.2^\circ$$

果然是一個略大於 $90^\circ$ 的逆行傾角。這跟現實中 Landsat、Sentinel 系列、福衛五號等地球觀測衛星的傾角（多落在 $97^\circ$–$99^\circ$）完全吻合。一條看似抽象的擾動公式，直接決定了人類拍攝地球的方式。

大氣阻力與軌道衰減：低軌衛星的「壽命時鐘」

對 LEO 衛星，最終決定它能活多久的往往是大氣阻力。阻力加速度為

$$\vec{a}_{\text{drag}} = -\frac{1}{2}\,\rho\,\frac{C_D A}{m}\,v_{\text{rel}}^2\,\hat{v}_{\text{rel}}$$

其中 $\rho$ 是當地大氣密度，$C_D$ 是阻力係數（典型約 $2.2$），$A$ 是迎風截面積，$m$ 是質量。比值 $B = C_D A / m$ 稱為彈道係數（ballistic coefficient）的倒數相關量；面積質量比 $A/m$ 越大，阻力影響越嚴重。

阻力很微妙地與直覺相反：它讓軌道衰減（decay），但會讓衛星加速。原因是阻力消耗的是軌道能量 $E = -\mu/(2a)$，能量降低意味著 $a$ 變小、軌道下降，而由活力公式（vis-viva）

$$v^2 = \mu\left(\frac{2}{r} - \frac{1}{a}\right)$$

可知 $a$ 變小時速度反而增加。換言之，阻力把衛星「往下推」進入更低、更快的軌道，但那裡大氣更稠密、阻力更大，於是進入正回饋——軌道衰減是加速崩潰的過程。這也是為什麼衛星再入大氣的最後幾天會急遽下墜。

關鍵在於大氣密度 $\rho$ 隨高度近似指數衰減：

$$\rho(h) \approx \rho_0 \exp\!\left(-\frac{h - h_0}{H}\right)$$

其中 $H$ 是「尺度高度（scale height）」，在 $400\ \mathrm{km}$ 附近約 $50$–$60\ \mathrm{km}$。更麻煩的是 $\rho_0$ 本身會隨太陽活動（11 年週期）變化達一個數量級以上：太陽活躍期高層大氣受紫外線加熱膨脹，密度大增，衛星壽命可能腰斬。這正是 2022 年 SpaceX 一批星鏈衛星在剛部署時遇上地磁風暴、大氣膨脹而當場損失數十顆的物理機制。

從太空垃圾治理的角度，這條密度曲線也決定了國際規範：通常要求 LEO 衛星任務結束後 25 年內自然再入（近年趨向收緊為 5 年）。設計者要嘛把衛星放在阻力夠大、會自己掉下來的高度，要嘛保留燃料主動降軌（deorbit burn），把它推進濃密大氣燒毀。

通訊鏈路預算：從三萬六千公里外聽見耳語

衛星的價值最終要靠通訊送回地面。從 GEO（$\sim 35\,786\ \mathrm{km}$）傳訊號到地面，會遇到一個殘酷的物理事實：自由空間路徑損耗（free-space path loss）。電磁波從點源向外擴散，功率密度隨距離平方衰減，接收天線只能截到極小一部分。路徑損耗為

$$L_{\text{fs}} = \left(\frac{4\pi d}{\lambda}\right)^2 = \left(\frac{4\pi d f}{c}\right)^2$$

鏈路預算把所有增益與損耗用分貝（dB）相加，核心是接收訊雜比 $C/N_0$，但先看接收功率：

$$P_r = P_t + G_t + G_r - L_{\text{fs}} - L_{\text{other}}$$

（單位皆為 dB / dBW / dBi）。

看一個例子：GEO 下行鏈路

設一顆 GEO 衛星以 $f = 12\ \mathrm{GHz}$（Ku 頻段）下行，發射功率 $P_t = 100\ \mathrm{W} = 20\ \mathrm{dBW}$，衛星天線增益 $G_t = 35\ \mathrm{dBi}$，地面接收碟天線增益 $G_r = 40\ \mathrm{dBi}$，距離 $d = 3.78\times10^7\ \mathrm{m}$。

先算波長：$\lambda = c/f = (3\times10^8)/(12\times10^9) = 0.025\ \mathrm{m}$。

路徑損耗：

$$L_{\text{fs}} = \left(\frac{4\pi \times 3.78\times10^7}{0.025}\right)^2$$

括號內 $= 1.900\times10^{10}$，平方後 $= 3.61\times10^{20}$。換成分貝：

$$L_{\text{fs}} = 10\log_{10}(3.61\times10^{20}) \approx 205.6\ \mathrm{dB}$$

這是一個驚人的數字——訊號被衰減了大約 $10^{20}$ 倍。接收功率：

$$P_r = 20 + 35 + 40 - 205.6 = -110.6\ \mathrm{dBW}$$

即約 $8.7\times10^{-12}\ \mathrm{W}$，不到 $9\ \mathrm{pW}$。這比你手機接收基地台的訊號還弱數百倍。但只要這個功率仍高於接收系統的熱雜訊地板（noise floor），加上現代通道編碼（如 LDPC、Turbo 碼）的糾錯能力，就能可靠地把資料解出來。整個衛星通訊產業，本質上就是在這條鏈路預算的紅線上做極限工程。

這也解釋了為什麼天線增益這麼關鍵：$G_t$ 與 $G_r$ 各貢獻數十 dB，是把訊號從「聽不見」拉回「聽得見」的主力。高增益意味著窄波束，於是又回到指向控制（pointing）的難題——GEO 衛星必須把波束精準對準地球上一塊區域，這又需要靠前面那些擾動修正來維持姿態與位置。所有環節在此扣連起來。

重點回顧

• 真實軌道不是克卜勒橢圓：地球扁率（$J_2$）、大氣阻力、日月引力、太陽輻射壓四種擾動持續改變軌道，量級從 $10^{-3}$ 到 $10^{-7}$ 不等，任務設計就是管理這些擾動。
• $J_2$ 造成軌道面歲差 $\dot\Omega \propto \cos i$，選 $i \approx 98^\circ$ 的逆行軌道可讓軌道面以每天 $0.9856^\circ$ 跟著太陽轉，造就免燃料的太陽同步軌道，地球觀測衛星的標準配置。
• 大氣阻力讓軌道衰減卻使衛星加速：能量耗損使 $a$ 變小，由活力公式速度反增，且更低處密度更大，形成加速崩潰的正回饋；大氣密度隨太陽活動劇烈起伏，直接決定衛星壽命與離軌規範。
• 鏈路預算用 dB 相加：自由空間路徑損耗隨 $d^2$ 與 $f^2$ 增長，GEO 下行可達 $\sim 205\ \mathrm{dB}$，接收功率僅數 pW，靠高增益天線與通道編碼維持可靠通訊。
• 各環節彼此扣連：擾動修正維持指向，指向決定窄波束能否命中地面，鏈路預算決定資料能否傳回——衛星工程是一個耦合系統，不能孤立看任一塊。

深入探討（研究所視角）

把擾動力學做嚴謹的工具，是變分方程（variational equations）與平均化理論（averaging theory）。直接積分含擾動的笛卡兒運動方程在數值上很差，因為主項與擾動項相差數個數量級，誤差累積快。較好的做法是改用軌道根數（orbital elements）作為狀態變數，寫出 Lagrange 行星方程（Lagrange planetary equations） 或對非保守力使用 Gauss 形式的變分方程。它們把每個根數的時間導數表達為擾動位能（或擾動加速度）的函數，例如對升交點赤經

$$\frac{d\Omega}{dt} = \frac{1}{n a^2 \sqrt{1-e^2}\,\sin i}\,\frac{\partial R}{\partial i}$$

其中 $R$ 是擾動位能。把 $J_2$ 的位能代入後對一圈軌道做平均（消去快變的真近點角），就得到前面那條 $\dot\Omega$ 的長期（secular）項。這種「分離快慢變數」的多尺度方法（method of multiple scales），是天體力學與一般非線性動力系統共通的語言。

值得區分三類擾動效應：長期項（secular）隨時間線性累積（如 $J_2$ 的 $\Omega$、$\omega$ 漂移）、長週期項（long-period）以遠長於軌道週期的尺度振盪、短週期項（short-period）在一圈內起伏。任務分析通常只關心長期與長週期項，因為短週期項會自我抵消。高精度應用（如 GNSS 定軌）則需要保留到 $J_4$、$J_6$ 甚至完整的球諧引力場（spherical harmonic gravity field，如 EGM2008 展開到上千階），並納入固體潮、海潮、相對論修正。

最後一個常被忽略卻深刻的環節是相對論時間效應對導航衛星的影響。GPS 衛星上的原子鐘同時受兩種相對論修正：狹義相對論的時間膨脹（衛星高速運動，鐘走慢）約 $-7\ \mu\mathrm{s}/\mathrm{day}$，與廣義相對論的引力紅移（衛星在較弱引力位能處，鐘走快）約 $+45\ \mu\mathrm{s}/\mathrm{day}$，淨效應約 $+38\ \mu\mathrm{s}/\mathrm{day}$。若不修正，定位誤差每天會累積到約 $11\ \mathrm{km}$，導航系統幾小時內就完全失效。工程師的做法是在發射前就把衛星鐘的頻率調慢一個對應量，並在偏心軌道上加入週期性修正項。這是相對論從抽象理論變成每天保障數十億人定位的最直接證據，也是衛星工程裡物理深度的縮影。