徐亞楠，劉舜

(1.廣州番禺職業(yè)技術學院，廣州 511483；2.陜西鐵路工程職業(yè)技術學院，陜西 渭南 714099)

0 引 言

全球?qū)Ш叫l(wèi)星系統(tǒng)(GNSS)以北斗衛(wèi)星導航系統(tǒng)(BDS)、GPS、Galileo 和GLONASS 四大系統(tǒng)為主，當前每個系統(tǒng)都能單獨提供全球定位服務，且部分衛(wèi)星能播發(fā)三頻甚至三頻以上的信號，為GNSS 精密數(shù)據(jù)處理提供了新思路[1-3].Galileo 系統(tǒng)是歐盟建設和研發(fā)的衛(wèi)星導航定位系統(tǒng)，于2016 年12 月15 日正式聲明公開服務，其建設計劃最早可追溯到20 世紀90 年代，該系統(tǒng)計劃由30 顆衛(wèi)星組成，且能播發(fā)E1、E5a、E5b 等多頻信號[4-7].雖然Galileo 系統(tǒng)正式開通服務時間較短，但國內(nèi)很多學者對其定位性能進行了研究.文獻[8]基于多個MGEX(Multi-GNSS Experiment)跟蹤站數(shù)據(jù)分析校正天線參數(shù)對Galileo系統(tǒng)雙頻精密單點定位(PPP)的影響，研究結果表明，與消電離層定位結果相比，水平定位精度相當，高程方向定位精度提升較大.文獻[9]推導了適用于多頻非組合(UC)觀測值的GNSS 衛(wèi)星精密定軌模型，并且基于多測站Galileo 雙頻和三頻數(shù)據(jù)進行定軌，發(fā)現(xiàn)雙頻UC 模型與IF 模型定軌精度相當，三頻UC模型和IF 模型浮點解精度較E1/E5b 雙頻組合結果有一定的提升.文獻[10]對比分析了Galileo 系統(tǒng)雙頻無電離層組合、雙頻非差非組合和三頻非差非組合PPP 精度以及收斂時間，發(fā)現(xiàn)Galileo 三頻PPP精度較優(yōu)，且相比雙頻PPP 精度有明顯提升，但較收斂時間的提升并不明顯.文獻[11]利用非組合模型解算了Galileo 系統(tǒng)三頻組合PPP 精度，發(fā)現(xiàn)三頻PPP精度與收斂時間優(yōu)于傳統(tǒng)雙頻PPP，尤其是在觀測數(shù)據(jù)異常時效果更為明顯.文獻[12]分析了GPS/Galileo組合PPP 精度，發(fā)現(xiàn)雙系統(tǒng)組合PPP 具有較好的定位結果，且定位精度與收斂時間較單系統(tǒng)具有一定的改善.綜合當前各學者對Galileo 系統(tǒng)研究的內(nèi)容與方法，本文基于自測試9 km 短基線，提出利用一種UC 定位模型解算Galileo 雙頻與三頻短基線數(shù)據(jù)，并且進一步分析了不同頻率組合下定位精度.

1 數(shù)學模型

在短基線定位中，雙差函數(shù)模型是常用的解算模型.該模型能有效消除接收機鐘差、衛(wèi)星鐘差，可以最大程度地削弱電離層延遲、對流層延遲和衛(wèi)星軌道誤差的影響，有利于模糊度的固定[13].雙差載波相位與偽距觀測方程[14]一般表示如下：

在利用UC 模型定位時，電離層誤差無法消除，因此將電離層誤差當作未知參數(shù)進行估計，在式(1)的基礎上可以進一步得到三頻非組合短基線定位模型[11,15-16]，表示如下：

式中：I FB 為頻間偏差；γ 為不同頻率間的電離層因子；其余符號表示含義與之前相同.

在式(1)、(2)的基礎上進一步利用卡爾曼濾波進行參數(shù)估計，利用LAMBDA 算法進行模糊度固定，最終解算得到接收機坐標.

2 數(shù)據(jù)解算分析

為詳細評估非組合模型下Galileo 系統(tǒng)雙頻與三頻短基線解算性能，試驗設計了一條約9 km長的短基線，數(shù)據(jù)采集時高度截止角設置為12°，采集時間為2020 年3 月12 日全天，采集間隔為30 s，接收機類型為TRIMBLE ALLOY，天線類型為TRM59800.00.

數(shù)據(jù)解算軟件為根據(jù)RTKLIB 編譯的開源程序，以GAMIT 軟件解算得到的靜態(tài)坐標作為參考坐標，對流層改正模型采用Saastamoinen 模型，其他誤差通過對應模型進行改正.在進行數(shù)據(jù)解算時，首先解算得到Galileo 系統(tǒng)E1/E5a、E1/E5b 以及E5a/E5b 雙頻組合短基線定位結果，然后解算得到E1/E5a/E5b 三頻組合短基線定位結果，最后根據(jù)解算得到的定位結果與參考坐標計算得到定位偏差與定位精度.

如圖1 所示，給出了Galileo 系統(tǒng)的衛(wèi)星可用數(shù)與位置精度因子(PDOP)值隨時間變化情況.由圖1可知，在整個觀測期內(nèi)，Galileo 系統(tǒng)的衛(wèi)星可用數(shù)為4～9 顆，平均衛(wèi)星可用數(shù)為6 顆.Galileo 系統(tǒng)的PDOP 值除個別歷元外，PDOP 值在4 cm 以內(nèi)，平均PDOP 值為2.13 cm.如圖2 所示，給出了Galileo 系統(tǒng)雙頻與三頻解算歷元率，由圖2 可知，Galileo 系統(tǒng)雙頻與三頻短基線歷元解算率都在95%以上，三種雙頻歷元解算率相當，而三頻歷元解算率較雙頻有一定提升.

圖1 衛(wèi)星可用數(shù)與PDOP 值

圖2 歷元解算率

根據(jù)解算得到坐標與參考坐標，計算得到水平方向與高程方向的定位誤差，Galileo 系統(tǒng)雙頻與三頻水平向與高程向定位誤差，如圖3～圖6 所示.

由圖3～圖5 可知，E1/E5a 和E5a/E5b 組合短基線定位E 方向定位誤差小于0.4 m，N 方向定位誤差小于0.3 m，U 方向定位誤差小于0.6 m，E1/E5b 組合短基線定位E 方向定位誤差小于0.3 m，N 方向定位誤差小于0.3 m，U 方向定位誤差小于0.6 m.由圖6可知，Galileo 三頻短基線定位誤差較雙頻有所減小，E 方向定位誤差小于0.3 m，N 方向定位誤差小于0.3 m，U 方向定位誤差小于0.4 m.

圖3 E1/E5a 組合短基線水平與高程向定位誤差

圖4 E1/E5b 組合短基線水平與高程向定位誤差

圖5 E5a/E5b 組合短基線水平與高程向定位誤差

圖6 E1/E5a/E5b 組合短基線水平與高程向定位誤差

為進一步直觀的表示Galileo 雙頻、三頻定位精度以及三頻較雙頻短基線定位精度的提升，統(tǒng)計了不同頻率組合東(E)、北(N)、天(U)三個方向的定位精度以及三頻較雙頻短基線定位精度的提升量，如表1所示.

表1 Galileo 雙頻與三頻短基線基線解算精度均方根(RMS)及三頻較雙頻定位精度提升統(tǒng)計

由表1 可知，E1/E5a 和E5a/E5b 短基線E 方向和U 方向定位精度相當，E 方向定位精度優(yōu)于3.5 cm，U 方向定位精度優(yōu)于9 cm，E1/E5a 組合N 方向定位精度優(yōu)于5 cm，E5a/E5b 組合N 方向定位精度優(yōu)于4 cm，E1/E5b 組合定位精度優(yōu)于E1/E5a 和E5a/E5b，E 方向定位精度優(yōu)于2 cm，N 方向定位精度優(yōu)于4.5 cm，U 方向定位精度優(yōu)于5 cm.E1/E5a/E5b 三頻組合短基線定位精度優(yōu)于任一雙頻組合，E 方向定位精度優(yōu)于1.5 cm，N 方向定位精度優(yōu)于2.5 cm，U 方向定位精度優(yōu)于3.5 cm.同時發(fā)現(xiàn)，三頻組合短基線定位精度較雙頻定位精度有明顯提升，提升量在30%～60%.

3 結束語

本文基于9 km 短基線實測數(shù)據(jù)，分析了非組合模型下Galileo 雙頻與三頻短基線解算精度，結果表明，Galileo 衛(wèi)星可用數(shù)和PDOP 值較優(yōu)，雙頻歷元解算率相當，三頻歷元解算率較雙頻有明顯提升.雙頻短基線定位中，E1/E5b 組合定位精度要優(yōu)于另外兩種雙頻組合定位精度，但是定位精度都在cm 級，水平定位精度優(yōu)于5 cm，高程精度優(yōu)于10 cm.三頻組合短基線定位精度較雙頻組合有明顯提升，E 方向定位精度優(yōu)于1.5 cm，N 方向定位精度優(yōu)于2.5 cm，U 方向定位精度優(yōu)于3.5 cm.