3. 星间链路传播时延建模

星间链路的时间同步精度，最终取决于对信号单向传播时延的估计精度。若把两颗卫星之间的传播路径视为一条“时变、色散、非均匀”的信道，那么传播时延 τ(t) 可以拆成四类可分别建模的贡献：几何距离引起的真空传播时延、相对运动引起的时延变化率、多普勒频移对载波相位与码相位的影响，以及大气与电离层折射带来的附加时延。本节按“先建立几何骨架，再叠加动态与介质修正”的顺序展开，并给出可直接用于链路预算的示例参数。

一、几何距离：时延的基准项

设两颗卫星在惯性系（如 J2000）中的位置矢量为 r₁(t) 与 r₂(t)，则瞬时几何距离为 d(t)=|r₂(t)−r₁(t)|，对应的真空传播时延为 τ_geo(t)=d(t)/c，其中 c 为真空光速。这是所有后续修正的基准。以低轨（LEO）与地球静止轨道（GEO）之间的链路为例，若某时刻 d=3.6×10⁷ m，则 τ_geo≈120 ms；若为两颗同轨面 LEO 卫星相距 2000 km，则 τ_geo≈6.67 ms。可见星间链路的时延跨度从毫秒级到百毫秒级，建模时必须先明确轨道类型与几何构型。

实际计算中，位置矢量通常由星历（如广播星历或精密星历）给出，需注意坐标系与时间系统的统一：位置在惯性系，时延在信号发射时刻 t_tx 与接收时刻 t_rx 之间定义。由于卫星在信号传播期间仍在运动，严格求解需迭代：先以 t_tx 时刻位置估算 d，得到 τ，再以 t_rx=t_tx+τ 时刻位置修正，通常 2～3 次迭代即可收敛到亚毫米级。

二、相对运动：时延变化率与一阶修正

星间链路的两颗卫星存在相对速度 v_rel，其沿视线方向的分量 v_los 决定了几何时延的变化率：dτ/dt = v_los/c。若 v_los=3 km/s，则时延变化率约为 10 μs/s，即每秒钟时延变化 10 微秒。对于需要亚纳秒级同步的链路，这一项必须在时间同步算法中实时补偿，否则会引入与时间成正比的相位漂移。

更精确地，接收时刻的时延可展开为 τ(t_rx)≈τ(t_tx)+ (v_los/c)·τ(t_tx)，其中第二项即一阶相对运动修正。若相对加速度较大（如高动态编队），还需引入二阶项。工程上常把 τ(t) 建模为多项式或采用卡尔曼滤波实时估计，以跟踪其变化。

三、多普勒效应：频域与时域的双重影响

相对运动不仅改变时延，还产生多普勒频移 f_d = −(v_los/c)·f_c，其中 f_c 为载波频率。以 Ka 频段 26 GHz、v_los=3 km/s 为例，f_d≈−260 kHz。多普勒对时间同步的影响体现在两方面：一是载波相位测量中引入与时间相关的相位旋转，若未补偿会直接污染相位式测距；二是码相位测量中，伪码的码率被“压缩”或“拉伸”，导致码相位时延估计出现偏差。

处理方式通常是在接收端利用已知的星历与相对速度先验，计算并剥离多普勒频移，再对残余频偏进行跟踪。对于双向时间比对（如 TWSTFT 或星间双向测距），多普勒在两条相反方向上符号相反，可在差分中部分抵消，但几何不对称时仍会残留。

四、大气与电离层：介质附加时延

星间链路若穿过地球大气层或电离层，信号传播速度不再是 c，而是 v=c/n，其中 n 为折射指数。附加时延可写为 Δτ=∫(n−1)/c dl。对电离层，n−1≈−40.3·N_e/f²，其中 N_e 为电子密度（el/m³），f 为频率。以 f=26 GHz、总电子含量 TEC=50 TECU（1 TECU=10¹⁶ el/m²）为例，电离层附加时延约为 0.3 ns；若频率降至 2 GHz，同样 TEC 下时延可达 50 ns 量级。可见高频链路对电离层不敏感，但低频链路必须修正。

对流层（中性大气）的附加时延主要来自干分量与湿分量，天顶方向约 2.3 m（干）与 0～0.4 m（湿），对应时延约 7.7 ns 与 0～1.3 ns。星间链路若视线高度角较低，需按映射函数放大。工程上常用模型（如 Saastamoinen、GPT2w）结合地面气象参数或数值天气模型进行修正。

五、建模流程与示例参数

把上述四项串起来，一个可执行的建模流程如下：第一步，读取两颗卫星在 t_tx 时刻的惯性系位置与速度；第二步，计算 d 与 τ_geo，并迭代修正到 t_rx；第三步，计算 v_los 与 dτ/dt，得到一阶运动修正；第四步，按 f_d=−(v_los/c)·f_c 计算多普勒，用于载波与码相位补偿；第五步，判断链路是否穿越电离层或对流层，若穿越则按 TEC 或气象模型计算 Δτ_media；第六步，将四项相加得到总时延，并输出不确定度。

下表给出一组教学示例参数（非实测，仅用于说明量级与计算方式）：

最后，建模精度取决于星历精度、时间系统统一程度以及介质模型适用性。对于亚纳秒级星间时间同步，几何项与运动项通常占主导，多普勒项需在载波环中实时跟踪，介质项则按链路几何与频段决定是否保留。把四类贡献分开建模、分别验证，再叠加成总时延，是工程上最稳妥的做法。