靖守讓，吳鵬，劉文祥，孫廣富

?

基于半長軸補償的廣播歷書參數改進計算方法

靖守讓，吳鵬，劉文祥，孫廣富

(國防科技大學電子科學與工程學院，湖南長沙，410073)

基于全球定位系統(tǒng)(global positioning system, GPS)歷書參數忽略了衛(wèi)星軌道近地點角的變化，無法準確描述中圓地球軌道(medium earth orbit, MEO)衛(wèi)星軌道的長期變化規(guī)律，直接使用歷書參數進行最小二乘擬合的傳統(tǒng)方法無法取得較高精度等問題，針對非歷書長期改正項對位置計算的影響進行分析，提出一種基于半長軸補償的廣播歷書改進計算方法。對擬合歷書增加角速度變化量，并利用其對半長軸平方根進行補償從而獲得廣播歷書參數。研究結果表明：通過此方法擬合6 d衛(wèi)星軌道數據，位置誤差為數千米，用戶測距誤差(user range error,UR)小于900 m，比傳統(tǒng)直接歷書擬合方法精度提高1倍以上；改進法獲得的歷書外推8 d的用戶測距誤差小于3 km，采用傳統(tǒng)方法獲得的歷書用戶測距誤差達8 km。

MEO衛(wèi)星；擴展歷書參數集；半長軸補償；歷書擬合；用戶測距誤差

星歷和歷書參數是導航電文的重要組成部分，使用星歷參數能夠計算衛(wèi)星精密的位置和速度，而歷書參數只能提供衛(wèi)星的粗略位置和速度。前者主要用于導航定位，后者則可用于衛(wèi)星可見性預報、信號傳播時延估計、多普勒頻率偏移估計和星座性能估計 等[1]。導航衛(wèi)星的歷書參數包括精度較低的鐘差參數和星歷參數[2]。本文所討論的歷書參數僅指用于計算衛(wèi)星位置的低精度星歷參數，為了敘述方便，直接將歷書中的低精度星歷參數稱為歷書參數。GPS接口控制文件(interface control document, ICD)[3]僅指出星歷參數是通過最小二乘擬合的方法獲得，而對歷書參數的獲取方法尚未報道。近年來，國內學者對星歷參數擬合進行了研究[4?9]，但對于歷書參數計算的研究較少。崔先強等[10]指出歷書參數計算的2種方法：一種是直接利用歷書參數表達式擬合衛(wèi)星歷書參數，簡稱直接歷書擬合法；另一種是利用星歷參數表達式計算衛(wèi)星星歷參數，然后取出相應的歷書參數，簡稱基于星歷擬合的歷書生成法，同時將Givens變換應用于擬合過程以提高穩(wěn)定性。莫中秋等[11]提出“兩步法”的改進算法，通過提高擬合初值精度的方式提高擬合成功率和擬合精度，并指出歷書擬合結果不應該選取最終收斂時的歷書參數，而應選擇在中間迭代過程中出現最小擬合誤差時的歷書參數。與GPS接口控制文件給出的衛(wèi)星在正常運行下歷書參數的用戶距離誤差(user range error,UR)900 m相比，文獻[10?11]提出的方法并未獲取更好的擬合結果(文中用距離誤差均為1，約為67%)。陳劉成等[12]給出了幾種改進的歷書擬合算法，同時也得到了較好的擬合結果，但幾種改進算法重新設計了歷書參數集并相應修改了歷書算法，擬合得到的歷書參數無法直接應用于導航電文生成。本文通過分析歷書參數忽略的長周期改正項對位置計算的影響，提出包括廣播歷書參數和角速度變化量的擴展擬合歷書參數的概念，并給出擴展歷書參數到廣播歷書的轉換方法。最后將得到的歷書參數與傳統(tǒng)直接歷書擬合得到的歷書參數以及GPS公布的廣播歷書參數精度進行比較，以檢驗此方法獲得歷書的性能。

1 直接歷書擬合方法

1.1 GPS歷書組成

表1所示為GPS廣播星歷參數符號及其含義，其中有上標“*”的參量為歷書參數項(參考時刻的軌道傾角0可以看作是歷書項參考時刻的軌道傾角改正量的等價項)。GPS星歷參數可以分為參考時刻的軌道根數()、長周期項改正參數(,dot,)、短周期改正項振幅(rs,-rc,us,uc,is,ic) 3種不同類型的項[13]，而歷書參數僅包含近似平根數項以及長周期改正參數中的升交點赤經變化率項。

表1 GPS星歷參數

1.2 直接歷書擬合精度

GPS接口控制文件給出了不同運行狀態(tài)時歷書參數的用戶測距誤差(UR)，如表2所示。

表2 GPS歷書EUR指標

統(tǒng)計精度時，在不考慮衛(wèi)星鐘差誤差的情況下，MEO衛(wèi)星軌道UR計算公式為

其中：為徑向軌道誤差；為軌道延跡誤差；為軌道法向誤差。

使用直接歷書擬合方法得到GPS衛(wèi)星1歷書精度如圖1所示。擬合衛(wèi)星軌道使用IGS(international GNSS service)發(fā)布的2011?12?25—12?31的精密軌道數據。

1—用戶測距誤差；2—位置誤差

從圖1可以看出：使用直接歷書擬合位置誤差達數萬米，UR達數千米。

2 基于半長軸補償的改進方法

與星歷參數相比，歷書參數忽略了角速度變化量?、軌道傾角變化率dot和短周期改正項。對于使用長弧段衛(wèi)星軌道的歷書擬合而言，短周期改正項影響可以忽略，因此，主要分析?和dot對位置計算的影響。

2.1 非歷書長周期項影響

?在擬合過程中吸收了0和的長期項和長周期項，主要是二階帶諧項引起的的長期漂移，也包括了日、月引力攝動力和太陽光壓攝動力，從物理規(guī)律上角速度變化量?應該是常值或長期變化值。軌道傾角在幾個主要攝動力作用下無長期項，軌道傾角變化率dot僅吸收了軌道傾角的長周期項[14]。通過逐步增加擬合步長的方法發(fā)現擬合得到的星歷參數?為10?9數量級，擬合得到dot比?小1個數量級。不同擬合弧段?和dot擬合結果見圖2。

(a) ?擬合結果；(b)dot擬合結果

1—衛(wèi)星1；2—衛(wèi)星2

圖2 不同擬合弧段?和dot擬合結果

Fig.2 Fitting result of ?anddotby different curve fit interval

(3)

式(2)和(3)所示為衛(wèi)星位置對?和dot的偏導數。其中：為長半軸長度；t為計算時刻與參考時刻之差；E為偏近點角；r為修正的半徑；u為修正的升交角距；i為修正的軌道傾角；為修正的升交點經度。

由于偏心率是1個較小量，其影響可忽略不計，因此，?對衛(wèi)星位置的影響主要體現在式(2)的第2項上。對導航衛(wèi)星而言，r為107數量級，使用軌道數據進行擬合通常?在10?9數量級；對歷書參數而言，t可能達到105數量級，因此，?項對衛(wèi)星位置影響可達數萬米；對軌道傾角變化率dot項，擬合結果一般在10?10數量級，其對衛(wèi)星位置的影響要比?小1個數量級，其對衛(wèi)星位置計算影響僅數千米。

為了驗證上述定性分析的正確性，通過將擬合得到的星歷參數中元素置零的方法驗證其對衛(wèi)星位置計算的影響，驗證結果如圖3所示。

1—設置?n為0 rad/s；2—設置idot為0 rad/s

經對比發(fā)現：將?置零后，對使用星歷參數計算衛(wèi)星位置影響達到數萬米，而dot影響比?影響小1個數量級。

2.2 半長軸補償方法及補償帶來的精度降低

為了得到較好的歷書擬合結果，需要考慮近地點角和平近點角的長期變化，即在擬合過程中考慮角速度變化量?。本文提出一種基于半長軸補償的廣播歷書計算方法，即使用一套包括廣播歷書和角速度變化量的擴展歷書參數()進行最小二乘擬合，然后根據?與的關系將擴展歷書參數轉換為廣播歷書參數。該方法簡稱為改進方法。

利用使用星歷用戶算法計算衛(wèi)星位置時，?項僅用于計算校正平均角速度，滿足

(5)

進一步分析使用式(5)進行半長軸補償對衛(wèi)星位置計算的影響。由長半軸與平均角速度的關系，可以推導出長半軸對平均角速度的導數為，而對于GPS等導航系統(tǒng)中的MEO衛(wèi)星，平均角速度約為地球自轉速度的2倍，長半軸對平均角速度的導數約為?1.2×1011，而典型的?為10?9數量級，此將?對半長軸的補償值僅為100 m數量級。結合計算衛(wèi)星位置用戶算法可知，在后續(xù)計算中半長軸也僅在計算改正后的軌道半徑時使用。因此，補償操作對衛(wèi)星位置計算的影響也在100 m數量級，遠低于忽略角速度變化量對衛(wèi)星位置計算的影響。

3 算例分析及性能驗證

仿真驗證數據使用IGS發(fā)布的2011?12?25—2012?01?07共14 d的精密星歷數據以及celestrak網站發(fā)布的參考時刻為644周，周內時間為233 472 s的GPS系統(tǒng)Yuma格式歷書數據。

GPS ICD指出歷書參數由控制站至少每6 d更新1次，并指出Block II和Block IIA衛(wèi)星使用3組歷書保證60 d跨度使用，且3組歷書均基于6 d軌道數據擬合；Block IIR/IIR-M，Block IIF和Block III衛(wèi)星使用5組歷書，前3組基于6 d數據擬合，另外2組基于32 d數據擬合計算。因此，本文仿真將前6 d精密星歷作為擬合弧段計算參考周內秒為259 200 s的歷書，其余8 d數據作為歷書外推精度的參考標準。由于具體擬合弧段未知，對GPS歷書而言，擬合弧段和外推弧段僅指示精密星歷對應時間段。

圖4所示為使用傳統(tǒng)直接歷書擬合法和改進法獲得GPS衛(wèi)星1歷書參數在擬合弧段和外推弧段的位置精度。擬合弧段和外推弧段直接擬合法得到歷書、改進法得到歷書以及GPS歷書的UR如表3所示。

從圖4和表3可以看出：

1)在擬合弧段使用改進法得到歷書精度明顯高于直接擬合法得到的歷書精度，精度提高1倍以上。

2)改進法在擬合弧段內位置誤差均在數千米以內，URE小于900 m，與GPS歷書精度相當。

3) 在外推弧段，直接歷書擬合得到的歷書參數發(fā)散較快，而改進法獲得歷書發(fā)散較小。

4) 改進法獲得的歷書外推8 d，UR統(tǒng)計結果小于3 km，與GPS歷書精度相當。

誤差比較：(a) X軸；(b) Y軸；(c) Z軸；(d) 三維

表3 3種歷書EUR精度統(tǒng)計結果

4 結論

1)在進行歷書擬合中必須考慮近地點角和平近點角的長期變化即角速度變化量。

2)基于半長軸補償的廣播歷書計算方法在位置精度和外推性能上明顯優(yōu)于傳統(tǒng)直接歷書擬合法。擬合6 d GPS軌道得到的歷書參數與GPS發(fā)布歷書參數精度相當，平均位置誤差為1 km數量級，UR滿足GPS正常運行狀態(tài)指標要求，并且使用改進法獲得的歷書參數具有良好的外推性能。在GPS歷書生成方法未公開的情況下，此方法可以作為一種替代方法生成廣播歷書，為信號源或新組建導航系統(tǒng)的廣播歷書參數生成提供參考。

參考文獻：

[1] 王夢麗, 陳華明, 王飛雪. GPS歷書數據的有效齡期[J]. 遙測遙控, 2007, 28(3): 31?35. WANG Mengli, CHEN Huaming, WANG Feixue. Effective age of GPS almanac data[J]. Journal of Telemetry, Tracking and Command, 2007, 28(3): 31?35.

[2] Kaplan E D, Hegarty C J. Understanding GPS: Principles and applications[M]. 2nd ed. Norwood, MA: Artech House, 2006: 67?69.

[3] Revision H. Navstar GPS space segment/navigation user segment interfaces, global positioning systems directorate systems engineering & integration[EB/OL]. [2013?09?24]. http://www.gps.gov/technical/icwg.

[4] 謝小剛, 曾大治, 龍騰, 等. 北斗GEO衛(wèi)星廣播星歷的直接擬合算法研究[J]. 華中科技大學學報(自然科學版), 2014, 42(3): 14?18. XIE Xiaogang, ZENG Dazhi, LONG Teng, et al. Research on broadcast ephemeris direct fitting algorithm of BeiDou GEO satellite[J]. Journal of Huazhong University of Science & Technology (Natural Science Edition), 2014, 42(3): 14?18.

[5] 黃勇, 胡小工, 王小亞, 等. 中高軌衛(wèi)星廣播星歷精度分析[J]. 天文學進展, 2006, 24(1): 81?88. HUANG Yong, HU Xiaogong, WANG Xiaoya, et al. Precision analysis of broadcast ephemeris for medium and high orbit satellites[J]. Progress in Astronomy, 2006, 24(1): 81?88.

[6] 劉光明, 廖瑛, 文援蘭, 等. 導航衛(wèi)星廣播星歷參數擬合算法研究[J]. 國防科技大學學學報, 2008, 30(3): 100?104. LIU Guangming, LIAO Ying, WEN Yuanlan, et al. Research on the fitting algorithm of broadcast ephemeris parameters[J]. Journal of National University of Defense Technology, 2008, 30(3): 100?104.

[7] 陳劉成, 李靜, 馬瑞, 等. 工程化廣播星歷參數擬合算法與接口設計[J]. 武漢大學學報(信息科學版), 2011, 36(1): 18?23. CHEN Liucheng, LI Jing, MA Rui, et al. The engineering design for broadcast ephemeris parameters fitting arithmetic and their interfaces[J]. Geomatics and Information Science of Wuhan University, 2011, 36(1): 18?23.

[8] 陳忠貴, 劉光明, 廖瑛, 等. 廣播星歷參數星上自主擬合算法[J]. 國防科技大學學報, 2011, 33(3): 1?4. CHEN Zhonggui, LIU Guangming, LIAO Ying, et al. Autonomously updated broadcast ephemeris algorithm[J]. Journal of National University of Defense Technology, 2011, 33(3): 1?4.

[9] 何峰, 王剛, 劉利, 等. 地球靜止軌道衛(wèi)星廣播星歷參數擬合與試驗分析[J]. 測繪學報, 2011, 40(增刊): 52?58. HE Feng, WANG Gang, LIU Li, et al. Ephemeris fitting and experiments analysis of GEO satellite[J]. Acta Geodaetica et Cartographica Sinica, 2011, 40(Suppl): 52?58.

[10] 崔先強, 焦海文, 秦顯平. 導航衛(wèi)星的歷書參數及其擬合算法[J]. 測繪科學技術學報, 2006, 23(1): 23?25. CUI Xianqiang, JIAO Haiwen, QIN Xianping. Ephemeris parameters of navigation satellite and their fitting arithmetic[J]. Journal of Geomatics Science and Technology, 2006, 23(1): 23?25.

[11] 莫中秋, 陳劉成, 董恩強, 等. 導航衛(wèi)星歷書擬合初值的選取與誤差分析[J]. 測繪科學技術學報, 2008, 25(2): 104?107. MO Zhouqiu, CHEN Liucheng, DONG Enqiang, et al. Improvement on initial values selection and errors analysis of almanac fitting of navigation satellites[J]. Journal of Geomatics Science and Technology, 2008, 25(2): 104?107.

[12] 陳劉成, 胡小工, 韓春好, 等. 導航衛(wèi)星歷書參數擬合算法研究[J]. 天文學報, 2008, 49(3): 288?296. CHEN Liucheng, HU Xiaogong, HAN Chunhao, et al. Research on almanac parameters fitting algorithm for navigation satellites[J]. Acta Astronomica Sinica, 2008, 49(3): 288?296.

[13] 黃華. 導航衛(wèi)星廣播星歷參數模型及擬合算法研究[D]. 南京: 南京大學天體測量與天體力學, 2012: 31?44. HUANG Hua. Research on the Broadcast ephemeris parameters model and its fitting Algorithm[D]. Nanjing: Nanjing University. Astrometry and Celestial Mechanics, 2012: 31?44.

[14] 徐紹銓, 張華海, 楊志強, 等. GPS測量原理及應用[M]. 3版. 武漢: 武漢大學出版社, 2008: 32?41. XU Shaoquan, ZHANG Huahai, YANG Zhiqiang, et al. GPS measurement principle and application[M]. 3rd ed. Wuhan: Wuhan University Press, 2008: 32?41.

Fitting method of improved almanac parameters with semi-major axis compensation

JING Shourang, WU Peng, LIU Wenxiang, SUN Guangfu

(College of Electronic Science and Engineering, National University of Defense Technology, Changsha 410073, China)

Considering that the medium earth orbit (MEO) satellite orbit is ineffectively described by the global positioning system (GPS) almanac without considering the change of the argument of perigee, the traditional algorithm using almanac parameters can’t directly obtain good effect via the least squares curve fit, an almanac calculation method based on semi-major axis compensation was proposed by analyzing the effect of elements for calculating satellite position, which was neglected by almanac parameters. An extended fitting almanac considering mean motion difference was provided. The broadcast almanac parameters were acquired by compensating the effect of mean motion difference on semi-major axis. The simulation results show that the satellite position error of almanac generated by the method is just several kilometers in 6 d, and the user range error is less than 900 m, whose precision doubles than that of the method fitted using almanac parameters.URof almanac generated by the proposed method is less than 3 km in 8 d extrapolation, while that of the traditional method is more than 8 km.

medium earth orbit (MEO); extended almanac parameters set; semi-major axis compensation; almanac fitting; user range error (URE)

10.11817/j.issn.1672-7207.2015.07.017

TN967.1

A

1672?7207(2015)07?2504?06

2014?07?11；

2014?09?22

國家自然科學基金資助項目(61403413) (Project(61403413) supported by the National Natural Science Foundation of China)

靖守讓，博士研究生，從事衛(wèi)星導航信息處理算法研究；E-mail: hanchongjsr@163.com

(編輯 陳燦華)