實時動態(tài)GPS數(shù)據(jù)處理的研究

易恒　樊東昊

【摘 要】針對GPS實時動態(tài)數(shù)據(jù)的特點，本文利用離散卡爾曼濾波算法進(jìn)行解算， 選用了重力船載GPS數(shù)據(jù)進(jìn)行實驗，實驗結(jié)果證明，對于船載短基線的解算，該算法可以獲得較高精度的定位結(jié)果，從而為解算重力場模型提供可靠的定位信息。

【關(guān)鍵詞】實時動態(tài)；GPS數(shù)據(jù)解算；短基線

【Abstract】According to the characteristics of GPS real-time dynamic data，this paper uses the discrete Kalman filter algorithm to solve the problem.The GPS data of gravity ship is used to experiment.The experimental results show that the algorithm can obtain high accuracy Positioning results，so as to solve the gravity field model to provide reliable positioning information.

【Key words】Real-time dynamic；GPS data calculation；Short baseline

0 引言

動態(tài)相對定位就是在固定站點上固定一個接收機，在運動的載體上架上觀測接收機。在運動的過程中，所有的接收機進(jìn)行同步觀測，用來確定運動載體相對于基準(zhǔn)站的瞬時位置。動態(tài)相對定位的特點就是實時確定載體的位置。

差分動態(tài)定位是利用安置在一個運動載體上的接收機和安置在地面上的一個或多個基準(zhǔn)站的接收機聯(lián)合測得改運動載體的位置，所以差分動態(tài)定位也稱為相對動態(tài)定位。根據(jù)實時性要求不同，差分動態(tài)定位分為實時差分動態(tài)定位和后處理差分動態(tài)定位。實時差分動態(tài)定位需要建立無線電數(shù)據(jù)傳輸，在觀測的同時解算出載體的位置；后處理差分動態(tài)不需要實時傳輸數(shù)據(jù)，而是在觀測結(jié)束后進(jìn)行處理。

針對GPS實時動態(tài)數(shù)據(jù)的特點，本文利用離散卡爾曼濾波算法進(jìn)行解算， 選用了重力船載GPS數(shù)據(jù)進(jìn)行實驗，實驗結(jié)果證明，對于船載短基線的解算，該算法可以獲得較高精度的定位結(jié)果，從而為解算重力場模型提供可靠的定位信息。

1 GPS動態(tài)定位中的離散卡爾曼濾波算法

如果測量誤差彼此之間是完全獨立的，那么應(yīng)用經(jīng)典的最小二乘理論是非常有效的，但是在處理動態(tài)數(shù)據(jù)時，其中一些誤差和時間是相關(guān)的，在這種情況下，卡爾曼濾波理論作為一種重要的最優(yōu)化估計理論被廣泛應(yīng)用于各種動態(tài)的數(shù)據(jù)處理中，尤其是在GPS單歷元定位中，近十年來得到了深入的發(fā)展和廣泛的應(yīng)用。Track就是采用離散的卡爾曼濾波算法。

在GPS數(shù)據(jù)處理時，觀測方程都是經(jīng)過線性化的，所以我們只討論線性系統(tǒng)的卡爾曼濾波算法。

1.1 觀測方程

1.3 計算步驟

整個卡爾曼濾波包括了3個步驟

2 實驗與分析

該數(shù)據(jù)為測定臺灣某海域內(nèi)海水重力，使用GPS觀測得到相關(guān)數(shù)據(jù)后進(jìn)行解算得到重力資料。對于精密星歷，可以從SOPAC網(wǎng)站下載了觀測日兩天內(nèi)的星歷。在本文中只對該數(shù)據(jù)的精度進(jìn)行解算。在海研一號船上布設(shè)了四臺GPS接收機作為觀測站，在岸邊布設(shè)了一臺GPS接收機作為基準(zhǔn)站。處理的數(shù)據(jù)時間段為一天。

船上四臺接收機分別為A、B、C、D，其分布及相關(guān)位置如圖：

該圖中顯示的坐標(biāo)為工程坐標(biāo)，是以船中線為坐標(biāo)軸。不管任意時間段內(nèi)，A-B A-D B-C C-D 的距離是不變的，可以利用這一點作為評定數(shù)據(jù)的精度。

從上圖可以得出各個站點之間的距離。

使用Track解算出各個站點的X、Y、Z坐標(biāo)后，然都算出各個歷元相鄰站點間的距離可得如圖：

從圖中可以看到解算出來的基線長與實際長度的差距在大部分幾厘米之內(nèi)，其中紅色線表示為實際長度。可看出解算GPS動態(tài)數(shù)據(jù)時的精度可達(dá)cm級。在圖中，基線A-B與實際長度有一定的差距，可能是由于GPS接收機擺放位置的原因?qū)е掠^測數(shù)據(jù)較差，從而使得基線解算出現(xiàn)較大的偏差。

從以上可以看出在海上，解算出的定位精度是比較高的，而且還可以提供航船的軌跡以及潮汐高程變化等信息。

3 總結(jié)

本文利用離散卡爾曼濾波算法進(jìn)行解算船載動態(tài)GPS數(shù)據(jù)，通過解算船載4站GPS位置信息，得到短基線的精度可達(dá)厘米，可以為重力船獲得重力數(shù)據(jù)的解算提供精度較高的定位信息。

【參考文獻(xiàn)】

[1]李征航，黃勁松.GPS測量與數(shù)據(jù)處理[M].武漢大學(xué)出版社，2010.

[2]蘇小寧.高頻GPS單歷元定位方法及其在地震學(xué)中的應(yīng)用與研究[J].2011.

[3]韓曉東.GPS單歷元動態(tài)定位的研究[J].2006.

[4]李林紅.何秀鳳.利用動態(tài)定位模塊Track進(jìn)行GPS單歷元定位研究[J].全球定位系統(tǒng)，2011.

[5]Thomas Herring.track short course.MIT，2012.

[6]GAMIT-GLOBKHomepage[EB/OL].http：//www-gpsg.mit.edu/～simon/gtgk/.

[7]蘇小寧，孟國杰，等.基于TRACK進(jìn)行GPS單歷元定位[J].大地測量與地球動力學(xué)，2009.

[8]楊志強，王利，等.GPS 監(jiān)測青藏板塊運動數(shù)據(jù)的GAMIT/GLOBK處理方法及有關(guān)問題[J].西安工程學(xué)院學(xué)報，1999.

[9]馬洪斌.賀黎明.新版GAM IT軟件的功能特點與應(yīng)用實例分析[J].礦山測量，2008.

[10]TRACK.GPS differential phase kinematic positioning program[EB/OL].http：//geoweb.mit.edu/～simon/gtgk/help/track.hlp.htm.2011-07-28.

[11]喻國榮.單歷元模糊度解算問題[J].測繪通報，2003，11：6-7.

[12]孫亞峰.單歷元定位技術(shù)及其在的地震觀測中的應(yīng)用[D].上海：同濟大學(xué)，2009.

[13]周星.基于GAMIT TRACK和Bernese PPP的地震監(jiān)測比較[J].測繪信息與工程，2010.

[責(zé)任編輯：田吉捷]