这篇由香港理工大学的研究团队发表的学术论文，题为《Multipath Parameter Estimation Using GNSS Receiver Doppler Measurements》，是一篇针对全球导航卫星系统（GNSS）在城市环境中面临的关键挑战——多径干扰——所进行的原创性研究。该研究首次深入探讨并建模了多径效应对接收机频率锁定环（FLL, Frequency-Locked Loop）和最终多普勒测量的影响，并在此基础上提出了一种新颖、低成本的多径信号参数估计与检测方法。该研究有望为提升低成本商业GNSS接收机在复杂城市环境中的定位精度与鲁棒性提供新的解决方案。
本研究的主要作者为香港理工大学的Jingxiaotao Fang、Bing Xu（通讯作者）、Xikun Liu和Li-Ta Hsu（资深会员）。研究被收录于《IEEE Transactions on Aerospace and Electronic Systems》期刊，虽然原文显示接收日期为0000年，但根据引用信息判断，其最终版本预计于2026年发表。该工作得到了国家自然科学基金和香港理工大学启动基金的支持。
研究的科学领域属于全球导航卫星系统（GNSS）信号处理与高精度定位，具体聚焦于多径误差的检测与缓解技术。研究动机源于一个显著的知识缺口：以往大量研究集中在多径对伪距（pseudorange）测量的影响上，并建立了广泛使用的多径误差包络（Multipath Error Envelope, MEE），但对多径在多普勒频率和伪距率（pseudorange rate）测量上影响的认知却存在矛盾。一方面，许多研究观测到在静态、动态及不同几何条件下，视距（Line-of-Sight, LOS）信号与反射信号之间存在不可忽视的多普勒频率差；另一方面，许多现有模型却假设此频率差可忽略不计。这种不一致性阻碍了基于多普勒测量的高级定位算法（如因子图优化、多普勒平滑、瞬时RTK等）的精度提升。因此，本研究旨在从理论和实验上澄清多径对多普勒测量的影响，并利用这一影响特征，开发一种基于标准接收机输出（无需修改硬件架构）的多径参数估计与缓解方法，特别服务于低成本的独立接收机用户。
本研究遵循从理论建模、算法提出、仿真验证到实际实验评估的完整流程。
1. 理论建模与特征分析： 首先，研究在作者先前工作的基础上，对接收机FLL在单视距信号加单反射信号的多径场景下进行了精确建模。该模型从接收到的复合信号出发，推导出理想提示（prompt）相关器的I/Q输出，进而得到复合信号的相位误差和频率鉴别器输出。通过一系列三角函数变换和近似，最终推导出接收机多普勒测量（f0(n)）的完整表达式（论文中式14）。该模型清晰地表明，多普勒测量误差（fmp(n)）是反射信号相对于视距信号的多普勒频移（fshift）、初始载波相位差（∆ϕshift）以及信号幅度比（α）的函数，并呈现出随时间振荡的特征。通过模拟分析，研究确认了该振荡的频率等于反射信号的多普勒频移，而其振幅则受频移和幅度比共同影响。此外，研究还绘制了伪距率误差包络图，并与经典的伪距误差包络进行对比，指出短延迟、强功率且伴有显著多普勒频移的反射信号会对多普勒测量产生更显著的影响，这为后续方法的应用场景提供了理论依据。
2. 算法设计与开发——多普勒非线性最小二乘估计器： 基于上述理论模型，本研究提出了核心创新方法：多普勒非线性最小二乘（D-NLS, Doppler-based Nonlinear Least Squares）估计器。算法的目标是利用一段观测窗口内的接收机多普勒测量序列，反演出四个未知信号参数：视距信号多普勒频率（f̃los）、反射信号的多普勒频移（f̃shift）、初始载波相位差（∆ϕ̃shift）以及幅度比（α̃）。 - 成本函数与优化： 算法通过最小化实际接收到的多普勒测量值与基于模型（式19）估计出的多普勒序列之间的残差二范数来求解参数。这是一个有界约束的非线性最小二乘优化问题，采用MATLAB的lsqnonlin求解器，结合信任域反射（trust-region-reflective）算法进行迭代优化。 - 初始化策略： 为解决非线性优化对初始值敏感的问题，研究设计了基于多普勒振荡特征的初始化方法：f̃los初始化为观测窗口内的多普勒测量均值；f̃shift初始化为观测到的振荡频率的倒数；α̃则利用推导出的多普勒振荡极值公式（式15）进行反推。 - 多径检测流程： 完整的算法流程如Algorithm 2所示。对每个卫星通道的每个估计窗口，先进行参数初始化，然后运行D-NLS优化。检测多径的依据是：如果优化后重建出的估计多普勒序列的振荡范围（最大值减最小值）超过预设阈值（研究中设为4 Hz），则判定该窗口存在多径。检测到多径后，可利用估计出的f̂los对接收机多普勒测量进行校正，并将更新后的测量值用于导航解算。
3. 仿真评估与性能边界分析： 在应用于真实数据前，研究通过蒙特卡洛仿真系统地评估了D-NLS估计器的性能，并推导了克拉美 罗下界（CRLB, Cramér–Rao Lower Bound）作为理论精度基准。评估主要关注两个关键因素：信噪比（C/N0）和测量样本大小（即估计窗口时长）。 - 仿真在固定参数（flos=0 Hz, fshift=15 Hz, ∆ϕshift=0°, α=2）下进行，通过添加不同水平的噪声（模拟不同C/N0）或改变样本数量，运行5000次蒙特卡洛试验计算参数估计的均方根误差（RMSE）。 - 结果表明，在C/N0较高（>30 dB-Hz）且样本大小适中（约300毫秒）时，估计器性能良好，对f̃los的估计接近CRLB。f̃shift、∆ϕ̃shift和α̃的估计误差相对更大，但对多径检测而言精度已足够。样本过大（>300ms）会因模型非线性增强和累积噪声过多而导致性能下降。
4. 真实场景实验与验证： 研究设计了行人和车载两种典型的城市动态场景实验，在香港的城区进行。 - 数据采集系统： 使用LabSat 3W系统记录GPS L1中频（IF）数据供后处理软件接收机（SDR）使用，同时使用NovAtel SPAN-CPT系统（通过分路器连接同一天线）提供厘米级精度的地面真值轨迹。 - 基准工具——光线追踪： 作为性能评估的基准，研究采用了光线追踪（Ray-Tracing）技术。该技术结合3D建筑模型、卫星星历和用户地面真值数据，计算出从卫星到接收机最可能的信号传播路径（包括反射），从而模拟出视距和反射信号的“真实”多普勒频率、路径延迟等参数。 - 实验设置： 场景1（行人）： 测试者垂直于建筑物表面往复行走，产生几何关系明确、反射信号多普勒频移变化的可控多径环境，重点关注PRN 22卫星。场景2（车载）： 车辆在密集的城市峡谷中行驶，卫星-接收机几何关系快速变化，多径短暂而剧烈，重点关注PRN 15卫星。SDR的PLL环路带宽分别设置为20 Hz（行人）和40 Hz（车辆），D-NLS处理的样本窗口统一为300 ms。
1. D-NLS估计器在实际数据中的表现： - 场景1（行人）： 对于受多径影响的PRN 22，接收机多普勒测量清晰地显示出理论预测的高频振荡特征。D-NLS估计器成功地重建了该振荡曲线，估计出的视距和反射信号多普勒频率与光线追踪模拟的“真实”值在整体趋势和多数时段内吻合良好。计算单差多普勒（以LOS卫星PRN 20为参考）后的RMSE，视距多普勒估计约为2 Hz，反射多普勒估计约为6 Hz，验证了方法的有效性。对于纯视距卫星PRN 17，估计器输出平稳，视距与反射多普勒估计值几乎重合，符合预期。 - 场景2（车载）： 在高速动态环境下，PRN 15的多普勒测量表现出更剧烈、周期更短的振荡（约45 ms周期，对应25 Hz频移）。D-NLS在光线追踪判定为多径的时段（约31-45秒）内，能够估计出与模拟值一致的反射信号多普勒。在非多径时段，两者出现不一致，分析表明这可能源于光线追踪模型的局限性（如建筑模型精度不足、未考虑漫反射等）、D-NLS将动态跟踪噪声误判为多径，或真实环境中存在弱多径而光线追踪未检出。
2. 与现有技术的对比分析： - 与MEDLL对比： 研究将D-NLS与经典的多径估计延迟锁定环（MEDLL）方法在行人场景下进行了直接比较。为公平起见，实现了高性能版本的MEDLL（使用161个相关器）。结果表明：在用户静态或低速阶段，MEDLL能稳定估计出合理的反射信号幅度比（2-4，相当于反射信号衰减6-12 dB）和码延迟（约0.2 chips），而D-NLS因缺乏显著振荡特征而性能有限。在用户动态阶段，情况反转：D-NLS能给出稳定、合理的幅度比估计（2-5），而MEDLL的码延迟估计出现剧烈波动，表明其在高动态噪声环境下失效。这凸显了两种方法优势的互补性：MEDLL基于相关器，在静态/准静态下分辨率高；D-NLS基于多普勒动力学特征，在动态场景下更鲁棒。 - 与SQM指标对比： 研究还将D-NLS的多径检测率与传统信号质量监测（SQM）指标（Delta测试和Ratio测试）进行了对比。在所有观测卫星中，D-NLS对多径卫星PRN 22的检测率约为50.8%（若只考虑动态时段，升至64%），同时对其他视距卫星保持了极低的虚警率（普遍低于5%）。相比之下，Ratio测试虽然对PRN 22也有44%的检测率，但对视距卫星PRN 17的虚警率高达93.6%，因其对相关函数（ACF）的任何畸变都过于敏感。Delta测试性能居中。这表明D-NLS在保持高检测灵敏度的同时，能更可靠地区分多径和噪声/动态波动。
3. 计算效率： 研究记录了算法处理时间。在所述硬件配置下，处理1秒的GNSS数据（6颗卫星），D-NLS平均耗时约1.38秒，而高性能MEDLL则需约440秒。这体现了D-NLS在计算效率上的显著优势，更有利于实时或近实时应用。
本研究系统地建立并验证了多径干扰下GNSS接收机多普勒测量的理论模型，揭示其特征性的振荡现象，并创新性地提出了基于此现象的D-NLS多径参数估计与检测方法。主要结论如下： 1. 多径效应确实会在接收机多普勒测量中引入不可忽视的、具有特定振荡模式的误差，该模式由反射信号相对于视距信号的参数（频移、相位、幅度）决定。 2. 所提出的D-NLS方法能够有效地从标准接收机多普勒输出中估计出这些信号参数，实现多径检测，并在动态用户场景下展现出优于传统MEDLL和SQM方法的性能，兼具较高的检测率、较低的虚警率和更好的计算效率。 3. 通过与光线追踪基准的比较，验证了该方法估计信号参数（尤其是多普勒频率）的准确性。
研究的价值体现在： - 科学价值： 填补了多径对GNSS多普勒测量影响机理研究的知识空白，提供了严谨的理论分析和模型，解决了该领域长期存在的观点矛盾。 - 应用价值： 提出了一种无需修改接收机硬件、仅利用其标准多普勒输出的低成本多径缓解方案，特别适用于资源受限的低成本商业接收机。校正后的多普勒测量可直接用于提升各类多普勒辅助的定位、测速和模糊度固定算法的性能。 - 工程价值： 研究指出D-NLS（对动态多径敏感）和MEDLL（对静态多径敏感）的互补特性，为未来开发自适应、混合型的多径缓解系统指明了方向。
研究在讨论中也坦诚指出了当前方法的局限性，并为未来工作提供了清晰的方向： - 局限性： 目前模型和算法主要针对单一主导反射信号且视径信号占优的场景；在极高噪声或用户动态极其剧烈时，性能会下降；未在模型中显式考虑噪声分布。 - 未来展望： 作者建议可探索更鲁棒的估计算法或机器学习技术以提升抗噪能力和处理多反射场景的能力；研究将校正后的多普勒频率反馈至接收机跟踪环路以直接改善跟踪性能；开发结合D-NLS和MEDLL优势的混合型自适应多径缓解系统。这些都为后续研究提供了富有潜力的切入点。