BDS衛星精密星歷中不同類型軌道的插值分析

李曉光，程鵬飛，成英燕

(1.遼寧工程技術大學，遼寧 阜新　123000；2.中國測繪科學研究院，北京　100830)

?

BDS衛星精密星歷中不同類型軌道的插值分析

李曉光1，2，程鵬飛2，成英燕2

(1.遼寧工程技術大學，遼寧 阜新123000；2.中國測繪科學研究院，北京100830)

為了進一步研究使用精密星歷獲取BDS衛星精確的、更高采樣率的軌道位置的問題，提出應用拉格朗日和牛頓插值對3類BDS衛星軌道進行內插的方法，分析了2種內插對3類衛星的適用性。結果表明，要想得到最佳收斂效果應用的收斂階數是不同的，GEO最佳收斂在6～7階，IGSO最佳收斂在8～10階，MEO最佳收斂在10～11階。

BDS；精密星歷；拉格朗日插值；牛頓插值；GEO

0　引言

在北斗衛星導航系統(BeiDou navigation satellite system，BDS)精密星歷的應用中，正確獲取BDS衛星精確的軌道位置是需要解決的基礎問題。衛星精密星歷采用sp3格式，其存儲方式為ASCII文本文件，內容包括表頭信息以及文件體[1-2]，文件體中每隔15 min給出一個衛星的位置[3]，有時還給出衛星的速度和鐘差，主要提供衛星精確的軌道位置。而BDS接收機的采樣率一般為30 s或者15 s，甚至更密；因此要想利用某一時刻的衛星位置，就必須對精密星歷進行高精度的插值。

本文利用武漢大學發布的2015年第74天廣播星歷采用拉格朗日和牛頓插值方法進行插值，并進行比較分析。

1　Lagrange插值方法原理

L1(xk)=yk,L1(xk+1)=yk+1。

(1)

則其線形插值多項式為

L1(x)=yklk(x)+yk+1lk+1(x)。

(2)

其中

式(2)稱為Lagrange線形插值基函數。

同理，相應的Lagrange插值多項式為

(3)

其中

式(3)是n次插值基函數，滿足條件

(4)

對于n階插值，即有n+1個已知點，是內插位于這n+1個點之間的任意位置的函數值。

拉格朗日多項式插值模型簡單，是經典的插值算法；但當精度不夠而需要增加新的插值節點時，原來的插值多項式包括連續乘積項都不能使用，必須重新構造1個插值多項式。

2　Newton插值方法原理

Newton插值將插值基函數[4]定義為

(5)

其多項式的系數是各階均差，定義為：

1)1階均差

(6)

2)2階均差

(7)

3)k階均差

(8)

根據均差定義，高階均差是低1階均差的均差，因此可以用表 1 中遞推關系來計算各階均差。

4)Newton均差插值多項式

(9)

從上述定義可以看出：基函數的系數就是表1中主對角線元素，于是在實際計算中，可以先將表1列出。當增加節點時，只需在表的最下面增加1行，右面增加1列高階均差即可，原來計算全部有效。這樣就避免了增減節點引起插值基函數的改變，導致整個公式的變化。

表1　Newton各階均差表

3　BDS精密星歷內插

本文采用2015年第74天 05∶00∶00至10∶00∶00的3類衛星的精密星歷，運用MATLAB[5]對07∶30∶00時刻衛星軌道進行2種方法的插值，并對原精密星歷數據進行對比。

1)GEO衛星。BDS在軌工作衛星有5顆地球靜止軌道(geostationary Earth orbit，GEO)衛星[6](如表2所示)。選取C01號衛星為案例進行插值分析。由于GEO衛星是靜止軌道衛星，因此這類衛星坐標基本不變[7]；但在3個方向也會呈現周期變化如圖1～圖3，取7 d的精密星歷坐標畫出。

表2　C01衛星插值結果　mm

圖1　X方向坐標變化

圖2　Y方向坐標變化

圖3　Z方向坐標變化

由圖1～圖3可以看出GEO衛星在X、Y方向上的變化存在1 d和0.5 d周期項，在Z方向上存在1 d周期項，對于它的周期變化可能是由衛星機動[8]、太陽光壓等因素造成的。

對于低階插值，2種方法都能夠迅速收斂。對于X、Z方向上的殘差，在6階插值時達到最低的mm級，階數大于6以后，用增加階數的方法提高精度不顯著，存在震蕩現象；Y方向上7階達到最低的mm級。在效率方面，由于只解算1組數據，耗時都很少，但Newton插值方法要比Lagrange插值方法快2/10。對于大數據解算來說，會產生較大時間差。

對于其他GEO衛星進行插值都能迅速收斂，最佳收斂階數大多分布在6～7階。

2)IGSO衛星。BDS在軌工作衛星有5 顆傾斜地球同步軌道(inclined geo-synchronous orbits，IGSO)衛星[6]。選取C06號衛星為案例進行插值分析，如表3所示。

對于IGSO衛星，2種插值也能迅速收斂，但較GEO衛星來說相對慢些：Y、Z坐標方向上在8階插值時達到最佳插值效果；而X方向收斂較慢，在10階插值時達到最佳插值效果。在計算效率上，Newton插值方法比Lagrange插值方法節約了大約1/3的時間。

對于其他IGSO衛星，通過計算得知最佳收斂階數大多分布在8～10階。

3)MEO衛星。BDS在軌工作衛星有4 顆中圓地球軌道(medium Earth orbit，MEO)衛星[6]。選取C11號衛星為案例進行插值分析，如表4所示。

表3　C06衛星插值結果　mm

表4　C11衛星插值結果　mm

MEO衛星是3類衛星中收斂速度最慢的衛星。3個方向最佳收斂階數產生較大差異，X方向在11階達到最佳，Y方向在10階達到最佳，而Z方向收斂較快，在8階達到最佳。階數的增加對效果沒有產生更好的提高，多在mm級震蕩。MEO衛星和GPS衛星軌道相近，但與眾多關于GPS衛星軌道插值的文獻比較發現[9-11]，GPS衛星軌道插值最優收斂階在8～9階。在計算效率上，Newton插值方法更節約時間。

其他MEO衛星收斂較慢，最佳收斂階數在10～11階，高階后產生震蕩。

4　結束語

通過利用2種插值方法對3類BDS衛星精密星歷進行內插，數據結果分析得出如下結論：

1)對BDS精密星歷進行內插，充分利用待插節點前后的數據信息能夠充分體現待插節點附近衛星運動的基本規律，插值精度較高。

2)BDS不同類型軌道的衛星要達到最優收斂效果，應用的內插階數不同。GEO衛星收斂最快，最佳收斂階數大多分布在6～7階；IGSO衛星收斂較慢，最佳收斂階數大多分布在8～10階；MEO衛星收斂最慢，最佳收斂階數在10～11階。Lagrange方法隨著插值階數的增大，舍入誤差會造成龍格現象的出現。因此插值階數越高，插值效果不一定更好。

3)對于2種插值方法，從插值精度上來說2者插值效果幾乎一樣；在編程實現上，Lagrange插值方法原理簡單易懂，更易于編程；從運算效率上來說Newton插值方法更節約時間，對于大數據計算而言，建議采用Newton插值方法，這樣也避免了龍格現象。

[1]楊元喜.BDS衛星導航系統的進展、貢獻與挑戰[J].測繪學報，2010，39(1)：1-6.

[2]李征航，黃勁松.GPS測量與數據處理[M].武漢：武漢大學出版社，2005：17-28.

[3]魏二虎，柴華.GPS精密星歷插值方法的比較研究[J].全球定位系統，2006，31(5)：13-15.

[4]李慶楊，關治，白峰衫.數值計算原理[M].北京：清華大學出版社，2000：5-20.

[5]龔純，王正林.MATLAB語言常用算法程序集[M].電子工業出版社，2010：30-40.

[6]施闖，趙齊樂，李敏，等.北斗衛星導航系統的精密定軌與定位研究[J].中國科學：地球科學，2012，42(6)：854-861.

[7]ASHKENAZI V，CHEN W，HILL C J，et al.Real-time autonomous orbit determination of LEO satellites using GPS[C]//The Institute of Navigation(ION).Proceedings of the 10th International Technical Meeting of the Satellite Division of the Institute of Navigation，Proceedings of -97.Kansas City，Missouri：ION，1997：16-19.

[8]徐天河，賀凱飛.顧及系統誤差的GEO衛星幾何法定軌[J].大地測量與地球動力學，2009，29(1)：64-69.

[9]蔡艷輝，程鵬飛，李夕銀.衛星坐標的內插和擬合[J].全球定位系統，2003，28(3)：10-13.

[10]何玉晶，楊九.基于拉格朗日插值方法的GPS IGS精密星歷值分析[J].測繪工程，2011，20(5)：60-62.

[11]洪櫻，歐吉坤，彭碧波.GPS衛星精密星歷和鐘差三種內插方法比較[J].武漢大學學報·信息科學版，2006，31(6)：516-518.

Analysis on BDS satellite precise ephemeris interpolation in different orbits

LI Xiaoguang1，2，CHENG Pengfei2，CHENG Yingyan2

(1.School of Geomatics，Liaoning Technical University，Fuxin，Liaoning 123000，China；2.Chinese Academy of Surveying and Mapping，Beijing 100830，China)

In order to further study on using BDS satellite ephemeris to obtain accurate orbital position with higher sampling rate，the paper used Lagrange and Newton interpolation to interpolate three types of BDS satellite orbits in this paper，and analyzed the feasibility of the two algorithms on the three types of satellites.Result showed that for getting the best effect of convergence，the convergence orders used with the satellites should be different，as the GEO optimal convergence in order 6～7，IGSO in 8～10，and MEO in 10～11 orders.

BDS；precise ephemeris；Lagrange interpolation；Newton interpolation；GEO

2016-04-21

國家自然科學基金項目(41374014)；中國測繪科學研究院基本科研業務費支持項目(7771405)。

李曉光(1991—)，男，山東濱州人，研究研究生，研究方向為北斗雙差等科研和應用。

10.16547/j.cnki.10-1096.20160307.

P228

A

2095-4999(2016)03-0031-04

引文格式：李曉光，程鵬飛，成英燕.BDS衛星精密星歷中不同類型軌道的插值分析[J].導航定位學報，2016，4(3)：31-34.(LI Xiaoguang，CHENG Pengfei，CHENG Yingyan.Analysis on BDS satellite precise ephemeris interpolation in different orbits[J].Journal of Navigation and Positioning，2016，4(3)：31-34.)