廣義延拓插值法在衛星精密星歷中的應用

申俊飛，鄭 沖，陳金春

（北京衛星導航中心，北京 100094）

1 引言

在衛星導航定位中需要不斷地計算出導航衛星的瞬時坐標和速度，而衛星坐標和速度是根據接收到的廣播星歷來計算的，而每次計算都需要占用大量內存，影響計算速度。衛星在軌道上的位置和速度用星歷表示，具體形式可以是衛星位置和速度的時間列表，也可以是一組以時間為引數的軌道參數。衛星通過發布星歷參數為用戶提供實時定位、精確導航以及授時等服務；用戶根據提供的星歷參數由用戶終端設備按固定公式計算出的衛星位置和速度稱之為導航衛星廣播星歷。當然廣播星歷還包括其他一些如衛星鐘等參數，這里只討論星歷中發播的軌道參數［1］。

廣播星歷是通過衛星導航系統的測軌分系統外推軌道得到，隨著時間的推移，預報軌道會偏離實際軌道，導航和定位的精度降低，所以只能保證一定時間段內的精度要求。隨著衛星導航系統的不斷發展，對廣播星歷的精度和實時性提出了更高的要求，因此，如何快速、精確、方便地為用戶提供導航衛星的位置和速度信息是衛星導航定位應用中的一個極其重要的問題［2］。

2 新插值方法

考慮到提供給用戶的廣播星歷的精確性和實時性要求，廣播星歷應是盡可能精確的提供衛星軌道預報值，即考慮了衛星所受的地心引力和各種攝動力；廣播星歷參數能很好的擬合這些值，計算誤差應盡量小，用戶利用廣播星歷參數計算衛星位置的過程應盡量簡單。

為此，將導航衛星星歷表示為某種特定的多項式（即將一定時間段內的衛星星歷擬合于選定的多項式），然后將多項式系數作為星歷參數廣播給用戶是一種可行的方式。目前在精密星歷數據處理中，應用較多的是拉格朗日多項式插值，切比雪夫多項式在軌道擬合標準化中也由于其良好的插值精度有著廣泛的應用［3－8］。

文獻 ［9］提出了廣義延拓插值法的基本原理與方法，并對其精度進行評價，但沒有探討3個參數的取值對插值結果的影響。本文采用廣義延拓插值法進行導航衛星廣播星歷的插值，探討了3個參數的取值與插值精度之間的關系。

3 廣義延拓插值原理

廣義延拓插值是在插值法和擬合法的基礎上發展起來的一種插值方法。它一方面在分片邊界點上滿足插值條件，使得分片之間的變化具有一定的協調性；另一方面又利用分片插值區域的周圍節點的信息，實現分片區域內部的最佳擬合。

設函數u有一組離散數據

欲在區域Ω上構造一個近似函數U滿足

將Ω劃分為m個互不重疊的子區域（e＝1，2，…，m），稱作單元域Ωe，單元域共有r個節點，其相應的節點坐標為xi（i＝1，2，…，r），節點函數值為ui（i＝1，2，…，r）。將單元域拓展，得單元延拓域Ω′e。

設Ω′e內有s個節點，其中屬于Ωe的有r個節點（r＜s）。利用延拓域內信息來構造單元域內逼近函數。

式（3）中，g1（x），g2（x），…，gt（x）為Ω′e上的一組基函數；ak（k＝1，2，…，t）為待定系數；t為逼近函數項數，且r＜t＜s。

利用最佳平方逼近，求a1，a2，…，at，使

同時，使Ue（x）滿足Ωe上的插值條件：

將兩式合并寫成為

由式（6）解出逼近函數式（3）中的待定系數，便可確定Ωe上的插值函數。式（6）稱為廣義延拓逼近方法的數學模型。

4 結果分析及應用

本文中用美國國家大地測量機構（national geodetic survey，NGS）提供的精密星歷來進行廣義延拓插值計算，然后對結果進行分析。軌道擬合的效果受到多項式階數、擬合時段長度等參數的影響，分析這些參數的影響以及參數間的關系對于更好地使用廣播星歷有一定的意義。

目前NGS精密星歷采樣間隔為15min 本文采用 NGS跟蹤站提供的2002－01－01PRN1衛星15min等間隔的GPS精密星歷，利用廣義延拓插值法將15min間隔的數據插值到5min間隔，依次改變單元域節點數r、延拓域節點數s和逼近函數的項數t等3個參數的取值，組成多種組合，比較在取不同參數時插值的精度，分析參數和精度之間的關系，并與GPS精密星歷中的衛星坐標值進行比較，計算結果見表1。

表1 插值結果與GPS精密坐標比較

表1中給出了3個參數取得不同值時，廣義延拓插值結果與精密星歷坐標的最大差值。可以看出，3個參數的取值和插值精度之間存在著一定的規律，當3個參數取值較小時，插值誤差較大，r＝3，s＝5，t＝4時最大誤差可達到90m以上。隨著r，s，t等3個參數的增大，插值的精度不斷提高，在一定的范圍內，精度保持穩定，具有較高的精度。從表中數據可以看出，插值精度與參數s的關系最為密切，可以看到當s大于10時才能得到較高的插值精度；而當s的取值固定時，插值精度隨著t的變化而改變，t和s越接近，插值精度就越高，當時，精度最高。一般地，當s大于10后，再增加3個參數的值，插值精度基本保持穩定，不會有明顯的變化，可以認為r＝7，s＝11，t＝10即滿足需要。X、Y、Z的最大誤差為2.7mm2.1mm1.8mm 具有很高的精 度，達到毫米級。參數繼續增加時，計算量增加，插值精度反而不斷降低。

由以上可見用廣義延拓插值法是可以滿足導航定位精度的要求。對于精密廣播星歷來說，如果采用r、s、t三個參數分別為7、11、10時，在廣播星歷播送時，只需要知道X、Y、Z三個方向上的系數，即t＝10時共計30個參數即可。

并且這30個參數可由主控站將衛星在在地球協議坐標系中的軌道直接擬合得到。用戶收到這30個精密廣播星歷參數后，應用廣義延拓逼近多項式，可直接計算衛星在地球協議坐標系中對應時刻的位置和速度，無需進行坐標轉換，物理意義明確，計算簡單，速度快。

5 結論

通過以上數據計算統計的結果和分析，可以初步得出以下結論：

（1）用項數為10的廣義延拓逼近多項式逼近擬合衛星軌道，將多項式系數作為導航衛星精密星歷參數廣播，這是一種新的廣播星歷形式，導航定位終端利用這組星歷參數解算實時導航定位時的所需衛星軌道，方法簡單、快捷。

（2）利用這組星歷參數解算衛星軌道，計算過程簡潔，計算速度快，占用計算機內存空間少、能滿足精度要求等優點，但同時可能造成廣播星歷中信息的損失，還需要進一步研究；同時對廣義延拓逼近多項式系數可根據實際精度需要進行合理選擇。該方法在GPS定位測量中為實現衛星軌道標準化提供了一種較好的、行之有效的方法。

（3）這種方法是否適用于具有多種軌道類型的導航衛星系統還待研究，這對于我國具有多種軌道類型共存的北斗衛星導航系統具有一定的參考價值。

［1］ 朱俊，文援蘭，廖瑛.一類適用于各種軌道類型的導航衛星廣播星歷研究［J］.航天控制，2005（6）：9－11.

［2］ 陳春花，高成發.衛星軌道標準化及其算法實現［J］.現代測繪，2008，31（3）：12－15.

［3］ 嚴麗，李萌.切比雪夫多項式擬合衛星軌道與鐘差的精度分析［J］.測繪科學，2013，38（3）：59－62.

［4］ 王義峰，文鴻雁，劉立龍.用切比雪夫多項式標準化 GPS衛星軌道［J］.測繪與空間地理信息，2008，31（1）：21－22.

［5］ 李明峰，江國焰，張凱.IGS精密星歷內插與擬合法精度的比較［J］.大地測量與地球動力學，2008，28（2）：77－80.

［6］ 魏二虎，柴華.GPS精密星歷插值方法的比較研究［J］.全球定位系統，2006，31（5）：13－15.

［7］ 孫華麗，張政治.全球定位系統衛星位置計算插值方法精度分析［J］.導航定位學報，2013，1（3）：89－93.

［8］ 劉偉平，郝金明.一種新的IGS精密星歷插值算法［J］.武漢大學學報：信息科學版，2011，36（11）：1320－1324.

［9］ 施滸立，顏毅華，徐國華.工程科學中的廣義延拓逼近法［M］.北京：科學出版社，2005.