基于北斗三頻寬巷組合的RTK 單歷元定位方法

基于北斗三頻寬巷組合的RTK單歷元定位方法

張超1,聞道秋1,潘樹國2,高旺1,張小國2

(1.東南大學交通學院，江蘇南京210096；2.東南大學儀器科學與工程學院，江蘇南京210096)

摘要：為發揮北斗三頻超寬巷或寬巷模糊度易于固定的優勢，提出一種利用寬巷觀測值的北斗單歷元定位方法。首先使用載波偽距組合求解最容易固定的(0,-1,1)組合模糊度，其次基于模糊度已固定的(0,-1,1)組合觀測值通過CIR方法單歷元求解第二個寬巷模糊度(1,-1,0)組合模糊度。(1,-1,0)組合觀測值噪聲和電離層延遲綜合影響最優，因此基于(1,-1,0)組合觀測值進行單歷元坐標解算。實驗結果表明，文中提出的單歷元定位方法，能夠依據單歷元數據最終成功解算出地面點坐標，統計精度為平面2.2 cm，高程5.4 cm，在滿足部分領域應用要求的情況下提高效率和抗干擾能力。

關鍵詞：北斗三頻；寬巷；單歷元定位；RTK；短基線

中圖分類號：P228.4

收稿日期：2014-09-25

基金項目：國家科技支撐計劃課題(2013BAJ13B01)

作者簡介：張超(1989-)，男，碩士研究生.

Single epoch algorithm of short baseline RTK based on triple-frequency wide-lane combination of BDS

ZHANG Chao1,WEN Dao-qiu1,PAN Shu-guo2,GAO Wang1

(1.School of Transportation，Southeast University，Nanjing 210096，China;2.School of Instrument Science and Engineering，Southeast University，Nanjing 210096，China)

Abstract：As for the advantage that triple-frequency extra-wide-lane or wide-line ambiguities are easily fixed，a single epoch algorithm of short baseline RTK based on triple-frequency wide-lane combination is proposed.First the extra-wide-lane ambiguities(0,-1,1) are resolved by carrier-pseudo combination method，then the wide-lane ambiguities(1,-1,0) by CIR(Cascade Integer Resolution) with the fixed ambiguities.The influence of noise and ionosphere delay of the wide-lane ambiguities(0,-1,1) are smallest，so this combination is chozen to resolve the coordinates based on single epoch data.Experimental results show that the single epoch algorithm proposed is feasible to calculate the coordinates，with the statistical positioning error in plane 2 cm，and the error in up direction about 5cm.

Key words：triple frequency of BDS；wide-lane；single epoch algorithm；RTK；short baseline

北斗衛星導航系統(BeiDou Navigation Satellite System，BDS)是我國獨立發展、自主運行并與其他衛星導航系統兼容的衛星導航定位系統，并于2012-12-27正式向亞太大部分地區提供連續無源定位、導航、授時等服務[1]。BDS系統衛星全部播發三頻信號，是目前唯一的全三頻工作的衛星導航系統。多個頻點的觀測數據具有獨特的優勢，通過不同頻率的線性組合，可以得到一系列波長較長的超寬巷或寬巷組合。結合偽距觀測值組合，可以快速解算超寬巷或寬巷組合觀測值模糊度[4-6]。模糊度已固定的超寬巷或寬巷觀測值即可視為高精度的“偽距”觀測量，并應用于高精度快速定位。

實時動態差分技術(Real-Time Kinematic，RTK)，是一種常用的GNSS測量方法。采用實時載波相位差分技術，可以在十秒左右甚至更短的時間內固定模糊度，為用戶提供厘米級的定位精度。但是，在強遮擋、可視衛星自身分布均勻性差以及電離層活躍等惡劣觀測條件下，RTK需要較長的時間才能固定，甚至不能固定，嚴重影響RTK的作業效率。另外，傳統的RTK算法采用載波、偽距列出的誤差方程組一般無法單歷元將模糊度固定，因此其定位時間無法達到真正的瞬時?？紤]到在一些領域，如水深測量的平面定位、高精度交通導航、道路斷面測量等并不要求常規RTK平面1～2 cm、高程優于3～4 cm的定位精度，如果能適當放寬定位精度，來換取定位所需時間的大幅縮短，不僅能夠滿足工程應用，還可以提高工作效率、適應比較惡劣的觀測環境。

1北斗觀測值及其線性組合

考慮到雙差電離層、對流層殘差的影響，北斗雙差載波觀方程和雙差偽距觀測方程[3-4]：

(1)

(2)

經過組合的載波和偽距雙差觀測方程：

(3)

(4)

其中，i,j,k為任意整數；α,β,γ為任意實數；Δφi(i=1,2,3)為第i個頻率上的雙差載波相位觀測值；ΔPi(i=1,2,3)為第i個頻率上的雙差偽距觀測值；λi(i=1,2,3)和ΔNi分別為頻率i上的載波波長和雙差模糊度；Δρ，ΔT，ΔI分別為雙差衛地距、雙差對流層延遲值和北斗B1頻點上的雙差電離層延遲值； f1，f2，f3分別為北斗系統的三個載波頻率，頻率為f1=1 561.098 MHz，f2=1 207.140 MHz，f3=1 268.520 MHz。

組合的載波觀測值與組合的偽距觀測值精度：

(5)

其中，μ(i,j,k)，μ(α,β,γ)為組合載波觀測值噪聲放大因子。對北斗三個頻率的觀測值進行線性組合，能夠構造出一系列具有長波特性的寬巷、超寬巷組合，其模糊度能夠得到快速固定[3]。

2三頻寬巷組合單歷元定位模型

首先，選擇一個最容易固定的北斗超寬巷組合，實現單歷元模糊度固定；然后，將上一步固定模糊度的超寬巷組合觀測值作為精度較高的觀測量，與第二個寬巷觀測值構成無幾何模型，實現第二個寬巷/超寬巷模糊度的單歷元固定；最后，采用噪聲較小的第二個寬巷觀測方程進行單歷元定位[7]。

2.1超寬巷/寬巷模糊度固定

1)求解寬巷組合(0,-1,1)模糊度，一般按式(6)進行求解[8]：

(6)

式中：ΔN(0,-1,1)為寬巷組合模糊度，假設載波噪聲σΔφ =0.005 m，σΔP =0.6 m，計算的解算中誤差為0.092周，固定成功率為100%，表明可以實現單歷元固定。但是由于(0,-1,1)組合觀測值噪聲較大(見表1)[3]，不適合直接用于定位，因此需要解算第二個對應觀測值噪聲較小的寬巷模糊度。

表1　兩種北斗寬巷載波組合

2)參照文獻[10]，選擇第二個噪聲放大因子較小的(1，-1，0)觀測值組合(見表1)與(0,-1,1)組成無幾何模型，根據式(7)求解寬巷組合(1,-1,0)模糊度，并用于其后的單歷元定位解算。

(7)

式中：η(i,j,k)，Δεφ(i,j,k)分別為載波組合觀測值的電離層延遲因子和組合雙差載波觀測值噪聲。

(8)

根據式(8)分析得知，當電離層延遲為10cm，載波噪聲為0.5cm時，模糊度解算成功率為99.81%，因此使用CIR方法求解第二個寬巷模糊度也基本可以實現單歷元固定[10]。

2.2基于寬巷組合(1,-1,0)的坐標解算定位解算方程

(9)通過網絡RTK區域建模內插得到雙差對流層和電離層延遲值，根據衛星高度角對方程進行定權[11-12]，

(10)

其中Pφ，Ppi為載波、偽距觀測方程權陣，ER，EP分別為參考星和非參考星高度角。式(9)所采用的寬巷組合其觀測噪聲約為L1的5.58倍，導致定位精度比常規RTK略低，但是，其帶來的優勢也是明顯的，即寬巷模糊度基本可以達到單歷元固定，從而實現單歷元定位。

3解算實驗

本文采用江蘇省北斗地基增強系統一期工程中的NJLH和NJLT兩個站點組成的基線數據進行實驗，基線長度約27.7km，取其中2014-05-06UTC時20:00～21:00的觀測數據進行實驗解算。該時段可見北斗衛星7顆，以時段內平均高度角最大的C06號衛星為參考衛星。

3.1模糊度解算結果

根據文中所提方法，對寬巷(0,-1,1)進行解算，每顆非參考衛星的解算結果如圖1、圖2所示。

圖1　(0,-1,1)組合超寬巷模糊度單歷元解算結果

圖2　(1,-1,0)組合寬巷模糊度單歷元解算結果

從圖1看出，依據單歷元觀測數據解算所得寬巷 (0,-1,1)模糊度，偏差一般在0.1周以內，最多不超過0.2周，因此，單歷元固定寬巷模糊度的方法是可行的。然后，利用CIR方法求解第二個寬巷模糊度(1,-1,0)，解算結果如圖2所示。

從圖2看出，依據CIR方法解算所得寬巷(1,-1,0)模糊度，偏差一般在0.3周以內，最大不超過0.4周，因此，可以實現第二個寬巷模糊度的固定。

3.2定位解算結果

對于3 600個歷元分別進行單歷元解算，具體解算結果見圖3，從圖3中可以看出，單歷元解算結果，在N和E方向誤差基本保持在4cm之內，U方向誤差基本保持在10cm之內。對結果進行統計，中誤差為σN=0.016 m,σE=0.014 m,σU=0.054 m。由于北斗衛星南北分布均勻性比較差，導致N方向定位精度略低于E方向；而U方向的定位精度顯著低于N，E方向，這主要是高程方向幾何強度弱于平面方向以及受殘余大氣誤差的影響。

圖3　寬巷組合單歷元定位偏差(N、E、U方向)

4結論

本文提出的基于寬巷觀測值的北斗單歷元定位方法，成功的利用了超寬巷組合(0,-1,1)模糊度容易固定的優勢實現了超寬巷組合的單歷元固定。通過選用觀測值噪聲和電離層延遲綜合影響最優的(1,-1,0)組合，并基于已固定模糊度的組合觀測(0,-1,1)值，采用用CIR方法，實現了單歷元解算寬巷模糊度組合(1,-1,0)模糊度。最后基于(1,-1,0)組合觀測值成功的進行了單歷元坐標解算。實驗結果表明，提出的單歷元定位方法，能夠依據單歷元數據最終成功解算出地面點坐標，統計精度為平面2.2cm，高程5.4cm。此定位精度略低于傳統RTK定位精度，但基于瞬時定位、無初始化時間；且模糊度基于單星解算，對高遮擋等部分觀測條件惡劣的環境下有較強的適應能力。在某些對精度要求不是很高，對效率十分敏感的領域，如水深測量的平面定位、高精度交通導航，有較重要的應用價值。

參考文獻：

[1]王家偉，范曉偉.中國北斗系統正式服務亞太地區[J].地理教學，2013(4)：4-6.

[2]李博峰，沈云中，周澤波.中長基線三頻GNSS模糊度的快速算法[J].測繪學報，2009，38(4)：296-301.

[3]李金龍.GNSS三頻精密定位數據處理方法研究[D].鄭州:解放軍信息工程大學,2011.

[4]RICHERTT，El-SHEIMYN.Optimallinearcombinationsoftriplefrequencycarrierphasedatafromfutureglobalnavigationsatellitesystems[J].GPSSolutions，2007，11(1)：11-19.

[5]CHUFY,YANGM，WUJ.Triple-frequencyGPSambiguityresolutionforpostprocessedmedium-rangebaselinedetermination：aphase-onlymethod[J].JournaloftheChineseInstituteofEngineers，2014，37(3):279-287.

[6]FENGYM，RIZOSC，HIGGINSM.MultiplecarrierambiguityresolutionandperformancebenefitsforRTKandPPPpositioningservicesinregionalareas[C]//ProceedingsofIONGNSS20thinternationaltechnicalmeetingofthesatellitedivision.2007：25-28.

[7]五岳，邱蕾.網絡RTK模式下多頻載波相位觀測值解算整周模糊度[J].測繪工程，2013,22(4):1-4.

[8]WANGK，ROTHACHERM.Ambiguityresolutionfortriple-frequencygeometry-freeandionosphere-freecombinationtestedwithrealdata[J].JournalofGeodesy，2013，87(6)：539-553.

[9]FORSSELLB，MARTIN-NEIRAM，HARRISZRA.CarrierphaseambiguityresolutioninGNSS-2[C]//Proceedingsofthe10thInternationalTechnicalMeetingoftheSatelliteDivisionofTheInstituteofNavigation(IONGPS1997).1997：1727-1736.

[10]王東會,彭競,王飛雪.北斗信號體制下三頻CIR法模糊度解算方法研究[J].全球定位系統，2009(6):13-16.

[11]高成發，趙毅，萬德鈞.GPS載波定位中雙差觀測值權的合理確定[J].測繪科學，2005,30(3):28-32.

[12]甘雨，隋立芬，張鶴.利用接收機鐘差實時濾波提高GNSS精度[J].測繪工程，2014,23(5):6-8.家偉，范曉偉.中國北斗系統正式服務亞太地區[J].地理教學，2013(4)：4-6.

[2]李博峰，沈云中，周澤波.中長基線三頻GNSS模糊度的快速算法[J].測繪學報，2009，38(4)：296-301.

[3]李金龍.GNSS三頻精密定位數據處理方法研究[D].鄭州:解放軍信息工程大學,2011.

[4]RICHERTT，El-SHEIMYN.Optimallinearcombinationsoftriplefrequencycarrierphasedatafromfutureglobalnavigationsatellitesystems[J].GPSSolutions，2007，11(1)：11-19.

[5]CHUFY,YANGM，WUJ.Triple-frequencyGPSambiguityresolutionforpostprocessedmedium-rangebaselinedetermination：aphase-onlymethod[J].JournaloftheChineseInstituteofEngineers，2014，37(3):279-287.

[6]FENGYM，RIZOSC，HIGGINSM.MultiplecarrierambiguityresolutionandperformancebenefitsforRTKandPPPpositioningservicesinregionalareas[C]//ProceedingsofIONGNSS20thinternationaltechnicalmeetingofthesatellitedivision.2007：25-28.

[7]五岳，邱蕾.網絡RTK模式下多頻載波相位觀測值解算整周模糊度[J].測繪工程，2013,22(4):1-4.

[8]WANGK，ROTHACHERM.Ambiguityresolutionfortriple-frequencygeometry-freeandionosphere-freecombinationtestedwithrealdata[J].JournalofGeodesy，2013，87(6)：539-553.

[9]FORSSELLB，MARTIN-NEIRAM，HARRISZRA.CarrierphaseambiguityresolutioninGNSS-2[C]//Proceedingsofthe10thInternationalTechnicalMeetingoftheSatelliteDivisionofTheInstituteofNavigation(IONGPS1997).1997：1727-1736.

[10]王東會,彭競,王飛雪.北斗信號體制下三頻CIR法模糊度解算方法研究[J].全球定位系統，2009(6):13-16.

[11]高成發，趙毅，萬德鈞.GPS載波定位中雙差觀測值權的合理確定[J].測繪科學，2005,30(3):28-32.

[12]甘雨，隋立芬，張鶴.利用接收機鐘差實時濾波提高GNSS精度[J].測繪工程，2014,23(5):6-8.

[責任編輯：李銘娜]