邱 蕾,劉正平,伍 岳,王海軍
(1.深圳市地籍測繪大隊,廣東 深圳 518000;2.湖北省地質環境總站,湖北 武漢 430034;3.三峽大學 土木與建筑學院,湖北宜昌443002;4.二炮士官學校作戰保障教研室,山東 青州262500)
美國從1999年實施GPS現代化措施,其中一個主要的措施是在L2載波上增加新的民用碼L2C碼。截止目前,美國已經發射了9顆現代化的GPS Bl ock IIR-M衛星,該衛星上調制了L2C碼,9顆衛星 分 別 是 PRN01、PRN05、PRN07、PRN12、PRN15、PRN17、PRN25、PRN29 和 PRN31。為了說明新增加的L2C碼以及L2載波的優勢,本文利用實際觀測到的9顆衛星的觀測數據,先分別從L2C碼的信號結構、多路徑效應和信噪比進行分析,然后從求解整周模糊度方面綜合進行分析。
GPS現代化后,GPS偽距觀測值主要有兩種,分別是調制在L1載波上的C/A碼和調制在L2載波上的L2C碼。C/A碼的碼元長度為1 023碼元,周期為1 ms,碼率為1.023 MHz,而L2C碼采用新的設計結構,包含了民用中等長度碼(CM)和民用長度碼(CL)碼。CM碼的長度為10 230碼元,周期為20 ms,并且調制有速率為25 bit/s的導航數據。CL碼的長度為767 250碼元,周期為1 500 ms,不含調制數據,碼率為1.023 MHz。
將調制了導航數據的CM碼與未調制數據的CL碼通過時分復用(TDM)的方式,就可以獲得新的L2C碼。正是因為使用了TDM方式,L2C碼的數據解調門限與C/A碼相比,降低了5 d B,這種性能使得L2C碼在信號較微弱的情況下,也可以捕捉到較好的信號。
多路徑效應檢核可采用式(1)和式(2)進行。

式中,C/A和L2C是雙頻偽距觀測值,φ1和φ2是L1載波和L2載波相位觀測值,MPCA和MPL2C是C A和L2C偽距的多路徑效應,為L1載波和L2載波的頻率的平方比。
為了分析新增加的L2C碼的數據質量,本文利用式(1)和式(2)對實際采集到的觀測數據進行了計算并且將結果繪制成表1和表2。數據采樣時間為1 s,數據采集時間為2012年083日—090日連續8 d的觀測數據。表1為測站連續8 d的多路徑效應每天平均中誤差。表2列出了8 d的平均中誤差。

表1 測站1連續8 d的多路徑統計結果m

表2 連續8 d的多路徑效應平均中誤差m
從表1中可以看出,PRN01衛星的觀測質量整體較差,PRN15的多路徑效應在C/A碼和L2C碼上最為接近,PRN01的多路徑效應在C/A碼和L2C碼上的多路徑效應相差最大為0.17 m。從表2可以看出,采取現代化措施后的GPS衛星,C/A碼的多路徑效應平均中誤差為0.31 m,L2C碼的多路徑效應平均中誤差為0.33 m,平均中誤差較為接近。
因此,利用GPS現代化后的衛星得到的L2C碼的多路徑效應與C/A碼的多路徑效應比較接近,兩者的數據觀測質量非常接近。
信噪比(SNR)是指接收機的載波信號強度與噪聲強度的比值,單位為d B-Hz,如式(3)所示。在GPS數據觀測質量中,信噪比可以用來分析接收機接收到的數據質量。一般信噪比越高,觀測數據質量越好。

式中:k為玻茲曼常數,S為信號平均功率,T為開氏溫度。
選取2012年年積日083日的全天觀測數據,數據采樣率為1 s,共86 400個歷元。圖1左圖為調制了 L2C碼的PRN25、PRN31、PRN1、PRN17及PRN1衛星的L1載波和L2載波的信噪比變化;圖1右圖為調制了L2C碼的PRN12、PRN29、PRN15和PRN17衛星的L1載波和L2載波的信噪比變化。圖1中,左縱坐標代表信噪比,單位為d B-Hz,右縱坐標為高度角變化,單位為(°)。圖中曲線變化為衛星的信噪比變化,虛線變化為衛星的高度角變化。
圖2為未調制有L2C信號的上述GPS衛星的 信噪比隨高度角變化圖。圖中單位與圖1一致。

圖1 調制L2C信號的GPS衛星信噪比圖

圖2 未調制L2C信號的GPS衛星信噪比圖
圖1 中L2載波的信噪比值變化范圍為20~40 d B-Hz,圖2中L2載波的信噪比值變化范圍為12~40 d B-Hz。所有衛星的L1載波的信噪比值變化范圍為23~53 d B-Hz。L1載波和L2載波信噪比的變化趨勢與高度角密切相關。高度角越大,信噪比越高,高度角越低,信噪比越低。
圖1中,L1載波和L2載波信噪比差值最大的為PRN29,其L1載波的信噪比變化范圍為25~45 d B-Hz,L2載波的信噪比變化范圍為20~30 d BHz,其差值在15 d B-Hz內。圖2中,L1載波和L2載波的差值明顯比圖1中的差值要大。同時,在高度角變化較快的地方,加載了L2C碼的L2載波信噪比的變化比未加載L2C碼的L2載波信噪比的變化更緩慢。
上述實驗結果表明,加載了L2C碼的L2載波信噪比比未加載L2C碼的L2載波以及L1載波的信噪比明顯提高。隨著高度角的變化,其信噪比變化較緩慢,加載了L2C碼的L2載波的抗干擾能力明顯得到提高。
載波相位觀測數據質量的好壞直接影響到整周模糊度固定的效率。在GPS現代化之前,由于L1載波的觀測數據質量較好,在單個頻率解算時通常采用L1載波固定整周模糊度而未采用L2載波數據。通過上面分析可知,GPS現代化后L2載波的觀測質量明顯提高。本實驗利用現代化后的L2載波與L1載波在固定整周模糊度效率上進行對比來說明L2載波的觀測質量得到了明顯提高。
在GPS基線解算中,通常采用式(4)進行雙差整周模糊度浮點解的解算,可以看出,雙差整周模糊度浮點解的解算受到多種因素的影響,包括多路徑效應和噪聲。

式中:▽ΔN為雙差整周模糊度,▽Δρ為雙差偽距觀測值,▽Δφ為雙差載波相位觀測值,▽Δs為雙差后的剩余誤差。
實驗中采用兩個測站的觀測數據,觀測時間為24 h,采樣率為1 s。兩個測站的坐標精確已知,接收機類型一致。將所有的數據進行處理后,只利用具有L2C觀測值的GPS衛星對按照式(4)計算雙差整周模糊度的浮點解。選取其中的3組衛星對5-15、5-29和5-12計算L1載波和L2載波的單歷元雙差整周模糊度浮點解,將雙差整周模糊度浮點解與其固定解真值的差值變化見圖3、圖4和圖5,紅色為L1載波,綠色為L2載波。各圖的橫坐標為觀測歷元,單位為s,左縱坐標為差值,單位為周,所繪制曲線為L1載波和L2載波雙差整周模糊度的浮點解與其真值的差值;右縱坐標為衛星的高度角,單位為度,直線為衛星的高度角。

圖3 衛星對5-15 L1及L2載波的雙差模糊度浮點解與真值差值

圖4 衛星對5-29L1及L2載波的雙差模糊度浮點解與真值差值
從各圖可以看出,衛星對5-15、衛星對5-29和衛星對5-12的L1載波和L2載波的模糊度浮點解與真值的差值均在(-0.2,0.4)周內。當衛星高度角在20°~40°內,L1載波與L2載波的模糊度浮點解與真值的差值僅在0.2周內。利用L1載波和L2載波計算模糊度浮點解的效果基本相同,都可以較好地計算出模糊度的浮點解,為進一步固定整周模糊度縮小搜索空間。因此,加載了L2C碼的L2載波的數據觀測質量與L1載波的數據觀測質量較為一致,模糊度浮點解與真值較為接近,能夠較好地提高計算整周模糊度的浮點解的精度,縮小模糊度的搜索空間。

圖5 衛星對5-12L1及L2載波的雙差模糊度浮點解與真值差值
隨著GPS現代化的推進,新增加的民用碼L2C碼對用戶將具有較大的優勢。從本文的實驗中可以看出,L2C碼的多路徑誤差和C/A碼的多路徑誤差較一致;加載了L2C碼的L2載波的信噪比比原L2載波的信噪比明顯提高,同時,隨著高度角的變化,信噪比變化較緩慢。通過單歷元整周模糊度解算實驗,說明了加載了L2C的L2載波在計算模糊度時與L1效果較為一致,能達到L1載波的解算效果。
[1] 王澤民,邱蕾,孫偉.GPS現代化后L2載波的定位精度研究[J].武漢大學學報:信息科學版,2008,33(8):779-782.
[2] 伍岳,邱蕾.網絡RTK模式下多頻載波相位觀測值解算整周模糊度[J].測繪工程,2013,22(4):1-4.
[3] 陸明泉,馮振明,張祺.GPSL2捕獲算法研究及性能分析[J].電子與信息學報,2012,32(2):2.
[4] LILIáNA SüKEOVá,MARCELO C S,RICHARD B,et al.GPS L2C Signal Quality Analysis[C].Pr oceedings of the 63rd Annual Meeting,2007,Cambridge,Massachusetts:232-241.
[5] 張笑凡,何宏,陳建峰.GPS技術中對流層延遲影響分析[J].測繪與空間地理信息,2014,37(2):187-188.
[6] 羅伙軍,李程,岳國棟.GPS數據質量檢查軟件的可視化開 發[J].測 繪 與 空 間 地 理 信 息,2014,37(2):162-164.
[7] 劉基余.GPS衛星導航定位原理與方法[M].北京:科學出版社,2003.