基于偽距相位和STPIR組合的北斗三頻周跳探測與修復

李 迪,柴洪洲,潘宗鵬,王華潤

(信息工程大學，河南 鄭州 450001)

基于偽距相位和STPIR組合的北斗三頻周跳探測與修復

李 迪,柴洪洲,潘宗鵬,王華潤

(信息工程大學，河南 鄭州 450001)

GNSS周跳探測中，電離層殘差法的適用性受數據采樣間隔的影響較大，同時聯合其它組合觀測量進行周跳修復時，周跳修復方程組易出現病態解。針對這些問題，文中提出一種可靠的北斗三頻周跳探測與修復算法，通過構造北斗三頻電離層殘差組合觀測量，進行二階歷元間差分，基于三頻偽距相位組合優選理論，選取適用于北斗三頻數據的偽距相位組合，結合兩種組合觀測量，優選條件數較小的組合系數矩陣進行周跳修復，最后通過北斗三頻實測數據驗證，結果表明：在數據采樣間隔較大的情況下，利用構建的三個組合觀測量可以探測出北斗三頻原始數據中的所有周跳，具有很好的修復效果。

周跳探測與修復；北斗三頻;STPIR;偽距相位組合

北斗衛星導航系統(BeiDou Navigation Satellite System,BDS)是我國正在實施的自主發展、獨立運行的全球衛星導航系統，計劃2020年左右建成覆蓋全球的衛星導航系統，屆時將為全球北斗用戶提供定位、導航和授時一體化服務。同時，BDS是整星座播發三頻信號的衛星導航系統，更多頻點的信號不僅增加了觀測量，為用戶提供了更多的選擇，而且為數據預處理，特別是周跳的探測和可靠修復帶來了更多的機遇。

周跳探測與修復是衛星導航定位數據預處理的核心組成部分，準確的探測與修復周跳可以提高導航定位的精度和可靠性[1]。電離層殘差法[2-4]和偽距相位組合法[5-6]均是較為常用的周跳探測方法。綜合電離層殘差法和MW組合，可以實現對三頻數據的周跳探測和修復[3]，但當采樣間隔較大時，對小周跳的探測能力將受到限制；范麗紅對電離層殘差組合進行二階歷元間差分，解決了采樣間隔過大對電離層殘差法的影響[4]，然而其所選組合觀測量僅能對雙頻數據進行周跳探測和修復；三頻偽距相位組合法在特定的采樣間隔下對周跳的探測和修復效果較好[5]，然而其組合觀測量選取標準、探測閾值條件和修復準測仍需進一步研究；在聯合三頻偽距相位組合法和三頻無幾何相位組合法進行周跳修復時[6]，往往由于系數矩陣條件數過大，導致方程組出現病態解問題。

本文根據電離層殘差法和偽距相位組合法的原理和特點，提出一種可靠的北斗三頻周跳探測與修復方法，通過構造基于北斗三頻的STPIR(second-order time-difference phase ionospheric residual)組合觀測量，并推導其探測閾值條件和修復準則，聯合兩個偽距相位組合(0,-1,1)和(-3,6,-2)進行周跳探測與修復，通過實例驗證，可以有效地探測和修復北斗三頻數據周跳。

1 周跳探測與修復的原理

1.1 三頻偽距相位組合

三頻偽距相位組合是基于三頻偽距和載波相位組合觀測值進行差分構造的組合觀測量，對于原始觀測值，忽略多路徑效應和觀測噪聲等影響的偽距和載波相位的非差觀測方程表示為[7]

(1)

(2)

式中：i代表不同頻率，ρ為包含對流層延遲和衛星鐘差的站星間幾何距離，I為載波頻率f1對應的電離層延遲誤差，Pi和φi分別為對應頻率的偽距和相位觀測值，λi為對應頻率的載波波長，qi為對應頻率相對于載波頻率f1上的電離層延遲放大系數，dPi和dφi分別為偽距和相位觀測量中的硬件延遲，Ni為整周模糊度。

根據多頻理論[7]，三頻偽距和載波相位的組合觀測方程表示為

(3)

(4)

式(3)和式(4)作差可得

(5)

再對Nαβγ進行歷元間差分，因為硬件延遲隨時間變化比較緩慢，可以忽略：

(6)

當歷元間的電離層延遲變化及其系數均較小的時候，式子右端第三項可以忽略。因此，偽距相位組合觀測量周跳估值及其標準差為

(7)

(8)

式中：σφ和σP分別代表載波相位和偽距觀測值的觀測精度。本文采用多數文獻設定的標準，取σφ=0.01周，σP=0.3m[7]。以4倍標準差(99.99%的置信水平)為周跳探測閾值，即偽距相位組合周跳探測的條件為

(9)

(10)

1.2 三頻STPIR組合

結合北斗三頻信號的實際情況，構造北斗三頻STPIR組合，三頻載波相位電離層殘差組合觀測量(phaseionosphericresidual)為

(11)

對式(11)進行歷元間的差分，即得到電離層殘差法的周跳檢驗量：

(12)

(13)

在二階歷元間差分后的周跳檢驗量中，電離層殘差影響明顯小于一階歷元間差分值，始終在0附近波動[4]，從而更有利于周跳的探測。

假設三個頻點的載波相位觀測中誤差均為mφ=σφ，根據誤差傳播定律，三頻STPIR組合周跳檢驗量的中誤差可以表示為

(14)

以4倍檢驗量中誤差為周跳探測閾值，即三頻STPIR組合探測周跳的條件表達式為

(15)

由式(13)可知，STPIR組合周跳值與各頻點的周跳值的關系式為

(16)

1.3 周跳探測組合觀測量的選取

對于三頻偽距相位組合，較優的組合系數應該滿足[9]：①組合波長較長以減少歷元間多路徑變化的影響；②ΔI的系數較小以減小歷元間電離層延遲變化的影響；③周跳估值標準差較小以確保組合觀測量對周跳確定具有較高成功率。

對于條件2而言，式(6)中記電離層殘差ΔI的系數為κ，則：

(17)

表1 較優的偽距相位組合

由于有多種偽距相位組合方式，在探測出周跳進行修復時，需要考慮兩個偽距相位組合系數與STPIR組合系數構成的系數矩陣的條件數。如果條件數過大，會引起方程組病態，組合量的微小變化將會引起解的不穩定[10]。對不同偽距相位組合方式與STPIR組合系數構成的系數矩陣2-范數條件數如表2所示(部分)。

表2 不同組合系數矩陣2-范數條件數

從表2可以看出，由偽距相位組合(0,-1,1)和(-3,6,-2)同電離層殘差組合構成的系數矩陣條件數最小，所以本文選取這三個組合系數構建周跳探測方程式。同時根據組合觀測量歷元間的變化聯立方程組，可以解算三個載波上分別產生的原始周跳值[8,10]：

AX=L.

(18)

由于系數矩陣A可逆，可由X=A-1L唯一確定原始三個載波相位觀測值中的周跳值。

2 數據實驗與分析

為了進一步分析本文方法的有效性，利用一組實測的北斗三頻數據進行實驗，實驗數據為2016-05-05 cutb站的北斗三頻靜態數據，采樣間隔為 30 s，選擇C03號衛星，共觀測2880個歷元。原始數據中沒有周跳和粗差，使用本文選取的三個組合觀測量進行數據檢驗，結果如圖1所示，三個組合觀測量均未超過周跳探測閾值(圖中虛線代表閾值)。

圖1 三個組合量周跳數據檢驗

2.1 周跳探測

為了研究上述三個組合觀測量的周跳探測能力及修復的效果，在實驗數據的第歷元間每隔200個歷元分別加入不敏感周跳組合(1,-1,-1)、(2,2,3)、(1,3,2)，探測結果如圖2～圖4所示。采用相同的方法，在實驗數據中添加一般性的周跳組合(10,-5,8),探測結果如圖5所示。

圖2 不敏感周跳組合(1,-1,-1)

圖3 不敏感周跳組合(2,1,2)

圖4 不敏感周跳組合(1,3,2)

圖5 一般性周跳組合(10,-5,8)

由圖2～圖4可知，偽距相位組合(0,-1,1)對在第二、三頻點上發生的相同大小的周跳組合(1,-1,-1)不敏感，但是偽距相位組合(-3,6,-2)和STPIR組合對該周跳組合均可以有效探測；偽距相位組合(-3,6,-2)對滿足-3ΔN1+6ΔN2-2ΔN3=0的周跳組合(2,2,3)不敏感，但是偽距相位組合(0,-1,1)和STPIR組合對該周跳組合均可以有效探測；STPIR組合對滿足ΔN1+λ2/λ1ΔN2-2λ3/λ1ΔN3≈0的周跳組合(1,3,2)不敏感，但兩個偽距相位組合對該周跳組合均可以有效探測。從圖5可以發現，對于一般性的周跳組合(10,-5,8),3個組合觀測量均可有效探測。

由上述分析可知，對于某一組合觀測量存在不敏感性質的周跳組合，均可被其他兩個組合觀測量有效探測，由于不存在對3個組合觀測量同時不敏感的周跳組合，所以可以保證任意周跳組合在3個組合觀測量的結合使用中均可被有效探測。

2.2 周跳修復

為了驗證所選3個組合觀測量的修復周跳能力，本文分別在實驗數據的第500，1 000，1 500和2 000歷元加入模擬周跳組合(0,-1,1)，(2,2,3)，(1,3,2)和(10,-5,8)。根據其所在歷元的各組合觀測值，由式(18)進行計算各自的原始周跳組合值，結果如表3所示。

表3 周跳組合修復結果

從表3可以看出，直接從線性方程組解算的周跳值與實際周跳值存在一定的差異，其中ΔN2和ΔN3存在的差值相同，最大不超過0.2周，ΔN1存在的差值最大，但不超過0.3周。所以對解算結果直接取整，可以獲取正確的周跳修復值。

3 結束語

1)電離層殘差法是一種常用的周跳探測方法，具有對小周跳敏感的特性，但是當采樣間隔大于15 s時，電離層延遲變化較大，常規的電離層殘差法不再適用于周跳的探測工作。由于相位電離層殘差二階歷元間差分組合有效地消除了電離層高階項的影響，在采樣間隔較大時，依舊可以有效地進行周跳探測。本文根據三頻無幾何相位組合(1,1,-2)，提出了一種組合系數為的三頻STPIR組合，可以更好地用于三頻數據的周跳探測與修復。

2)在進行周跳修復時，為了防止方程組出現病態的問題，篩選系數矩陣條件數較小的偽距相位組合(0,-1,1)和(-3,6,-2)，與三頻STPIR組合構成北斗三頻數據的周跳探測與修復組合。該組合能夠消除電離層的影響，有效探測出北斗三頻數據中包括不敏感周跳在內的所有周跳組合，且修復后的周跳值與實際周跳值之間的差異不大于0.3周，通過直接取整可以獲取正確的周跳修復值。

[1] 姚一飛,高井祥,王堅,等.北斗三頻載波觀測值的周跳實時探測與修復[J].中國礦業大學學報,2014,43(6):1140-1148.

[2] 趙曉陽,黃張裕,魯卓康,等.基于電離層殘差法探測與修復周跳的特性分析[J].海洋測繪,2015,35(2):65-68.

[3] 陶庭葉,何偉,高飛,等.綜合電離層殘差和超寬巷探測和修復北斗周跳[J].中國慣性技術學報,2015,23(1):54-58.

[4] 范麗紅,王利,張明,等.基于MW與STPIR組合的周跳探測與修復方法研究[J].武漢大學學報(信息科學版),2015,40(6):790-794.

[5] 謝愷,柴洪洲,王敏.不同采樣間隔下的三頻周跳探測與修復算法研究[J].大地測量與地球動力學,2014,34(1):139-143.

[6] 黃令勇,宋力杰,王琰,等.北斗三頻無幾何相位組合周跳探測與修復[J].測繪學報,2012,41(5):763-768.

[7] 李林陽,呂志平,崔陽,等.偽距相位和無幾何相位組合探測與修復多頻周跳的比較[J].大地測量與地球動力學,2015,35(3):396-400.

[8] CAI Changsheng,LIU Zhizhao,XIA Pengfei,et al.Cycle Slip Detection and Repair for Undifferenced GPS Observations Under High Ionospheric Activity[J].GPS Solution,2012,17(2):247-260.

[9] 李金龍,楊元喜,徐軍毅,等.基于偽距相位組合實時探測與修復GNSS三頻非差觀測數據周跳[J].測繪學報,2011,40(6):717-722.

[10] 王華潤,柴洪洲,謝愷.北斗三頻無幾何、消電離層組合周跳探測方法研究[J].大地測量與地球動力學,2015,35(3):406-410.

[責任編輯：李銘娜]

BDS triple-frequency cycle-slip detection and repair based on Code-Phase and STPIR

LI Di,CHAI Hongzhou,PAN Zongpeng,WANG Huarun

(Information Engineering University,Zhengzhou 450001,China)

The applicability of the ionosphere residual method is influenced by the observation interval in the GNSS cycle-slip detection.The ill-posed solution often occurs in the cycle-slip repair equations combined with other combinations.To solve those problems,a reliable cycle-slip detection method is proposed for BDS triple-frequency measurement.The method constructs a second-order time-difference phase ionosphere residual combination measurement,and then,selects a code-phase combination suitable for BDS triple-frequency.Combining with two kinds of combination,this paper repairs the cycle-slip with combination coefficient matrix which has small condition number.Finally,BDS triple-frequency data have been used to validate the algorithm,and the experimental result indicates that the method proposed in this paper can effectively detect and repair all cycle-slip even in the large observation interval.

cycle-slip detection and repair;BDS triple-frequency;STPIR;code-phase combination

引用著錄：李迪,柴洪洲,潘宗鵬，等.基于偽距相位和STPIR組合的北斗三頻周跳探測與修復[J].測繪工程，2017，26(5)：71-75，80.

10.19349/j.cnki.issn1006-7949.2017.05.015

2016-06-27

國家自然科學基金資助項目(41274045;41574010;41604013)

李 迪(1991-)，男，碩士研究生.

P228

A

1006-7949(2017)05-0071-05