6.2 环球对时与卫星校时：原子钟实验与 GPS

6.2 环球对时与卫星校时：原子钟实验与 GPS

本站摘要：1971 年原子钟随民航机环球飞行，向东飞回的钟比地面慢了 59 纳秒、向西飞的快了 273 纳秒，与相对论预言吻合。今天的 GPS 卫星钟每天比地面快 38 微秒，不修正则定位每天漂移十公里量级。本节把这两张工程账单逐行算清。

6.1 验完μ子的票，本站把时间膨胀搬进人造钟表与导航系统。这是全册离日常生活最近的一节——你手机里的定位就在用这些数字。

环球对时

**背景。**1971 年，黑菲与基廷把四台铯原子钟装上商用客机，先向东环球、再向西环球，出发前与美国海军天文台的参考钟对齐。目的只有一个：看飞过的钟回到地面后读数差多少纳秒。

**操作与账目结构。**飞行钟的账单有两栏。狭义栏：钟相对惯性系运动，按γ变慢；地球在自转，东飞与西飞的"总速度"不同——东飞顺着自转，总速度更大、钟更慢；西飞逆着自转，总速度小到低于地面钟，飞钟反而快。引力栏（广义相对论的贡献，6.4 预告）：飞机高、引力弱，钟走得略快。两栏代数和才是预言。

# 演算 1：狭义栏的量级估算（取中纬度典型值做常数近似）
import math
c = 299_792_458

v_plane = 250.0                                 # 相对地面的typical航速 m/s
v_rot   = 463.0 * math.cos(math.radians(40))    # 北纬40度的自转线速度
T_east, T_west = 41.2*3600, 48.6*3600           # 东/西向累计飞行时间

v_east  = v_rot + v_plane                       # 东飞：顺自转，总速度更大
v_west  = abs(v_rot - v_plane)                  # 西飞：逆自转，总速度小于地面钟
dt_east = -((v_east**2 - v_rot**2) / (2*c**2)) * T_east
dt_west = -((v_west**2 - v_rot**2) / (2*c**2)) * T_west
print(f"自转线速度 = {v_rot:.0f} m/s, 东飞总速度 = {v_east:.0f}, 西飞总速度 = {v_west:.0f}")
# 输出: 自转线速度 = 355 m/s, 东飞总速度 = 605, 西飞总速度 = 105
print(f"狭义栏：东飞 ≈ {dt_east*1e9:.0f} 纳秒, 西飞 ≈ {dt_west*1e9:+.0f} 纳秒")
# 输出: 东飞 ≈ -198 纳秒, 西飞 ≈ +113 纳秒
# 正式分析按逐段航迹积分得 -184 与 +96 纳秒——常数近似即可对上量级与符号。
​

**结果。**正式分析（含逐段航迹积分与引力栏）的预言与实测：

**解读。**东飞钟净亏、西飞钟净赚，方向与大小都落在预言区间内。妙处在于两栏符号相反：单看狭义栏或引力栏都不对账，必须两本账一起算——相对论不是"修正项美学"，是两栏都精确到纳秒的复式记账。变式：2014 年有人用单座氢原子钟在国内航班复测，2016 年国际空间站计划把这类实验搬到轨道——账目结构不变，γ更大。

卫星校时：GPS 的每日账单

GPS 卫星在两万公里高的轨道上以约 3.9 公里每秒飞行，星上原子钟的账单同样两栏：

# 演算 2：GPS 星钟的每日漂移账单
import math
c = 299_792_458
GM = 3.986004418e14            # 地球引力参数（用于引力栏，公式来自广义相对论）
R_E, r_orbit = 6.371e6, 2.66e7 # 地球半径与卫星轨道半径（20200 km 高度）
v_sat = math.sqrt(GM / r_orbit) # 圆轨道速度
day = 86400.0

# 狭义栏：星钟运动快，每天慢 v²/2c² × 1天
sr = -(v_sat**2 / (2*c**2)) * day
# 引力栏：轨道处引力势更高，每天快 GM/c² × (1/R_E − 1/r) × 1天
gr = (GM/c**2) * (1/R_E - 1/r_orbit) * day
net = sr + gr
print(f"轨道速度 = {v_sat/1e3:.3f} km/s")
# 输出: 3.871 km/s
print(f"狭义栏（钟变慢）= {sr*1e6:+.1f} 微秒/天")   # 输出: -7.2 微秒/天
print(f"引力栏（钟变快）= {gr*1e6:+.1f} 微秒/天")    # 输出: +45.7 微秒/天
print(f"净值 = {net*1e6:+.1f} 微秒/天")              # 输出: +38.5 微秒/天
print(f"不修正的定位漂移 = {abs(net)*c/1e3:.1f} 公里/天")  # 输出: 11.6 公里/天
print(f"折合每分钟 = {net*1e9/1440:.1f} 纳秒")        # 输出: 26.8 纳秒
​

账单要点：引力栏是狭义栏的六倍多，两栏相抵后净剩每天快 38 微秒。38 微秒听着少，乘上光速就是每天 11.6 公里的定位漂移——导航系统出厂前必须预调星钟频率（把 10.23 兆赫的钟调低约万分之四）并让地面持续监控。相对论在这里不是学术装饰，是系统设计的输入参数。

卫星对时还有一个专业细节：环球与星地钟差计算用的是"同向约定"的对钟（光信号单向校准），而地球在自转，绕圈回来的信号会累积一个与回路面积成正比的时差（萨格纳克效应）。环球飞行实验的数据处理必须扣掉这项（量级上百纳秒），GPS 的地面站网也在实时修正它。这不是相对论出错，而是"旋转平台上的几何"必须单独记账——它提醒我们：本册的换算表建立在惯性系上，旋转的地球是"足够接近惯性"的近似，账房先生要心里有数。

图 GPS 星钟的两栏账

⚠️ 两个常见追问。其一："星钟快是广义相对论，不是你这本教程的功劳？"——狭义栏 −7.2 微秒完全是本册内容（v 已知就能算），引力栏只借了广义相对论的一个公式，6.4 会给出接口；没有狭义栏，账根本不平。其二："卫星速度 3.87 km/s 算出来 γ 才 1 加 8 乘 10 的负 11 次方，至于吗？"——单看比例确实微小，乘以一天 86400 秒再乘光速，公里级的漂移就是从这个"微小"里长出来的。相对论工程的全部戏剧性都在于：尾数 × 大基数 = 整数。

💡 一个反直觉的收获：GPS 反过来也成了相对论的永久实验台——星地钟差数据持续被学界拿来检验洛伦兹对称性，精度一年比一年高。一个理论最好的应用，往往是它最严格的考场。

本节要点回顾

• 环球实验：东飞实测 −59 纳秒、西飞 +273 纳秒，与两栏合计的预言吻合
• 账目结构：狭义栏看相对惯性系的速度（自转参与加减），引力栏看高度
• 东慢西快：东飞顺自转总速度大、钟更慢；西飞总速度小、钟反超地面
• GPS 每日账：狭义 −7.2、引力 +45.7、净 +38.5 微秒，漂移 11.6 公里每天
• 工程动作：星钟出厂预调频率约万分之四，地面持续监控
• 方法论：尾数乘大基数——相对论工程的一切"神奇"都是这个算式

钟的票验完，下一站验质量的票：太阳每秒烧掉四百万吨，账单就写在 5.2 的公式里。