王 樂 賈小林

1)長安大學地質(zhì)工程與測繪學院，西安 710064

2)西安測繪研究所，西安 710054

北斗廣播星歷參數(shù)擬合及精度分析*

王 樂1)賈小林2)

1)長安大學地質(zhì)工程與測繪學院，西安 710064

2)西安測繪研究所，西安 710054

根據(jù)GPS ICD文檔中定義的一組18參數(shù)廣播星歷模型及其用戶算法推導了三維坐標對各參數(shù)的偏導數(shù)表達式，建立了18參數(shù)和衛(wèi)星位置矢量之間的線性函數(shù)模型，利用最小二乘平差原理計算北斗衛(wèi)星18參數(shù)廣播星歷，分別通過北斗衛(wèi)星定軌數(shù)據(jù)對算法進行測試，結果驗證了算法的正確性，同時證明在4小時之內(nèi)，擬合誤差優(yōu)于0.01 m。

北斗衛(wèi)星;用戶算法;擬合分析;GPS ICD文檔;18參數(shù)廣播星歷模型

1 引言

眾所周知，獲取衛(wèi)星的位置和速度是進行導航定位的前提條件，而衛(wèi)星的位置和速度又是由衛(wèi)星的廣播星歷參數(shù)來表征的。因此，研究導航衛(wèi)星廣播星歷參數(shù)的生成算法，對于提高衛(wèi)星導航定位的精度具有非常重要的應用價值。

通常所用的廣播星歷模型是以美國的GPS為代表的16參數(shù)模型，已經(jīng)有很多學者對GPS 16參數(shù)模型的擬合算法的效率和精度進行了深入研究［1－3］，然而在最新發(fā)布的 GPS ICD 文檔中定義了一組最新的18參數(shù)廣播星歷，對廣播星歷18參數(shù)擬合算法進行分析研究和實際數(shù)據(jù)測試，能夠檢驗這種新型星歷模型的有效性。

2 北斗導航衛(wèi)星廣播星歷18參數(shù)擬合算法

18 參數(shù)廣播星歷以開普勒軌道根數(shù)加攝動參數(shù)的形式給出，與GPS廣播星歷16參數(shù)比較起來在長半軸、平均角速度、升交點赤經(jīng)的計算方法存在差異。

2.1 18參數(shù)廣播星歷偏導數(shù)的推導

利用18 參數(shù)的用戶算法［4，5］計算tk時刻地面點在地固直角坐標系中的坐標為:

其余16個參數(shù)的偏導數(shù)表達式可參考文獻［5］。上述過程表明在推導廣播星歷18參數(shù)偏導數(shù)的過程中忽略了調(diào)和項對偏導數(shù)的影響，經(jīng)算例驗證這種忽略調(diào)和項的算法是切實可行的。

2.2 廣播星歷18參數(shù)的最小二乘擬合

根據(jù)已知的廣播星歷用戶算法，衛(wèi)星坐標與星歷參數(shù)可表示成

將式(12)在toe時刻處進行泰勒級數(shù)展開，舍掉二階和二階以上的小量后得到的線性化表達式為:

設現(xiàn)在有K(K≥6)個歷元，則誤差方程為

根據(jù)最小二乘原理有:

為了提高參數(shù)擬合的精度，通常需進行循環(huán)迭代，迭代時將上一次計算結果作為下一次迭代的初值，則第i次計算結果為:

3 算例及分析

采用北斗衛(wèi)星 2012－09－03—06日的采樣間隔為5分鐘的精密軌道數(shù)據(jù)(三維坐標值和三維速度值)分別對GEO、MEO和IGSO衛(wèi)星的18參數(shù)廣播星歷進行擬合。

該組精密軌道數(shù)據(jù)是利用長春、烏魯木齊、喀什、臨潼和昆明5個國內(nèi)跟蹤站以及8個境外站的接收機，連續(xù)對各階段入網(wǎng)的11顆衛(wèi)星進行觀測，獲取偽距、相位數(shù)據(jù)，通過事后精密定軌得到的在軌衛(wèi)星精密星歷，通過重疊弧段統(tǒng)計其精度為:徑向平均誤差為0.46 m，切向平均誤差為2.62 m，法向平均誤差為1.39 m。

為了分析擬合精度，我們分別統(tǒng)計了用6個至61個已知點擬合17參數(shù)的精度值，本文采用的精度指標為用戶測距誤差(URE):

圖1～3分別表示MEO、GEO、IGSO衛(wèi)星18參數(shù)擬合精度隨擬合點數(shù)增加的變化情況。從圖中可以明顯看出三種衛(wèi)星整體擬合URE隨著擬合已知點的增加而增大，但是當已知點個數(shù)在50個(約4小時)之內(nèi)時擬合誤差增速較慢，而已知點個數(shù)大于50個(約4小時)時，擬合精度發(fā)散速度急劇增大。同時在已知點個數(shù)小于50個時，對于GEO和IGSO衛(wèi)星的參數(shù)擬合精度變化呈現(xiàn)出約1小時的周期性。

表1～3分別羅列了MEO、GEO、IGSO前10組不同時長的擬合結果。從表中可以看出，三種衛(wèi)星的參數(shù)擬合精度在擬合時長為2小時和3小時時最為均勻，保持在3 mm左右，隨著擬合時長的增加擬合精度下降速度增大，當擬合時長為4小時，擬合誤差最大為7 mm，當擬合時長為5小時，擬合誤差最大達3 cm。

圖1 MEO衛(wèi)星擬合URE隨擬合點數(shù)變化圖Fig.1 Fitting URE of MEO satellite changed with the fitting point numbers

圖2 GEO衛(wèi)星擬合URE隨擬合點數(shù)變化圖Fig.2 Fitting URE of GEO satellite changed with the fitting point numbers

圖3 IGSO衛(wèi)星擬合URE隨擬合點數(shù)變化圖Fig.3 Fitting URE of IGSO satellite changed with the fitting point numbers

表1 MEO衛(wèi)星不同擬合時長擬合URE表Tab.1 Fitting URE of MEO satellite with different length of time

表2 GEO衛(wèi)星不同擬合時長擬合URETab.2 Fitting URE of GEO satellite with different length of time

表3 IGSO衛(wèi)星不同擬合時長擬合URETab.3 Fitting URE of IGSO satellite with different length of time

4 結論

1)算例驗證了本文推導的廣播星歷參數(shù)擬合算法以及偏導數(shù)公式是正確的，且這種新型的廣播星歷18參數(shù)模型同樣適用于北斗衛(wèi)星導航定位系統(tǒng);

2)當擬合時間不超過4小時，各類衛(wèi)星擬合精度優(yōu)于0.01米，完全可以滿足實時導航定位的精度要求。但需要注意的是當擬合時段如果小于1小時，在迭代計算時部分參數(shù)擬合誤差增大，擬合迭代次數(shù)增加，甚至出現(xiàn)迭代不收斂情況，因此擬合時段需大于1小時結果才有效;

3)此種擬合算法的內(nèi)部精度穩(wěn)定，且隨著擬合時間的增長擬合精度逐步下降。

1 崔先強，焦文海，賈小林.GPS廣播星歷參數(shù)擬合算法［J］.測繪學院學報，2004，21(4):244 －246.(Cui Xianqiang，Jiao Wenhai and Jia Xiaolin.GPS broadcast ephemeris parameter fitting algorithm［J］.Journal of Institute of Surveying and Mapping，2004，21(4):244 －246)

2 戴曉蕾，等.GPS兩種廣播星歷參數(shù)擬合及外推精度分析［J］.大地測量與地球動力學，2013.(1):53－58.(Dai Xiaolei，et cl.Accuracy analysis of two GPS broadcast epheemeris fitting and extrapotation algorithm［J］.Journal of Geodesy and Geodynamics，2013，(1):53 －58)

3 陳留成，唐波.參考系選擇對Kepler廣播星歷參數(shù)擬合精度的影響［J］.飛行器測控學報，2006，25(4):19－25.(Chen Liucheng and Tang Bo.Impact of reference system choice on Kepler fitting of broadcast ephemeris earameters precision［J］.Journal of Aircraft Control，2006，25(4):19－25)

8 劉基余.GPS衛(wèi)星導航定位原理與方法［M］.北京:科學出版社2008.(Liu Jiyu.GPS satellite navigation principle and method［M］.Beijing:Science Press，2008)

5 崔先強，等.兩種GPS廣播星歷參數(shù)算法的比較［J］.空間科學學報，2006，26(5):382 －387.(Cui Xianqiang，et al.Comparison of two kinds of algorithm of GPS broadcast ephemeris parameters［J］.Journal of Space Science，2006，26(5):382－387)

6 呂志偉，易維勇，曾志林.GPS廣播星歷參數(shù)擬合算法及其分析［J］.測繪科學技術學報，2010，27(2):83－85.(Lü Zhiwei，Yi Weiyong and Zeng Zhilin.GPS broadcast ephemeris parameter fitting algorithm and analysis［J］.Geomatics Science and Technology，2010，27(2):83 －85)

BROADCAST EPHEMERIS PARAMETER FITTING AND PRECISION ANALYSIS OF BEIDOU SATELLITE

Wang Le1)and Jia Xiaolin2)
1)School of Geology Engineering and Geomatics，Chang’an University，Xi’an710054
2)Xi’an Research Institute of Surverying and Mapping，Xi’an710054

The GPS ICD document defined a new group of 18 parameters broadcast ephemeris model，according to the user algorithm of this new parameter model，the partial derivative expression of each parameter is derived，and the linear function model is established between the 18 parameters and satellite position vector.In order to test this algorithm，the Beidou satellite orbit determination data are used to calculate the 18 parameters of Beidou satellite broadcast ephemeris by the principle of least－squares adjustment.The result verified the correctness of the algorithm，and simultaneously proved 4 hours fitting error is better than 0.01 m.

Beidou satellite;user algorithmbroadcast ephemeris parameter;fitting analysis;GPS ICD document;18 parameters broadcast ephemeris model

P207

A

1671－5942(2013)05－0071－05

2013－04－08

王樂，女，1988年生，碩士研究生，研究方向為衛(wèi)星定軌.E－mail:815557438@qq.com