GPS定位这件事很多人天天在用但真正拆开来看它其实就是一句话卫星把“我在哪里、现在几点”用无线电波告诉你你的手机把这些信息收齐再算一算自己在哪里。这个过程中GPS信号负责把信息从两万公里外送下来导航电文负责把“卫星的位置和时间”说清楚。这篇笔记主要想聊透这两件事——信号长什么样电文里装了些什么以及接收机到底怎么把它们变成坐标。GSP系列的笔记如果只讲原理不讲链路总觉得隔了一层。所以我这篇会顺着“信号生成—播发—接收—解调—电文解读”这条线往下走中间穿插一些我在实际解调信号时踩过的坑和总结的参数规律。内容适合刚接触卫星导航、想搞懂接收机到底在做什么的人也适合做定位开发、做信号处理但一直没时间系统梳理底层的朋友。看完你应该能回答这几个问题C/A码和P码是干嘛用的导航电文的1500个比特是怎么组织的接收机又是怎么从一堆噪声里把卫星数据抠出来的1. GPS信号从卫星到手机的物理链路1.1 频率选择背后的工程逻辑GPS卫星播发的信号核心频率有三个L1是1575.42MHzL2是1227.60MHzL5是1176.45MHz。民用最常见的还是L1因为C/A码就调制在这上面。为什么选L波段这个频段有天然的工程优势——电离层延迟小、波长适中、大气吸收弱、天线尺寸合理。如果频率太高信号容易被雨衰和大气吸收吃掉太低的话天线又大得离谱而且电离层折射会严重到没法用。L1频率1575.42MHz是怎么来的它其实是基准频率10.23MHz的154倍。这个基准频率还派生出一系列关键参数C/A码码率1.023MHzP码码率10.23MHz都是基准频率的整数倍。整个GPS系统像一个精密的时钟工厂所有信号的时序都咬合在同一套时钟体系里这给接收机做时间测量带来了极大的便利。你可以把基准频率理解成一副齿轮的主轴其他频率全部挂着它的倍数转谁也不会乱。实际接收的时候卫星发射功率约50W但经过两万公里自由空间衰减到达地面的信号功率大约只有-160dBW量级。这个数字有多小差不多是手机收到基站信号的百亿分之一比热噪声底还低。GPS接收机要从噪声里把信号捞出来靠的就是扩频增益——这个下面细讲。1.2 C/A码与P码扩频通信的“暗号”GPS信号里最核心的扩频码包括C/A码和P码。C/A码是粗捕获码码率1.023Mchip/s码长1023个码片周期正好1毫秒。P码是精密码码率10.23Mchip/s周期长达约7天主要给军用和高精度用户使用一般民用接收机只用到C/A码。为什么需要扩频码这里用生活类比最好理解想象一个巨大的房间里几百个人同时在说话。如果大家都用同一种语言谁也听不清谁。但如果每个人用不同的、互不相关的“暗语”听的人只要对准其中一个暗语就能把对应的那个人说的话从嘈杂声中分离出来。C/A码就是卫星和接收机之间约定的暗语。GPS卫星虽然各自播发相同频率的信号但每颗卫星分配不同的C/A码——这叫码分多址。C/A码的具体生成方式是Gold码由两个10级移位寄存器生成每个码周期1023个码片不同卫星用不同的相位抽头组合互相之间的相关性极低。接收机本地复制出对应卫星的C/A码与接收信号做相关运算只有当两个码相位严格对齐时相关峰才会出现。这个峰就是捕捉到卫星的标志同时也是测量伪距的基础。P码的码速率是C/A码的10倍码片更窄所以理论测距精度也高约10倍。但它结构复杂还有加密的Y码版民用领域基本接触不到。我对P码的态度是理解它的存在价值但不必深挖序列结构因为实际工程中C/A码已经能满足绝大多数定位需求了。1.3 调制方式与信号功率分配GPS信号在L1上同时调制了C/A码和P码两者共用载波但相位正交——一个是I支路一个是Q支路互不干扰。调制方式用的BPSK二进制相移键控原理简单说就是码片是0就发0度相位码片是1就发180度相位。接收机解调时靠“翻转”来还原码序列。L1信号的总功率在那50W里要怎么分典型分配是C/A码约占50%P码约占50%两者功率相近。这样做的好处是军用和民用都能获得稳定信号强度不至于谁抢了谁的资源。L2上早期只播发P码后来增加了L2C民用信号L5则是新一代民用信号带宽更大、抗干扰更强。一个容易被忽略但非常关键的点是卫星上所有信号共用同一个原子钟。无论是L1还是L2无论是C/A码还是导航电文它们的时间基准完全一致。接收机测到C/A码的相位本质就是在测量信号从卫星到接收机的传播时间这个传播时间乘以光速就是伪距。没有统一的时钟体系整个测距过程就无从谈起。2. 信号捕获与跟踪定位前的必经之路2.1 多普勒效应与二维搜索接收机要解调导航电文第一步不是直接读比特流而是先完成“捕获”。捕获要解决两个问题一是找到卫星C/A码的码相位二是估计当前信号的多普勒频移。为什么会有多普勒因为卫星不是静止的它以约3.87km/s的速度绕地球飞行。卫星靠近你时收到的信号频率偏高远离时频率偏低。L1频率上最大的多普勒偏移大约在±5kHz左右如果接收机在高速移动的载体上还要叠加载体运动带来的频移总范围可能扩展到±10kHz。捕获过程本质上是一个二维搜索过程。C/A码一个周期有1023个码相位多普勒频率需要按一定步进扫描——一般步进取500Hz。频率范围±10kHz就是40个频点再乘以1023个码相位接收机要搜索约四万个“码相位-频率”组合单元。每颗卫星搜一遍整个搜索空间相当可观。现代接收机靠的是并行相关器才能把搜索速度压下来。相关器数量从几十到上千不等。每个相关器可以同时检查一个频率单元下所有码相位的相关性。我在实际测试某跨平台系统时发现入门级接收机冷启动搜索时间一般在30秒到几分钟不等而高端接收机因为相关器多、算法优化好可以把冷启动压缩到几秒。这个差异本质上就是硬件并行度和软件策略差异的体现。2.2 跟踪环路载波环与码环的配合捕获完成后接收机已经知道大概的码相位和多普勒频率但这只是粗略估计。接下来要进入“跟踪”阶段精确锁定载波相位和码相位。跟踪环节有两条环路载波跟踪环和码跟踪环。载波环负责把多普勒频移去除把信号搬到基带码环负责精确对齐本地C/A码与接收信号让相关峰值最大化。码环通常用的是延迟锁定环DLL里面有早、即时、晚三路相关器。早路提前半个码片晚路延后半个码片即时路正好对齐。当早路和晚路的相关值相等时说明本地码精确对齐如果不等环路滤波器就会产生一个误差信号驱动本地码发生器朝正确方向调整。这个过程每毫秒都在跑保证码相位误差被控制在可接受范围。载波环分两种锁相环PLL和锁频环FLL。PLL精度高但对动态应力敏感FLL抗动态性能好但精度差。实际接收机通常先用FLL粗锁频率再切到PLL精锁相位。数据解调时导航电文比特就隐藏在即时路相关器输出的符号翻转里——比特翻转一次相位变化180度接收机通过检测这种翻转来恢复电文。2.3 捕获跟踪过程中的常见坑我在做信号采集和基带处理实验时遇到过几个非常典型的坑这里整理一下采样率必须是C/A码码率的整数倍吗不一定但采用整数倍会极大简化码相位对齐逻辑。我用过非整数倍的采样率结果捕获后的码相位和变频导致的残留频偏纠缠在一起排查了很久才发现是采样率非整数倍导致本地码发生器频率与真实码率略有偏差。捕获阈值设太低会把噪声误判成卫星。早期的调试程序里我把相关峰的门限设成固定值结果在城市峡谷环境里频繁误捕获定位结果乱跳。后来改成自适应门限以当前噪声底为参考动态调整误捕获率大幅下降。环路带宽不能拍脑袋定。带宽太宽噪声大带宽太窄动态响应慢。静态接收机可以接受1Hz以下的带宽车载动态场景至少需要几Hz。我通常是做成可变带宽信号刚捕获时用宽带宽锁定稳定后再收窄。3. 导航电文定位所需的“数字大脑”3.1 电文结构帧、子帧与字的组织导航电文是GPS信号里真正“说话”的内容。它告诉接收机卫星当前在哪、它的钟准不准、卫星是否健康、整个星座的大致分布等信息。没有这些信息就算把C/A码跟踪得再好也换不出坐标。GPS导航电文的基本结构是一个超帧包含25个帧每一帧1500比特播发速率50bps所以一帧恰好播30秒。25帧完整播完需要12.5分钟。每一个帧分成5个子帧每个子帧300比特播发6秒。每个子帧又分成10个字每字30比特播发0.6秒。前三个子帧实时播发的是当前卫星自己的星历和时钟修正参数每30秒重复一次接收机只要连续接收30秒就能拿到完整定位所需参数。第四子帧是电离层修正参数、UTC参数等第五子帧则用来循环播发25颗卫星的历书——历书是粗略的轨道参数用于快速搜星和预判可见卫星精度远不如星历。接收机要收完整整25帧才拼出全星座历书这就是为什么冷启动搜星比热启动慢那么多——热启动时接收机已经把历书存进内存了。3.2 子帧内容拆解星钟参数与开普勒星历导航电文里对定位精度影响最大的是子帧1里的时钟修正参数和子帧2、3里的星历参数。卫星时钟修正采用二阶多项式模型。电文里播发三个系数af0、af1、af2分别对应钟差、钟速和钟漂变化率还附带一个参考时间toc。接收机算出当前时刻t与toc的差值Δt卫星钟差就等于af0 af1×Δt af2×Δt²。这里的单位是秒算出来要乘上光速才是距离误差。GPS卫星上的原子钟虽然极其精准但依然会有微小的漂移如果不修正每毫秒的钟差就对应300公里的测距误差。星历则是描述卫星椭圆轨道的参数集合。传统开普勒根数有六个轨道长半轴的平方根、偏心率、轨道倾角、升交点赤经、近地点角距、平近点角。GPS广播星历在这基础上还增加了轨道摄动修正参数比如纬度辐角修正项、轨道半径修正项、倾角修正项以及轨道平面进动率等。接收机拿到这组参数后通过开普勒方程迭代求解卫星在轨道面上的位置再经过坐标旋转到地心地固坐标系就得到卫星在某一时刻的精确位置。整个计算过程在接收机里通常几百微秒就能完成但逻辑链路很长任何一个参数解错定位结果都会偏到离谱。3.3 电文的时间同步与比特同步导航电文的每个比特持续20毫秒50bps就是50个比特每秒。接收机在完成码跟踪之后还要做比特同步——找到每个比特的边界。因为C/A码周期是1毫秒一个比特正好跨20个码周期。接收机通过检测相关值符号的翻转位置来对齐比特边界。这一步如果没做对后续解帧就会全乱。在电文层次还有比特同步到字同步、字同步到子帧同步的逐级过程。每个子帧的第一个字是遥测字前8位是固定的前导码接收机靠搜索这段时间充当前导码的固定模式来实现字同步。第二个字是交接字里面包含周内秒和子帧号。周内秒是GPS时间体系里从每周日零点开始数的秒数范围0到604799秒。有了这个信息接收机才知道当前电文处在哪一周的第几秒整个系统才能对齐到GPS时间上。关于时间有一点特别重要GPS时间不是协调世界时。GPS时是连续的原子时标不插入闰秒。协调世界时会因为地球自转逐渐变慢而不定期插入闰秒导致两者差出整数秒。导航电文里有UTC参数接收机拿到后可以做转换。我在第一次做时间转换时忘了加闰秒结果定位坐标差了好几百米——因为时间错了卫星位置算出来就是错的。后来我习惯在任何时间处理逻辑里都显式检查一下闰秒表而不是用硬编码。3.4 电文里那些容易被忽略的辅助参数除了星历和星钟电文里还有好几个容易被忽略但关键时刻能救命的参数电离层延迟修正参数一组8个系数alpha和beta各4个接收机可以用它们代入Klobuchar模型估算电离层延迟。单频用户靠这个补偿能修正掉约50%以上的电离层误差。卫星健康状态包含在子帧1和子帧4/5里接收机应跳过不健康卫星。我在调试测试软件时见过有人没检查健康标志结果把一颗标记为不健康的卫星数据也拿去做定位解算坐标虽然能用但精度明显变差。URA用户测距精度指数这个值告诉接收机这颗卫星的测距精度预期。高精度应用可以把URA差的卫星降权或剔除简单定位时可能无所谓但做精密定位时这是个好用的筛选条件。4. 实操解析从信号到定位坐标4.1 搭建一套信号解析验证环境要说我从零到一验证GPS信号链路的过程其实没想象中复杂。我建议想深入理解的读者都试一遍工具链完全开源免费你需要的只是一台能接收GPS信号的软件接收机可以用通用软件无线电外设一台电脑加上一套开源的GPS基带处理代码。大致流程是这样的先用软件无线电设备在L1频段采集一段原始采样数据采样率建议至少4MHz。然后把数据交给开源GPS软件接收机处理。这个接收机会先做捕获找出可见卫星的码相位和多普勒频率接着进入跟踪输出每颗卫星的I/Q相关值再做比特同步和帧同步解出导航电文最后解析星历参数并计算卫星位置和接收机位置。我实测过用2.4MHz采样率解调C/A码勉强可行但捕获灵敏度不理想城里树多的地方容易丢星。升到4MHz以上才稳定很多。采样率的选择本质上是在数据量、运算量和灵敏度之间找平衡。你采集一秒的4MHz采样数据16位量化的话就是8MB处理起来不多不少正好适合在普通电脑上实时跑。4.2 电文解析中的关键细节解调出导航电文比特流之后下一步就是做“帧解析”。这里有几个可以直接抄的细节前导码搜索不要死等。实际解码时比特错误不可避免如果你严格要求一段连续比特与前导码完全一致可能会漏掉好几帧。建议允许前导码有一定容错比如8位里允许1位翻转同时用奇偶校验码做二次确认。GPS每个字的最后6位是奇偶校验位专门用来验证前30位数据的有效性。星历参数的解析必须注意“尺度因子”。导航电文为了节省比特所有星历参数都用定点数编码每个参数都有特定的尺度因子比如偏心率E的尺度因子是2的负21次方。如果你忘乘尺度因子算出来的卫星位置会错得离谱。我见过不少新手在这个地方翻车排查了很久才发现是尺度因子没乘。子帧号必须校验。交接字里包含了子帧号如果你解到的子帧号与当前正在解析的数据位置对不上说明帧同步出了问题直接丢弃重新搜索前导码是明智的做法别强行去拼帧。这一整套下来你会发现一个很有意思的点真正决定定位精度的往往不是捕获和跟踪而是导航电文里那些参数的解析精度。信号处理只要做到“能解调、不错位”剩下的就是数据解读的活了。4.3 常见问题与排查技巧实录这里把我在调试中最常遇到的疑问和解决办法整理成一张速查表直接对照着查就行现象 | 可能原因 | 排查方向 捕获不到卫星信号 | 采样率过低、信号功率太弱、多普勒估计步进偏大 | 提高采样率、检查天线位置、缩小频率搜索步进到250Hz 能捕获但跟踪不稳定 | 环路带宽不合适、动态应力过大 | 载波环带宽调整为5Hz左右码环带宽1Hz左右或加入3阶环路 电文帧同步失败 | 比特同步错误、前导码容错不足 | 检查比特边界增大前导码容错位数确认奇偶校验逻辑正确 输出坐标漂移严重 | 某颗不健康卫星参与解算、电离层模型未启用 | 过滤不健康卫星启用Klobuchar模型修正 时间转换后位置偏差大 | 闰秒未更新、UTC参数未解析 | 更新闰秒表检查UTC参数解析逻辑在实际调试中我最大的一次教训是以为数据采集好了就万事大吉结果软件接收机处理出来的卫星位置始终差了十几公里。查了两天才发现采集设备的分频器配置错误导致采样率实际是3.996MHz而不是标称的4MHz。接收机本地码发生器按4MHz来假设码速率结果码相位不断缓慢漂移跟踪环一路在补偿这个偏差。这个案例提醒我采集设备的实际时钟必须校准不能盲目相信配置。另一个经验是在城市峡谷环境里多径效应会让跟踪环路产生偏置。这时候你从相关峰形状能看出门道——正常情况相关峰是对称的有强多径时会看到副峰或拖尾。我之前是直接看早路和晚路的差值是否稳定一旦出现持续性偏差就判定该颗星可能存在多径将其从解算中剔除。位置精度能从几十米提高到十米以内。如果你打算自己做电文解析一个小建议是先用模拟信号源生成测试场景而不是一上来就接真实天线。模拟信号的参数是已知的方便你逐项验证解析逻辑。等解析结果和模拟源完全一致了再接天线做实测。这个流程能帮你把“软件问题”和“环境问题”分开排查省掉无数调试时间。4.4 定位解算层面的衔接电文解析完成后定位解算其实是一个成熟的数学问题。接收机至少需要4颗卫星的伪距测量值因为未知数有四个接收机三维坐标加接收机钟差。伪距等于几何距离加上接收机钟差引起的光速偏差再经过电离层和对流层修正。实际计算中接收机会先做一个粗略的位置估计然后用牛顿迭代法对非线性观测方程逐次线性化求解直到位置修正量小于某个阈值。解算过程中要用到加权最小二乘或卡尔曼滤波权重取决于前面提到的URA和仰角。低仰角卫星的电离层误差和多径误差都大权重应该降低。这是导航解算的基本素养很多新手容易忽略仰角权重结果定位误差被低仰角卫星拖大。迭代求解时还有一个容易被忽略的坑初始位置不能随便给。如果初始位置离真实位置太远牛顿迭代可能发散。我通常先用伪距最小二乘的闭式解做个粗定位把结果作为牛顿迭代的初值这样即使城市峡谷环境也能稳定收敛。5. 关于GPS信号学习路线的个人建议如果你看完笔记想进一步深入我给一条自认为高效的实践路线先从C/A码捕获开始写一个最小实现只处理一个卫星、一个时刻的相关结果不用管跟踪。然后再逐步增加多通道跟踪、电文解析、位置解算每加一层都拿模拟信号验证。这条路我走过虽然前期会比较枯燥但每一步的反馈都清晰可见比直接铺开整套接收机代码要容易坚持得多。我自己的体会是GPS信号和导航电文这两块表面上是通信和编码问题本质上其实是时间测量和几何解算问题。你把“时间戳怎么打”“卫星时间怎么对”这两件事想透了整个GPS原理就贯通了一大半。下一篇笔记我打算写GPS定位中的误差源和差分定位到时候就可以把这篇讲到的电文参数、伪距测量继续往前推一步。
