做GNSS接收机或者导航定位应用的朋友十有八九会遇到一个现象终端同时收着GPS和北斗定位解算出来的坐标明明在一个地方可两个系统的伪距观测值放在一起比对总差着一截。这个差不是随机噪声更像一个稳定的偏移。前段时间我正好整理了一批L1/B1频段的实测数据把GPS L1和北斗B1的伪距偏差从原理到数据处理完整过了一遍有些结论和排查方法值得单独拿出来说说。先说清楚这东西是什么、能解决什么问题。GPS有一个经典民用频点L1中心频率1575.42 MHz北斗这边对应的是B1频段最常用的B1I中心频率1561.098 MHz。追新信号的还会接触B1C1575.42 MHz但本篇重点讨论的是B1I因为市面上大量北斗模块和板卡依旧以B1I为主。伪距偏差简单理解就是接收机输出的伪距观测值家庭条件同等的情况下跟“理想真值”之间存在的一个系统性偏移。偏了多少、为什么偏、怎么把它测出来再消掉直接影响双系统组合定位的精度和PPP精密单点定位的收敛速度。适合正在做多系统接收机算法、双频天线设计以及从事高精度定位测试的朋友参考。1. 为什么要把GPS L1和北斗B1放在一起测1.1 两个频点听上去像实际差得不少GPS L1和北斗B1I都是民用导航频点码率也都是1.023 McpsC/A码和B1I码的周期都是1毫秒码长都是1023个码片。频率差多少1575.42 MHz减去1561.098 MHz约14.32 MHz。这个数值看着不大但放在射频链路里已经是不可忽略的差别了。简单算一下相对带宽14.32 MHz除以1575.42 MHz约0.9%。对于一个带宽十几兆的接收机来说这两个频点分别处在滤波器通带的不同位置群延迟响应不太可能一样。我再举个例子帮助理解。射频前端的带通滤波器、声表面波滤波器SAW、低噪声放大器LNA都是有相位频率特性的信号的群延迟是频率的函数。如果一个滤波器通带中心设计在1575 MHz附近那么GPS L1信号经过它时落在群延迟曲线比较平坦的区域而B1I信号落在边缘或过渡带延迟就可能差出十几纳秒甚至几十纳秒。光是一个纳秒换算成伪距就是0.299米。前端处理不当L1/B1之间差出数米伪距偏差一点都不稀奇。1.2 伪距偏差到底是什么偏差伪距观测方程里接收机输出的伪距包括了几何距离、卫星钟差、接收机钟差、电离层延迟、对流层延迟、硬件延迟和多径误差。伪距偏差核心其实就是硬件延迟差。同一个接收机处理不同系统信号时码生成器、相关器、滤波器路径不一致会引入一个相对固定的额外延迟。这个额外延迟在业内有两个容易混的术语一个是码间偏差Inter-Signal Bias另一个是频间偏差Inter-Frequency Bias。GPS L1和北斗B1频点不同、信号结构也不同所以偏差同时具备“频间”和“码间”属性。如果放在同一个接收机内部来看它的典型范围可以从几厘米到十几米不等完全取决于接收机射频方案、基带配置和天线设计。最关键的一点是这种偏差不能用几何模型直接消除只能靠标定。标定的过程就是实测分析。这也是为什么我坚持要把GPS和北斗放到同一套测试环境里做对比的原因。2. 频点差异到偏差的传导链路2.1 射频前端是最大“贡献者”接收机从天线收到信号经过带通滤波器、LNA、混频器到ADC采样。每一个环节都会对信号引入幅度和相位变化。其中带通滤波器的群延迟特性对伪距影响最直接。我用一个常见场景说明。很多低成本接收机为了抗带外干扰会选用带宽较窄的SAW滤波器比如16 MHz带宽。这种滤波器在通带边缘的群延迟变化特别剧烈。GPS L1和北斗B1I相隔14.32 MHz如果滤波器中心频率没有精确对准某一个频点另一个频点就相当于工作在通带边缘。结果就是两个频点的信号到达基带处理器的时间基准不一样伪距直接产生常量偏移。更麻烦的是群延迟还会随温度漂移。设备从冷启动到稳定工作内部温度变化十几度滤波器的群延迟可能挪动零点几纳秒到几纳秒。对应伪距偏移就是分米级到米级的慢慢变化。做高精度测量时这种温度相关性会体现在伪距残差的低频漂移上很容易被误判成钟差变化。2.2 相关器和码环里的“隐性延迟”基带部分同样有坑。接收机对卫星信号的伪距测量本质是本地码和接收码做相关找到相关峰最大值的位置。这个位置要对应到接收机的时间基准上一般是时钟的某一拍。但从相关器输出到伪距最终输出中间还要经过码环滤波器、载波环辅助、观测量生成等步骤。GPS C/A码和北斗B1I码虽然码率相同但码发生器的实现方式、电文调制方式不同导致接收机在跟踪环路里对其中的码相位估计会有细微差别。尤其是在低载噪比环境下码环滤波器带宽不同码相位估计的收敛特性不同伪距输出会出现一个缓慢变化的偏移。在强信号场景下这种偏移可以忽略但在遮挡环境或弱信号条件下伪距偏差会明显增大。如果接收机支持载波相位测量通常会用载波环辅助码环来平滑伪距。这个辅助权重在不同系统之间也不一定相同所以即便射频前端校准得很干净基带环路配置差异仍然会留下残余偏差。2.3 电离层延迟的频间差也不是小数目L1和B1I虽然频率接近但电离层延迟是频率的函数延迟量与频率平方成反比。对比一下假如当前电离层对GPS L1的延迟是3米那么对B1I的延迟约为3乘以1575.42 / 1561.098的平方约等于3.055米。这中间差了5到6厘米。理论上电离层延迟可以用双频组合或者模型修正掉但在单频接收机里如果直接用两个系统的原始伪距做差这5到6厘米就会混在系统间偏差里面。所以做实测分析时不要把原始伪距直接相减然后宣布这就是硬件偏差。必须先扣除电离层延迟的频间差异否则结果既包含硬件偏差也包含电离层残余后续用这个值去做系统间偏差补偿就没法干净。3. 实测方案设计3.1 测试环境与设备选择伪距偏差实测最怕环境不可控多径、天线相位中心、馈线差异都会把结论搞脏。我这次采用的方案是零基线对比测试。天线使用一只双频测量型天线同时支持GPS L1和北斗B1相位中心对两个频点的一致性良好。馈电天线输出后接一个高隔离度功分器分成两路。接收机一台多系统高精度接收机能同时输出GPS L1 C/A和北斗B1I的伪距/载波观测值。这样做的好处是两台接收机的时基不同步问题不存在所有偏差都归到同一个接收机内部路径上。天线架设位置我选了一个楼顶的开阔观测点四周没有高反射物仰角5度以上的天空基本干净。如果手上只有单系统接收机也可以用两台接收机共享一根天线但它们之间需要精确的时频同步通常用PPS脉冲对齐操作复杂度高不少。建议有条件的优先用一台多系统接收机做同机对比数据质量更可控。3.2 数据采集与配置采集时长建议至少24小时覆盖不同卫星几何分布和电离层变化周期。设置采样间隔为30秒仰角掩角设为5度。关键配置项是确认接收机输出的是原始观测值而不是已经做了系统间偏差修正的“净化”数据。关闭接收机内置的低精度多径抑制算法如果可配置或者在处理端记录算法开关状态。记录载噪比和卫星仰角用于事后筛选。保存星历用于解算卫星位置和钟差。同时记录设备温度后面分析偏差漂移时温度数据非常有价值。3.3 伪距偏差的提取流程拿到原始数据后我按以下步骤处理用广播星历或精密星历计算每颗卫星的三维位置。已知测站精确坐标计算出每颗卫星到测站的几何距离。对每颗卫星、每个历元的伪距观测值做标准修正扣除卫星钟差、相对论效应、对流层模型延迟、地球自转效应。计算单系统单星的伪距残差。这个残差里包含了接收机钟差、硬件延迟、电离层残余和多径误差。做星间单差或者系统间单差来消除接收机钟差。采用同历元下同一颗卫星的GPS和北斗观测值进行系统间单差效果非常直观。差值的均值就是两个系统的伪距偏差差值的抖动则反映了多径和噪声水平。4. 实测结果与分析4.1 静态观测的偏差统计以我这次实测的某型双频接收机为例。选取同一颗同时被GPS和北斗跟踪的卫星实际是两颗分别对应GPS和北斗的卫星用系统间单差处理24小时数据后B1I伪距相对L1伪距的偏差大约稳定在2.31米。连续24小时内的标准差是0.09米最大波动0.42米。这个0.42米的波动主要来自低仰角时段的多径高仰角卫星时段偏差变化通常小于0.2米。也就是说对这台接收机而言L1/B1伪距偏差是一个相当稳定的常数完全可以通过标定后扣除。如果接收机和天线是分离购买的偏差数值可能完全不同。比如我之前遇过一套组合天线是宽带的接收机射频前端中心频率按GPS L1优化结果B1I伪距偏差达到7.8米。所以切勿把别人设备标定出的偏差值直接搬到自己系统上。4.2 不同仰角和多径环境下的表现把数据按仰角分成5到30度、30到60度、60到90度三档统计能看到明显规律。仰角范围伪距偏差均值米抖动标准差米说明5° - 30°2.450.21低仰角多径影响明显30° - 60°2.330.08偏差趋稳60° - 90°2.280.06接近硬件真实偏差低仰角数据段的偏差均值比高仰角段高了约0.17米这个差异主要是多径叠加在伪距上形成的正向偏移。做系统间偏差标定时建议只采用30度仰角以上的数据否则会把多径误差一并拟合进去。另外我也做了遮挡环境下的快速检查把天线挪到靠墙位置只留半边天空同样的接收机硬件B1I相对L1的伪距偏差均值上升了0.6米左右抖动增大到0.6米以上。这说明多径环境对系统间偏差的短期计算有显著干扰实际使用中尽量在开阔环境标定。4.3 频间偏差IFB的提取处理过程中有一个关键细节系统间单差之后还要扣除电离层频间差异。具体做法是先用双频观测值如果接收机支持或模型估算L1和B1I各自的电离层延迟然后把B1I伪距统一折算到L1频率上再与L1伪距做差。折算公式的本质是利用电离层延迟反比于频率平方的关系将B1I伪距值对应的电离层误差修正到L1等价的水平避免把电离层残余当作硬件偏差。本次静态实测中电离层当天属于中等活动水平。如果不做频间修正直接做差偏差数值比修正后多出约5厘米在电离层活跃日这个差异可以达到10到20厘米。因此高精度应用中单频条件下的IFB标定必须结合电离层模型进行不能做“裸差”。5. 常见问题与排查技巧实录5.1 天线相位中心与馈线长度不一致实测中最常被忽略的是天线相位中心误差。不同天线对不同频段的相位中心位置并不完全重合可能差出几个毫米到几个厘米。对伪距偏差来说这相当于一个常量偏移很难从数据里自动分离。馈线的问题更隐蔽。如果两路接收机各自接一根不同长度的馈线或者使用不同材质的线缆信号到达接收机的时间基准就会不一样引入额外伪距差。一组经验数据室1米标准RG58馈线在L1频段大约造成1.4纳秒的延迟折合0.42米伪距差B1频段也有类似量级。零基线测试必须保证馈线长度一致最好使用同一根天线通过功分器连接。5.2 AGC自动增益控制引发的偏差变化如果接收机前端AGC调节范围过大会在强信号和弱信号场景下改变射频链路的增益分配从而影响不同频点的群延迟。典型表现是户外开阔地标定时偏差是2.3米把接收机拿到室内窗边后同一套硬件的偏差变成了2.9米而且发射不稳定。遇到这种情况先用信号模拟器在可控载噪比下测试逐档调节信号功率记录伪距偏差变化。如果偏差随功率明显摆动说明AGC设置需要调整例如给AGC设定一个合理的静态工作点。少数专业接收机允许手动锁定增益这类测试环境建议直接锁定。5.3 卫星钟差与星历误差混入还有一个踩过多次的坑处理星历的时候直接用广播星历解算卫星位置和钟差但GPS和北斗广播星历的钟差精度和参考时间基准不同。北斗广播星历的钟差精度通常略低于GPS短期变化也更明显。如果处理窗口内正好遇到卫星钟跳变单差值会产生一个持续时间较短的脉冲容易污染统计结果。所以我处理时对每个单差序列都做了中位数滤波同时剔除了偏离中位数超过3倍标准差的历元。对于实时测试如果偏差值出现突然跳动优先查看当前历元参与解算的卫星数、卫星健康状态以及星历龄期。5.4 排查流程速查表现象优先排查方向快速验证方法偏差整体偏大超过10米射频滤波器群延迟、功分器隔离度换用宽频参考天线/功分器对比偏差随时间漂移设备温度变化、AGC状态记录温度曲线强制冷却后重测偏差随卫星仰角变化明显多径、天线相位中心只统计30度以上数据移动天线位置偏差出现短时脉冲卫星钟差跳变、星历切换检查卫星健康状态剔除异常历元不同接收机偏差不一致基带码环配置、滤波器带宽用同一射频前端分别接两台接收机后处理6. 实测之外几个容易被带偏的细节6.1 陶瓷天线设计对双频偏差的影响最近不少做嵌入式导航模块的朋友问L1/B1双频陶瓷天线怎么选。陶瓷天线的介电常数和尺寸决定谐振频率一般单频L1陶瓷天线尺寸是25毫米乘以25毫米中心频率对准1575.42 MHz。想让B1I频点也工作良好最简单的办法是选择宽频陶瓷天线但带宽拓宽后天线增益在整个通带内不会太均匀边缘频点增益下降可能达到1到2 dB直接影响载噪比。更重要的是天线匹配网络。双频陶瓷天线在馈电点通常需要增加匹配电路来兼顾两个频点如果匹配不好反射系数在其中的某个频点会变大等效改变信号幅度和相位响应这对伪距偏差的稳定性相当不利。实测中我遇到过一次双频陶瓷天线在B1I频点的回波损耗只有-6 dB比L1差了3 dB以上接收机输出的B1I伪距噪声明显变大偏差抖动也因此增大。6.2 处理软件选择的坑处理伪距偏差数据有人习惯用商业软件一键出结果但我更建议用开源或自编脚本分步处理原因很简单你要清楚每一步改了什么。市面上一些后处理软件会在后台自动做系统间偏差改正或者电离层模型修正这会导致你拿到的“测试结果”其实是软件处理后的结果而不是原始硬件特性。做标定的人一旦把这种产品化的偏差回代到自研定位算法里就会两头打架。我现在的做法是先导出一个历元原始观测值自己写脚本重新走一遍解算流程确保每个修正项都有据可查。数据量大的时候再用自动化脚本批量处理并在输出中添加中间变量方便后期查验。6.3 接收机固件升级后记得重新标定这个细节非常容易忽视。某型接收机在一次固件更新后官方优化了B1I码环滤波器参数整机对北斗B1I的通道延迟发生了变化。因为我及时对比了升级前后的零基线数据才发现偏差变了约0.8米。如果不重新标定系统间偏差表就是旧的用在新固件上组合定位的融合结果会莫名其妙变差。所以我的建议是每次固件或配置参数变更哪怕只改了环路带宽、跟踪策略都要做一次快速零基线对比测试。成本不高但能帮你避免很多在线定位异常问题。7. 后续的延展方向上面这些工作集中在B1I频点。如果设备已经支持B1C信号别忘了B1C和GPS L1的中心频率完全一致都是1575.42 MHz但调制方式不同一个是B1C的BOC或QMBOC一个是C/A码的BPSK-R它们在带宽和群延迟上的偏差特征又是另一套故事。频点相同并不意味着偏差为零码型不同、处理方式不同仍会出现分米级偏差。要深入扩展下一步可以做三件事第一在消音室或微波暗室里用模拟信号源扫频测出射频前端在不同频点上的群延迟曲线把硬件偏差建模成频率和温度的函数第二把伪距偏差和载波相位偏差联合估计借助载波相位的低噪声特性提高偏差标定精度第三在多系统组合PPP和RTK解算中引入实时估计的频间偏差看看对浮点解收敛时间有多少实际改善。8. 写在最后的个人经验做这类实测分析我最大的感触是伪距偏差里藏的东西比想象中多但也比想象中稳定。大部分时间里只要射频链路和天线没有大的设计缺陷GPS L1和北斗B1I的伪距偏差就是一个固定常数标定一次能用很久。真正让偏差失效的从来不是原理不明而是环境变化和外部因素干扰比如温度、AGC状态、天线周围多径还有某些软件后处理层面的“好心修正”。所以在自己的项目里我会建议测试人员先把固定偏差标定流程固化下来形成标准作业文档后续每次换硬件、改固件、动天线都跑一遍同样的流程。这样积累起来的偏差档案才是多系统高精度定位里最可靠的一笔资产。
