周 巍 郝金明 賈小林

(1)解放軍信息工程大學(xué)測繪學(xué)院，鄭州 450052 2)中國人民解放軍63883部隊，洛陽 471003 3)西安測繪研究所，西安710054)

不同數(shù)據(jù)預(yù)處理方法對北斗衛(wèi)星精密定軌精度的影響*

周 巍1，2)郝金明1)賈小林3)

(1)解放軍信息工程大學(xué)測繪學(xué)院，鄭州 450052 2)中國人民解放軍63883部隊，洛陽 471003 3)西安測繪研究所，西安710054)

提出綜合利用高次差法和多項式擬合法對COMPASS-MEO衛(wèi)星進(jìn)行周跳探測的原理及實現(xiàn)方法，并采用COMPASS監(jiān)測站實測數(shù)據(jù)使用不同的數(shù)據(jù)預(yù)處理策略對COMPASS-MEO衛(wèi)星進(jìn)行精密定軌試驗。結(jié)果表明，根據(jù)實測數(shù)據(jù)情況，合理選擇周跳探測方法和預(yù)處理策略，可以改善定軌精度。

COMPASS;周跳;φLC-PLC組合;多項式擬合;定軌

1 引言

衛(wèi)星導(dǎo)航系統(tǒng)將導(dǎo)航衛(wèi)星作為動態(tài)已知點，實現(xiàn)精密定軌是首先需要解決的問題，而觀測數(shù)據(jù)質(zhì)量是影響衛(wèi)星定軌精度的重要因素之一。對載波相位觀測數(shù)據(jù)而言，未被檢測出的周跳將被位置參數(shù)的改正數(shù)所吸收，從而導(dǎo)致參數(shù)估計的有偏。長期以來學(xué)者們對GPS周跳探測與修復(fù)問題做了深入的研究，對COMPASS系統(tǒng)具有重要的借鑒作用。與GPS相比，COMPASS提供更多類型的觀測數(shù)據(jù)，對于各類數(shù)據(jù)以前并沒有使用的經(jīng)驗，傳統(tǒng)GPS數(shù)據(jù)預(yù)處理中周跳探測方法是否適用于COMPASS，是否能夠改善定軌精度，都有待進(jìn)一步研究。

本文提出了綜合利用高次差法和多項式擬合法探測周跳的方法，并利用COMPASS監(jiān)測站實測數(shù)據(jù)，分析了不同周跳探測方法和數(shù)據(jù)預(yù)處理方案對COMPASS-MEO衛(wèi)星定軌結(jié)果的影響。

2 周跳探測原理

2.1 利用φLC－PLC組合探測周跳

周跳探測原理基于粗差理論，構(gòu)造能夠精確反映周跳變化的檢驗序列，使周跳在該檢驗序列中以粗差的形式體現(xiàn)，進(jìn)而確定周跳的位置和大小。

構(gòu)造雙頻消電離層組合為：

φLC-PLC組合消除了幾何距離、電離層、對流層、鐘差的影響，對該組合相鄰歷元間求差，若無周跳發(fā)生，殘差應(yīng)該在0值上下隨機波動，而周跳會以粗差的形式突顯于該檢驗序列中。由于φLC-PLC組合放大了觀測噪聲(約為P1噪聲的3倍［1］)，直接影響周跳探測的靈敏度，實驗發(fā)現(xiàn)該方法對10周以內(nèi)的周跳不敏感，且利用φLC-PLC組合探測周跳只能反映周跳位置，并不能分離出某一個頻率上的周跳幅值。

2.2 綜合利用高次差法和多項式擬合法探測周跳

多項式擬合法根據(jù)m個無周跳的載波相位觀測值擬合一個n階多項式，根據(jù)多項式預(yù)估下一個觀測值并與實測值比較，從而發(fā)現(xiàn)周跳并修正整周計數(shù)［2］。

其數(shù)學(xué)模型為：

如圖1和表1所示，對于觀測質(zhì)量較好的大采樣率數(shù)據(jù)而言，該方法對周跳的探測十分敏感，且周跳估值與預(yù)設(shè)值非常接近。

圖1 多項式擬合法探測周跳Fig.1 Cycle slip detection by polynomial fitting

表1 多項式擬合法周跳檢測結(jié)果(單位：周)Tab.1 Results of cycle slip detection by polynomial fitting (unit：cycle)

采用多項式擬合法探測周跳存在以下問題：每弧段啟動第一個窗口(擬合窗寬度m)擬合系數(shù)并對第m+1歷元進(jìn)行外推時，若前m個歷元就已發(fā)生周跳，會引起之后一系列誤報。而通過對大量實測數(shù)據(jù)分析發(fā)現(xiàn)，這種現(xiàn)象比較常見，是必須要解決的問題。

采用 Remondi提出利用差分法進(jìn)行周跳探測［4］。采用原始觀測數(shù)據(jù)或者其線性組合值歷元間求差，可以消除一些誤差，如接收機鐘差、衛(wèi)星鐘差、大氣折射延遲等，當(dāng)求至5次差時，如果沒有周跳發(fā)生，該高次差序列呈現(xiàn)偶然誤差特性，而周跳則會破壞這種規(guī)律性且幅值被放大，這類似一個高通濾波器，既阻尼低頻又消除其常數(shù)部分，高頻部分(如周跳)則被放大，根據(jù)該特性可以確定周跳發(fā)生的位置［5］。

聯(lián)合使用高次差法和多項式擬合法探測周跳，在多項式擬合法啟動第一個窗口內(nèi)，首先對該m個數(shù)據(jù)使用高次差法進(jìn)行探測，如果發(fā)生周跳，確定其發(fā)生時刻，從該歷元起重新開窗，直到有m個“干凈”的觀測數(shù)據(jù)起再啟動多項式擬合法探測周跳。

3 基于Hatch濾波的COMPASS雙頻載波相位平滑偽距

鑒于載波相位測量精度高，能夠精確反映偽距的變化，且多路徑誤差影響大約是測碼偽距觀測值的1/10等［6］原因，利用載波相位信息輔助碼偽距測量，以獲得比單獨采用碼偽距測量更高的精度［7］，基本公式為：

式(3)中包含歷元ti的電離層折射誤差相對于前i個歷元平均值的變化。當(dāng)歷元間的電離層折射發(fā)生變化時，偽距平滑值中包含2倍電離層延遲變化影響，當(dāng)電離層延遲發(fā)生突然抖動時，單頻偽距相位平滑的結(jié)果中會包含較大的偏差［8］。因此，單頻偽距進(jìn)行多歷元數(shù)據(jù)的平滑處理時，如果所采用弧段的長度不適當(dāng)，可能會發(fā)生單頻偽距相位平滑值的離散情況。

通過對大量實測數(shù)據(jù)的分析，本文選取的B1、B3頻率窄相關(guān)類型偽距和載波相位用于定軌運算。如式(1)構(gòu)造雙頻電離層組合，將式(3)轉(zhuǎn)化為：采用Hatch濾波方法［9］，雙頻相位平滑偽距逐歷元遞推計算公式表示為：

選取2008年9月5日COMPASS-MEO臨潼站觀測數(shù)據(jù)，如(5)式獲得B1、B3雙頻相位平滑偽距值，與B1頻率窄相關(guān)偽距對比，分析基于Hatch濾波的COMPASS-MEO雙頻相位平滑偽距值對精度的改善，圖2、圖3分別為偽距和載波相位平滑偽距擬合噪聲殘差圖，計算相應(yīng)RMS值統(tǒng)計如表2。

結(jié)果表明，B1頻率窄相關(guān)偽距觀測噪聲約為0.200 m，載波相位平滑偽距后為毫米級，大大降低了噪聲水平。故采用雙頻載波相位平滑偽距數(shù)據(jù)進(jìn)行定軌可以獲得更高的精度。

4 定軌結(jié)果分析

選取2008年3月4、5、6日，上海站、長春站、臨潼站、昆明站、烏魯木齊站、南極站觀測數(shù)據(jù)，對COMPASS-MEO衛(wèi)星進(jìn)行動力學(xué)法定軌，比較不同周跳探測方法、不同預(yù)處理策略對定軌結(jié)果的影響。

圖2 偽距擬合殘差圖Fig.2 Residual of pseudorange polynomial fitting

圖3 載波相位平滑偽距擬合殘差圖Fig.3 Residual of pseudorange polynomial fitting of carrrier phase

表2 臨潼站數(shù)據(jù)擬合噪聲統(tǒng)計表(單位：m)Tab.2 Polynomial fitting noise of Lintong station(unit：m)

4.1 數(shù)據(jù)預(yù)處理

以2008年3月4臨潼站相位觀測數(shù)據(jù)為例，分別使用方案2和方案3對觀測數(shù)據(jù)進(jìn)行預(yù)處理，比較其結(jié)果。

如圖4所示，該時間段內(nèi)相位數(shù)據(jù)發(fā)生了較大的周跳，兩種方法都能夠?qū)⑵涑晒μ綔y。多項式擬合法可以確定某一頻率發(fā)生周跳的位置和大小，并將其成功修復(fù);而無論哪一個頻率相位數(shù)據(jù)發(fā)生周跳，組合歷元間求差都會以粗差形式顯現(xiàn)，但由于該段載波相位數(shù)據(jù)頻繁發(fā)生周跳，按照方案2的數(shù)據(jù)分段思想，會使得段內(nèi)數(shù)據(jù)量較少，可能會影響定軌結(jié)果。

圖4 周跳探測結(jié)果Fig.4 Results of cycle slip detection

4.2 精密定軌

采用3種不同策略對觀測數(shù)據(jù)進(jìn)行預(yù)處理后進(jìn)行精密定軌。

方案1：采用B1頻率窄相關(guān)偽距觀測量定軌;

方案2：利用φLC-PLC對B1、B3頻率組合進(jìn)行周跳探測，標(biāo)記周跳發(fā)生的位置將數(shù)據(jù)分段，段內(nèi)相位平滑偽距后定軌;

方案3：綜合利用高次差法和多項式擬合法分別對B1、B3頻率相位數(shù)據(jù)進(jìn)行周跳探測和修復(fù)，使用修復(fù)后的全弧段相位平滑偽距數(shù)據(jù)定軌(圖5)。

采用定軌殘差的RMS作為評定精度的指標(biāo)，統(tǒng)計結(jié)果如圖6和表3所示。

圖5 數(shù)據(jù)處理流程圖Fig.5 Flow chart of data processing

圖6 各站定軌殘差RMS統(tǒng)計Fig.6 Residual RMS statistics graph of orbit determination at each station

表3 各站定軌殘差RMS統(tǒng)計(單位：m)Tab.3 Residual RMS statistics of orbit determination at each station(unit：m)

結(jié)果表明：使用相位平滑偽踞數(shù)據(jù)可以獲得比使用偽踞觀測值更高的定軌精度;對于幅值較大的周跳，方案2、3中的兩種方法均能夠成功探測，此時方案2、3定軌精度基本相當(dāng)，如表2中長春站、昆明站等;對于頻繁發(fā)生周跳的觀測數(shù)據(jù)(如表2中上海站、臨潼站)，聯(lián)合使用高次差法和多項式擬合法進(jìn)行周跳的探測并將其修復(fù)，能夠改善定軌精度。對于較頻繁發(fā)生周跳的載波相位數(shù)據(jù)，采用修復(fù)周跳全弧段平滑的策略較分段平滑能夠獲得更高的定軌精度，分析其原因，對于式(4)應(yīng)用誤差傳播定律得到雙頻相位平滑偽距值方差為：

另外，對于漸進(jìn)穩(wěn)定的Hatch濾波器，當(dāng)濾波時間充分長后，其最優(yōu)濾波值將漸進(jìn)地不依賴于濾波初值的選?。?0］，而當(dāng)數(shù)據(jù)量較小時，初值中偽距的因素會對平滑偽距值產(chǎn)生影響。

5 結(jié)論

本文將φLC-PLC組合法、綜合利用高次差法和多項式擬合法應(yīng)用于COMPASS-MEO衛(wèi)星載波相位數(shù)據(jù)周跳探測，利用精密定軌結(jié)果對比分析不同周跳探測方法和預(yù)處理策略對精度的影響，在COMPASS數(shù)據(jù)處理過程中，應(yīng)根據(jù)實際觀測數(shù)據(jù)情況，綜合考慮周跳修復(fù)難度、精度要求、運算效率等因素，合理選擇周跳處理方法及預(yù)處理策略。

1 常志巧.利用非差觀測量確定GPS衛(wèi)星導(dǎo)航軌道的方法研究［D］.解放軍信息工程大學(xué)，2006.(Chang Zhiqiao.A study of the orbit determination of GPS satellite using zerodifference observation［D］.Information Engineering University，2006)

2 魏子卿，葛茂榮.GPS相對定位的數(shù)學(xué)模型［M］.北京：測繪出版社，1998.(Wei Ziqing and Ge Maorong.The mathematical model of GPS relative positioning［M］.Beijing：Publishing House of Surveying and Mapping，1998)

3 羅峰，姚宜斌，宋偉偉.綜合利用多項式擬合和載波相位變化率探測單頻GPS周跳［J］.全球定位系統(tǒng)，2007，32(5)：9-13.(Luo Feng，Yao Yibin and Song Weiwei.Detection for cycle slips using polynomial fitting and carrier phase rate method［J］.GNSS World of China，2007，32 (5)：9-13)

4 劉基余.GPS衛(wèi)星導(dǎo)航定位原理與方法［M］.北京：科學(xué)出版社，2006.(Liu Jiyu.GPS satellite navigation and positioning principles and methods of navigation［M］.Beijing：Science Press，2006)

5 嚴(yán)新生，王一強，白征東.聯(lián)合使用高次差法和 T urboEdit法自動探測、修復(fù)周跳［J］.測繪通報，2007，53 (9)：23-26.(Yan Xinsheng，Wang Yiqiang and Bai Zhengdong.Automatically detect and repair cycle slips using all iance of high difference and turboedit［J］.Bulletin of Surveying and Mapping，2007，53(9)：23-26)

6 Misra Pratap and Enge Per.Global positioning system：signals，measurements，and performance［M］.Lincoln，Mass：Ganga-Jamuna Press，2001.

7 唐衛(wèi)明，劉智敏.GPS載波相位平滑偽距精度分析與應(yīng)用探討［J］.測繪信息與工程，2005，30(3)：37-39.(Tang Weiming and Liu Zhimin.Precision evaluating and applications of GPS code pseudorange smoothing algorithm based on carrier observables［J］.Journal of Geomatics，2005，30(3)：37-39)

8 范士杰.GPS精密單點定位的算法實現(xiàn)與精度分析［D］.武漢大學(xué)，2004.(Fan Shijie.The algorithm and accuracy analysis of GPS precise point positioning［D］.Wu Han University，2004)

9 Hatch R.The synergism of GPS code and carrier measurements［A］.In：Proceedings of the Third International Geodetic Symposium on Satellite Doppler Positioning［C］.Las Cruces，New Mexico.，USA，1982，1 213-1 232.

10 付夢印，鄧印紅，張繼偉.Kalman濾波理論及其在導(dǎo)航系統(tǒng)中的應(yīng)用［M］.北京：科學(xué)出版社，2003.(Fu Mengyin，Deng Yinhong and Zhang Jiwei.Dynamic Kalman filtering theory and its application in navigation system［M］.Beijing：Science Press，2003)

INFLUENCE ON ACCURACY OF COMPASS ORBIT PRECISION DETERMINATION WITH DIFFERENT DATA PREPROCESSING METHODS

Zhou Wei1，2)，Hao Jinming1)and Jia Xiaolin3)

1)Institute of Surveying and Mapping，Information Engineering University，Zhengzhou 450052 2)63880 Troops，Luoyang 471003 3)Xi’an Research Institute of Surveying and Mapping，Xi’an710054

How to detect the cycle slip is always the key to data preprocessing for satellite orbit determination.The method of LC-PC combination，Multi-difference combined Polynomial fitting for cycle slip detecting is analyzed，and then the test of orbit precision determination using COMPASS observations is made.It shows that，as for the observations of COMPASS receiver，the accuracy of orbit determination will be improved by choosing the methods of cycle slip detecting and data preprocessing rationally.

COMPASS;cycle slip;φLC-PLCcombination;polynomial fitting;orbit determination

1671-5942(2011)06-0084-05

2011-06-07

國家自然科學(xué)基金(40604003)

周巍，女，1981年生，博士研究生，主要從事GNSS理論研究與數(shù)據(jù)處理.E-mail：zhouwei_0611@163.com

P227

A