龐晉煒, 陳 瀟

(1.中國科學院大學,北京 100089；2.中國科學院空天信息創(chuàng)新研究院,北京 100089)

衛(wèi)星定位技術(shù)先后經(jīng)歷了傳統(tǒng)單點定位、經(jīng)典相對定位和精密單點定位，與傳統(tǒng)單點定位和經(jīng)典相對定位相比，精密單點定位(precise point positioning, PPP)技術(shù)具有可獨立作業(yè)，不需要在已知點架設基準站，不受距離限制，數(shù)據(jù)處理簡單等優(yōu)勢，其應用前景不可限量[1]。網(wǎng)絡輔助導航是一種有效的導航增強技術(shù)，可明顯縮短衛(wèi)星導航的首次定位時間，提升終端定位性能[2]。隨著5G技術(shù)的發(fā)展，地面移動通信網(wǎng)絡將迎來新的時代。因此，如何利用地面網(wǎng)絡輔助增強衛(wèi)星精密單點定位是一個重要的研究課題。

目前，一些相關(guān)文獻已經(jīng)對網(wǎng)絡輔助的信息格式及傳輸方式展開研究。文獻[3]研究了國際海運事業(yè)無線電技術(shù)委員會(Radio Technical Commission for Maritime Services, RTCM)協(xié)議第三版本中3.1版本和3.2版本兩種格式之間的差異，并指出RTCM V3.2格式中輔助GNSS技術(shù)將主要應用于GNSS高精度實時定位領(lǐng)域。文獻[4]研究了A-GPS基于無線資源位置服務協(xié)議(radio resource location service protocol, RRLP)通過全球移動通信系統(tǒng)(global system for mobile communications, GSM)網(wǎng)絡利用信令控制方式下發(fā)輔助信息的方法。泛在無線信號指的是非專門用于導航定位的無線電信號，包括移動蜂窩信號、無線上網(wǎng)信號等。RTCM SC-104格式是國際海事無線電技術(shù)委員會104特別委員會(Radio Technical Commission for Maritime Services, Special Committee No.104)制定的DGPS數(shù)據(jù)通用格式。NMEA 0183是美國國家海洋電子協(xié)會(National Marine Electronics Association)為海用電子設備制定的標準格式。文獻[5]通過對RTCM SC-104和NMEA 0183這兩種國際標準格式進行擴展，使其包含泛在無線信號定位相關(guān)的輔助數(shù)據(jù)信息，實現(xiàn)輔助數(shù)據(jù)格式的標準化。但是，上述文獻中均未涉及5G網(wǎng)絡輔助精密單點定位的狀態(tài)空間表達式(state space representation, SSR)信息，可見SSR輔助數(shù)據(jù)在5G網(wǎng)絡中的傳輸格式問題有待解決。

根據(jù)第三代合作伙伴計劃(3rd generation partnership project, 3GPP)組織的規(guī)定[6]，5G網(wǎng)絡將使用長期演進定位協(xié)議(long term evolution positioning protocol, LPP)提供輔助增強功能，于是，基于LPP定位協(xié)議[7]提出了一種用于5G網(wǎng)絡傳輸?shù)腟SR輔助數(shù)據(jù)消息格式，并對其定位結(jié)果進行仿真分析。本文首先設計了SSR數(shù)據(jù)的LPP消息格式，然后針對該格式給出了解碼方案，最后通過仿真，利用SSR數(shù)據(jù)進行事后靜態(tài)精密單點定位，驗證消息格式的有效性，并對GPS系統(tǒng)PPP和GPS+GLONASS雙系統(tǒng)PPP定位結(jié)果進行比較分析。

1 LPP協(xié)議

LPP是用于位置服務器與目標設備間點對點通信的定位協(xié)議，該協(xié)議支持的定位方式有可觀察到達時間差定位法(observed time difference of arrival，OTDOA)、網(wǎng)絡增強全球?qū)Ш叫l(wèi)星系統(tǒng)(assisted-global navigation satellite system, A-GNSS)、增強型小區(qū)標識定位法(enhancedCell-ID，E-CID)、地面信標系統(tǒng)定位法(terrestrial beacon system，TBS)等。一次LPP會話中，位置服務器或目標設備可獲取相關(guān)的測量值或輔助數(shù)據(jù)等信息。一次LPP會話又包含一個或多個LPP事務，一個事務只能進行單一操作，比如只能進行輔助數(shù)據(jù)的傳輸，或只能進行位置信息的傳輸。在LPP協(xié)議層，不同的事務利用事務ID進行區(qū)分。每個LPP事務都會涉及一個或多個LPP消息，每個LPP消息包含一個或多個定位方式的相關(guān)信息。

LPP消息包括事務ID、事務終止標志符、序列號、確認位、LPP消息主體，其中，事務ID用于區(qū)分不同的事務，事務終止標志符指示事務是否結(jié)束，序列號用于在接收器上檢測重復的LPP消息，當LPP消息在控制平面上傳輸時，確認位字段包含在LPP確認和請求確認的任何LPP消息中。LPP消息主體用于標識消息類型并包含特定消息類型的信息，包括“請求功能”類型、“提供功能”類型、“請求輔助數(shù)據(jù)”類型、“提供輔助數(shù)據(jù)”類型、“請求位置信息”類型、“提供位置信息”類型等。在LPP消息主體為“提供輔助數(shù)據(jù)”類型時，其內(nèi)容支持A-GNSS、OTDOA等定位方式。本文主要討論“提供輔助數(shù)據(jù)”類型下，定位方式為A-GNSS，支持PPP的輔助數(shù)據(jù)類型的LPP消息格式，輔助數(shù)據(jù)類型包括SSR軌道改正參數(shù)、SSR鐘差改正參數(shù)和SSR碼偏差參數(shù)[7]。

狀態(tài)空間概念是為了提供單一GNSS誤差源的狀態(tài)信息而提出的，SSR信息可以分別對GNSS的各項誤差源進行描述。相對于觀測空間表示方法(observation space representation, OSR)，SSR更適合于GNSS單點定位的誤差改正[8]。GNSS的SSR向量包括下列幾個參數(shù)：衛(wèi)星軌道改正、衛(wèi)星鐘差改正、衛(wèi)星信號偏差(包括衛(wèi)星硬件或軟件碼偏差和載波相位偏差)電離層延遲參數(shù)、對流層延遲參數(shù)、狀態(tài)質(zhì)量指示參數(shù)[9]。文獻[6]給出的SSR輔助數(shù)據(jù)信息如表1～表3。

表1 SSR軌道改正參數(shù)

表2 SSR鐘差改正參數(shù)

表3 SSR碼偏差

2 誤差模型

2.1 觀測方程

偽距和載波相位觀測方程分別為

(1)

(2)

2.2 SSR改正

(3)

接下來進行坐標系轉(zhuǎn)換，將軌道改正量ΔO轉(zhuǎn)換到地心地固坐標系下，得到:

(4)

Xorbit=Xbroadcast-ΔX

(5)

精密鐘差由廣播星歷和SSR鐘差改正參數(shù)計算得出，設SSR數(shù)據(jù)參考時刻為t0，由廣播星歷得到的衛(wèi)星鐘差記為Tbroadcast，由SSR鐘差改正參數(shù)得到的鐘差改正量記為ΔT，精密鐘差記為Tclock。SSR鐘差改正參數(shù)分別記為C0、C1和C2，則時刻t的鐘差改正量為

(6)

式(6)中：Vc是真空中的光速。則精密鐘差為

Tclock=Tbroadcast-ΔT

(7)

衛(wèi)星和接收機終端在發(fā)射和接收信號時，不同頻率以及同一頻率不同測距信號通道間都存在時延偏差，不同信號通道之間的時延偏差稱為差分碼偏差。硬件延遲是GNSS誤差源之一，差分碼偏差可用于消除硬件延遲[9-10]。

2.3 改正參數(shù)及模型選取

選取無電離層組合相位和偽距觀測值，截止高度角為15°，對流層采用Saastamonien模型改正，利用雙頻接收信號消除電離層影響，硬件延遲改正數(shù)、軌道改正數(shù)和鐘差改正數(shù)由NTRIP服務器提供，同時采用相對論效應、天線相位中心、固體潮等較為成熟的改正模型[11-13]。定位流程如圖1所示。

圖1 定位流程圖Fig.1 Positioning flow chart

3 關(guān)鍵模塊及傳輸流程

3.1 LPP消息格式設計

針對GPS和GLONASS系統(tǒng)進行設計，SSR輔助數(shù)據(jù)的LPP消息是位串，由圖2表示，其中第一行的最左位是第一位以及最高位，最后一行的最右位是最后一位以及最低位。該消息格式字節(jié)對齊(即8位的倍數(shù))，各字段比特數(shù)如表4所示。

圖2 SSR數(shù)據(jù)LPP消息格式Fig.2 SSR data LPP message format

表4 LPP消息各字段比特數(shù)

LPP消息分為消息頭和消息主體,LPP消息頭包括事務ID(TsID)、保留位、事務結(jié)束標志位(F)和序列號(SN)。事務ID由消息類型指示符(type)(8比特)和消息來源指示符(re)(1比特)組成。規(guī)定消息類型指示符為“00001111”時，傳輸支持精密單點定位的SSR輔助數(shù)據(jù)，消息來源符指示消息來自位置服務器還是終端。規(guī)定消息來源指示符為1時，LPP消息來自位置服務器，消息來源指示符為0時，LPP消息來自終端。保留位供今后設計使用，占6比特。事務結(jié)束符占1比特，指示當前事務是否結(jié)束。序列號占8比特，序列號用于在接收器上檢測重復接收的LPP消息。

消息主體包含歷元(e)、SSR軌道改正參數(shù)更新周期(UI1)、SSR鐘差改正參數(shù)更新周期(UI2)、SSR碼偏差參數(shù)更新周期(UI3)、衛(wèi)星參考基線(RD)、SSR改正數(shù)IOD(iodssr)、保留位、軌道改正支持衛(wèi)星數(shù)量指示符(OCsN)、SSR軌道改正參數(shù)集(L1)、鐘差改正支持衛(wèi)星數(shù)量指示符(CCsN)、SSR鐘差改正參數(shù)集(L2)、碼偏差支持衛(wèi)星數(shù)量指示符(CBsN)和碼偏差參數(shù)集(L3)。其中，歷元時間范圍為0～604 799 s；UI1、UI2和UI3分別占4比特，分別對應不同的更新時間，對應關(guān)系如表5所示，消息的發(fā)送周期取三者的最大公約數(shù)；SSR改正數(shù)IOD用于表示SSR信息和廣播星歷之間的對應關(guān)系；支持衛(wèi)星數(shù)量指示符表示該消息體支持改正的衛(wèi)星顆數(shù)，最大支持64顆衛(wèi)星；改正參數(shù)集包含改正用到的有效參數(shù)。

表5 UI與SSR更新周期對照表

圖3 E1、E2和E3結(jié)構(gòu)圖Fig.3 Structure of E1, E2 and E3

表6 E1、E2和E3字節(jié)數(shù)

3.2 LPP解碼模塊

首先識別消息前8位，獲取該LPP消息傳輸信息類型，接著獲取第2、3字節(jié)，了解LPP消息的來源、事務是否結(jié)束等信息，解碼流程如圖4所示。在解析E1、E2和E3參數(shù)包的過程中，原始SSR參數(shù)是帶符號位的二進制比特序列，需要先轉(zhuǎn)換為十進制，再乘以相應的精度因子，最終得到實際的SSR參數(shù)。

圖4 LPP消息解碼流程圖Fig.4 LPP message decoder

3.3 SSR數(shù)據(jù)通信仿真流程

本文使用RTKLIB進行仿真試驗，選取合適的掛載點，接收并保存RTCM3格式SSR數(shù)據(jù)，通過RTCM轉(zhuǎn)LPP格式模塊處理RTCM3格式的SSR數(shù)據(jù)，并生成LPP格式消息。參考文獻[14]使用MATLAB語言仿真5G物理廣播信道，并在其中傳輸LPP格式消息，接收端按照解碼方案解碼獲取SSR數(shù)據(jù)。圖5為SSR數(shù)據(jù)通信仿真流程圖。

圖5 SSR數(shù)據(jù)通信仿真流程圖Fig.5 SSR data transmittion simulation flow chart

4 性能分析

4.1 信息量分析

SSR消息長度隨更新周期、衛(wèi)星數(shù)和信號類型數(shù)變化而變化，在仿真試驗中，取軌道改正數(shù)更新周期為5 s，鐘差改正數(shù)更新周期為5 s，碼偏差改正數(shù)更新周期為5 s，衛(wèi)星數(shù)OCsN=CCsN=CBsN=x，支持y種信號碼偏差改正，那么，在60 s的時間里，SSR消息總字節(jié)數(shù)記為B1，計算公式為

B1=12[8+(1+17x)+(1+10x)+

(1+5xy)]

(8)

從圖6中可以看出，支持的信號類型數(shù)一定時，總字節(jié)數(shù)隨衛(wèi)星數(shù)的增加而線性增加。當衛(wèi)星數(shù)為10，信號類型數(shù)為2時，SSR消息字節(jié)數(shù)為1 980 bytes；當衛(wèi)星數(shù)為18，信號類型數(shù)為2時，SSR消息字節(jié)數(shù)為2 940 bytes。

圖6 60 s內(nèi)總字節(jié)數(shù)Fig.6 The size in bytes within 60 seconds

4.2 定位結(jié)果與分析

選取CLK93掛載點，接收并保存RTCM3格式的SSR信息，靜態(tài)情況下利用NOVATEL接收機接收實際的衛(wèi)星信號，并將觀測值和星歷保存為renix文件。按照定位流程圖，利用這些源數(shù)據(jù)分別進行GPS系統(tǒng)事后PPP定位和GPS+GLONASS雙系統(tǒng)事后PPP定位，定位結(jié)果如圖7和圖8所示。

圖7 GPS系統(tǒng)PPP結(jié)果Fig.7 PPP results of GPS single system

圖8 GPS+GLONASS雙系統(tǒng)PPP結(jié)果Fig.8 PPP results of GPS+GLONASS dual-system

從圖7和圖8中可以看出，GPS+GLONASS雙系統(tǒng)PPP整體收斂至0.2 m約需40 min， GPS系統(tǒng)PPP整體收斂至0.2 m約需50 min，可見，在收斂時間上，GPS+GLONASS雙系統(tǒng)比GPS系統(tǒng)更快，并且從圖中可以看出，雙系統(tǒng)更加穩(wěn)定。

本文將觀測數(shù)據(jù)按照5、10、15、20、25、30、40、50、60、90、120 min不同的觀測時長，分別統(tǒng)計了兩種定位方式下PPP所達到的定位精度，如表6所示。可以看到在第5 min的時候，GPS系統(tǒng)在E、N和U方向上的均方根分別為0.97、0.78、0.55 m，GPS+GLONASS雙系統(tǒng)在E、N和U方向上的均方根分別為0.42、0.11和0.90 m，均在1 m以內(nèi)。相比GPS系統(tǒng)，GPS+GLONASS雙系統(tǒng)在E和N方向上有明顯的改善，分別有40%和64%的提升，而在U方向上，兩者精度相當。在整個觀測過程中，兩種定位方式的可見星數(shù)量統(tǒng)計表如表7，再結(jié)合圖6，可以發(fā)現(xiàn)，雙系統(tǒng)以1 040 bytes/min的額外信息換取了上述精度的提升。

表7 衛(wèi)星數(shù)統(tǒng)計表

表6 PPP精度(RMS)統(tǒng)計表

5 結(jié)論

利用MATLAB仿真5G物理廣播控制信道，模擬了LPP消息從生成到解碼的過程，并利用解碼得到的SSR信息進行了事后靜態(tài)PPP定位。定位結(jié)果表明，不論是GPS系統(tǒng)PPP定位，還是GPS+GLONASS雙系統(tǒng)PPP定位，定位精度均能夠達到分米級，在E、N方向上，GPS+GLONASS雙系統(tǒng)精度明顯優(yōu)于GPS系統(tǒng)，以每分鐘1 040 bytes的額外信息換取了40%和64%的精度提升，在U方向上，GPS系統(tǒng)與GPS+GLONASS雙系統(tǒng)定位精度相當，得到以下結(jié)論：

(1)消息格式可用于事后靜態(tài)PPP定位。

(2)隨著消息中衛(wèi)星信息量的增大，定位精度會相應增大。