2. 卫星轨道与地面站几何基础

要判断一颗卫星何时能被某个地面站“看见”,必须先回答三个问题:卫星此刻在哪里、地面站此刻在哪里、两者之间的方向如何随时间变化。这三个问题分别对应轨道六根数、星下点轨迹以及地心惯性系与地固系之间的转换。本节建立这套时空基准,后续所有可见性计算都建立在它之上。

一、轨道六根数:用六个参数锁定轨道

在二体问题假设下,卫星绕地球的运动由六个独立参数完全确定,称为轨道六根数。它们分为两组:一组描述轨道形状与大小,一组描述轨道在空间中的姿态与卫星在轨道上的位置。

  • 半长轴 a:决定轨道尺度,直接关联轨道周期。圆轨道时 a 即轨道半径。
  • 偏心率 e:描述轨道扁平程度,e=0 为圆轨道,0<e<1 为椭圆轨道。
  • 轨道倾角 i:轨道平面与赤道面的夹角,0°为赤道轨道,90°为极轨道。
  • 升交点赤经 Ω:从春分点方向量到升交点方向的赤经角,决定轨道平面绕地轴的旋转位置。
  • 近地点幅角 ω:在轨道平面内,从升交点量到近地点方向的角度。
  • 真近点角 ν:卫星当前位置相对近地点的角度,是随时间变化最快的量。

前五个参数在无摄动时基本固定,第六个参数 ν 随时间推进。工程上常把 ν 换成平近点角 M,通过开普勒方程 M = E − e·sinE 解出偏近点角 E,再换算到 ν,从而得到任意时刻的卫星位置。这一步是后续所有几何计算的起点。

记忆要点:a、e 管“轨道多大、多扁”;i、Ω、ω 管“轨道平面怎么摆、近地点朝哪”;ν(或 M)管“卫星现在跑到哪”。

二、星下点轨迹:把三维轨道投影到地面

星下点是卫星与地心连线与地球表面的交点,其经纬度随时间变化形成的曲线就是星下点轨迹。它把抽象的轨道参数转化为地面站工程师能直接使用的“卫星在地图上的走位”。

星下点纬度由轨道倾角和卫星在轨道上的位置共同决定,最大值等于轨道倾角(对顺行轨道)。星下点经度则受两个因素影响:一是卫星绕地运动带来的经度变化,二是地球自转使地面参考系相对惯性空间转动。因此星下点轨迹通常呈“S”形或波浪形,且相邻两圈轨迹在经度上向西偏移,偏移量约为每圈 360°/圈周期×地球自转周期。

对地面站可见性而言,星下点轨迹的价值在于:地面站与星下点的地心角距离决定了卫星是否在地平线以上。若地面站坐标为 (φs, λs),星下点坐标为 (φu, λu),则地心角 γ 满足:

cos γ = sin φ_s · sin φ_u + cos φ_s · cos φ_u · cos(λ_s − λ_u)

当 γ 小于可见性门限角时,卫星对地面站可见。门限角由地面站最低仰角要求决定,例如要求仰角 ≥ 5°,则门限地心角约为 70° 左右(随轨道高度变化)。

三、地心惯性系与地固系:两套坐标系的转换

轨道六根数描述的是卫星在惯性空间中的运动,通常定义在地心惯性系(ECI)中:原点在地心,Z 轴指向北极,X 轴指向春分点,Y 轴按右手定则确定。而地面站位置、星下点经纬度、天线指向等都在地固系(ECEF)中描述:原点在地心,Z 轴指向北极,X 轴指向本初子午线与赤道交点,随地球一起转动。

两套坐标系之间的转换核心是地球自转角。设某时刻为 t,从 J2000 历元起算的格林尼治平恒星时为 θG(t),则 ECI 到 ECEF 的旋转矩阵为绕 Z 轴旋转 θG:

r_ECEF = R_z(θ_G) · r_ECI
R_z(θ_G) = [ cosθ_G  sinθ_G  0
            -sinθ_G  cosθ_G  0
             0       0       1 ]

实际工程中还需考虑岁差、章动、极移等修正,但在教学与初步可见性分析中,仅用 θG 旋转即可获得足够精度。转换完成后,卫星位置与地面站位置都统一在 ECEF 中,两者相减得到站星矢量,再转换到地面站本地坐标系(东-北-天)即可计算仰角与方位角。

地心惯性系 ECI (不随地球转动,轨道六根数定义于此) Z (北极) X (春分点) Y 地球 轨道面 卫星 升交点 r 地固系 ECEF (随地球转动,地面站与星下点定义于此) Z (北极) X (本初子午线) Y 地球 地面站 星下点 站星矢量 γ R_z(θ_G) 旋转 θ_G ECI → ECEF 转换后,卫星与地面站统一在同一坐标系,方可计算仰角、方位角与可见性

四、从坐标系到可见性判据

把上述三步串起来,就得到可见性计算的标准链路:由轨道六根数在 ECI 中递推卫星位置 → 用 θG 旋转到 ECEF → 与地面站 ECEF 坐标求差得到站星矢量 → 转换到本地东-北-天坐标系 → 计算仰角 el 与方位角 az。当 el 大于等于地面站最低仰角门限时,判定为可见。

本地坐标系下,若站星矢量为 (E, N, U),则仰角 el = arcsin(U / |r|),方位角 az = arctan2(E, N)。这两个量直接决定天线指向,也是后续可见时间窗口计算的输入。

常见错误:把星下点经纬度直接当作卫星位置代入站星矢量计算。星下点只是地心连线与地面的交点,卫星在星下点正上方高度 h 处,必须加上高度分量才能得到真实卫星位置。

五、教学示例:一组可复算的参数

设某地面站位于北纬 30°、东经 120°,海拔 0 km。某卫星轨道六根数为:a = 7000 km,e = 0.001,i = 60°,Ω = 45°,ω = 0°,ν = 30°。在 ECI 中计算卫星位置后,取 θG = 100° 旋转到 ECEF,得到卫星 ECEF 坐标约为 (−3200, 5400, 2800) km(示例近似值)。地面站 ECEF 坐标约为 (−2700, 4700, 3200) km。两者相减得站星矢量,转换到本地坐标系后,可算得仰角约 35°、方位角约 120°。若地面站最低仰角门限为 5°,则该时刻卫星可见。