4. 卫星轨道与可见性窗口计算
测控任务计划的第一道“硬约束”来自几何：卫星什么时候飞到地面站“看得见”的位置。要回答这个问题，必须把轨道预报、升交点、测控弧段、仰角门限和窗口生成串成一条完整的计算链。本节按“先定义参照、再算几何、后生成窗口”的顺序展开，所有示例参数均为教学假设值，用于说明方法而非实测结果。
一、参照约定与基本量
本文统一采用地心惯性系（ECI，J2000）描述轨道，采用地固系（ECEF）描述地面站位置，两者通过地球自转和岁差章动模型转换。卫星位置矢量记作 r_sat，地面站位置矢量记作 r_gs，则站星视线矢量为 ρ = r_sat − r_gs。仰角定义为视线与当地水平面的夹角：
El = arcsin( (ρ · u_up) / |ρ| )
u_up = r_gs / |r_gs| （球模型下的当地垂线方向）方位角 Az 由视线在当地东、北方向上的投影确定，用于判断天线指向和过顶方向。仰角门限 El_min 是地面站能否建立测控链路的最低仰角，典型值取 5°、7° 或 10°，取决于站址遮蔽、频段和大气衰减要求。
二、轨道预报：从两行根数到位置序列
工程上最常用的输入是 TLE（两行根数）。TLE 给出历元时刻的平运动 n、偏心率 e、倾角 i、升交点赤经 Ω、近地点幅角 ω 和平近点角 M。预报流程为：
- 由
n求半长轴a = (μ / n²)^(1/3)，μ 取 398600.4418 km³/s²； - 用开普勒方程
M = E − e·sinE迭代求偏近点角E，再得真近点角ν和轨道半径r； - 在轨道面内得到位置，再经
ω、i、Ω三次旋转转到 ECI； - 考虑 J2 摄动对
Ω、ω、M的长期变化率，按预报时长外推； - 按固定步长（如 10 s 或 30 s）输出 ECI 位置速度序列，再转换到 ECEF。
预报精度取决于 TLE 龄期。对低轨卫星，TLE 龄期超过 3～5 天，位置误差可达数公里量级，会直接平移窗口边界，因此任务前应使用最新 TLE 或精密星历。
三、升交点与轨道面几何
升交点是卫星由南向北穿越赤道面的点，其赤经即 Ω。升交点赤经决定了轨道面在惯性空间中的朝向，而地球在惯性系下自转，因此同一颗卫星相邻两圈的地面轨迹会向西平移约 360° × T_orbit / 86164 s。对周期约 95 min 的低轨卫星，这个西移量约为 24°，是决定地面站每天可见圈次分布的关键量。降交点则用于太阳同步轨道的地方时控制，与测控窗口的昼夜分布直接相关。
四、测控弧段与仰角门限
把预报得到的地固系卫星位置与地面站位置逐点求仰角，仰角从低于门限升到高于门限再降回门限，中间连续满足 El ≥ El_min 的时间段就是一个测控弧段。弧段特征量包括：
- 起始时刻与结束时刻：仰角穿越门限的时刻，用线性插值细化到秒级；
- 最大仰角
El_max：决定链路余量和可用速率； - 弧段时长：低轨单站典型 5～12 min；
- 过顶方位：由方位角变化判断是东侧过境还是西侧过境。
门限抬高会同时压缩弧段时长并推迟起始时刻。以 500 km 圆轨道、10° 门限为例，单站弧段约 8 min；若门限提高到 15°，弧段可能缩短到 5 min 左右。因此门限不是越小越好，需与链路预算和天线跟踪能力匹配。
五、窗口生成方法与冲突消解
窗口生成是把所有地面站、所有可见弧段按时间轴排列，再叠加卫星能源、测控资源、数传需求和任务优先级约束。基本步骤为：
- 对每个站—星组合计算全部弧段，形成候选窗口集合；
- 按
El_max、弧段时长、任务优先级打分排序； - 同一副天线同一时刻只能服务一颗星，做资源冲突检测；
- 对冲突窗口按优先级保留，低优先级窗口回退到次优弧段或次日；
- 输出窗口表，含站名、卫星、起止时刻、最大仰角、可用时长和分配任务。
下面这张图给出从 TLE 到窗口表的完整数据流与关键判断节点，可作为实现时的结构参照。
六、一个可复算的教学示例
设某低轨卫星轨道高度 500 km、倾角 97.4°、升交点赤经 45°，地面站位于北纬 40°、东经 116°、海拔 50 m，仰角门限取 10°。用 30 s 步长粗扫一天，得到该站可见 5～6 圈，单圈弧段约 7～9 min，最大仰角在 40°～85° 之间波动。若把门限改为 5°，同一圈弧段可延长约 1.5～2 min，起始时刻提前约 40 s。这些数值随 TLE 历元和站址变化，实际任务中应以最新星历重算。
七、工程实现要点
实现时建议把“预报—转换—仰角—弧段—窗口”拆成独立函数，便于替换星历模型和门限策略。弧段边界用二分法或线性插值细化到 1 s 以内；对多站多星场景，先用粗步长生成候选区间，再对候选区间细化，可把计算量降低一个数量级。窗口表输出后，还应与卫星能源约束、数传速率和测控任务优先级做二次校验，确认没有把不可执行的弧段写入计划。