布金偉，左小清，金立新，李海強，周羽佳，常 軍

(1. 昆明理工大學國土資源工程學院，云南 昆明 650093； 2. 中鐵第一勘察設計院集團有限公司，陜西 西安 710043； 3. 甘肅鐵道綜合工程勘察院有限公司，甘肅 蘭州 730000； 4. 國家測繪地理信息局第一大地測量隊，陜西 西安 710054)

偽距單點定位具有定位速度快且不存在整周模糊度問題，原理簡單易于編程實現(xiàn)等優(yōu)點。目前，一些學者對GPS+BeiDou、GPS+GLONASS、GPS+Galileo、GPS+GLONASS+Galileo和GPS+GLONASS+BeiDou多系統(tǒng)組合偽距單點定位進行了相關的研究[1-5]，研究表明：多系統(tǒng)組合偽距單點定位可以提高定位的穩(wěn)定性和精度[6-8]。

然而，對于BeiDou+GLONASS+Galileo多系統(tǒng)組合SPP的研究相對較少。因此，我國BDS對多系統(tǒng)組合SPP解算性能的貢獻有待進一步研究。鑒于此，本文在現(xiàn)有研究的基礎上，分析BeiDou+GLONASS+Galileo之間的時間和坐標系統(tǒng)差異，并推導了BeiDou+GLONASS+Galileo多系統(tǒng)組合SPP的數(shù)學模型；然后從可見衛(wèi)星數(shù)、DOP值、不同高度角下定位性能和定位精度幾個方面分析其不同組合模式下(單系統(tǒng)：BeiDou、GLONASS、Galileo；雙系統(tǒng)：BeiDou+GLONASS、BeiDou+Galileo、GLONASS+Galileo；三系統(tǒng)：BeiDou+GLONASS+Galileo)的SPP解算性能，最終得出有益結論。

1 BeiDou+GLONASS+Galileo多系統(tǒng)組合SPP基礎理論

1.1 時間基準和坐標基準的統(tǒng)一

為解算測站的坐標，必須首先實現(xiàn)BeiDou+GLONASS+Galileo三系統(tǒng)的時間系統(tǒng)和坐標系統(tǒng)的統(tǒng)一。文獻[9]和文獻[10]指出：BeiDou時間基準采用BDST，GLONASS時間基準采用GLONASST，Galileo時間基準采用GST。對于3個系統(tǒng)時間基準之間的差異和聯(lián)系，文獻[11]已經(jīng)給出比較詳細的敘述和轉換公式，這里不再闡述。

文獻[12—14]指出，BeiDou坐標基準采用CGCS2000坐標系，GLONASS坐標基準采用PZ-90坐標系，Galileo坐標基準采用GTRF坐標系。鑒于WGS-84坐標系目前已經(jīng)較為成熟且被廣泛使用，本文考慮將各框架按照赫爾墨特變換模型解算出的轉換參數(shù)統(tǒng)一至WGS-84坐標系下。在滿足精度的前提下，CGCS2000、WGS-84和GTRF坐標系統(tǒng)的差異較小，三者之間的差異為cm級。對于單點定位而言，其定位精度在m級，因此在數(shù)據(jù)處理過程中，可不考慮三者坐標系統(tǒng)差異引起的結果偏差。而PZ-90與WGS-84坐標系的差異能夠達到20 m[11]，故需考慮PZ-90和WGS-84坐標之間的轉換，目前世界公認精度最高的轉換參數(shù)為俄羅斯飛行控制中心(Russian Mission Control Center,RMCC)利用全球激光跟蹤站測算得到的轉換參數(shù)[9,11]。

1.2 BeiDou+GLONASS+Galileo多系統(tǒng)組合SPP定位模型

在多系統(tǒng)組合定位時，在某歷元t可以同時觀測到多個系統(tǒng)的衛(wèi)星信號，因此對于BeiDou+GLONASS+Galileo組合SPP的觀測方程可表示為[15]

(1)

(2)

(3)

在觀測方程式(1)—式(3)的基礎上，設測站近似坐標為(x0,y0,z0)，并記δx=x-x0，δy=y-y0，δz=z-z0，然后在(x0,y0,z0)處用泰勒級數(shù)展開，可得到誤差方程

(4)

(5)

三系統(tǒng)組合偽距單點定位待估參數(shù)中含有3個坐標分量參數(shù)和3個接收機鐘差，當觀測的衛(wèi)星數(shù)大于偽距定位的待估參數(shù)個數(shù)時，可采用最小二乘法原理求解，此時必要觀測衛(wèi)星數(shù)為6顆。由最小二乘估計可得

(6)

式中，P為觀測值權陣。

從而求出測站的坐標為

[xyz]T=[x0y0z0]T+[δxδyδz]T

2 試驗分析

為了分析BeiDou+GLONASS+Galileo組合偽距單點定位解算性能，試驗數(shù)據(jù)選取MGEX跟蹤站的部分測站2016年9月6日和10月25日共2 d的觀測數(shù)據(jù)，數(shù)據(jù)采樣率為30 s。試驗方案采用：單系統(tǒng)(BeiDou、GLONASS、Galileo，簡寫為C、R、E)、雙系統(tǒng)(BeiDou+GLONASS、BeiDou+Galileo、GLONASS+Galileo，簡寫為CR、CE、RE)與三系統(tǒng)(BeiDou+GLONASS+Galileo，簡寫為CRE)組合共7種定位模式。

2.1 衛(wèi)星可見數(shù)和DOP值分析

為分析BeiDou+GLONASS+Galileo組合的衛(wèi)星可見數(shù)和DOP值，本文選取測站ANMG進行試驗，圖1給出了BeiDou+GLONASS+Galileo組合在不同截止高度角(高度角分別設為15°、30°、40°)下的衛(wèi)星可見數(shù)和DOP值比較。圖2給出了1 d內BeiDou+GLONASS+Galileo組合模式在不同截止高度角情況下的平均DOP值和平均可見衛(wèi)星數(shù)統(tǒng)計。從圖1和圖2可知，當截止高度角為15°時，1 d內BeiDou+GLONASS+Galileo三系統(tǒng)組合可見衛(wèi)星數(shù)均在10顆以上，平均可見衛(wèi)星數(shù)在14顆左右。DOP值較小且波動平緩穩(wěn)定，平均DOP值小于2，說明衛(wèi)星的空間幾何結構分布較好；當截止高度角為30°時，BeiDou+GLONASS+Galileo三系統(tǒng)組合衛(wèi)星可用性仍然較高，全天內平均可見衛(wèi)星數(shù)為9.68顆。平均HDOP、VDOP、PDOP、GDOP值均小于5，而此時在1 d 內大約有2 h DOP值波動較大；當截止高度角大于40°時，所有組合模式的DOP值波動都較大，相比較而言，BeiDou+GLONASS+Galileo三系統(tǒng)組合全天內DOP波動情況比其他組合模式稍好，主要是由于三系統(tǒng)組合在任何時間段內可見衛(wèi)星數(shù)都多于其他組合模式，提高了衛(wèi)星的可用性和可靠性。結合圖2 (e)分析可知，7種定位模式下，從1 d內能觀測衛(wèi)星的時間和衛(wèi)星數(shù)來看可以得出關系：BeiDou+GLONASS+Galileo>BeiDou+GLONASS>BeiDou+Galileo>BeiDou>GLONASS+Galileo>GLONASS>Galileo。由于文中所選的測站在觀測1 d內GLONASS和Galileo可見衛(wèi)星數(shù)較少，無法獲得可靠的DOP值，故針對BeiDou、BeiDou+GLONASS、BeiDou+Galileo、GLONASS+Galileo、BeiDou+GLONASS+Galileo 5種模式可以得出DOP值大小比較關系：BeiDou+GLONASS+Galileo

圖1 多種組合模式在不同截止高度角下的衛(wèi)星可見數(shù)和DOP值比較

圖2 平均DOP值和平均可見衛(wèi)星數(shù)

為了分析BeiDou+GLONASS+Galileo三系統(tǒng)組合的DOP值和可見衛(wèi)星數(shù)相比其他組合模式的改善情況，表1—表5分別給出了CRE組合較其他模式的DOP值改善統(tǒng)計結果和可見衛(wèi)星數(shù)增加百分比統(tǒng)計。由表1—表4可知，當截止高度角小于30°時，HDOP、VDOP、PDOP、GDOP三系統(tǒng)組合較其他組合模式分別改善10%～61%、12%～59%、11%～58%、13%～55%；當截止高度角為40°時，CRE三系統(tǒng)組合的HDOP、VDOP、PDOP、GDOP較CR和CE組合分別改善13.33%和27.78%、15.56%和25.49%、14.89%和27.27%、15.00%和27.14%；較C和RE組合改善率均大于45%；當截止高度角為45°時，CRE三系統(tǒng)組合的HDOP和VDOP較C組合分別改善55.81%和42.86%，CRE三系統(tǒng)組合垂直方向上的改善效果稍差于水平方向。由表5可知，無論何種情況的截止高度角，CRE三系統(tǒng)組合的可見衛(wèi)星數(shù)都多于單系統(tǒng)和雙系統(tǒng)組合，特別是在截止高度角大于30°時，全天內的可見衛(wèi)星數(shù)能保持在6顆以上，提高了定位的可靠性和穩(wěn)定性。

2.2 不同高度角下定位穩(wěn)定性分析

圖3給出了測站GMSD在單系統(tǒng)、雙系統(tǒng)和三系統(tǒng)7種模式下單點定位結果(X、Y、Z方向的偏差為DX、DY、DZ，點位偏差為DP，圖3(a)截止高度角為15°，圖3(b)截止高度角為25°，圖3(c)截止高度角為30°)。從圖3可以看出：①對于同一種定位模式，隨著高度角的增加，定位穩(wěn)定性有所降低；②在不同截止高度角下，BeiDou+GLONASS+Galileo三系統(tǒng)組合單點定位穩(wěn)定性都最優(yōu)，盡管截止高度角增加時(30°時)，X、Z方向的偏差和點位偏差均小于10 m，Y方向的偏差小于12 m，點位偏差小于10 m。

表1 HDOP值改善統(tǒng)計結果 (%)

表2 VDOP值改善統(tǒng)計結果 (%)

表3 PDOP值改善統(tǒng)計結果 (%)

表4 GDOP值改善統(tǒng)計結果 (%)

2.3 定位精度分析

為分析單系統(tǒng)(C)、雙系統(tǒng)(CR、CE)和三系統(tǒng)(CRE)4種模式單點定位的定位精度情況，圖4分別給出了1 d內2個測站(其中圖4(a)、(c)、(e)、(g)為GMSD測站統(tǒng)計結果，圖4(b)、(d)、(f)、(h)為ANMG測站統(tǒng)計結果)不同模式的單點定位精度直方圖。如圖4所示，橫軸表示定位的平面誤差，縱軸表示分布密度，圖中的曲線為密度曲線，從圖中可以明顯看出，對于單系統(tǒng)而言，GMSD測站BeiDou系統(tǒng)(如圖4(a)所示)的平面定位精度為12 m，并且定位精度主要分布在0～12 m范圍內，密度曲線的波峰位于3～7 m范圍內，ANMG測站BeiDou系統(tǒng)(如圖4(b)所示)的平面定位精度為8 m，定位精度主要分布在0～8 m范圍內，密度曲線的波峰位于2～5 m范圍內；對于雙系統(tǒng)而言，GMSD測站CR、CE組合系統(tǒng)(如圖4(c)、(e)所示)的平面定位精度分別為10 m、12 m，定位精度分別主要分布在0～10 m、1～12 m范圍內，密度曲線的波峰分別位于3～6 m、3～7 m范圍內，ANMG測站CR、CE組合系統(tǒng)(如圖4(d)、(f)所示)的平面定位精度分別為8 m、8 m，定位精度分別主要分布在0～8 m、0～8 m范圍內，密度曲線的波峰分別位于2～5 m、2～4 m范圍內；然而對于三系統(tǒng)而言，GMSD測站CRE系統(tǒng)(如圖4(g)所示)的平面定位精度為10 m，平面定位精度主要分布在0～8 m范圍內，密度曲線的波峰位于3～5 m范圍內，ANMG測站CRE系統(tǒng)(如圖4(h)所示)的平面定位精度為8 m，定位精度主要分布在0～8 m，密度曲線的波峰位于2～4 m范圍內；通過以上2個測站的分析，相比單系統(tǒng)和雙系統(tǒng)組合模式，三系統(tǒng)組合定位在1 d內的定位精度優(yōu)于10 m，定位精度較高且定位結果相對集中。

圖3 多種組合模式單點定位結果在X、Y、Z方向的偏差和點位偏差

圖4 1 d內4種模式單點定位精度直方圖

高度角C/顆CR/顆CE/顆RE/顆CRE/顆CRE/C/(%)CRE/CR/(%)CRE/CE/(%)CRE/RE/(%)15°6.8812.469.728.4115.29122.2022.7657.3781.8420°6.7011.269.197.0513.75105.1322.1349.5595.1225°5.939.727.815.6811.6195.8819.4248.55104.3030°5.208.246.654.499.6886.3317.5845.73115.8440°4.216.255.062.907.1168.7913.6740.32145.3645°3.455.054.092.235.6964.6712.6638.96154.64

3 結 論

本文通過對單系統(tǒng)(C、R、E)、雙系統(tǒng)(CR、CE、RE)與三系統(tǒng)(CRE)共7種定位模式進行單點定位解算和結果分析，得出如下結論：

(1) 針對文中所選的測站而言，當高度角達到40°時，在可見衛(wèi)星數(shù)方面，無論單系統(tǒng)、雙系統(tǒng)還是三系統(tǒng)組合可見衛(wèi)星數(shù)均小于10顆，1 d內大部分時間都無法觀測到衛(wèi)星，而BeiDou+GLONASS+Galileo三系統(tǒng)組合大于18 h能觀測平均6顆以上，BeiDou+GLONASS雙系統(tǒng)組合次之；在DOP值方面，CRE三系統(tǒng)組合的HDOP、VDOP、PDOP、GDOP較其他組合模式改善大約13%～45%，1 d內DOP值變化較其他組合模式波動小。

(2) 單系統(tǒng)模式中BDS系統(tǒng)定位穩(wěn)定性最優(yōu)，Galileo系統(tǒng)最差，雙系統(tǒng)模式中BeiDou+GLONASS、BeiDou+Galileo組合系統(tǒng)定位穩(wěn)定性最優(yōu)，GLONASS+Galileo系統(tǒng)最差，而BeiDou+GLONASS+Galileo組合系統(tǒng)定位穩(wěn)定性優(yōu)于其他6種模式。

(3) 從定位精度來看，相比單系統(tǒng)和雙系統(tǒng)組合模式，BeiDou+GLONASS+Galileo三系統(tǒng)組合定位在1 d內的定位精度較優(yōu)且定位結果相對集中。