周 雨，李英冰

（武漢大學測繪學院，湖北 武漢430079）

0 引 言

仿真GPS觀測數(shù)據(jù)的系統(tǒng)的誤差和隨機誤差都是可以人為控制的，使用的各種模型參數(shù)都是精確已知的，使得GPS仿真數(shù)據(jù)在測試GPS數(shù)據(jù)解算軟件的正確性、基于GPS觀測值的定軌方案論證分析、評論各種GPS算法性能等方面，都比實際觀測數(shù)據(jù)更有優(yōu)勢［1－3］。

GPS的仿真數(shù)據(jù)一般由模擬生成器或仿真軟件實現(xiàn)，如GPSLab、瑞士的伯爾尼大學研制的精密GPS定位軟件中的GPS仿真模塊等［4－6］。但GPS仿真生成器的價格較高，又無法獲得其他國外仿真軟件的源代碼［4］，因此GPS仿真對我們研究衛(wèi)星導航，如北斗衛(wèi)星導航定位系統(tǒng)等，都有著非常重要的現(xiàn)實意義。

從單點定位的原理出發(fā)，考慮電離層延遲、對流層延遲、衛(wèi)星鐘差、接收機鐘差、多路徑效應、相對論效應、地球自轉(zhuǎn)改正、測量噪聲等影響因素［6－8］，自主設計面向用戶的GPS觀測值仿真軟件，模擬偽距和載波的生成，并以IGS跟蹤站的實際觀測數(shù)據(jù)為參考，對模擬值與真實值進行定量比較，并定性的分析出現(xiàn)誤差的原因，進一步明確模擬觀測的可行性和穩(wěn)定性。

1 仿真觀測值的生成原理

GPS觀測值的仿真，是利用用戶任意給出的測站三維坐標來模擬該測站觀測的GPS信號。GPS衛(wèi)星發(fā)射的信號由載波、測距碼和導航電文三部分組成，主要對偽距和載波進行模擬，對他們建立如下的觀測方程［7－8，3］。

式中：i為衛(wèi)星編號；t為時刻；ρi為接收機至觀測碼偽距；XiYiZi為衛(wèi)星的WGS－84坐標；XYZ為接收機的WGS－84坐標；c為光速；VtR為接收機鐘差；為衛(wèi)星鐘差；為電離層延遲改正；為對流層延遲改正；δρi為地球自轉(zhuǎn)改正；δρimul為多路徑效應改正；ei，εi分別為觀測噪聲；λ為相應載波的波長；φi為載波相位觀測值；Ni為整周模糊度。

在式（1）、（2）中，衛(wèi)星坐標XiYiZi和衛(wèi)星鐘差Vits可以根據(jù)IGS發(fā)布的精密星歷進行拉格朗日插值得到，接收機坐標XYZ可以由模擬條件直接給出，電離層延遲改正Viion、對流層延遲改正Vitrop、地球自轉(zhuǎn)改正δρi、多路徑效應δρimul可以采用相應模型改正，測量噪聲eiεi可以通過隨機模型直接生成。

2 仿真觀測值軟件的設計

GPS觀測值仿真的具體過程為：首先，根據(jù)模擬的觀測時段讀取對應的IGS發(fā)布的間隔15min的精密星歷，并將其進行拉格朗日插值得到與模擬采樣率間隔一致的精密星歷；其次，從該插值的精密星歷文件中獲取衛(wèi)星坐標、衛(wèi)星速度和衛(wèi)星鐘差，根據(jù)衛(wèi)星坐標與模擬的測站坐標判斷衛(wèi)星的可視情況，并計算出測站到衛(wèi)星的幾何距離；復次，根據(jù)模型公式對電離層延遲和對流層延遲進行計算，按照相應模型公式計算地球自轉(zhuǎn)延遲、相對論效應延遲、多路徑效應延遲，并對接收機鐘差和測量噪聲進行隨機模擬；再次，將各項延遲改正按照用戶需要加入測站到衛(wèi)星的幾何距離中，得到偽距的仿真值，再加入相位觀測的初始模糊度和隨機周跳值，得到模擬的載波相位觀測值；最后將數(shù)據(jù)按照Rinex文件格式進行整理，生成用戶需求的觀測值文件。其具體流程圖如圖1所示，程序的運行的界面如圖2所示。

圖1 程序的流程圖

3 仿真實驗與結(jié)果

圖2 軟件運行截圖

利用上述的仿真觀測軟件，選擇現(xiàn)有的IGS全球跟蹤站作為參考位置，下載實際觀測值與仿真觀測值進行對比研究。僅選取了上海站（shao）的模擬觀測結(jié)果與真實觀測結(jié)果進行對比分析。在整個仿真過程中，參數(shù)設置依次為：觀測時段為2012年2月10日0時0分至23時59分，采樣率為30s，截至高度角為5°，電離層改正采用Klobuchar模型，對流層改正選用Hopfield模型，多路徑效應改正、接收機鐘差、地球自轉(zhuǎn)改正、相對論效應改正、測量噪聲均采用相應的模型。

從圖3可以看出虛擬觀測到的衛(wèi)星與實際觀測的衛(wèi)星基本一致，由于虛擬觀測設定的截止高度角為5°，與上海站（shao）的實際觀測情況有出入，所以出現(xiàn)了部分時段虛擬觀測的衛(wèi)星與實際觀測的衛(wèi)星數(shù)不符。圖4示出了2號衛(wèi)星的真實C1、L1、L2值與虛擬C1、L1、L2值的差值，在200個觀測歷元中，C1差值變化不大且都在0到10 m的范圍中，L1和L2差值都在－30周至30周的范圍內(nèi)，基本滿足GPS單點定位的精度要求。圖5示出了對所有衛(wèi)星的真實C1值與虛擬C1值的差值統(tǒng)計情況，差值基本出現(xiàn)在0到10m的范圍內(nèi)，仿真的觀測值與實際觀測值相差不大，主要誤差可能是模型改正不精確造成的，還有其他的一些誤差可能來自實際觀測的外部環(huán)境、IGS精密星歷文件。

圖3 可視衛(wèi)星總數(shù)對比

4 結(jié) 論

基于GPS觀測值的生成原理，采用面向?qū)ο缶幊谭椒?，開發(fā)了GPS觀測值的仿真軟件，該軟件經(jīng)過初步的試用，得到了較高精度的GPS仿真觀測值，基本滿足單點定位的精度要求，可以為其他相關(guān)科學的研究與應用提供參考。

隨著計算機硬件和軟件技術(shù)的發(fā)展，仿真的成本會越來越低，仿真的結(jié)果會越來越真實，仿真系統(tǒng)在GNSS導航的各個領域都將有更加廣泛的應用。

［1］ 匡翠林，鄒 璇，李 陶.利用IGS數(shù)據(jù)產(chǎn)品進行高精度GPS觀測數(shù)據(jù)仿真［J］.系統(tǒng)仿真學報，2007，19（12）：2857－2859.

［2］ BARTLETT S，LEGRAND R，STATION，E.A－nalysis and simulation of discrete Kalman filtering applied to a GPS receiver［J］.Naval Engineers Journal，1996，108（1）：29－36.

［3］ ELSA M，MAURICIO G，CLAUDIO B.SiGOG：simulated GPS observation generator［J］.GPS Solution，2005（9）：250－254.

［4］ 范國清，王 威，郗曉寧.高精度GPS觀測數(shù)據(jù)的實時仿真研究［J］.國防科技大學學報，31（3）：60－64.

［5］ 郝 明，崔篤信，王慶良，等.高動態(tài)GPS觀測數(shù)據(jù)的仿真研究［J］.測繪科學，2008，33（6）：60－62.

［6］ 朱偉剛，陶春燕.GPS觀測數(shù)據(jù)的建模與仿真［J］.海洋測繪，2007，27（5）：26－29.

［7］ 裴 霄，王解先.GPS偽距單點定位的精度分析及改進［J］.海洋測繪，2012，32（1）：5－7.

［8］ 廖 華.GPS偽距單點定位算法的綜合比較［J］.測繪科學，2011，36（1）：20－21.