王亦欣

(湖南省第一測(cè)繪院，湖南衡陽(yáng) 421001)

1 概述

在GPS后處理中，經(jīng)過(guò)三維無(wú)約束平差(最小約束平差)后，一般要對(duì)基線殘差進(jìn)行分析以剔除粗差。無(wú)約束平差是以獨(dú)立基線向量及其方差－協(xié)方差陣為觀測(cè)量。平差后一般輸出基線向量各分量改正數(shù)(殘差)、基線邊長(zhǎng)、方位和點(diǎn)位精度信息等。常用的對(duì)基線殘差探測(cè)分析方法有X2檢驗(yàn)、τ檢驗(yàn)、環(huán)閉合差和重復(fù)基線分析等。但τ檢驗(yàn)對(duì)粗差探測(cè)往往具有不確定性，以下筆者試對(duì)其詳細(xì)分析。

2 τ檢驗(yàn)

觀測(cè)值誤差是服從正態(tài)分布的隨機(jī)變量，故殘差/標(biāo)準(zhǔn)殘差的期望是可預(yù)計(jì)的。τ檢驗(yàn)就是用標(biāo)準(zhǔn)殘差去統(tǒng)計(jì)測(cè)定觀測(cè)值殘差是否位于期望界限之內(nèi)，再計(jì)算出一個(gè)閥值來(lái)檢驗(yàn)每個(gè)標(biāo)準(zhǔn)殘差。每個(gè)被檢驗(yàn)的標(biāo)準(zhǔn)差可能有兩種結(jié)果:

①τ檢驗(yàn)通過(guò)，標(biāo)準(zhǔn)殘差不大于閥值，通常表明觀測(cè)值中無(wú)粗差存在。

②τ檢驗(yàn)失敗，標(biāo)準(zhǔn)殘差大于閥值，應(yīng)進(jìn)一步查驗(yàn)粗差情況。

τ檢驗(yàn)的嚴(yán)密公式如下:

構(gòu)造統(tǒng)計(jì)量:

其中:

vi為第i個(gè)觀測(cè)值的殘差;^σ0為單位權(quán)中誤差估值;qi為第i個(gè)觀測(cè)值的協(xié)因數(shù);f為自由度，等于多余觀測(cè)值的個(gè)數(shù)。

那么有:

α為顯著性水平。因t(f)呈對(duì)稱分布，如圖1所示。

圖1 τ檢驗(yàn)置信區(qū)間圖

通常選擇顯著性水平α=0.05，置信水平為95%時(shí)，若:

(tα/2(f)可在t分布表中查得)即Vi大于其中誤差的2倍時(shí)，認(rèn)為與它相應(yīng)的觀測(cè)量含有粗差。

為了減少棄真概率，根據(jù)“3σ”原則，選擇置信水平接近99.7%時(shí)，即Vi大于其中誤差的3倍時(shí)，則認(rèn)為其觀測(cè)量含有粗差。

需要注意的是，殘差未通過(guò)統(tǒng)計(jì)檢驗(yàn)，并不意味著觀測(cè)值中必然存在粗差，只是應(yīng)進(jìn)一步查驗(yàn)粗差情況，對(duì)可能存在粗差的基線作出標(biāo)記。一個(gè)觀測(cè)值的粗差往往會(huì)影響到其他觀測(cè)值殘差的推算。因此，這種標(biāo)記可能會(huì)張冠李戴，如果有一個(gè)甚至多個(gè)觀測(cè)值做上了這樣的標(biāo)記，可逐一檢測(cè)是否真正存在粗差。

3 實(shí)例分析

以中海達(dá)Hi－TargetTMGPS數(shù)據(jù)處理軟件算例進(jìn)行分析。該軟件實(shí)際提供的τ值為檢驗(yàn)值與τ值的比值，如果該值小于±1.0，則說(shuō)明該觀測(cè)值不含粗差，如它大于±1.0，則意味著該觀測(cè)值含有粗差。

例如:有GPS控制網(wǎng)如圖2所示。

經(jīng)三維無(wú)約束平差后進(jìn)行τ檢驗(yàn)，如表1所示。

表1和圖3表明有3條基線含有粗差。而這3條基線向量各分量的改正數(shù)如表2所示。

τ檢驗(yàn)列表 表1

圖3中，背景曲線為理論上的τ值分布曲線，直方圖為實(shí)際曲線。

圖3 τ檢驗(yàn)直方圖

基線分量改正數(shù) 表2

對(duì)基線向量各分量的改正數(shù)可以利用以下簡(jiǎn)單公式進(jìn)行檢驗(yàn)(見《全球定位系統(tǒng)(GPS)測(cè)量規(guī)范》GB/T18314－2001)。

式中a為固定誤差，mm;b為比例誤差系數(shù)，ppm;d為相鄰點(diǎn)間的距離，km;對(duì) GPS E級(jí)網(wǎng)，取 a=10 mm，b=10 ppm。

經(jīng)過(guò)計(jì)算，基線4001→4015.3165標(biāo)準(zhǔn)差 σ為4.40 cm;基線 4001→ZBXX.2971標(biāo)準(zhǔn)差 σ為4.97 cm;基線4002→4015.3164標(biāo)準(zhǔn)差 σ為3.18 cm。顯而易見，上述3條基線向量各分量的改正數(shù)都遠(yuǎn)遠(yuǎn)小于各自的3σ，并不存在粗差。

接下來(lái)，我們刪除4004號(hào)站點(diǎn)(該點(diǎn)離這3條基線并不近)，變成如圖4所示網(wǎng)形。

重新平差，再進(jìn)行τ檢驗(yàn)，如表3、圖5所示。

圖4 改變后的GPS控制網(wǎng)圖

τ檢驗(yàn)列表 表3

圖5 τ檢驗(yàn)直方圖

結(jié)果顯示，基線4001→4015.3165還含有粗差，而另外 2條基線 4001→ZBXX.2971和基線 4002→4015.3164則不再含有粗差了。觀察這3條基線向量各分量的改正數(shù)與前次平差的改正數(shù)相比并沒(méi)有多大的改動(dòng)。如表4所示。

基線分量改正數(shù) 表4

4 結(jié)語(yǔ)

從以上的算例分析可以看出，τ檢驗(yàn)對(duì)粗差的探測(cè)具有不確定性。如果殘差τ檢驗(yàn)通過(guò)則說(shuō)明觀測(cè)值不包含粗差。如果τ檢驗(yàn)失敗，粗差不一定存在，需要進(jìn)一步探查。可采用最有效的粗差探測(cè)方法:獨(dú)立環(huán)閉合差檢驗(yàn)。也可采用式(6)檢驗(yàn)，如果 V△X、V△Y、V△Z大于3σ，那么可以肯定該基線或該基線附近存在粗差，需要剔除粗差基線。因此，在 GPS網(wǎng)數(shù)據(jù)處理中，對(duì)平差結(jié)果的分析，過(guò)于看重τ檢驗(yàn)的結(jié)果是沒(méi)有必要的。

［1］李征航，黃勁松.GPS測(cè)量與數(shù)據(jù)處理［M］.武漢:武漢大學(xué)出版社，2005

［2］王國(guó)祥.GPS基線解算與質(zhì)量控制［J］.鐵道勘察，2005(6):7～9

［3］陶本藻，姚宜斌.可靠性分析與數(shù)據(jù)探測(cè)［J］.武漢大學(xué)學(xué)報(bào)(信息科學(xué)版)，2002(6):607

［4］GB/T18314－2009.全球定位系統(tǒng)(GPS)測(cè)量規(guī)范［S］.

［5］盛驟，謝式千，潘承毅.概率論與數(shù)理統(tǒng)計(jì)［M］.北京:高等教育出版社，2008