李國梁, 錢雨陽

(1.中南勘察設(shè)計院集團有限公司,湖北 武漢 430071 2. 山東科技大學(xué) 測繪科學(xué)與工程學(xué)院,山東 青島 266590)

0 引 言

北斗衛(wèi)星導(dǎo)航系統(tǒng)(BDS)是我國根據(jù)我國基本國情以及基本需求而自主設(shè)計研發(fā)的導(dǎo)航系統(tǒng),分三步建設(shè)完成[1-4]:試驗系統(tǒng)-區(qū)域系統(tǒng)-全球系統(tǒng)．自從1983年提出建設(shè)BDS以來,我國BDS經(jīng)歷了翻天覆地的變化,2000年我國發(fā)射第一課BDS衛(wèi)星,開始我國BDS建設(shè)之旅,到2007年完成第一階段的建設(shè),共四顆衛(wèi)星組成北斗一號(BDS-1)．到2012年底,完成北斗二號(BDS-2)的建設(shè),共14顆衛(wèi)星組網(wǎng),BDS-2可以向亞太地區(qū)提供導(dǎo)航與服務(wù)．自2015年發(fā)射第一顆北斗三號(BDS-3)試驗衛(wèi)星,標志著我國進入第三階段的建設(shè),預(yù)計2020完成BDS-3的建設(shè),屆時將有35顆衛(wèi)星組網(wǎng),向全球提供導(dǎo)航定位服務(wù)．

本文利用含有BDS-3的IGS跟蹤站數(shù)據(jù)進行標準單點定位分析,評估了BDS-3的數(shù)據(jù)質(zhì)量,并且其定位精度與BDS-2和GPS單系統(tǒng)定位精度進行對比分析．

1 BDS-3衛(wèi)星情況

本文選取的2018年第200天的數(shù)據(jù),根據(jù)觀測數(shù)據(jù)分析以及BDS-3發(fā)射列表,總共可以接收到8顆BDS-3衛(wèi)星,不算之前發(fā)射的試驗衛(wèi)星,BDS-3衛(wèi)星軌道類型都為中地球軌道(MEO),具體發(fā)射時間以及類型如表1所示.

表1 BDS-3衛(wèi)星發(fā)射一覽表

相比于BDS-2衛(wèi)星,BDS-3衛(wèi)星增加了3個民用頻點,新衛(wèi)星頻點的增加使BDS將會廣泛存在多系統(tǒng)頻點重合,這為實現(xiàn)不同系統(tǒng)間的同頻組合定位提供了機會．BDS-3頻段的具體情況如表2所示.

表2 BDS-3衛(wèi)星信號一覽表

2 BDS單點定位原理

BDS與GPS不同,BDS由三種衛(wèi)星星座構(gòu)成,GEO、MEO和IGSO,因此在計算BDS衛(wèi)星位置時也會與GPS衛(wèi)星位置計算方法表現(xiàn)出差異性,其中BDS的MEO和IGSO的位置計算方法與GPS一致,只是在計算過程中要注意BDS與GPS引力常數(shù)不同,以及計算出來的衛(wèi)星位置所在坐標系不一致,而GEO則與其他兩個星座的計算方法存在差異[5-6]．MEO和IGSO的衛(wèi)星位置計算公式如下[7]:

(1)

式中：(Xk,Yk,Zk)為衛(wèi)星位置坐標;Ωk為升交點赤經(jīng);(xk,yk)為衛(wèi)星平面坐標;ik為軌道傾角;

GEO衛(wèi)星位置的計算如下:

GEO衛(wèi)星在自定義坐標系統(tǒng)中的坐標為

(2)

GEO衛(wèi)星在CGCS2000坐標系中的坐標為

(3)

在計算出衛(wèi)星的軌道坐標之后,讀取觀測值文件,對獲取的偽距觀測值進行預(yù)處理,初始化數(shù)據(jù):

(4)

接下來計算衛(wèi)星發(fā)射信號初始時刻,經(jīng)過迭代之后,重新算出某一歷元的偽距觀測方程:

(5)

假設(shè)接收機起始近似坐標為(X0,Y0,Z0),計算得到接收機位置坐標改正量為(δx,δy,δz),所以接收機坐標為

(6)

在計算得到坐標位置之后,需要對其精度進行分析,一般只對其外符合精度進行分析,計算公式如下:

(7)

(8)

3 實驗與結(jié)果分析

3.1 數(shù)據(jù)質(zhì)量分析

BDS數(shù)據(jù)質(zhì)量是后續(xù)導(dǎo)航定位和授時(PNT)服務(wù)的基礎(chǔ),因此對BDS數(shù)據(jù)質(zhì)量進行分析是非常重要的內(nèi)容．常規(guī)數(shù)據(jù)質(zhì)量的評估指標有數(shù)據(jù)完整率、信噪比、多路徑和周跳比等[8]．為了對BDS-3數(shù)據(jù)質(zhì)量進行全面分析,本文選取iGMAS發(fā)布的ZHON、PETH和DWIN等3個全球跟蹤站,2018年第200天的數(shù)據(jù),采樣率為30 s,評估C20、C22、C27、C30這四顆衛(wèi)星的數(shù)據(jù)質(zhì)量,主要從信噪比和多路徑兩方面進行分析．對于數(shù)據(jù)處理,則選取3個站其中一個站的信噪比成圖．

圖1 C20信噪比 圖2 C22信噪比

圖3 C27信噪比 圖4 C30信噪比

如圖1～4所示,所選取衛(wèi)星的信噪比都隨著高度角的增加而增加,起始的衛(wèi)星信噪比基本都大于30 dB-Hz,當(dāng)高度角達到60°時,信噪比趨于平穩(wěn),并且每顆衛(wèi)星的信噪比趨于平穩(wěn)后都大于50 dB-Hz．

圖5 C20多路徑 圖6 C22多路徑

圖7 C27多路徑 圖8 C30多路徑

圖9 BDS-2多路徑

如圖5～9所示,每顆衛(wèi)星的多路徑都在2 m以內(nèi),隨著高度角的增加而減小,并且在BDS-2衛(wèi)星中觀測到的系統(tǒng)誤差在BDS-3衛(wèi)星中并未觀測到,這對提高偽距單點定位精度有著非常重要的意義．

3.2 BDS-3標準單點定位精度分析

為了對BDS-3衛(wèi)星的標準單點定位精度進行分析,本文選取上述3個站中的PETH站作為單點定位精度解算分析站,利用GAMIT對該站進行網(wǎng)解算,水平精度可以達到mm,豎直精度可以達到cm.將解算結(jié)果作為真值,以便分析標準單點定位結(jié)果的外符合精度．為了排除BDS-2衛(wèi)星對定位結(jié)果的影響,在計算過程中將BDS-2衛(wèi)星剔除,只計算BDS-3衛(wèi)星標準單點定位結(jié)果．相反將BDS-2標準單點定位結(jié)果作為對比精度,在計算時也需要將BDS-3的衛(wèi)星剔除．在選取測站的全天數(shù)據(jù)中,截取其中20個小時的數(shù)據(jù)進行分析．

如圖10所示,BDS-2的衛(wèi)星數(shù)要多于BDS-3的衛(wèi)星數(shù),因為目前BDS-2發(fā)射的衛(wèi)星要比BDS-3多6顆．

圖10 衛(wèi)星數(shù)

如圖11和表3所示,BDS-2和BDS-3三個方向的定位誤差都在3 m以內(nèi),三個方向的定位精度是E方向的最好,U方向的最差,平面精度明顯優(yōu)于高程方向的精度．BDS-3相對BDS-2,每個方向的精度都有所提高,E方向精度提高了5%,N方向提高了14%,U方向提高了8%．雖然衛(wèi)星數(shù)影響SPP的精度,但是在進行數(shù)據(jù)質(zhì)量分析時發(fā)現(xiàn),在BDS-3衛(wèi)星中不存在系統(tǒng)偏差,此外,BDS-3衛(wèi)星的性能在其他方面相對于BDS-2都有很大的改進,導(dǎo)致BDS-3定位精度更高．相信隨著BDS-3衛(wèi)星的增多,BDS-3的SPP將會更優(yōu)于BDS-2．

圖11 北斗二號和北斗三號標準單點定位(SPP)定位誤差

表3 BDS-2和BDS-3的SPP定位誤差RMS統(tǒng)計

4 結(jié) 論

本文闡述了我國BDS-3衛(wèi)星的現(xiàn)狀,結(jié)合IGS連續(xù)跟蹤站的數(shù)據(jù),對BDS-2和BDS-3的數(shù)據(jù)質(zhì)量以及SPP精度進行分析,結(jié)果表明:

1)BDS-3衛(wèi)星數(shù)據(jù)質(zhì)量優(yōu)良,并且沒有觀測到存在于BDS-2衛(wèi)星中的系統(tǒng)偏差;

2)單系統(tǒng)實驗結(jié)果表明,BDS-3的定位精度優(yōu)于BDS-2的定位精度,在三個方向上都有所提升．