吝志勇，藍(lán)悅明

(1.保定市城市設(shè)計院，河北 保定071000;2.武漢大學(xué) 測繪學(xué)院，湖北 武漢430079)

如何在盡可能短的時間內(nèi)，獲得盡可能高精度的結(jié)果，是GNSS 應(yīng)用研究領(lǐng)域的熱點(diǎn)問題。目前的研究多是基于大樣本數(shù)據(jù)考慮的，對小樣本數(shù)據(jù)的情況考慮較少，而灰色系統(tǒng)理論的特色就是研究“小樣本”、“貧信息”的不確定性問題?；疑到y(tǒng)理論認(rèn)為，盡管客觀系統(tǒng)的表象復(fù)雜，數(shù)據(jù)凌亂，但它們總有自身的整體功能，必然蘊(yùn)藏某種內(nèi)在規(guī)律，關(guān)鍵是如何選擇適當(dāng)?shù)姆椒ㄈネ诰蚝屠盟?］。

1 灰色理論序列算子構(gòu)造與數(shù)據(jù)重生成

1.1 序列算子

對于GNSS 測量而言，由于觀測過程中受到諸如與衛(wèi)星有關(guān)的誤差、信號傳播誤差、與接收機(jī)有關(guān)的誤差，以及其他誤差等的影響，相應(yīng)的觀測值序列是一個沖擊擾動序列，要從沖擊擾動序列出發(fā)實(shí)現(xiàn)對真實(shí)行為序列系統(tǒng)變化規(guī)律的正確把握和認(rèn)識，必須首先設(shè)法排除系統(tǒng)行為數(shù)據(jù)所受到的沖擊波干擾，還數(shù)據(jù)以本來面目。設(shè)：

式(1)為系統(tǒng)真實(shí)行為序列，而觀測到的系統(tǒng)行為數(shù)據(jù)序列為：

式中：ε=(ε1，ε2，…，εn)為沖擊擾動項，X為沖擊擾動序列。

設(shè)D為作用于X的算子，X經(jīng)算子D作用后所得序列記為：

其中，XD為一階算子作用序列。

序列算子的作用可以進(jìn)行多次，相應(yīng)地，有二階、三階算子作用序列。

1.2 實(shí)用緩沖算子的構(gòu)造

設(shè)原始觀測值序列為：

則當(dāng)X為單調(diào)增長序列、單調(diào)衰減序列或震蕩序列時，D皆為弱化算子。

2 方程病態(tài)性及其改善

2.1 方程病態(tài)性討論

利用少數(shù)幾個歷元的觀測數(shù)據(jù)進(jìn)行基線解算時，由誤差方程組成的法方程是一個病態(tài)方程［2-3］。由于法矩陣條件數(shù)較大，法矩陣求逆不穩(wěn)定，觀測值微小的變化會導(dǎo)致解不穩(wěn)定。

Gauss-Markov 模型為：

式中：A為系數(shù)矩陣;X為待求參數(shù)向量;L為觀測值;Δ 是觀測噪聲;P為相應(yīng)觀測值的權(quán)陣。相應(yīng)參數(shù)X的最小二乘解為：

其中：N=ATPA;W=ATPL。當(dāng)法矩陣非奇異時，參數(shù)解存在且唯一;假設(shè)N、W有微小的擾動δN、δW，引起估計量有偏差δ，可表示為：

顯然，δN、δW對δ的影響與N-1密切相關(guān)。如果N的條件數(shù)比較大，當(dāng)N、W有一個微小的變化時，參數(shù)解相應(yīng)將產(chǎn)生較大的變化。

2.2 方程病態(tài)性改善

有偏估計中的嶺估計是克服方程病態(tài)性的常用方法，而嶺參數(shù)的選擇又是其中的關(guān)鍵和難點(diǎn)。根據(jù)TIKHONOV 正規(guī)化原理，對應(yīng)式(5)嶺估計的估計準(zhǔn)則為：

相應(yīng)嶺估計的解為：

考慮到少數(shù)幾個歷元數(shù)據(jù)組成的法矩陣特征值較大，選取正規(guī)化矩陣R代替式(9)中的單位矩陣I，相應(yīng)的方法稱為修正的TIKHONOV 正規(guī)化方法，估計準(zhǔn)則為：)

經(jīng)大量實(shí)際算例研究證明，用L 曲線法選取正規(guī)化參數(shù)α =1 時效果最好，正規(guī)化矩陣R的選取方式有很多［4-7］，筆者采用如下的形式：

式中：j為組成法方程的歷元數(shù);I為3 ×3 的單位矩陣;β 為權(quán)比系數(shù)，可以根據(jù)待定點(diǎn)近似坐標(biāo)的偏差和相位觀測值的精度比來確定，經(jīng)多個實(shí)例驗(yàn)證，其取值范圍為10-6～10-4。于是，正規(guī)化后的參數(shù)解為：

與最小二乘中的法矩陣相比，式(12)中的法矩陣中增加了R項。該措施的實(shí)質(zhì)是在基線分量對應(yīng)的部分附加微小的約束，并讓正規(guī)化矩陣隨著歷元數(shù)的變化而變化，正是由于R的參與，法矩陣的病態(tài)性得到抑制。

3 實(shí)例分析

灰色理論緩沖算子的引入可以減弱觀測數(shù)據(jù)中某些不確定性誤差的影響，改善觀測值的質(zhì)量。從定性分析的角度來看，利用緩沖處理后的數(shù)據(jù)進(jìn)行參數(shù)解算應(yīng)能取得較好的結(jié)果。為了體現(xiàn)小樣本數(shù)據(jù)的特點(diǎn)，實(shí)驗(yàn)中均選取了5 個歷元的觀測值。觀測值選取時，盡量選在衛(wèi)星觀測信號連續(xù)且穩(wěn)定的時間段內(nèi)，以避免粗差或周跳等帶來的影響。實(shí)驗(yàn)中，首先利用軟件TGO 對整個時段(1.5 h)內(nèi)的觀測值進(jìn)行平差計算，作為待求基線解的真值。實(shí)驗(yàn)中筆者設(shè)計了3 種不同的方案：

(1)利用原始觀測數(shù)據(jù)組成同步觀測值，建立單差、雙差模型，在雙差模式下求解模糊度浮動解，再根據(jù)模糊度浮動解及方差協(xié)方差陣，進(jìn)行模糊度搜索與固定，最終求得基線參數(shù)的固定解。

(2)利用原始觀測值φi進(jìn)行序列算子構(gòu)造并對原始觀測值序列φ 進(jìn)行緩沖處理，得到新的觀測值序列φ'(虛擬觀測值序列)后再建立雙差模型進(jìn)行參數(shù)求解計算，參數(shù)解算同實(shí)驗(yàn)方案(1)。

(3)利用兩個不同測站的原始觀測值φi、φj，組成雙差觀測值φ▽Δ，然后用雙差觀測值來構(gòu)造緩沖序列算子，進(jìn)行雙差觀測值序列的重生成，最后用新生成的雙差觀測值φ'▽Δ進(jìn)行參數(shù)的求解計算。3 種實(shí)驗(yàn)方案解算過程如圖1 所示。

圖1 基線解算處理流程圖

為了驗(yàn)證上述3 種方案的解算效果，筆者選取了3 條不同長度的基線觀測數(shù)據(jù)進(jìn)行解算，解算結(jié)果如表1 所示。

表1 基線解算結(jié)果對比表

從表1 可以看出：對于第一條基線來說，方案(2)、方案(3)的解算效果要優(yōu)于方案(1)，利用緩沖后的數(shù)據(jù)求得的結(jié)果更接近真值;對于后兩條基線來說，方案(2)的效果最好，解算的結(jié)果更接近真值。

4 結(jié)論

對于給定不同長度的基線觀測數(shù)據(jù)，尤其是短基線觀測數(shù)據(jù)來說，利用灰色理論緩沖算子對少數(shù)幾個歷元的觀測數(shù)據(jù)進(jìn)行數(shù)據(jù)重生成，然后再進(jìn)行參數(shù)求解計算，其結(jié)果比數(shù)據(jù)直接解算的效果要好，說明灰色理論在短基線解算中具有一定的優(yōu)越性。

［1］ 鄧聚龍.灰色理論基礎(chǔ)［M］.武漢：華中科技大學(xué)出版社，2002：153 -192.

［2］ 郭建峰. 測量平差系統(tǒng)病態(tài)性的診斷與處理［D］.鄭州：中國人民解放軍信息工程大學(xué)，2002.

［3］ 林勝良.病態(tài)線性方程組解法研究［D］.杭州：浙江大學(xué)，2005.

［4］ 崔希璋，於宗儔，陶本藻，等.廣義測量平差［M］.武漢：武漢測繪科技大學(xué)出版社，2001：163 -170.

［5］ 王振杰.大地測量中不適定問題的正規(guī)化解法研究［D］. 武漢：中國科學(xué)院測量與地球物理研究所，2003.

［6］ 陽仁貴，歐吉坤.提高GPS 單頻相位模糊度解算成功率的新方法［J］. 南京航空航天大學(xué)學(xué)報，2005，37(3)：279 -283.

［7］ 王振杰，歐吉坤，田愛奎.改進(jìn)的ARCE 方法及其在單頻GPS 快速定位中的應(yīng)用［J］. 武漢大學(xué)學(xué)報(信息科學(xué)版)，2005，30(6)：488 -491.