李 慧， 趙 琳， 李 亮， 丁繼成

(哈爾濱工程大學(xué) 自動(dòng)化學(xué)院， 哈爾濱 150001)

一種雙頻非組合實(shí)時(shí)精密定位方法

李 慧， 趙 琳， 李 亮， 丁繼成

(哈爾濱工程大學(xué) 自動(dòng)化學(xué)院， 哈爾濱 150001)

為滿足單系統(tǒng)單基線雙頻數(shù)據(jù)條件下的實(shí)時(shí)精密定位需求，提出一種雙頻非組合實(shí)時(shí)精密定位技術(shù)，基于站間-星間載波相位及偽距觀測量雙差觀測模型，實(shí)現(xiàn)單系統(tǒng)單基線雙頻非組合RTK(Real Time Kinematic). 通過分析雙差模型觀測量冗余度，確立模型殘余誤差處理策略，設(shè)定狀態(tài)向量，推導(dǎo)并建立狀態(tài)預(yù)測方程及測量方程，實(shí)時(shí)更新狀態(tài)向量變換矩陣，根據(jù)隨機(jī)模型調(diào)整兩種觀測量數(shù)據(jù)的權(quán)重，最后利用擴(kuò)展卡爾曼濾波器技術(shù)得到實(shí)時(shí)定位結(jié)果. 文中基于幾組中長基線實(shí)驗(yàn)，通過考察定位結(jié)果的三維定位誤差及整周模糊度成功固定率，驗(yàn)證該方法的有效性. 實(shí)驗(yàn)結(jié)果表明，在中長基線條件下進(jìn)行實(shí)時(shí)定位，該方法精度可以達(dá)到厘米級，基線長度為135.6 km時(shí)，整周模糊度固定成功率為97.3%，東向、北向、天向的CEP95定位誤差分別為1.35 cm、1.84 cm、7.08 cm. 雙頻非組合技術(shù)充分利用差分技術(shù)的優(yōu)勢消除與距離無關(guān)的相關(guān)誤差，并有效地避免了觀測值組合過程所引起的觀測噪聲，可以實(shí)現(xiàn)中長基線條件下的厘米級實(shí)時(shí)定位.

實(shí)時(shí)精密定位；雙頻；非組合；中長基線；擴(kuò)展卡爾曼

載波相位觀測量因其高精度被廣泛應(yīng)用在全球衛(wèi)星導(dǎo)航系統(tǒng)(Global Navigation Satellite System, GNSS)中，但只有當(dāng)準(zhǔn)確解算出整周模糊度，利用載波相位觀測量進(jìn)行定位才變得有意義. 差分定位原理主要依靠載波相位和偽距雙差可以消除相關(guān)誤差，但當(dāng)接收機(jī)間距離超過10 km時(shí)，該誤差處理策略失效[1]. 隨距離增加的大氣層延遲誤差殘差及缺乏足夠的數(shù)據(jù)冗余度去解決整周模糊度是制約其基線長度的一個(gè)重要因素[2]. 對于中長基線，電離層與對流層在經(jīng)過站間、星間雙差后的殘差都應(yīng)該作為一未知參數(shù)被考慮[3]，但這會(huì)加重整周模糊度求解模型的病態(tài)性，造成模糊度固定失敗. 卡爾加里大學(xué)利用幾種觀測量線性組合包括雙差寬巷組合和電離層無關(guān)線性組合及卡爾加里大學(xué)提出的模型，消除電離層和對流層誤差對中長基線下整周模糊度求解的影響[4]；部分學(xué)者利用三頻觀測數(shù)據(jù)，通過電離層無關(guān)組合消除電離層延遲誤差對整周模糊度求解的影響[5-7]；除了上述方法還可以利用Thai電離層地圖解決RTK中模型中電離層誤差問題[8]；對于對流層誤差，最常用的方法是通過經(jīng)驗(yàn)?zāi)Ｐ颓蠼猓チ_茨瓦夫大學(xué)針對波蘭地區(qū)研究了兩種對流層延遲誤差模型，并對其在單基線和多基線場景下的精度進(jìn)行了研究分析[9]；對于快速精密定位常見的有兩種對流層延遲估計(jì)模型[10]；近些年來，隨著CORS系統(tǒng)等地面基站的建設(shè)，多基站輔助解決大氣誤差策略應(yīng)運(yùn)而生，利用網(wǎng)絡(luò)基站通過插值得到網(wǎng)絡(luò)電離層延遲修正量[11]，但插值過程的誤差是不可避免的，因此發(fā)展了一系列方法，如對天頂對流層參數(shù)限制提出一種正則化方法[12]，利用經(jīng)驗(yàn)對流層模型GT2生成一個(gè)高精度的有先驗(yàn)參數(shù)的對流層延遲，改善傳統(tǒng)網(wǎng)絡(luò)RTK模型中整周模糊度求解方程的病態(tài)性；利用2個(gè)基站對用戶進(jìn)行相對精密定位[13]，其中一個(gè)基站離用戶基站非常近，可以看做零基線，其中零基線的基站跟用戶站有著相似的電離層及對流層延遲. 將已知基線長度、不同基線之間整周模糊度的關(guān)系以及相似的電離層對流層延遲等已知條件充分應(yīng)用到相對精密定位中，改善整周模糊度的固定情況以及最終的定位精度；科廷大學(xué)提出了利用多系統(tǒng)提高中長基線整周模糊度固定率的方法[14-15]，更多利用多系統(tǒng)組合數(shù)據(jù)提高精密定位精度的研究見文獻(xiàn)[16-18].

為滿足單一系統(tǒng)，沒有足夠的觀測網(wǎng)絡(luò)數(shù)據(jù)情況時(shí)的實(shí)時(shí)定位需求，本文提出了單系統(tǒng)、單基線條件下的一種雙頻非組合實(shí)時(shí)精密定位技術(shù)，采用雙頻雙差模型，消除衛(wèi)星鐘差和接收機(jī)鐘差等與距離無關(guān)的相關(guān)誤差的影響，同時(shí)減弱諸如軌道誤差、大氣折射誤差等系統(tǒng)性誤差的影響，對于對流層延遲誤差，本文利用Saastamoinen模型估計(jì)后對原始觀測量進(jìn)行補(bǔ)償，對于影響整周模糊度解算的主要誤差電離層雙差誤差則設(shè)為未知參數(shù)，最后通過擴(kuò)展卡爾曼濾波器進(jìn)行用戶位置求解. 文中單獨(dú)處理兩頻率上的觀測量數(shù)據(jù)，有效避免了觀測值組合過程所引起的觀測噪聲.

1 數(shù)學(xué)模型

1.1雙頻雙差基礎(chǔ)模型

(1)

1.2誤差處理策略

高精度RTK定位主要依據(jù)載波相位觀測量，為提高整周模糊度解算的準(zhǔn)確性，必須設(shè)法解決式(1)中影響定位結(jié)果的各個(gè)誤差. 對于短基線而言，在經(jīng)過站間及星間雙差后，認(rèn)為消除或近似消除了衛(wèi)星鐘差、接收機(jī)鐘差、軌道誤差、對流層延遲誤差、電離層延遲誤差等，但如果用戶和基站間基線長度超過10 km，有些誤差經(jīng)過雙差后的殘差需要分析、考慮，尤其是與距離相關(guān)的一些誤差，如電離層誤差、對流層誤差及衛(wèi)星軌道誤差，本文的誤差處理策略主要考慮了這三種誤差的處理方法.

1.2.1 衛(wèi)星軌道誤差

衛(wèi)星軌道誤差對最終定位結(jié)果的影響隨基線長度的增加而增加見下式.

.

(2)

式中: dρ為衛(wèi)星軌道誤差，ρ為衛(wèi)星與接收機(jī)間距離，db為因衛(wèi)星軌道誤差對用戶定位結(jié)果造成的誤差，b為基線長度. 衛(wèi)星到接收機(jī)間的距離大約為20 000 km，而中長基線的基線長度從10～300 km間不等. 對于衛(wèi)星軌道誤差，有些單位提供了不同精度的軌道產(chǎn)品，對于IGS產(chǎn)品的精度如表1[19].

表1 IGS產(chǎn)品精度表

1.2.2 對流層延遲誤差

對流層天頂延遲總量由90%的干延遲分量和10%的濕延遲分量構(gòu)成[20]，對流層延遲誤差通常表示為天頂方向的對流層折射量與仰角的映射函數(shù)之積，如式

T(e)=mh(e)×zh+mw(e)×zw.

(3)

式中:e為天頂角，zh,zw分別為對流層的干濕分量，mh,mw分別為干濕分量的仰角映射函數(shù). 關(guān)于對流層天頂延遲量相關(guān)模型，見參考文獻(xiàn)[21-23].

對于基線長度超過50 km的用戶，應(yīng)將雙差對流層殘余作為未知參數(shù)考慮，但這會(huì)加重整周模糊度求解模型的病態(tài)性[12]. 因此對于對流層延遲本文采用模型估計(jì)方法，盡量降低對流層延遲的誤差干擾. 對流層天頂延遲量采用Saastamoinen模型估計(jì)[24]，得到對流層天頂延遲的干延遲分量和由水蒸氣引起的濕延遲分量，映射函數(shù)采用GMF[25]. 對于對流層延遲誤差，本文直接在原始觀測量上進(jìn)行補(bǔ)償，從而降低定位模型中的未知參數(shù)個(gè)數(shù)，增加模型方程觀測量冗余度.

1.2.3 電離層延遲誤差

電離層延遲誤差是影響用戶定位的主要誤差源之一，不同測站的電離層延遲誤差不同，同一測站不同觀測方向上的電離層延遲誤差也不同. 一般對于短基線用戶而言，其電離層誤差與基站誤差近似認(rèn)為相等，因此可以通過差分技術(shù)消除，但當(dāng)用戶基線超過10 km后，電離層差分殘余誤差將對定位結(jié)果造成較大誤差.

為了消除中長基線的電離層延遲誤差，許多學(xué)者提出了利用電離層無關(guān)組合消除電離層延遲誤差，但雙頻電離層無關(guān)組合會(huì)使得測量噪聲增加，由誤差傳播定律，偽距觀測量的無電離層組合觀測值的觀測噪聲為單頻率上雙差偽距噪聲的2.978倍[26]，并且因?yàn)榻M合觀測量的波長太短，造成短時(shí)間內(nèi)模糊度固定困難[27]，因此該方法大部分基于網(wǎng)絡(luò)基站或三頻數(shù)據(jù)條件下[4-7]. 針對單基線單系統(tǒng)雙頻數(shù)據(jù)條件下中長基線的實(shí)時(shí)精密定位，本文將電離層延遲站間-星間雙差殘差值設(shè)為未知參數(shù)，假設(shè)用戶和接收機(jī)在頻率fn(n=1,2)上同時(shí)觀測到m顆衛(wèi)星，則選定參考星后，對電離層延遲誤差做站間-星間差分處理，得到m-1個(gè)雙差電離層延遲量作為擴(kuò)展卡爾曼的未知狀態(tài)向量求解.

1.2.4 多徑誤差

由公式(1)，除了衛(wèi)星軌道誤差、對流層延遲誤差、電離層延遲誤差，用戶在進(jìn)行了站間-星間雙差后，還存在多徑誤差. 本文中對多徑誤差不作重點(diǎn)研究，假設(shè)接收機(jī)處于空曠區(qū)域，因此對于文中所提的雙頻非組合實(shí)時(shí)精密定位方法，從建模到利用實(shí)際數(shù)據(jù)進(jìn)行驗(yàn)證，均沒有考慮多徑誤差.

1.3單基線雙頻非組合模型

1.3.1 EKF建模及求解

EKF是通過線性化處理來實(shí)現(xiàn)非線性濾波估計(jì)，是用高斯分布來逼近系統(tǒng)狀態(tài)的后驗(yàn)概率密度. 相對于處理非線性濾波的粒子波波器等，在計(jì)算速度上，EKF具有明顯的優(yōu)勢. 本文采用EKF對式(1)中的雙頻雙差模型進(jìn)行求解，首先利用泰勒展開式對非線性函數(shù)進(jìn)行線性化處理，然后再在標(biāo)準(zhǔn)卡爾曼濾波框架下進(jìn)行遞歸濾波，包括測量更新和時(shí)間更新兩個(gè)過程. 線性化后的雙差非組合定位的狀態(tài)預(yù)測方程和測量方程如式：

Xk=Fk,k-1Xk-1+Wk,

(4)

yk=HkXk+ek.

(5)

式中:Xk為k時(shí)刻狀態(tài)向量，yk為k時(shí)刻觀測值，F(xiàn)k,k-1為第k-1時(shí)刻到k時(shí)刻的狀態(tài)轉(zhuǎn)移矩陣，Wk為狀態(tài)模型噪聲，ek為測量噪聲，假設(shè)狀態(tài)模型噪聲Wk和測量噪聲ek互為不相關(guān)的白噪聲序列，其協(xié)方差矩陣分別為Qw和Re. 假設(shè)1號(hào)衛(wèi)星為選定的參考衛(wèi)星，則未知狀態(tài)向量Xk為

(6)

yk=(T.

(7)

其中：

(8)

測量模型向量b(x)、觀測矩陣Hg為：

(9)

(10)

(11)

其中:Ss為星間單差矩陣，Gg為用戶-衛(wèi)星方向余弦向量，E0為單位矩陣，具體如式(12)到式(15).

(12)

(13)

(14)

(15)

測量誤差協(xié)方差矩陣Re描述測量過程中的誤差，具體如式

(16)

其中:D為站間-星間雙差矩陣，RΦ,n為頻率n上的非差相位觀測值的協(xié)方差矩陣，RP,n為頻率n上的非差偽距觀測值的協(xié)方差矩陣.

利用第k-1個(gè)歷元的得到的狀態(tài)向量對應(yīng)的方差協(xié)方差矩陣P及過程噪聲協(xié)方差矩陣計(jì)算第k個(gè)歷元的狀態(tài)向量對應(yīng)的協(xié)方差矩陣，計(jì)算求得卡爾曼濾波器增益K，進(jìn)而得到狀態(tài)向量在第k個(gè)歷元的的預(yù)測結(jié)果，如式：

Pk=(E-KkHg,k)Pk,k-1

(17)

式(17)為標(biāo)準(zhǔn)卡爾曼濾波器的測量更新過程. 再按式(18)進(jìn)行時(shí)間更新過程，為第k-1個(gè)歷元求解做準(zhǔn)備.

(18)

其中:

(19)

(20)

式中：τ為相鄰兩歷元間的時(shí)間間隔，qr、qν、qI、qN分別為位置、速度、電離層延遲及整周模糊度的過程噪聲協(xié)方差矩陣，共同構(gòu)成矩陣過程噪聲方差協(xié)方差矩陣Qw. 為了避免因過度依賴模型而造成的濾波結(jié)果發(fā)散現(xiàn)象，利用Qw矩陣調(diào)節(jié)狀態(tài)協(xié)方差矩陣P，避免P變得太小或者為0而導(dǎo)致的濾波器失效. 根據(jù)卡爾曼濾波原理，最終的濾波結(jié)果是由模型狀態(tài)預(yù)測值及觀測值進(jìn)行加權(quán)更新得到的，如果觀測精度較高，則增加其權(quán)重，使得濾波結(jié)果傾向于觀測結(jié)果，如果觀測精度較低，則增加狀態(tài)預(yù)測值對濾波結(jié)果的貢獻(xiàn).

在利用EKF得到L1及L2兩個(gè)頻率上的雙差整周模糊度的浮點(diǎn)解及其方差協(xié)方差矩陣后，利用LAMBDA算法進(jìn)行雙差整周模糊度的固定，得到用戶最后定位結(jié)果，關(guān)于LAMBDA算法具體見文獻(xiàn)[28-29].

1.3.2 狀態(tài)向量變換矩陣

1)參考衛(wèi)星不變，其它共視衛(wèi)星增加或減少. 如果當(dāng)前歷元丟失非參考星，則將該衛(wèi)星對應(yīng)的整周模糊度項(xiàng)及電離層延遲項(xiàng)去掉，同時(shí)需去掉對應(yīng)狀態(tài)向量方差協(xié)方差矩陣P的相關(guān)行和列；如果當(dāng)前歷元增加非參考星，則構(gòu)建該非參考星的站間-星間雙差載波相位觀測量，利用原有的已知雙差模糊度向量及相應(yīng)衛(wèi)星的觀測量，求出該時(shí)刻的基線向量，進(jìn)而求得新增衛(wèi)星的雙差模糊度向量，如公式(21)，或者直接利用其偽距和載波相位觀測量，初始化其雙差模糊度向量，如公式(22)，對于雙差電離層延遲項(xiàng)，則直接賦予其0的初值，同時(shí)將方差協(xié)方差矩陣P相關(guān)行和列元素置為0，如式(23).

Nnew=Hnewbur-yΦ,new,

(21)

(22)

(23)

2)參考星改變. 新歷元數(shù)據(jù)在解算前，需要檢測參考衛(wèi)星較上個(gè)歷元是否發(fā)生變換，若依據(jù)高度角最大原則，假設(shè)第k+1個(gè)歷元參考星發(fā)生了變化，則需要變換矩陣Tk,N將第k個(gè)歷元中以參考星p得到的雙差模糊度解向量映射至新參考衛(wèi)星pnew，具體如公式(24)，對于雙差電離層延遲項(xiàng)，利用變換矩陣Tk,I將第k個(gè)歷元中以參考星p得到的雙差電離層延遲量映射至新參考衛(wèi)星pnew，如式(25)，其中，Tk,N和Tk,I為單位陣對應(yīng)新參考衛(wèi)星的列統(tǒng)一減去1得到，由對應(yīng)狀態(tài)向量中的位置、速度形成的單位陣，結(jié)合Tk,N和Tk,I共同構(gòu)成從第k個(gè)歷元到第k+1個(gè)歷元的狀態(tài)解向量的變換矩陣Tk, 根據(jù)Tk，利用對第k+1個(gè)歷元對應(yīng)的狀態(tài)向量方差協(xié)方差矩陣Pk+1,N,new進(jìn)行更新，具體如式子(26).

Nk+1,new=Tk,NNk,

(24)

(25)

(26)

2 實(shí)驗(yàn)分析

為驗(yàn)證本文提出的單基線雙頻非組合實(shí)時(shí)精密定位方法的性能，利用4組不同基線長度的實(shí)際數(shù)據(jù)對其進(jìn)行了考察. 本次實(shí)驗(yàn)所采用的數(shù)據(jù)為2016年5月1日UNAVCO-PBO機(jī)構(gòu)位于美國加利福尼亞州圣阿爾多、圣米格爾、科林加、Huasna及德拉諾的基站采集的數(shù)據(jù)，接收機(jī)采用TRIMBLE NETRS，天線類型為TRM29659.00，連續(xù)采集24 h數(shù)據(jù)，采樣間隔為30 s，共2 880個(gè)歷元，各組基線的具體情況如表2所示.

表2 實(shí)驗(yàn)數(shù)據(jù)

本次實(shí)驗(yàn)中的4組基線的基站接收機(jī)均位于空曠地帶，可以忽略多徑誤差，對流層延遲誤差已通過建模對原始觀測量進(jìn)行了補(bǔ)償，因此，雙差雙頻非組合定位模型中的主要誤差為電離層延遲誤差. 4組實(shí)驗(yàn)中，在利用擴(kuò)展卡爾曼濾波器對進(jìn)行狀態(tài)參數(shù)估計(jì)時(shí)，狀態(tài)向量中的位置參數(shù)初始值取單點(diǎn)定位得到的結(jié)果，對應(yīng)的初始標(biāo)準(zhǔn)差設(shè)置為10 m，電離層延遲參數(shù)初始值設(shè)為0，因文中處理的最長基線為135.6 km，對應(yīng)的初始標(biāo)準(zhǔn)差值統(tǒng)一設(shè)置為0.1 m，當(dāng)基線長度更長時(shí)，需要適當(dāng)調(diào)整初始標(biāo)準(zhǔn)差值，否則會(huì)影響模糊度快速固定，雙差整周模糊度初始值可由雙差載波相位觀測值和雙差偽距觀測值得到，對應(yīng)的標(biāo)準(zhǔn)差設(shè)為20周. 文中同時(shí)利用了載波相位觀測量和偽距觀測量，這兩種觀測量的精度差別較大，因此采用隨機(jī)加權(quán)模型進(jìn)行調(diào)整，其中載波相位觀測量對應(yīng)的觀測噪聲標(biāo)準(zhǔn)差設(shè)為0.009 m，偽距觀測量對應(yīng)的觀測噪聲標(biāo)準(zhǔn)差設(shè)為0.9 m，兩觀測量內(nèi)部的不同衛(wèi)星-基站雙差觀測量采用等權(quán)模式.

圖1～4為實(shí)驗(yàn)A、B、C、D根據(jù)本文所建立的實(shí)時(shí)定位模型估計(jì)出的雙差電離層延遲誤差，其中，圖中橫坐標(biāo)為觀測時(shí)間，縱坐標(biāo)為L1載波對應(yīng)的雙差電離層延遲誤差.

圖1 A組數(shù)據(jù)估計(jì)出的L1載波對應(yīng)的雙差電離層誤差

Fig.1 Residual double difference ionospheric delay on L1 of data A

圖2 B組數(shù)據(jù)估計(jì)出的L1載波對應(yīng)的雙差電離層誤差

Fig.2 Residual double difference ionospheric delay on L1 of data B

圖3 C組數(shù)據(jù)估計(jì)出的L1載波對應(yīng)的雙差電離層誤差

Fig.3 Residual double difference ionospheric delay on L1 of data C

圖4 D組數(shù)據(jù)估計(jì)出的L1載波對應(yīng)的雙差電離層誤差

Fig.4 Residual double difference ionospheric delay on L1 of data D

不同的顏色代表不同共視衛(wèi)星對的雙差電離層延遲，橫坐標(biāo)為觀測時(shí)間，以小時(shí)為單位，對應(yīng)UTC(GMT)時(shí)間0點(diǎn)到24點(diǎn)，即本地時(shí)間16點(diǎn)到第二天16點(diǎn)，圖中后半段對應(yīng)白天中午的電離層延遲，因此波動(dòng)較大. 從圖中可以看出，當(dāng)基線長度僅為20.7 km時(shí)，雙差后的電離層延遲殘差高的能達(dá)到10 cm，這對于厘米級的精密定位是不能容忍的，并且隨著基線長度的增加，雙差電離層延遲殘差越來越大，當(dāng)基線長度為135.6 km時(shí)，達(dá)到50 cm.

圖5～8為實(shí)驗(yàn)A、B、C、D的三維定位誤差，其中橫坐標(biāo)為橫坐標(biāo)為觀測時(shí)間，以小時(shí)為單位，對應(yīng)UTC(GMT)時(shí)間0點(diǎn)到24點(diǎn)，即本地時(shí)間16點(diǎn)到第二天16點(diǎn)，縱坐標(biāo)為定位誤差，以米為單位，圖中不同的顏色及線型分別代表東向、北向及天向誤差.

表3給出了本次實(shí)驗(yàn)中4組基線數(shù)據(jù)的詳細(xì)定位結(jié)果. 當(dāng)用戶與基站間的基線長度為20.7 km時(shí)，利用本文提出的實(shí)時(shí)精密定位方法，整周模糊度固定率可達(dá)98.8%，每次整周模糊度固定時(shí)間開銷約為0.015 s，在利用CEP95準(zhǔn)則統(tǒng)計(jì)下的東北天定位誤差均為厘米級，RMS誤差東向和北向達(dá)毫米級，天向?yàn)槔迕准?；隨著用戶與基站間基線長度增加，整周模糊度固定率逐漸下降，當(dāng)基線長度為135.6 km時(shí)，整周模糊度固定率為97.3%，此時(shí)，CEP95準(zhǔn)則下東北天三個(gè)方向上的誤差仍為厘米級，RMS誤差為厘米級.

圖5 A組數(shù)據(jù)三維定位誤差

圖6 B組數(shù)據(jù)三維定位誤差

圖7 C組數(shù)據(jù)三維定位誤差

圖8 D組數(shù)據(jù)三維定位誤差

實(shí)驗(yàn)基線∕km固定率∕%CEP95定位誤差∕cm東向北向天向RMS誤差∕cm東向北向天向A20.798.81.291.503.620.690.712.13B49.097.92.222.014.691.262.042.35C105.897.71.172.134.332.300.892.21D135.697.31.351.847.082.001.354.00

3 結(jié) 論

本文在單系統(tǒng)、單基線及雙頻數(shù)據(jù)條件下，提出了一種適用于中長基線的雙頻非組合實(shí)時(shí)精密定位方法. 該方法以載波相位、偽距的站間-星間雙差為基礎(chǔ)模型，通過分析中長基線條件下主要雙差殘余誤差量，提出相應(yīng)誤差處理策略，建立精密定位解算模型. 為了避免引起觀測矩陣秩虧，將主要誤差源雙差電離層延遲作為未知量，而對對流層延遲則利用Saastamoinen模型在原始觀測量上進(jìn)行補(bǔ)償. 該方法充分利用了載波相位和偽距雙頻觀測量數(shù)據(jù)，通過隨機(jī)模型調(diào)整兩種觀測量數(shù)據(jù)的權(quán)重，以降低低精度偽距觀測量的應(yīng)用帶來的定位誤差，最后利用擴(kuò)展卡爾曼濾波器實(shí)現(xiàn)定位，并通過實(shí)時(shí)建立當(dāng)前歷元與上個(gè)歷元的狀態(tài)變換矩陣，保證擴(kuò)展卡爾曼濾波的連續(xù)性和傳遞的準(zhǔn)確性. 通過4組基線實(shí)驗(yàn)驗(yàn)證了文中基于擴(kuò)展卡爾曼濾波的雙頻非組合實(shí)時(shí)精密定位方法的有效性，基線長度達(dá)135.6 km時(shí)，該方法仍能實(shí)現(xiàn)厘米級的實(shí)時(shí)定位. 該方法單獨(dú)處理兩頻率上的觀測量，有效避免了觀測值組合過程所引起的觀測噪聲，是中長基線條件下單基線實(shí)時(shí)精密定位的另一種思路.

[1] PAZIEWSKI J, WIELGOSZ P. Assessment of GPS+ Galileo and multi-frequency Galileo single-epoch precise positioning with network corrections[J]. GPS Solutions, 2014, 18(4):571-579. DOI:10.1007/s10291-013-0355-3.

[2] KASHANI I, GREJNER D, WIELGOSZ P. Toward instantaneous network based real time kinematic GPS over 100 km distance[J]. Navigation, 2005, 52(4):239-245. DOI:10.1002/j.2161-4296.2005.tb00366.x.

[3] WIELGOSZ P, KASHANI I, GREJNER D. Analysis of long-range network RTK during a severe ionospheric storm[J], Journal of Geodesy, 2005, 79(9):524-531.DOI:10.1007/s00190-005-0003-y.

[4] WANG D, GAO C, PAN S. Single-epoch integer ambiguity resolution for long-baseline RTK with ionosphere and troposphere estimation[C]// SUN J, JIAO W, WU H, et al. Lecture Notes in Electrical Engineering. Berlin Heidelberg: Springer, 2013: 125-137. DOI:10.1007/978-3-642-37398-5_12.

[5] ODIJK D, TEUNISSEN P, TIBERIUS C. Triple-frequency ionosphere-free phase combinations for ambiguity resolution[C]// ASKILDT M. Proc. of the ENC-GNSS 2002. Copenhagen, Denmark: The European Navigation Conference press, 2002: 27-30.

[6] WANG K, ROTHACHER M. Ambiguity resolution for triple-frequency geometry-free and ionosphere-free combination tested with real data[J]. Journal of Geodesy, 2013, 87(6): 539-553. DOI:10.1007/s00190-013-0630-7.

[7] ODIJK D. Ionosphere-free phase combinations for modernized GPS[J]. Journal of surveying engineering, 2003, 129(4): 165-173. DOI:10.1061/(ASCE)0733-9453(2003)129:4(165).

[8] CHAROENKALUNYUTA T, SATIRAPOD C. Effect of Thai ionospheric Maps (THIM) model on the performance of network based RTK GPS in Thailand[J]. Survey Review, 2014, 46(334): 1-6. DOI: 10.1179/1752270613Y.0000000055.

[9] HADAS T, KAPLON J, BOSY J, et al. Near-real-time regional troposphere models for the GNSS precise point positioning technique[J]. Measurement Science and Technology, 2013, 24(5): 055003. DOI:10.1088/0957-0233/24/5/055003.

[10]WIELGOSZ P, KRUKOWSKA M, PAZIEWSKI J, et al. Performance of ZTD models derived in near real-time from GBAS and meteorological data in GPS fast-static positioning[J]. Measurement Science and Technology, 2013, 24(12): 125802. DOI:10.1088/0957-0233/24/12/125802.

[11]HU G, ABBEY D A, CASTLEDEN N, et al. An approach for instantaneous ambiguity resolution for medium-to long-range multiple reference station networks[J]. GPS Solutions, 2005, 9(1): 1-11. DOI:10.1007/s10291-004-0120-8.

[12]PAN S G, MENG X, WANG S L, et al. Ambiguity resolution with double troposphere parameter restriction for long range reference stations in NRTK System[J]. Survey Review, 2015, 47(345): 429-437. DOI: 10.1179/1752270614Y.0000000144.

[13]PAZIEWSKI J. Precise GNSS single epoch positioning with multiple receiver configuration for medium-length baselines: methodology and performance analysis[J]. Measurement Science and Technology, 2015, 26(3): 035002. DOI: 10.1088/0957-0233/26/3/035002.

[14]ODIJK D, ARORA B S, TEUNISSEN P J G. Predicting the success rate of long-baseline GPS+ Galileo (partial) ambiguity resolution[J]. Journal of Navigation, 2014, 67(3): 385-401. DOI:10.1017/S037346331400006X.

[15]ODOLINSKI R, TEUNISSEN P J G, ODIJK D. Combined GPS+ BDS+ Galileo+ QZSS for long baseline RTK positioning[J]. Measurement Science and Technology, 2015, 26(4): 045801.

[16]PAZIEWSKI J, STEPNIAK K. New on-line system for automatic postprocessing of fast-static and kinematic GNSS data[C]// CYGAS D, TOLLAZZI T. International conference on environmental engineering (ICEE) selected papers. Vilnius, Lithuania: Vilnius Gediminas Technical University Press Technika, 2014:1-7. DOI: 10.3846/enviro.2014.235.

[17]KHODABANDEH A, TEUNISSEN P J G. Array-based satellite phase bias sensing: theory and GPS/BeiDou/QZSS results[J]. Measurement Science and Technology, 2014, 25(9): 095801. DOI:10.1088/0957-0233/25/9/095801.

[18]PAZIEWSKI J, WIELGOSZ P. Accounting for Galileo-GPS inter-system biases in precise satellite positioning[J]. Journal of Geodesy, 2015, 89(1): 81-93. DOI:10.1007/s00190-014-0763-3.

[19]International GNSS Service. GPS satellite ephemerides / Satellite & Station Clocks[OL]. [2017-04-06]. http://www.igs.org/productsl.

[20]WIELGOSZ P, PAZIEWSKI J, KRANKOWSKI A, et al. Results of the application of tropospheric corrections from different troposphere models for precise GPS rapid static positioning[J]. Acta Geophysica, 2012, 60(4): 1236-1257. DOI:10.2478/s11600-011-0078-1.

[21]BOHM J, MOLLER G, SCHINDELEGGER M, et al. Development of an improved empirical model for slant delays in the troposphere (GPT2w)[J]. GPS Solutions, 2015, 19(3): 433-441. DOI: 10.1007/s10291-014-0403-7.

[22]LI W, YUAN Y B, OU J K, et al. A new global zenith tropospheric delay model IGG trop for GNSS applications[J]. Chinese Science Bulletin, 2012, 57(17): 2132. DOI: 10.1007/s11434-012-5010-9.

[23]JIN S G, LUO O, REN C. Effects of physical correlations on long-distance GPS positioning and zenith tropospheric delay estimates[J]. Advances in Space Research, 2010, 46(2): 190-195. DOI:10.1016/j.asr.2010.01.017.

[24]SAASTAMOINEN J. Contributions to the theory of atmospheric refraction[J]. Bulletin Géodésique (1946-1975), 1973, 107(1): 13-34. DOI:10.1007/BF02522083.

[25]BOHM J, NIELL A, TREGONING P, et al. Global Mapping Function (GMF): A new empirical mapping function based on numerical weather model data[J]. Geophysical Research Letters, 2006, 33(7): 1-4. DOI: 10.1029/2005GL025546.

[26]唐衛(wèi)明, 劉經(jīng)南, 施闖, 等. 三步法確定網(wǎng)絡(luò) RTK 基準(zhǔn)站雙差模糊度[J]. 武漢大學(xué)學(xué)報(bào) (信息科學(xué)版), 2007, 32(4): 305-308. DOI:10.3321/j.issn:1671-8860.2007.04.007.

TANG W M, LIU J N, SHI C, et al. Threesteps 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. DOI:10.3321/j.issn:1671-8860.2007.04.007.

[27]BETZ J W. Engineering satellite-based navigation and timing: global navigation satellite systems, Signals, and Receivers[M]. Hoboken, New Jersey: John Wiley & Sons, 2015: 520-529.

[28]TEUNISSEN P J G. The least-squares ambiguity decorrelation adjustment: a method for fast GPS integer ambiguity estimation[J]. Journal of geodesy, 1995, 70(1-2): 65-82. DOI:10.1007/BF00863419.

[29]WANG J, STEWART M P, TSAKIRI M. A discrimination test procedure for ambiguity resolution on-the-fly[J]. Journal of Geodesy, 1998, 72(11): 644-653. DOI:10.1007/s001900050204.

Real-timeprecisepositioningusingdual-frequencyuncombineddata

LI Hui, ZHAO Lin, LI Liang, DING Jicheng

(College of Automation, Harbin Engineering University, Harbin 150001, China)

To satisfy the real-time precise positioning requirement of single-system single-baseline with dual-frequency data, a novel approach based on dual-frequency uncombined data is proposed. The single-system single-baseline and dual-frequency uncombined Real Time Kinematic (RTK) is realized based on inter-station-inter-satellite carrier phase and pseudorange measurement model. Based on the analysis of the residual errors of the double difference model, the state vector is established, and then the state prediction equation and the measurement equation are derived. The real-time transformation matrix of the state vector is established at each epoch, and the weights of observations including pseudorange and carrier phase are adjusted by a random model. Then extended Kalman filter is used to estimate the position of the user receiver. The approach is tested with several groups of medium-range and long-range single-baseline data which are the actual collection of satellites in the experiment. The effectiveness of the approach is verified by examining the three-dimensional errors of the positioning result as well as the success rate of the ambiguities. It is shown that when the baseline is 135.6 km, the positioning accuracy of the proposed approach can reach the centimeter level with the 97.3% success rate of the ambiguities, and the CEP95 positioning errors of east and north and sky directions are 1.35 cm, 1.84 cm, 7.08 cm respectively. This approach makes full use of the advantages of differential technology to eliminate the relative errors independent of distance and effectively avoid the noise caused by the linear combination of the multiple frequency measurements and can satisfy real-time precise positioning for medium and long single-baseline.

real-time precise positioning; dual-frequency; uncombined; medium and long baseline; extended Kalman filter

10.11918/j.issn.0367-6234.201703023

P223

A

0367-6234(2017)11-0122-08

2017-03-06

國家自然科學(xué)基金重點(diǎn)項(xiàng)目(61633008)、國家自然科學(xué)基金青年科學(xué)基金(61304235、61401114)資助的課題.

李 慧(1987—)，女，講師，博士研究生；趙 琳(1968—)，男，教授，博士生導(dǎo)師

李 慧，E-mail: lihuiheu@hotmail.com

(編輯苗秀芝)