馬卓希,楊力,賈小林,趙海山

(1.信息工程大學(xué) 導(dǎo)航與空天目標工程學(xué)院,鄭州 450001;2 西安測繪研究所,西安 710054)

基于多項式模型的改進衛(wèi)星鐘差預(yù)報方法

馬卓希1,楊力1,賈小林2,趙海山1

(1.信息工程大學(xué) 導(dǎo)航與空天目標工程學(xué)院,鄭州 450001;2 西安測繪研究所,西安 710054)

摘要：根據(jù)衛(wèi)星鐘固有的特性,提出了一種通用的衛(wèi)星鐘差預(yù)報方法。該方法在多項式模型的基礎(chǔ)上,對多項式模型的選擇、周期項修正、殘差建模等方面進行了改進。首先利用擬合數(shù)據(jù)的殘差標準差指標進行多項式模型的自適應(yīng)選擇,繼而利用頻譜分析的方法對鐘差序列的周期項進行了建模,最后利用AR模型對消除趨勢項的隨機項進行建模。試驗結(jié)果表明，改進的鐘差預(yù)報模型能夠顯著地提高鐘差的預(yù)報精度,對于穩(wěn)定度較差的衛(wèi)星的預(yù)報效果尤為顯著。改進模型的衛(wèi)星鐘預(yù)報6 h預(yù)報精度總體優(yōu)于1 ns,24 h預(yù)報精度優(yōu)于2 ns.

關(guān)鍵詞：GPS;星載原子鐘;多項式模型;鐘差預(yù)報

0引言

衛(wèi)星軌道及精密鐘差的預(yù)報一直是實時高精度定位領(lǐng)域的重要研究內(nèi)容。IGS組織向用戶發(fā)布事后、快速及超快速的衛(wèi)星軌道和鐘差產(chǎn)品。GPS的超快速的預(yù)報軌道和鐘差產(chǎn)品的精度分別為5 cm和3 ns左右[1]。近年來,精密軌道的預(yù)報精度得到了較大幅度的改進,然而衛(wèi)星鐘差的預(yù)報精度相對較差,難以滿足亞分米級的實時精密單點定位的需求[2]。

為了提高衛(wèi)星預(yù)報精度,國內(nèi)外學(xué)者提出了多種衛(wèi)星鐘預(yù)報的模型和方法,常見的衛(wèi)星鐘預(yù)報模型有:多項式預(yù)報模型、灰色模型、求和自回歸滑動平均模型等[3-7]。鄭作亞等對不同預(yù)報模型的預(yù)報精度進行了對比分析[3];黃觀文等利用頻譜分析的方法,對鐘差序列的周期項進行建模和預(yù)報,并取得了較高的預(yù)報精度[4];范旭亮等基于ARIMA模型對衛(wèi)星鐘差的短期預(yù)報性能進行了分析[5]。盡管國內(nèi)外學(xué)者進行了大量的研究工作,但仍然存在預(yù)報算法難以克服粗差、相位跳變等數(shù)據(jù)異常、預(yù)報模型多樣化但普適性不強、預(yù)報模型沒有充分結(jié)合原子鐘的物理特性等諸多問題,上述問題將極大程度地制約衛(wèi)星導(dǎo)航系統(tǒng)的服務(wù)性能。因此,研究高精度的衛(wèi)星鐘差預(yù)報算法,進一步提高衛(wèi)星鐘

差的預(yù)報精度和適應(yīng)性具有重要的意義。

本文在多項式預(yù)報模型的基礎(chǔ)上,提出了一種改進的衛(wèi)星鐘差預(yù)報方法。先詳細闡述了其相關(guān)理論及實現(xiàn)方法,然后結(jié)合算例,利用IGS實測的鐘差數(shù)據(jù)對算法的短期預(yù)報效果進行了論證。

1改進的衛(wèi)星鐘差預(yù)報方法

由于多項式模型基本能反映衛(wèi)星鐘鐘速和鐘漂的物理特性,且算法的復(fù)雜程度較低,因此,本文采用的預(yù)報模型基于多項式模型進行改進?？傮w而言,算法的實現(xiàn)方法共包含三個步驟: 1) 進行數(shù)據(jù)預(yù)處理。排除粗差、相位跳變、數(shù)據(jù)中斷等異常對后續(xù)處理的影響; 2) 在多項式模型的基礎(chǔ)上,對多項式模型的選擇進行改進,并對周期項進行建模; 3 )利用AR模型對消除趨勢項后的殘差進行建模[8]。算法的實現(xiàn)流程圖如圖1所示。本文重點對模型的改進策略及建模方法進行闡述。

圖1　改進的衛(wèi)星鐘差預(yù)報算法流程圖

1.1多項式模型的自適應(yīng)選擇

對于衛(wèi)星鐘差預(yù)報而言,選擇多項式模型的擬合階次主要取決于衛(wèi)星鐘的頻漂特性是否明顯。由于星載原子鐘的變化具有不穩(wěn)定性。不同衛(wèi)星會表現(xiàn)出不同的物理特性;對于同一顆衛(wèi)星而言,在不同的時間段的物理特性也會發(fā)生變化。因此,在實際應(yīng)用中,不能籠統(tǒng)地將其歸類至單一的模型,需要利用擬合數(shù)據(jù)段的殘差標準差指標進行多項式模型的自適應(yīng)選擇。模型的具體表達式及判斷準則為

(1)

式中: x(tj)表示第j個歷元的衛(wèi)星鐘差值; tj表示第j個歷元的時刻; v(j)為該衛(wèi)星鐘第j個歷元鐘差的殘差; RMSLine、RMSQuadratic分別為利用一次多項式及二次多項式擬合的殘差標準差。

1.2周期項建模

利用二次多項式模型對GPS所有在軌衛(wèi)星的1天數(shù)據(jù)進行擬合,其擬合的殘差序列如圖2所示。從圖2中可以發(fā)現(xiàn),幾乎所有的在軌衛(wèi)星的擬合殘差序列均表現(xiàn)出不同幅度的周期特性,其周期特性將對鐘差序列的擬合及后續(xù)的預(yù)報工作產(chǎn)生較大的負面影響。受星載原子鐘的穩(wěn)定度因素影響,其中G04、G10、G29的擬合殘差較大。

圖2　GPS在軌衛(wèi)星擬合殘差序列圖

采用頻譜分析的方法對鐘差序列的周期項進行建模。頻譜分析主要利用數(shù)據(jù)形式呈波形的性質(zhì),通過對數(shù)據(jù)波形的幅值進行分析,對頻譜圖中具有較大幅值特性的顯著周期項進行提取,達到消除系統(tǒng)誤差的目的。對于鐘差殘差序列中的系統(tǒng)性誤差,通過快速傅立葉級數(shù)展開式計算序列中各數(shù)值的頻譜值。傅立葉變換的展開式

(2)

式中,X(k)為k時段的頻譜值,x(n)為時間序列的值。其具體的理論及實現(xiàn)見參考文獻[9-10]。在公式1的基礎(chǔ)上,顧及周期項改正的多項式模型為

(3)

式中: x0,y0,β,Am和Bm分別為待求的模型參數(shù),m=1,2,…n,m的大小由功率譜決定; ωk表示該衛(wèi)星鐘的頻率變化特征參數(shù);利用頻譜分析進行周期項建模的具體步驟為

1) 采用多項式模型對觀測數(shù)據(jù)進行擬合,得到其擬合殘差序列vi.

2) 對擬合殘差vi進行降噪處理,利用快速傅立葉變換進行處理,通過對其功率譜進行分析,確定其顯著周期項,在程序設(shè)計中只選取了前6個最大的功率譜對應(yīng)的頻率。

3) 利用最小二乘的參數(shù)方法求解待估參數(shù),并利用求解的參數(shù)進行后續(xù)的鐘差預(yù)報。

1.3隨機項建模

通過周期項建模使得鐘差序列的趨勢項得到較大幅度的削弱,利用AR模型對扣除多項式項及周期項后的殘差v(j)進行建模。AR模型,即自回歸模型,由ARIMA模型簡化而來,使用該模型的一個重要前提是需要分析的時間序列是平穩(wěn)的。在實際應(yīng)用中,通常通過差分的方法消除數(shù)據(jù)的趨勢項,從而滿足模型對數(shù)據(jù)平穩(wěn)性的要求。基于差分處理的ARIMA模型通常記為{Xt}～ARIMA(p,d,q)。其中,{Xt}為時間序列,p,q為模型階數(shù),d為差分的次數(shù)。當(dāng)d為0時,其數(shù)學(xué)表達式為

Xt=φ1Xt-1+φ2Xt-2+…+φpXt-p-

θ1at-1-θ2at-2-…-θqat-q+at,

(4)

式中: φ,θ為待估參數(shù); {at}為白噪聲序列,當(dāng)p,q分別為0時,模型可相應(yīng)地變化為MA模型和AR模型。

為了確定更為準確的模型參數(shù),可通過AIC準則確定模型參數(shù),對于AR模型而言,其AIC判別函數(shù)

(5)

2試驗和結(jié)果

本文所提出的鐘差模型的預(yù)報精度,利用IGS采樣率為30 s的精密鐘差數(shù)據(jù)進行試驗分析。采用的數(shù)據(jù)的時段為2015年3月1日至3月14日,共兩周的數(shù)據(jù),并以IGS精密鐘差值作為基準,計算預(yù)報鐘差與鐘差基準值的殘差和均方根誤差(RMS),并以此作為評估依據(jù),評估衛(wèi)星鐘差的預(yù)報精度。由于G26衛(wèi)星在采用的數(shù)據(jù)時段內(nèi)未提供服務(wù),因此,本文及后續(xù)實驗中未涉及對該顆衛(wèi)星的結(jié)果分析。

為了分析改進方法利用不同擬合時長的預(yù)報精度,試驗分別采用6 h、12 h、18 h和24 h的數(shù)據(jù)作為歷史數(shù)據(jù)進行數(shù)據(jù)擬合,表1～表3分別示出了不同類型的衛(wèi)星利用不同模型預(yù)報的均方根誤差。其中BLOCK IIA衛(wèi)星采用一次多項式模型的預(yù)報結(jié)果與之對比,其余衛(wèi)星采用二次多項式模型進行對比分析[11-12]。表中的預(yù)報精度為利用2周數(shù)據(jù)計算的平均統(tǒng)計結(jié)果。

表1　不同類型衛(wèi)星6 h預(yù)報精度(單位:ns)

表2　不同類型衛(wèi)星12 h預(yù)報精度(單位:ns)

表3　不同類型衛(wèi)星24 h預(yù)報精度(單位:ns)

從表1～表3中可以得到,總體而言,采用改進的預(yù)報算法的預(yù)報精度優(yōu)于多項式模型,能夠較為顯著地提高不同類型衛(wèi)星鐘差的預(yù)報精度。改進算法的衛(wèi)星鐘預(yù)報6 h預(yù)報精度總體優(yōu)于1 ns,24 h預(yù)報精度優(yōu)于2 ns.

從數(shù)據(jù)擬合長度與預(yù)報精度的關(guān)系方面分析,并非擬合弧段越長,預(yù)報的精度越高。其原因可歸結(jié)為擬合時段距離預(yù)測時段越遠,其對預(yù)報的相關(guān)性越弱。但當(dāng)采用的擬合數(shù)據(jù)量較少時,由于對少量數(shù)據(jù)擬合并不能完全反映星載原子鐘物理特性的變化規(guī)律,因此預(yù)報精度較低。通過表中數(shù)據(jù)可以得到看出,當(dāng)擬合弧段為6 h時,采用改進模型的預(yù)報精度明顯低于采用更長擬合弧段的預(yù)報精度;當(dāng)擬合弧段為12 h時,改進的預(yù)報模型能夠達到較優(yōu)的預(yù)報效果,當(dāng)擬合數(shù)據(jù)量繼續(xù)增多時,預(yù)報精度稍有降低,但變化并不顯著。

為了更加直觀地分析不同星載原子鐘的預(yù)報效果,利用2015年3月13日的精密鐘差數(shù)據(jù)進行擬合并預(yù)報后續(xù)24 h的衛(wèi)星鐘差。對不同類型星載原子鐘的預(yù)報鐘差與IGS鐘差基準的殘差進行統(tǒng)計,其預(yù)報殘差隨時間序列的變化關(guān)系如圖3～圖6所示。

圖4　BLOCK IIR-M在軌衛(wèi)星預(yù)報殘差序列 　　　(a) 預(yù)報模型； (b)改進的預(yù)報模型

圖5　BLOCK IIA在軌衛(wèi)星預(yù)報殘差序列 　　　(a) 預(yù)報模型； (b)改進的預(yù)報模型

圖6　BLOCK IIR在軌衛(wèi)星預(yù)報殘差序列 　　　(a) 預(yù)報模型； (b)改進的預(yù)報模型

從圖3～圖6中可以得到,衛(wèi)星鐘差的預(yù)報精度與星載原子鐘的穩(wěn)定度相關(guān),星載原子鐘的穩(wěn)定度越差,其預(yù)報的精度越低。同其他類型的衛(wèi)星相比,BLOCK IIA衛(wèi)星的預(yù)報精度最低,且預(yù)報序列呈現(xiàn)不規(guī)則的變化特性。由于改進的模型充分對周期項、隨機項進行了建模,同多項式模型相比,改進模型的誤差積累的程度有較大幅度的降低,對于穩(wěn)定度較差的衛(wèi)星預(yù)報的改進效果尤為顯著。

3結(jié)束語

本文在多項式模型研究的基礎(chǔ)上,提出了一種改進的衛(wèi)星鐘差預(yù)報方法,該模型充分考慮了星載原子鐘的物理特性,并重點對多項式模型的選擇、周期項修正、殘差隨機項修正等方面進行了改進。試驗表明,改進的鐘差預(yù)報模型能夠顯著地提高短期鐘差的預(yù)報精度,對于穩(wěn)定度較差的衛(wèi)星的預(yù)報效果尤為顯著。改進算法的衛(wèi)星鐘預(yù)報6 h預(yù)報精度總體優(yōu)于1 ns,24 h預(yù)報精度優(yōu)于2 ns.在實際應(yīng)用中,可依據(jù)具體情況對模型進行擴充或簡化。

文中僅研究了改進算法對GPS衛(wèi)星的短期預(yù)報效果。在后續(xù)工作中,探索衛(wèi)星鐘差的中長期變化規(guī)律,并將算法應(yīng)用于BDS衛(wèi)星鐘差的預(yù)報工作,是下一步的研究重點。

參考文獻

[1] IGS. IGS產(chǎn)品[EB/OL].http://igsws.unavco.org/components/prods.html.

[2] 付文舉.GNSS在軌衛(wèi)星鐘特性分析及鐘差預(yù)報研究[D].西安:長安大學(xué),2014.

[3] 鄭作亞,盧秀山.幾種GPS衛(wèi)星鐘差預(yù)報方法比較及精度分析[J].山東科技大學(xué)學(xué)報·自然科學(xué)版,2008,27(4):6-11.

[4] 黃觀文.基于頻譜分析的IGS精密星歷衛(wèi)星鐘差精度分析研究[J]. 武漢大學(xué)學(xué)報·信息科學(xué)版,2008,33(5):497-499.

[5] 范旭亮,王曉紅，張顯云，等.基于ARIMA模型的衛(wèi)星鐘差短期預(yù)報研究[J].測繪與空間地理信息,2015,38(1):105-110.

[6] 崔先強, 焦文海. 灰色系統(tǒng)模型在衛(wèi)星鐘差預(yù)報中的應(yīng)用[J]. 武漢大學(xué)學(xué)報·信息科學(xué)版,2005, 30(5):447-450.

[7] 鄭作亞,黨亞民，盧秀山，等.附有周期項的預(yù)報模型及其在GPS衛(wèi)星鐘差預(yù)報中的應(yīng)用研究[J].天文學(xué)報，2010,51(1):95-102.

[8] 唐桂芬,許雪晴，曹紀東，等?；谕ㄓ苗姴钅Ｐ偷谋倍风婎A(yù)報精度分析[J].中國科學(xué):物理學(xué) 力學(xué) 天文學(xué)，2015,45(7):1-6.

[9] GREENGARD L,LEE J. Accelerating the nonuniform fast Fourier transform[J].Siam Review,2004, 46(3):443-454.

[10]LEE J-Y,GREENGARF L. The type 3 nonuniform FFT and its application[J].Journal of Computational Physics, 2005,206(1):1-5.

[11]HUANG Guanwen, ZHANG Qin,XU Guuochang,etal. Real-time clock offset prediction with an improved model[J].GPS Solutions,2014,18(1):95-104.

[12]IGS.IGS_CB分析中心[EB/OL].http://igscb.jpl.nasa.gov/pipermail/igsmail/2000/003036.html.

馬卓希(1991-),男,碩士生,主要從事GNSS精密鐘差解算及其預(yù)報方面的研究。

楊力(1965-),男,博士,教授,主要從事衛(wèi)星導(dǎo)航定位等教學(xué)研究。

賈小林(1972-),男,博士,研究員,主要從事導(dǎo)航衛(wèi)星定軌方面的研究。

趙海山(1991-),男,碩士生,主要從事GNSS理論與數(shù)據(jù)處理方面的研究。

The Optimization Method for the Satellite Clock Error Prediction Based on Polynomial Model

MA Zhuoxi1,YANG Li1,JIA Xiaolin2,ZHAO Haishan1

(1.InstituteofNavigationandSpaceTargetEngineering,InformationEngineeringUniversity,Zhengzhou450001,China;2.Xi'anInstituteofSurveyingandmapping,Xi'an710054,China)

Abstract: According to the physical characteristics of satellite clock, this paper adopts a universal clock offset model based on polynomial prediction model. Multiple improvements have been verified,which contains the self-adaptive selection of polynomial model, modeling for periodic term and modeling for the residual term. Firstly, the adaptive selection of polynomial fitting model is judged by the standard deviation of residual. Then the periodical term of clock offset is modeled using spectrum analysis. Finally the stochastic model after removing the periodical term is constructed using AR model. The experiment results show that the proposed model can greatly improve the short-term prediction accuracy, especially for less stable satellite clock. The overall accuracy of 6 h prediction can be better than 1 ns and 24 h prediction accuracy better than 2 ns.

Keywords:GPS; satellite atomic clock; Polynomial model; clock offset prediction

doi:10.13442/j.gnss.1008-9268.2016.02.005

收稿日期：2015-11-02

中圖分類號：P228.4

文獻標志碼：A

文章編號：1008-9268(2016)02-0027-06

作者簡介

聯(lián)系人： 馬卓希 E-mail: mzhx-1991@163.com