王耀鼎,劉文祥,王飛雪

(國防科學(xué)技術(shù)大學(xué) 電子科學(xué)與工程學(xué)院,長沙 410073)

?

基于多項式擬合法北斗三頻周跳探測改進技術(shù)研究

王耀鼎,劉文祥,王飛雪

(國防科學(xué)技術(shù)大學(xué) 電子科學(xué)與工程學(xué)院,長沙 410073)

傳統(tǒng)三頻線性組合的周跳檢測量一般通過歷元間作差產(chǎn)生,這樣要求組合后必須保證不含幾何距離,使得線性組合的方式存在一個強約束,很多具有優(yōu)良特性的單純載波相位組合無法使用。而將多項式擬合法應(yīng)用到三頻線性組合載波后,則無需限制組合無幾何約束。靜態(tài)場景下,本文提出了一種基于多項式擬合法的三頻周跳探測改進算法,通過定義電離層放大因子、鐘差放大因子以及組合噪聲/波長放大因子,選擇出三組線性組合,以此來探測周跳。實測結(jié)果表明,本算法可以有效的探測出1周的小周跳。

北斗系統(tǒng);周跳探測;三頻線性組合;多項式擬合法;放大因子

0 引 言

利用載波相位測量值進行高精度導(dǎo)航定位的前提是對周跳進行預(yù)處理。北斗三頻數(shù)據(jù)的出現(xiàn),為周跳的探測與修復(fù)提供了更多具有電離層延遲小、波長長、組合噪聲小的組合[1]。國內(nèi)外很多專家學(xué)者均對三頻周跳探測與修復(fù)技術(shù)展開了研究。

文獻[2]和文獻[3]通過檢查大周跳、小周跳、特殊周跳三級檢測法最終將所有周跳組合檢測出來。文獻[4]利用三頻偽距/載波組合探究了周跳探測算法在GEO、IGSO、MEO三種不同衛(wèi)星下的探測效果。文獻[5]利用空間搜索法探測周跳。文獻[7]通過定義波長參數(shù)、電離層參數(shù)、噪聲參數(shù),并以此為目標(biāo)函數(shù)進行周跳探測。文獻[8]和文獻[9]通過電離層兩次歷元間作差,削弱了電離層延遲對周跳探測精度的影響。文獻[10]分析了不同采樣間隔下無幾何相位組合和偽距/相位組合的探測性能。文獻[11]推導(dǎo)了組合觀測值電離層比例因子和噪聲比例因子的關(guān)系。文獻[12]提出將觀測數(shù)據(jù)進行平滑之后再進行周跳探測。文獻[13]提取了偽距多路徑誤差,改善探測精度。文獻[14]探測方法對采樣間隔較大、電離層較活躍的場景下有較強的適應(yīng)性。

靜態(tài)場景下,可以考慮用多項式擬合星地距離。通過線性組合,擴大組合后波長,降低了接收機鐘差對探測結(jié)果的影響。實測結(jié)果表明,本算法可以有效的探測出1周的小周跳。

1 多項式擬合法原理

多項式擬合法基于觀測量平滑變化,此時對于觀測量可以用一個多項式進行擬合。由于星地距離4次差分后殘差已接近于白噪聲,因此可利用4階多項式對觀測量進行擬合。其擬合模型為

φi= a0+a1(ti-t0)+a2(ti-t0)2+

a3(ti-t0)3+a4(ti-t0)4

(1)

式中: φi為時刻的載波相位； ai,i=0,…,4為多項式擬合系數(shù)； t0為擬合窗口起始時間； ti為擬合窗口內(nèi)i時刻。目前文獻一般利用多項式擬合法直接對單頻原始載波相位進行擬合。對于北斗三頻載波相位,進行線性組合后,其組合結(jié)果并不破壞原始載波相位的平滑性。因此,本文提出將多項式擬合法應(yīng)用于三頻線性組合后的載波相位,對組合后載波相位進行多項式擬合。

2 實數(shù)線性組合約束

首先寫出原始載波相位觀測值表達式:

φ= λφ=r+c(δtu-δts)-I+

T+εφ-λN,

(2)

φi,j,k=iφ1+jφ2+kφ3.

(3)

波長為

(4)

定義頻率無關(guān)量為

s=r+c(δtu-δts)+T.

(5)

則將式(2)展開,得:

∑tempρ·(-Ntemp),

tempρ=i,j,k; temp=1,2,3.

(6)

可以通過限制組合系數(shù)(i,j,k),使得式(6)中的某一項置零。

考慮周跳探測的誤差來源,提出如下無電離層、無鐘差、最小組合噪聲三種約束。

2.1 無電離層約束

限制組合后電離層系數(shù)為0,經(jīng)過化簡,得無電離層組合約束為

(7)

對于北斗系統(tǒng),可得:

i+1.293j+1.231k=0.

(8)

滿足上式的實數(shù)解即可滿足無電離層約束。

定義式(8)為電離層放大因子αi,j,k,其表達式為

αi,j,k=i+1.293j+1.231k.

(9)

其實際上構(gòu)成了一個以(i,j,k)為自變量的空間中的平面,則空間中越靠近該平面的點,其對應(yīng)組合后的電離層放大因子越小。

2.2 無鐘差約束

限制組合后鐘差系數(shù)為0,即非頻率部分系數(shù)為0,得無鐘差約束為

if1+jf2+kf3=0,

(10)

代入數(shù)值,得:

763i+590j+620k=0.

(11)

滿足上式的(i,j,k)可以消除鐘差,其同樣定義了一個平面,空間中越靠近該平面的點其鐘差影響越小。

定義鐘差放大因子為βi,j,k,其值為上式相對于f1的放大倍數(shù)

(12)

代入數(shù)值,得

βi,j,k=i+0.7733j+0.8126k.

(13)

該值越小,鐘差影響越小。

2.3 最小組合噪聲約束

定義以周為單位的組合噪聲放大因子為

(14)

由上式可知,(i,j,k)在空間中越靠近原點,噪聲放大因子越小。

(15)

上式值越小,說明噪聲放大程度越小于波長放大程度,則性能越好。

3 最優(yōu)整數(shù)組合

為保證組合后的模糊度為整數(shù),載波組合系數(shù)必須為整數(shù)。因此需要在整數(shù)范圍內(nèi)搜索能夠最大限度接近第2節(jié)描述的三種約束條件的整數(shù)組合。同時,為便于周跳探測,組合后的波長應(yīng)盡可能長。將整數(shù)組合進行分類,定義組合系數(shù)之和為M,即:

M=i+j+k.

(16)

同時,定義Mx為滿足i+j+k=±x的(i,j,k)的類,由此可以得到多類Mx,由文獻[1]可知,M0類和M1類具有優(yōu)良的特性,而在進行周跳探測時,只需要構(gòu)造三組線性無關(guān)的組合即可。因此,研究M0類和M1類中的組合已經(jīng)滿足要求,故不需再考慮其他類。

3.1 M0類

此時i+j+k=0,由于噪聲放大因子不可過大,因此,限制(i,j,k)∈[-5,5]進行搜索,挑選其中的較優(yōu)組合如表1所示.

表1 M0類組合

為提高多項式擬合法精度,應(yīng)選擇電離層放大因子、鐘差放大因子及噪聲放大因子均較小的量。

由表1可知,(0,-1,1)是最佳選擇,因此,將其確定為第一個線性組合。對于多項式擬合法,接收機鐘差對其影響較大,因此需選擇鐘差放大因子較小的組合。考察組合(1,4,-5),其鐘差放大因子在剩余組合中最小,雖然其噪聲放大因子較大,但其噪聲/波長放大因子γi,j,k較小,因此其性能仍然較好,故選擇其為第二個線性組合。

由于限制i+j+k=0,因此M0類中最多有2個線性獨立的組合,故第三組組合需要從M1類中進行選擇。

3.2 M1類

此時i+j+k=±1,同樣,為防止噪聲放大因子過大,限制(i,j,k)∈[-5,5]進行搜索,挑選其中的較優(yōu)組合如表2所示。

表2 M1類組合

由上表可知,噪聲放大因子最小的組合是(1,0,0)(0,1,0)(0,0,1),但是其噪聲/波長放大因子γi,j,k較大,故不采用。對于組合(-4,2,3),其鐘差放大因子及噪聲/波長放大因子均最小,雖然電離層放大因子較大,但是由于在產(chǎn)生檢測量時用多項式擬合法,因此對于電離層變化不太大的情況,其電離層延遲變化可以很好的通過多項式擬合法消弱其影響。綜上,第三組線性組合選擇為(-4,2,3)。

4 改進算法流程

由第三節(jié)分析可知,選擇的載波相位組合系數(shù)為(0,-1,1),(1,4,-5),(-4,2,3),假設(shè)其組合后的檢測量為L,則其周跳組合表達式為

(17)

式中:ΔNi為各頻點的周跳值; li為各線性組合的周跳檢測量。由于系數(shù)矩陣行列式為1,故矩陣是可容許的,則由上式可以解算出各個頻點的周跳值,其表達式為

(18)

以第一組線性組合為例,介紹改進算法周跳檢測量的產(chǎn)生過程。

根據(jù)式(1)和式(3),對組合后載波相位φi,j,k進行多項式擬合,窗口寬度取6,得到多項式系數(shù)。

預(yù)測第7個歷元的組合觀測值,其預(yù)測值為

a3(t7-t0)3+a4(t7-t0)4,

(19)

φ0,-1,1=0·φ1+(-1)·φ2+1·φ3,

(20)

則周跳檢測量l1為

(21)

另外兩組線性組合的周跳檢測量產(chǎn)生過程與之完全相同,不再贅述。

對li直接取整,即為周跳檢測值。

5 實測數(shù)據(jù)驗證

利用北斗一期監(jiān)測接收機2015年12月22日9時至15時的觀測數(shù)據(jù),采樣間隔為1 s.選取C09衛(wèi)星進行研究。圖1示出了三種組合的周跳檢測量情況。

為測試改進算法對特殊周跳(1,1,1)的檢測性能,在三個頻點400歷元處均加入1周小周跳,結(jié)果如圖2所示。

圖1 三種組合的周跳檢測量

圖2 改進算法三種組合的周跳檢測量

理論上三種組合在400歷元處的檢測量為0,0,1,實際結(jié)果為-0.0079,0.0231,1.0470,取整后結(jié)果符合理論值。

隨機在某些歷元加入模擬周跳,其周跳探測情況如表3所示。

對應(yīng)的周跳探測結(jié)果圖，如圖3所示。

圖3 模擬周跳三種組合的檢測量

表3 模擬周跳探測與修復(fù)結(jié)果

6 結(jié)束語

在靜態(tài)場景下,星地距離可以很好的用多項式進行擬合。利用多項式擬合法對組合后載波進行周跳探測,使得組合后觀測值即保持了以周為單位的載波的組合系數(shù)的整數(shù)特性,又可以不再限制幾何距離無關(guān),這為探討新的具有優(yōu)良特性的線性組合方式開辟了道路。通過限制鐘差放大因子和噪聲/波長放大因子,選擇了三組載波幾何相關(guān)組合,取代了偽距/相位組合,充分利用載波測量的高精度特性,提高了探測精度。通過加入模擬周跳,驗證本算法可以探測小周跳、特殊周跳、大周跳等周跳類型。該方法可應(yīng)用于北斗三頻非差觀測數(shù)據(jù)中各頻點周跳的探測。對于動態(tài)場景,由于很難用多項式進行擬合,故如何提高動態(tài)場景下的探測精度有待進一步研究。

[1] 張小紅,何錫揚.北斗三頻相位觀測值線性組合模型及特性研究[J].中國科學(xué),2015,45(5):601-10.

[2] 曹新運,王堅.GPS三頻非差觀測值探測與修復(fù)周跳[J].武漢大學(xué)學(xué)報(信息科學(xué)版),2014,39(4):450-457.

[3] CLARA DE LACY M, REGUZZONI M,SANSO F. Real-time cycle slip detection in triple-frequency GNSS[J].GPS Solutions,2012,16(3):353-362.

[4] 劉俊,何秀鳳,劉炎雄.北斗三頻非差觀測數(shù)據(jù)的周跳探測與修復(fù)[J].大地測量與地球動力學(xué),2014,34(4):117-22.

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

[6] 常志巧,劉利,何海波.基于最優(yōu)觀測組合的三頻周跳探測與修復(fù)[C]//CSNC2010第一屆中國衛(wèi)星導(dǎo)航學(xué)術(shù)年會論文集,2010.

[7] 韓紹偉.GPS組合觀測值理論及應(yīng)用[J].測繪學(xué)報,1995,24(2):8-13.

[8] 黃令勇,翟國君,歐陽永忠,等.削弱電離層影響的三頻TurboEdit周跳處理方法[J].測繪學(xué)報,2015,44(8):840-847.

[9] 黃令勇,翟國君,歐陽永忠,等.三頻GNSS電離層周跳處理[J].測繪學(xué)報,2015,44(7):717-725.

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

[11] 孫保琪,歐吉坤,盛傳貞,等.一種適于Compass周跳探測的三頻數(shù)據(jù)優(yōu)化組合[J].武漢大學(xué)學(xué)報(信息科學(xué)版),2010,35(10):1157-1160.

[12] 范建軍,王飛雪,郭桂蓉.GPS三頻非差觀測數(shù)據(jù)周跳的自動探測與改正研究[J].測繪科學(xué),2006,31(5):24-26.

[13] 肖國銳,隋立芬,戚國賓,等.大地測量與地球動力學(xué)[J].大地測量與地球動力學(xué),2015,35(4):671-675.

[14] 王興,劉文祥,李柏渝,等.北斗三頻非差觀測數(shù)據(jù)周跳探測與修復(fù)[J].國防科技大學(xué)學(xué)報,2015:1-7.

A New Triple-Frequency Cycle Slip Detection Algorithm Based on Polynomial

WANG Yaoding,LIU Wenxiang,WANG Feixue

(CollegeofElectronicScienceandEngineering,NationalUniversityofDefenseTechnology,Changsha410073,China)

Traditional methods of Tri-frequency cycle slip detection are based on the difference of epochs, which requires that the combination do not contain geometric distance. And that causes a strong constraint in the combination. When the polynomial fitting applies to the Tri-frequency cycle slip detection, the constraint is relieved. In static case, this novel applies the polynomial to Tri-frequency cycle slip, and defines ionosphere amplification factor、clock bias amplification factor and noise amplification. The results shows that this method can detect at least 1 slip.

BeiDou; cycle slip; tri-frequency Linear combination; polynomial; amplification factor

10.13442/j.gnss.1008-9268.2016.05.016

2016-07-11

P228.4

A

1008-9268(2016)05-0079-05

王耀鼎 (1992-),男,山東濱州人,碩士,研究方向為星基導(dǎo)航與定位技術(shù)。

劉文祥 (1981-),男,江西宜春人,博士,講師,研究方向為星基導(dǎo)航與定位技術(shù)。

王飛雪 (1971-),男,福建長汀人,博士,教授,博士生導(dǎo)師,主要研究方向為衛(wèi)星導(dǎo)航定位、擴頻信號處理、電子系統(tǒng)抗干擾。

聯(lián)系人： 王耀鼎 E-mail: wangyaodingsdu@126.com