基于方差分量估計的不同截止高度角下的組合單點定位

馮 彪，柴洪洲，王 敏，潘宗鵬，李 迪

(信息工程大學(xué)，河南 鄭州 450000)

基于方差分量估計的不同截止高度角下的組合單點定位

馮 彪，柴洪洲，王 敏，潘宗鵬，李 迪

(信息工程大學(xué)，河南 鄭州 450000)

BDS/GPS/GLONASS組合系統(tǒng)定位時，由于系統(tǒng)間衛(wèi)星測距精度的差異性，需要合理確定衛(wèi)星間權(quán)比，Helmert方差分量估計常被用于確定不同類觀測值間權(quán)比；而當(dāng)觀測值含有粗差時，Helmert方差分量估計定位結(jié)果容易被粗差污染或收斂失真，出現(xiàn)大的偏差。文中基于Helmert方差分量估計，引入等價權(quán)因子IGGIII函數(shù)，建立抗差Helmert方差分量估計權(quán)函數(shù)模型，對比分析其在低截止高度角10°、15°和20°下，在BDS/GPS/GLONASS組合系統(tǒng)定位中的應(yīng)用及定位精度，并討論分析在高截止高度角30°和40°下，組合系統(tǒng)和單系統(tǒng)BDS的定位精度。實驗結(jié)果表明：當(dāng)觀測值無明顯粗差時，Helmert方差分量估計和抗差分量估計的定位精度相當(dāng)，略低于高度角權(quán)函數(shù)的定位結(jié)果，點位精度RMS優(yōu)于2.5 m；含粗差時，抗差解定位精度最高；當(dāng)截止高度角為30°時，BDS單系統(tǒng)定位精度RMS優(yōu)于5 m，而組合系統(tǒng)RMS接近3 m；為40°時，組合系統(tǒng)平面精度RMS優(yōu)于2 m，三維精度RMS優(yōu)于6 m，而單系統(tǒng)不能定位。

BDS/GPS/GLONASS組合系統(tǒng)定位；截止高度角；Helmert方差分量估計；IGGIII等價權(quán)函數(shù)；抗差估計；定位精度

北斗衛(wèi)星導(dǎo)航系統(tǒng)(BDS)于2012年正式向我國及周邊地區(qū)提供區(qū)域服務(wù)，預(yù)示著全球衛(wèi)星導(dǎo)航系統(tǒng)(GNSS)進入多系統(tǒng)兼容、合作的發(fā)展時代。單一系統(tǒng)往往由于用戶位置和觀測條件受限，不能提供足夠的衛(wèi)星數(shù)以獲取較高精度的位置信息，滿足不了用戶需求，多系統(tǒng)聯(lián)合導(dǎo)航定位逐漸成為GNSS發(fā)展研究的熱點。但由于各衛(wèi)星導(dǎo)航系統(tǒng)星座和衛(wèi)星軌道分布不同，衛(wèi)星的測距精度存在差異以及系統(tǒng)的定位能力的差別，GNSS聯(lián)合定位時，再根據(jù)先驗精度確定觀測值間權(quán)比不再準(zhǔn)確可靠。因此，對于組合導(dǎo)航系統(tǒng)，各導(dǎo)航系統(tǒng)間和衛(wèi)星間的權(quán)重因子優(yōu)化和其導(dǎo)航定位性能是研究多系統(tǒng)組合導(dǎo)航定位的重要方面。

1942年，一種驗后方差估計的方法被提出，即Helmert方差分量估計。它是通過迭代計算自適應(yīng)的調(diào)整不同類或不等精度觀測值之間的權(quán)比，使權(quán)的分配更具合理性[1-2]。但由于基于最小二乘的方差分量估計不具有抗差性，若觀測值含有粗差，粗差會污染隨機模型，導(dǎo)致迭代收斂失真或不收斂，定位結(jié)果也偏離真值[3-4]。為了減弱粗差對Helmert方差分量估計結(jié)果的影響，可以應(yīng)用抗差Helmert方差分量估計的方法來合理調(diào)整3類觀測值之間的權(quán)比，降低粗差觀測值的權(quán)，削弱其對定位結(jié)果的影響[5-6]。關(guān)于組合系統(tǒng)的導(dǎo)航定位性能[7-11]，相關(guān)的分析比較都是在低截止高度角下10°或15°的，相對單系統(tǒng)BDS和GLONASS，其導(dǎo)航定位性能提升明顯，而GPS不明顯。但對于截止高度角較大的情況，如城市、峽谷地區(qū)等的應(yīng)用，組合系統(tǒng)的導(dǎo)航定位性能及與單系統(tǒng)的差別分析較少，特別對正處于建設(shè)中的北斗系統(tǒng)，目前的應(yīng)用能否滿足上述條件的導(dǎo)航定位服務(wù)。因而，本文主要研究分析基于方差分量估計的BDS/GPS/GLONASS組合系統(tǒng)在低截止高度角10°、15°和20°下，3種不同權(quán)函數(shù)的定位精度和抵御粗差的能力以及探討在高截止高度角30°和40°下，組合系統(tǒng)和單系統(tǒng)BDS的導(dǎo)航定位性能。

1 BDS/GPS/GLONASS組合系統(tǒng)理論

1.1 時空基準(zhǔn)統(tǒng)一

GPS、BDS、GLONASS 3個導(dǎo)航系統(tǒng)都是采用獨立的時間和坐標(biāo)系統(tǒng)。GPS采用GPS時GPST，起算歷元為協(xié)調(diào)世界時UTC的1980年1月6日0時；BDS采用北斗時BDT，起算歷元為UTC的2006年1月1日0時。GPST和BDT都是以原子時為時間基準(zhǔn)，是一個連續(xù)的時間系統(tǒng)，不存在跳秒。但由于閏秒和兩個維持時間系統(tǒng)維持的差異，GPST與BDT存在整14 s和微小同步誤差。GLONASS采用GLONASS時GLONASST，GLONASST是采用莫斯科地區(qū)的協(xié)調(diào)世界時UTC，與協(xié)調(diào)世界時UTC時存在3 h偏差。GLONASST也存在跳秒，且與UTC保持一致。

GPS采用WGS84坐標(biāo)系統(tǒng)，BDS采用CGCS2000坐標(biāo)系統(tǒng)，二者定義一致，僅框架實現(xiàn)有差異，同一點在兩個坐標(biāo)系統(tǒng)下的坐標(biāo)差異在cm級[12]。GLONASS采用PZ-90.02，與ITRF2000只存在原點的平移，3個軸的定向一致，兩坐標(biāo)系偏差維持在1～2 dm內(nèi)[13]。因此，鑒于偽距定位的精度低于cm級，組合系統(tǒng)單點定位不需進行坐標(biāo)轉(zhuǎn)換。

1.2 觀測方程

GPS、GLONASS和BDS組合單點定位的觀測方程為

2 組合系統(tǒng)觀測值隨機函數(shù)模型

2.1 Helmert方差分量估計

將各衛(wèi)星系統(tǒng)偽距觀測值分為3類觀測值，根據(jù)Helmert方差分量估計原理，其計算過程[14]如下：

1)第一次最小二乘平差時，GPS、GLONASS、BDS觀測值根據(jù)衛(wèi)星高度角先驗定權(quán)P1，P2，P3，先驗權(quán)比為

(2)式中，eli為衛(wèi)星i的高度角。

3)按下式進行方差分量估計

(3)式中:

N=ATPA=A1TP1A1+A2TP2A2+A3TP3A3=N1+N2+N3.

(4)

其中:n1為GPS觀測值的個數(shù),n2為GLONASS觀測值的個數(shù)，n3為BDS觀測值的個數(shù)。

4)重新定權(quán)

(5)

針對上述過程中易出現(xiàn)迭代不收斂或收斂失真現(xiàn)象，分析原因一是可能在于各類觀測值的個數(shù)不均衡，如果一類觀測值個數(shù)太少，Helmert方差分量估計確定的權(quán)比就易失真；二是在觀測數(shù)據(jù)預(yù)處理過程中，如果各類觀測值殘差相差太大，迭代過程可能就出現(xiàn)不收斂或迭代次數(shù)太多，S矩陣就為奇異陣，無法定位解算。因此，解算過程中需注意歷元中各類衛(wèi)星數(shù)，以決定是否采用Helmert方差分量估計；在遇到迭代過程發(fā)散，需再對各衛(wèi)星殘差進行分析，重新篩選衛(wèi)星或適當(dāng)降低一些衛(wèi)星的權(quán)重，如抗差，再進行解算。

2.2 基于IGGIII的抗差Helmert方差分量估計

由最小二乘原理知，當(dāng)觀測值中含有粗差時，觀測值的殘差易受粗差的影響，從而導(dǎo)致Helmert方差分量估計出現(xiàn)震蕩現(xiàn)象，迭代結(jié)果失真[15]。為了避免粗差對Helmert方差分量估計結(jié)果的影響，可以采用等價權(quán)的方法對含有粗差的數(shù)據(jù)降權(quán)處理[2]。在BDS、GPS、GLONASS組合定位中，系統(tǒng)間觀測值相互獨立，相應(yīng)的等價權(quán)為

(6)

式中:Pi為原始權(quán)陣，wi為自適應(yīng)權(quán)因子，采用IGGⅢ函數(shù)[16]，如下式：

(7)

式中：k1，k2為常數(shù)，一般取k1=1.5～2.0，k2=3.0～8.0。結(jié)合Helmert方差分量估計，抗差Helmert方差分量估計的主要步驟是在步驟3)后，引入等價權(quán)因子重新定權(quán)，再進行方差分量估計。

3 實驗分析

實驗數(shù)據(jù)采用MGEX nnor站2015年12月12日24 h的觀測數(shù)據(jù)，采樣率為30 s，共計2 880個歷元。定位結(jié)果精度分析采用將各歷元解算的坐標(biāo)與IGS發(fā)布的測站參考真值間的偏差轉(zhuǎn)換為站心之間坐標(biāo)E、N、U 3個方向，分別統(tǒng)計3個方向分量的內(nèi)外符合精度STD和RMS。為了對比分析抗差Helmert方差分量估計、Helmert方差分量估計以及常規(guī)高度角定權(quán)3種權(quán)函數(shù)模型對BDS/GPS/GLONASS組合系統(tǒng)定位精度的影響和抵御粗差的能力以及組合系統(tǒng)的導(dǎo)航定位精度，進行以下兩種方案實驗：

方案1：截止高度角為10°、15°和20°時，3種不同權(quán)函數(shù)下，組合系統(tǒng)的導(dǎo)航定位精度分析；

方案2：截止高度角為30°和40°時，組合系統(tǒng)和單系統(tǒng)BDS的定位精度分析。

實驗方案1中，衛(wèi)星截止高度角為10°的定位解算結(jié)果如圖1～圖7所示。圖1、圖2分別為GPS、BDS單系統(tǒng)定位結(jié)算結(jié)果；圖3、圖4分別為GPS/BDS組合定位Helmert方差解和抗差解；圖5、圖6、圖7分別為BDS/GPS/GLONASS組合定位3種權(quán)函數(shù)模型的定位結(jié)果。

在截止高度角為10°下，不同衛(wèi)星系統(tǒng)和權(quán)函數(shù)模型的定位結(jié)果統(tǒng)計如表1所示，統(tǒng)計歷元數(shù)都大于2 800，接近總歷元數(shù)。

從圖1～圖7和表1可知，GPS定位精度優(yōu)于BDS，兩系統(tǒng)定位結(jié)果在U方向表現(xiàn)一致，都存在周期性波動，GPS更顯著；而在E、N方向，BDS變化曲線穩(wěn)定性明顯差于GPS。組合系統(tǒng)定位精度相對GPS單系統(tǒng)提升較少，15%左右，而對BDS提升約40%，顯著提高。但在E方向上，組合系統(tǒng)定位精度低于GPS單系統(tǒng)，N方向上提升也較小，U方向上改善顯著，說明衛(wèi)星數(shù)增加對提高GPS單系統(tǒng)U方向精度作用最大。組合系統(tǒng)BDS/GPS/GLONASS相對BDS/GPS的定位結(jié)果，兩者定位精度相當(dāng)，可能在于觀測到的GLONASS衛(wèi)星數(shù)較少，BDS/GPS組合定位精度已較高，衛(wèi)星數(shù)再增加對精度提升作用較小。組合系統(tǒng)BDS/GPS或BDS/GPS/GLONASS，抗差Helmert方差分量估計定位解與Helmert方差分量估計定位解精度一樣，稍差于高度角權(quán)函數(shù)模型，可能因為測站觀測數(shù)據(jù)質(zhì)量較好，無明顯粗差，反而方差分量估計迭代過程中，不同系統(tǒng)間衛(wèi)星數(shù)差異和各衛(wèi)星后驗殘差相差較大，致確定的系統(tǒng)間權(quán)比出現(xiàn)大的偏差，因而迭代解算結(jié)果更差，如圖中的“尖刺”點。也說明當(dāng)多余觀測數(shù)不充足時或組合系統(tǒng)中某個系統(tǒng)可見星較少時，采用Helmert方差分量估計調(diào)整系統(tǒng)間權(quán)比不再適用。

圖1 GPS單系統(tǒng)定位結(jié)果在E、N、U方向的坐標(biāo)差值

圖2 BDS單系統(tǒng)定位結(jié)果在E、N、U方向的坐標(biāo)差值

圖3 BDS/GPS組合定位Helmert方差解在E、N、U方向的坐標(biāo)差值

圖4 BDS/GPS組合定位抗差解在E、N、U方向的坐標(biāo)差值

圖5 BDS/GPS/GLONASS組合高度角定權(quán)函數(shù)模型定位解算結(jié)果

圖6 BDS/GPS/GLONASS組合Helmert方差分量估計權(quán)函數(shù)模型定位解算結(jié)果

圖7 BDS/GPS/GLONASS組合抗差分量估計權(quán)函數(shù)模型定位解算結(jié)果

表1 不同衛(wèi)星系統(tǒng)定位結(jié)果 m

方案1中，衛(wèi)星截止高度角為15°和20°時，不同模型定位結(jié)果統(tǒng)計如表2、表3所示。

表2 15°下不同模型定位結(jié)果統(tǒng)計 m

表3 20°下不同模型定位結(jié)果統(tǒng)計 m

從表1～表3的統(tǒng)計結(jié)果可知，隨衛(wèi)星截止高度角增大，導(dǎo)航系統(tǒng)定位精度逐漸降低。對GPS單系統(tǒng)，定位精度下降趨勢明顯，主要是衛(wèi)星數(shù)減少。而BDS，10°和15°情況下，系統(tǒng)定位精度及N、E、U 3方向上差異很小，20°時定位精度下降較低，體現(xiàn)了BDS異構(gòu)星座的特點。兩系統(tǒng)都呈現(xiàn)出在E、N方向上，定位精度受衛(wèi)星數(shù)變化的影響小于U方向。同樣，對于3種權(quán)函數(shù)模型，高度角權(quán)函數(shù)模型定位精度最高，Helmert方差分量估計和抗差分量估計對提高定位精度的作用不是很顯著。組合系統(tǒng)BDS/GPS和BDS/GPS/GLONASS定位精度都優(yōu)于3 m，精度相當(dāng)；相對單系統(tǒng)，精度都有顯著提高。

針對上述試驗中3種權(quán)函數(shù)定位結(jié)果精度相當(dāng)，Helmert方差分量估計和抗差分量估計對提升定位精度的作用并不明顯的問題，分析原因可能主要在于實驗中采用的nnor站觀測數(shù)據(jù)質(zhì)量較好。而在實際觀測過程中，因為觀測條件和儀器的差異，GNSS觀測數(shù)據(jù)質(zhì)量可能會稍差，含部分粗差。因此，為進一步分析高度角權(quán)函數(shù)、Helmert方差分量估計和抗差分量估計抵御粗差的能力，在原始觀測數(shù)據(jù)中隨機加入10～20 m不等的偏差，3種權(quán)函數(shù)的定位結(jié)果如圖8～圖10所示。

圖8 含粗差下BDS/GPS/GLONASS組合高度角定權(quán)函數(shù)模型定位解算結(jié)果

圖9 含粗差下BDS/GPS/GLONASS組合Helmert方差分量估計權(quán)函數(shù)模型定位解算結(jié)果

圖10 含粗差下BDS/GPS/GLONASS組合抗差分量估計權(quán)函數(shù)模型定位解算結(jié)果

從圖8～圖10中明顯可以看出，在觀測值含粗差時，高度角權(quán)函數(shù)定位結(jié)果受粗差影響較大，發(fā)生了大的偏差；Helmert方差分量估計定位結(jié)果也同樣受到粗差的影響，定位解算過程中方差分量估計迭代失真和迭代過程不收斂，無法定位，即圖9中“尖刺”點和折線連接部分；而抗差分量估計能較好抵御粗差的影響，改善粗差對定位結(jié)果的“污染”，提高定位精度。

針對組合導(dǎo)航系統(tǒng)衛(wèi)星數(shù)成倍增加和更優(yōu)的衛(wèi)星空間幾何構(gòu)型的特點，進一步探討分析其在截止高度角較大的情況下的導(dǎo)航定位性能，即實驗方案2。圖11為組合系統(tǒng)和單系統(tǒng)在衛(wèi)星截止高度角為30°時的可視衛(wèi)星數(shù)，圖12～圖16分別為GPS、BDS、GLONASS 3種系統(tǒng)不同組合模式在高度角權(quán)函數(shù)下，各系統(tǒng)的定位結(jié)果在E、N、U 3方向的坐標(biāo)差值。30°不同模型定位結(jié)果統(tǒng)計見表4。

從圖11可知，在高度截止角為30°時，BDS單系統(tǒng)整天內(nèi)可見衛(wèi)星數(shù)都大于5顆，為6～10顆，相對GPS、GLONASS可觀測時間更長，說明BDS星座分布更具優(yōu)勢。組合系統(tǒng)相對任一單系統(tǒng)可見衛(wèi)星數(shù)更多，G/R/C可見衛(wèi)星數(shù)最多，絕大部分觀測時段可見衛(wèi)星數(shù)大于15顆；G/C可見衛(wèi)星數(shù)在13顆左右，C/R為11顆左右；G/R最少，可見衛(wèi)星數(shù)與觀測時段有關(guān)，變化較大。

圖11 導(dǎo)航系統(tǒng)的可見衛(wèi)星數(shù)

圖12 BDS定位結(jié)果

圖13 BDS/GPS組合系統(tǒng)定位結(jié)果

圖14 GPS/GLONASS組合系統(tǒng)定位結(jié)果

圖15 BDS/GLONASS組合系統(tǒng)定位結(jié)果

圖16 BDS/GPS/GLONASS組合系統(tǒng)定位結(jié)果

表4 30°不同模型定位結(jié)果統(tǒng)計 m

從圖12～圖16和表4可以看出，在衛(wèi)星高度截止角較大的情況下，組合系統(tǒng)相對單系統(tǒng)的導(dǎo)航定位優(yōu)勢更明顯。單系統(tǒng)GPS和GLONASS由于可見衛(wèi)星數(shù)不足或可觀測時間較短，不能提供導(dǎo)航定位服務(wù)，而組合系統(tǒng)全天都能用于導(dǎo)航定位。特別，BDS單系統(tǒng)整天可視衛(wèi)星數(shù)也能滿足定位，定位精度也相對較高，為4.459 m，與BDS/GLONASS精度相當(dāng)。組合系統(tǒng)BDS/GPS/GLONASS定位精度最高，點位精度為3.098 m；BDS/GPS組合定位精度優(yōu)于GPS/GLONASS和BDS/GLONASS，為3.220 m；BDS/GLONASS組合定位精度為4.311 m，與BDS單系統(tǒng)相差無幾；GPS/GLONASS組合系統(tǒng)由于衛(wèi)星數(shù)較少，且不穩(wěn)定，定位精度最差。

隨著衛(wèi)星截止高度角升高，組合系統(tǒng)可見星也逐漸減少。在截止高度角為40°時，單系統(tǒng)BDS和GPS/GLONASS、BDS/GLONASS組合系統(tǒng)大部分時段可觀測衛(wèi)星數(shù)較少，不能定位或精度太差；而BDS/GPS和BDS/GPS/GLONASS組合系統(tǒng)整天可見星數(shù)依舊大于6顆，為7～14顆，其定位結(jié)果如圖17、圖18所示。40°不同模型定位結(jié)果統(tǒng)計見表5。

從圖中可以看出，組合系統(tǒng)定位結(jié)果在E、N方向上較穩(wěn)定，平面位置精度仍相對較高，在3 m以內(nèi)；而U方向上波動較大，絕大部分觀測時段在-10～+10 m內(nèi)波動，個別時刻達到20 m左右；BDS/GPS/GLONASS組合系統(tǒng)相對BDS/GPS系統(tǒng)，增加的衛(wèi)星數(shù)對改善U方向的精度作用很小，甚至可能降低，如圖中10～11點間所示。

圖17 BDS/GPS組合系統(tǒng)定位結(jié)果

圖18 BDS/GPS/GLONASS組合定位結(jié)果

表5 40°不同模型定位結(jié)果統(tǒng)計 m

從表5可知，在高度截止角為40°時，BDS/GPS組合定位精度為5.821 m，平面定位精度優(yōu)于3 m，高程精度稍差，為5.474 m；BDS/GPS/GLONASS組合系統(tǒng)定位精度為5.864 m，平面定位精度與BDS/GPS一樣，高程精度略差于BDS/GPS，兩者導(dǎo)航精度都小于6 m。這也充分說明組合系統(tǒng)的導(dǎo)航定位優(yōu)勢，在高截止衛(wèi)星高度角下，如城市、山區(qū)和峽谷地區(qū)等，可見星數(shù)仍足夠定位，且導(dǎo)航定位精度較高，也能連續(xù)提供服務(wù)。

4 結(jié) 論

本文針對BDS/GPS/GLONASS組合系統(tǒng)聯(lián)合定位時，不同權(quán)函數(shù)模型——高度角定權(quán)、Helmert方差分量估計和抗差分量估計的定位精度和組合導(dǎo)航系統(tǒng)及單系統(tǒng)BDS導(dǎo)航定位性能進行分析，并探討了3種權(quán)函數(shù)對粗差的抵御能力。實驗結(jié)果表明，當(dāng)觀測值含有粗差時，抗差分量估計定位結(jié)果最好。

在低截止衛(wèi)星高度角下，BDS/GPS/GLONASS組合系統(tǒng)定位精度與BDS/GPS相當(dāng)，較GPS單系統(tǒng)定位精度提升15%，BDS提升40%，隨衛(wèi)星截止高度角增大，GPS精度提升增大到23%，而BDS不變；特別，當(dāng)衛(wèi)星截止高度角為30°時，GPS單系統(tǒng)可見衛(wèi)星數(shù)低于4顆，不能定位；而BDS系統(tǒng)由于其異構(gòu)星座分布的特點，可視衛(wèi)星數(shù)為6～9顆，其導(dǎo)航定位精度優(yōu)于5 m；且隨著我國北斗的建設(shè)，衛(wèi)星數(shù)將不斷增多，其城市高截止高度角地區(qū)的應(yīng)用將不再是問題。

BDS/GPS/GLONASS組合系統(tǒng)在城市或峽谷地區(qū)的應(yīng)用，即當(dāng)衛(wèi)星截止高度角較大時，等于30°或40°，組合系統(tǒng)相對單系統(tǒng)的導(dǎo)航定位優(yōu)勢更為突出，單系統(tǒng)往往不能定位或?qū)Ш蕉ㄎ环?wù)不連續(xù)，而組合系統(tǒng)能提供連續(xù)的導(dǎo)航定位，且精度優(yōu)于6 m。

上述結(jié)論僅只是相對于中國區(qū)域或周邊部分地區(qū)。

[1] GRAFAREND E.W,Schaffrin B.Variance-Covariance-Component Estimation of Helmet Type[J],Surveying and Mapping,1979.No.1.

[2] CAI C,GAO L P Y.A Precise Weighting Approach with Application to Combined L1/B1 GPS/BeiDou Positioning[J].Journal of Navigation,2014,67(5):911-925.

[3] 劉長建，馬高峰.抗差Helmert方差分量估計及其應(yīng)用[J].北京測繪，2002(1)：16-18.

[4] 秦顯平，楊元喜.抗差方差分量估計在衛(wèi)星定軌中的應(yīng)用[J].大地測量與地球動力學(xué)，2003,23(4)：40-43.

[5] 高曉,戴吾蛟.抗差Helmert方差分量估計在GPS/BDS組合定位中的應(yīng)用[J].大地測量與地球動力學(xué),2014,34(1):173-176.

[6] 張勇,田林亞,徐君民,等.抗差方差分量估計在GPS/GLONASS組合導(dǎo)航中的應(yīng)用[J].測繪科學(xué)技術(shù)學(xué)報,2013,30(2):132-135.

[7] 李鶴峰,黨亞民,秘金鐘,等.BDS與GPS、GLONASS多模融合導(dǎo)航定位時空統(tǒng)一[J].大地測量與地球動力學(xué),2013,33(4):73-78.

[8] 何俊,袁小玲,曾琪,等.GPS/BDS/GLONASS組合單點定位研究[J].測繪科學(xué),2014,39(8):124-128.

[9] 徐龍威,劉暉,張明,等.不同截止高度角多模GNSS組合單點定位性能分析[J].大地測量與地球動力學(xué),2015(6):987-991.

[10] SANTERRE R,PAN L,CAI C,et al.Single Point Positioning Using GPS,GLONASS and BeiDou Satellites[J].Positioning,2014,5(4):107-114.

[11] MAN X,SUN F,LIU S,et al.Analysis of Positioning Performance on Combined BDS/GPS/GLONASS[M].China Satellite Navigation Conference (CSNC) 2015 Proceedings: Volume II.Springer Berlin Heidelberg,2015:399-409.

[12] 程鵬飛,文漢江,成英燕,等.2000國家大地坐標(biāo)系橢球參數(shù)與GRS 80和WGS 84的比較[J].測繪學(xué)報,2009,38(3):189-194.

[13] 郭際明,孟祥廣,李宗華,等.GLONASS衛(wèi)星廣播星歷精度分析[J].大地測量與地球動力學(xué),2011,31(1):68-71.

[14] 劉長建,馬高峰.Helmert方差分量估計的粗差檢驗與抗差解[J].測繪信息與工程,2002,27(6):5-7.

[15] 黃維彬.近代平差理論及其應(yīng)用[M].北京：解放軍出版社，1992.

[16] 楊元喜.自適應(yīng)動態(tài)導(dǎo)航定位 [M].北京：測繪出版社，2006.

[責(zé)任編輯：劉文霞]

Combined single point positioning with different cut-off elevation angle using variance component estimation

FENG Biao, CHAI Hongzhou, WAN Min, PAN Zongpeng, LI Di

(Information Engineering University, Zhengzhou 450000, China)

When the combination of BDS/GPS/GLONASS system is positioned, due to the difference of the accuracy of satellite ranging between systems, a reasonable determination of weight ratio between satellites is needed.Helmert variance component estimation is often used to determine weight ratio of different observations. However, when the observations contain gross errors, the positioning results of Helmert variance component estimation are easy to be contaminated by gross errors or convergence distortion. This paper, based on Helmert variance component estimation, introduces the equivalent weight factor IGGIII function, to establish weight function model of robust Helmert variance component estimation, comparing and analyzing its application and positioning accuracy in combined of BDS/GPS/GLONASS system positioning under low cut-off angles of 10°, 15° and 20°. Then it discusses and analyzes the combination system and BDS positioning accuracy with the high cut-off elevation angle of 30°and 40°. Experimental results show that: when observations have no significant gross errors, the positioning accuracy of Helmert variance component estimation and robust component is quite similar, and slightly lower than that of elevation angle weigh function, of which the positioning accuracy is better than 2.5 m; while containing outliers, the positioning accuracy of robust solution is the highest; when the cut-off elevation angle is 30°,the positioning accuracy of the BDS single system is better than 5 m, while that of the combination system is nearly 3 m; when is 40°, its horizontal accuracy is better than 2 m, and 3D accuracy is better than 6 m, while the single system can’t be positioned.

combined BDS/GPS/GLONASS positioning; cut-off elevation angle; Helmert variance component estimation; IGIII equivalent weight function; robust estimation; positioning accuracy

引用著錄：馮彪，柴洪洲，王敏，等.基于方差分量估計的不同截止高度角下的組合單點定位[J].測繪工程，2017，26(4)：32-42.

10.19349/j.cnki.issn1006-7949.2017.04.007

2016-07-20

國家自然科學(xué)基金資助項目(41574010；41274045)

馮 彪(1990-)，男，碩士研究生.

P228.4

A

1006-7949(2017)04-0032-11