2.4 星历与过顶：从两行轨道数据算可见窗口

2.4 星历与过顶：从两行轨道数据算可见窗口

为什么业余爱好者用一台笔记本就能精确预告卫星过顶，误差以秒计？秘密全在两行文本里。本节介绍卫星世界的"时刻表"——星历（TLE，两行轨道数据），用一个简化几何模型亲手算出过顶窗口与多普勒变化，把前三节的轨道、时延、频段参数串成一件趁手的工具。这节是全册动手密度最高的一节，也是理解第 4 章切换预测的前置。

先读懂时刻表：TLE 的两行长什么样

任何公开追踪平台发布的卫星轨道数据长这样：以 "1 " 开头的一行与以 "2 " 开头的一行，每行 69 个字符，每个字符位固定含义。对通信工程真正要紧的字段只有几个：第二行的倾角（决定覆盖纬度带）、升交点赤经（决定轨道面朝向）、平均运动（每天绕地球圈数，倒数即周期）、偏心率（近圆轨道接近 0）；第一行的历元时间（这组参数代表的时刻）与阻力项（用于修正大气拖拽）。专业轨道预报用 SGP4 模型解析这两行数据，精度对过顶预报而言足够到秒级。本节不实现完整 SGP4——那是几百行经过验证的专用代码，任何轨道计算库都有现成实现——而是抽出它背后的圆轨道几何，用一个五十行的简化模型理解过顶的物理骨架，误差在分钟量级，但足以讲清原理并支撑工程估算。

一个必要的实践提醒：TLE 有时效性。低轨卫星受大气阻力与轨道维持机动影响，几周前的 TLE 预报当下过顶可能偏差数十秒到数分钟，做演示足够，做信关站调度要用最新数据并配专业模型。星上侧同理：卫星把自身星历广播给终端（3GPP NTN 里叫星历辅助信息），终端据此预判过顶——这个机制是第 4 章"切换三重难点"的第一根支柱。

图：一次过顶的完整时间窗——仰角与多普勒的同步变化

亲手算窗口：五十行的圆轨道模型

下面的模型做三步：按开普勒第三定律算周期；在时间轴上推进卫星星下点；对每个时刻算观测者的仰角，仰角越过阈值的时段即过顶窗口。为聚焦几何，忽略地球自转（中纬度、十分钟窗口内自转只移动几度，对窗口长度影响是分钟以下的），这在估算层面可接受，正式工具请用专业库的 SGP4 实现。

import math

Re, MU = 6371.0, 398600.0          # 地球半径 km、引力常数

def pass_window(alt_km, incl_deg, raan_deg, lat_deg, lon_deg,
                t0_min=0.0, span_min=90.0, min_elev_deg=25.0, step_s=10):
    """圆轨道近似过顶预报：返回仰角超过阈值的时段"""
    a = Re + alt_km
    T = 2*math.pi*math.sqrt(a**3/MU)                 # 周期(秒)
    inc, raan, lat = map(math.radians, (incl_deg, raan_deg, lat_deg))
    lon0, wins, above = math.radians(lon_deg), [], None
    t = t0_min*60
    while t <= (t0_min+span_min)*60:
        u = 2*math.pi*t/T                            # 相位
        # 轨道坐标系 -> 地心坐标（忽略自转的简化）
        x = math.cos(u)*math.cos(raan) - math.sin(u)*math.sin(raan)*math.cos(inc)
        y = math.cos(u)*math.sin(raan) + math.sin(u)*math.cos(raan)*math.cos(inc)
        z = math.sin(u)*math.sin(inc)
        lat_s, lon_s = math.asin(z), math.atan2(y, x)
        # 星下点与观测者的球面距离 -> 地心角 -> 仰角
        dlon = lon_s - lon0
        gc = math.acos(max(-1,min(1, math.sin(lat)*math.sin(lat_s) +
                 math.cos(lat)*math.cos(lat_s)*math.cos(dlon))))
        rho = math.sqrt(a**2 + Re**2 - 2*a*Re*math.cos(gc))
        # 仰角：asine 公式，由星-地几何直接推出
        elev = math.degrees(math.asin(max(-1, min(1, (a*math.cos(gc) - Re) / rho))))
        if elev >= min_elev_deg and above is None: above = t
        if elev < min_elev_deg and above is not None:
            wins.append((above/60, t/60)); above = None
        t += step_s
    return T/60, wins

T, wins = pass_window(alt_km=550, incl_deg=53, raan_deg=40, lat_deg=30, lon_deg=116)
print(f"轨道周期约 {T:.1f} 分钟")
for s, e in wins:
    print(f"过顶窗口: 第 {s:.1f} ~ {e:.1f} 分钟, 持续约 {e-s:.1f} 分钟")
​

对 550 公里、倾角 53 度的典型参数，这段代码输出的窗口长度在几分钟量级，与 2.3 节"单星可见时间只有几分钟"的说法互相印证。把高度改成 1200 公里重跑，窗口明显变长——高轨星"赖在天上"的时间更多，这正是 OneWeb 每颗星服务时间窗更充裕的几何原因，也是它可以用更少星数的另一半原因。

从窗口到频率：顺手把多普勒也算掉

过顶全程多普勒频移单调地从正穿过零到负，峰值出现在窗口边缘。粗估式为相对速度除以波长：550 公里轨道速度约 7.6 公里每秒，20 GHz 下行波长 1.5 厘米，峰值频移约 500 千赫——链路两端不做预补偿，解调器直接失锁。3GPP NTN 的做法正是本节模型的工程化：终端从星历推算窗口内的多普勒斜率与定时提前斜率，提前装定预补偿参数（第 4.2 节展开）。你在这节写的五十行，就是那套机制的雏形。

常见问题

为什么不直接用现成的轨道计算库？

应该用，正式工程里没有任何理由手写 SGP4。本节手搓简化模型的目的是让机制透明：当你调用库函数拿到过顶窗口时，知道里面是轨道几何、仰角阈值与时间步进三件事的组合。工具越顺滑，理解机制的机会越少——手搓一次是为了以后用得明白。

手机上也能算过顶吗？

能。天文与卫星追踪应用每天都在做同样的事，手机算力绰绰有余。真正要花心思的是信号级应用：过顶预报给出的是"什么时候看得见"，通信终端还需要把窗口细化为"什么时候连得上"——把 4.2 节的参数补偿曲线叠加到窗口上，才是完整的接入计划。

本节要点回顾

• TLE 是行业通用时刻表：倾角、升交点赤经、平均运动、历元四个字段撑起过顶预报。
• 圆轨道模型五十行可写：精度分钟级，足够理解原理与做工程估算，正式调度用 SGP4 与最新星历。
• 窗口长度随高度增长：550 公里约几分钟，1200 公里明显更长，这是星座选高的另一重原因。
• 多普勒从窗口几何直接推：20 GHz 下行峰值约 500 千赫，预补偿参数就藏在你刚算的曲线里。

单颗星的窗口算清了，但用户要的是"永不掉线"——星与星、星与地之间怎么接力？下一章进入星间链路。