陳 杰，黃智剛，李 銳，韓春陽

(1.北京航空航天大學(xué) 電子信息工程學(xué)院，北京 100191；2.北京衛(wèi)星導(dǎo)航中心，北京 100094)

1 引言

目前，全球衛(wèi)星導(dǎo)航系統(tǒng)(global positioning system，GPS) 的開放服務(wù)定位精度為10 m左右[1]，尚不能滿足對精度和可靠性要求高的用戶需求，系統(tǒng)性能的提升技術(shù)已經(jīng)成為衛(wèi)星導(dǎo)航領(lǐng)域的研究熱點之一，其中差分增強技術(shù)作為有效提升全球衛(wèi)星導(dǎo)航系統(tǒng)性能的手段之一，大大推動了全球衛(wèi)星導(dǎo)航系統(tǒng)在不同領(lǐng)域的應(yīng)用。差分增強技術(shù)根據(jù)服務(wù)區(qū)域的范圍分類，可分為廣域和局域，根據(jù)觀測值分類，可分為基于碼偽距的差分增強和基于載波相位的精密定位，典型代表分別為廣域增強系統(tǒng)(wide area augmentation system， WAAS)、全球差分GPS系統(tǒng)(global differential GPS， GDGPS)、局域增強系統(tǒng)(local area augmentation system， LAAS)和連續(xù)運行參考站(continuously operating reference stations， CORS)。

為了便于測試分析和縮短研發(fā)周期，科研單位和商業(yè)機構(gòu)開發(fā)了多種仿真平臺，按照輸出方式可以分為射頻信號輸出和純數(shù)字輸出兩類。輸出射頻信號的典型仿真平臺是美國思博倫(Spirent)公司的GSS8000系統(tǒng)，該系統(tǒng)在控制器計算機上運行SimGEN軟件并通過信號發(fā)生器提供射頻輸出信號，可用于單個接收機的測試。受限于計算機硬件，組網(wǎng)的多臺機箱無法實現(xiàn)時鐘的納秒級同步，無法實現(xiàn)系統(tǒng)級仿真。純數(shù)字輸出的仿真平臺能夠克服時鐘同步的問題，歐洲航天局(European space agency，ESA)開發(fā)的伽利略系統(tǒng)仿真工具(Galileo system simulation facility， GSSF)[2]可配置星座、監(jiān)測站和用戶，具備衛(wèi)星可視性、監(jiān)測站覆蓋深度等初步分析功能。然而GSSF僅能按照仰角配置用戶等效測距誤差，無法單獨配置各個偽距誤差影響因子，此外用戶段只支持靜態(tài)用戶或者區(qū)域網(wǎng)格，未提供用戶接收機動態(tài)模擬功能，限制了用戶定制仿真場景的深度。用戶獲得ESA授權(quán)才能使用GSSF，進一步限制了其應(yīng)用范圍。

本文設(shè)計了面向差分和完好性研究需求的全球衛(wèi)星導(dǎo)航系統(tǒng)綜合仿真平臺，能夠靈活深入配置仿真參數(shù)，特有完好性故障仿真功能，具備用戶接收機多類型運動能力，采用接收機自主交換格式[3](receiver independent exchange format， RINEX)和標(biāo)準(zhǔn)產(chǎn)品3軌道格式[4](extended standard product 3 orbit format， SP3-c)輸出，向用戶插件開放數(shù)據(jù)接口，支持差分增強技術(shù)的系統(tǒng)級和用戶級仿真數(shù)據(jù)生成與分析。

2 總體設(shè)計

全球衛(wèi)星導(dǎo)航系統(tǒng)綜合仿真平臺主要用于衛(wèi)星導(dǎo)航系統(tǒng)的數(shù)字仿真，為差分增強技術(shù)研究提供高逼真度的數(shù)據(jù)源。該平臺允許用戶通過文本和圖形界面兩種方式靈活配置星座、偽距誤差、監(jiān)測站和用戶接收機等參數(shù)，具備完好性故障仿真的功能，按照普遍采用的RINEX和SP3格式輸出仿真數(shù)據(jù)，開放數(shù)據(jù)接口提高平臺可拓展能力。

該仿真平臺采用模塊化設(shè)計，分為五大模塊：控制模塊、交互模塊、數(shù)據(jù)庫模塊、運算模塊和輸出模塊，五大模塊的內(nèi)部結(jié)構(gòu)及模塊間關(guān)系如圖1所示。

圖1 全球衛(wèi)星導(dǎo)航系統(tǒng)綜合仿真平臺

交互模塊包括文本文件讀取模塊和圖形界面交互模塊，前者負責(zé)從文本讀入星座和監(jiān)測站配置信息，后者用于實現(xiàn)選擇衛(wèi)導(dǎo)系統(tǒng)和查看衛(wèi)星參數(shù)、查看和編輯監(jiān)測站參數(shù)、設(shè)定用戶接收機類型和運動狀態(tài)參數(shù)、配置偽距誤差模型參數(shù)、配置完好性故障參數(shù)、設(shè)定仿真時間等功能。

數(shù)據(jù)庫模塊用于存儲、讀取和修改衛(wèi)星數(shù)據(jù)庫、偽距誤差模型數(shù)據(jù)庫、監(jiān)測站數(shù)據(jù)庫、用戶接收機數(shù)據(jù)庫和完好性故障數(shù)據(jù)庫。衛(wèi)星數(shù)據(jù)庫包含星座、名稱、軌道、時鐘等參數(shù)，偽距誤差模型數(shù)據(jù)庫包含星歷、星鐘、電離層、對流層、接收機鐘差、熱噪聲及多徑誤差等模型參數(shù)，監(jiān)測站數(shù)據(jù)庫包含名稱、位置、時鐘、觀測高度截止角等參數(shù)，用戶數(shù)據(jù)庫包含位置、運動狀態(tài)、接收機觀測類型等參數(shù)，完好性故障數(shù)據(jù)庫包含星歷故障、星鐘故障、電離層風(fēng)暴、監(jiān)測站故障的起止時間和故障模型參數(shù)。

運算模塊包括偽距誤差模型、完好性故障模型和用戶接收機動態(tài)模型。偽距誤差模型分為星鐘誤差模型、星歷誤差模型、電離層誤差模型、對流層誤差模型、接收機鐘差模型和熱噪聲及多徑誤差模型，用于模擬信號從產(chǎn)生到接收運算過程中的關(guān)鍵影響因素。完好性故障模型分為星歷故障模型、星鐘故障模型、電離層風(fēng)暴模型、監(jiān)測站時鐘故障模型和接收機單通道故障模型。用戶接收機動態(tài)模型用于生成用戶接收機的實時位置，并計算多普勒頻移對觀測值的影響。

輸出模塊負責(zé)輸出仿真結(jié)果。為了便于數(shù)據(jù)交換，該平臺采用標(biāo)準(zhǔn)格式輸出數(shù)據(jù)：根據(jù)RINEX 2.10格式輸出導(dǎo)航電文和觀測數(shù)據(jù)，其中每個監(jiān)測站包含3臺接收機，經(jīng)過一致性校驗后輸出2臺接收機的觀測數(shù)據(jù)，觀測類型包括多個頻點的碼偽距和載波相位；根據(jù)SP3-c格式輸出精密星歷。此外能夠以動畫形式演示星下點軌跡。

控制模塊實現(xiàn)仿真平臺整體控制，提供接口用于各模塊間數(shù)據(jù)交換。此外，控制模塊具備調(diào)用用戶插件的功能，預(yù)置了格網(wǎng)電離層垂直誤差(grid ionospheric vertical error，GIVE)、接收機自主完好性監(jiān)測 (receiver autonomous integrity monitoring，RAIM)、用戶測距精度(user range accuracy， URA)、用戶差分距離誤差(user differential range error， UDRE)和監(jiān)測站覆蓋深度等插件的按鈕。

圖2 仿真平臺界面(BDS和GPS雙系統(tǒng))

3 關(guān)鍵技術(shù)

3.1 參數(shù)配置靈活

差分增強系統(tǒng)一般針對特定的衛(wèi)星導(dǎo)航系統(tǒng)向一定區(qū)域內(nèi)用戶提供服務(wù)，大氣條件、接收環(huán)境根據(jù)區(qū)域不同而存在差異。為了滿足用戶深入便捷定制仿真場景的需求，該平臺向用戶提供了文本和圖形化界面兩種方式，允許用戶配置星座、偽距誤差、監(jiān)測站和用戶接收機等參數(shù)。

星座參數(shù)可通過圖形界面和文本兩種方式配置，圖形界面提供衛(wèi)星系統(tǒng)選擇、衛(wèi)星信息查閱功能，文本包含衛(wèi)星的名稱、健康狀況和軌道參數(shù)。該平臺預(yù)置我國的北斗衛(wèi)星導(dǎo)航系統(tǒng)(BeiDou navigation satellite system， BDS)、美國的全球定位系統(tǒng)(global positioning system，GPS)、俄羅斯的格洛納斯衛(wèi)星導(dǎo)航系統(tǒng)(global navigation satellite system， GLONASS)和歐盟的伽利略衛(wèi)星導(dǎo)航系統(tǒng)(Galileo navigation satellite system，Galileo)的星座配置文本。

偽距誤差模型參數(shù)可通過圖形界面配置。星鐘誤差模型采用高斯馬爾科夫模型，可配置參數(shù)為天穩(wěn)定度、秒穩(wěn)定度、初始頻偏、初始鐘偏；星歷誤差模型采用高斯馬爾科夫模型，可配置參數(shù)為相關(guān)時間、采樣間隔和星歷誤差標(biāo)準(zhǔn)差；電離層誤差模型采用Klobuchar模型[5]，可配置參數(shù)分為生成參數(shù)和廣播參數(shù)兩部分，預(yù)置太陽黑子平靜期和活躍期兩組參數(shù)；對流層誤差模型采用WAAS的對流層模型[6]，可配置不同緯度區(qū)域的對流層參數(shù)；接收機鐘差模型采用原子鐘模型[7]，可配置參數(shù)為初始頻偏、初始鐘偏和采樣間隔；熱噪聲及多徑誤差模型采用WAAS的熱噪聲和多徑模型[6]。該平臺以實測數(shù)據(jù)為依據(jù)預(yù)設(shè)各個模型參數(shù)。

監(jiān)測站參數(shù)配置方式分為文本和圖形界面兩種方式，文本僅包含監(jiān)測站的名稱、類別、健康狀況和經(jīng)緯高，圖形界面顯示的監(jiān)測站信息除了文本文件所含內(nèi)容外，還包括監(jiān)測站時鐘初始鐘差和頻差、對流層延遲干分量和濕分量，用戶還可通過直接點擊主界面地圖實現(xiàn)新建、刪除和編輯監(jiān)測站信息。該平臺預(yù)置中國區(qū)域27個監(jiān)測站配置文本。

用戶接收機參數(shù)通過圖形界面配置，可配置參數(shù)包括初始經(jīng)緯高、單頻或雙頻接收機、接收機高度截止角和用戶運動信息，用戶運動類型可分為靜止、勻速直線和勻速圓周運動三類。

3.2 完好性故障仿真

差分增強系統(tǒng)通過播發(fā)完好性信息來提升衛(wèi)星導(dǎo)航系統(tǒng)的可靠性，當(dāng)衛(wèi)星導(dǎo)航系統(tǒng)出現(xiàn)故障或某項運行參數(shù)偏差過大，使系統(tǒng)提供的定位結(jié)果超出規(guī)定限值時，增強系統(tǒng)需及時發(fā)現(xiàn)并及時通告用戶。

為了滿足用戶分析差分增強技術(shù)對完好性風(fēng)險處理能力的需求，該平臺設(shè)計了完好性故障仿真功能。根據(jù)實測數(shù)據(jù)分析可知，星鐘誤差、星歷誤差和電離層誤差是導(dǎo)致完好性風(fēng)險的主要因素，此外接收機鐘差也會導(dǎo)致定位結(jié)果超差。因此將完好性故障分為衛(wèi)星故障、電離層風(fēng)暴和監(jiān)測站故障三類，可單獨設(shè)定每一類故障的起止時間。

(1)

Ik=IKlobuchar+IError

(2)

監(jiān)測站故障分為監(jiān)測站時鐘故障和接收機單通道故障，用戶可對每一個監(jiān)測站單獨設(shè)定時鐘故障etk和接收機通道故障。在監(jiān)測站故障有效期間，由于監(jiān)測站的3臺接收機時鐘共源，時鐘故障疊加到每一臺接收機的觀測量。

tk=tk|k-1+etk

(3)

接收機天線間距離較近，因此接收機單通道故障可由一致性檢驗算法檢出，僅輸出無故障的2臺接收機數(shù)據(jù)。

3.3 平臺可擴展

為了提高該平臺的可擴展性，便于用戶在該平臺基礎(chǔ)上開展后續(xù)研究工作，用戶自定義插件可通過該平臺開放的數(shù)據(jù)接口，獲取衛(wèi)星數(shù)據(jù)、偽距誤差模型數(shù)據(jù)、監(jiān)測站數(shù)據(jù)、用戶接收機數(shù)據(jù)和完好性故障數(shù)據(jù)，可按照RINEX格式和SP3格式讀取該平臺輸出數(shù)據(jù)。程序預(yù)留了GIVE、RAIM、URA、UDER和監(jiān)測站覆蓋深度等5個插件按鈕，用戶可直接點擊調(diào)用相應(yīng)的插件。

以“監(jiān)測站覆蓋深度分析”插件為示例，該插件通過數(shù)據(jù)接口獲取該平臺預(yù)置的GPS星座和27個監(jiān)測站，完成監(jiān)測站覆蓋深度分析，分析結(jié)果如圖4所示。

圖4 監(jiān)測站覆蓋深度插件示例

4 性能測試

該平臺旨在為差分和完好性研究提供規(guī)范逼真的數(shù)據(jù)源，因此需要對其進行性能測試，以確保該平臺輸出數(shù)據(jù)達到接收機實測數(shù)據(jù)的性能，并可滿足差分增強技術(shù)仿真需求。對于用戶接收機端完好性算法而言，并不關(guān)心各誤差源分量的特性，只關(guān)心偽距殘差的整體數(shù)學(xué)特征，并且性能測試以定位誤差為評估標(biāo)準(zhǔn)。GrafNav是美國NovAtel公司開發(fā)的全球衛(wèi)星導(dǎo)航系統(tǒng)(global navigation satellite system， GNSS)后處理軟件，具備高精度定位和長距離實時動態(tài)差分法(real-time kinematic，RTK)解算能力，本節(jié)利用GrafNav 8.40在定位域上對該平臺輸出數(shù)據(jù)進行性能分析。

4.1 靜態(tài)單點定位

為了驗證該平臺輸出數(shù)據(jù)符合標(biāo)準(zhǔn)、分析輸出數(shù)據(jù)的定位精度，進行了靜態(tài)定位測試。以BJFS站在2012-01-01的觀測數(shù)據(jù)為對比，在仿真平臺中添加與BJFS站相同坐標(biāo)的監(jiān)測站BEIJ，仿真時長為1 d，采樣間隔為30 s，均為雙頻接收機。

首先，利用GrafNav進行基本定位，如圖5所示，左側(cè)為仿真數(shù)據(jù)定位結(jié)果，右側(cè)為實測數(shù)據(jù)定位結(jié)果，可以得到如下結(jié)論：

其一，定位誤差相近：仿真數(shù)據(jù)定位落在半徑為6 m的圓內(nèi)(99%)，實測數(shù)據(jù)定位落在半徑為6 m的圓內(nèi)(99%)。

其二，定位分布存在差異：仿真數(shù)據(jù)定位呈條狀，實測數(shù)據(jù)定位呈散點狀。星歷誤差和星鐘誤差是定位的主要誤差源，該平臺采用馬爾科夫模型，前后時隙相關(guān)性較強，難以逼真地模擬真實誤差。

圖5 基本定位對比

隨后，利用GrafNav進行精密單點定位(precise point positioning，PPP )，PPP利用精密星歷消除軌道誤差和星鐘誤差對定位的影響，電離層延遲、對流層延遲、接收機鐘差等因素對定位結(jié)果起主導(dǎo)作用。如圖6所示，左側(cè)為仿真數(shù)據(jù)精密定位結(jié)果，落在半徑為0.5 m的圓內(nèi)(99%)，右側(cè)為實測數(shù)據(jù)精密定位結(jié)果，落在半徑為0.3 m的圓內(nèi)(99%)。由此可知，仿真數(shù)據(jù)定位結(jié)果略差于實測數(shù)據(jù)定位結(jié)果相比，但都顯示出條狀的分布特性。

圖6 精密單點定位對比

通過分析基本定位和精密單點定位的結(jié)果可知，該平臺接收機輸出數(shù)據(jù)規(guī)范有效，在靜態(tài)條件下與實測數(shù)據(jù)的定位性能相當(dāng)。

4.2 動態(tài)差分增強

為了驗證該平臺可支持差分增強技術(shù)研究，進行了動態(tài)差分增強仿真。模擬民航飛機航路區(qū)飛行進行分析，以監(jiān)測站BEIJ為基準(zhǔn)站，飛機起始位置為(115.0°E，39.6°N)，飛行高度為10 000 m，航行速度為東向200 m/s，仿真時長為15 min，采樣間隔為1 s，均為雙頻接收機。圖7為GrafNav解算得到的飛行軌跡，用戶接收機由西向東經(jīng)過BEIJ站上空，地面投影的航跡長約140 km。

圖7 運動軌跡

為了分析差分增強效果，在圖8繪制了實曲線和虛曲線，分別表示差分前和差分后的水平定位誤差。定位誤差(均方根)從2.30 m降到0.11 m，定位精度大幅提高，與實測數(shù)據(jù)性能相當(dāng)[9]，可見該平臺能夠支持差分增強技術(shù)的仿真。用戶可通過布設(shè)多個監(jiān)測站，實現(xiàn)差分增強技術(shù)和精密定位技術(shù)的系統(tǒng)級分析研究。

圖8 水平差分增強效果

5 結(jié)束語

本文設(shè)計了以差分和完好性研究需求為導(dǎo)向的全球衛(wèi)星導(dǎo)航系統(tǒng)綜合仿真平臺，具備參數(shù)靈活配置、數(shù)據(jù)標(biāo)準(zhǔn)化輸出和平臺可擴展的特點，特有完好性故障仿真功能。性能測試結(jié)果表明，該平臺模擬接收機性能達到真實接收機的性能，可以同時支持WAAS、GDGPS、LAAS和CORS等差分增強技術(shù)與精密定位技術(shù)的用戶級和系統(tǒng)級仿真研究。

[1] U.S. Department of Defense. Global Positioning System Standard Positioning Service Performance Standard[EB/OL]. [2013-09-07].http：//www.gps.gov/technical/ps/2008-SPS-performance-standard.pdf.

[2] ZIMMERMANN F， HAAK T， HILL C. Galileo System Simulation Facility[EB/OL]. [2013-09-07].http：//www.gssf.info/Documents/SESP%202004%20-%20Galileo%20System%20Simulation%20Facility%20-%20Issue%201%20-%20Revision%201.pdf .

[3] GURTNER W. RINEX：The Receiver Independent Exchange Format Version 2.10[EB/OL].[2013-09-07].http：//igscb.jpl.nasa.gov/igscb/data/format/rinex210.txt.

[4] HILLA S. The Extended Standard Product 3 Orbit Format (SP3-c)[EB/OL].[2013-09-07].http：//igscb.jpl.nasa.gov/igscb/data/format/sp3c.txt.

[5] KLOBUCHAR J A. Ionospheric Time-delay Algorithm for Single-frequency GPS Users [EB/OL]. [2013-09-07].http：//www.dtic.mil/dtic/tr/fulltext/u2/a229968.pdf.

[6] RTCA. DO-229D， Minimum Operational Performance Standards for Global Positioning System/Wide Area Augmentation System Airborne Equipment [S].

[7] SHIN M Y， PARK C， LEE S J. Atomic Clock Error Modeling for GNSS Software Platform[C]//Proceedings of Position， Location and Navigation Symposium， 2008 IEEE/ION. Monterey， CA：[s.n.]，2008：71-76.

[8] LUO Ming， PULLEN S， DENNIS J， et al. LAAS Ionosphere Spatial Gradient Threat Model and Impact of LGF and Airborne Monitoring[EB/OL]. [2013-09-07].http：//gps.stanford.edu/papers/LuoIONGPS03.pdf.

[9] 黑志堅， 周秋生， 曲建光， 等. GPS RTK 測量成果的精度估計及應(yīng)用探討[J]. 哈爾濱工業(yè)大學(xué)學(xué)報， 2006， 38(8)：1295-1298.