邵 晨，曾慶化,2,3，邱文旗，許 睿

(1. 南京航空航天大學導航研究中心，南京 211106；2.先進飛行器導航、控制與健康管理工業(yè)和信息化部重點實驗室，南京 211106；3.衛(wèi)星通信與導航江蘇高校協(xié)同創(chuàng)新中心, 南京 211106)

0 引言

全球?qū)Ш叫l(wèi)星系統(tǒng)(Global Navigation Satellite System，GNSS)多徑效應是指GNSS信號由于傳輸過程中存在遮擋，使接收信號中混有多徑信號，造成信號幅值、載波相位和偽碼延遲變化，引入定位解算誤差，從而影響接收機定位精度的現(xiàn)象。根據(jù)偽隨機噪聲碼(Pseudo Random Noise code，PRN)的自相關(guān)特性，當多徑信號相對于直達信號的延遲在1個碼片之內(nèi)時，接收機接收的合成信號和本地產(chǎn)生信號之間的相關(guān)函數(shù)發(fā)生畸變，產(chǎn)生解算誤差，這種由多徑效應引起的誤差稱為多徑誤差。

隨著GNSS現(xiàn)代化進程的推進，用戶對于定位精度的要求也越來越高。影響GNSS定位精度的誤差有很多,部分如鐘差、對流層誤差、電離層誤差等可以通過改善模型、差分技術(shù)等方法進行抑制;而多徑誤差由于沒有空間相關(guān)性，很難通過上述方法削弱，在誤差源中的影響越來越大。因此，研究多徑抑制技術(shù)對提高衛(wèi)星導航定位精度具有重要意義。

由圖1可知，多徑誤差產(chǎn)生的環(huán)境分為視距(Line-of-Sight，LOS)環(huán)境和非視距(Non-Line-of-Sight，NLOS)環(huán)境，視距環(huán)境下接收信號中含有直射信號，而非視距環(huán)境下的接收信號中沒有直射信號。由于GNSS多徑抑制基帶處理算法僅適用于視距環(huán)境下多徑誤差的抑制，本文對非視距環(huán)境下的多徑誤差抑制技術(shù)不再詳述。

(a) 非視距環(huán)境

(b) 視距環(huán)境圖1 多徑誤差生成過程Fig.1 Generation of multipath errors

視距環(huán)境下,多徑誤差抑制技術(shù)分為天線增強技術(shù)、基帶處理算法和后處理技術(shù)三類，如圖2所示。

圖2 多徑抑制技術(shù)Fig.2 Multipath mitigation technique

天線增強技術(shù)通過改進天線，使方向圖主瓣對準直射信號方向，減少反射多徑信號的接收，該方法主要在硬件層面通過天線選址、特殊天線、安裝扼流圈等方式，抑制來自地面或低仰角的多徑信號，但難以消除來自空中的多徑信號。

基帶處理算法主要通過分析信號結(jié)構(gòu)，改進跟蹤環(huán),減小或消除多徑信號對于自相關(guān)函數(shù)的影響，從而減小多徑誤差。

后處理技術(shù)通過對接收機測量數(shù)據(jù)進行后期處理以減小多徑誤差。如利用高斯牛頓方法跟蹤多徑增益及延遲參數(shù)，利用后期建筑信息協(xié)助多徑判斷進行選星，該類方法需要其他輔助信息和大量高精度測量數(shù)據(jù)，多徑誤差抑制效果好但實時性較差。

上述多徑抑制技術(shù)中，基帶處理算法受信號體制影響較大，衛(wèi)星播發(fā)信號的不斷改善和發(fā)展給多徑抑制基帶處理算法的研究帶來了更大的機遇和挑戰(zhàn)。本文從信號體制角度出發(fā)，對GNSS多徑抑制基帶處理算法的研究現(xiàn)狀，典型算法的基本原理、優(yōu)缺點和發(fā)展趨勢做了較為詳細的論述，其中部分典型算法如表1所示。

由表1可知，多徑抑制基帶算法目前的相關(guān)研究分為兩類：第一類是對傳統(tǒng)多徑抑制算法進行研究，提高傳統(tǒng)算法的跟蹤精度和抗干擾性能，典型算法有窄相關(guān)技術(shù)、ELS、CCRW和MEDLL等；第二類是對新體制信號如北斗B1C、B2a信號等進行研究，充分利用多峰、數(shù)據(jù)導頻分離、主瓣帶寬大等信號結(jié)構(gòu)特點來提高環(huán)路跟蹤性能，如基于去模糊的多徑抑制算法，能有效抑制或消除BOC信號的多峰特性從而避免誤鎖，同時獲得更好的抗多徑性能，典型算法有BPSK-like、Bump-Jump、ASPeCT、PCF等，下文對此進行詳細討論。

表1 GNSS多徑抑制基帶處理典型算法

1 傳統(tǒng)多徑抑制算法

GNSS多徑信號對接收機信號處理的影響主要體現(xiàn)在兩方面：

一方面,合成信號本身引入了載波和偽距測量值誤差，直接影響跟蹤環(huán)路對偽距等重要導航信息的解算。

另一方面,碼跟蹤環(huán)通過接收信號與本地參考信號的相關(guān)函數(shù)進行碼相位識別，如圖3所示，若接收信號中含有多徑信號，會使得本地信號和接收信號的相關(guān)函數(shù)發(fā)生畸變，鑒相器在判斷峰值位置時出現(xiàn)偏差，從而導致跟蹤誤差。

圖3 合成信號和直達信號歸一化相關(guān)函數(shù)Fig.3 Normalized correlation function of synthetic signal and direct signal

多徑信號也會對載波跟蹤產(chǎn)生影響，但由于多徑信號對碼環(huán)的影響遠大于對載波環(huán)的影響，所以多徑抑制算法研究多在碼環(huán)。傳統(tǒng)多徑抑制算法從原理上主要分為參量式和非參量式，對于中長時延多徑的抑制效果較好。

1.1 非參量式

非參量算法通過改進接收機的環(huán)路結(jié)構(gòu)和鑒相器函數(shù)以降低對多徑信號的敏感程度。1992年，文獻[6]提出了窄相關(guān)技術(shù)，通過壓縮碼片間隔提高抗多徑性能。1994年，文獻[7]提出了ELS(Early Late Slope)技術(shù)，通過對相關(guān)函數(shù)兩側(cè)斜率進行估計,判斷相關(guān)函數(shù)的畸變程度，跟蹤精度比窄相關(guān)技術(shù)更高。1996年，文獻[8]提出了Strobe相關(guān)器技術(shù)，文獻[10]提出了脈沖間隙相關(guān)器(Pulse Aperture Correlator，PAC)技術(shù)。2002年，文獻[11]提出了碼相關(guān)參考波形(Code Correlation Reference Wavefo-rms，CCRW)技術(shù)，通過改變本地碼波形獲得理想的鑒相函數(shù)以抑制多徑。2006年，文獻[12]在CCRW的基礎(chǔ)上提出了一種參考波形的設(shè)計方法，能夠獲得期望的相關(guān)函數(shù)。2009年，文獻[13]提出了高分辨率相關(guān)器(High Resolution Correlator，HRC)。2012年，文獻[14]提出了一種基于小波變換的多徑抑制方法，利用小波變換判斷偽碼相關(guān)函數(shù)的奇異點，從而實現(xiàn)多徑抑制。2013年，文獻[15]在CCRW技術(shù)的基礎(chǔ)上進一步研究了參考波形和閘寬對于多徑抑制的影響。2018年，文獻[16]提出了一種新的參考波形，對于高階BOC信號具有良好的多徑抑制效果。2020年，文獻[17]提出了一種可滿足多目標約束條件的參考波形設(shè)計方法。其中，Strobe相關(guān)器、PAC和HRC統(tǒng)稱為Double-Delta技術(shù)，相關(guān)器間距相同時，該類技術(shù)的跟蹤精度較窄相關(guān)技術(shù)更高。CCRW技術(shù)通過設(shè)計不同的參考波形可以達到本地相關(guān)器線性組合的效果，相較于Double-Delta具備更為優(yōu)異的抗多徑性能。非參量式算法的跟蹤環(huán)路結(jié)構(gòu)簡單，跟蹤速度較快，但多徑性能受以下限制：1)受鑒相器性能和相關(guān)器間距的寬度限制等影響，不能完全消除多徑影響；2)受前端帶寬影響，前端帶寬越窄時，相關(guān)峰越平，鑒相函數(shù)多徑分辨率越低，需要增加帶寬和信號采樣頻率，但會導致環(huán)路抗干擾和跟蹤性能下降。

1.1.1 窄相關(guān)技術(shù)

衛(wèi)星導航發(fā)展的早期，接收機碼跟蹤環(huán)的超前滯后相關(guān)器間距(Code Spacing)設(shè)置較大，一般為1個碼片，能夠減小初始捕獲時間，但是會導致碼環(huán)鑒相器函數(shù)單調(diào)范圍較寬，跟蹤精度較低，無法滿足高精度定位需求。1992年，NovAtel公司提出了窄相關(guān)技術(shù)以提高跟蹤精度。窄相關(guān)的碼跟蹤鑒相器是超前碼與滯后碼相減得到的，該技術(shù)通過減小相關(guān)器間距，提高鑒相器對峰值點的靈敏度，降低多徑信號帶來的誤差，但是不能完全消除多徑誤差，具有一定局限性。一般而言，相關(guān)器間距越小，多徑誤差包絡(luò)面積越小，抗多徑性能越好，但隨著相關(guān)器間距的減小，環(huán)路穩(wěn)定性能變差，容易造成失鎖等現(xiàn)象。

1.1.2 Double-Delta技術(shù)

Double-Delta技術(shù)通過兩對相關(guān)器構(gòu)建鑒相函數(shù)，主要包括HRC技術(shù)、Strobe相關(guān)器技術(shù)等。鑒相函數(shù)表達式為

(1)

其中，、、、為兩組超前滯后相關(guān)值，相距2，通過選取不同的值達到不同的效果。HRC技術(shù)選取的值為1，Strobe相關(guān)器則為2，兩者抗多徑性能相近。

圖4中，表示偽碼延遲，表示相關(guān)函數(shù)。由圖4可知，相關(guān)器間距越小，Double-Delta的抗多徑能力越強，相關(guān)器間距相同時，Double-Delta技術(shù)的抗多徑性能要優(yōu)于窄相關(guān)技術(shù)。

圖4 Double-Delta技術(shù)Fig.4 Double-Delta technology

1.1.3 ELS技術(shù)

1994年，文獻[7]提出了ELS技術(shù)，相較于Dou-ble-Delta技術(shù)，ELS技術(shù)同樣是使用兩組相關(guān)器，但是鑒相器設(shè)計的原理不同。ELS技術(shù)利用多徑信號導致碼相關(guān)函數(shù)在峰值兩側(cè)主瓣的坡度不同這一特性，對峰值兩側(cè)坡度進行估計，從而準確得出峰值的位置。

由圖5可知，、、、為兩組超前滯后相關(guān)值，ELS技術(shù)計算得到的多徑誤差為峰值兩側(cè)直線交點的橫坐標，、、分別為相鄰相關(guān)器之間的間隔，一般取==12，此時鑒相器表達式為

(2)

與窄相關(guān)技術(shù)類似，相關(guān)器間隔越小，抗多徑性能越好。ELS技術(shù)依賴于兩側(cè)相關(guān)函數(shù)斜率估計的準確度，抗多徑性能與窄相關(guān)技術(shù)相比有明顯的提升，但是在信噪比較低時，斜率的估計往往不夠準確，會造成較大的誤差。

圖5 ELS技術(shù)Fig.5 ELS technology

1.1.4 CCRW技術(shù)

2002年，文獻[11]提出了CCRW技術(shù)，經(jīng)過多次改進已較為成熟。該方法的原理是通過設(shè)計特殊的參考波形作為本地碼，將其與接收信號相關(guān)后得到不同的鑒相曲線，抗多徑性能與參考波形設(shè)計相關(guān)?；緟⒖疾ㄐ沃饕袃深悾旱谝活悆H在碼片變化處產(chǎn)生參考波形，第二類在每個碼片邊緣處產(chǎn)生并隨之變化，部分波形的產(chǎn)生方式如圖6所示。

圖6 CCRW參考波形Fig.6 Reference waveform of CCRW

不同參考波形可以選擇相同的鑒相函數(shù)，歸一化表達式為

(3)

其中，、分別為同相、正交分量的即時相關(guān)值；、分別為同相、正交分量的參考波相關(guān)值。當參考波形為矩形碼時，CCRW技術(shù)與窄相關(guān)技術(shù)等效，參考波形為時，CCRW技術(shù)與Double-Delta技術(shù)等效。通過選擇不同的參考波，CCRW可以達到與其他非參量式算法類似的效果，無需增加額外相關(guān)器，環(huán)路簡單。參考波的閘寬越小，多徑抑制性能越好，參考波形為和時,抗多徑性能比Double-Delta技術(shù)更優(yōu)。

1.2 參量式

參量式算法的原理是通過大量的相關(guān)器獲取多徑信號的特征參數(shù)觀測值，再通過估計算法得到多徑參數(shù)最優(yōu)估計值，最后剝離出直射信號。1995年，文獻[20]提出了多徑消除延遲鎖相環(huán)(Multipath Estimation Delay Locked Loop，MEDLL)技術(shù)，通過估計多徑參數(shù)有效抑制中長時延多徑信號，但該方法計算量大，實時性較差。1997年，文獻[21]提出了多徑消除技術(shù)(Multipath Mitigation Technology，MMT)，在MEDLL的基礎(chǔ)上通過預先假設(shè)多徑個數(shù),簡化多徑模型，降低了計算量。2009年，文獻[22]提出了Vision相關(guān)器(Vision Correlator)技術(shù)，通過精確測量碼片過渡時刻的射頻信號性能，提取清晰的碼片過渡時刻波形，再剝離多徑信號，多徑性能相較于MMT得到進一步提高和優(yōu)化。這類方法對中長時延多徑誤差的抑制效果好，跟蹤精度高，但是由于需要估計多徑信號并進行剝離，所以需要較多的相關(guān)器，跟蹤環(huán)路結(jié)構(gòu)較為復雜，計算量較大。除此之外，隨著濾波和智能算法研究的深入，一些改進算法相繼被提出。2017年，文獻[23]采用遺傳粒子濾波對多徑信號進行估計,提高了多徑估計的精度。2018年，文獻[24]提出了基于差分進化(Differential Evolution，DE)改進粒子濾波的多徑估計算法，提高了非高斯噪聲下的估計精度。2019年，文獻[25]提出了利用擴展卡爾曼濾波(Extended Kalman Filter，EKF)提高MEDLL的多徑估計精度。同年，文獻[26]利用信息論學習實現(xiàn)多徑估計，該算法無需假設(shè)多徑數(shù)目,能夠快速實現(xiàn)參數(shù)估計。其他還有利用Newton法優(yōu)化求解，加入TK算子降低計算復雜度等針對最大似然估計算法的改進。參量式算法需要估計多徑信號的特征信息，跟蹤精度相對于非參量式算法更高，但環(huán)路設(shè)計復雜，計算量較大。

接下來以MEDLL技術(shù)為例介紹參量式算法的實現(xiàn)過程，該技術(shù)是一種基于多徑模型的參數(shù)估計方法,通過最大似然估計算法估計多徑偽碼延遲、載波延遲和幅值。

圖7 MEDLL跟蹤環(huán)Fig.7 Tracking loop of MEDLL

由圖7可知，中頻信號分別與本地載波sin()、cos()相乘分為(同相)分量和(正交)分量，之后分別進入2+1個相關(guān)器，每個相關(guān)器間隔為個碼片。整個跟蹤過程中含有大量待估參數(shù)，一般使用迭代求解。迭代過程分為參數(shù)估計、相關(guān)峰恢復和檢測流程控制。其中，參數(shù)估計過程為假設(shè)輸入僅有一路信號，估計該信號幅值、載波相位和偽碼延遲；相關(guān)峰恢復過程為根據(jù)給定的幅值、載波相位和偽碼延遲，輸出估計的相關(guān)峰波形；檢測流程為根據(jù)估計波形的殘差預設(shè)多徑路數(shù)，判斷是否進行下一輪迭代。針對檢測流程繁瑣的問題，文獻[29]提出了一種改進的MEDLL技術(shù)，該算法通過設(shè)定接收信號中最大路徑數(shù)進行余量估計，能夠有效提高直射信號的估計效率，但是依舊沒有解決MEDLL需要大量相關(guān)器、計算量大的主要問題。MEDLL算法是目前最有效的多徑檢測和估計算法之一，跟蹤精度較高，實時性較差，多徑抑制性能主要受噪聲影響。

1.3 傳統(tǒng)多徑抑制算法應用挑戰(zhàn)：短時延多徑

一般而言，以上方法可以在多數(shù)場合發(fā)揮較好的多徑抑制作用，然而它們都有一個共同的缺陷，當面對短時延多徑信號(多徑延遲小于0.5碼片)時，上述方法的抗多徑性能會有一定的下降。短時延多徑信號是多徑信號中的一類特殊信號，兩者具有相似的特性，但短多徑信號的延遲更小。短多徑抑制技術(shù)的難點在于短多徑中直達信號的剝離和延遲的估計。關(guān)于短多徑抑制算法的研究如下，2009年，文獻[30]在HRC相關(guān)器中增加了修正因子以減小短多徑誤差。2010年，文獻[31]提出了一種基于斜率的多徑估計(the Slope-Based Multipath Estimation，SBME)方法，這種方法能夠很好地抑制延遲小于0.35碼片的短多徑信號。2011年，文獻[32]提出了一種新型碼環(huán)路鑒相器(Short Multipath Insensitive Code Loop Discriminator，SMICLD)，這種方法碼環(huán)路簡單，對短時延多徑有較好的抑制效果。2012年，文獻[33]在SMICLD的基礎(chǔ)上進行改進，提出了MSMICD鑒相器，該方法計算量小，處理同相多徑信號時性能優(yōu)于SMICLD。2014年，文獻[34]提出了改進窄相關(guān)的短多徑抑制方法，在傳統(tǒng)窄相關(guān)器的基礎(chǔ)上，額外增加一個相關(guān)器用于對短時延多徑誤差進行定量估計并實時修正。同年，文獻[1]提出了在Strobe相關(guān)器算法中引入修正因子，對碼環(huán)鑒相器進行修正以抑制短多徑。2019年，文獻[35]提出了一種改進的對短多徑不敏感的Strobe鑒相(Modified Short Multipath Insensitive Strobe Discriminator，MSMISD)算法，通過構(gòu)造延遲相關(guān)值和設(shè)計反相多徑補償因子抑制短多徑，相較于傳統(tǒng)Strobe相關(guān)器抗多徑性能更好。短延時多徑抑制基帶處理算法研究相對較少，大都是在傳統(tǒng)長時延多徑抑制算法的基礎(chǔ)上進行改進，精度和抗干擾性能還有待提高。

2 新體制信號多徑抑制算法

隨著衛(wèi)星導航定位系統(tǒng)的發(fā)展，現(xiàn)代化衛(wèi)星導航系統(tǒng)通過采用二進制偏移載波(Binary Offset Carrier，BOC)調(diào)制和混合二進制偏移載波(Multiplexed Binary Offset Carrier，MBOC)調(diào)制解決頻段擁擠的問題。BOC信號的數(shù)學模型表示如下

()=(-)cos(2π+)

(4)

=sign(sin(2π+))

(5)

其中，為信號幅值；為PRN碼；為偽碼延遲(單位是碼片)；為載波頻率；為載波相位；為方波子載波；為子載波頻率；決定子載波類型，取0時為sin-boc，取 12π時為cos-boc。

新體制信號從信號設(shè)計的角度減小多徑誤差，帶來更高性能的同時提升了接收機設(shè)備的設(shè)計復雜度?，F(xiàn)代化GNSS部分信號已采用BOC調(diào)制或在其基礎(chǔ)上進行變形的調(diào)制方式，如北斗的B1C信號采用QMBOC調(diào)制，Galileo系統(tǒng)中E1信號采用CBOC調(diào)制和GPS中L1C采用TMBOC調(diào)制等。

新體制信號多徑抑制算法針對于新體制信號設(shè)計，主要分為兩類：一類是基于去模糊的多徑抑制算法，另一類是針對不同信號特性設(shè)計的處理方法。

2.1 基于去模糊的多徑抑制算法

BOC、MBOC信號雖然有著良好的抗多徑性能，但是由于其自相關(guān)函數(shù)的多峰特性，容易導致跟蹤環(huán)路發(fā)生誤鎖現(xiàn)象，在使用Strobe相關(guān)器、HRC等常用多徑抑制算法時，在±0.5碼片附近存在模糊點。為解決該問題，一般采用基于去模糊的多徑抑制算法，這類方法由去模糊算法和傳統(tǒng)多徑抑制算法組成。首先通過去模糊算法解決BOC信號的模糊性問題，降低多徑敏感度，再通過傳統(tǒng)多徑抑制算法進一步提高多徑抑制性能。由于傳統(tǒng)多徑抑制算法上文已提及，接下來具體闡述去模糊算法的發(fā)展歷程。

1999年，文獻[38]提出了峰跳法(Bump-Jump)，在傳統(tǒng)BPSK跟蹤方法的基礎(chǔ)上增加了遠超前和遠滯后相關(guān)器，監(jiān)測主峰兩側(cè)相關(guān)峰值，判斷是否鎖定在主峰，該方法需要積累較長的時間來完成失鎖檢測。2000年，文獻[39]提出了非相干邊帶算法，通過對BOC信號進行單邊帶濾波消除多峰特性，但是該算法會造成信號能量的浪費，導致跟蹤精度降低。2003年，文獻[40]提出了一種基于子載波相關(guān)的去模糊算法。2004年，文獻[41]提出了一種亞載波相位消除的去模糊算法，這類方法通過消除矩形波的方式獲得單峰相關(guān)函數(shù)。2006年，文獻[42]提出了BPSK-like算法，相較于非相干邊帶法的優(yōu)點在于能夠同時利用兩個邊帶能量,提高接收機處理的衛(wèi)星信號能量。2007年，文獻[43]提出了信號相關(guān)邊峰消除(Autocorrelation Side-Peak Cancellation Technique， ASPeCT)算法。同年，文獻[44]提出了一種基于偽相關(guān)函數(shù)(Pseudo Correlation Function，PCF)的算法。2013年，文獻[48]提出了一種MBOC信號去模糊算法，設(shè)計了雙超前減滯后鑒相器去除MBOC信號的多峰性。2014年，文獻[49]提出了一種基于PCF算法的CBOC信號去模糊算法，跟蹤性能更優(yōu)。2018年，文獻[50]提出了一種基于CCRW算法的QMBOC信號去模糊算法，能有效消除多峰性，但由于閘波設(shè)計不便，具有一定局限性。

采用去模糊算法消除BOC信號的多峰性是后續(xù)使用鑒相器改進、參數(shù)估計等傳統(tǒng)多徑抑制算法的前提，同時可以通過構(gòu)造特殊相關(guān)函數(shù)進一步提高整個算法的多徑抑制性能。接下來對ASPeCT、PCF等典型去模糊算法進行介紹，這類方法通過多個互相關(guān)函數(shù)的非線性組合對原相關(guān)函數(shù)進行改進，重構(gòu)的相關(guān)函數(shù)的多峰性大大減弱，從而減小了跟蹤環(huán)路誤鎖的可能性。

2.1.1 ASPeCT算法

2007年，文獻[43]提出了ASPeCT算法，經(jīng)過多次改進已較為成熟。該算法是一種基于自相關(guān)邊峰消除技術(shù)的跟蹤算法，利用BOC信號與PRN碼的互相關(guān)函數(shù)和BOC信號的自相關(guān)函數(shù)，對碼環(huán)路相關(guān)結(jié)果進行重構(gòu)，重構(gòu)之后的相關(guān)函數(shù)為

(6)

其中，()是BOC信號的自相關(guān)函數(shù)；()是BOC信號和本地PRN序列的互相關(guān)函數(shù)；是調(diào)節(jié)因子，其取值影響ASPeCT相關(guān)函數(shù)中兩個較小邊鋒的大小，需要根據(jù)應用的不同需求進行調(diào)整。

由圖8可知，BOC(,)信號與本地偽碼序列的互相關(guān)函數(shù)曲線關(guān)于原點中心對稱，且存在兩個峰值。這兩個峰值的碼相位延遲與BOC(,)信號自相關(guān)函數(shù)的兩個邊峰所在的碼相位延遲相同，位于±0.5碼片處。由于兩側(cè)峰峰值符號相反，ASPeCT算法利用平方運算消除符號對相關(guān)峰的影響，從而削弱側(cè)峰對整個偽碼跟蹤環(huán)路的影響。ASPeCT算法構(gòu)建的相關(guān)函數(shù)沒有削弱主峰的幅值和寬度，保持了BOC(,)信號原有的特性，跟蹤性能上優(yōu)于傳統(tǒng)的BPSK跟蹤方法;但是該方法并不能完全消除自相關(guān)邊峰，在主峰的兩側(cè)仍然留有兩個較小的邊峰，初始相位較大時可能會導致跟蹤不穩(wěn)定。針對該問題，文獻[51]提出了改進ASPeCT算法相關(guān)函數(shù)，構(gòu)建一個無邊峰的相關(guān)函數(shù)，新相關(guān)函數(shù)的相關(guān)主峰相對于改進前更寬，這會導致環(huán)路的跟蹤精度降低，但是新相關(guān)函數(shù)只有一個主峰，完全消除了邊峰的影響，穩(wěn)定性更好。ASPeCT算法構(gòu)建的相關(guān)函數(shù)只能消除BOC(,)的邊峰影響，對于BOC(,)信號和MBOC信號不適用。

圖8 ASPeCT歸一化相關(guān)函數(shù)Fig.8 Normalized correlation function of ASPeCT

2.1.2 PCF算法

2007年，文獻[44]提出了PCF算法，經(jīng)過多次改進已較為成熟，是一種通過構(gòu)建偽相關(guān)函數(shù)實現(xiàn)無模糊跟蹤的算法。PCF算法的主要原理為：在中頻信號處理過程中構(gòu)建兩組本地參考信號()和()，分別與接收信號進行相關(guān)運算，相關(guān)函數(shù)為()和()，最后利用相關(guān)結(jié)果構(gòu)建無相關(guān)邊峰的偽相關(guān)函數(shù)，跟蹤性能受寬度因子影響。相關(guān)函數(shù)為

()=|()|+|()|-
|()+()|

(7)

相關(guān)函數(shù)圖像如圖9所示。

圖9 PCF歸一化相關(guān)函數(shù)Fig.9 Normalized correlation function of PCF

由圖9可知，PCF算法能夠徹底消除BOC信號的多峰性，同時利用合適的寬度因子獲得更窄的相關(guān)主峰和更好的抗多徑性能。2021年，文獻[52]提出了一種子函數(shù)組合相關(guān)的無模糊跟蹤算法，有效地適用于對BOC(,)和CBOC(6,1,1/11)信號的跟蹤，但是與PCF算法相似，需要設(shè)計特殊的參考波形，組合函數(shù)在實現(xiàn)無模糊的同時會損失一部分跟蹤性能。

2.2 其他處理方法

新體制信號多徑抑制算法除了去模糊算法與傳統(tǒng)多徑抑制算法組合的方式之外，還可以根據(jù)信號結(jié)構(gòu)特征設(shè)計特殊算法，提高信號利用功率和多徑抑制性能。該類方法已實現(xiàn)工程應用的有聯(lián)合跟蹤、TM61跟蹤算法和雙重估計技術(shù)等。

2.2.1 聯(lián)合跟蹤

部分新體制信號為數(shù)據(jù)、導頻分離的結(jié)構(gòu)，數(shù)據(jù)分量包含電文信息，用于獲取偽距和導航信息，導頻分量僅有擴頻碼，不包含電文信息，從而增加了相干積分時間，精度更高。傳統(tǒng)的多徑抑制算法僅針對單路信號進行觀測量的估計，這會造成有用功率的浪費，從而導致跟蹤精度的降低，針對該問題，可以使用聯(lián)合跟蹤算法解決。2002年，文獻[54]針對L5信號提出了鑒相器和濾波器層面的聯(lián)合跟蹤方法，并研究了該方法的適用性。2010年，文獻[55]針對E1信號提出了數(shù)據(jù)/導頻聯(lián)合捕獲跟蹤算法，提高了捕獲跟蹤的穩(wěn)定性。2011年，文獻[56]提出了不同權(quán)重的非相干和相干數(shù)據(jù)/導頻聯(lián)合方法，提高了信號的利用率。同年，文獻[57]針對L5信號提出了三種復雜度不同的聯(lián)合跟蹤方法，分別為相關(guān)器輸出聯(lián)合、鑒相器輸出聯(lián)合和濾波器輸出聯(lián)合。2013年,文獻[15]提出了聯(lián)合多徑抑制算法，將CBOC信號的導頻和數(shù)據(jù)通道利用碼參考波形技術(shù)實現(xiàn)多徑抑制，再進行非相干合并消除BOC(6,1)分量的影響，達到了該信號的多徑抑制極限。2018年，文獻[58]針對B1C信號，采用導頻輔助數(shù)據(jù)聯(lián)合跟蹤和數(shù)據(jù)/導頻幅值聯(lián)合跟蹤，提高了跟蹤的精度和穩(wěn)定性。2019年，文獻[51]分析比較了三種聯(lián)合跟蹤算法的跟蹤精度和抗多徑性能，為以后的研究提供參考。2020年，文獻[59]提出了一種B1C信號的無模糊捕獲算法，采用數(shù)據(jù)/導頻非相干組合的捕獲策略，實現(xiàn)了B1C信號的無模糊捕獲。同年，文獻[60]提出了一種基于EKF的數(shù)據(jù)/導頻聯(lián)合跟蹤方法，在聯(lián)合跟蹤模型的基礎(chǔ)上引入卡爾曼濾波器，進一步減小跟蹤誤差。

由圖10可知，聯(lián)合跟蹤方式分為相關(guān)器、鑒相器和濾波器三類。相關(guān)器聯(lián)合通過數(shù)據(jù)和導頻相關(guān)積分進行組合，需要判斷電文與二次碼符號，環(huán)路結(jié)構(gòu)最簡單，但是在低載噪比時符號判斷受影響，從而導致精度下降。鑒相器聯(lián)合通過信號發(fā)射功率比對數(shù)據(jù)和導頻鑒相器輸出進行加權(quán)組合，環(huán)路結(jié)構(gòu)相比相關(guān)器聯(lián)合更復雜，在低載噪比環(huán)境下精度更高。濾波器聯(lián)合通過信號發(fā)射功率比對數(shù)據(jù)和導頻環(huán)路濾波器輸出進行加權(quán)組合，復雜度最高，在低載噪比環(huán)境下精度最高。

圖10 聯(lián)合跟蹤結(jié)構(gòu)Fig.10 Combined tracking structure

2.2.2 TM61跟蹤算法

TM61跟蹤最早是由ENAC的信號處理實驗室針對于TMBOC(6,1,1/11)信號提出的。2007年，文獻[61]將其用于跟蹤CBOC(6,1,1/11)信號。2019年，文獻[51]研究了TM61跟蹤在QMBOC調(diào)制信號上的應用。

由圖11可知，TM61算法中延遲鎖定環(huán)P(即時)支路將BOC(1,1)分量與接收信號進行相關(guān)運算，利用BOC(1,1)信號的高功率特性獲取接收信號的數(shù)據(jù)信息;E(超前)和L(滯后)支路采用接收信號與本地BOC(6,1)信號進行互相關(guān)，利用BOC(6,1)信號的高精度特性對接收信號進行跟蹤。該算法充分利用了BOC(1,1)信號和BOC(6,1)信號的特性，抗多徑能力強。但是該算法有兩個缺點，一是TM61跟蹤是僅適用于BOC(1,1)信號和BOC(6,1)信號復用調(diào)制的信號，具有一定局限性；二是由于該方法并沒有充分利用MBOC信號的能量，跟蹤精度低于使用復用信號的相關(guān)結(jié)果。

圖11 TM61算法跟蹤環(huán)路Fig.11 TM61 tracking loop structure

2.2.3 雙重估計技術(shù)

由式(4)、式(5)可知，BOC調(diào)制的原理是在原有BPSK調(diào)制的基礎(chǔ)上，加上一個二進制副載波對BPSK信號進行二次擴頻。針對該類信號，2008年，文獻[66]提出了雙重估計技術(shù)(Dual Estimate Technology，DET),分別使用碼環(huán)和副載波環(huán)對偽碼和副載波進行跟蹤。2018年，文獻[67]提出了一種基于改善DET的多徑抑制算法，設(shè)計了一種DET雙環(huán)交互互相關(guān)函數(shù)，多徑抑制性能高于傳統(tǒng)DET算法。

DET方法的原理是將本地偽碼信號和本地副載波信號分離，分別進行相關(guān)。從偽碼層面看，相關(guān)函數(shù)與BPSK偽碼相關(guān)函數(shù)形狀類似，具有單峰性，但是精度較低。從副載波層面看，精度較高，但是具有模糊度。通過兩者的非線性組合可以取得精度較高且無模糊度的結(jié)果，DET的性能由副載波跟蹤精度決定。目前,副載波跟蹤方式可以分為兩種：一種是延遲鎖定環(huán)跟蹤(Subcarrier Delay Loc-ked Loops，SDLL)方式，該方式計算復雜度較低，要求的硬件資源少，考慮到應用實現(xiàn)難度，一般用于偽碼速率較低或BOC階數(shù)較低的信號;另一種是鎖相環(huán)跟蹤(Subcarrier Phase Locked loops，SPLL)方式，對軟硬件要求較高，在跟蹤精度和抗多徑性能上相對于SDLL具有優(yōu)勢，一般適用于偽碼速率較高或BOC階數(shù)較高的信號。

3 結(jié)論與展望

目前,GNSS多徑抑制基帶處理算法相關(guān)研究如圖12所示，按信號體制可以分為傳統(tǒng)多徑抑制算法和新體制信號多徑抑制算法。傳統(tǒng)多徑抑制算法主要針對于BPSK信號，分為非參量式和參量式，其中非參量式算法通過改進鑒相器，降低其對于多徑的敏感性，提高抗多徑性能,該類方法環(huán)路簡單，計算量較?。粎⒘渴剿惴ㄐ枰烙嫸鄰叫盘柕奶卣?，跟蹤精度較高，計算量較大。新體制信號多徑抑制算法針對于BOC、MBOC信號，按實現(xiàn)方式分為兩類：第一類是基于去模糊的多徑抑制算法，這類方法適用范圍較廣，首先去除信號多峰性，然后改善傳統(tǒng)算法的多徑抑制性能，進一步提高跟蹤精度；第二類是針對特殊信號結(jié)構(gòu)的處理算法，此類方法根據(jù)不同的信號結(jié)構(gòu)進行設(shè)計，提高信號在跟蹤過程中的利用功率，從而提高跟蹤精度，相較于第一類方法精度更高，但適用范圍更窄。

圖12 多徑抑制基帶處理算法框架Fig.12 Framework of multipath mitigation baseband processing algorithm

由上述分析可知，GNSS多徑抑制基帶處理算法相關(guān)研究面對的挑戰(zhàn)主要有兩方面：一方面是越來越復雜的應用場景對多徑抑制算法魯棒性、快速性、準確性提出更高的要求;另一方面是不斷發(fā)展的衛(wèi)星信號體制使得多徑抑制算法需要向針對性、高效性的方向發(fā)展。隨著衛(wèi)星信號設(shè)計和信號處理技術(shù)的不斷發(fā)展，多徑抑制基帶處理算法也需要隨之不斷進步，從而進一步提升衛(wèi)星導航系統(tǒng)服務性能和服務質(zhì)量。在此基礎(chǔ)上，對GNSS多徑抑制基帶處理算法研究前景進行展望：

1)改善算法性能：

①提高多徑抑制算法的性能，提高Double-Delta、MEDLL、ELS等算法的跟蹤精度和計算效率，合理改進算法使其滿足應用需求。對于新體制信號,一方面需要加強去模糊算法的研究，使其在達到無模糊的同時保留原信號的良好特性,如主峰寬度、信號強度等，從而進一步提升多徑抑制性能;另一方面,針對信號結(jié)構(gòu)進行新型跟蹤算法的研究，如數(shù)據(jù)/導頻聯(lián)合跟蹤、TM61跟蹤算法等，提高新體制信號的能量利用率，發(fā)揮新體制信號的優(yōu)秀性能。

②改善特殊多徑的抑制性能，一是加強短多徑抑制算法研究，隨著衛(wèi)星應用越來越廣泛，特別在復雜城市環(huán)境中，多徑信號往往是近距離和短時延的，使得短多徑抑制技術(shù)具有較高的研究價值；二是加強高階BOC信號多徑抑制算法研究，目前BOC多徑抑制技術(shù)多集中于低階BOC信號，高階BOC信號的研究相對較少，高階BOC信號的帶寬較大，其自相關(guān)函數(shù)的相關(guān)峰窄而且非常密集，有的甚至存在一定的非對稱性，給跟蹤和多徑抑制帶來挑戰(zhàn)，如Galileo E1頻點和北斗B3頻點的BOC(12,2.5)等。

2)多源信息融合：結(jié)合衛(wèi)星導航的發(fā)展趨勢，將信息融合技術(shù)應用于基帶處理算法解決多徑問題，例如加強矢量跟蹤技術(shù)的研究，矢量跟蹤算法的原理是將各通道的鑒相器誤差及其他輔助信息(如慣導)輸入到導航濾波器，再生成環(huán)路控制量調(diào)節(jié)本地NCO，在高動態(tài)、遮擋和空曠的環(huán)境下，相較于標量跟蹤算法具有更高的跟蹤精度和多徑抑制能力，目前應用主要受結(jié)構(gòu)復雜和計算量大的限制，可以從濾波算法角度出發(fā)進行改善。

3)算法應用優(yōu)化：加強新播發(fā)信號多徑抑制算法的應用研究，通過改進與新播發(fā)信號結(jié)構(gòu)類似的多徑抑制算法，快速地實現(xiàn)新信號的精確跟蹤，提高算法的適用性，如早期北斗B1C信號的多徑抑制算法的研究可以借鑒GPS的TMBOC信號和Galileo的CBOC信號，加快產(chǎn)品化進程。