張 鑫， 崔曉偉， 馮振明

(清華大學 電子工程系，北京 100084)

基于統計質量控制的測距故障檢測方法

張 鑫， 崔曉偉， 馮振明

(清華大學 電子工程系，北京 100084)

衛(wèi)星導航接收機在城區(qū)環(huán)境使用時，多徑信號的存在使得碼自相關函數發(fā)生畸變，從而導致測距結果與真實距離間出現較大偏差?，F有的接收機基帶多徑消除算法復雜度高，實現較為困難。而接收機自主正直性監(jiān)測(RAIM)算法的故障檢測能力有限，難以同時檢測多個測距故障。提出使用統計質量控制方法來實現對測距故障的檢測，根據質量控制圖的輸出結果調整測距信息在定位解算中的權重，從而提高定位精度。該方法能夠對每顆衛(wèi)星的碼跟蹤性能進行單獨監(jiān)測，算法實現的復雜度低。仿真結果表明：使用該方法能夠有效指示多徑信號導致的測距誤差，從而為改善城區(qū)環(huán)境下的定位精度提供了新思路。

多徑； 測距誤差； 故障檢測； 統計質量控制

0 引 言

在城區(qū)環(huán)境，來自某些方位的衛(wèi)星信號與建筑物反射信號進行疊加而形成多徑信號。接收機在跟蹤多徑信號時碼相關函數會發(fā)生畸變，導致測距出現誤差,利用含有測距誤差的觀測量進行定位會導致定位結果發(fā)生偏移。

現有的削弱多徑信號對定位結果影響的方法可劃分如下:其一是使用高精度抗多徑天線，抑制多徑信號的接收[1,2],這種方法從源頭上扼制了多徑信號接收，但使得接收機體積增大，成本增高。其二是在接收機基帶跟蹤中添加多徑信號參數估計與多徑信號消除算法[3～5],該類算法實現復雜度較高。在接收機同時跟蹤多顆衛(wèi)星信號時，如果對每路信號都執(zhí)行多徑消除算法，則計算實時性難以保證。另一種是接收機自主正直性監(jiān)測(RAIM)算法，通過對定位解算后的殘差進行一致性檢測來排除故障觀測量[6,7]。但RAIM算法只利用了偽距殘差信息，而沒有使用在接收機對導航信號進行跟蹤時獲得的信息,這使得該算法在多個測距故障同時發(fā)生時無法進行有效檢測且檢測性能對可見衛(wèi)星的星座分布十分敏感[7]。

事實上，隨著多個衛(wèi)星導航系統的建設與完善，可見衛(wèi)星數目顯著增加。在這種情況下，接收機可以直接剔除含有測距誤差的觀測量，而不會引起定位精度因子(DOP)的惡化(DOP值由可見衛(wèi)星的幾何分布決定，其值越小，表示定位性能越好)，因此，無需對多徑信號參數進行估計。另一方面，在檢測到故障觀測量后，也可不直接剔除，而是通過降低該觀測量在定位解算中的權重以削弱其對定位結果的影響，這樣則不會影響DOP值。

為此，本文提出一種新的測距故障檢測方法，使其能夠充分利用信號跟蹤過程中的信息，與RAIM算法相比有效提升了故障檢測能力。而相比于多徑消除算法,則大幅降低計算復雜度，能夠在可見星較多時滿足實時性要求。該方法從質量控制角度出發(fā)[8,9]，借助質量控制領域中的統計監(jiān)測方法實現對信號跟蹤過程的質量監(jiān)測，以判別是否存在測距故障。

1 實現思路

算法的設計思路如圖1所示。

圖1 算法基本思路Fig 1 Basic idea of algorithm

如圖1所示，本文把測距信息的提取看作一個生產過程，直達信號為正常的輸入信息,噪聲和多徑信號為不可控的輸入信息。通過對信號跟蹤內部信息監(jiān)測(提取質量控制觀測量)，使用統計質量控制方法判別測距結果是否存在異常。如果在質量控制過程中檢測到測距異常，則將其以加權的形式反饋到定位解算模塊中，即降低異常觀測量在定位中的權重(權重為零時表示舍棄該觀測量，這與RAIM算法中隔離故障觀測量的思路是一致的)。

統計質量控制方法中的常用工具為控制圖,選取合適的控制圖可以對質量控制觀測量的統計特性(均值、方差等)的變化進行快速有效的檢測。如果質量控制觀測量能夠較為準確地反映系統運行的真實情況，則利用該觀測量構造的控制圖可以對系統運行質量即測距信息提取的質量進行評估?？梢姡疚乃崴惴ǖ暮诵脑谟谠O計質量控制觀測量和建立控制圖，其實現方法將在第4章中詳細說明。

2 測距誤差模型

2.1 多徑信號模型

包含有多徑信號的接收衛(wèi)星信號表達式為

(1)

其中,A為中頻信號幅度，C(t)表示衛(wèi)星的偽隨機碼信號，fd為多普勒頻率，φ0為初始相位，αk表示第k個多徑信號相對于直達信號的衰減幅度，τk表示第k個多徑信號的延遲，φk(t)表示第k個多徑信號的相位變化函數,n(t)為噪聲項。通常假定在積分時間T內，多徑信號的相位呈線性變化，即

φk(t)=φk,0+2πfkt.

(2)

其中，fk表示多徑信號相位的變化頻率，稱為衰落頻率,衰落頻率滿足fk?1/T。

2.2 多徑信號測距誤差

使用超前—滯后碼鑒相器[6]，記ye,yl分別為超前路和滯后路的積分結果，d為相關器的間距，則鑒相器輸出為

(3)

在信號穩(wěn)態(tài)跟蹤后，鑒相函數的過零點位置指示了測距誤差。為研究不同多徑信號參數下的測距誤差大小，將多徑信號與直達信號的幅度比固定為0.5，對不同多徑參數設置下的碼跟蹤誤差進行仿真，仿真結果如圖2所示。

圖2 測距誤差仿真結果Fig 2 Simulation results of ranging error

3 統計監(jiān)測方法

控制圖是統計質量控制領域的常用工具。使用控制圖能夠檢測觀測量中均值和方差的變化。其中指數移動加權平均(EWMA)控制圖實現簡單，對檢測量的分布模型不敏感[10]，因而，本文選取EWMA控制圖作為統計監(jiān)測方法，對測距誤差進行連續(xù)監(jiān)測，通過將控制圖的輸出與控制門限進行比較實現測距異常告警。

3.1 測距誤差觀測量構造

在使用控制圖方法對接收信號質量進行統計監(jiān)測之前，需要首先構造觀測量，使其與圖2所示的測距誤差仿真結果呈現正相關，即構造的觀測量要能夠反映出由多徑信號引入的測距誤差的大小變化。

考慮到多徑信號的幅度通常小于直達信號幅度，因而碼相關函數的最大值所在位置可以作為直達信號碼相位的參考點。因此，本文通過增加相關器的個數來構造誤差觀測量。

當環(huán)路進入穩(wěn)定跟蹤狀態(tài)時，準時支路的碼相位記為p，其相關值為yp。此時超前支路的相關值記為ye=yp-0.5，滯后支路的相關值為yl=yp+0.5。引入碼相位位于p±0.1,p±0.2,p±0.3,p±0.4處的相關函數值計算模塊，將這些相關值計算結果與yp一起進行比較，找到最大相關值所對應的碼相位，將其作為誤差觀測量，記為T,數學表達式如下

T=Δx.

(4)

其中,Δx滿足yp-Δx=max{yp-Δx},Δx=0,±0.1,±0.2,±0.3,±0.4。

進一步采用仿真方式對誤差觀測量的合理性進行驗證。圖3給出了在不同參數(與圖2仿真中參數變化方式相同)下，使用式(4)所述方法構造誤差觀測量所得到的觀測結果。對比圖3和圖2，可以看到，該觀測量在不同多徑參數設置下的分布與真實碼跟蹤誤差的分布具有較好的相關性。

圖3 誤差觀測量仿真結果Fig 3 Simulation results of error observations

進一步可通過計算誤差觀測量與真實誤差之間的相關系數來評估該觀測量構造是否合理。相關系數計算方法如下

(5)

其中,x和y分別表示不同仿真樣點下，測距真實誤差輸出和誤差觀測量輸出所構成的向量。對于不同的多徑信號幅度，表1給出了相關系數的計算結果。

表1 觀測量與真實值的相關系數計算Tab 1 Correlation coefficient calculation between observations and true values

從表1可見，觀測量與真實誤差之間呈高度線性相關，隨著多徑信號幅度增大，相關性逐漸增強。事實上，當多徑信號幅度較小時，其造成的測距誤差也較小。因此，式(4)給出的觀測量可用于對多徑信號導致的測距誤差進行檢測。

3.2 EWMA控制圖建立

在第i個觀測時刻，EWMA控制圖的統計量Zi表示為[8]

Zi=λTi+(1-λ)Zi-1.

(6)

其中，λ為常數且0<λ≤1，Ti為第i個觀測時刻的觀測值，Z0通常選為監(jiān)測量無偏移時的目標值μ。對于多徑誤差監(jiān)測，當相關函數沒有發(fā)生畸變時，其觀測值應為0，故Z0=μ=0。

將式(6)展開可得

(7)

可見，EWMA統計量為當前觀測量和所有歷史觀測值的加權和，且權重隨著觀測值的滯后呈指數下降。若各個觀測值之間相互獨立且方差恒為σ2，則檢驗統計量Zi的方差可以寫為

(8)

根據式(8)可設定控制圖的判決門限,即控制上限(UCL)和控制下限(LCL)

(9)

(10)

參數L的大小表示質量控制過程對觀測值偏移大小的容忍程度,參數λ和L共同決定了EWMA控制圖方法對測距誤差的檢測性能。

4 仿真結果

首先對直達信號進行仿真分析，根據觀測量的統計結果設定控制圖參數。仿真信號采樣率為5 MHz，信噪比為-23 dB。根據誤差觀測量的統計結果，設定控制圖參數為μ=0，σ=0.049 8，λ=0.02，L=6。

向直達信號中加入兩路多徑信號，以測試控制圖方法能否對測距誤差進行有效監(jiān)測。加入的多徑信號參數如表2。

表2 多徑信號參數設置Tab 2 Multipath signal parameters settings

由于信號為仿真生成，其每個采樣點的真實碼相位是已知的。在使用環(huán)路對信號跟蹤過程中，記錄每個采樣點跟蹤得到的碼相位信息。將其與真實碼相位做差后，得到碼跟蹤誤差。碼跟蹤誤差曲線如圖4(a)所示,圖4(b)給出了檢驗統計量的控制圖輸出。其中，UCL,LCL為控制門限,當檢驗統計量超過控制門限時，表明當前接收信號存在異常。

圖4 跟蹤誤差與控制圖輸出Fig 4 Tracking error and control chart output

對比圖4(a)與(b)可以看出：控制圖輸出能夠對碼跟蹤異常進行準確指示。對于導航接收機，可以使用本文提出的方法對每個衛(wèi)星信號的碼跟蹤結果進行檢驗，若某顆衛(wèi)星的控制圖輸出超出控制門限，則表明使用該顆衛(wèi)星信號獲取的測距觀測量中可能存在誤差。此時在定位解算中，可通過剔除該顆衛(wèi)星的觀測量或者降低該顆衛(wèi)星觀測量在定位結果中的權重來提升定位精度。

5 結 論

在城區(qū)環(huán)境下的衛(wèi)星導航應用中，多徑信號的存在是導致定位精度降低的主要原因。本文使用統計質量控制方法，通過構造合適的質量控制觀測量對每顆衛(wèi)星信號的跟蹤狀態(tài)分別進行監(jiān)測，從而判別該顆衛(wèi)星的測距值是否存在異常。該方法實現簡單，復雜度低。仿真表明:使用該方法能夠對多徑信號導致的碼跟蹤誤差進行有效指示，從而為降低城區(qū)環(huán)境下由多徑信號引入的定位性能下降提供了新思路。

[1] Maqsood M,Gao S,Brown T,et al.A compact multipath mitigating ground plane for multiband GNSS antennas[J].Antennas and Propagation,2013,61(5):2775-2782.

[2] Nayak R A,Cannon M E,Wilson C,et al.Analysis of multiple GPS antennas for multipath mitigation in vehicular navigation[C]∥Institute of Navigation National Technical Meeting/Anaheim,CA,2000:26-28.

[3] Closas P,Fernandez-Prades C,Fernandez-Rubio J A.A Bayesian approach to multipath mitigation in GNSS receivers[J].Selected Topics in Signal Processing,2009,3(4):695-706.

[4] Townsend B,Fenton P,Van Dierendonck K,et al.L1 carrier phase multipath error reduction using MEDLL technology[C]∥Procee-dings of ION GPS,CA,1995:1539-1544.

[5] Phan Q H,Tan S L,McLoughlin I.GPS multipath mitigation:A nonlinear regression approach[J].GPS Solutions,2012,17(3):371-380.

[6] Parkinson B W.Global positioning system[M].Washington DC,USA:American Institute of Aeronautics and Astronautics,1996.

[7] Macabiau C,Gerfault B,Nikiforov I,et al.RAIM performance in presence of multiple range failures[C]∥Proceedings of the National Technical Meeting of the Institute of Navigation,2005:779-791.

[8] Montgomery D C.Introduction to statistical quality control[M].Hoboken:Wiley,2007.

[9] Reynolds,Marion R,Lou J Y.An evaluation of a GLR control chart for monitoring the process mean[J].Journal of Quality Technology,2010,42(3):287-310.

[10] Ingolfsson,Armann,Emanuel Sachs.Stability and sensitivity of an EWMA controller[J].Journal of Quality Technology,1993,25(4):271-287.

Ranging fault detection method based on statistical quality control

ZHANG Xin， CUI Xiao-wei， FENG Zhen-ming

(Department of Electronic Engineering，Tsinghua University，Beijing 100084，China)

When satellite navigation receiver is used in urban area,existence of multipath signal cause distortion of code autocorrelation function,resulting in a large bias between ranging result and real distance.The baseband multipath cancellation algorithm of reciever is hard to implement due to its high computational complexity.While it is difficult for detecting ability of receiver autonomous integrity monitoring (RAIM) algorithm is limited to detect multiple failures.Propose a new method to detect ranging fault based on the idea of statistical quality control.According to output of quality control chart,adjust the weight of ranging information in positioning calculation,thus increase precision of positioning.The character of code tracking for each satellite is monitored separately,complexity of algorithm implementation is low.The simulation results demonstrate that the method can effectively indicate ranging error by multipath signal,which provides a new idea to improve location precision in urban environment.

multipath; ranging error; fault detection; statistical quality control

2014—01—23

TN 967.1

A

1000—9787(2014)04—0150—04

張 鑫(1987-)，男，北京人，博士研究生，主要研究領域為衛(wèi)星導航技術。