3. 星地链路几何与坐标系

要分析星地链路的天线增益与指向，第一步不是算增益，而是把“谁在哪里、朝哪个方向”说清楚。卫星在轨道上高速运动，地面站随地球自转，天线又安装在地球表面的某个具体位置，如果坐标系选错或混用，方位角、俯仰角、离轴角都会整体偏移，增益计算随之失真。本节把星地链路中常用的四类坐标系串起来，并给出天线指向角的严格定义。

一、四类坐标系的分工

地心惯性系（ECI）是描述卫星轨道运动的自然坐标系。原点在地心，Z 轴指向北极，X 轴指向春分点，三轴在惯性空间中不随地球转动。卫星的轨道六根数、位置速度矢量通常先在 ECI 中给出，因为开普勒方程和摄动模型都建立在惯性系上。它的特点是“不转”，适合描述卫星自身的运动规律。

地固系（ECEF）随地球一起自转，原点同样在地心，Z 轴指向北极，X 轴指向本初子午线与赤道的交点。地面站坐标、地球模型、地图数据都建立在 ECEF 中。ECI 与 ECEF 之间通过绕 Z 轴旋转一个格林尼治恒星时角（GAST）相互转换，这个角度随时间变化，是星地几何计算中最关键的时变量之一。

站心系（ENU/SEZ）以地面站为原点，东向为 E 轴、北向为 N 轴、天顶为 U 轴，构成右手直角坐标系。它把“地心视角”转换为“站在地面抬头看”的视角，是天线指向角最直接的参考系。卫星在 ECEF 中的位置减去地面站位置，再经一次旋转矩阵投影到 ENU，就得到卫星相对地面站的东、北、天三向分量。

天线本体坐标系则进一步以天线口径面为基准，定义天线电轴的指向。天线增益随离轴角变化，而离轴角正是天线电轴与“地面站到卫星视线”之间的夹角。因此，前三个坐标系解决“视线方向”，本体坐标系解决“天线朝向”，两者之差决定实际增益。

二、从 ECI 到 ENU 的转换链

实际计算中，卫星位置往往以 ECI 给出，而天线指向需要 ENU 下的视线矢量，中间要经过两次关键变换。设卫星在 ECI 中的位置为 r_eci，地面站在 ECEF 中的位置为 r_sta_ecef，则流程如下：

1) ECI → ECEF：
   r_sat_ecef = Rz(GAST) · r_eci
   其中 Rz(θ) 为绕 Z 轴旋转 θ 的矩阵

2) ECEF → ENU（站心）：
   Δr = r_sat_ecef − r_sta_ecef
   [E, N, U]ᵀ = R_enu · Δr
   R_enu 由站址纬度 φ、经度 λ 决定：
   E = −sinλ·Δx + cosλ·Δy
   N = −sinφ·cosλ·Δx − sinφ·sinλ·Δy + cosφ·Δz
   U =  cosφ·cosλ·Δx + cosφ·sinλ·Δy + sinφ·Δz

得到 ENU 分量后，俯仰角和方位角即可直接写出：

El = arcsin( U / |Δr| )
Az = atan2( E, N )   （从北向东为正，范围 0°～360°）
斜距 ρ = |Δr| = sqrt(E² + N² + U²)

这里必须强调符号约定：方位角以正北为 0°、顺时针为正，俯仰角以本地水平面为 0°、天顶为 90°。若采用“以南为 0°”或“逆时针为正”的约定，所有角度都会镜像翻转，后续增益查表必然出错。全文与图示统一采用“北 0°、东 90°、天顶 90°”的约定。

三、天线指向角的定义

天线指向角不是一个角，而是一组角。最核心的是离轴角 θ_off，即天线电轴方向与视线方向之间的夹角。若天线电轴在 ENU 中的单位矢量为 a，视线单位矢量为 u，则：

θ_off = arccos( a · u )

天线增益 G(θ_off) 随离轴角增大而下降，主瓣内近似抛物线，旁瓣则出现起伏。工程上常把“指向误差”定义为电轴与视线的夹角，而把“跟踪误差”定义为电轴与卫星真实方向的夹角。两者在理想情况下相等，但在有轨道预报误差、标定残差或机械回差时会出现差异。

此外还有两个辅助角：方位角 Az 和 俯仰角 El，它们描述视线在站心系中的方向，是天线座架转动的直接指令。对于方位-俯仰座架，Az 和 El 就是两个转轴的目标角；对于 X-Y 座架或极轴座架，则需要再做一次座架坐标变换。无论哪种座架，最终都要回到 θ_off 来评估增益损失。

四、一个可复算的示例

设某地面站位于北纬 30°、东经 120°，海拔 0 km。某时刻卫星在 ECEF 中的位置为（单位：km）：

r_sat_ecef = [ -4500, 5200, 3200 ]
r_sta_ecef ≈ [ -2140, 3706, 3170 ]   （由纬度30°、经度120°、海拔0算得）

则站星矢量 Δr = [ -2360, 1494, 30 ] km。代入 ENU 旋转公式（φ=30°，λ=120°）：

E = −sin120°·(−2360) + cos120°·1494 ≈ 2044 − 747 = 1297 km
N = −sin30°·cos120°·(−2360) − sin30°·sin120°·1494 + cos30°·30
  ≈ −590 − 647 + 26 = −1211 km
U = cos30°·cos120°·(−2360) + cos30°·sin120°·1494 + sin30°·30
  ≈ 1770 + 1940 + 15 = 3725 km

于是斜距 ρ ≈ 4140 km，俯仰角 El = arcsin(3725/4140) ≈ 64.1°，方位角 Az = atan2(1297, −1211) ≈ 133°。若天线电轴恰好对准视线，θ_off = 0；若电轴在方位上偏 0.5°、俯仰上偏 0.3°，则 θ_off ≈ 0.58°，可按天线方向图查得对应的增益下降量。这个示例说明：坐标系转换的每一步都会直接影响最终角度，进而影响增益评估。

五、工程实现中的注意事项

• 时间系统要统一。GAST 依赖 UT1，而轨道预报常用 UTC 或 TAI，转换前必须统一到同一时间尺度，否则自转角误差会直接变成站星几何误差。
• 地球模型要一致。地面站 ECEF 坐标若用 WGS-84 椭球，卫星位置也应使用同一椭球参数，避免高程基准不一致引入几十米量级的偏差。
• 角度单位要显式。三角函数库有弧度与角度之分，建议内部统一用弧度计算，仅在输入输出时转换。
• 奇异点要处理。当卫星接近天顶时，方位角对微小位置误差极其敏感，此时应优先用 θ_off 评估指向，而不是依赖 Az 的数值稳定性。