高 曉 戴吾蛟

1)中南大學(xué)地球科學(xué)與信息物理學(xué)院，長沙 410083

2)湖南省精密工程測量與形變?yōu)暮ΡO(jiān)測重點(diǎn)實(shí)驗(yàn)室，長沙410083

1 引言

由于BDS 與GPS 均采用碼分多址通信模式，GPS/BDS 組合定位優(yōu)勢更為明顯。GPS 與BDS是兩個(gè)獨(dú)立的衛(wèi)星系統(tǒng)，組合定位時(shí)須合理確定兩類觀測值權(quán)比，即確定觀測值的方差-協(xié)方差陣，其數(shù)值取決于觀測值隨機(jī)噪聲水平［1］。隨著衛(wèi)星不斷運(yùn)動，隨機(jī)噪聲發(fā)生改變，采用事先給定的固定權(quán)比進(jìn)行組合定位顯然是不合理的。Helmert 方差分量估計(jì)是一種通過迭代計(jì)算自適應(yīng)確定不同類觀測值權(quán)比的驗(yàn)后方差分量估計(jì)方法［2］，在數(shù)據(jù)處理領(lǐng)域得到了廣泛的應(yīng)用，并取得了良好的效果［3－5］。基于最小二乘原理的驗(yàn)后方差分量估計(jì)模型對粗差觀測值很敏感，當(dāng)觀測值包含粗差時(shí)，粗差會“污染”隨機(jī)模型，導(dǎo)致迭代收斂失真［6］，定位結(jié)果出現(xiàn)較大偏差?？共頗elmert 方差分量估計(jì)基于相關(guān)等價(jià)權(quán)原理確定自適應(yīng)降權(quán)因子，通過降低粗差觀測值的權(quán)重削弱其對定位結(jié)果的影響，提高定位精度。本文擬將該方法應(yīng)用于GPS/BDS 組合偽距定位以及載波相位相對定位中以檢驗(yàn)該方法的有效性與可行性。

2 抗差Helmert 方差分量估計(jì)模型

2.1 Helmert 方差分量估計(jì)模型

嚴(yán)密的Helmert 方差分量估計(jì)公式為［3］:

其中，n1、n2分別為兩類觀測值個(gè)數(shù)，N 為觀測值設(shè)計(jì)矩陣，V 為觀測值殘差矩陣。解得與之后，重新調(diào)整兩類觀測值的權(quán):

式中，C 為常數(shù)。

2.2 基于IGGIII 模型的抗差Helmert 方差分量估計(jì)模型

楊元喜［7］根據(jù)穩(wěn)健估計(jì)思想構(gòu)建觀測向量抗差等價(jià)權(quán)，其中IGGⅢ模型為:

式中，k1、k2為常數(shù)，參照文獻(xiàn)［8］，k1=1.5 ～2.0，k2=3.0 ～8.5。δ 為方差因子，其取值為:

由于標(biāo)準(zhǔn)化殘差向量v 可以較好地反映觀測向量中的異常信息，當(dāng)觀測值帶有粗差時(shí)，該異常信息將及時(shí)反映到殘差向量中，利用該模型可以降低粗差觀測值的影響，避免粗差污染定位結(jié)果，保證定位的可靠性。由于觀測值精度不同導(dǎo)致不同類觀測值組成的殘差向量也存在系統(tǒng)性差異，采用自適應(yīng)因子降權(quán)處理時(shí)無法準(zhǔn)確分辨粗差與不同類觀測值正確殘差信息。因此，應(yīng)先根據(jù)方差分量估計(jì)思想統(tǒng)一不同類觀測值的中誤差，確定不同類觀測值合理權(quán)比，再利用上述模型抵抗粗差影響，從而提高定位精度與可靠性。

2.3 抗差Helmert 方差分量估計(jì)算法流程

1)根據(jù)高度角模型進(jìn)行初始定權(quán)［9］，即

其中，e 為觀測衛(wèi)星高度角，a、b、σ 根據(jù)經(jīng)驗(yàn)設(shè)定;

2)利用初始權(quán)陣進(jìn)行首次平差計(jì)算，獲取觀測值殘差信息;

3)按照式(1)計(jì)算GPS 與BDS 兩類觀測值方差估值)，重新定權(quán);

4)按照式(3)計(jì)算自適應(yīng)抗差因子，再次定權(quán)，式(3)中方差因子σ 可以按照式(4)確定;

5)基于抗差Helmert 方差分量估計(jì)模型重新進(jìn)行平差計(jì)算，獲取新的觀測值殘差信息;

3 算例分析

3.1 實(shí)驗(yàn)數(shù)據(jù)

采用GPS/BDS 雙系統(tǒng)接收機(jī)于2012年11月28日在中南大學(xué)校本部采集4 小時(shí)數(shù)據(jù)，采樣間隔設(shè)定為1 秒，共計(jì)14 400 個(gè)觀測歷元。統(tǒng)計(jì)測段內(nèi)GPS 與BDS 觀測衛(wèi)星數(shù)如圖1?？梢?，GPS 可視衛(wèi)星數(shù)為6 ～10 顆。

圖1 GPS/BDS 觀測衛(wèi)星數(shù)Fig.1 Number of available satellites of GPS and BDS

圖2 為測段內(nèi)GPS/BDS 組合系統(tǒng)GDOP 值分時(shí)統(tǒng)計(jì)圖。由圖2 可知，當(dāng)天組合系統(tǒng)衛(wèi)星分布強(qiáng)度較高(GDOP 最小值為1.45，最大值為2.48，均值為1.82)，各歷元GDOP 值均小于閾值(4.0)，測段內(nèi)觀測數(shù)據(jù)無異常。

3.2 解算結(jié)果分析

高度角模型與Helmert 方差分量估計(jì)模型不足之處在于無法抵御粗差影響。抗差Helmert 方差分量估計(jì)由于加入抗差模型能夠較好地抵制粗差。為體現(xiàn)這一特點(diǎn)，隨機(jī)選取C01 衛(wèi)星8 個(gè)歷元的C1碼觀測值，人為添加三倍中誤差(30 m)的粗差，分別基于高度角模型、Helmert 方差分量估計(jì)模型以及抗差Helmert 方差分量估計(jì)模型進(jìn)行組合偽距單點(diǎn)定位計(jì)算，結(jié)果如表1 與圖3 所示。

從表1 和圖3 可以看出，添加粗差的8 個(gè)歷元基于不同模型定位結(jié)果差異顯著。高度角模型無法抵御粗差的影響，定位結(jié)果出現(xiàn)明顯偏差，水平方向與高程方向均呈現(xiàn)明顯的“尖刺”狀。當(dāng)觀測值包含粗差時(shí)，Helmert 方差分量模型明顯收斂失真，特別是東西方向定位偏差更加明顯?？共頗elmert 方差分量估計(jì)模型在三個(gè)方向均無明顯震蕩現(xiàn)象，真誤差曲線較之其他模型更加平緩，可見抗差Helmert方差分量估計(jì)模型不僅可以合理確定GPS 與BDS組合偽距觀測值權(quán)比，同時(shí)可以有效抵制粗差影響，提高定位精度與可靠性。

圖2 GPS/BDS 衛(wèi)星GDOP 值分時(shí)統(tǒng)計(jì)Fig.2 Statistics of time-share of GDOP for GPS and BDS

表1 不同定權(quán)方案偽距單點(diǎn)定位結(jié)果對比(單位:m)Tab.1 Comparison of results of pseudorange single point position with different stochastic models(unit:m)

圖3 不同定權(quán)方案組合偽距定位真誤差序列Fig.3 Compareation of real errors for different models in combined GPS and BDS pseudorange position

圖4 不同定權(quán)方案組合相對定位真誤差序列Fig.4 Compareation of real errors for different models in combined GPS and BDS pseudorange position

為檢驗(yàn)抗差Helmert 方差分量估計(jì)模型對于載波相位觀測值的有效性，隨機(jī)選取C01 衛(wèi)星7 個(gè)歷元L1觀測值，人為添加三倍中誤差(0.6 周)的粗差，基于不同隨機(jī)模型進(jìn)行GPS/BDS 組合動態(tài)相對定位計(jì)算，結(jié)果如圖4。同時(shí)，將添加粗差的7 個(gè)歷元定位結(jié)果進(jìn)行統(tǒng)計(jì)，結(jié)果如表2。

由圖4 和表2 可以看出，高度角模型無法抵御粗差影響，7 個(gè)添加粗差的歷元定位結(jié)果存在明顯偏差，東西方向更加明顯。Helmert 方差分量估計(jì)模型通過合理確定GPS 與BDS 兩類觀測值權(quán)比提高了定位精度，但是該模型存在收斂失真現(xiàn)象，特別是測段后期定位結(jié)果出現(xiàn)了明顯的震蕩現(xiàn)象，個(gè)別歷元定位結(jié)果偏差甚至大于高度角模型，可見粗差已經(jīng)“污染”隨機(jī)模型導(dǎo)致收斂失真，降低了定位可靠性。抗差Helmert 方差估計(jì)模型較好地抑制粗差影響，7 個(gè)帶有粗差觀測值的歷元定位結(jié)果與真值偏差較小，可見該模型有助于提高定位精度與可靠性。

4 結(jié)論

由于GPS 與BDS 存在系統(tǒng)性差異，且觀測值隨機(jī)噪聲水平隨衛(wèi)星運(yùn)動而發(fā)生改變，組合定位時(shí)依據(jù)固定權(quán)比進(jìn)行定位是不合理的。Helmert 方差分量估計(jì)方法可以合理確定GPS 與BDS 兩類觀測值權(quán)比，提高定位精度。但當(dāng)觀測值存在粗差時(shí)，方差分量估計(jì)方法會出現(xiàn)收斂失真現(xiàn)象，導(dǎo)致定位結(jié)果偏差較大。基于IGGIII 模型的抗差Helmert 方差分量估計(jì)模型采用等價(jià)權(quán)的方法對粗差觀測值進(jìn)行降權(quán)處理，削弱甚至消除粗差觀測值的影響，同時(shí)可以動態(tài)確定兩類觀測值權(quán)比，定位精度與可靠性均有提高。

1 何海波，楊元喜.GPS 觀測量先驗(yàn)方差-協(xié)方差矩陣實(shí)時(shí)估計(jì)［J］.測繪學(xué)報(bào)，2001，30(1):42－47.(He Haibo and Yang Yuanxi.Real-time estimation of a prior variance-covariance for GPS observations［J］.Acta Geodaetica et Cartographica Sinica，2001，30(1):42－47)

2 高曉，戴吾蛟.基于方差分量估計(jì)確定GPS/BDS 組合定位先驗(yàn)權(quán)比［J］.大地測量與地球動力學(xué)，2013，(2):136－138.(Gao Xiao and Dai Wujiao.Determination of prior weight ratio in combination of GPS-BDS position based on variance component estimation，2013，(2):136－138)

3 段舉舉，沈云中.基于方差分量估計(jì)的GPS/GLONASS 組合點(diǎn)定位［J］，測繪通報(bào)，2011，(4):4－6.(Duan Juju and Shen Yunzhong.GPS/GLONASS combined position using variance component estimation［J］.Bulletin of Surveying and Mapping，2011，(4):4－6)

4 范百興，李廣云.Helmert 方差分量估計(jì)在跨河水準(zhǔn)測量中的應(yīng)用［J］.測繪學(xué)院學(xué)報(bào)，2004(3)，21(1):27－29.(Fan Baixing and Li Guanyun.Application of Helmert-WF variance component estimation to river-crossing leveling［J］.Journal of Institute of Surveying and Mapping，2004(3)，21(1):27－29)

5 高俊強(qiáng)，胡燦.Helmert 方差分量估計(jì)在隧道貫通中控制導(dǎo)線平差的研究［J］.中國礦業(yè)大學(xué)學(xué)報(bào)，2006(1)，35(1):125－129.(Gao Junqiang and Hu Can.Study on application of Helmert variance estimation to tunnel runthrough control traverse difference［J］.Journal of China University of Mining ＆ Technology，2006(1)，35(1):125－129)

6 劉長建，馬高峰.Helmert 方差分量估計(jì)的粗差檢驗(yàn)與抗差解［J］.測繪信息與工程，2002，27(6):5－7.(Liu Changjian and Mu Gaofeng.Blunder test and robust solution of variance component estimation for Helmert type［J］.Journal of Geomatics，2002，27(6):5－7)

7 Yang Y X.Robust estimation of geodetic datum transformation［J］.Journal of Geodesy，1999，73:345－349.

8 楊元喜，宋力杰，徐天河.大地測量相關(guān)觀測抗差估計(jì)理論［J］.測繪學(xué)報(bào)，2002，31(2):95－99.(Yang Yuanxi，Song Lijie and Xu Tianhe.Robust parameter estimation for geodetic correlated observations［J］.Acta Geodaetica et Cartographica Sinica，2002，31(2):95－99)

9 戴吾蛟，丁曉利，朱建軍.基于觀測值質(zhì)量指標(biāo)的GPS 觀測量隨機(jī)模型分析［J］.武漢大學(xué)學(xué)報(bào)(信息科學(xué)版)，2008，33(7):718－722.(Dai Wujiao，Ding Xiaoli and Zhu Jianjun.Comparation of GPS stochastic models based on observation quality indices［J］.Geomatics and Information Science of Wuhan University，2008，33(7):718－722)