GPS衛(wèi)星星歷仿真

崔立魯

(成都大學 建筑與土木工程學院，四川 成都 610106)

0 引 言

隨著衛(wèi)星導航系統(tǒng)及相關應用越來越廣泛，人們對衛(wèi)星導航定位的要求也越來越高.一個可靠性高、適用性強、計算精度高的數(shù)據(jù)處理軟件直接決定了導航定位的解算精度.而檢驗數(shù)據(jù)處理軟件最有效的方法是利用衛(wèi)星仿真數(shù)據(jù)進行驗算，其中衛(wèi)星星歷的仿真極為重要.目前，國外已有相關研究機構開展了這方面的研究［1－2］，而國內也有一些科研機構開始對GPS 仿真進行研究［3－5］，但相關GPS 仿真方法的研究和文獻還比較少.鑒于此，本研究簡要地闡述了GPS 衛(wèi)星廣播星歷的仿真原理，利用衛(wèi)星參數(shù)數(shù)據(jù)對該衛(wèi)星廣播星歷數(shù)據(jù)實現(xiàn)仿真模擬，并與真實的廣播星歷進行比較，驗證了該仿真原理的可靠性和適用性，實現(xiàn)了衛(wèi)星仿真的關鍵性步驟.

1 GPS 廣播星歷仿真

1.1 廣播星歷參數(shù)

GPS 廣播星歷由一系列參數(shù)組成，這些參數(shù)以開普勒軌道參數(shù)為基礎，再加上一些表示軌道蛇洞的調和系數(shù)及隨時間的變率，共有16 個:一個參考時刻，6 個相應參考時刻的開普勒軌道參數(shù)，9 個反映攝動力影響的參數(shù).這些變率在仿真廣播星歷時以RINEX 格式輸出，具體參數(shù)如表1 所示.

表1 中，各參數(shù)的具體含義為:TOC 為時間參數(shù)，參考時刻的年、月、日、時、分、秒;a0為衛(wèi)星鐘誤差的常數(shù)項，單位s;a1為衛(wèi)星鐘誤差的變化速率，單位s/s;a2為a1的變化速率，單位s/s2;IODE 為星歷數(shù)據(jù)的年齡，單位s;Crs為在軌徑向方向上周期改正正弦項的振幅，單位m;Δn 為平均角速度的長期變化率，單位rad/s;M0為參考時刻的平近點角，單位rad;Cuc為在軌道延跡方向上周期改正余弦項的振幅，單位rad;e 為衛(wèi)星軌道的離心率;Cus為在軌道延跡方向上周期改正正弦項的振幅，單位rad為衛(wèi)星軌道長半軸的平方根，單位In V2;t0為星歷參考時刻所在星期中的秒數(shù)，單位s;Cic為軌道傾角周期改正余弦項的振幅，單位rad;Ω0為參考時刻升交點赤徑的主項，單位rad;Cis為軌道傾角周期改正正弦項的振幅，單位rad;i0參考時刻軌道傾角，單位rad;Crc為在軌道徑向方向上周期改正正弦項的振幅，單位m;ω 為近地點角距，單位rad;Ω·為升交點赤徑在赤道平面中的長期變化，單位rad/s;i 為軌道傾角變化率，單位rad/s;IODC 為衛(wèi)星鐘數(shù)據(jù)年齡，單位s.

表1 廣播星歷軌道及各軌道參數(shù)

1.2 星歷參數(shù)仿真模型

1.2.1 時間參數(shù)生成模型.

星歷文件中的時間參數(shù)包括GPST 周計數(shù)weekno、獲得星歷參數(shù)時刻t0、星歷數(shù)據(jù)年齡IODE、衛(wèi)星鐘數(shù)據(jù)年齡IODC 等.其中，weekno 和t0均需要根據(jù)UTC 來生成.令UTC 時間為tUTC，它與GPST 起始時刻所對應的起始時刻之間相隔的星期數(shù)為WUTC，同時考慮模擬所在時刻的時區(qū)為Tzone，則，

其中，n 為UTC 與GPST 之間的跳秒數(shù)，toc是衛(wèi)星導航電文中第一子幀的參考時刻，t1是計算時刻改正數(shù)所作測量的最后觀測時刻.由于在本研究的仿真過程中主要基于軟件方法的實現(xiàn)，所以后2 個參數(shù)均置零.

1.2.2 時鐘改正數(shù)生成模型.

本研究在仿真計算時，對星歷文件中的3 個時鐘改正參數(shù)，鐘偏a0、鐘漂a1、鐘頻a2，按照下述關系進行計算，設定參考時刻的3 個參數(shù)為，則有，

在本研究中，這3 個參數(shù)設置為零.

1.2.3 軌道參數(shù)的生成.

根據(jù)求得的慣性坐標系下衛(wèi)星的位置和速度求解衛(wèi)星初始時刻的6 個開普勒軌道根數(shù).

1)根數(shù)a，e，M 的計算.

其中，

2)3 個定向i，Ω，ω 根數(shù)的計算.

Pz= sinisinω，

Qz= sinicosω，

Rx= sinisinΩ，

Ry= sinicosΩ，

Rz= cosi

根據(jù)上面5 個式子可解出，

2 仿真數(shù)據(jù)計算和分析

2.1 計算方法

開普勒軌道根數(shù)擬合是一迭代收斂的過程，初值的選取對迭代結果具有一定的影響.首先計算參考時刻衛(wèi)星對應的慣性系下的位置速度，然后以二體問題的方法求解初始6 個開普勒軌道根數(shù)，并將軌道根數(shù)轉化到需要的地固坐標系中作為初值，將其余9 個攝動調和參數(shù)初值設為0，其具體計算方法如下:①根據(jù)讀取的精密星歷插值求出任意時刻衛(wèi)星的位置和速度;②將地固系下的toe時刻衛(wèi)星狀態(tài)矢量，通過極移、歲差和章動逆向變換到J2000 坐標系，求得在J2000 下的開普勒軌道根數(shù);③將升交點赤經Ω 與參考時刻格林尼治恒星時CAST 組差，得到地固系下的初始軌道根數(shù);④由衛(wèi)星在慣性系下的坐標與速度計算衛(wèi)星在初始時刻的6 個開普勒軌道根數(shù);⑤計算出的6 個開普勒軌道根數(shù)及9 個攝動項參數(shù)(初值設為0)，計算衛(wèi)星的近似坐標及誤差方程系數(shù)陣(通過求導得到)，并列出誤差方程;⑥由最小二乘法迭代求解15 個參數(shù).

2.2 試驗數(shù)據(jù)比較

本研究計算了3 號衛(wèi)星從1月31日開始到2月12日連續(xù)13 d 的每歷元時刻的瞬時軌道參數(shù).表2為3 號衛(wèi)星第一組平差計算所得第一個歷元時刻的瞬時軌道參數(shù)，表3 為3 號衛(wèi)星廣播星歷中第一個歷元時刻的瞬時軌道參數(shù)，從這兩個表中數(shù)據(jù)對比可以看出求解的6 個開普勒軌道根數(shù)與9 個攝動項參數(shù)的差異.

表2 3 號衛(wèi)星第一組平差計算所得第一歷元時刻的瞬時軌道參數(shù)

表3 3 號衛(wèi)星廣播星歷中第一歷元時刻的瞬時軌道參數(shù)

從求出的第一歷元的廣播星歷參數(shù)與實際接收機所接收到的廣播星歷數(shù)據(jù)比較可以看出，各參數(shù)之間差異很小，對個別攝動參數(shù)差別較大的數(shù)據(jù)進行替換計算，其對衛(wèi)星位置的影響在分米級以內.

表4 是根據(jù)求出的參數(shù)計算出衛(wèi)星位置與精密星歷中衛(wèi)星位置求差所得的數(shù)值.

表4 精密星歷與仿真星歷計算出的衛(wèi)星位置差值

表4 數(shù)據(jù)表明，按上述方法求解共15 個參數(shù)，由9 個歷元計算得到的軌道曲線基本上通過每個歷元的衛(wèi)星位置，計算得到的軌道與衛(wèi)星真實軌道基本一致.

圖1、圖2 給出了仿真星歷的3 個參數(shù)與接收機接收到的廣播星歷中對應的3 個參數(shù)的差值，這2顆衛(wèi)星分別為PRN3、PRN8.圖1、圖2 中，縱坐標的單位為rad，時間區(qū)間為1月31日到4月12日，數(shù)據(jù)對應的時刻為每日的02:00 時.

圖1 PRN3 模擬的3 個參數(shù)與實際參數(shù)差值曲線圖

圖2 PRN8 模擬的3 個參數(shù)與實際參數(shù)差值曲線圖

從圖1、2 可知，仿真星歷的3 個參數(shù)M0、Ω、ω 值與接收機接收到的廣播星歷中的M0、Ω、ω 值之差在一定范圍內波動，但在所列圖中個歷元參數(shù)計算出的衛(wèi)星位置與精密星歷中衛(wèi)星位置差值都在厘米級.由此可知，仿真星歷能夠滿足檢驗數(shù)據(jù)軟件算法的需要.

3 結 論

GPS 衛(wèi)星星歷仿真是GPS 衛(wèi)星仿真中比較重要的一個環(huán)節(jié)，也是驗證衛(wèi)星數(shù)據(jù)處理軟件的重要依據(jù).本研究主要介紹了GPS 衛(wèi)星廣播星歷仿真的基本原理，并將仿真數(shù)據(jù)與真實數(shù)據(jù)進行了比較，利用仿真數(shù)據(jù)計算出相應的衛(wèi)星位置數(shù)據(jù)并與真實衛(wèi)星位置數(shù)據(jù)進行了比較.數(shù)據(jù)比較的結果表明:利用本研究提出的衛(wèi)星星歷仿真方法計算得到的衛(wèi)星位置數(shù)據(jù)能夠較好地反應真實衛(wèi)星的情況，有助于衛(wèi)星數(shù)據(jù)處理軟件的驗證.

［1］Kouba J，Héroux P.Precise point positioning using IGS orbit and lock products［J］.GPS Sol，2002，5(2):101－106

［2］許其鳳.GPS 衛(wèi)星導航與精密定軌［M］.北京:解放軍出版社，1994.

［3］崔立魯，何秀鳳，羅志才，等.Galileo 系統(tǒng)定位性能仿真模擬分析［J］.測繪信息與工程，2007，32(4):10－12.

［4］崔立魯，羅志才，鐘波，等.觀測噪聲對頻域輸入輸出法數(shù)據(jù)融合的影響［J］.大地測量與地球動力學，2009，29(1):79－82.

［5］崔立魯.重力衛(wèi)星加速度數(shù)據(jù)的插值算法研究［J］.測繪工程，2012，21(3):21－25.