ＧＰＳ相位觀測值中周跳的探測與修復

摘 要：對于初始化載波相位整周模糊度的GPS定位而言，如果在整周模糊度初始化過程中發(fā)生載波相位周跳，將嚴重影響整周模糊度的確定。針對GPS載波相位定位中出現(xiàn)的周跳問題，提出了一種新的周跳的檢測方法，首先用多項式模型的卡爾曼濾波對載波相位測量值進行預處理，然后再用奇偶矢量法探測和修復周跳。通過某OEM板輸出的數(shù)據(jù)計算結(jié)果仿真，驗證了新方法可以有效地探測出一周周跳。

關(guān)鍵詞：GPS；周跳探測；載波相位；卡爾曼濾波；奇偶矢量法

中圖分類號：P228-4 文獻標識碼：B

文章編號：1004373X(2008)0501603

Cycle Slips Detection and Recovery for GPS Carrier Phase

SHI Xiaoli1，2，CHEN Xiaoping ZHONG Erjie2

(1.Research Institute of Electronic Science and Technology，University of Electronic Science and Technology of China，Chengdu，610054，China;

2.Applied Mathematics College，University of Electronic Science and Technology of China，Chengdu，610054，China)

Abstract：In GPS positioning，the cycle slip occurred in the period of integer ambiguity initializating must be carefully removed.This paper discusses cycle slip problem and its solving method.New methods for enhancing detection capability are presented.Firstly，using Kalman filter in the model of the polynomial to pretreat the carrier phase，and then a parity vector test method is used to detect cycle slips.A test is carried out and implies that the new method can detect and repair small cycle slip.

Keywords：GPS;cycle slips;carrier phase;Kalman filter;parity vector

1 引 言

在利用GPS載波相位進行定位時，要達到高精度，必須保證觀測數(shù)據(jù)中無載波相位周跳，因此正確檢測與修復周跳是GPS數(shù)據(jù)處理中的一個關(guān)鍵問題^[1]。周跳探測與修復的基本思路是通過合理的數(shù)據(jù)處理，得出精確反映周跳變化的檢測量序列，從檢測量序列中探測出周跳發(fā)生的位置和大小并修復，最后將周跳改正后的數(shù)據(jù)序列參與解算。目前已有的檢測周跳的方法主要有差分法、多普勒頻移法、偽距相位組合法^[1]、電離層殘差法^[2]、卡爾曼濾波法^[3]和站間或站星間差分法。這幾種方法各有特點，但也都有其局限性：差分法主要用于探測較大的周跳，因為差分法在放大周跳的同時也放大了噪聲；多普勒頻移探測周跳時，必須保證至少前四個時元的載波相位觀測值沒有周跳，繼而用他們來探測和修復第五個時元的載波相位觀測值的周跳，另外以上兩種方法在動態(tài)情況下不能探測出小周跳；偽距相位組合法周跳的探測能力倚賴于偽距觀測量的精度，一般需要接收機能夠輸出P碼偽距，對單頻接收機不適用；電離層殘差法是目前效果最好的周跳探測方法之一，但也不能應用于單頻接收機^[5]；卡爾曼濾波法在當運動載體機動加速度較大時不能有效探測和修復周跳。

本文提出一種新方法，基本思想是：首先采用三階多項式模型的卡爾曼濾波對載波相位測量值進行預處理，找出一組不含有周跳的觀測值，然后再用奇偶矢量法^[4]探測下一觀測值是否含有周跳，依次循環(huán)，完成周跳探測與修復。

2 載波相位測量值的預處理

在進行探測和修復周跳時，一些方法必須保證至少前四個歷元的觀測值不含有周跳，從而用來探測和修復第五個觀測值的周跳。本文首先利用三階多項式模型的卡爾曼濾波進行預處理，用殘差平方和檢驗前N個歷元的觀測值有無周跳，選出不含有周跳的N個歷元，然后用奇偶矢量法探測第N⁺¹個歷元是否含有周跳，依次循環(huán)。

首先建立GPS載波相位測量的動態(tài)模型和測量方程分別為：

3 奇偶矢量法探測周跳原理

3.1 奇偶空間矢量的形成

假設沒有測量誤差的情況下，GPS載波相位的測量模型為：

3.2 周跳的探測

由于奇偶矢量[WTHX]p[WTBX]直接反映了觀測誤差的信息，基于奇偶矢量可構(gòu)造檢驗統(tǒng)計量，進行周跳的檢測。假設在上述中判斷出了一組載波相位觀測值在N個歷元內(nèi)沒有周跳，下面來判斷第N+1個歷元是否含有周跳。

根據(jù)這N+1個歷元所組成的載波相位測量模型計算出奇偶空間矩陣[WTHX]Q[WTBX]2，將奇偶空間矢量投影到

則檢測到第k個歷元含有周跳。

3.3 周跳的修復

探測出周跳后，去除測量模型中發(fā)生周跳的測量值，得到新的測量模型：

4 算法步驟及算例仿真實現(xiàn)

4.1 算法步驟

為了更好地說明該方法，由以上理論設計出周跳探測和修復的具體算法為：

第一步：根據(jù)建立的動態(tài)方程和測量方程用N點進行三階多項式平滑(在此取N=7)，求出狀態(tài)估計值；

第二步：計算殘差平方和C和δ，如果滿足式(6)，判斷這N歷元內(nèi)沒有周跳，如果不成立，則取下N個歷元并跳到第一步，直到式(6)滿足為止；

4.2 仿真結(jié)果

為驗證該方法的正確性，我們利用一組無周跳數(shù)據(jù)人為地增加跳值后進行了周跳的探測和修復，進行下列模擬實驗，所用數(shù)據(jù)(表1)為某OEM板輸出的數(shù)據(jù)，數(shù)據(jù)采樣間隔為1 s，表2中為人為添加跳值時探測的結(jié)果。

表1 相位觀測值(無周跳)

經(jīng)過載波相位預處理得知，根據(jù)建立的動態(tài)方程和測量方程用7點進行三階多項式平滑，判斷出第1~7個歷元中不含有周跳，然后用這7個歷元來探測下面12個歷元，如表2所示，在第8個歷元中加入1周，在第18個歷元中加入1周，可以看出，L(8)和L(18)發(fā)生了周跳，其他歷元均沒有發(fā)生周跳。利用式(15)得到周跳的估值后即可對載波相位測量值進行修復。

5 結(jié) 語

討論，并提出一種新方法，首先用多項式模型的卡爾曼濾波對載波相位測量值進行預處理，然后采用奇偶矢量法進行探測周跳。通過對某OEM板上輸出的數(shù)據(jù)進行仿真分析，結(jié)果表明該方法可以探測出小周跳，從而驗證了該方法的有效性。

參考文獻

［1］張守信.GPS衛(wèi)星測量定位理論與應用[M].長沙:國防科技大學出版社，1996.

［2］賈沛璋，吳連大.單頻GPS周跳檢測與估計算法[J].天文學報，2001，42(2):192-197.

［3］廖向前，黃順吉.GPS載波相位的周跳檢測方法[J].電子科技大學學報，1997，26(6):27-31.

［4］陳金平.GPS完善性增強研究[D].鄭州:解放軍信息工程大學測繪學院，2001.

［5］Chai Yanju，Ou Jikun .Method for Detecting and Repairing Cycle Slips in GPS Navigation[J].Transactions of Nanjing University of Aeronautics，2005:121-124.

［6］肖飛，羅斌風.GPS接收機射頻模塊的設計[J].現(xiàn)代電子技術(shù)，2006，29(12):120-121，125.

注：“本文中所涉及到的圖表、注解、公式等內(nèi)容請以PDF格式閱讀原文。”