ASHBY N. Relativity in the Global Positioning SystemJ/OL. Living Reviews in Relativity, 2003, 6: 1. 2026-07-19. https://www.livingreviews.org/lrr-2003-1.
全球定位系统中的相对论效应
《相对论综述》，2003年第6卷第1篇
作者：尼尔·阿什比
美国科罗拉多大学博尔德分校物理系，邮编80309--0390
邮箱：Neil.Ashby@colorado.edu
个人主页：http://www.colorado.edu/physics/Web/directory/faculty/ashby_n.html
收稿日期：2003年1月8日
发表日期：2003年1月28日
摘要
全球定位系统（GPS）依托卫星与地面高精度原子钟实现全球定位与授时。这类时钟存在引力与运动带来的频率偏移，若不精细处理各类相对论效应，整套系统将无法正常工作。本文基于狭义、广义相对论阐述卫星导航理论基础，梳理实际工程中必须纳入的相对论原理与效应：光速不变原理、等效原理、萨格纳克效应、时间膨胀、引力频移、同时性的相对性。本文还将基于TOPEX/POSEIDON卫星搭载GPS接收机开展的实验，介绍相对论相关实测验证；同时解释卫星轨道调整引发的频率跳变这类新近发现的相对论现象，并介绍GPS各类应用。
本文采用德国3.0版知识共享署名-非商业性使用-禁止演绎协议授权。
出版说明与使用条款
《相对论综述》由马克斯·普朗克引力物理研究所出版，是同行评审开放获取期刊，地址：德国波茨坦米伦贝格大街1号，邮编14476，国际标准连续出版物号ISSN 1433-8351。
本文采用德国3.0版知识共享署名-非商业性使用-禁止演绎协议授权：http://creativecommons.org/licenses/by-nc-nd/3.0/de/
由于《相对论综述》论文会持续更新，推荐引用格式如下：
Neil Ashby，"全球定位系统中的相对论效应"，《相对论综述》，2003，6，1. 在线文献：引用日期<日期>，http://www.livingreviews.org/lrr-2003-1
标注引用日期可唯一区分你阅读的论文版本。
论文修订记录
《相对论综述》提供两种更新机制：
快速修订：作者可补充短期研究成果、行业趋势或重要文献短评，由责任编辑同行评审；若本文完成快速修订，此处会列出修改摘要。
重大更新：包含大量新增、改写内容，需完整外审，发布时赋予新文献编号。
如需查阅完整版本更迭记录，请访问论文在线历史页面：http://www.livingreviews.org/lrr-2003-1
2007年6月21日修订说明：文末新增增强系统、全球导航系统两节；多处文本通用更新（删除"近期"等过时表述）；结论段落重新措辞；新增参考文献25；排版适配期刊新版式；第13页末项符号修正；第17页补充语句；第21页增加图4引用；第26页WGS-84(837)替换旧编号；第37、38页新增独立小节。
目录
1 引言 第5页
2 参考坐标系与萨格纳克效应 第7页
3 GPS坐标时与国际原子时（TAI） 第10页
4 坐标时的工程实现 第13页
5 卫星时钟的相对论效应 第15页
6 TOPEX/POSEIDON卫星相对论验证实验 第19页
7 多普勒效应 第24页
8 星间测距 第25页
9 轨道调整引发时钟频移 第26页
10 次级相对论效应 第35页
11 卫星导航增强系统 第37页
12 全球卫星导航系统 第38页
13 应用场景 第39页
14 结论 第40页
参考文献 第41页
文章目录
-
- 全球定位系统中的相对论效应
- [1 引言](#1 引言)
- [2 参考坐标系与萨格纳克效应](#2 参考坐标系与萨格纳克效应)
- [3 GPS坐标时与国际原子时TAI](#3 GPS坐标时与国际原子时TAI)
- [4 坐标时的工程实现](#4 坐标时的工程实现)
- [5 卫星时钟的相对论效应](#5 卫星时钟的相对论效应)
- [6 TOPEX/POSEIDON卫星相对论验证实验](#6 TOPEX/POSEIDON卫星相对论验证实验)
- [7 多普勒效应](#7 多普勒效应)
- [8 星间测距](#8 星间测距)
- [9 轨道调整引发时钟频移](#9 轨道调整引发时钟频移)
- [10 次级相对论效应](#10 次级相对论效应)
-
- [1. 夏皮罗传播延迟](#1. 夏皮罗传播延迟)
- [2. 空间曲率测地距离修正](#2. 空间曲率测地距离修正)
- [3 相位缠绕](#3 相位缠绕)
- [4 太阳系天体潮汐引力红移](#4 太阳系天体潮汐引力红移)
- [11 卫星导航增强系统](#11 卫星导航增强系统)
- [12 全球卫星导航系统（GNSS）](#12 全球卫星导航系统（GNSS）)
- [13 应用场景](#13 应用场景)
- [14 结论](#14 结论)
- 参考文献
1 引言
全球定位系统（GPS）可划分为三大核心分系统：空间段、控制段、用户段。
空间段核心为24颗搭载原子钟的卫星（含备用卫星与星载备用时钟）；6个轨道平面，每个平面布设4颗卫星，轨道倾角相对地球赤道面55°。轨道排布保证地表任意位置几乎随时能观测到4颗及以上卫星。每颗卫星向外广播绑定计时信号，这类信号可视为时空事件序列，携带发射位置、发射时刻信息；导航电文同时下发卫星星历、卫星健康状态，以及美国海军天文台维持的协调世界时UTC(USNO)解算参数。
控制段由多处地面监测站组成，持续采集卫星观测数据，传输至科罗拉多斯普林斯主控站 ；主控分析全星座轨道与时钟预报，将校正参数上注卫星，由卫星转发给地面用户。
用户段为所有接收卫星信号、解算自身位置、速度与本地时钟时刻的设备。
GPS由美国国防部运营，存在涉密技术细节 ；多家机构独立监测GPS信号，可获取卫星星历与钟差数据 。定位精度常规可达5--10厘米 ，载波相位测量精度普遍优于1毫米。
GPS地面接收两路载波信号：L1（154×10.23 MHz）、L2（120×10.23 MHz）。L1载波调制两类伪随机噪声码(pseudorandom noise codes)：1.023 MHz粗捕获码（C/A码）、10.23 MHz加密精码（P码）。精码接收机可同时接收两路载波，完成电离层延迟校正；民用用户仅能使用C/A码，对应两种定位服务：精定位服务（P码）、标准定位服务（C/A码）。美国国防部曾具备选择可用性（SA）干扰机制，人为抖动信号频率，将民用C/A码定位误差限制在±100米；克林顿总统于2000年5月下令永久关闭SA。
他国自研卫星导航体系加速成型（俄罗斯 GLONASS 逐步完善），就算美国 GPS 被干扰，对手可直接用俄导航，限制作用大幅削弱；
现代武器低成本惯性导航、光电制导、地形匹配完全可以脱离卫星定位，SA 拦不住军用打击手段；
美军自身作战也高度依赖民用 GPS 配套基础设施（民航、后勤运输），SA 全球干扰反而拖累己方后勤。
GPS的技术根基是超高稳定原子钟。图1给出高性能铯 钟的艾伦偏差随平均采样时间 τ \tau τ变化曲线。若一组时钟初始同步，经过时长 T T T比对后，艾伦偏差可表征时钟固有噪声带来的相对频率均方根偏差 ；频偏、频漂属于额外系统误差 ，需单独修正。图中同时给出GPS接收机常用石英振荡器艾伦偏差曲线：石英短期稳定性优于 铯钟，但100秒后长期稳定性能远弱于铯钟；实际时钟艾伦偏差会在标称曲线附近浮动。
图1 铯原子钟与石英振荡器典型艾伦偏差随平均时长 T T T变化曲线
铯钟曲线代表GPS在轨最优时钟性能：一台铯钟初始化静置一日后，走时误差约 5 × 10 − 14 5×10^{-14} 5×10−14量级，对应4纳秒 ；而相对论效应带来的时误差远大于该数值。
本文旨在解释GPS系统中各类相对论效应的校正方法。尽管近地卫星速度、地球引力场均为弱场近似 ，但相对论修正量不可忽略；相关效应包含相对运动一阶、二阶多普勒频移、引力频移、地球自转引发萨格纳克效应。若未正确修正，定位、授时将产生无法接受的巨大误差。GPS是狭义、广义相对论基础原理大量落地的工程实例，具备极高教学研究价值；同时可依托GPS开展相对论实验，只是测量精度一般不超越传统引力测试手段。
GPS定位与授时原理可简洁表述：设4台同步原子钟 分别在时空点 r j \boldsymbol{r}_j rj、时刻 t j t_j tj发射窄脉冲（下标 j = 1 , 2 , 3 , 4 j=1,2,3,4 j=1,2,3,4区分4颗卫星）；假设地面接收机在同一瞬时 t t t接收四路信号 ，根据光速不变原理可得方程组：
c 2 ( t − t j ) 2 = ∣ r − r j ∣ 2 , j = 1 , 2 , 3 , 4 (1) c^{2}(t-t_{j})^{2}=|\boldsymbol{r}-\boldsymbol{r}_{j}|^{2},\quad j=1,2,3,4 \tag{1} c2(t−tj)2=∣r−rj∣2,j=1,2,3,4(1)
其中真空光速定义值 c = 299792458 m/s c=299792458\ \text{m/s} c=299792458 m/s。联立四式可解接收事件四维时空坐标 ( r , t ) (\boldsymbol{r},t) (r,t)，光速不变原理是GPS底层核心理论。1纳秒计时误差将造成约30厘米定位偏差；同时必须严格统一卫星时钟同步所用参考系，否则式(1)失效。
接收机三维坐标(x,y,z) + 接收机本地钟差δt，共 4 个未知量，至少需要 4 颗卫星列 4 个独立方程才能唯一求解；
3 颗卫星只能解空间位置、消不掉时钟带来的测距偏差，这个关键点讲得很清楚。
计时脉冲可理解为圆偏振电磁波波列内特定相位翻转位置（不依赖于观察者惯性系的通用时间标记），电磁场张量在此处过零，理论上可供相对运动观测者统一识别时序事件（同时是具有相对性的，但是一个事件是客观发生的，"事件本身"------即在时空中某一点发生了"场张量过零"这一物理事实------是绝对的，不依赖于任何参考系。所有观测者都必须承认这个事件发生了。）
2 参考坐标系与萨格纳克效应
绝大多数GPS用户固定在旋转地球表面，或地表低速移动。因此早期系统设计决定，卫星星历 采用地心地固坐标系（ECEF，WGS-84(G873)）广播；该坐标系模型地球以固定角速度 ω E = 7.2921151467 × 10 − 5 rad/s \omega_{E}=7.2921151467 ×10^{-5}\ \text{rad/s} ωE=7.2921151467×10−5 rad/s绕固定转轴自转。相对论分析时，ECEF具体版本无本质影响；真实地球自转轴、自转角速度存在微小变化，本文不展开讨论，仅将GPS的ECEF视为巴黎国际计量局（BIPM）定义的国际地面参考系衍生坐标系。
ECEF坐标系本身随着地球一起旋转，原点在地心，z轴指向北极，x轴指向本初子午线与赤道的交点，y轴与之垂直。
需重点说明：导航电文仅提供发射时刻对应的卫星 ECEF位置函数；各卫星发射时刻 t j t_j tj互不相等，每次观测对应的坐标系空间指向均不同。为代入式(1)计算，接收机必须将每一组卫星观测数据分别旋转至同一惯性系求解传播延迟；解算完成后再旋转回ECEF得到地面位置。该流程逻辑极其复杂，技术文档10有完整详细推导。
虽然导航以ECEF为主，但电磁波传播等物理过程在惯性系中描述更简洁；式(1)仅能在惯性系成立，直接在地固系使用会出现严重误差。行业常讨论的标准局部惯性系原点与地心重合，随地球在太阳系引力场自由下落， z z z轴对齐J2000历元地球角动量轴；该局部惯性系可通过坐标变换关联国际天球参考系ICRF（以约500颗射电星坐标定义的太阳系质心惯性系）。
物理定律（如光速不变）仅在惯性系中具有最简洁的形式。因此GPS在计算时，不会直接在旋转的ECEF坐标系中构建方程，而是先将卫星轨道参数在ICRF惯性系中解算，再将结果通过坐标变换转换到ECEF，并修正地球自转带来的Sagnac效应。
历元 = 固定标准时间戳：J2000.0 = 2000 年 1 月 1 日 12:00 地球力学时（TT）
航天、GPS、卫星轨道计算统一锁定这个时刻做静态基准，因为地球自转轴、春分点一直在缓慢漂移（岁差、章动），必须固定一个瞬间当标准，否则坐标系会持续变形。
地固 WGS84 坐标系 Z 轴，随极移持续变动，不能用于惯性轨道 / 相对论计算。
GPS地固坐标系的时间单位为国际单位制秒（美国海军原子钟组复现），长度单位为国际单位制米；这一设计令GPS局部观测不受太阳系天体随时间变化的尺度效应干扰。（米、秒基准不依赖太阳 / 行星）
先忽略引力势，推导惯性系至旋转坐标系的时空变换。惯性系柱坐标下 度规（闵氏时空）：
− d s 2 = − ( c d t ) 2 + d r 2 + r 2 d ϕ 2 + d z 2 (2) -d s^{2}=-(c d t)^{2}+d r^{2}+r^{2} d \phi^{2}+d z^{2} \tag{2} −ds2=−(cdt)2+dr2+r2dϕ2+dz2(2)
以恒定角速度 ω E \omega_E ωE旋转的坐标系 ( t ′ , r ′ , ϕ ′ , z ′ ) (t',r',\phi',z') (t′,r′,ϕ′,z′)变换关系：
t = t ′ , r = r ′ , ϕ = ϕ ′ + ω E t ′ , z = z ′ (3) t=t',\ r=r',\ \phi=\phi'+\omega_{E} t',\ z=z' \tag{3} t=t′, r=r′, ϕ=ϕ′+ωEt′, z=z′(3)
代入得到朗之万旋转度规 ：
− d s 2 = − ( 1 − ω E 2 r ′ 2 c 2 ) ( c d t ′ ) 2 + 2 ω E r ′ 2 d ϕ ′ d t ′ + ( d σ ′ ) 2 (4) -d s^{2}=-\left(1-\frac{\omega_{E}^{2} r^{\prime 2}}{c^{2}}\right)\left(c d t'\right)^{2}+2 \omega_{E} r^{\prime 2} d \phi' d t'+\left(d \sigma'\right)^{2} \tag{4} −ds2=−(1−c2ωE2r′2)(cdt′)2+2ωEr′2dϕ′dt′+(dσ′)2(4)
其中坐标距离元简写 ( d σ ′ ) 2 = ( d r ′ ) 2 + ( r ′ d ϕ ′ ) 2 + ( d z ′ ) 2 (d \sigma')^{2}=(d r')^{2}+(r' d \phi')^{2}+(d z')^{2} (dσ′)2=(dr′)2+(r′dϕ′)2+(dz′)2。
式(3)中 t = t ′ t=t' t=t′看似简单，实则旋转系的时间变量 t ′ t' t′底层依托惯性系坐标时定义，GPS采用完全相同的设计思路。
现假设地面观测者在旋转系 内，试图用爱因斯坦同步法（光速不变原理）搭建全局同步时钟网络。光沿类光世界线传播，令 d s 2 = 0 d s^{2}=0 ds2=0；仅保留小量 ω E r ′ / c \omega_{E} r'/c ωEr′/c一阶项化简：
( c d t ′ ) 2 − 2 ω E r ′ 2 d ϕ ′ ( c d t ′ ) c − ( d σ ′ ) 2 = 0 (5) (c d t')^{2}-\frac{2 \omega_{E} r^{\prime 2} d \phi'\left(c d t'\right)}{c}-\left(d \sigma'\right)^{2}=0 \tag{5} (cdt′)2−c2ωEr′2dϕ′(cdt′)−(dσ′)2=0(5)
整理解出 c d t ′ c d t' cdt′：
c d t ′ = d σ ′ + ω E r ′ 2 d ϕ ′ c (6) c d t'=d \sigma'+\frac{\omega_{E} r^{\prime 2} d \phi'}{c} \tag{6} cdt′=dσ′+cωEr′2dϕ′(6)
式中 r ′ 2 d ϕ ′ / 2 r'^{2} d \phi'/2 r′2dϕ′/2为旋转系内矢量自转轴扫出、投影至赤道平面无穷小面积元 d A z ′ dA_z' dAz′。光沿任意路径传播的总坐标时积分：
∫ p a t h d t ′ = ∫ p a t h d σ ′ c + 2 ω E c 2 ∫ p a t h d A z ′ 光传播 (7) \int_{path} d t'=\int_{path} \frac{d \sigma'}{c}+\frac{2 \omega_{E}}{c^{2}} \int_{path} d A_{z}' \quad 光传播 \tag{7} ∫pathdt′=∫pathcdσ′+c22ωE∫pathdAz′光传播(7)
若地面观测者忽略地球自转，仅用 ∫ d σ ′ / c \int d \sigma'/c ∫dσ′/c同步时钟，会产生与光路包围投影面积成正比 的自洽偏差。举例：沿赤道向东同步一圈，地球参数 2 ω E / c 2 = 1.6227 × 10 − 21 s/m 2 2 \omega_{E}/c^{2}=1.6227 ×10^{-21}\ \text{s/m}^2 2ωE/c2=1.6227×10−21 s/m2，赤道半径 a 1 = 6378137 m a_1=6378137\ \text{m} a1=6378137 m，赤道包围面积 π a 1 2 = 1.27802 × 10 14 m 2 \pi a_1^{2}=1.27802 ×10^{14}\ \text{m}^2 πa12=1.27802×1014 m2，代入第二项得总修正：
2 ω E c 2 ∫ p a t h d A z ′ = 207.4 ns (8) \frac{2 \omega_{E}}{c^{2}} \int_{path} d A_{z}'=207.4\ \text{ns} \tag{8} c22ωE∫pathdAz′=207.4 ns(8)
从惯性系视角解释：地球向东自转，光线需要额外传播距离追上移动参考点 ；仅在地固系套用光速同步会产生巨大误差。向东绕行一圈同步链路末端时钟滞后首钟207.4纳秒，向西绕行则超前207.4纳秒。
惯性系中可携带慢动钟实现时间传递，移动钟时间膨胀效应可忽略 ；但旋转系内移动标准钟的固有时会受地球自转影响。对式(4)因式分解，运动钟固有时增量 d τ d\tau dτ：
( d τ ) 2 = ( d s / c ) 2 = d t ′ 2 1 − ( ω E r ′ c ) 2 − 2 ω E r ′ 2 d ϕ ′ c 2 d t ′ − ( d σ ′ c d t ′ ) 2 (9) (d \tau)^{2}=(d s/c)^{2}=d t^{\prime 2}\left1-\\left(\\frac{\\omega_{E} r\^{\\prime}}{c}\\right)\^{2}-\\frac{2 \\omega_{E} r\^{\\prime 2} d \\phi'}{c\^{2} d t'}-\\left(\\frac{d \\sigma'}{c d t'}\\right)\^{2}\\right \tag{9} (dτ)2=(ds/c)2=dt′21−(cωEr′)2−c2dt′2ωEr′2dϕ′−(cdt′dσ′)2(9)
低速下 ( d σ ′ / c d t ′ ) 2 ≪ 1 (d \sigma'/c d t')^{2}\ll1 (dσ′/cdt′)2≪1可舍去，仅保留 ω E r ′ / c \omega_{E} r'/c ωEr′/c一阶项化简：
d τ = d t ′ − ω E r ′ 2 d ϕ ′ c 2 (10) d \tau=d t'-\frac{\omega_{E} r^{\prime 2} d \phi'}{c^{2}} \tag{10} dτ=dt′−c2ωEr′2dϕ′(10)
积分得到移动钟对应的坐标时关系：
∫ p a t h d t ′ = ∫ p a t h d τ + 2 ω E c 2 ∫ p a t h d A z ′ 便携时钟授时 (11) \int_{path} d t'=\int_{path} d \tau+\frac{2 \omega_{E}}{c^{2}} \int_{path} d A_{z}' \quad 便携时钟授时 \tag{11} ∫pathdt′=∫pathdτ+c22ωE∫pathdAz′便携时钟授时(11)
对比式(7)：无论用光同步还是移动原子钟，旋转系必然产生路径相关时差；而底层惯性系内同步自洽。式(7)(11)说明：完成同步后叠加 + 2 ω E / c 2 ∫ d A z ′ +2 \omega_{E}/c^{2} \int dA_z' +2ωE/c2∫dAz′修正，即可在地固系复现统一坐标时。萨格纳克sagnac效应存在多重物理解释：可视为向旋转地表局部惯性系洛伦兹变换带来同时性相对性；也可理解为便携钟固有时与地面基准钟的差值。
1980年代初，秒定义咨询委员会、国际无线电咨询委员会正式采纳该修正用于全球时间标准比对；对应GPS设计：所有地面、在轨原子钟统一在地球中心局部惯性系（ECI）内同步6。
单颗卫星同时覆盖两处地面基准站时，可采用共视法比对两地原子钟；萨格纳克修正可达数百纳秒，是比对核心项。1984年利用3、4、6、8号GPS卫星同步开展环球萨格纳克闭环实验：实验站点为美国国家标准局（NBS）、德国布伦瑞克物理技术研究院（PTB）、东京天文台（TAO），修正量区间240--350纳秒，完成90组独立环球测量；三组时差相加残差均值仅5纳秒，不足理论萨格纳克效应的2%4。
3 GPS坐标时与国际原子时TAI
GPS的坐标时 t ′ = t t'=t t′=t在地固旋转系内落地，全部同步操作依托自洽惯性系完成。本节引入地球自身引力场模型，假设地球质量静态，建立原点在地心、随地球自由下落无旋转局部惯性系，写出各向同性坐标下爱因斯坦场方程近似解：
− d s 2 = − ( 1 + 2 V c 2 ) ( c d t ) 2 + ( 1 − 2 V c 2 ) ( d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2 ) (12) -d s^{2}=-\left(1+\frac{2 V}{c^{2}}\right)(c d t)^{2}+\left(1-\frac{2 V}{c^{2}}\right)\left(d r^{2}+r^{2} d \theta^{2}+r^{2} \sin ^{2} \theta d \phi^{2}\right) \tag{12} −ds2=−(1+c22V)(cdt)2+(1−c22V)(dr2+r2dθ2+r2sin2θdϕ2)(12)
式中 r , θ , ϕ r,\theta,\phi r,θ,ϕ为球极坐标， V V V为牛顿地球引力势近似：
V = − G M E r 1 − J 2 ( a 1 r ) 2 P 2 ( cos θ ) (13) V=-\frac{G M_{E}}{r}\left1-J_{2}\\left(\\frac{a_{1}}{r}\\right)\^{2} P_{2}(\\cos \\theta)\\right \tag{13} V=−rGME1−J2(ra1)2P2(cosθ)(13)
参数说明： G M E = 3.986004418 × 10 14 m 3 / s 2 G M_{E}=3.986004418 ×10^{14}\ \text{m}^3/\text{s}^2 GME=3.986004418×1014 m3/s2（地球引力常数乘积）； J 2 = 1.0826300 × 10 − 3 J_{2}=1.0826300 ×10^{-3} J2=1.0826300×10−3地球二阶带谐系数； a 1 = 6.3781370 × 10 6 m a_{1}=6.3781370 ×10^{6}\ \text{m} a1=6.3781370×106 m地球赤道半径； θ \theta θ为极角（自旋转对称轴向下计量）； P 2 P_2 P2为二阶勒让德多项式。式(12)仅保留 V / c 2 V/c^2 V/c2一阶项，更高阶多极矩对GPS相对论计算影响极小。
将式(12)通过式(3)类变换转为地固旋转球坐标，仅保留 1 / c 2 1/c^2 1/c2量级项：
− d s 2 = − 1 + 2 V c 2 − ( ω E r ′ sin θ ′ c ) 2 ( c d t ′ ) 2 + 2 ω E r ′ 2 sin 2 θ ′ d ϕ ′ d t ′ + ( 1 − 2 V c 2 ) ( d r ′ 2 + r ′ 2 d θ ′ 2 + r ′ 2 sin 2 θ ′ d ϕ ′ 2 ) (15) \begin{aligned} -d s^{2}= & -\left1+\\frac{2 V}{c\^{2}}-\\left(\\frac{\\omega_{E} r' \\sin \\theta'}{c}\\right)\^{2}\\right\left(c d t'\right)^{2}+2 \omega_{E} r^{\prime 2} \sin ^{2} \theta' d \phi' d t' \\ & +\left(1-\frac{2 V}{c^{2}}\right)\left(d r^{\prime 2}+r^{\prime 2} d \theta^{\prime 2}+r^{\prime 2} \sin ^{2} \theta' d \phi^{\prime 2}\right) \end{aligned} \tag{15} −ds2=−1+c22V−(cωEr′sinθ′)2(cdt′)2+2ωEr′2sin2θ′dϕ′dt′+(1−c22V)(dr′2+r′2dθ′2+r′2sin2θ′dϕ′2)(15)
旋转系度规00分量定义等效引力势 Φ \Phi Φ：
g 00 ′ = − 1 + 2 V c 2 − ( ω E r ′ sin θ ′ c ) 2 ≡ − ( 1 + 2 Φ c 2 ) (16) g_{00}'=-\left1+\\frac{2 V}{c\^{2}}-\\left(\\frac{\\omega_{E} r' \\sin \\theta'}{c}\\right)\^{2}\\right \equiv-\left(1+\frac{2\Phi}{c^{2}}\right) \tag{16} g00′=−1+c22V−(cωEr′sinθ′)2≡−(1+c22Φ)(16)
Φ \Phi Φ包含地球静态引力势与旋转离心势两项叠加。
大地水准面
式(12)(15)中坐标时速率以无穷远静止时钟为基准；GPS坐标时紧密关联国际原子时TAI。TAI由巴黎国际计量局基于全球数百台一级原子钟、氢钟加权计算生成，全部原始固有时统一归算至地固系平均海平面等势面（大地水准面）得到统一时标。
协调世界时UTC与TAI相差整数跳秒；定期插入闰秒保证UTC贴合地球自转。各授时机构维护独立本地时标，美国海军天文台时标记为UTC(USNO)；GPS时与UTC(USNO)绑定，除闰秒差值外偏差控制在100纳秒内，日常差值小于40纳秒。GPS设备无法承受闰秒跳变，会造成信号失锁、瞬态误差。
本节推导大地水准面静止原子钟走时速率；地球椭球令地面钟离地心距离不同，同时地球自转带来运动时差。使用ECEF度规式(15)，地面静止时钟简化度规：
− d s 2 = − ( 1 + 2 V c 2 − ω E 2 r ′ 2 sin 2 θ ′ c 2 ) ( c d t ′ ) 2 (17) -d s^{2}=-\left(1+\frac{2 V}{c^{2}}-\frac{\omega_{E}^{2} r^{\prime 2} \sin ^{2} \theta'}{c^{2}}\right)\left(c d t'\right)^{2} \tag{17} −ds2=−(1+c22V−c2ωE2r′2sin2θ′)(cdt′)2(17)
代入赤道处 θ ′ = π / 2 , r ′ = a 1 \theta'=\pi/2,\ r'=a_1 θ′=π/2, r′=a1求解大地水准面等效势常数 Φ 0 / c 2 \Phi_0/c^2 Φ0/c2：
Φ 0 c 2 = − G M E a 1 c 2 − G M E J 2 2 a 1 c 2 − ω E 2 a 1 2 2 c 2 = − 6.95348 × 10 − 10 − 3.764 × 10 − 13 − 1.203 × 10 − 12 = − 6.96927 × 10 − 10 (18) \begin{aligned} \frac{\Phi_{0}}{c^{2}} & =-\frac{G M_{E}}{a_{1} c^{2}}-\frac{G M_{E} J_{2}}{2 a_{1} c^{2}}-\frac{\omega_{E}^{2} a_{1}^{2}}{2 c^{2}} \\ & =-6.95348 × 10^{-10}-3.764 × 10^{-13}-1.203 × 10^{-12} \\ & =-6.96927 × 10^{-10} \tag{18} \end{aligned} c2Φ0=−a1c2GME−2a1c2GMEJ2−2c2ωE2a12=−6.95348×10−10−3.764×10−13−1.203×10−12=−6.96927×10−10(18)
等效势三部分贡献：地球质量单极引力项（主导）、二阶带谐椭球修正、旋转离心项；离心修正比主项小500倍，椭球修正小2000倍。国际天文学联合会最新决议定义标准值 Φ 0 / c 2 = 6.969290134 × 10 − 10 \Phi_0/c^2=6.969290134 ×10^{-10} Φ0/c2=6.969290134×10−10，本文计算值精度满足GPS需求。
大地水准面上静止钟固有时与坐标时关系：
d τ = d s / c = d t ′ ( 1 + Φ 0 c 2 ) (19) d \tau=d s/c=d t'\left(1+\frac{\Phi_{0}}{c^{2}}\right) \tag{19} dτ=ds/c=dt′(1+c2Φ0)(19)
结论：旋转大地水准面上原子钟相对无穷远钟慢约 7 × 10 − 10 7×10^{-10} 7×10−10，该效应比高性能铯钟本征稳定度大一万倍。令 Φ = Φ 0 \Phi=\Phi_0 Φ=Φ0可由式(16)反演大地水准面半径 r ′ r' r′关于极角 θ ′ \theta' θ′的幂级数展开：
r ′ = ( 6356742.025 + 21353.642 x ′ 2 + 39.832 x ′ 4 + 0.798 x ′ 6 + 0.003 x ′ 8 ) m , x ′ = sin θ ′ (20) r'=\left(6356742.025+21353.642\, x^{\prime 2}+39.832\, x^{\prime 4}+0.798\, x^{\prime 6}+0.003\, x^{\prime 8}\right)\ \text{m},\quad x'=\sin\theta' \tag{20} r′=(6356742.025+21353.642x′2+39.832x′4+0.798x′6+0.003x′8) m,x′=sinθ′(20)
该模型为参考椭球，真实大地水准面与模型差值最高百米，高阶多极矩对时钟计时影响微弱，仅改变 Φ 0 \Phi_0 Φ0第六位有效数字。
大地水准面为旋转系等势面，地表理想原子钟走时速率完全一致。高纬度站点离地心更近，引力红移更强，但自转线速度更低，二阶多普勒蓝移抵消引力效应，两者在参考椭球上严格平衡。
式(12)坐标时以无穷远钟为基准，现重新定义以大地水准面原子钟速率为基准的新坐标时 t ′ ′ t'' t′′：
t ′ ′ = ( 1 + Φ 0 / c 2 ) t ′ = ( 1 + Φ 0 / c 2 ) t (21) t''=\left(1+\Phi_{0} / c^{2}\right) t'=\left(1+\Phi_{0} / c^{2}\right) t \tag{21} t′′=(1+Φ0/c2)t′=(1+Φ0/c2)t(21)
修正系数约 7 × 10 − 10 7×10^{-10} 7×10−10。替换后地固系完整度规（仅保留 1 / c 2 1/c^2 1/c2项）：
− d s 2 = − 1 + 2 ( Φ − Φ 0 ) c 2 ( c d t ′ ′ ) 2 + 2 ω E r ′ 2 sin 2 θ ′ d ϕ ′ d t ′ ′ + ( 1 − 2 V c 2 ) ( d r ′ 2 + r ′ 2 d θ ′ 2 + r ′ 2 sin 2 θ ′ d ϕ ′ 2 ) (22) \begin{aligned} -d s^{2}= & -\left1+\\frac{2\\left(\\Phi-\\Phi_{0}\\right)}{c\^{2}}\\right\left(c d t''\right)^{2}+2 \omega_{E} r^{\prime 2} \sin ^{2} \theta' d \phi' d t'' \\ & +\left(1-\frac{2 V}{c^{2}}\right)\left(d r^{\prime 2}+r^{\prime 2} d \theta^{\prime 2}+r^{\prime 2} \sin ^{2} \theta' d \phi^{\prime 2}\right) \end{aligned} \tag{22} −ds2=−1+c22(Φ−Φ0)(cdt′′)2+2ωEr′2sin2θ′dϕ′dt′′+(1−c22V)(dr′2+r′2dθ′2+r′2sin2θ′dϕ′2)(22)
同理无旋转ECI惯性系度规更新为：
− d s 2 = − ( 1 + 2 ( V − Φ 0 ) c 2 ) ( c d t ) 2 + ( 1 − 2 V c 2 ) ( d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2 ) (23) -d s^{2}=-\left(1+\frac{2\left(V-\Phi_{0}\right)}{c^{2}}\right)(c d t)^{2}+\left(1-\frac{2 V}{c^{2}}\right)\left(d r^{2}+r^{2} d \theta^{2}+r^{2} \sin ^{2} \theta d \phi^{2}\right) \tag{23} −ds2=−(1+c22(V−Φ0))(cdt)2+(1−c22V)(dr2+r2dθ2+r2sin2θdϕ2)(23)
式(22)(23)表明大地水准面静止钟固有时等于新坐标时 t ′ ′ t'' t′′，完美贴合地面标准原子钟基准。
4 坐标时的工程实现
下文省略 t ′ ′ t'' t′′撇号，统一记为 t t t，该时标基准为大地水准面UTC(USNO)，同步基准为地心局部惯性系ECI，度规简化为式(24)：
− d s 2 = − ( 1 + 2 ( V − Φ 0 ) c 2 ) ( c d t ) 2 + ( 1 − 2 V c 2 ) ( d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2 ) (24) -d s^{2}=-\left(1+\frac{2\left(V-\Phi_{0}\right)}{c^{2}}\right)(c d t)^{2}+\left(1-\frac{2 V}{c^{2}}\right)\left(d r^{2}+r^{2} d \theta^{2}+r^{2} \sin ^{2} \theta d \phi^{2}\right) \tag{24} −ds2=−(1+c22(V−Φ0))(cdt)2+(1−c22V)(dr2+r2dθ2+r2sin2θdϕ2)(24)
式中 ( V − Φ 0 ) (V-\Phi_0) (V−Φ0)项源于ECI惯性系内引力势随空间变化，直接决定运动钟固有时与坐标时差值。广义相对论框架下式(24)在覆盖地球与GPS星座的坐标邻域内单值有效（规范自由度不唯一，但本文选用最简规范），适合作为全系统同步基准。
对式(24)变形，提取 ( c d t ) 2 (c d t)^2 (cdt)2公因子，定义ECI系卫星速度 v 2 = d r 2 + r 2 d θ 2 + r 2 sin 2 θ d t 2 d t 2 v^2=\frac{d r^{2}+r^{2} d \theta^{2}+r^{2} \sin ^{2} \theta d t^{2}}{d t^{2}} v2=dt2dr2+r2dθ2+r2sin2θdt2，仅保留 1 / c 2 1/c^2 1/c2一阶项化简，运动原子钟固有时增量：
d τ = d s / c = 1 + ( V − Φ 0 ) c 2 − v 2 2 c 2 d t (27) d \tau=d s / c=\left1+\\frac{\\left(V-\\Phi_{0}\\right)}{c\^{2}}-\\frac{v\^{2}}{2 c\^{2}}\\right d t \tag{27} dτ=ds/c=1+c2(V−Φ0)−2c2v2dt(27)
沿卫星轨道积分得到坐标时修正关系：
∫ p a t h d t = ∫ p a t h d τ 1 − ( V − Φ 0 ) c 2 + v 2 2 c 2 (28) \int_{path} d t=\int_{path} d \tau\left1-\\frac{\\left(V-\\Phi_{0}\\right)}{c\^{2}}+\\frac{v\^{2}}{2 c\^{2}}\\right \tag{28} ∫pathdt=∫pathdτ1−c2(V−Φ0)+2c2v2(28)
无引力场直观模型：构建无旋转地心局部惯性系，全域架设一组通过爱因斯坦同步的虚拟标准钟，读数为坐标 t t t；地表各类实体原子钟需结合位置、速度套用式(28)修正，使本地钟时刻与瞬时重合惯性虚拟钟时刻对齐。叠加地球引力势后，额外增加引力红移修正项，即式(28)第一项。
5 卫星时钟的相对论效应
分析卫星钟优先选用ECI地心惯性系，萨格纳克效应仅针对地面运动接收机，卫星层面无需考虑 。联立引力频移、二阶多普勒效应，依托式(28)推导卫星钟总修正； Φ 0 \Phi_0 Φ0包含大地水准面尺度修正，地球二阶带谐项贡献 -3.76 ×10\^{-13} 相对频偏， G P S 必须计入。地球二阶带谐引力势随距离快速衰减，此前常忽略，后文讨论其微小影响；暂时仅保留单极引力项 相对频偏，GPS必须计入。地球二阶带谐引力势随距离快速衰减，此前常忽略，后文讨论其微小影响；暂时仅保留单极引力项 相对频偏，GPS必须计入。地球二阶带谐引力势随距离快速衰减，此前常忽略，后文讨论其微小影响；暂时仅保留单极引力项V=-G M_E/r。
假设卫星遵循开普勒轨道，忽略高阶多极矩，牛顿轨道能量守恒可简化式(28)。 a a a轨道半长轴， e e e偏心率，EC系卫星地心距：
r = a ( 1 − e 2 ) 1 + e cos f (29) r=\frac{a\left(1-e^{2}\right)}{1+e \cos f} \tag{29} r=1+ecosfa(1−e2)(29)
f f f真近点角（自近地点沿轨道计量），真近点角与偏近点角 E E E转换关系：
{ cos f = cos E − e 1 − e cos E sin f = 1 − e 2 sin E 1 − e cos E (30) \begin{cases} \cos f=\frac{\cos E-e}{1-e \cos E} \\ \sin f=\sqrt{1-e^{2}} \frac{\sin E}{1-e \cos E} \end{cases} \tag{30} {cosf=1−ecosEcosE−esinf=1−e2 1−ecosEsinE(30)
地心距另一种表达：
r = a ( 1 − e cos E ) (31) r=a(1-e \cos E) \tag{31} r=a(1−ecosE)(31)
偏近点角超越方程（开普勒方程）：
E − e sin E = G M E a 3 ( t − t p ) (32) E-e \sin E=\sqrt{\frac{G M_{E}}{a^{3}}}\left(t-t_{p}\right) \tag{32} E−esinE=a3GME (t−tp)(32)
t p t_p tp为过近地点坐标时。单位质量轨道总机械能守恒：
1 2 v 2 − G M E r = − G M E 2 a (33) \frac{1}{2} v^{2}-\frac{G M_{E}}{r}=-\frac{G M_{E}}{2 a} \tag{33} 21v2−rGME=−2aGME(33)
将 v 2 v^2 v2代入式(28)，沿轨道积分卫星钟坐标时增量：
Δ t = ∫ p a t h d τ 1 + 3 G M E 2 a c 2 + Φ 0 c 2 − 2 G M E c 2 ( 1 a − 1 r ) (34) \Delta t=\int_{path} d\tau \left 1+\\frac{3GM_{E}}{2ac\^{2}}+\\frac{\\Phi_{0}}{c\^{2}}-\\frac{2GM_{E}}{c\^{2}}\\left( \\frac{1}{a}-\\frac{1}{r} \\right) \\right \tag{34} Δt=∫pathdτ1+2ac23GME+c2Φ0−c22GME(a1−r1)(34)
恒定频偏两项数值：
3 G M E 2 a c 2 + Φ 0 c 2 = + 2.5046 × 10 − 10 − 6.9693 × 10 − 10 = − 4.4647 × 10 − 10 (35) \frac{3 G M_{E}}{2 a c^{2}}+\frac{\Phi_{0}}{c^{2}}=+2.5046 × 10^{-10}-6.9693 × 10^{-10}=-4.4647 × 10^{-10} \tag{35} 2ac23GME+c2Φ0=+2.5046×10−10−6.9693×10−10=−4.4647×10−10(35)
负号代表在轨钟固有走时快于地面大地水准面钟，主导原因为引力蓝移。为让地面观测到基准10.23 MHz，卫星出厂预先压低固有频率：
1 − 4.4647 × 10 − 10 × 10.23 MHz = 10.2299999543 MHz (36) \left1-4.4647 × 10\^{-10}\\right × 10.23\ \text{MHz}=10.2299999543\ \text{MHz} \tag{36} 1−4.4647×10−10×10.23 MHz=10.2299999543 MHz(36)
该偏移在地面出厂调试完成。
图2 圆轨道卫星时钟总相对频移随轨道径向距离变化曲线
近地航天飞机：速度极高时间膨胀红移占主导；GPS、格洛纳斯、地球同步轨道卫星引力蓝移更强；两者抵消临界轨道半径约9545千米；伽利略轨道半径约30000千米。
历史背景：1977年6月23日NTS-2卫星搭载首台在轨铯钟，当时学界对相对论修正符号、量级存在争议，甚至有人质疑相对论效应必须修正。卫星时钟内置频率合成器，若在轨实测频偏符合广义相对论预测，可在轨切换补偿。卫星铯钟开机20天实测相对地面频偏 + 442.5 × 10 − 12 +442.5×10^{-12} +442.5×10−12，理论预测 + 446.5 × 10 − 12 +446.5×10^{-12} +446.5×10−12，误差在在轨钟精度范围内，在4.2倍地球半径轨道完成约1%精度相对论验证。
卫星钟额外频偏来源：钟漂、环境扰动、入轨半长轴偏差；导航电文下发钟差多项式补偿。早期GPS出厂频偏计算遗漏地球离心势项，1980年代初发现，时隔8年更新系统规范；铷钟发射易产生不可预测频跳，一般不预设固定出厂偏移，入轨实测频偏后写入导航电文。
轨道偏心周期性修正
对式(34)随时间项积分，偏近点角时间导数：
d E d t = G M E / a 3 1 − e cos E (37) \frac{d E}{d t}=\frac{\sqrt{G M_{E} / a^{3}}}{1-e \cos E} \tag{37} dtdE=1−ecosEGME/a3 (37)
低阶近似 d s / c ≈ d t ds/c \approx dt ds/c≈dt，积分周期项：
∫ 2 G M E c 2 ( 1 r − 1 a ) d s c ≃ 2 G M E c 2 ∫ ( 1 r − 1 a ) d t = 2 G M E a c 2 ∫ e cos E 1 − e cos E d t = 2 G M E a c 2 e ( sin E − sin E 0 ) = + 2 G M E a c 2 e sin E + 常数 (38) \begin{aligned} \int\left\\frac{2 G M_{E}}{c\^{2}}\\left(\\frac{1}{r}-\\frac{1}{a}\\right)\\right \frac{d s}{c} & \simeq \frac{2 G M_{E}}{c^{2}} \int\left(\frac{1}{r}-\frac{1}{a}\right) d t \\ & =\frac{2 G M_{E}}{a c^{2}} \int \frac{e \cos E}{1-e \cos E} d t \\ & =\frac{2 \sqrt{G M_{E}} a}{c^{2}} e\left(\sin E-\sin E_{0}\right) \\ & =+\frac{2 \sqrt{G M_{E} a}}{c^{2}} e \sin E+\text{常数} \end{aligned} \tag{38} ∫c22GME(r1−a1)cds≃c22GME∫(r1−a1)dt=ac22GME∫1−ecosEecosEdt=c22GME ae(sinE−sinE0)=+c22GMEa esinE+常数(38)
积分常数并入整体钟差滤波，仅保留周期性相对论钟差修正：
Δ t r = 4.4428 × 10 − 10 e a sin E 单位：秒 / 米 (39) \Delta t_{r}=4.4428 × 10^{-10}\ e\sqrt{a}\ \sin E \quad \\text{单位：秒}/\\sqrt{\\text{米}} \tag{39} Δtr=4.4428×10−10 ea sinE单位：秒/米 (39)
该修正由接收机软件实时计算；偏心率 e = 0.01 e=0.01 e=0.01卫星最大修正量23纳秒，由引力频移、二阶多普勒随轨道距离周期性耦合产生。等价矢量形式（工程软件通用）：
Δ r e l = − 2 r s a t ⋅ v s a t c 2 (40) \Delta_{rel}=-\frac{2 \boldsymbol{r}^{sat} \cdot \boldsymbol{v}^{sat}}{c^{2}} \tag{40} Δrel=−c22rsat⋅vsat(40)
r s a t , v s a t \boldsymbol{r}^{sat},\boldsymbol{v}^{sat} rsat,vsat为发射时刻ECI卫星位置、速度矢量，可由开普勒轨道公式严格证明。
GPS将偏心修正交由接收机处理；格洛纳斯直接将该修正编入卫星下发钟差参数 τ n , γ n \tau_n,\gamma_n τn,γn，接收机无需额外计算。早期卫星算力有限是GPS设计根源，如今数十家接收机厂商设备固化该逻辑，难以反向修改；但若采集原始接收机观测数据，可直接测量该相对论效应。
6 TOPEX/POSEIDON卫星相对论验证实验
航空航天公司报告声称高偏心率动态接收机（低轨卫星）无法套用式(38)(39)偏心修正，存在概念错误，若属实GPS二代、三代卫星自主导航重构成本极高。为此航空航天公司William Feess提出利用TOPEX卫星六通道GPS接收机实测该周期性相对论效应，这是唯一在轨接收机实测引力+时间膨胀耦合周期效应实验。
接收机伪距测量原理：本地复现与接收码多普勒匹配复制码，相关锁定后由码相位翻转点识别卫星发射坐标时；本地自由振荡钟存在固有偏差，需4颗卫星联立解位置+本地钟差。记卫星下发原始钟时 t j ′ t_j' tj′，叠加偏心相对论修正得到真实GPS发射时 t j = t j ′ + ( Δ t r ) j − b j t_j=t_j'+(\Delta t_r)_j-b_j tj=tj′+(Δtr)j−bj， b j b_j bj为卫星固有钟偏（含SA选择可用性扰动）；接收机本地实测接收时 t R ′ t_R' tR′，真实GPS接收时 t R = t R ′ − b R t_R=t_R'-b_R tR=tR′−bR， b R b_R bR本地钟偏。伪距定义：
ρ j = c ( t R ′ − t j ′ ) = c ( t R + b R ) − ( t j + b j − ( Δ t r ) j ) (41) \rho {j}=c(t{R}^{\prime }-t_{j}^{\prime })=c\left (t_{R}+b_{R})-(t_{j}+b_{j}-(\\Delta t_{r})_{j}) \\right \tag{41} ρj=c(tR′−tj′)=c(tR+bR)−(tj+bj−(Δtr)j)(41)
卫星接收机真实几何距离：
∣ r R ( t R ) − r j ( t j ) ∣ = c ( t R − t j ) (42) |\boldsymbol{r}{R}(t{R})-\boldsymbol{r}{j}(t{j})|=c\left(t_{R}-t_{j}\right) \tag{42} ∣rR(tR)−rj(tj)∣=c(tR−tj)(42)
联立得到实验观测方程：
∣ r R ( t R ) ∣ − ρ j + c b R − c b j + c ( Δ t r ) j = 0 (43) |\boldsymbol{r}{R}(t{R})|-\rho_{j}+c b_{R}-c b_{j}+c\left(\Delta t_{r}\right)_{j}=0 \tag{43} ∣rR(tR)∣−ρj+cbR−cbj+c(Δtr)j=0(43)
TOPEX卫星核心任务测海面高度，搭载高精度六通道双频GPS接收机与雷达高度计；轨道半径7714 km，轨道周期约6745秒，倾角66.06°，偏心率仅0.000057，近圆轨道非常适合实验。卫星、GPS星座轨道由喷气推进实验室（JPL）百站多普勒独立精密定轨，不依赖星载钟。
双频L1/L2载波相位可消除电离层延迟；JPL Yoaz Bar-Sever提供1995年10月22--23日完整观测数据集：
- 25颗卫星J2000 ECI位置速度，采样15分钟，径向精度10 cm、水平30 cm；
- TOPEX天线相位中心ECI位置速度，每分钟采样，径向3 cm、水平10 cm；
- 基于地面监测站卫星钟参数，5分钟更新，测距精度5--10 cm；
- TOPEX接收机伪距、载波相位，10秒采样，单时刻最多跟踪6颗卫星。
1995年GPS时比UTC快10秒，UTC闰秒会造成GPS时持续领先，轨道文件高次插值匹配钟文件300秒时间戳；本人与Feess独立插值，卫星位置毫米级吻合。数据预处理：载波相位平滑伪距、补偿SA干扰、修正卫星天线质心偏移。
变形式(43)求解TOPEX本地钟偏：
∣ r R ( t R ) − r j ( t j ) ∣ − ρ j − c b j + c Δ t r = − c b R (44) |\boldsymbol{r}{R}(t{R})-\boldsymbol{r}{j}(t{j})|-\rho_{j}-c b_{j}+c \Delta t_{r}=-c b_{R} \tag{44} ∣rR(tR)−rj(tj)∣−ρj−cbj+cΔtr=−cbR(44)
每一观测时刻同时观测5--6颗卫星，冗余观测精确求解本地钟偏与噪声均方误差。
图3 基于1571组观测得到的TOPEX卫星钟偏（单位：米）
钟偏整体线性漂移叠加轨道周期大尺度周期波动，该曲线为本地钟偏最优估计；变形方程分离相对论偏心项：
∣ r R ( t R ) − r j ( t j ) ∣ − ρ j − c b j + c b R = − c Δ t r (45) |\boldsymbol{r}{R}(t{R})-\boldsymbol{r}{j}(t{j})|-\rho_{j}-c b_{j}+c b_{R}=-c \Delta t_{r} \tag{45} ∣rR(tR)−rj(tj)∣−ρj−cbj+cbR=−cΔtr(45)
单独分析某颗卫星偏心效应时需剔除该星观测，否则均方误差由29 cm增至35 cm。
图4 TOPEX钟偏均值均方误差随时间变化曲线
图5 13号卫星实测相对论偏心效应（散点）与理论预测曲线对比
该卫星偏心率 e = 0.01486 e=0.01486 e=0.01486，理论与实测残差均方根22 cm，为最大效应10.2 m的2.2%。
对五颗高偏心率卫星归一（除以 e e e、时间原点平移至近地点）合并绘图：
图6 五颗高偏心率卫星归一化通用偏心效应曲线
单圈观测存在系统性通道偏差，整体理论与实测偏差约2.5%。
7 多普勒效应
卫星钟出厂已预校正恒定频偏，偏心周期修正交由接收机，地面用户仍需处理一阶纵向多普勒频移；静态引力场中电磁波传播坐标频率守恒，仅考虑地球单极、四极引力势即可推导接收频率。设卫星引力势 V j V_j Vj、速度 v j \boldsymbol{v}_j vj，发射基准频率 f 0 = 10.23 MHz f_0=10.23\ \text{MHz} f0=10.23 MHz，接收机速度 v R \boldsymbol{v}_R vR、引力势 V R V_R VR，传播单位矢量 N \boldsymbol{N} N，接收频率：
f R = f 0 1 + − V R + v R 2 / 2 + Φ 0 + 2 G M E / a + 2 V j c 2 ( 1 − N ⋅ v R / c ) ( 1 − N ⋅ v j / c ) (46) f_{R}=f_{0}\left 1+\\frac {-V_{R}+v_{R}\^{2}/2+\\Phi _{0}+2GM_{E}/a+2V_{j}}{c\^{2}}\\right \frac{\left( 1-\boldsymbol{N}\cdot \boldsymbol{v}{R}/c \right)}{\left( 1-\boldsymbol{N}\cdot \boldsymbol{v}{j}/c \right)} \tag{46} fR=f01+c2−VR+vR2/2+Φ0+2GME/a+2Vj(1−N⋅vj/c)(1−N⋅vR/c)(46)
地面大地水准静止接收机简化形式：
f R = f 0 1 + 2 G M E c 2 ( 1 a − 1 r ) ( 1 − N ⋅ v R / c ) ( 1 − N ⋅ v j / c ) (47) f_{R}=f_{0}\left 1+\\frac{2 G M_{E}}{c\^{2}}\\left( \\frac{1}{a}-\\frac{1}{r} \\right) \\right \frac{\left(1-\boldsymbol{N}\cdot \boldsymbol{v}{R} / c\right)}{\left(1-\boldsymbol{N}\cdot \boldsymbol{v}{j} / c\right)} \tag{47} fR=f01+c22GME(a1−r1)(1−N⋅vj/c)(1−N⋅vR/c)(47)
中括号项对应偏心周期相对论修正；一阶多普勒项量级 10 − 5 10^{-5} 10−5，远大于 e × 10 − 10 e×10^{-10} e×10−10量级周期修正。
8 星间测距
自主导航（Autonav）需卫星间直接信号传递坐标时；卫星预校正恒定频偏，叠加偏心修正得到坐标时后发射测距信号，接收端反向修正得到本机原子时。下标 S S S代表卫星校正后坐标时， i i i发射星、 j j j接收星，发射坐标时：
T ( i ) = T S ( i ) + 2 G M a i c 2 e i sin E i (48) T^{(i)}=T_{S}^{(i)}+\frac{2\sqrt{G M a_{i}}}{c^{2}} e_{i} \sin E_{i} \tag{48} T(i)=TS(i)+c22GMai eisinEi(48)
信号传播坐标时近似 l / c l/c l/c， l l l为发射瞬间卫星间距；忽略夏皮罗延迟，接收坐标时： T ( j ) = T ( i ) + l / c T^{(j)}=T^{(i)}+l/c T(j)=T(i)+l/c，反向得到接收卫星本机坐标时：
T S ( j ) = T ( j ) − 2 G M a j c 2 e j sin E j (49) T_{S}^{(j)}=T^{(j)}-\frac{2 \sqrt{G M a_{j}}}{c^{2}} e_{j} \sin E_{j} \tag{49} TS(j)=T(j)−c22GMaj ejsinEj(49)
联立完整星间时传公式：
T S ( j ) = T S ( i ) + l c − 2 G M a j c 2 e j sin E j + 2 G M a i c 2 e i sin E i (50) T_{S}^{(j)}=T_{S}^{(i)}+\frac{l}{c}-\frac{2 \sqrt{G M a_{j}}}{c^{2}} e_{j} \sin E_{j}+\frac{2 \sqrt{G M a_{i}}}{c^{2}} e_{i} \sin E_{i} \tag{50} TS(j)=TS(i)+cl−c22GMaj ejsinEj+c22GMai eisinEi(50)
l l l为信号传播真实路径长；改用发射时刻卫星间距 ∣ Δ r ∣ = ∣ r j ( t i ) − r i ( t i ) ∣ |\Delta r|=|\boldsymbol{r}_j(t_i)-\boldsymbol{r}_i(t_i)| ∣Δr∣=∣rj(ti)−ri(ti)∣表达：
l = ∣ Δ r ∣ + Δ r ⋅ v j c (51) l=|\Delta \boldsymbol{r}|+\frac{\Delta \boldsymbol{r} \cdot \boldsymbol{v}_{j}}{c} \tag{51} l=∣Δr∣+cΔr⋅vj(51)
代入完整星间时传方程：
T S ( j ) = T S ( i ) + ∣ Δ r ∣ c − 2 G M a j c 2 e j sin E j + 2 G M a i c 2 e i sin E i + Δ r ⋅ v j c 2 (52) T_{S}^{(j)}=T_{S}^{(i)}+\frac{|\Delta \boldsymbol{r}|}{c}-\frac{2 \sqrt{G M a_{j}}}{c^{2}} e_{j} \sin E_{j}+\frac{2 \sqrt{G M a_{i}}}{c^{2}} e_{i} \sin E_{i}+\frac{\Delta \boldsymbol{r} \cdot \boldsymbol{v}_{j}}{c^{2}} \tag{52} TS(j)=TS(i)+c∣Δr∣−c22GMaj ejsinEj+c22GMai eisinEi+c2Δr⋅vj(52)
末项与地面接收机偏心修正式(40)无关，为星间传播独有时差项。
9 轨道调整引发时钟频移
轨道机动带来半长轴变化会产生可测相对论频跳，此前长期被忽略；2000年7月SV43卫星轨道转移实测验证该效应。半长轴变化 δ a \delta a δa对应的相对频偏：
δ f f = + 3 G M E δ a 2 c 2 a 2 (55) \frac{\delta f}{f}=+\frac{3 G M_{E} \delta a}{2 c^{2} a^{2}} \tag{55} fδf=+2c2a23GMEδa(55)
半长轴减小，卫星平均引力势更低、轨道速度更大，总频偏负向；SV43机动前后轨道半长轴：
机动前2000.07.22： a = 2.656139556 × 10 7 m , v = 3873.947951 m/s a=2.656139556×10^7\ \text{m},\ v=3873.947951\ \text{m/s} a=2.656139556×107 m, v=3873.947951 m/s
机动后2000.07.30： a = 2.654267359 × 10 7 m , v = 3875.239113 m a=2.654267359×10^7\ \text{m},\ v=3875.239113\ \text{m} a=2.654267359×107 m, v=3875.239113 m
理论预测频偏 δ f / f = − 1.734 × 10 − 13 \delta f/f=-1.734×10^{-13} δf/f=−1.734×10−13，实测 − 1.85 × 10 − 13 -1.85×10^{-13} −1.85×10−13，偏差在测量误差内。2000年10月二次机动预测 + 1.75 × 10 − 13 +1.75×10^{-13} +1.75×10−13，无实测数据。
地球二阶带谐 J 2 J_2 J2带来轨道半长轴周期振荡，振幅约1658 m，周期两倍轨道周期；对应频偏振幅 0.695 × 10 − 14 0.695×10^{-14} 0.695×10−14，与轨道机动频偏同量级，高精度钟需计入。利用拉格朗日摄动方程推导椭球引力对轨道、钟频的扰动；无摄开普勒轨道单位能量守恒 v 2 2 − G M E r = − G M E 2 a \frac{v^2}{2}-\frac{GM_E}{r}=-\frac{GM_E}{2a} 2v2−rGME=−2aGME，代入频偏方程：
Δ f f = − 1 2 v 2 c 2 − G M E r c 2 − Φ 0 c 2 (53) \frac{\Delta f}{f}=-\frac{1}{2} \frac{v^{2}}{c^{2}}-\frac{G M_{E}}{r c^{2}}-\frac{\Phi_{0}}{c^{2}} \tag{53} fΔf=−21c2v2−rc2GME−c2Φ0(53)
替换速度项得到分恒定、周期两部分频偏：
Δ f f = − 3 G M E 2 a c 2 − Φ 0 c 2 + 2 G M E c 2 1 r − 1 a (54) \frac{\Delta f}{f}=-\frac{3 G M_{E}}{2 a c^{2}}-\frac{\Phi_{0}}{c^{2}}+\frac{2 G M_{E}}{c^{2}}\left\\frac{1}{r}-\\frac{1}{a}\\right \tag{54} fΔf=−2ac23GME−c2Φ0+c22GMEr1−a1(54)
首两项为出厂恒定频偏，末项积分得到 e sin E e\sin E esinE周期修正； J 2 J_2 J2引入额外周期时偏：
δ t J 2 = − G M E a 0 3 J 2 a 1 2 sin 2 i 0 2 c 2 sin ( 2 ω 0 + 2 n t ) (89) \delta t_{J_{2}}=-\sqrt{\frac{G M_{E}}{a_{0}^{3}}} \frac{J_{2} a_{1}^{2} \sin ^{2} i_{0}}{2 c^{2}} \sin \left(2 \omega_{0}+2 n t\right) \tag{89} δtJ2=−a03GME 2c2J2a12sin2i0sin(2ω0+2nt)(89)
接收机需补偿项取负号，峰值对定位误差约1.43 cm。轨道半长轴平均可精确计算轨道平均频移 Δ f / f = 3 ϵ / c 2 \Delta f/f=3\epsilon/c^2 Δf/f=3ϵ/c2， ϵ \epsilon ϵ单位轨道能量。轨道机动频移公式预测精度优于3.3%。
10 次级相对论效应
厘米级定位需额外微小相对论修正，毫米级建模必须纳入：
1. 夏皮罗传播延迟
类光世界线 d s 2 = 0 ds^2=0 ds2=0积分得到传播时滞：
Δ t d e l a y = Φ 0 c 2 l c + 2 G M E c 3 ln r 1 + r 2 + l r 1 + r 2 − l (96) \Delta t_{delay}=\frac{\Phi_{0}}{c^{2}} \frac{l}{c}+\frac{2 G M_{E}}{c^{3}} \ln \left \\frac{r_{1}+r_{2}+l}{r_{1}+r_{2}-l} \\right \tag{96} Δtdelay=c2Φ0cl+c32GMElnr1+r2−lr1+r2+l(96)
第二项经典夏皮罗延迟，第一项因坐标时基准抵消部分贡献，地星链路总修正小于2 cm，常规定位忽略。
2. 空间曲率测地距离修正
径向两点固有距离与坐标距离差值毫米量级：
∫ r 1 r 2 d r 1 + G M E c 2 r = r 2 − r 1 + G M E c 2 ln ( r 2 r 1 ) (97) \int_{r_{1}}^{r_{2}} d r\left1+\\frac{G M_{E}}{c\^{2} r}\\right=r_{2}-r_{1}+\frac{G M_{E}}{c^{2}} \ln \left(\frac{r_{2}}{r_{1}}\right) \tag{97} ∫r1r2dr1+c2rGME=r2−r1+c2GMEln(r1r2)(97)
GPS轨道与地表差值约6.3毫米，激光卫星测距比对需计入。
3 相位缠绕
GPS右旋圆极化信号，接收机自旋产生额外相位累积，等效多普勒偏移，四卫星定位时不影响解算，但高速旋转接收机需建模。
4 太阳系天体潮汐引力红移
日、月引力仅产生潮汐二阶微小频移，量级 10 − 16 10^{-16} 10−16，现有GPS系统全部忽略。
11 卫星导航增强系统
GPS信号功率低，易干扰、室内峡谷遮挡，且系统受控于美军，各国开发广域增强系统：
- WAAS（美国广域增强系统）、EGNOS（欧洲静地增强系统）、日本MTSAT：静地卫星下发电离层、钟差校正，无新增相对论项；
- QZSS日本准天顶系统：静地大偏心率轨道，相对论恒定频偏 − 5.39 × 10 − 10 -5.39×10^{-10} −5.39×10−10，偏心修正振幅290纳秒，仅覆盖亚太属增强系统。
12 全球卫星导航系统（GNSS）
- GLONASS俄罗斯格洛纳斯：轨道略低于GPS，出厂恒定频偏小，卫星下发已包含偏心相对论修正，接收机无需额外计算；
- GALILEO欧洲伽利略：卫星不预调出厂频偏，所有相对论修正交由接收机软件，多模接收机需两套相对论计算模块；
- BEIDOU北斗：早期静地转发式系统主要修正萨格纳克效应，远期独立全球兼容伽利略。
13 应用场景
GPS民用规模远超军用，覆盖电力同步、深空测控、脉冲星观测、重力实验、交通调度、测绘、农业、海洋考古、气象探测、建筑形变、地震地壳监测、搜救等；1974年双脉冲星相对论验证依托GPS共视授时；全球时间标准比对、灾害预警、自动驾驶均高度依赖厘米级定位；全球导航市场规模超300亿美元，伽利略预计年收入百亿欧元。
14 结论
GPS是狭义、广义相对论落地最完备的工程试验平台，提供大量教学实例；格洛纳斯、伽利略、北斗均依托地心局部惯性系同步理论，GPS数十年稳定运行充分验证该相对论框架。未来国际空间站搭载冷原子钟（稳定度 5 × 10 − 14 / τ 5×10^{-14}/\sqrt{\tau} 5×10−14/τ ），将开展更高精度相对论测试，进一步提升导航系统精度。
参考文献
1 "IAU Resolutions adopted at the 24th General Assembly (Manchester, August 2000)", project homepage, L'Observatoire de Paris: Le Bureau National de M´etrologie - Syst`emes de R´ef´erence Temps Espace (BNM - SYRTE), (2000). URL (cited on 20 December 2002): http://syrte.obspm.fr/IAU_resolutions/Resol-UAI.htm. 3
2 NAVSTAR GPS Space Segment and Navigation User Interfaces, Interface Control Document, ICD-GPS-200C, (ARINC Research Corporation, Fountain Valley, CA, 1993). Related online version (cited on 11 June 2007): http://www.navcen.uscg.gov/pubs/gps/icd200/. 5, 9
3 Department of Defense World Geodetic System 1984 -- Its Definition and Relationships with Local Geodetic Systems, NIMA Technical Report, TR8350.2, (National Imagery and Mapping Agency, Bethesda, MD, 2004), 3rd corr. edition. Related online version (cited on 11 June 2007): http://earth-info.nga.mil/GandG/publications/tr8350.2/tr8350_2.html. 2
4 Allan, D.W., Weiss, M., and Ashby, N., "Around-the-World Relativistic Sagnac Experiment", Science, 228, 69--70, (1985). DOI. 2
5 Alley, C., "Proper time experiments in gravitational fields with atomic clocks, aircraft, and laser light pulses", in Meystre, P., and Scully, M.O., eds., Quantum Optics, Experimental Gravitation, and Measurement Theory, Proceedings of the NATO Advanced Study Institute on Quantum Optics and Experimental General Relativity, August 1981, Bad Windsheim, Germany, NATO Science Series: B, vol. 94, pp. 363--427, (Plenum Press, New York, 1983). 5
6 Ashby, N., An Earth-Based Coordinate Clock Network, NBS Technical Note, 659; S.D. Catalog # C13:46:659, (U.S. Dept. of Commerce, U.S. Government Printing Office, Washington, DC, 1975). 2
7 Ashby, N., and Allan, D.W., "Practical Implications of Relativity for a Global Coordinate Time Scale", Radio Sci., 14, 649--669, (1979). DOI. 4
8 Ashby, N., and Bertotti, B., "Relativistic Effects in Local Inertial Frames", Phys. Rev. D, 34, 2246--2258, (1986). DOI. 10
9 Ashby, N., and Spilker Jr, J.J., "Introduction to Relativistic Effects on the Global Positioning System", in Parkinson, B.W., and Spilker Jr, J.J., eds., Global Positioning System: Theory and Applications, Vol. 1, Progress in Astronautics and Aeronautics, vol. 163, 18, pp. 623--697, (American Institute of Aeronautics and Astronautics, Washington, DC, 1996). 3, 9
10 Ashby, N., and Weiss, M., Global Positioning System Receivers and Relativity, NIST Technical Note, TN 1385, (National Institute of Standards and Technology, Boulder, CO, 1999). 2
11 Buisson, J.A., Easton, R.L., and McCaskill, T.B., "Initial Results of the NAVSTAR GPS NTS2 Satellite", in Rueger, L.J. et al., ed., 9th Annual Precise Time and Time Interval (PTTI) Applications and Planning Meeting, Proceedings of the meeting, held at NASA Goddard Space Flight Center, November 29 -- December 1, 1977, pp. 177--200, (Technical Information and Administrative Support Division, Goddard Space Flight Center, Greenbelt, MD, 1978). 5
12 Epstein, M., Stoll, E., and Fine, J., "Observable Relativistic Frequency Steps Induced by GPS Orbit Changes", in Breakiron, L.A., ed., 33rd Annual Precise Time and Time Interval (PTTI) Systems and Applications Meeting, Proceedings of a meeting held at Long Beach, California, 27 -- 29 November 2001, (U.S. Naval Observatory, Washington, DC, 2002). Related online version (cited on 11 June 2007): http://tycho.usno.navy.mil/ptti/index1.html. 9, 9, 9, 9, 9
13 Fitzpatrick, P., The Principles of Celestial Mechanics, (Academic Press, New York, 1970). 5, 9, 9, 9
14 Fliegel, H.F., and DiEsposti, R.S., GPS and Relativity: an Engineering Overview, The Aerospace Corporation Report, ATR-97(3389)-1, (The Aerospace Corporation, El Segundo, CA, 1996). 6
15 Garland, G.D., The Earth's Shape and Gravity, (Pergamon Press, New York, 1965). 3
16 Hulse, R.A., and Taylor, J.H., "Discovery of a pulsar in a binary system", Astrophys. J., 195, L51--L53, (1975). DOI, ADS. 13
17 Kraus, J.D., Antennas, (McGraw-Hill, New York, 1988), 2nd edition. 10
18 Lambeck, K., in Geophysical Geodesy: The Slow Deformations of the Earth, pp. 13--18, Oxford Science Publications, (Clarendon Press, Oxford; New York, 1988). 3
19 Malys, S., and Slater, J., "Maintenance and Enhancement of the World Geodetic System 1984", in Proceedings of the 7th International Technical Meeting of The Satellite Division of The Institute of Navigation (ION GPS-94), September 20 -- 23, 1994, Salt Palace Convention Center - Salt Lake City, UT, pp. 17--24, (Institute of Navigation, Fairfax, VA, 1994). 2
20 Mashhoon, B., "Nonlocal Theory of Accelerated Observers", Phys. Rev. A, 47, 4498--4501, (1993). DOI. 10
21 Mashhoon, B., "On the Coupling of Intrinsic Spin with the Rotation of the Earth", Phys. Lett. A, 198, 9--13, (1995). DOI. 10
22 Parkinson, B.W., and Spilker Jr, J.J., "Overview of GPS Operation and Design", in Parkinson, B.W., and Spilker Jr, J.J., eds., Global Positioning System: Theory and Applications, Vol. 1, Progress in Astronautics and Aeronautics, vol. 2, pp. 29--56, (American Institute of Aeronautics and Astronautics, Washington, DC, 1996). 6
23 Taylor, J.H., "Binary Pulsars and Relativity", Rev. Mod. Phys., 66, 711--719, (1994). DOI. 13
24 Tetewsky, A.K., and Mullen, F.E., "Effects of Platform Rotation on GPS with Implications for GPS Simulators", in Proceedings of the 9th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GPS-96), September 17 -- 20, 1996, Kansas City Convention Center, Kansas City, Missouri, pp. 1917--1925, (Institute of Navigation, Fairfax, VA, 10
25 Galileo Open Service Signal In Space Interface Control Document (SIS ICD), (European Space Agency / Galileo Joint Undertaking, Brussels, 2006). Related online version (cited on 6 November 2009): http://www.gsa.europa.eu/go/galileo/os-sis-icd. (document), 12