基于移動基準站的艦船系統(tǒng)試驗數(shù)據(jù)處理方法研究

彭曉剛，馮 欣，李景巖

(1.91404部隊，河北 秦皇島 066000；2.66444部隊,北京 100042)

基于移動基準站的艦船系統(tǒng)試驗數(shù)據(jù)處理方法研究

彭曉剛1，馮 欣2，李景巖1

(1.91404部隊，河北 秦皇島 066000；2.66444部隊,北京 100042)

針對艦船系統(tǒng)試驗基準數(shù)據(jù)處理利用GNSS定位數(shù)據(jù)通過固定基準站差分解算后再獲得相對關系解算的不足，分析該方法面臨的難題及難以滿足遠距離大航區(qū)試驗的現(xiàn)實問題，提出基于載波相位實時動態(tài)差分處理的移動基準站解算方法。在研究GNSS差分定位原理基礎上，結合經(jīng)典載波相位雙差方法，通過單歷元實時差分，建立無固定基準站的相對數(shù)據(jù)解算模型。基于采集的實測數(shù)據(jù)，分析基線解算精度，驗證所提出方法在艦船系統(tǒng)試驗基準數(shù)據(jù)處理中的可行性和質(zhì)量效果。

艦船系統(tǒng)試驗；基準數(shù)據(jù)處理；移動基準站；差分處理

艦船系統(tǒng)試驗數(shù)據(jù)處理主要利用全球?qū)Ш叫l(wèi)星系統(tǒng)(Global Navigation Satellite System，GNSS)測量數(shù)據(jù)計算目標平臺相對于被試裝備的空間幾何關系(距離、方位、高度或俯仰等)，并把該數(shù)據(jù)作為檢驗被試裝備性能的基準數(shù)據(jù)，以判斷被試裝備所獲取數(shù)據(jù)與基準數(shù)據(jù)的差異，進而分析被試裝備精度參數(shù)，以檢驗、鑒定裝備的性能指標[1]。

以往基準數(shù)據(jù)測量通過架設一個或多個固定基準站加活動站的模式，采用商業(yè)差分軟件對被試裝備架設平臺及目標平臺接收機原始數(shù)據(jù)進行處理，先分別獲得各平臺位置信息(通常是經(jīng)度、緯度、高度)后，再通過空間模型轉(zhuǎn)換，得到目標相對于被試裝備的相對關系。該處理模式存在諸多不足，如：需要基準站觀測數(shù)據(jù)，處理過程相對復雜，動用測控裝備多；試驗中航路設計、平臺布設一定程度上受到基準站架設條件制約；某些地區(qū)突發(fā)性緊急項目獲得固定基準站坐標難度較大；基準站需要布置在陸地、海島上，遠離陸地的動態(tài)試驗數(shù)據(jù)處理如果仍舊采用固定基準站，移動站與基準站的距離將會變得很長等等[2-4]。

針對上述不足，本文在分析傳統(tǒng)GNSS差分處理技術的基礎上，針對艦船系統(tǒng)試驗特點，探索建立單歷元載波相位動態(tài)差分相對定位模型，以實現(xiàn)無基準站的相對數(shù)據(jù)解算。最后通過實例驗證算法的可靠性和有效性，探討處理精度對試驗影響。

1 差分觀測模型

試驗中，構成解算基線的目標平臺與被試平臺距離動態(tài)變化，因此需要構造不受距離約束且不考慮運動狀態(tài)的幾何觀測量[2,5-7]。對原始觀測數(shù)據(jù)進行周跳的探測與修復[6]，并對其濾波以消除隨機噪聲影響；再進行站間和星間雙差處理，以消除部分公共誤差；建立基于雙差模型的觀測方程，通過平差，獲得模糊度浮點解的精確估計；最后，應用算法搜索整周模糊度的固定解[6-9]。關于周跳探測與修復方法可以參看作者在文獻[6]的論述，本文側(cè)重探討其他相關處理。

以BDS接收機數(shù)據(jù)為例，在載波相位差分定位中，接收機i觀測衛(wèi)星r的載波相位觀測方程為[9-11]

(1)

(2)

進行雙差處理[1-2,5]，得到

(3)

(4)

2 整周模糊度求解

(5)

式中：A為系數(shù)陣；X由bur向量分量組成；N由ΔNi,j向量組成；ε是觀測噪聲向量。記投影算子為

(6)

式中：I為單位陣；P為觀測值權陣；R為載波相位雙差殘差靈敏度矩陣。

用R左乘以式(5)的兩邊，令V=Rε(V為殘差向量)，考慮RA=0，則有

(7)

得模糊度向量的標準差為(n為可視衛(wèi)星的個數(shù))

(8)

求解正確模糊度的目標函數(shù)可變?yōu)椋篤TPV=min。設初始相位殘差向量為

(9)

欲使殘差向量的值最小，即等價于

(10)

則

(11)

式中：b0為初始相位殘差向量，單位為周。在兩邊左乘以R矩陣，則：

(12)

由于R是冪等陣，且其秩為n-3，經(jīng)過變換，可得

(13)

式中，R1=[R11R21]T，R11為R的分解，其大小為(n-1)×(n-3)，R21為R的分解，其大小為3×(n-3)，V1=ΔN1-ΔN1f。由結果結合其他解算處理，即可確定唯一正確的整周模糊度，進而再得到基線解b[6-7,14]。

3 模型驗證及結果分析

為了驗證算法的精度、正確性，分別進行靜態(tài)和動態(tài)實驗，對實際采集數(shù)據(jù)用上述算法進行解算，并將所得結果與標準值進行比較以驗證算法的可靠性。

3.1 靜態(tài)驗證

利用某地CORS站HLD1模擬測試平臺，某臨時架設測點為目標平臺，采用2016-05-02上午的觀測數(shù)據(jù)，并用另一CORS站HLD2數(shù)據(jù)對臨時點差分后的結果解算作為已知值(測試平臺和目標平臺距離：1 569.767 4 m，方位：244.163 0°，俯仰：0.304 9°)，采用上述方法進行處理，得到的基線向量與已知值做差。解算后，模糊度固定解8 430個歷元(占89%)，浮點解1 015個歷元(占10.7%)，單點解24個歷元(0.3%)。得到解算結果和誤差曲線如圖1～圖4所示。

圖1 解算結果

圖2 解算基線誤差分布

圖3 解算方位誤差分布

圖4 解算俯仰誤差分布

由結果剔除個別異常數(shù)據(jù)可知，基線誤差在0.2 m附近，方位誤差在±0.02°范圍內(nèi)，俯仰誤差在±0.03°范圍內(nèi)。

3.2 動態(tài)驗證

利用2015-10-29上午實際動態(tài)試驗數(shù)據(jù)(航路見圖5)，以某基準站觀測數(shù)據(jù)用商業(yè)軟件waypoint差分完的結果解算兩個船之間的相對關系作為已知值，采用上述方法進行處理，得到的基線向量與已知值做差。解算后，固定解15 305個歷元(占57.3%)，浮點解11 399個歷元(占42.7%)。得到解算結果和誤差曲線如圖5～圖10所示。

圖5 實際航路

圖6 解算結果

圖7 距離變化趨勢

圖8 距離誤差解算結果

圖9 方位誤差解算結果

圖10 俯仰誤差解算結果

由結果剔除個別異常數(shù)據(jù)可知，距離誤差在±0.25 m附近，方位誤差在±0.01°范圍內(nèi)，俯仰誤差在±0.015°范圍內(nèi)。

3.3 結果分析

驗證一，由測試環(huán)境可知，CORS站架設條件好，周圍無遮擋、無干擾，數(shù)據(jù)接收穩(wěn)定。所用活動站接收機為動態(tài)測量模式，架設的附近有高樓和電線等影響，數(shù)據(jù)質(zhì)量不如CORS站設備，解算質(zhì)量有所下降。

驗證二，船上動態(tài)試驗設備架設環(huán)境復雜，周圍存在復雜電磁環(huán)境，且有甲板、海面、桅桿的多路徑信號，導致個別歷元解算模糊度不固定。所處理數(shù)據(jù)與商業(yè)差分軟件比對發(fā)現(xiàn)，在模糊度浮點解算的區(qū)間，商業(yè)軟件也存在浮點解算現(xiàn)象，說明是由于當時測量環(huán)境導致數(shù)據(jù)質(zhì)量下降造成的。浮點解導致的距離誤差基本在1倍波長附近，約20 cm左右，導致的方位誤差與基線長度近似成反比，基線越長，方位、俯仰精度相對越高。

艦船系統(tǒng)試驗中，一般的高精度警戒探測設備，距離精度在米級，方位精度在0.2°左右，俯仰精度在0.2°左右。由上述驗證結果可知，能夠滿足試驗需要。但是，對于更高精度的跟蹤、控制設備試驗，其精度要求高，上述算法暫不具備該功能，需要進一步深化研究。

4 結束語

在基準數(shù)據(jù)處理模式上進行探索，提出的方法具有可操作性，能夠?qū)崿F(xiàn)整周模糊度的快速解算，可以準確獲取較高精度相對定位信息。通過靜態(tài)和動態(tài)驗證分析表明，提出的新方法具有可行性，靜態(tài)和動態(tài)精度較高，驗證了算法的正確性和有效性。

算法解算的結果能夠滿足一定的試驗需求，可以節(jié)約試驗成本，提高處理效率，滿足遠離陸地的試驗環(huán)境需要。但在工程應用及對更高精度設備試驗方面還需做進一步深化研究。

[1] 彭曉剛.GPS姿態(tài)測量及其在海軍靶場試驗數(shù)據(jù)處理中的應用研究[D].鄭州:信息工程大學,2011.

[2] 劉立龍,鴻雁,唐詩華.基于移動基準站DGPS整周模糊度快速求解的研究[J]. 測繪信息與工程,2007,32(2):12-13.

[3] 段榮,趙修斌,龐春雷，等.一種GPS移動基準站精密相對定位新算法[J]. 四川大學學報(工程科學版), 2015,47(3):130-136.

[4] YANG Yunchun, SHARPE R T, HATCH R R. A Fast Ambiguity Resolution Technique for RTK Embedded Within a GPS Receiver [C]. The 15th International Technical Meeting of Satellite Division of the Institute of Navigation, Portland, Oregon, 2002.

[5] 汪捷,徐冠楠.基于GNSS動態(tài)差分相對定位方法的研究與探討[J].現(xiàn)代導航,2015(3):250-256.

[6] 彭曉剛,李景巖. 北斗系統(tǒng)姿態(tài)測量中周跳探測與修復方法[J]. 測繪科學, 2016, 41(9): 29-32.

[7] 任小偉.載波相位差分相對定位的模糊度求解[J]. 導航定位學報,2014,2(1):20-22.

[8] 王艷麗,雷繼兆,趙笛. 一種GPS載波相位差分相對定位算法[J]. 計算機測量與控制,2015,23(7):2447-2449.

[9] 劉玥.支持在航載波相位差分定位方法研究[D].哈爾濱：哈爾濱工程大學,2013.

[10] 趙偉,萬德均,劉建業(yè)，等. GPS載波相位整周模糊度的在航快速算法研究[J].中國慣性技術學報,2004,12(1):33-36.

[11] 馮超,田蔚風,金志華. GPS載波相位雙差整周模糊度在航求解[J]. 中國慣性技術學報, 2004,12(5):33-37.

[12] 鄒進貴,王海城,隗劍秋. RTK技術在航跡、航速測量中的應用與研究[J]. 測繪信息與工程, 2001(1):4-6.

[13] 唐衛(wèi)明,鄧辰龍,高麗峰．北斗單歷元基線解算算法研究及初步結果[J]．武漢大學學報(信息科學版),2013,38( 8):897-901.

[14] 劉立龍.動態(tài)對動態(tài)GPS高精度定位理論及其應用研究[D].武漢:武漢大學,2005.

[責任編輯：張德福]

A data processing method in warship system test based on moving base station

PENG Xiaogang1, FENG Xin2, LI Jingyan1

(1. Troops 91404，Qinhuangdao 066000，China;2.Troops 66444, Beijing 100042，China)

Aiming at the insufficiency of benchmark data processing in warship system test based on fixed base station differential calculation through GNSS positioning data, this paper analyses the difficulty that conventional method is facing with, which can not fit the test of long-range and great navigation area. It proposes the data processing method of carrier phase differential calculation based on moving base station. Based on these analyses of GNSS differential calculation principle, combined with classic carrier phase double-difference, it builds a differential calculation model based on moving base station through real-time single epoch processing. Using the data collected in real-time situation, it gets the baseline calculating precision, and validates the feasibility and result quality in warship system test benchmark data processing.

warship system test; benchmark data processing; moving base station; difference processing

引用著錄：彭曉剛，馮欣，李景巖.基于移動基準站的艦船系統(tǒng)試驗數(shù)據(jù)處理方法研究[J].測繪工程，2017，26(4)：12-15.

10.19349/j.cnki.issn1006-7949.2017.04.003

2016-07-24

彭曉剛(1978-)，男，高級工程師，博士.

P228

A

1006-7949(2017)04-0012-04