基于WRELAX的小运算量GNSS多径干扰抑制算法

贾琼琼 吴仁彪 王文益 卢 丹 王 璐

(中国民航大学天津市智能信号与图像处理重点实验室，天津，300300)

引 言

在全球导航卫星系统逐步发展的过程中，对其性能的需求在不断提升。然而，由于卫星导航信号的脆弱性，导致其易受到干扰影响。其中多径干扰是指来自卫星的信号,被接收机天线周围物体反射后，与直达卫星信号一起进入天线，从而对接收机形成干扰。多径干扰导致的直接后果是产生跟踪和测量误差,降低定位精确度。相关研究指出，多径干扰导致的伪距误差可以达到米级甚至百米级，从而严重危及到系统的可靠性和定位精度。因此,多径干扰抑制技术一直是卫星导航接收机设计领域的研究热点[1-5]。

窄带相关技术通过减小早晚相关器间隔来降低多径对码跟踪环路的影响，是目前较为常用的一种多径干扰抑制技术。Michael S B对窄带相关器模型进行了分析[6]，Cannon M E等分析了窄带相关器技术的性能[7]，Novatel公司的Dierendonck等对其理论和性能进行了较为详细的阐述分析[8]。理论上，相关器间隔越小，对多径干扰环境下的误差改善越好。然而，这是基于接收机信号通道为无限带宽，而这个假设条件在实际中无法得到满足。当相关器间隔小于通道等效带宽的倒数时，码相位跟踪误差将趋于一个常数，因此不能通过无限降低相关器间隔的方法来减小码相位跟踪误差。

多径估计延迟锁定环技术(Multipath estimating delay lock loop, MEDLL) 是NovAtel公司1992年提出的多相关器接收机技术[9]。MEDLL的核心思想是其接收机中的多个相关器获得相关函数的多个采样值，进而根据最大似然准则进行迭代优化，这使得较之于窄带相关技术，MEDLL的复杂度更高。本团队提出的基于WRELAX算法是通过最小化非线性最小二乘(None-linear least square, NLS)代价函数的估计方法，它逐个估计出直达和多径信号的时延信息，然而该方法同样存在计算量大的问题[10]。为了进一步减小算法运算量，本文在WRELAX算法的基础上，结合卫星导航信号相关函数的特点，对接收数据与本地参考信号的互相关函数、本地参考信号的自相关函数进行加窗截取，取出相关函数中信息量较大的部分，即相关函数主峰值附近的范围，从而大大减小了WRELAX参数估计中迭代处理的数据长度，也就大大减少了算法的运算量，便于工程化。

1 数据模型

以GPS卫星信号为例，其结论可直接推广到其他GNSS系统。卫星发射信号可表示为

x(t)=d(t)c(t)cos(ωct)

(1)

式中：d(t)为导航电文，c(t)为C/A码。考虑接收机可收到P条反射的多径干扰，此时接收到的信号模型可以表示为

(2)

(3)

由于卫星距地面较远，可认为多径与直达信号具有相同的多普勒频率ωd0[10]，进一步化简式(3)得

(4)

(5)

因此，本文处理的目的是从接收到的混合数据中得到直达卫星信号的准确参数估计。

2 本文算法

基于非线性最小二乘准则的WRELAX算法[11-14]能够将多维非线性寻优搜索问题转化为一系列一维搜索，从而大大减少了算法的运算量。与现有的时延估计算法相比，WRELAX更为系统和有效，且对信号波形没有限制。下面介绍基于WRELAX的GNSS多径干扰抑制算法[10]。若对载波解调后的式(5)进一步作希尔伯特变换，得到对应的复数据可表示为

(6)

进一步对式(6)数据以Ts为间隔进行采样，可得

(7)

(8)

通过准确估计上述模型中的时延信息来抑制多径干扰。

为了得到直达信号和多径信号的初始码相位，引入相关函数

r(nTs)=corr(c(nTs),y(nTs))

(9)

式中corr(·)代表相关运算。因为信号和噪声之间是不相关的，所以式(9)可以进一步表示为

(10)

其中rc(nTs)为参考信号c(nTs)的自相关函数，w(nTs)=corr(c(nTs),e(nTs))为噪声与参考信号互相关的结果。

对于GPS卫星信号来说，C/A码的波形如图1所示，在整个C/A码的持续时间内，信号幅值不断变化，也就是C/A码的信息等概率分布在整个持续时间内。C/A码的相关函数如图2所示，可以发现其相关函数的信息量集中在±1个码片Tc范围内。因此，考虑对接收信号与本地参考信号的相关函数以及参考相关函数进行加窗截取处理，即只截取相关函数信息量较大的部分。由于当多径干扰与直达信号延迟超过1.5个码片时将不对码跟踪环路产生影响，故在截取相关函数时截取范围取±2个码片。

图1 C/A码波形 图2 C/A码的相关函数 Fig.1 Waveform of C/A code Fig.2 Correlation function of C/A code

对式(10)加窗处理后可进一步表示为

(11)

(12)

定义NLS代价函数为

(13)

(14)

(15)

(16)

那么式(13)的代价函数可以重新表示为

(17)

(18)

将式(18)代入上面的代价函数，可以得到

(19)

(20)

且

(21)

图 3 本文算法处理流程Fig.3 Flow chart of proposed algorithm

由于GPS C/A码的总长度为1 023个码片，而对相关函数加窗后其数据长度与窗函数的长度有关，若选取以相关函数峰值为中心的±2个码片范围，那么后续参数估计过程中的数据长度将会变成原相关域WRELAX算法的4/1 023，运算量也必然减少为原来的4/1 023，因此本节介绍的加窗后处理的WRELAX算法运算量大大减小，从而为工程实现提供了有利条件。

3 仿真实验

3.1 加窗的影响分析实验

仿真产生码延迟为0的卫星信号，采样频率为5 MHz，计算该卫星信号与本地参考信号之间的互相关函数如图4所示。可以看出，经过加窗后，相关函数的长度大大减小，但是在窗函数范围内的相关函数主峰值形状保持不变，窗函数之外的相关函数其余部分会有信息损失，但是这些对于参数估计不起决定作用。

加窗后前后WRELAX算法的代价函数如图5所示，可以看出，加窗前后代价函数的主瓣相当，即加窗本身对参数估计不会产生严重影响，同上面的结论一致，虽然加窗使数据长度缩短，但是其含有信息的部分却完全保留。

图4 加窗对相关函数的影响 图5 加窗前后的代价函数Fig.4 Impact of window on correlation function Fig.5 Impact of window on cost function

3.2 各种算法性能对比实验

用卫星导航软件接收机处理仿真生成卫星信号和一路多径信号。设置直达信号与多径干扰的幅值比(Direct to multipath ratio, DMR)分别取0.1,0.3,0.5,0.7，相对延迟τ2-τ1从0～1.5个码片，无噪声、无限带宽条件下，分别使用传统相关器、窄相关器、MEDLL和本文所提两种基于信号分离的时延估计算法得到的多径误差包络如图6所示。

图6 不同DMR下各方法的误差包络曲线Fig.6 Comparison of code delay estimation error envelope in different DMR

从图6的实验结果可以看出，相关器间隔d=1时传统相关器存在较大的时延估计误差，随着相关器间隔的缩小，时延估计误差也逐渐减小。MEDLL技术、相关模型WRELAX及加窗后处理的算法作为高分辨率的时延估计方法能提供比窄带相关器更好的性能。为了进一步观察这3种基于时延估计方法的性能，对上述实验结果进行放大得到图7。从图7中可以看出，当DMR=0.1～0.3时，这些方法的性能相当，也就是说，多径干扰比较弱时(实际情况中经过反射后的多径干扰都会产生严重的衰减)，加窗后处理的方法在保证能以较低运算量实现的前提下也能保证其性能不受较大损失；当多径干扰逐渐增强，如DMR=0.5时，加窗后处理的方法性能会逐渐劣于相关域WRELAX算法和MEDLL，但是其性能仍优于窄带相关技术。加窗后处理的方法由于没有考虑主峰值之外的部分，因此参数估计的迭代过程会导致这部分的信息有所损失，最终导致其估计性能较之原始的WRELAX算法有所下降。

图7 不同DMR下各方法的误差包络曲线(细节放大)Fig.7 Comparison of code delay estimation error envelope in different DMR (Detail amplification)

本文所提算法是从减少运算量的角度，通过缩减参数搜索范围对文献[15]所提WRELAX算法的改进，文献[15]给出了噪声情况下WRELAX的性能分析，证明了相对于窄相关和MEDLL等常用的多径抑制技术，WRELAX算法具有良好的性能。通过理论分析和仿真实验证明了所提降低运算量的算法与原来的WRELAX算法性能相当，由于新算法只是巧妙地利用卫星信号自身良好的自相关特性来降低参数估计的搜索区间，因此不会影响算法的抗噪声性能。

4 结束语

本文根据GNSS信号相关函数的特点，提出了一种加窗处理的相关域WRELAX多径干扰抑制算法。该算法通过对接收数据与本地参考信号的相关函数以及本地参考信号的自相关函数加窗处理，使得参与参数估计的相关函数长度大大减小，而由于相关函数的信息量主要集中在相关峰值附近的范围内，因此加窗后相关函数的信息量不会有明显损失。最后，通过仿真实验进一步验证了本文算法能够在以较低运算量实现的前提下，同时保证性能不会受到太大损失，尤其是多径干扰较弱时。

参考文献:

[1] Kalyanaraman S K, Braasch M S, Kelly J M. Code tracking architecture influence on GPS carrier multipath[J]. IEEE Transactions on Aerospace and Electronic Systems, 2006,42(2):548-561.

[2] Kos T, Markezic I, Pokrajcic J. Effects of multipath reception on GPS positioning performance[C]// ELMAR, 2010 Proceedings. [S.l.]: IEEE, 2010:399-402.

[3] 贺成艳,卢晓春.用于GNSS基带信号分析的相干累加算法[J]. 数据采集与处理,2009，23(S1):124-127 .

He Chengyan , Lu Xiaochun. New coherent accumulation algorithm for further property analysis of GNSS baseband signals[J]. Journal of Data Acquisition and Processing, 2009，23(S1):124-127.

[4] 付龙飞,张水莲,李世凯. 一种新的码片内多径时延估计方法[J]. 数据采集与处理,2012,26(2):189-195.

Fu Longfei, Zhang Shuilian, Li Shikai. New subchip multipath time-delay estimation method[J]. Journal of Data Acquisition and Processing, 2012,26(2):189-195.

[5] Alfred R L. GPS landing system reference antenna[J]. IEEE Antennas and Propagation Society, 2010,52(1):104-113.

[6] Michael S B. GPS multipath model validation[C]//IEEE Position Location and Navigation Symposium Conference Proceedings. Atlanta: GA, 1996:672-678.

[7] Cannon M E，Lachapelle G，Qiu W，et al. Performance analysis of a narrow correlator spacing receiver for precise static GPS positioning[C]// IEEE Position Location and Navigation Symposium. Las Vegas: IEEE, 1994:2994-2999.

[8] Van Dierendonck A J，Fenton P，Ford T. Theory and performance of narrow correlator spacing in a GPS receiver[J]. Journal of The Institute of Navigation, 1992,39(3):265-283.

[9] Townsend B R, Fenton P C, Dierendonck K J V, et al. Performance evaluation of the multipath estimating delay lock loop[J]. Navigation, 1995, 42(3):502-514.

[10] 李杰，吴仁彪，卢丹，等. 基于信号分离估计理论的GPS多径抑制算法 [J].信号处理，2011,27(12):1884-1888.

Li Jie, Wu Renbiao, Lu Dan, et al. GPS multipath mitigation algorithm based on signal separation estimation theory[J]. Chinese Signal Processing，2011,27(12):1884-1888.

[11] Li Jian, Wu Renbiao. An efficient algorithm for time delay estimation[J]. IEEE Transactions on Signal Processing, 1998,46(8):2231-2235.

[12] Li Jian, Zheng Dunmin, Stoica P. Angle and waveform estimation via relax[J]. IEEE Transactions on Aerospace and Electronic Systems, 1997,33(3):1077-1087.

[13] Wu Renbiao, Li Jian, Liu Zhengshe. Super resolution time delay estimation via MODE-WRELAX [J]. IEEE Transactions on Aerospace and Electronic Systems, 1999,35(1):294-307.

[14] Wu Renbiao, Li Jian. Time-delay estimation via optimizing highly oscillatory cost functions[J]. IEEE Journal of Oceanic Engineering, 1998,23(3):235-244.

[15] Jia Qiongqiong, Wu Renbiao, Wang Wenyi, et al. Multipath interference mitigation in GNSS via WRELAX[J]. GPS Solutions, 2017, 21(2):487-498.