基于多通道SQM指标联合的矢量接收机多径干扰检测方法.pdf

DOI：10.12265/j.gnss.2023037基于多通道 SQM 指标联合的矢量接收机多径干扰检测方法贾琼琼，朱传胜(中国民航大学 天津市智能信号与图像处理重点实验室,天津 300300)摘要：矢量接收机利用各跟踪通道间信息共享达到通道间相互辅助的作用，因此相对于传统的标量接收机能在复杂环境下提供更好的性能.然而在城市峡谷等复杂环境下，多径传输会严重限制矢量接收机的性能，若能及时发现多径并加以处理则能消除多径干扰的影响.本文从多径对矢量接收机的影响机理出发，研究基于信号质量监测(SQM)的多径检测问题.在分析现有的针对标量接收机信号 SQM 指标的基础上，结合矢量接收机通道耦合的特点，提出一种联合多通道 SQM 指标的矢量接收机多径检测方法.最后，通过仿真实验验证了所提方法能够有效地检测与直达信号幅度比小于 0.5、相对码延时在 0.30.6 码片的多径干扰.关键词：全球导航卫星系统(GNSS)；多径干扰；矢量跟踪；信号质量监测；多通道联合检测中图分类号：P228.4文献标志码：A文章编号：1008-9268（2023）03-0110-100引言全球导航卫星系统(GNSS)能为其覆盖范围内的用户提供高精度的位置和时间信息，因此被广泛地应用于多个领域.矢量接收机通过充分利用各跟踪通道间信息共享能够达到通道间相互辅助的作用，因此在复杂环境下具有比传统标量接收机更优的性能.在城市峡谷等环境下，建筑物或其他障碍物会对卫星信号产生反射或散射，这些信号和直达信号一起进入接收机的现象称为多径传播.研究表明，多径干扰引起的伪距误差可以达到米级甚至百米级1，这足以影响系统的可靠性和定位精度.尽管在复杂环境下，矢量接收机较之标量接收机具有性能优势，但是，多径仍然会对矢量接收机引入跟踪误差，导致其导航解偏离真实值.因此，如果能够及时检测到多径并加以处理，则有望降低多径对矢量接收机的影响，从而提高其在城市峡谷环境下的鲁棒性.信号质量监测(SQM)技术最初是用来监测卫星异常的，它通过对接收机跟踪通道的相关值等输出进行监测以判断是否有干扰出现.SQM 已被广泛用于检测标量接收机中的多径干扰2-8.文献 3-4 利用早码、即时码和晚码相关器的输出定义了对称和非对称 Delta 和 Ratio 指标，以监测卫星信号相关峰值中是否存在由多径引起的畸变.文献 5 利用多径干扰下跟踪环路相关曲线两侧不对称，提出了斜率不对称(SAM)指标来监测卫星信号是否受到多径影响.文献 6 通过定义 DoubleDelta 指标研究了 SQM 技术在 GNSS 多径检测和缓解中的性能，提出了SQM 变化包络，并将其与相应的多径误差包络进行比较，以分析 SQM 指标对 GNSS 多径检测和缓解的敏感性和有效性.考虑到多径也会影响卫星信号的载波相位，Mubarak 等7-8提出了早-晚码相位(ELP)技术，利用早-晚相关器输出之间的相位差来检测多径干扰.目前，基于 SQM 实现矢量接收机多径干扰检测的研究还较少.文献 9 通过推导模型描述了矢量接收机下码和载波跟踪误差与多径延迟、相位和衰落频率的关系，并结合广义极大似然比检验和多相关器技术提出了一种多径检测方法，验证了其检测多径的有效性，但是没有与常用的 SQM 指标进行对比分析多径干扰下的检测性能；文献 10 将早减晚相关器输收稿日期：2023-03-08资助项目：中国民航大学民航航班广域监视与安全管控技术重点实验室开放基金(202202)通信作者：贾琼琼E-mail：第 48卷第3期全球定位系统Vol.48，No.32023年6月GNSS World of ChinaJune，2023C/N0出和估计器结合用作信号质量指标来检测多径，通过排除多径污染卫星提高了定位精度以及矢量跟踪性能.然而，现有的研究主要是基于矢量接收机各通道相关器输出进行检测，没有考虑不同 SQM 指标在标量接收机和矢量接收机的区别.事实上，对于矢量接收机，多径干扰存在时不仅受干扰信号 SQM指标出现异常，而且正常信号的指标值也会受到影响.基于上述原因，本文提出基于多通道 SQM 指标联合的矢量接收机多径干扰检测方法，该方法利用多径影响在矢量接收机各跟踪通道之间传递的特点，综合各通道的 SQM 指标以使多径的影响得以积累，从而更易于发现多径干扰；在此基础上，再利用 DoubleDelta 指标对实际接收到多径的通道加以识别.本文方法较之现有方法能更快速地发现多径，这有利于及时采取相应的措施对其进行处理，从而降低甚至消除多径对矢量接收机的影响.本文最后通过仿真实验验证了所提方法的有效性.本文其他章节安排如下：第 1 节给出了多径信号模型，并讨论了矢量跟踪结构和标量跟踪结构的区别，进而得到了多径对矢量接收机影响的特点；第2 节在分析常用 SQM 指标在矢量接收机中多径监测特点的基础上，提出了基于多通道 SQM 指标联合的矢量接收机多径检测方法；第 3 节通过仿真实验验证了本文所提多径监测方法的有效性；第 4 节进行了总结.1信号模型及问题描述 1.1 信号模型ni本文以 GPSL1C/A 信号为例来分析多径的影响，所得结论可推广至其他卫星信号.假设当前时刻可见卫星数为，其中第颗卫星信号还存在一条反射路径的回波，此时包含多径干扰的中频信号可以表示为10s(t)=nj=1AjCj(tj)cos(2fD,jt+j)+MPiAiCi(tiMPi)cos(2fD,it+i+MPi)+(t).(1)AjjfD,jjjCj()jMPiiMPiMPi式中：、和分别为第颗卫星信号的幅度、码时延、多普勒频率和载波相位；为第颗卫星信号的伪随机噪声码(PRN)；为第颗卫星信号反射路径的幅度相对于直达信号幅度的衰减比；和分别为该路径信号相对于直达路径的传输时延和(t)20载波相位，由于反射点一般在接收机附近，对信号多普勒频率的影响较小，因此可假设反射路径与直达路径的多普勒频率相同11；为均值为零方差为的高斯白噪声.GNSS 接收机码跟踪环路通常配置三个相关器，分别为早码(Early,E)、即时码(Prompt,P)和晚码(Late,L)相关器，其目的是通过比较三个相关器的输出来调整本地参考信号的码相位控制量.当接收信号为式(1)中的包含反射路径的多径传输时，三个相关器的输出分别为：IE=AjR(j+d/2)cos(j)+MPiAiR(MPi+d/2)cos(MPi)QE=AjR(j+d/2)sin(j)+MPiAiR(MPi+d/2)sin(MPi);(2)IP=AjR(j)cos(j)+MPiAiR(MPi)cos(MPi)QP=AjR(j)sin(j)+MPiAiR(MPi)sin(MPi);(3)IL=AjR(jd/2)cos(j)+MPiAiR(MPid/2)cos(MPi)QL=AjR(jd/2)sin(j)+MPiAiR(MPid/2)sin(MPi).(4)IEIPILQEQPQLjjMPiMPidR()式中：、分别为同相支路的早码，即时码和晚码相关值；、分别为正交支路的早码，即时码和晚码相关值；和分别为直达信号与本地信号之间的码相位差和载波相位差；和分别为多径信号与本地信号之间的码相位差和载波相位差；为早码和晚码相关器之间的间隔；为 PRN码相关函数.R()=1|Tc,|Tc0,其它,(5)Tc为码宽.MPi=0对于不包含反射路径的卫星信号来说，其自相关函数的主峰值范围内可以看成对称三角形，如图 1 中黑色线条所示.对于包含反射路径的卫星信号来说，其相关函数会发生畸变，如图 1 中绿色线条所示，表示建设性()多径发生时的相关函数.第3期贾琼琼，等：基于多通道 SQM 指标联合的矢量接收机多径干扰检测方法111直达信号合成信号反射信号码延时/码片ACFPLE101MPi图 1 多径干扰下的自相关函数 1.2 矢量接收机原理简介传统的标量接收机各跟踪通道之间相互独立，如图 2 所示，因此当多径干扰出现时，只有存在多径干扰的跟踪通道的相关函数会发生畸变，导致相关器的输出产生偏离，而对于不存在多径干扰的通道相关器的输出则不受影响.对于如图 3 所示的基于矢量延迟锁定环(VDLL)的矢量接收机12来说，由于其各跟踪通道之间相互耦合，所有跟踪通道的本地信号码相位都是由反馈的导航解来控制的.由于导航解是各跟踪通道综合作用的结果，那么利用导航解来控制各通道本地信号的码相位就意味着各跟踪通道会受到其他跟踪通道的影响.为了讨论多径对矢量接收机的影响，下面首先简单介绍矢量接收机的工作原理.跟踪环路中频信号相关器码/载波NCO码/载波鉴别器环路滤波器通道1通道NPVT.导航处理器图 2 标量接收机跟踪结构中频信号跟踪环路相关器码NCO载波NCO通道1通道N.码鉴别器载波鉴别器环路滤波器EKF图 3 矢量接收机跟踪结构X=u,v,ctb,ctdTu=(ux,uy,uz)v=(vx,vy,vz)ctbtd定义为接收机的状态矢量，其中和分别为接收机在地心地固坐标系(ECEF)中的三维位置和三维速度；为光速；为接收机钟差；为钟漂.定义接收机状态误差矢量为X=u,v,ctb,ctdT.(6)u=(ux,uy,uz)v=(vx,vy,vz)tbtd式中：、分别为ECEF下的三维位置误差和三维速度误差；和分别为接收机的钟差误差和钟漂误差.矢量接收机一般选用扩展卡尔曼滤波器(EKF)作为导航滤波器，其本质是求解下列问题：Xk+1=Fk,k+1Xk+wkZk+1=Hk+1Xk+1+k.(7)Fk,k+1wkHk+1kZk+1式中：为状态预测矩阵；为过程噪声矢量；为观测矩阵；为观测噪声矢量13.矢量接收机的观测量表示为Zk+1=1,2,n,1,2,n.(8)j,k+1其中，前半部分为各通道的伪距误差观测量，其与码鉴别器输出的关系可以表示为j,k+1=j,k+1cfCA.(9)fCA j,k+1式中：为 PRN 的码频率(对 GPSL1 信号来说为1.023MHz).对于本文的 VDLL 矢量接收机来说，其伪距率误差通过下式计算：j,k+1=fd,j,k+1cfL1(vj,k+1vk+1)aTj,k+1ctd,j.(10)fd,j,k+1jfL1vj,k+1jaj,k+1=(ax,j,k+1,ay,j,k+1,az,j,k+1)j式中：为第颗卫星的多普勒频移频率；为载波频率；为第颗卫星的速度；为接收机和第颗卫星间的方向矢量.EKF 求解式(7)分为两个步骤：第一步，根据当前时刻的导航解预测下一历元的接收机状态，即Xk+1=Fk,k+1X+k.(11)Xk+1式中：上标“”为先验估计值；“+”为后验估计值.第一步为根据得到的状态预测值来更新本地参考信号的码相位，j,k+1=j,k+(pj,k,k+1+tk,k+1vk+1)(aj,k+1)T.(12)pj,k,k+1jtk,k+1式中，为第颗卫星在时间的位移矢量.Zk+1X+k+1k+1第二步为根据当前历元的观测量求解状态误差的估计值，并利用该值校正式(11)所得到的状态预测值，得到第个历元导航解13：X+k+1=Xk+1+X+k+1.(13)112全 球 定 位 系 统第48卷X+k+1同时利用对式(12)中所有跟踪通道的码相位进行校正：j,k+1=j,k+1+u+k+1(a+j,k+1)T+ctb,k+1.(14)u+k+1X+k+1a+j,k+1X+k+1式中：为中的位置分量；为由校正更新后的得到的方向矢量.1.3 多径对矢量接收机的影响ikii假设第颗卫星信号在第个历元受到多径干扰，此时接收信号为式(1)所示.跟踪通道的信号为直达和反射路径的和，致使该通道的相关器输出发生变化，进一步导致鉴别器输出发生偏离，最终导致该通道的伪距误差观测量出现偏差.同时，由于其他通道信号只有直达卫星信号，因此相应的伪距观测量不存在偏差.当第个通道存在偏差的伪距误差观测量和其他通道正常的伪距误差观测量送入 EKF 时，会导致误差状态的估计值出现偏差：X+k=X+k+Xb.(15)Xbi式中，偏差来自第个跟踪通道伪距误差观测量.此时利用式(13)进行导航解更新会导致导航解出现偏差：X+k=X+k+Xb.(16)X+k利用有偏差的导航解对各跟踪通道的码相位进行校正，会导致所有跟踪通道的码相位受到存在偏差的导航解影响而出现跟踪偏差.k+1i假设导航解中速度估计不受影响，那么两个历元间接收机状态增量不存在偏差，则第个历元预测的码相位仍然存在偏差.当利用这些码相位生成本地参考信号时，则所有跟踪通道的本地参考信号都会偏离该颗卫星直达信号的码相位.对于第颗卫星信号来说，接收到的信号仍然为直达信号和反射信号的合成信号，此时对应的相关函数如图 4(a)所示，而对于其他卫星信号来说，接收信号中仅包含直达信号，因此相关函数仍然为正常的三角形，但是本地参考信号的码相位偏差会导致相关函数的中心位置发生偏移，如图 4(b)所示.受干扰信号合成信号反射信号其他信号码延时/码片码延时/码片ACFACF110011MPij(a)存在多径干扰的跟踪通道(b)不存在多径干扰的跟踪通道图 4 多径干扰下矢量跟踪中不同跟踪通道信号相关函数2矢量接收机多径检测方法本节首先给出了常用的 SQM 指标，分析这些指标在标量接收机与矢量接收机中检测多径干扰的有效性，在此基础上结合矢量接收机跟踪误差在各通道之间的传递特性，提出了基于多通道 SQM 指标联合的矢量接收机多径检测方法.2.1 常用 SQM 指标dIL1IL2d1=1.8d2=0.2表 1 给出了常用的 SQM 指标.其中 Delta 指标是利用早码和晚码之间的差异来检测相关函数的不对称性.在没有干扰的情况下，指标值通常为零.然而，当存在多径时，相关函数曲线会发生畸变，峰值两侧不再对称，如图 1 所示.此时，Delta 指标将偏离零值.Ratio 指标与 Delta 指标类似，它是通过对一对早-晚相关器的两个值求和来检测多径干扰出现时相关函数的失真，主要用于检测相关峰出现异常尖锐或平坦的情况.在本文中，Ratio 和 Delta 指标中早晚相关器的间隔设置为 0.1 码片.DoubleDelta 指标是基于由即时相关器归一化的两对早-晚相关器之间的差来定义的，被广泛应用于 GNSS 失真检测6,14.一般，两对早-晚相关器的间隔分别为 0.2 和 1 码片6.ELP指标利用早-晚相关器输出之间的相位差来检测干扰，包含同相分量和正交分量，已被证实可用于检测多径干扰.多径存在时相关峰两侧会出现失真，通常延迟的反射信号会使相关曲线的右侧更加明显，因此 Slope 指标一般由晚码相关器组成15.Slope 指标中的和对应的相关器间隔分别设置为和码片，无干扰时 Slope 指标的均值为1.表 1 常用的 SQM 指标接收机类型指标名称指标计算表达式均值标量跟踪Delta3IEILIP0Ratio3IE+ILIP1.9DoubleDelta6(IE1IL1)(IE2IL2)IP0ELP8arctan(QEIE)arctan(QLIL)0Slope152(IL1IL2)IP(d1d2)1矢量跟踪文献16IPILIE+IL|IE|IL|IPIEIE+IL|IE|Tu+NummJ()(t)Tu+NummJ()(t)TlMMP.(22)PfaMclearPdMMP式中：虚警概率为不存在干扰时超过给定门限的联合指标数与样本总量之间的比率；检测概率为多径干扰存在时超过给定门限的联合指标数与样本总量之间的比率.通过预设虚警率，根据式(21)可以计算对应的检测门限，然后由式(22)得到多径干扰下指标的检测率.然而，多通道联合的方式却不能识别出存在多径干扰的跟踪通道，因此，本文将进一步利用 DoubleDelta 指标来判断干扰所在的跟踪通道，这是基于上2.1 节所得到的该指标只有在存在多径干扰的跟踪通道时会出现异常.图 8 给出了本文所提多径检测方法的流程图，具体步骤可以描述如下：步骤 1：根据各跟踪通道的相关器输出，分别计算出其对应的 SQM 指标值(这里为 Ratio 或 Delta).mJ()步骤 2：根据步骤 1 所得到的各跟踪通道的SQM 指标值，利用式(17)计算对应的联合检验统计量.mJ()步骤 3：判断是否超过检测门限，如果超过门限，则判定有信号受到多径干扰，进入下一步识别受干扰信号；否则，判定信号间不存在多径.第3期贾琼琼，等：基于多通道 SQM 指标联合的矢量接收机多径干扰检测方法115mDDmDD步骤 4：计算各跟踪通道的 DoubleDelta 指标，当超过检测门限，则表明该跟踪通道信号受到多径干扰；否则不存在多径干扰，从而识别出多径干扰所在的跟踪通道.相关器输出SQM指标计算联合指标否否是是mJ()=m,()|N=1判断mJ()是否超过门限计算各个跟踪通道Double Delta指标mDD第个通道信号的mDD是否超过门限信号存在多径多径不存在图 8 本文方法流程图3实验分析在本节中，我们通过分析单个指标和多通道联合指标来评估多径干扰检测性能.针对 GPSL1 波段信号进行数据仿真，将所提的检测方法分别应用在矢量跟踪上，在多径干扰的环境下进行检测实验和性能评估.仿真实验的相关参数设置如表 2 所示，其中 2 号卫星信号受到多径干扰，其他信号无干扰，干扰时间为 1020s.表 2 仿真实验参数设置相关参数取值接收机中频/MHz1.405接收机采样率/MHz5.714数据长度/s30多径干扰时间/s1020受干扰卫星号2首先对表 1 中总结的常用 SQM 指标进行多径干扰下的性能分析.为了得到准确的检测概率，采用恒虚警率方法进行处理.设定虚警概率为 10%，通过式(21)计算多径检测的阈值，检测概率是通过计算每 500ms 超过门限的样本数除以该时间间隔内的样本总数来获得的，根据式(22)检测概率的计算方法可以得到所求指标在整个数据时间内的检测率，如图 9 所示.可以看出在没有多径干扰的前 10s，所有指标的检测率都在 10%左右，符合虚警率为 10%的设定；在 1020s 多径干扰期间，所有指标的检测率都出现了变化，大多数指标的检测率在 50%上下波动，而 ELP 的检测率在 40%以下低于其他指标.图 10显示了图 9 中所有指标的接受者操作特征曲线(ROC)，ROC 曲线上的点越靠近左上角，检测的准确性就越高.检测窗口为 10s，即多径干扰存在的 1020s.可以看到，随着虚警率的提高所有指标的检测率都发生相应的变化，Ratio、Delta、DoubleDelta 和Slope 这 4 个 SQM 指标的检测性能比较接近.01020时间/s检测率/%30100806040200 ELPRatio 文献16Delta SlopeDouble Delta图 9 SQM 指标检测率随时间变化趋势116全 球 定 位 系 统第48卷020406080100虚警率/%检测率/%100806040200 ELPRatio 文献16Delta SlopeDouble Delta图 10 SQM 指标检测率随虚警率变化趋势为了进一步分析多径干扰下不同 SQM 指标的检测性能，分别研究了不同指标检测率随码延时和多径相对幅度变化趋势，虚警率同样设置为 10%.图 11给出了多径相对幅度为 0.3 时，指标检测率随码延时变化结果.当码延时增大时，只有 Slope 指标的检测率出现提高；并且当码延时大于 0.5 码片时，Slope 指标的检测率明显高于其他指标.图 12 给出了多径码延时为 0.3 码片时，指标检测率随多径相对幅度的变化结果.可以看到，随着相对幅度的增大，所有指标的检测率都出现不同程度的提高，DoubleDelta 指标检测率提升最为明显，当相对幅度为 0.8 时，DoubleDelta 指标的检测率能达到 98%.检测率/%10080604020000.20.40.60.81.0多径码延时/码片ELPRatio文献16DeltaSlopeDouble Delta图 11 SQM 指标检测率随多径码延时变化趋势检测率/%10080604020000.20.40.60.8多径相对幅度 ELPRatio文献16Delta SlopeDouble Delta图 12 SQM 指标检测率随多径相对幅度变化趋势接下来分析联合指标的多径干扰检测性能.对Ratio 和 Delta 指标进行多通道联合，并与 DoubleDelta、Slope 指标进行对比分析.图 13 给出了多通道联合指标检测率随时间变化趋势.可以看到，联合后的 Ratio 和 Delta 指标与没有经过联合的 Ratio 和Delta 指标相比，检测率都有不同程度的提高.其中联合 Ratio 指标检测率提升得更加明显，检测率可达到76%，且高于其他指标.进一步分析多通道联合指标和其他指标检测率随虚警率变化的结果，如图 14 所示.显然，联合 Ratio 指标的检测性能优于其他指标，而联合 Delta 指标与 Delta 指标相比，检测性能有一定的提升.检测率/%100806040200010时间/s2030Ratio联合RatioDelta联合DeltaSlopeDouble Delta图 13 多通道联合指标检测率随时间变化趋势检测率/%100806040200020406080100虚警率/%Ratio联合RatioDelta联合DeltaSlopeDouble Delta图 14 多通道联合指标检测率随虚警率变化趋势同样地，在多径干扰存在的时间内，分析多径码延时和相对幅度对多通道联合指标检测性能的影响.图 15 为多通道联合指标检测率随多径码延时变化的结果.可以看出：联合 Ratio 指标在码延时为 0.30.6码片时，检测率高于其他指标；当码延时为 0.4 码片时，联合 Ratio 指标检测率可达到 76%.图 16 显示了多通道联合指标检测率随多径相对幅度变化的结果，当多径相对幅度小于 0.5 时，联合 Ratio 指标的检测性能优于其他指标.第3期贾琼琼，等：基于多通道 SQM 指标联合的矢量接收机多径干扰检测方法117检测率/%10080604020000.20.40.60.81.0多径码延时/码片Ratio联合RatioDelta联合DeltaSlopeDouble Delta图 15 多通道联合指标检测率随码延时变化趋势00.20.40.60.8多径相对幅度检测率/%100806040200Ratio联合RatioDelta联合DeltaSlopeDouble Delta图 16 多通道联合指标检测率随多径相对幅度变化趋势多通道联合指标检测到多径干扰后，通过二次检测识别出受干扰的信号，采用 2.2 节结合 DoubleDelta 指标的方法进行干扰识别.根据图 8 对联合Ratio 指标进行二次检测，将二次检测后检测率随虚警率变化的结果与另外四个指标进行比较，如图 17所示.可以看出，经过二次检测后的联合 Ratio 检测性能优于文献 16 和 ELP 指标.020406080100虚警率/%检测率/%100806040200ELP联合Ratio文献16二次检测SlopeDouble Delta图 17 二次检测后多通道联合指标检测率随虚警率变化趋势4结束语本文研究了矢量接收机中的多径干扰检测问题，结合矢量接收机存在误差会在跟踪通道间传递的特性，提出了基于多通道 SQM 指标联合的矢量接收机多径检测方法.通过仿真实验验证了较之单一的Ratio 和 Delta 指标，本文所提出的联合 Ratio 和联合Delta 指标多径干扰检测性能均有不同程度的提升，其中联合 Ratio 指标在多径相对幅度小于 0.5 且码延时为 0.30.6 码片时，检测性能优于其他指标.参考文献谢钢.GPS原理与接收机设计M.北京:电子工业出版社,2009.1PIRSIAVASH A,BROUMANDAN A,LACHAPELLE G.CharacterizationofsignalqualitymonitoringtechniquesformultipathdetectioninGNSSapplicationsJ.Sensors,2017,17(7):1579.DOI:10.3390/s170715792IRSIGLER M.Multipath propagation,mitigation andmonitoring in the light of Galileo and the modernizedGPSD.Munich:BundeswehrUniversity,2008.3FANTINOM,MOLINOA,MULASSANOP,etal.Signalquality monitoring:correlation mask based on ratio testmetricsformultipathdetectionC/ProceedingsofIGNSSSymposium2009.4FRANCO-PATINODM,SECO-GRANADOSG,DOVISF.SignalqualitychecksformultipathdetectioninGNSSC/International Conference on Localization and GNSS(ICL-GNSS),2013.5PIRSIAVASH A,BROUMANDAN A,LACHAPELLE G.Performance evaluation of signal quality monitoringtechniques for GNSS multipath detection and mitigationC/International Technical Symposium on Navigation andTiming(ITSNT),2017.6MUBARAKO,DEMPSTERAG.PerformancecomparisonofELPandDELPformultipathdetectionC/Proceedingsofthe 22nd International Technical Meeting of the SatelliteDivisionofTheInstituteofNavigation(IONGNSS2009),2009:2276-2283.7MUBARAKOM,DEMPSTERAG.Analysisofearlylatephase in single-and dual-frequency GPS receivers formultipathdetectionJ.GPSsolutions,2010,14(1):381-388.DOI:10.1007/s10291-010-0162-z8AMANI E,DE BOER J R,VIGNEAU W,et al.GPSmultipathinducederrorsforthevectortrackingloop:insightintomultipathdetectionC/EuropeanNavigationConference(ENC2016).2016.9AMANIE,DJOUANIK,DEBOERJR,etal.Correlator-basedmultipathdetectiontechniqueforaglobalpositioning10118全 球 定 位 系 统第48卷system/GNSS receiverJ.IET radar,sonar and navigation,2018,12(7):783-793.DOI:10.1049/iet-rsn.2017.0379LUOY,WANGYQ,WUSL,etal.MultipatheffectsonvectortrackingalgorithmforGNSSsignalJ.Sciencechinainformation sciences,2014,10(57):1-13.DOI:10.1007/s11432-014-5153-111张欣然,李洪,杨春,等.欺骗干扰对GNSS矢量跟踪环路的影响J.清华大学学报(自然科学版),2022,62(1):163-171.12ZHAO S H,LU M Q,FENG Z M.Implementation andperformanceassessmentofavectortrackingmethodbasedona software GPS receiverJ.Journal of navigation,2011,64(S1):S151-S161.DOI:10.1017/S037346331100044013PIRSIAVASH A,BROUMANDAN A,LACHAPELLE G.How effective are signal quality monitoring techniques forGNSSmultipathdetectionJ/OL.(2018-01-01)2023-01-20.InsideGNSS,2018,13(1):40-47.14KHAN A M,IQBAL N,KHAN A A,et al.Detection ofintermediate spoofing attack on global GNSS receiverthroughslopebasedmetricsJ.Journalofnavigation,2020,73(5):1052-1068.DOI:10.1017/S037346332000016815ZHANG X R,LI H,YANG C,et al.Signal qualitymonitoringspoofingidentificationmethodsforGNSSvector16trackingstructureC/Proceedingsofthe2021InternationalTechnicalMeetingofTheInstituteofNavigation,2021.SUN C,CHEONG J W,DEMPSTER A G,et al.GNSSspoofing detection by means of signal quality monitoring(SQM)metriccombinationsJ.IEEEaccess,2018,99(6):1-1.DOI:10.1109/ACCESS.2018.287594817WIPERMP,GIRONFJ,PEWSEYA.Objectivebayesianinference for the half-normal and half-t distributionsJ.Communications in statistics theory and methods,2008,37(20):3165-3185.DOI:10.1080/0361092080210518418MANFREDINIEG,DOVISF,MOTELLAB.ValidationofasignalqualitymonitoringtechniqueoverasetofspoofedscenariosC/20147thESAWorkshoponSatelliteNavigationTechnologiesandEuropeanWorkshoponGNSSSignalsandSignalProcessing(NAVITEC),2014.19作者简介贾琼琼(1986)，女，硕士，副教授，研究方向为信号处理在卫星导航中的应用研究.朱传胜(1997)，男，硕士，研究方向为卫星导航与多径干扰监测.Multipath interference detection method for vector receivers based onjoint multi-channel SQM metricsJIAQiongqiong，ZHUChuansheng(Tianjin Key Laboratory for Advanced Signal Processing,Civil Aviation University of China,Tianjin 300300,China)Abstract：Vectorreceiversprovidebetterperformanceincomplexenvironmentsthantraditionalscalarreceivers by sharing information between tracking channels to assist each other.However,multipathinterferencecanseverelylimittheperformanceofvectorreceiversincomplexenvironmentssuchasurbancanyons,andtheeffectsofmultipathcanbeeliminatedifitisdetectedanddealtwithinatimelymanner.Thispaperresearchesmultipathdetectionbasedonsignalqualitymonitoring(SQM)fromtheperspectiveofthemechanismofmultipatheffectsonvectorreceivers.Firstly,thispaperproposesamultipathdetectionmethodforvectorreceiverswithjointmulti-channelSQMmetrics,basedontheanalysisofexistingSQMmetricsforscalarreceiversignals,combinedwiththecharacteristicsofvectorreceiverchannelcoupling.Finally,thispaperverifisthroughsimulationexperimentsthattheproposedmethodcaneffectivelydetectmultipathwithamplituderatiolessthan0.5andtherelativecodedelayintherangeof0.30.6chips.Keywords：GlobalNavigationSatelliteSystem(GNSS);multipathinterference;vectortracking;signalqualitymonitoring(SQM);multi-channelcombineddetection第3期贾琼琼，等：基于多通道 SQM 指标联合的矢量接收机多径干扰检测方法119