★ FEATURED ARTICLE 简介《GPS三频组合观测量的特征及应用研究》是国防科学技术大学学者吕志成等撰写的学术论文压缩包内为单篇PDF全文大小295KB适合GPS定位系统开发、卫星导航算法研究及精密测量专业人员作为参考文献。文章从GPS现代化引入L2C、L5民用信号出发介绍双差载波相位观测方程及卫星轨道误差、电离层时延、对流层时延、观测噪声等影响因素。随后推导三频组合观测值的一般表达式分析组合波长、模糊度及误差放大特性并重点研究电离层无关线性组合的构建方法与性能。针对多种特殊线性组合形式文中论证合理选择组合有助于改善整周模糊度解算成功率提升粗差探测与周跳检测修复的有效性。最后总结三频组合观测值在电离层误差改正、周跳探测与修复、整周模糊度解算、精密定位中的应用前景为高精度定位算法设计、系统优化与研究提供理论支撑与实践指导。已有77人学习下载。1. GPS定位系统中的三频组合这篇论文为什么值得细读做GPS定位系统开发的人在翻阅参考文献时大概率会撞上三频组合这个话题。以前处理长基线或高动态数据单频靠电离层模型硬扛双频消一阶项但残差仍然存在三频信号L1、L2C、L5全部可用之后载波相位观测值通过线性组合能同时调节波长、电离层放大倍数和噪声放大倍数——这是双频时代没有的自由度。这篇论文把三频组合观测量的通用表达式、误差影响因素、电离层无关组合的构建方法以及EWL/WL/ML/NL几种特殊组合的特征值完整梳理了一遍对专业指导周跳探测、整周模糊度解算和精密定位的数据处理策略很有参考价值。下文按论文的推导顺序拆开讲把参数怎么算、组合怎么选、坑在哪里一次说清。2. 双差载波相位观测方程为什么三频研究要从这一步开始2.1 双差算子到底消掉了什么观测方程逐项拆解GPS精密定位很少直接用原始载波相位观测值常规做法是先做双差。双差算子的作用是把接收机钟差和卫星钟差一次性消掉同时大幅削弱卫星轨道误差、电离层时延和对流层时延这些空间相关误差。论文给出的双差载波相位观测方程为∇Δφ_i (1/λ_i)(∇ΔR ∇ΔdR ∇ΔT) − (f_i²/f₁²)·(∇ΔI₁/λ_i) ∇ΔN_i ∇Δε_i/λ_i其中∇Δ表示双差算子φ_i是第i个频率的载波相位观测值单位周λ_i是对应波长单位mR是接收机与卫星间的几何距离dR是卫星轨道误差I₁是L1频率上的电离层时延单位mT是对流层时延N_i是整周模糊度ε_i包含接收机噪声和多径效应。这个式子的关键信息有三层。第一层卫星钟差和接收机钟差在双差过程中被完全消除所以方程里根本没有钟差项。第二层残余误差只剩轨道误差残差、电离层时延残差、对流层时延残差和观测噪声四类。第三层电离层时延在式子里乘了频率比f_i²/f₁²这说明电离层误差与频率的平方成反比——正是这个频率相关性让三频线性组合有了操作空间。论文特别强调双差观测方程是后面所有线性组合的基础。三频组合不是在原始单差观测值上做而是在双差观测值上做。这一点容易被忽略实际处理数据时如果输入的不是双差观测量后面的组合系数和误差放大因子全都不成立。2.2 差分消不掉的残差多径与接收机噪声的随机模型约定双差能消掉钟差、大幅削弱空间相关误差但有两类误差它无能为力多径效应和接收机噪声。原因是这两者在时间和空间上都不具备明显的相关性差分不仅消不掉在某些几何条件下还会被放大。论文的处理建议是把它们放到随机模型里考虑也就是给不同高度角、不同频率的观测值赋予不同的方差用加权的方式控制它们对解算的影响。轨道误差残差和电离层残差则要分场景看。短基线情况下两个测站看到的卫星轨道误差基本相同双差后残差很小可以忽略但长基线或者高动态场景下轨道残差和电离层残差会成为主导误差源。单频系统只能靠模型改正双频系统可以构建电离层无关组合消掉一阶项但消完一阶项之后残差依然存在。三频的意义在于多了一个频率后可以用三个观测值构造出同时满足长波长、低噪声、电离层无关等更苛刻条件的组合这是双频做不到的。下表总结了双差后各项残余误差的行为和处理方式。误差类型双差后表现典型处理手段卫星钟差、接收机钟差完全消除无需处理轨道误差残差短基线可忽略长基线显著精密星历或作为参数估计电离层时延残差短基线小长基线/高动态大三频组合消除或模型改正对流层时延残差与频率无关无法靠组合消除模型改正加参数估计多径效应时空无相关性差分无效随机模型加权、抗多径天线接收机噪声无法通过差分削弱随机模型加权、选低噪声接收机对流层延迟不依赖频率这是三频组合解决不了的问题只能靠模型加估计。所以三频组合针对的核心目标是电离层残差和观测噪声这两个频率相关项。3. 三频线性组合的一般形式波长、模糊度与三个放大因子3.1 组合观测量的数学定义系数必须取整才能保住整数模糊度三频组合的思想很直接把三个频率的载波相位观测值各乘一个系数后相加得到一个新的组合观测量。论文给出的一般表达式为φ_c k₁·φ₁ k₂·φ₂ k₃·φ₃对应的组合模糊度为N_c k₁·N₁ k₂·N₂ k₃·N₃要让组合模糊度仍然保持整数特性k₁、k₂、k₃在实际使用中必须取整数。这一点是后面所有模糊度解算方法能够成立的根基——如果系数带小数组合模糊度就不是整数取整法直接失效。组合观测量的波长由组合频率决定λ_c 1 / (k₁/λ₁ k₂/λ₂ k₃/λ₃)代入GPS三个民用频率L11575.42MHz、L21227.60MHz、L51176.45MHz就能算出每个组合对应的波长。比如超宽巷组合系数(0,1,-1)组合后的波长为5.86m远大于任何一个原始频率的波长。原理在于组合频率k₁f₁k₂f₂k₃f₃变得很小波长自然被拉长。这里有一个重要的因果关系组合波长的分母是k_i/λ_i之和这个和越小波长越长而几何相关误差的放大因子又正比于这个和。所以长波长组合天然具备压低几何误差的特性。论文将这一点在误差分析中单独展开。3.2 三个放大因子怎么算几何、电离层、噪声各有各的标尺把三频组合表达式代入双差观测方程后可以得到三个关键的误差放大因子这是选择组合系数时必须先算清楚的三组数。几何相关误差放大因子dG (k₁λ₁ k₂λ₂ k₃λ₃) / λ₁它表示组合观测量的轨道误差、对流层误差和几何距离误差相对于L1载波的放大倍数。当dG等于0时该组合对几何相关误差完全不敏感称为几何无关组合。电离层误差放大因子分两种表示方式。以米为单位时记为α_I以周为单位时记为β_I两者的差别只是波长换算。当α_I和β_I接近0时组合观测值基本不受电离层时延影响这就是电离层无关组合。电离层无关线性组合在应用中意义最大论文用了单独一章来讨论它的构建方法。测量噪声放大因子同样分m和周两种单位。假设各频率载波相位观测误差相互独立且以米为单位的观测噪声均方差等于相应波长的a倍a取经验值2.5%那么组合噪声放大因子为dσ(m) √(k₁²k₂²k₃²) · λ_c / λ₁dσ(周) √(k₁²k₂²k₃²)注意周单位噪声放大因子只与系数平方和有关与波长无关。这个特性让某些组合虽然波长很长但周单位噪声并不放大——这正是超宽巷组合能被用于周跳探测的根本原因。下表把三个放大因子的计算方式和关注点放在一起。放大因子计算公式理想值主要影响几何放大dG(k₁λ₁k₂λ₂k₃λ₃)/λ₁0几何无关轨道误差、对流层误差、几何距离电离层放大mα_I由频率比与系数导出0电离层无关电离层时延及残差噪声放大m√(k₁²k₂²k₃²)·λ_c/λ₁越小越好模糊度固定成功率与定位精度噪声放大周√(k₁²k₂²k₃²)越小越好周跳探测灵敏度三频组合选系数的过程本质上就是在这个三维空间里找折中点想要波长长几何误差就小但m单位噪声往往被放大想要噪声小波长通常被压短模糊度又不好固定。3.3 噪声约束从哪来Jung阈值与k₁²k₂²k₃²100选组合系数不是随心所欲的论文引用了Jung给出的一个实用门限为了保证取整法求解整周模糊度有效组合观测量的噪声必须小于组合波长的一半即σ_c λ_c/2。把噪声放大公式代入这个不等式整理后得到√(k₁²k₂²k₃²) 1/(2a)按经验值a2.5%即原始观测噪声约为波长的2.5%计算最终约束为k₁² k₂² k₃² 100这个约束看起来简单实际筛选组合时非常有用。后面构造电离层无关组合时它就是一道硬性过滤门槛满足电离层无关的组合很多但噪声超标的组合直接淘汰。做过一阵子三频数据处理的人应该会有同感这个阈值是少有的不用玄学就能定死的参数。4. 电离层无关组合构建一个带不等式约束的整数搜索问题4.1 电离层无关的约束条件一个方程把三个频率绑在一起构建电离层无关组合的目标是让组合观测值的电离层误差为0。根据前面的误差分析这意味着组合后的电离层放大因子必须满足k₁/f₁ k₂/f₂ k₃/f₃ 0这个式子把三个频率和三个系数绑在了一个方程里。三个未知数一个方程解不唯一需要加约束才能得到具体系数。论文给出的约束条件还包括系数必须是整数且满足噪声门限k₁²k₂²k₃²100。到这里问题变成一个带不等式约束的非满秩整数方程组求解。手算不现实论文的处理办法是枚举加取整。4.2 浮点解与定点解整数搜索的具体筛选步骤我按论文思路整理出一套可以直接照做的筛选步骤第一步固定k₁和k₃在[-10,10]区间内取整数k₂先当作实数处理。之所以把k₂留出来是因为三个频率里L2居中解出来的浮点值更容易落在合理区间。第二步由电离层无关方程解出k₂的浮点解k₂_float −f₂·(k₁/f₁ k₃/f₃)第三步对k₂_float做四舍五入得到定点解k₂_int。浮点解与定点解的偏差越小组合后残留的电离层误差越小这一步是筛选质量的关键。第四步用噪声约束k₁²k₂²k₃²100过滤候选组合。第五步计算每个候选组合的波长λ_c、电离层放大因子和噪声放大因子按浮点解与定点解偏差从大到小排序取前几组作为候选。用(1,-6,5)这个组合来验证这套流程。取k₁1k₃5算出的浮点解为k₂_float −1227.60 × (1/1575.42 5/1176.45) ≈ −5.9966四舍五入后k₂-6浮点解与定点解偏差只有0.0034非常接近理想值所以这个组合的电离层残差极小。同时k₁²k₂²k₃²1362562顺利通过噪声门限。这个流程的妙处在于把连续的浮点解转化为整数解再用偏差大小反推组合质量。实际写程序时这个搜索在毫秒级就能跑完比盲目穷举系数高效得多。4.3 五组候选组合怎么看波长、噪声、电离层放大系数的取舍论文按上述方法给出了五组满足电离层无关条件的候选组合及其特征值。这些组合有一个共同点除了(1,-6,5)之外其他组合的波长都小于任意单一载波的波长。短波长对整周模糊度解算不利但它们的电离层误差相对L1降低了大概一到两个数量级在基线解算中可以忽略不计。(1,-6,5)是三组中比较特殊的一个波长约3.22m远大于L1的19cm。波长长意味着几何相关误差小电离层残差也小这两个特性让它在理论上很好看。但代价同样明显——以米为单位的观测噪声放大了约两个数量级。如果拿它直接做精密定位坐标解会非常吵。这就是典型的三频组合困境没有一个组合是样样都好的。电离层无关组合适合作为中间观测量用来约束电离层误差或辅助模糊度搜索但不适合直接作为最终定位观测量。实际工程中通常是把它们嵌进序贯解算流程而不是单独拎出来用。5. 特殊组合选型避坑EWL/WL/ML/NL的参数权衡与常见翻车点5.1 四种特殊组合的参数表一个表看懂长波与低噪声的矛盾论文在电离层无关组合之外还专门分析了另一种选系数的策略不是追求某个误差为零而是选择能突出某一项特性的组合。四种有代表性的组合是超宽巷EWL、宽巷WL、中巷ML和窄巷NL参数对比如下。组合名称系数(k₁,k₂,k₃)组合波长m单位噪声放大周单位噪声放大特性EWL超宽巷(0,1,-1)586.1cm约43.6倍√2≈1.41倍波长极长电离层小m噪声大WL宽巷(1,-1,0)86.2cm约6.4倍√2≈1.41倍波长较长各项指标均衡ML中巷(1,0,-1)75.1cm约5.6倍√2≈1.41倍略短于宽巷噪声略低NL窄巷(1,1,0)10.7cm约0.8倍√2≈1.41倍波长短电离层大m噪声最小L1单频(1,0,0)19.0cm1倍基准1倍常规处理基准有一个细节值得注意这四种组合的周单位噪声放大因子都是√2因为它们的系数平方和都等于2。也就是说以周为单位看噪声这些组合几乎没有代价但以m为单位看EWL的噪声被放大了40多倍。很多人在这一步把单位搞混直接导致后面一系列误判。EWL的波长586.1cm最夸张接近6m电离层误差和几何误差都非常小这些特性对整周模糊度解算极为有利。但它的m单位噪声太大不适合直接用于精密定位。WL和ML波长在1m上下电离层和噪声都适中更适合作为模糊度解算的中间步骤。NL波长只有10.7cm电离层和几何误差大模糊度不好固定但m单位噪声小一旦模糊度固定它的定位精度最高。5.2 组合选型的一般原则序贯组合先长波后短波的原因观察EWL、WL、ML三个组合的特征值变化会发现一条清晰的规律波长从5.86m递减到0.75mm单位噪声放大系数从43.56一路降到5.58周单位噪声保持√2不变电离层时延略有增加但幅度有限。这条规律直接给出了三频模糊度解算的序贯思路先用长波组合固定模糊度再用短波组合提高精度。具体来说先处理EWL或WL组合把它们的模糊度确定下来然后用已经固定的模糊度作为约束逐步过渡到波长更短、噪声更小的组合最终得到L1或窄巷的模糊度。每一步都在上一步的基础上缩小搜索范围成功率比直接解L1单频模糊度高得多。论文中提到的TCAR、CIR、MCAR方法本质上都是这个思路的变体。这个序贯思想也解释了为什么(1,-6,5)这类电离层无关组合虽然不能直接用但依然有研究价值——它可以在序贯链路的某个环节充当电离层误差的约束条件帮助其他组合更快地收敛。5.3 避坑记录三频组合使用中最容易被绕进去的五个问题这一章列出几个实际处理三频数据时经常遇到的坑每一条都是从现象反推原因再给解决路径。提示下面的问题在论文的理论框架下都有明确解释遇到时先回头看一眼对应的放大因子比盲目调参数有效。现象一直接用EWL超宽巷组合输出定位结果坐标序列抖动明显甚至比单频还差。原因EWL的m单位噪声放大超过40倍。虽然以周为单位看噪声很小但换算到米后噪声被剧烈放大直接进入定位解会导致坐标噪声暴涨。解决EWL只做周跳探测、粗差定位和模糊度初值估计最终定位交给L1或窄巷组合。把它当作辅助工具而不是定位工具。现象二用(1,-6,5)做电离层无关组合在信号质量一般的场景下模糊度固定失败率很高。原因该组合m单位噪声被放大两个数量级噪声门限守不住。信号质量好时还能靠长波长优势压住一旦信号变差噪声立刻成为主导。解决只在数据质量好的静态场景使用(1,-6,5)动态或遮挡严重的场景改用WL或ML组合先固定宽巷模糊度再逐级过渡。现象三组合出来的模糊度残差总在0.3到0.5周之间游走取整法成功率很低。原因组合系数不是整数或者用了浮点系数近似导致组合模糊度不具备整数特性。另一个可能是浮点解与定点解偏差过大电离层残差超出了允许范围。解决检查k₁、k₂、k₃是否全部为整数用4.2节的浮点解与定点解偏差评估电离层残差偏差大于万分之一量级就先排除。现象四EWL周跳探测不报警WL却频繁报警分不清是真周跳还是电离层扰动。原因周跳探测和电离层扰动的响应在不同组合里表现不同。EWL对电离层不敏感周跳信号非常干净WL对电离层残差的灵敏度更高电离层活跃时容易出现虚警。解决两个组合同时做交叉验证。EWL报警可以确认周跳WL报警需要先看电离层活跃程度必要时用三频组合把电离层项单独估出来再判断。现象五写程序穷举组合系数搜出几百组算出来要么波长太短、要么噪声超标、要么电离层残差太大。原因只枚举了k₁和k₃的整数格点没有按浮点解约束k₂也没有计算真实波长、电离层放大因子和噪声放大因子做二次过滤。解决按4.2节的完整流程走一遍——固定k₁/k₃、解k₂浮点值、取整、算三个放大因子、按噪声门限过滤、按偏差排序。十几行代码就能从几千组候选里筛出最合理的五组不需要人工一条条试。6. 三频组合在解算中的落地用法从周跳修复到模糊度序贯解算6.1 周跳探测、模糊度解算与验证方法三频组合在数据处理里的三个主要应用场景是电离层误差修正、周跳探测与修复、整周模糊度解算。论文对这三个场景都有涉及落地时最实用的是周跳探测和序贯模糊度固定。周跳探测的推荐路径是先用EWL组合构建周跳检测量。因为EWL波长接近6m周单位噪声只有√2倍电离层误差极小周跳发生时会在组合残差序列里出现明显的台阶状跳变肉眼都能看出来。检测到周跳后用WL组合做二次确认两个组合同时跳变的位置基本可以判定为真实周跳。修复时先恢复EWL模糊度再逐级恢复WL和L1模糊度修复后的数据再进入解算流程。这一段我自己跑过的流程和论文描述完全对得上三频在周跳处理上的优势确实是实打实的。模糊度解算则按序贯链路推进第一步用EWL固定超宽巷模糊度第二步回代固定WL宽巷模糊度第三步通过ML过渡最后固定L1或NL窄巷模糊度。每一步用上一个组合已经固定的模糊度作为已知约束把搜索空间压缩到极小。验证方法建议拿一组静态三频数据分别用单频、双频和三频流程处理对比模糊度固定成功率和高程精度再画一张三条组合残差叠加图观察EWL那一条是否明显比WL和单频干净。这个对比做完对三频组合的价值会有直观感受。这里有一条某开发者用血泪换来的习惯每次选定一组组合系数先列出λ_c、α_I、dσ(m)、dσ(周)四个数再决定这个组合用于周跳、模糊度还是最终定位。从那以后我每次处理三频数据都强制走一遍这个检查基本没有再因为选错组合翻过车。这份论文PDF里表1和表2的原始数据适合做三频数据处理时当作案头参考按需取用即可希望帮到你。本文还有配套的精品资源点击获取
阅读完成 · 觉得有帮助？
