高 猛，徐愛功，祝會忠，葛茂榮,2，楊秋實

1. 遼寧工程技術(shù)大學測繪與地理科學學院，遼寧 阜新 123000； 2. 德國地學研究中心，波茨坦 14473

?

BDS網(wǎng)絡(luò)RTK參考站三頻整周模糊度解算方法

高 猛1，徐愛功1，祝會忠1，葛茂榮1,2，楊秋實1

1. 遼寧工程技術(shù)大學測繪與地理科學學院，遼寧 阜新 123000； 2. 德國地學研究中心，波茨坦 14473

北斗衛(wèi)星導(dǎo)航系統(tǒng)是目前唯一一個全星座提供三頻觀測數(shù)據(jù)的衛(wèi)星導(dǎo)航定位系統(tǒng)，三頻觀測值有助于載波相位整周模糊度的快速、準確固定。本文提出了一種BDS網(wǎng)絡(luò)RTK參考站三頻整周模糊度解算方法。首先利用B2、B3頻率的觀測值及嚴格的模糊度固定標準確定超寬巷整周模糊度，將固定的超寬巷整周模糊度與其他寬巷整周模糊度的線性關(guān)系作為約束條件，然后估計寬巷整周模糊度、相對天頂對流層延遲誤差和電離層延遲誤差，并搜索確定寬巷整周模糊度。利用固定的寬巷整周模糊度與三頻載波相位整周模糊度的整數(shù)線性關(guān)系，將線性關(guān)系加入載波相位整周模糊度參數(shù)估計觀測模型中，然后確定載波相位整周模糊度。使用實測的CORS網(wǎng)BDS三頻觀測數(shù)據(jù)進行算法驗證，結(jié)果表明，該方法可正確有效地實現(xiàn)參考站間三頻載波相位整周模糊度的快速解算。

BDS；網(wǎng)絡(luò)RTK；三頻；整周模糊度；相對天頂對流層延遲誤差；電離層延遲誤差

北斗衛(wèi)星導(dǎo)航系統(tǒng)(BeiDou navigation satellite system，BDS)是我國正在實施的自主發(fā)展、獨立運行的全球衛(wèi)星導(dǎo)航系統(tǒng)[1]。截至2016年4月初，5顆新一代北斗導(dǎo)航衛(wèi)星準確進入預(yù)定軌道，標志著北斗衛(wèi)星導(dǎo)航系統(tǒng)由區(qū)域運行向全球拓展的建設(shè)目標邁出了堅實的一步。作為全球首個全星座播發(fā)三頻衛(wèi)星導(dǎo)航信號的衛(wèi)星導(dǎo)航系統(tǒng)，多頻觀測信息極大地改善了北斗衛(wèi)星導(dǎo)航系統(tǒng)的觀測冗余度，提高了導(dǎo)航定位的精度和效率[2-3]。北斗衛(wèi)星導(dǎo)航系統(tǒng)提供的標準服務(wù)定位精度為10 m，測速精度為0.2 m/s，授時精度為20 ns[4]。BDS標準定位服務(wù)不能滿足用戶高精度定位的需求，BDS網(wǎng)絡(luò)RTK方法是提高BDS定位精度，實現(xiàn)高精度實時動態(tài)定位服務(wù)的重要手段。

BDS網(wǎng)絡(luò)RTK的核心問題之一是參考站間載波相位整周模糊度的準確確定，只有正確固定了參考站模糊度才能得到高精度的綜合誤差或是建立高精度的誤差模型[5]。由于網(wǎng)絡(luò)RTK參考站一般相距幾十千米以上，所以參考站間電離層延遲和對流層延遲等誤差對雙差觀測值的影響大于模糊度的半個波長，即使在使用雙頻觀測數(shù)據(jù)和參考站坐標已知的情況下，整周模糊度也難以與誤差分離。

近年來，國內(nèi)外學者對網(wǎng)絡(luò)RTK參考站間模糊度進行了大量研究。文獻[6]利用卡爾曼濾波算法進行了參考站間整周模糊度的動態(tài)解算。文獻[7—8]提出了網(wǎng)絡(luò)RTK參考站間單歷元模糊度搜索方法，解決了參考站間模糊度單歷元解算問題。文獻[9]研究了參考站間雙差整周模糊度快速解算的三步法。文獻[10]研究了不固定參考站坐標情況下利用GPS雙頻觀測數(shù)據(jù)仿三頻數(shù)據(jù)的三頻模糊度快速算法。文獻[11]提出了利用電離層約束模型的網(wǎng)絡(luò)RTK模糊度解算方法。上述研究工作主要針對GPS網(wǎng)絡(luò)RTK，BDS能夠向亞太地區(qū)提供連續(xù)無源的導(dǎo)航、定位和授時服務(wù)之后，學者們已利用BDS實測數(shù)據(jù)進行網(wǎng)絡(luò)RTK參考站模糊度的研究。文獻[12]使用綜合噪聲最小的弱電離層組合構(gòu)造了一種三頻幾何模糊度解算策略，由于其使用觀測量組合而不是原始觀測量，故不可避免地放大了觀測噪聲影響，影響了參數(shù)估計精度。文獻[13]研究了基于三頻無幾何模型的BDS網(wǎng)絡(luò)RTK方法，主要實現(xiàn)了超寬巷和寬巷模糊度的準確固定，并且具有較高的解算效率，但未對BDS原始頻率模糊度進行解算。文獻[14]利用B1、B2頻率載波相位整周模糊度的線性關(guān)系，通過雙差電離層延遲誤差模型進行B1、B2載波相位整周模糊度單歷元固定，但該方法適用于電離層平穩(wěn)的情況，電離層誤差空間相關(guān)模型化較難的條件下單歷元模糊度解算的可靠性得不到保證。相對于GPS，BDS寬巷組合觀測值的噪聲較大，模糊度解算的收斂速度慢，固定成功率不高，BDS雙頻單歷元方法由于以單顆衛(wèi)星為研究對象，沒有綜合利用所有衛(wèi)星的觀測值信息。

鑒于上述情況，應(yīng)充分利用BDS提供三頻觀測信息的優(yōu)勢，進行BDS網(wǎng)絡(luò)RTK參考站三頻整周模糊度解算方法的研究。本文利用BDS三頻超寬巷長波長及三頻載波相位整周模糊度之間的整數(shù)線性關(guān)系，提出了BDS網(wǎng)絡(luò)RTK參考站三頻整周模糊度解算方法。該方法利用B2、B3頻率的觀測值得到超寬巷整周模糊度，將利用多個模糊度固定標準固定的超寬巷整周模糊度與寬巷整周模糊度的線性關(guān)系作為約束條件，避免直接利用噪聲較大的超寬巷組合觀測值對參數(shù)估計的影響，然后估計寬巷整周模糊度、相對天頂對流層延遲誤差和電離層延遲誤差，并搜索確定寬巷整周模糊度。在進行原始頻率雙差整周模糊度確定過程中也不使用寬巷組合觀測值，而是利用固定的寬巷整周模糊度與三頻載波相位整周模糊度的整數(shù)線性關(guān)系，將線性關(guān)系加入載波相位整周模糊度參數(shù)估計觀測模型中，最終可正確有效地確定載波相位整周模糊度。

1 BDS網(wǎng)絡(luò)RTK參考站三頻整周模糊度解算

雙差組合觀測值具有可消除接收機鐘差和衛(wèi)星鐘差，大大削弱軌道誤差、電離層延遲誤差、對流層延遲誤差等誤差影響的優(yōu)點，雙差載波相位模糊度保持整周特性[15]，所以，本文對BDS參考站間的載波相位雙差整周模糊度進行解算。通常情況下，BDS網(wǎng)絡(luò)RTK參考站一般都設(shè)在比較開闊的地方，且硬件設(shè)備齊全，多采用抗多路徑天線，衛(wèi)星的多路徑效應(yīng)較小，可忽略多路徑的影響，則BDS參考站間的B1、B2和B3頻率上的雙差偽距和載波相位觀測方程可表示為

(2)

利用BDS 3個頻率的載波相位觀測值進行組合，可以得到一系列具有較長波長的寬巷和超寬巷組合觀測值[16-17]，其中B2、B3頻率的超寬巷組合觀測值對應(yīng)的波長約為4.884 m，其電離層延遲誤差的影響是B1頻率電離層延遲誤差的1.6倍；B1、B2和B1、B3頻率的寬巷組合觀測值對應(yīng)的波長約為0.847 m和1.025 m，二者的電離層延遲誤差影響較接近，其載波相位和偽距組合觀測值的觀測噪聲小于B2、B3超寬巷組合觀測值的觀測噪聲。

1.1 參考站間B2、B3超寬巷整周模糊度解算

參考站間的雙差電離層延遲誤差、對流層延遲誤差和軌道誤差等誤差的殘差影響隨著站間距離的增加逐漸增大，因此，采用不受電離層延遲誤差、對流層延遲誤差和衛(wèi)星軌道誤差等誤差和幾何距離影響的M-W組合計算參考站間B2、B3雙差超寬巷整周模糊度。如果參考站A、B同步觀測衛(wèi)星p、q，偽距和載波相位觀測值的M-W組合觀測值如式(3)所示，根據(jù)式(4)可計算B2、B3雙差超寬巷整周模糊度

(3)

(4)

式(3)、式(4)中，MW為M-W組合觀測值；N32為B2、B3超寬巷整周模糊度；上標表示衛(wèi)星，下標表示參考站。其他符號與式(1)和式(2)含義相同。雙差M-W組合觀測值消除了雙差電離層延遲誤差一階項、對流層延遲誤差、站星間幾何距離(包含了軌道誤差)，不受衛(wèi)星鐘差、接收機鐘差的影響。從式(3)和式(4)可以看出，M-W組合模型求解B2、B3超寬巷整周模糊度僅受偽距和載波相位觀測值噪聲的影響，與參考站間距離無關(guān)，適合進行長距離參考站間載波相位整周模糊度的解算。

為了保證式(3)、式(4)的計算為連續(xù)弧段的觀測數(shù)據(jù)，利用M-W組合觀測值對BDS觀測數(shù)據(jù)進行粗差和周跳探測[18]。由于M-W組合觀測值主要受偽距觀測值的觀測噪聲影響，通過取平均值的方法削弱觀測值噪聲的影響，相應(yīng)的B2、B3雙差超寬巷整周模糊度和方差分別為

(5)

(6)

(7)

式中，prob0為B2、B3雙差超寬巷整周模糊度固定為最近整數(shù)的概率；b為B2、B3雙差超寬巷實數(shù)模糊度；bround為B2、B3雙差超寬巷實數(shù)模糊度的就近整數(shù)；σ為B2、B3雙差超寬巷整周模糊度的中誤差；erfc(·)為誤差函數(shù)。

利用嚴格的模糊度固定標準確定B2、B3雙差超寬巷整周模糊度，若雙差超寬巷模糊度滿足準則條件，則雙差超寬巷模糊度為正確值。B2、B3雙差超寬巷整周模糊度固定為整數(shù)的準則為：

(1) 實數(shù)雙差超寬巷模糊度與其最近整數(shù)之差的絕對值小于0.25周。

(2) 實數(shù)雙差超寬巷模糊度就近取整成功率大于0.99，保證固定雙差超寬巷模糊度的可靠性。

(3) 雙差超寬巷模糊度的方差小于等于0.005。

(4) 對于任意兩個以上的參考站，雙差超寬巷模糊度的代數(shù)和在理論上為零。一般BDS網(wǎng)絡(luò)RTK用于計算區(qū)域誤差的參考站數(shù)量大于3，以參考站A、B和C為例，則

(8)

1.2 參考站間寬巷整周模糊度解算

參考站間的B1、B2和B1、B3雙差寬巷組合載波相位觀測方程可表示為

(9)

(10)

(11)

(12)

式中，mf(·)為GMF投影函數(shù)；MF為星間投影函數(shù)之差；Twet為對流層濕延遲誤差；ZTDwet為天頂對流層濕延遲誤差；RZTDwetAB=ZTDwetA-ZTDwetB為參考站A、B的相對天頂對流層濕延遲誤差。式(12)右端第2個等式最后一項量級較小，可忽略，因此采用一個相對天頂對流層濕延遲參數(shù)估計所有可視衛(wèi)星的雙差對流層濕延遲誤差。

由于參考站坐標精確已知，假定歷元i，參考站A、B同步觀測到s+1顆衛(wèi)星，聯(lián)合式(9)—式(11)可得到B1、B2和B1、B3雙差寬巷組合載波相位觀測值的觀測方程

Lw(i)=Hw(i)Xw(i)

(13)

式中

式(13)在多歷元數(shù)據(jù)處理中將雙差電離層延遲誤差作為歷元參數(shù)，相對天頂對流層濕延遲誤差作為分段常數(shù)，每個歷元利用固定的寬巷模糊度計算得到的雙差電離層延遲誤差僅作為下一歷元的初值，采用參數(shù)消去法消掉法方程中與雙差電離層延遲誤差相關(guān)的信息，法方程中不進行電離層延遲誤差信息的疊加，只進行可被連續(xù)跟蹤時段內(nèi)所有衛(wèi)星的寬巷整周模糊度和相對天頂對流層濕延遲誤差的法方程疊加。在參數(shù)估計中，采用衛(wèi)星高度角定權(quán)法對雙差寬巷組合載波相位觀測值進行定權(quán)，并將正確固定的超寬巷整周模糊度與寬巷整周模糊度的線性關(guān)系作為強約束條件，各雙差衛(wèi)星對應(yīng)式(11)的權(quán)給予比衛(wèi)星高度角最高的寬巷載波相位觀測值的權(quán)稍大的權(quán)值。根據(jù)最小二乘計算原理，可估計雙差電離層延遲誤差、相對天頂對流層濕延遲誤差及雙差寬巷整周模糊度，之后搜索并確定雙差寬巷整周模糊度，并利用類似于式(8)的準則進行整周模糊度閉合條件檢驗。

1.3 參考站間B1、B2和B3整周模糊度解算

參考站間的雙差寬巷整周模糊度確定之后，其與B1、B2及B3雙差整周模糊度具有以下整數(shù)線性關(guān)系

(14)

將式(14)代入式(2)中得到B1、B2及B3雙差載波相位觀測方程

(15)

假定歷元i，參考站A、B同步觀測到s+1顆衛(wèi)星，由式(15)可得B1、B2和B3雙差載波相位觀測方程

L(i)=A(i)X(i)

(16)

其中

2 算例與分析

利用西南地區(qū)某省CORS網(wǎng)的BDS實測三頻數(shù)據(jù)進行算法檢驗，采用1 d的觀測數(shù)據(jù)，采樣間隔為1 s，衛(wèi)星截止高度角為10°。為了驗證所提出的BDS參考站三頻整周模糊度解算的可行性和有效性，利用精密數(shù)據(jù)處理軟件解算的雙差模糊度作為準確值進行驗證。該試驗數(shù)據(jù)共有3個參考站，參考站分布如圖1所示，BASE_A到BASE_B距離為70 km，BASE_B到BASE_C距離為113 km，BASE_C到BASE_A距離為78 km。

算例驗證中分別以衛(wèi)星C02-C09、C08-C09和C13-C09雙差衛(wèi)星為例，其中C09為觀測時段內(nèi)選擇的基準衛(wèi)星。3組雙差衛(wèi)星的首個連續(xù)觀測弧段對應(yīng)的各參考站間B2、B3雙差超寬巷整周模糊度計算結(jié)果如圖2所示。

圖1 測站分布Fig.1 Stations distribution

從圖2可以看出，衛(wèi)星C02-C09、C08-C09和C13-C09對應(yīng)的實數(shù)雙差超寬巷整周模糊度與其最近整數(shù)之差的絕對值遠遠小于0.25周， 3組雙

差衛(wèi)星在BASE_A-BASE_B、BASE_B-BASE_C、BASE_C-BASE_A參考站間的雙差超寬巷整周模糊度固定解分別為(-8，-5，13)、(-22，-7，29)和(-79，-79，0)，雙差超寬巷整周模糊度滿足閉合條件式(8)，保證了超寬巷整周模糊度的準確性。由于解算超寬巷整周模糊度的數(shù)學模型與站間距離無關(guān)，3條基線中BASE_B-BASE_C對應(yīng)的模糊度變化最為平滑，原因可能是其他兩條基線中的共同參考站BASE_A的接收機質(zhì)量和觀測環(huán)境稍差導(dǎo)致的。利用超寬巷波長較長的優(yōu)勢和雙差超寬巷整周模糊度固定為整數(shù)的準則可快速確定超寬巷整周模糊度。衛(wèi)星C02-C09、C08-C09和C13-C09在BASE_A、BASE_B和BASE_C組成的閉合網(wǎng)中模糊度收斂需要的平均歷元個數(shù)分別為5、8和3。

圖2 B2、B3雙差超寬巷整周模糊度Fig.2 Double differenced extra-wide-lane integer ambiguity of B2、B3

B2、B3雙差超寬巷整周模糊度準確固定之后，采用最小二乘方法估計雙差寬巷整周模糊度，并搜索確定雙差寬巷整周模糊度。圖3給出了各參考站間首個完整觀測弧段內(nèi)利用LAMBDA算法解算雙差寬巷整周模糊度的ratio值。由圖3可以看出，由于解算寬巷整周模糊度的數(shù)學模型綜合利用了所有衛(wèi)星的觀測值進行模糊度解算，加之寬巷整周模糊度之間的線性關(guān)系對觀測方程的強約束，雙差寬巷模糊度解算的可靠性很高，在雙差超寬巷整周模糊度準確固定的條件下，雙差寬巷整周模糊度可迅速準確固定。

圖3 雙差寬巷整周模糊度的ratio值Fig.3 Ratio of double differenced wide-lane integer ambiguity

進一步利用固定的寬巷整周模糊度構(gòu)造三頻整周模糊度之間的整數(shù)線性關(guān)系，將三頻整周模糊度之間的整數(shù)線性關(guān)系加入三頻雙差載波相位觀測方程中，綜合利用所有衛(wèi)星3個頻率的觀測值進行最小二乘參數(shù)估計。各參考站間首個完整觀測弧段內(nèi)衛(wèi)星C02-C09、C08-C09和C13-C09解算得到的B1雙差整周模糊度實數(shù)值如圖4所示。精密數(shù)據(jù)處理軟件解算的3組雙差衛(wèi)星在各參考站間的B1雙差整周模糊度準確值分別為(15，90，-105)、(218，92，-310)和(-9，86，-77)。由圖4可知，B1雙差整周模糊度實數(shù)值與模糊度準確值相差最大不超過1周，參數(shù)估計得到的模糊度實數(shù)值精度較高，利用LAMBDA算法可在兩個歷元內(nèi)準確的實現(xiàn)B1雙差整周模糊度固定。

圖4 B1雙差整周模糊度Fig.4 Double differenced integer ambiguity of B1

利用LAMBDA算法搜索并確定B1雙差整周模糊度之后，為了進一步驗證本文算法所解算的B1雙差整周模糊度的正確性和可靠性，可通過B1和B2雙差載波相位觀測方程線性組合計算得到B1頻率的雙差電離層延遲誤差和雙差對流層延遲誤差，雙差電離層延遲誤差和雙差對流層延遲誤差隨時間的變化情況可對所解算的B1雙差整周模糊度做出可靠性判斷，固定正確的雙差整周模糊度對應(yīng)的雙差電離層誤差和雙差對流層延遲誤差會在合理的范圍內(nèi)進行有規(guī)律的連續(xù)變化，相鄰歷元間相差較小，不會出現(xiàn)突然增大或減小、不連續(xù)及突變的情況，而且滿足閉合基線雙差電離層延遲誤差和雙差對流層延遲誤差之和等于零。

衛(wèi)星C02-C09高度角和電離層穿刺點的經(jīng)緯度差值變化如圖5所示，圖5中左側(cè)Y軸表示衛(wèi)星高度角，右側(cè)Y軸表示電離層穿刺點的經(jīng)緯度差值。由于GEO和IGSO衛(wèi)星均為高軌衛(wèi)星，衛(wèi)星運行角速度較慢，特別是GEO衛(wèi)星，其相對參考站的運行角速度幾乎為零，因此C02的高度角變化很小，而IGSO衛(wèi)星的運動周期約為24h，C09高度角變化較明顯。由于GEO衛(wèi)星電離層穿刺點變化較小，所以C09的穿刺點經(jīng)緯度變化決定衛(wèi)星C02-C09電離層穿刺點的經(jīng)緯度差值的變化。由于IGSO衛(wèi)星的運行軌道星下點軌跡以子午線為中心，所以C09的電離層穿刺點經(jīng)度變化較小，緯度變化較大。因此衛(wèi)星C02-C09電離層穿刺點的經(jīng)緯度差值中，經(jīng)度方向的差值變化較慢，緯度方向的差值變化稍大。衛(wèi)星C02-C09對應(yīng)的B1頻率的雙差電離層延遲誤差和雙差對流層延遲誤差如圖6和圖7所示。

圖5 C02-C09高度角和穿刺點差值Fig.5 The elevations and IPP difference of C02-C09

圖6 C02-C09電離層延遲誤差Fig.6 Ionosphere delay error of C02-C09

圖7 C02-C09對流層延遲誤差Fig.7 Tropospheric delay error of C02-C09

衛(wèi)星C08-C09高度角和電離層穿刺點的經(jīng)緯度差值變化如圖8所示。C08和C09均為IGSO衛(wèi)星，運動周期約為24h，其高度角變化較大，尤其是C08變化較明顯。從圖8可看出，衛(wèi)星C08-C09電離層穿刺點的經(jīng)緯度差值變化相對于衛(wèi)星C02-C09較快。衛(wèi)星C08-C09電離層穿刺點的經(jīng)緯度差中，經(jīng)度方向的差值變化較緩慢，而緯度方向的差值變化較快。衛(wèi)星C08-C09對應(yīng)的B1頻率的雙差電離層延遲誤差和雙差對流層延遲誤差如圖9和圖10所示。

衛(wèi)星C13-C09高度角和電離層穿刺點的經(jīng)緯度差值變化如圖11所示。C13為MEO衛(wèi)星，相比于GEO和IGSO衛(wèi)星，其衛(wèi)星運行角速度較快，運動周期約為12h，其高度角變化較大。衛(wèi)星C13-C09電離層穿刺點的經(jīng)緯度差值變化相對于衛(wèi)星C02-C09和C08-C09較快。衛(wèi)星C13-C09對應(yīng)的B1頻率的雙差電離層延遲誤差和雙差對流層延遲誤差如圖12和圖13所示。

圖8 C08-C09高度角和穿刺點差值Fig.8 The elevations and IPP difference of C08-C09

圖9 C08-C09電離層延遲誤差Fig.9 Ionosphere delay error of C08-C09

圖10 C08-C09對流層延遲誤差Fig.10 Tropospheric delay error of C08-C09

從衛(wèi)星C02-C09、C08-C09和C13-C09的雙差電離層延遲誤差和雙差對流層延遲誤差的變化情況可以看出，在觀測歷元的大部分時段內(nèi)二者變化比較平緩，雙差電離層延遲誤差隨著雙差衛(wèi)星電離層穿刺點的經(jīng)緯度方向的差值進行有規(guī)律的變化，雙差對流層延遲誤差隨著衛(wèi)星高度角進行有規(guī)律的變化，相鄰歷元間相差較小，不存在不連續(xù)及突變的情況，并且雙差大氣延遲誤差滿足閉合基線大氣延遲誤差之和等于零，即BASE_A-BASE_B、BASE_B-BASE_C和BASE_C-BASE_A上同一對雙差衛(wèi)星的大氣延遲誤差之和為零。3組雙差衛(wèi)星的雙差大氣延遲誤差在合理范圍內(nèi)進行有規(guī)律的連續(xù)變化進一步驗證了本文算法所解算的B1雙差整周模糊度的準確性和可靠性。

圖11 C13-C09高度角和穿刺點差值Fig.11 The elevations and IPP difference of C13-C09

圖12 C13-C09電離層延遲誤差Fig.12 Ionosphere delay error of C13-C09

圖13 C13-C09對流層延遲誤差Fig.13 Tropospheric delay error C13-C09

圖14給出了BASE_A-BASE_B、BASE_B-BASE_C和BASE_C-BASE_A參考站間24h觀測數(shù)據(jù)的B1、B2雙差寬巷整周模糊度和B1雙差整周模糊度的固定成功率，以精密數(shù)據(jù)處理軟件解算的雙差整周模糊度作準確值，使用LAMBDA算法搜索確定整周模糊度(Ratio>3)。圖中B1B2(triple)和B1(triple)表示本文方法的B1、B2雙差寬巷整周模糊度和B1雙差整周模糊度的固定成功率，B1B2(dual)和B1(dual)表示在雙頻觀測值條件下利用MW組合和無電離層組合解算的寬巷整周模糊度和B1整周模糊度的固定成功率。

圖14 模糊度固定成功率Fig.14 The success rate of ambiguity fixing

從圖14中可以看出，由于BDS偽距觀測值的噪聲較大，消耗了較多的觀測歷元用于確定寬巷模糊度，導(dǎo)致很多歷元模糊度不能固定，在寬巷整周模糊度固定之后，利用無電離層組合模型進行B1整周模糊度解算的成功率也偏低。相比于基于雙頻觀測值的模糊度解算方法，本文方法對應(yīng)的3條基線的B1、B2雙差寬巷整周模糊度具有較高的固定成功率，均為99.99%，在B1、B2雙差寬巷整周模糊度固定之后，進一步確定B1雙差整周模糊度，其固定成功率分別為98.83%、99.24%和98.81%，B1雙差整周模糊度部分不能固定的原因是受大氣延遲誤差影響較大的低高度角衛(wèi)星導(dǎo)致的。新方法中B1、B2雙差寬巷整周模糊度和B1雙差整周模糊度的固定速度較常規(guī)的雙頻模糊度解算方法有明顯提高，經(jīng)統(tǒng)計，新方法中波長較長的超寬巷整周模糊度的平均固定速度為9.8個歷元，B1、B2雙差寬巷整周模糊度的平均固定速度為2.5個歷元，B1雙差整周模糊度的平均固定速度為1.5個歷元，而常規(guī)的基于雙頻觀測值的模糊度解算方法中B1、B2雙差寬巷整周模糊度的平均固定速度為303.7個歷元，B1雙差整周模糊度的平均固定速度為598.4個歷元。從B1、B2雙差寬巷整周模糊度和B1雙差整周模糊度的固定成功率和固定速度可以看出，相比于GPS，BDS的信號質(zhì)量和接收設(shè)備不夠完善，常規(guī)方法將使BDS網(wǎng)絡(luò)RTK需要較長的初始化時間，影響B(tài)DS網(wǎng)絡(luò)RTK的實時性，新方法能提高BDS參考站整周模糊度的固定速度，縮短參考站網(wǎng)啟動和提供誤差改正信息所需的時間，能快速地為流動站用戶提供高精度的定位服務(wù)。本文方法在進行寬巷整周模糊度和B1整周模糊度確定過程中綜合利用了所有衛(wèi)星的觀測值信息，在超寬巷和寬巷整周模糊度固定之后，僅將其構(gòu)造的模糊度之間的整數(shù)線性關(guān)系作為約束條件，不直接利用超寬巷和寬巷載波相位觀測值，也不進行觀測值間的線性組合，避免了放大噪聲的不利影響，保證了參數(shù)估計的精度，進而可正確有效地固定參考站間三頻載波相位整周模糊度，縮短參考站網(wǎng)初始化時間。

3 結(jié) 論

本文研究了一種BDS網(wǎng)絡(luò)RTK參考站三頻整周模糊度解算方法。該方法利用多歷元BDS三頻觀測值，同時考慮大氣延遲誤差及超寬巷整周模糊度與寬巷整周模糊度的線性約束條件，并將三頻載波相位整周模糊度之間的整數(shù)線性關(guān)系加入原始頻率雙差載波相位整周模糊度參數(shù)估計觀測模型中，避免了噪聲較大的超寬巷及寬巷組合觀測值對參數(shù)估計的影響，有助于載波相位整周模糊度的快速、準確固定。通過實測CORS網(wǎng)的BDS三頻觀測數(shù)據(jù)進行試驗表明，該算法能夠?qū)崿F(xiàn)參考站間整周模糊度的準確固定，且可靠性高。本文方法是針對多歷元觀測值的BDS網(wǎng)絡(luò)RTK參考站三頻整周模糊度解算，對BDS網(wǎng)絡(luò)RTK參考站間單歷元三頻整周模糊度解算方法將作進一步研究。

[1] 楊元喜, 李金龍, 徐君毅, 等. 中國北斗衛(wèi)星導(dǎo)航系統(tǒng)對全球PNT用戶的貢獻[J]. 科學通報, 2011, 56(21): 1734-1740. YANG Yuanxi, LI Jinlong, XU Junyi, et al. Contribution of the Compass Satellite Navigation System to Global PNT Users[J]. Chinese Science Bulletin, 2011, 56(26): 2813-2819.

[2] 施闖, 趙齊樂, 李敏, 等. 北斗衛(wèi)星導(dǎo)航系統(tǒng)的精密定軌與定位研究[J]. 中國科學: 地球科學, 2012, 42(6): 854-861. SHI Chuang, ZHAO Qile, LI Min, et al. Precise Orbit Determination of BeiDou Navigation Satellites with Precise Positioning[J]. Science China Earth Sciences, 2012, 55(7): 1079-1086.

[3] 楊元喜. 北斗衛(wèi)星導(dǎo)航系統(tǒng)的進展、貢獻與挑戰(zhàn)[J]. 測繪學報, 2010, 39(1): 1-6. YANG Yuanxi. Progress, Contribution and Challenges of Compass/BeiDou Satellite Navigation System[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(1): 1-6.

[4] MONTENBRUCK O, HAUSCHILD A, STEIGENBERGER P, et al. Initial Assessment of the Compass/BeiDou-2 Regional Navigation Satellite System[J]. GPS Solutions, 2013, 17(2): 211-222.

[5] 高星偉, 陳銳志, 趙春梅. 網(wǎng)絡(luò)RTK算法研究與實驗[J]. 武漢大學學報(信息科學版), 2009, 34(11): 1350-1353. GAO Xingwei, CHEN Ruizhi, ZHAO Chunmei. A Network RTK Algorithm and Its Test[J]. Geomatics and Information Science of Wuhan University, 2009, 34(11): 1350-1353.

[6] 周樂韜, 黃丁發(fā), 袁林果, 等. 網(wǎng)絡(luò)RTK參考站間模糊度動態(tài)解算的卡爾曼濾波算法研究[J]. 測繪學報, 2007, 36(1): 37-42. ZHOU Letao, HUANG Dingfa, YUAN Linguo, et al. A Kalman Filtering Algorithm for Online Integer Ambiguity Resolution in Reference Station Network[J]. Acta Geodaetica et Cartographica Sinica, 2007, 36(1): 37-42.

[7] 高星偉, 劉經(jīng)南, 葛茂榮. 網(wǎng)絡(luò)RTK基準站間基線單歷元模糊度搜索方法[J]. 測繪學報, 2002, 31(4): 305-309. GAO Xingwei, LIU Jingnan, GE Maorong. An Ambiguity Searching Method for Network RTK Baselines between Base Stations at Single Epoch[J]. Acta Geodaetica et Cartographica Sinica, 2002, 31(4): 305-309.

[8] 祝會忠, 劉經(jīng)南, 唐衛(wèi)明, 等. 長距離網(wǎng)絡(luò)RTK基準站間整周模糊度單歷元確定方法[J]. 測繪學報, 2012, 41(3): 359-365. ZHU Huizhong, LIU Jingnan, TANG Weiming, et al. The Algorithm of Single-epoch Integer Ambiguity Resolution between Long-range Network RTK Base Stations[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(3): 359-365.

[9] 唐衛(wèi)明, 劉經(jīng)南, 施闖, 等. 三步法確定網(wǎng)絡(luò)RTK基準站雙差模糊度[J]. 武漢大學學報(信息科學報), 2007, 32(4): 305-308. TANG Weiming, LIU Jingnan, SHI Chuang, et al. Three Steps Method to Determine Double Difference Ambiguities Resolution of Network RTK Reference Station[J]. Geomatics and Information Science of Wuhan University, 2007, 32(4): 305-308.

[10] 李博峰, 沈云中, 周澤波. 中長基線三頻GNSS模糊度的快速算法[J]. 測繪學報, 2009, 38(4): 296-301. LI Bofeng, SHEN Yunzhong, ZHOU Zebo. A New Method for Medium and Long Range Three Frequency GNSS Rapid Ambiguity Resolution[J]. Acta Geodaetica et Cartographica Sinica, 2009, 38(4): 296-301.

[11] LI Bofeng, SHEN Yunzhong, FENG Yanming, et al. GNSS Ambiguity Resolution with Controllable Failure Rate for Long Baseline Network RTK[J]. Journal of Geodesy, 2014, 88(2): 99-112.

[12] FENG Yanming. GNSS Three Carrier Ambiguity Resolution Using Ionosphere-reduced Virtual Signals[J]. Journal of Geodesy, 2008, 82(12): 847-862.

[13] 高旺, 高成發(fā), 潘樹國, 等. 北斗三頻寬巷組合網(wǎng)絡(luò)RTK單歷元定位方法[J]. 測繪學報, 2015, 44(6): 641-648. GAO Wang, GAO Chengfa, PAN Shuguo, et al. Single-epoch Positioning Method in Network RTK with BDS Triple-frequency Widelane Combinations[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(6): 641-648.

[14] 祝會忠, 徐愛功, 高猛, 等. BDS網(wǎng)絡(luò)RTK中距離參考站整周模糊度單歷元解算方法[J]. 測繪學報, 2016, 45(1): 50-57. ZHU Huizhong, XU Aigong, GAO Meng, et al. The Algorithm of Single-epoch Integer Ambiguity Resolution between Middle-range BDS Network RTK Reference Stations[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(1): 50-57.

[15] 魏子卿, 葛茂榮. GPS相對定位的數(shù)學模型[M]. 北京: 測繪出版社, 1998. WEI Ziqing, GE Maorong. Mathematical Model of GPS Relative Positioning[M]. Beijing: Surveying and Mapping Press, 1998.

[16] TANG Weiming, DENG Chenlong, SHI Chuang, et al. Triple-frequency Carrier Ambiguity Resolution for BeiDou Navigation Satellite System[J]. GPS Solutions, 2014, 18(3): 335-344.

[17] CHU Fengyu, YANG Ming, WU J. A New Approach to Modernized GPS Phase-only Ambiguity Resolution over Long Baselines[J]. Journal of Geodesy, 2016, 90(3): 241-254.

[18] 黃令勇, 翟國君, 歐陽永忠, 等. 削弱電離層影響的三頻TurboEdit周跳處理方法[J]. 測繪學報, 2015, 44(8): 840-847. DOI: 10.11947/j.AGCS.2015.20140380. HUANG Lingyong, ZHAI Guojun, OUYANG Yongzhong, et al. Triple-frequency Turboedit Cycle-slip Processing Method Weakening Ionospheric Activity[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(8): 840-847. DOI: 10.11947/j.AGCS.2015.20140380.

[19] GE M, GENDT G, ROTHACHER M, et al. Resolution of GPS Carrier-phase Ambiguities in Precise Point Positioning (PPP) with Daily Observations[J]. Journal of Geodesy, 2008, 82(7): 389-399.

[20] DONG Danan, BOCK Y. Global of GPS Positioning System Network Analysis with Phase Ambiguity Resolution Applied to Crustal Deformation Studies in California[J]. Journal of Geophysical Research: Solid Earth, 1989, 94(B4): 3949-3966.

[21] 潘宗鵬, 柴洪洲, 劉軍, 等. 基于部分整周模糊度固定的非差GPS精密單點定位[J]. 測繪學報, 2015, 44(11): 1210-1218. PAN Zongpeng, CHAI Hongzhou, LIU Jun, et al. GPS Partial Ambiguity Resolution Method for Zero-difference Precise Point Positioning[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(11): 1210-1218.

[22] 葛茂榮, 劉經(jīng)南. GPS定位中對流層折射估計研究[J]. 測繪學報, 1996, 25(4): 285-291. GE Maorong, LIU Jingnan. The Estimation Methods for Tropospheric Delays in Global Positioning System[J]. Acta Geodaetica et Cartographica Sinica, 1996, 25(4): 285-291.

(責任編輯：陳品馨)

The Algorithm of Triple-frequency Ambiguity Resolution between BDS Network RTK Reference Stations

GAO Meng1,XU Aigong1,ZHU Huizhong1,GE Maorong1,2,YANG Qiushi1

1. School of Geomatics，Liaoning Technical University，F(xiàn)uxin 123000，China；2. German Research Center for Geosciences，Potsdam 14473，Germany

BeiDou navigation satellite system is currently the only system which the whole constellation offers triple-frequency observation data, triple-frequency observation is conducive to carrier phase ambiguity fixing quickly and accurately. The algorithm of triple-frequency integer ambiguity resolution between BDS network RTK reference stations is proposed. The observation of B2 frequency and B3 frequency and the strict criterion of ambiguity fixing are used to determine extra-wide-lane integer ambiguity. The linear relationship between fixed extra-wide-lane integer ambiguity and other wide-lane integer ambiguity is selected as constraint conditions, then the wide-lane integer ambiguity, relative zenith tropospheric delay error and ionospheric delay error are estimated and the wide-lane integer ambiguity can be fixed by search mode. There is a kind of integer linear relationship between fixed wide-lane integer ambiguity and triple-frequency carrier phase integer ambiguity, this linear relationship is added to the observation model of carrier phase ambiguity parameter estimation and the integer ambiguity can be fixed. This algorithm is tested by the measured triple-frequency data of BDS CORS network. The results indicate that triple-frequency carrier phase ambiguity can be fixed accurately and effectively with the method.

BeiDou navigation satellite system；network RTK；triple-frequency；integer ambiguity；relative zenith tropospheric delay error；ionospheric delay error

The National Key Research and Development Program (No.2016YFC0803102)；The National Natural Science Foundation of China (Nos.41504010；41474020；41504030)；The National High-tech Research and Development Program of China (863 Program) (No.2014AA123101)；The Project for Liaoning High School Innovation Team (No.LT2015013)；The Program for Liaoning Doctoral Foundation (Nos.20141141；201501126)

XU Aigong

高猛，徐愛功，祝會忠，等.BDS網(wǎng)絡(luò)RTK參考站三頻整周模糊度解算方法[J].測繪學報，2017,46(4)：442-452.

10.11947/j.AGCS.2017.20160179. GAO Meng, XU Aigong, ZHU Huizhong, et al.The Algorithm of Triple-frequency Ambiguity Resolution between BDS Network RTK Reference Stations[J]. Acta Geodaetica et Cartographica Sinica,2017,46(4):442-452. DOI:10.11947/j.AGCS.2017.20160179.

P288

A

1001-1595(2017)04-0442-11

國家重點研發(fā)計劃(2016YFC0803102);國家自然科學基金(41504010；41474020；41504030);國家863計劃(2014AA123101);遼寧省高等學校創(chuàng)新團隊項目(LT2015013);遼寧省博士啟動基金(20141141；201501126)

2016-04-22

高猛(1988—)，男，博士生，研究方向為高精度GNSS定位算法與研究。First author： GAO Meng(1988—)，male，PhD candidate，majors in the algorithm and application of GNSS for high-accuracy dynamic positioning.

E-mail： gaomeng512@163.com

徐愛功

E-mail： xu_ag@126.com

修回日期： 2017-02-01