利用滑動式Lagrange插值方法擬合衛(wèi)星精密星歷

郭忠臣

宿州學(xué)院環(huán)境與測繪工程學(xué)院，宿州，234000

利用滑動式Lagrange插值方法擬合衛(wèi)星精密星歷

郭忠臣

宿州學(xué)院環(huán)境與測繪工程學(xué)院，宿州，234000

為了得到精確的衛(wèi)星三維坐標(biāo)，應(yīng)用滑動式Lagrange插值方法對GPS精密星歷內(nèi)插，給出衛(wèi)星位置插值公式。通過設(shè)置不同的插值階數(shù)，對插值精度統(tǒng)計分析。結(jié)果表明：插值精度隨著階數(shù)的增加而提高，當(dāng)階數(shù)達到11階時，插值精度較高，X、Y和Z三個方向的RMS分別達到0.378、0.514、0.306 mm，且均值偏差都在0.1 mm左右，精度略優(yōu)于其他階數(shù)，可滿足導(dǎo)航方面的需求。

滑動式Lagrange；插值；精密星歷；GPS

在高精度衛(wèi)星導(dǎo)航和深空探測等領(lǐng)域中，需使用精確的衛(wèi)星三維坐標(biāo)。當(dāng)前使用較多的為IGS(International GNSS Service)發(fā)布的精密星歷產(chǎn)品，但精密星歷是按照一定時間間隔提供的[1-2]，若想得到任意時刻的衛(wèi)星位置，必須對精密星歷進行擬合和插值。本文采用滑動式Lagrange插值算法對IGS提供的精密星歷進行處理，并對插值結(jié)果進行統(tǒng)計分析，為后續(xù)研究提供參考。

1 滑動式Lagrange插值算法原理

利用Lagrange多項式可內(nèi)插出節(jié)點之內(nèi)任意時刻的衛(wèi)星位置。Lagrange插值函數(shù)可表示為[3-5]：

(1)

滑動式Lagrange插值算法是在原有方法的基礎(chǔ)上構(gòu)建一個活動的“窗口”，以保證插值點的位置時刻保持在插值范圍的中心。當(dāng)選取的節(jié)點數(shù)為奇數(shù)時，最佳插值點有2個，選取的節(jié)點數(shù)為偶數(shù)時，最佳插值點即為中間時刻。如選取9個節(jié)點時，可構(gòu)造8階Lagrange擬合多項式，此時最佳插值點位于第4至第5和第5至第6個節(jié)點之間；選取8個節(jié)點時，可構(gòu)造7階Lagrange擬合多項式，此時最佳插值點位于第4和第5個節(jié)點之間，具體說明見圖1。

圖1 滑動式Lagrange插值示意圖

根據(jù)IGS發(fā)布的精密星歷可得到各給定時刻衛(wèi)星的三維坐標(biāo)X(t)，本文選取偶數(shù)個節(jié)點作為插值點[3]，依據(jù)上述原理，可得任意時刻t的衛(wèi)星位置的表達式為[6]:

(2)

其中，X(t)、X(ti)為對應(yīng)時刻t和ti的位置，i,j≠0，2n-1為插值階數(shù)。

2 實驗分析

本文采用WHU站提供的時間間隔為15 min、包含衛(wèi)星坐標(biāo)和鐘差的1 791周(GPST)的精密星歷數(shù)據(jù)，共672組三維坐標(biāo)，提取該周內(nèi)以1 800 s為時間間隔(0 s，1 800 s，3 600 s，5 400 s，7 200 s，…，604 800 s)的坐標(biāo)作為已知數(shù)據(jù)，對本周內(nèi)26 100～578 700 s之間每隔1 800 s的時刻位置進行插值擬合。選取PG04號衛(wèi)星進行Lagrange插值計算，并對插值結(jié)果進行統(tǒng)計分析。

表1和圖2列出了使用不同階數(shù)的Lagrange插值精度的統(tǒng)計結(jié)果，通過分析可知：(1)隨著階數(shù)的增加，插值精度明顯提高，但當(dāng)階數(shù)增加到一定程度時，精度趨于穩(wěn)定；(2)當(dāng)階數(shù)增加到11階后，插值精度較高，X、Y和Z三個方向的RMS分別達到0.378、0.514、0.306 mm，且均值偏差都在0.1 mm左右，精度略優(yōu)于其他階數(shù)；(3)隨著階數(shù)增加，插值過程中的計算量增加，對計算機的耗損增加，當(dāng)需要大量插值結(jié)果時，可從精度與耗時等方面綜合考慮，以確定最佳階數(shù)。

圖2 不同階數(shù)插值精度統(tǒng)計

表1 不同階數(shù)插值精度統(tǒng)計

表2列出了部分時刻的精密星歷和當(dāng)階數(shù)為11階時的插值結(jié)果，表3給出了相應(yīng)時刻的插值誤差。由表3可知，當(dāng)階數(shù)為11階時，X、Y和Z三個方向的插值誤差均優(yōu)于1 mm，能滿足導(dǎo)航方面的需求。

表2 部分時刻的精密星歷與插值結(jié)果

4 結(jié)束語

(1)滑動式Lagrange插值算法因算法簡單、插值精度高而廣泛應(yīng)用于各個領(lǐng)域，用它對精密星歷插值，可使插值誤差保持在1 mm以內(nèi)，基本滿足導(dǎo)航方面的需求。

(2)本文選取不同階數(shù)進行插值，插值精度隨著階數(shù)的增加而增加，但當(dāng)階數(shù)增加到一定程度時，插值精度趨于穩(wěn)定，甚至?xí)陀诘碗A的插值精度。

(3)文中選取了一周的精密星歷進行插值計算，但若想對該周的起始與結(jié)束時刻進行插值，需聯(lián)合使用前一周和后一周數(shù)據(jù)。

表3 部分時刻插值誤差

[1]王興,高井祥,王堅,等.利用滑動式切比雪夫多項式擬合衛(wèi)星精密坐標(biāo)和鐘差[J].測繪通報,2015(5):6-8

[2]何玉晶,楊力.基于拉格朗日插值方法的GPSIGS精密星歷插值分析[J].測繪工程,2011,20(5):60-62

[3]雷雨,趙丹寧,高玉平,等.基于滑動式Lagrange插值方法的GPS精密星歷內(nèi)插分析[J].測繪工程,2013,22(2):34-36

[4]向榮榮,蘭孝奇,劉鋒,等.拉格朗日和切比雪夫多項式在精密星歷插值中的應(yīng)用[J].測繪與空間地理信息,2015,38(11):131-134

[5]何玉晶,楊力.基于拉格朗日插值方法的GPSIGS精密星歷插值分析[J].測繪工程,2011,20(5):60-62

[6]郭忠臣,高井祥,曹新運.北斗廣播星歷軌道及鐘差精度分析[J].煤炭技術(shù),2015,34(7):113-116

(責(zé)任編輯：汪材印)

P228.4

：A

：1673-2006(2017)07-0106-03

10.3969/j.issn.1673-2006.2017.07.027

2017-03-09

郭忠臣(1992-)，安徽阜陽人，碩士，助教，研究方向：GNSS數(shù)據(jù)處理。