附加原子鐘物理模型的PPP時(shí)間傳遞算法

于合理，郝金明，劉偉平，田英國(guó)，鄧 科

1. 信息工程大學(xué)導(dǎo)航與空天目標(biāo)工程學(xué)院，河南 鄭州 450001； 2. 北斗導(dǎo)航應(yīng)用技術(shù)河南省協(xié)同創(chuàng)新中心，河南 鄭州 450001

?

附加原子鐘物理模型的PPP時(shí)間傳遞算法

于合理1,2，郝金明1,2，劉偉平1,2，田英國(guó)1,2，鄧 科1

1. 信息工程大學(xué)導(dǎo)航與空天目標(biāo)工程學(xué)院，河南 鄭州 450001； 2. 北斗導(dǎo)航應(yīng)用技術(shù)河南省協(xié)同創(chuàng)新中心，河南 鄭州 450001

傳統(tǒng)精密單點(diǎn)定位(PPP)時(shí)間傳遞算法通常把接收機(jī)鐘差當(dāng)作相互獨(dú)立的白噪聲逐歷元進(jìn)行估計(jì)，而忽略了鐘差參數(shù)歷元間的相關(guān)性。針對(duì)這一問(wèn)題，本文提出了一種附加原子鐘物理模型的PPP時(shí)間傳遞算法。該算法通過(guò)利用Kalman濾波對(duì)高穩(wěn)定度的原子鐘鐘差進(jìn)行建模，拓展傳統(tǒng)PPP時(shí)間傳遞模型中的接收機(jī)鐘差參數(shù)，并給出了Kalman濾波過(guò)程噪聲協(xié)方差和初始狀態(tài)向量的確定方法。試驗(yàn)結(jié)果表明：該算法可以有效避免傳統(tǒng)算法時(shí)間傳遞結(jié)果需要一定收斂時(shí)間的問(wèn)題，使解算結(jié)果更加符合原子鐘的物理特性，能夠顯著提高時(shí)間傳遞結(jié)果的精度和穩(wěn)定性，可將單站時(shí)間傳遞精度平均提高58%，站間時(shí)間傳遞精度平均提高51%。

原子鐘；物理模型；精密單點(diǎn)定位；時(shí)間傳遞

在高精度遠(yuǎn)程時(shí)間傳遞中衛(wèi)星雙向時(shí)間頻率傳遞(two way satellite time and frequency transfer，TWSTFT)和精密單點(diǎn)定位(precise point positioning，PPP)時(shí)間傳遞技術(shù)是精度最高的兩種方法[1-3]。TWSTFT需要專門的中介衛(wèi)星建立通信鏈路，且參與雙方都需要專門的發(fā)射接收設(shè)備，設(shè)備價(jià)格昂貴，目前世界上只有少數(shù)較大的實(shí)驗(yàn)室應(yīng)用。而PPP時(shí)間傳遞算法是單站作業(yè)，操作靈活方便，時(shí)間傳遞的精度可達(dá)亞納秒量級(jí)，其短期穩(wěn)定度甚至優(yōu)于TWSTFT技術(shù)[1-3]，已經(jīng)被廣泛應(yīng)用于時(shí)間傳遞領(lǐng)域。隨著國(guó)際GNSS服務(wù)組織(International GNSS Services,IGS)和國(guó)際時(shí)頻機(jī)構(gòu)的深入合作，許多實(shí)驗(yàn)室已經(jīng)將PPP方法用于時(shí)間傳遞和原子時(shí)的計(jì)算[3]，國(guó)際計(jì)量局還將其用于維護(hù)協(xié)調(diào)世界時(shí)的運(yùn)作[1]。文獻(xiàn)[1—5]對(duì)PPP時(shí)間傳遞技術(shù)進(jìn)行了研究，并取得了良好效果。文獻(xiàn)[6—7]還對(duì)非組合PPP時(shí)間傳遞方法進(jìn)行了研究，該方法能有效避免傳統(tǒng)PPP算法由于觀測(cè)值組合過(guò)程中觀測(cè)噪聲及多路徑效應(yīng)誤差被放大而造成的對(duì)接收機(jī)鐘差解算的不利影響，解算結(jié)果優(yōu)于傳統(tǒng)PPP算法。傳統(tǒng)的PPP時(shí)間傳遞算法通常是將各歷元接收機(jī)鐘差看作相互獨(dú)立的白噪聲逐歷元進(jìn)行估計(jì)，未充分發(fā)掘接收機(jī)鐘差參數(shù)在相鄰歷元間變化可能較為平穩(wěn)的約束信息。同時(shí)由于接收機(jī)鐘差、測(cè)站高程和對(duì)流層參數(shù)的高度相關(guān)，解算得到的鐘差序列往往會(huì)偏離接收機(jī)鐘差序列多項(xiàng)式的特性[8-9]，時(shí)間傳遞結(jié)果還會(huì)存在一定時(shí)間的收斂過(guò)程。為此，本文利用原子鐘的物理特性對(duì)接收機(jī)鐘差建模來(lái)進(jìn)一步提高PPP時(shí)間傳遞的精度和穩(wěn)定性。

利用接收機(jī)鐘物理特性對(duì)鐘差進(jìn)行建模的思想并不是最新提出的，美國(guó)噴氣推進(jìn)實(shí)驗(yàn)室(Jet Propulsion Laboratory，JPL)在20世紀(jì)80年代為了提高GNSS軌道等相關(guān)參數(shù)的精度就采用不同的隨機(jī)過(guò)程模型來(lái)對(duì)接收機(jī)鐘差和對(duì)流層進(jìn)行建模[10]。但由于早期的接收機(jī)一般采用穩(wěn)定度較差的石英鐘，即使采用復(fù)雜的數(shù)學(xué)模型也無(wú)法準(zhǔn)確描述石英鐘的運(yùn)行特征，鐘差建模對(duì)解算結(jié)果的影響微乎其微。而目前許多時(shí)間實(shí)驗(yàn)室都已經(jīng)配備了高穩(wěn)定度的原子鐘，這使通過(guò)鐘差物理建模提高時(shí)間傳遞精度變得可能。近年來(lái)眾多學(xué)者對(duì)鐘差建模進(jìn)行了大量的研究，取得了豐碩的成果。在衛(wèi)星原子鐘建模方面，文獻(xiàn)[11]在顧及原子鐘物理模型的基礎(chǔ)上提出了適用于原子鐘的Kalman濾波器；文獻(xiàn)[12—15]進(jìn)一步研究了原子鐘濾波過(guò)程中狀態(tài)矩陣和觀測(cè)噪聲矩陣的確定方法；文獻(xiàn)[16]利用原子鐘物理特性對(duì)鐘差建模實(shí)現(xiàn)了對(duì)銣原子鐘的狀態(tài)監(jiān)測(cè)。在接收機(jī)鐘差建模方面，文獻(xiàn)[17—19]對(duì)低軌衛(wèi)星星載接收機(jī)端的鐘差建模方法進(jìn)行了研究；文獻(xiàn)[9]對(duì)高穩(wěn)定性地面接收機(jī)鐘的隨機(jī)模型進(jìn)行了研究；文獻(xiàn)[8，20—21]對(duì)通過(guò)接收機(jī)鐘物理建模來(lái)提高定位精度和收斂速度進(jìn)行了研究，并取得了一定的效果。目前針對(duì)接收機(jī)鐘差建模提高時(shí)間傳遞性能方面的研究較少，且處于起步階段。文獻(xiàn)[22]通過(guò)接收機(jī)鐘差建模對(duì)接收機(jī)鐘差參數(shù)估計(jì)進(jìn)行了一些研究，但主要是針對(duì)石英晶振接收機(jī)，并不能有效提高時(shí)間傳遞精度。文獻(xiàn)[1]也對(duì)基于原子鐘鐘差模型的時(shí)間傳遞算法進(jìn)行了初步的研究。

因此，本文利用原子鐘的物理特性，采用Kalman濾波方法對(duì)接收機(jī)鐘差進(jìn)行建模，給出了Kalman濾波過(guò)程噪聲協(xié)方差矩陣、初始狀態(tài)向量及其協(xié)方差矩陣的確定方法，并重點(diǎn)分析了附加原子鐘物理模型對(duì)單站時(shí)間傳遞及站間時(shí)間同步結(jié)果的影響。

1 附加原子鐘物理模型的PPP時(shí)間傳遞算法

1.1 基本觀測(cè)模型

對(duì)于雙頻接收機(jī)，PPP時(shí)間傳遞算法一般采用消電離層組合觀測(cè)量[1]，觀測(cè)方程如下

(1)

(2)

1.2 附加原子鐘物理模型

原子鐘具有較高的頻率穩(wěn)定度，其鐘差通常采用簡(jiǎn)單的一次多項(xiàng)式或二次多項(xiàng)式進(jìn)行描述[23]。限于篇幅，本文僅以二次多項(xiàng)式為例進(jìn)行分析

x(t)=a0+a1(t-t0)+a2(t-t0)2+εx(t)

(3)

式中，x(t)為t時(shí)刻的接收機(jī)鐘差；a0、a1、a2為確定性分量，分別表示參考時(shí)刻t0的初始鐘差、鐘速、鐘漂，當(dāng)a2=0時(shí)，式(3)可簡(jiǎn)化為線性模型；εx(t)為隨機(jī)變化的不確定性分量。

PPP時(shí)間傳遞算法通常采用Kalman濾波進(jìn)行參數(shù)估計(jì)，一般是將接收機(jī)鐘差當(dāng)作獨(dú)立的白噪聲，逐歷元估計(jì)鐘差參數(shù)，未充分發(fā)掘接收機(jī)鐘差參數(shù)在相鄰歷元間變化可能較為平穩(wěn)的約束信息。若考慮接收機(jī)鐘差歷元間的關(guān)系，可附加原子鐘的物理模型進(jìn)行約束，拓展傳統(tǒng)PPP模型中的接收機(jī)鐘差參數(shù)，使解算的鐘差序列仍然體現(xiàn)原子鐘的物理特性。具體做法是在PPP時(shí)間傳遞基本觀測(cè)方程式(1)與式(2)估計(jì)接收機(jī)鐘差參數(shù)的基礎(chǔ)上，同時(shí)估計(jì)鐘速和鐘漂參數(shù)。此時(shí)觀測(cè)方程中未知參數(shù)為測(cè)站三維坐標(biāo)參數(shù)、鐘差、鐘速、鐘漂、對(duì)流層和模糊度參數(shù)。為了便于推導(dǎo)和直觀地顯示附加鐘差模型的方法，公式推導(dǎo)過(guò)程僅考慮鐘差模型參數(shù)，觀測(cè)方程簡(jiǎn)化為

Lk=AkXk+Δk

(4)

式中，Lk為觀測(cè)量；Ak=[100]為設(shè)計(jì)矩陣；Xk=[x(t+τ)y(t+τ)z(t+τ)]T為tk時(shí)刻的三維狀態(tài)向量，x(t)、y(t)、z(t)為接收機(jī)鐘差(時(shí)差)、鐘速(頻率)和鐘漂(頻漂)，τ為tk時(shí)刻與tk-1時(shí)刻的時(shí)間間隔；Δk為觀測(cè)噪聲。

Kalman濾波三維狀態(tài)空間方程為

(5)

令

(6)

可得狀態(tài)空間方程的簡(jiǎn)化形式

Xk=Φk,k-1Xk-1+Wk

(7)

式中，Φk,k-1為狀態(tài)轉(zhuǎn)移矩陣；Wk為三維過(guò)程噪聲向量，其協(xié)方差矩陣Q可以表示成過(guò)程噪聲參數(shù)的函數(shù)[13]

(8)

式中，q1、q2和q3分別可用調(diào)頻白噪聲、調(diào)頻隨機(jī)游走噪聲和調(diào)頻隨機(jī)奔跑噪聲來(lái)描述。

若原子鐘線性頻漂不明顯，可采用二維狀態(tài)模型來(lái)描述接收機(jī)的動(dòng)態(tài)過(guò)程，此時(shí)狀態(tài)轉(zhuǎn)移矩陣和協(xié)方差陣需要根據(jù)三維狀態(tài)模型進(jìn)行相應(yīng)的調(diào)整。

2 過(guò)程噪聲協(xié)方差矩陣的確定

系統(tǒng)過(guò)程噪聲協(xié)方差矩陣的確定是Kalman濾波的關(guān)鍵，可由原子鐘的濾波過(guò)程噪聲參數(shù)來(lái)表示。過(guò)程噪聲參數(shù)可以通過(guò)哈達(dá)瑪總方差或阿倫方差進(jìn)行估計(jì)。

2.1 基于哈達(dá)瑪總方差的過(guò)程噪聲估計(jì)

原子鐘頻漂較明顯時(shí)，可采用哈達(dá)瑪總方差來(lái)估計(jì)過(guò)程噪聲參數(shù)。基于時(shí)差數(shù)據(jù)的哈達(dá)瑪總方差計(jì)算公式為[5,24]

(9)

當(dāng)同時(shí)考慮調(diào)相白噪聲和狀態(tài)方程中存在的3種調(diào)頻噪聲，且假設(shè)調(diào)相白噪聲與其他3種噪聲不相關(guān)時(shí)，可推導(dǎo)出公式[5]

(10)

式中，q0為對(duì)應(yīng)于調(diào)相白噪聲的測(cè)量噪聲參數(shù)。

2.2 基于阿倫方差的過(guò)程噪聲估計(jì)

原子鐘頻漂不明顯時(shí)，可以采用阿倫方差來(lái)估計(jì)過(guò)程噪聲參數(shù)。基于時(shí)差數(shù)據(jù)阿倫方差計(jì)算公式為[5,24-25]

(11)

當(dāng)同時(shí)考慮調(diào)相白噪聲和狀態(tài)方程中存在的3種調(diào)頻噪聲時(shí)，可推導(dǎo)出公式[5]

(12)

設(shè)定不同的采樣間隔，就可以得到不同的哈達(dá)瑪總方差(或阿倫方差)公式，聯(lián)立各式利用最小二乘算法即可解得過(guò)程噪聲參數(shù)q0、q1、q2和q3，進(jìn)而確定系統(tǒng)過(guò)程噪聲協(xié)方差矩陣。

3 初始狀態(tài)向量及協(xié)方差確定

Kalman濾波算法需要給定初始狀態(tài)向量及其協(xié)方差矩陣。在無(wú)先驗(yàn)信息時(shí)，通常令初始狀態(tài)向量為零向量，協(xié)方差矩陣為aI，其中a為很大的正數(shù)，I為單位矩陣，以確保濾波結(jié)果的無(wú)偏性。但在沒(méi)有先驗(yàn)信息時(shí)，濾波器初始化會(huì)比較困難，且會(huì)影響濾波器的穩(wěn)定性。在進(jìn)行連續(xù)解算或已知先驗(yàn)信息時(shí)，若已知3個(gè)時(shí)差測(cè)量值l(0)、l(1)、l(2)，τ1、τ2為這3個(gè)觀測(cè)量的采樣間隔，則根據(jù)式(13)—式(21)即可確定狀態(tài)向量初值X(2)=[x(2)y(2)z(2)]以及協(xié)方差矩陣P(2)，其中協(xié)方差矩陣為對(duì)角陣[15-16]

x(2)=l(2)

(13)

(14)

(15)

P11(2)=q0

(16)

(17)

(18)

(19)

(20)

(21)

4 試驗(yàn)分析

由JPL發(fā)布的“l(fā)oglist.txt”文件可知，截至2016年6月2日，IGS已經(jīng)有超過(guò)130個(gè)跟蹤站配備了高頻穩(wěn)度的原子鐘，其中71個(gè)配備了氫原子鐘[26]，這為本文算法性能的驗(yàn)證提供了很大的便利。為了驗(yàn)證附加原子鐘物理模型的PPP時(shí)間傳遞算法的應(yīng)用效果，選取2015年1月5日配備了氫原子鐘的WAB2、CRO1、PTBB、NRC1和MDVJ 5個(gè)IGS跟蹤站的觀測(cè)數(shù)據(jù)進(jìn)行試驗(yàn)分析。觀測(cè)數(shù)據(jù)采樣間隔為30 s，采用偽距和載波相位消電離層組合觀測(cè)量，偽距與載波相位的權(quán)比為1/10 000，高度截止角設(shè)為10°。試驗(yàn)采用靜態(tài)模式進(jìn)行處理，利用卡爾曼濾波估計(jì)三維坐標(biāo)、接收機(jī)鐘差、鐘速、對(duì)流層延遲和模糊度參數(shù)。衛(wèi)星軌道和鐘差采用IGS提供的事后精密星歷和鐘差產(chǎn)品，為了避免衛(wèi)星和接收機(jī)硬件延遲偏差對(duì)接收機(jī)鐘差造成的系統(tǒng)性偏差影響，利用CODE發(fā)布的DCB產(chǎn)品對(duì)其進(jìn)行改正[27]。分別采用以下兩種方案進(jìn)行PPP解算。

方案1：將接收機(jī)鐘差當(dāng)作白噪聲，采用傳統(tǒng)的逐歷元估計(jì)一維鐘差參數(shù)的方法。

方案2：顧及接收機(jī)鐘差歷元間的相關(guān)性，采用附加原子鐘物理模型的方法進(jìn)行解算，過(guò)程噪聲由阿倫方差根據(jù)IGS事后精密產(chǎn)品解算得到。

目前IGS提供的事后接收機(jī)鐘差產(chǎn)品精度為75 ps，可將其作為真值進(jìn)行比對(duì)。由于IGS事后精密鐘差產(chǎn)品的采樣間隔為5 min，因此本文以5 min為采樣間隔對(duì)方案1和方案2的時(shí)間傳遞精度進(jìn)行分析。限于篇幅，僅給出了WAB2和CRO1兩跟蹤站采用兩種方案得到的結(jié)果與IGS時(shí)間序列對(duì)比，分別如圖1與圖2所示。同時(shí)將兩種方案的鐘差結(jié)果分別與IGS最終結(jié)果作差可以得到對(duì)應(yīng)的差值序列，圖3與圖4分別給出了WAB2和CRO1兩跟蹤站單站時(shí)間傳遞結(jié)果的差值序列。

從圖1—圖4可知：

(1) 方案2的收斂速度明顯快于方案1，單歷元就可以達(dá)到很高的精度。而方案1由于未使用原子鐘物理模型作為先驗(yàn)信息進(jìn)行約束，需要一定時(shí)間的收斂過(guò)程才能達(dá)到較高的精度，WAB2和CRO1兩測(cè)站分別需要22 min、35 min才能收斂。另外，對(duì)方案2得到的坐標(biāo)，對(duì)流層和模糊度結(jié)果進(jìn)行分析，其收斂速度均有一定程度的提高。

圖1 WAB2-GPST兩種方案結(jié)果與IGS時(shí)間序列對(duì)比Fig.1 WAB2-GPST results of the two schemes compared with IGS time series

圖2 CRO1-GPST兩種方案結(jié)果與IGS時(shí)間序列對(duì)比Fig.2 CRO1-GPST results of the two schemes compared with IGS time series

圖3 WAB2-GPST兩種方案時(shí)間傳遞結(jié)果與IGS差值序列Fig.3 Difference between WAB2-GPST results and IGS from two schemes

(2) 方案2由于附加了原子鐘物理模型，得到的接收機(jī)鐘差時(shí)間序列的連續(xù)性和穩(wěn)定性明顯優(yōu)于方案1，與IGS發(fā)布的事后精密鐘差結(jié)果的一致性較高，更能準(zhǔn)確地描述接收機(jī)鐘差序列的變化規(guī)律。

圖4 CRO1-GPST兩種方案時(shí)間傳遞結(jié)果與IGS差值序列Fig.4 Difference between CRO1-GPST results and IGS from two schemes

(3) 從圖1明顯可以看出，方案1解算的鐘差序列略微偏離了接收機(jī)鐘本身的多項(xiàng)式特性，這是由于方案1解算的鐘差序列之間是相互獨(dú)立的，解算的過(guò)程中吸收了許多不屬于鐘差的其他誤差源而造成。方案2的濾波結(jié)果更能準(zhǔn)確可靠地反映原子鐘的物理特性。

表1中給出了WAB2、CRO1、PTBB、NRC1和MDVJ5個(gè)跟蹤站相應(yīng)的誤差統(tǒng)計(jì)結(jié)果。表中SD為標(biāo)準(zhǔn)差(standarddeviation)，用來(lái)衡量鐘差誤差序列自身的離散程度，RMS為均方根誤差用來(lái)衡量鐘差估計(jì)值同真值之間的誤差。

表1 兩種方案結(jié)果時(shí)間傳遞精度

Tab.1 Accuracy of the time transferring results from different schemes ns

由表1可知，各測(cè)站方案2時(shí)間傳遞結(jié)果的SD和RMS均優(yōu)于方案1，方案1的鐘差解相對(duì)較為離散，方案2的鐘差解穩(wěn)定度較高，與IGS結(jié)果的符合程度更高。方案2附加原子鐘物理模型后，能夠顯著提高PPP時(shí)間傳遞的精度，可將WAB2、CRO1、PTBB、NRC1和MDVJ 5個(gè)跟蹤站的單站時(shí)間傳遞精度(RMS)分別提高71%、62%、73%、53%和33%，平均提高58%。該算法對(duì)于促進(jìn)PPP時(shí)間傳遞技術(shù)在國(guó)際時(shí)間實(shí)驗(yàn)室等高精度時(shí)間同步領(lǐng)域的應(yīng)用具有一定意義。

為了進(jìn)一步分析本文算法應(yīng)用于站間時(shí)間同步中的精度，將各站利用PPP獲得的原子鐘的鐘差結(jié)果相減即可得到站間時(shí)間同步結(jié)果，表2給出了兩種方案得到的不同時(shí)間鏈路的時(shí)間同步序列的精度統(tǒng)計(jì)結(jié)果，為更好地觀察時(shí)間同步序列變化情況，圖5以CRO1-WAB2時(shí)間鏈路為例給出了站間時(shí)間同步結(jié)果同IGS結(jié)果的差值序列圖。

表2 兩種方案時(shí)間同步結(jié)果精度

Tab.2 Accuracy of the time transferring results from different schemes ns

圖5 CRO1-WAB2鏈路兩方案時(shí)間同步結(jié)果與IGS結(jié)果差值序列Fig.5 Difference between CRO1-WAB2 results and IGS from two schemes

根據(jù)表2和圖5可知，方案2得到的站間時(shí)間同步結(jié)果的精度明顯優(yōu)于方案1，與IGS時(shí)間同步結(jié)果一致性較高，能夠可靠地描述測(cè)站間時(shí)間同步序列的變化規(guī)律。由圖5可知，方案1采用傳統(tǒng)算法CRO1-WAB2鏈路需要大約35 min才能收斂，而方案2附加原子鐘物理模型后幾乎不需要收斂時(shí)間即可達(dá)到很高的精度。方案2相比方案1，能夠?qū)RO1-WAB2、CRO1-PTBB、CRO1-NRC1和CRO1-MDVJ鏈路的站間時(shí)間同步精度分別提高51%、49%、29%和73%，平均提高51%。

5 結(jié) 論

鑒于當(dāng)前許多國(guó)際時(shí)頻實(shí)驗(yàn)室和IGS跟蹤站均配置有高性能原子鐘，本文通過(guò)附加原子鐘物理模型對(duì)現(xiàn)有PPP時(shí)間傳遞模型進(jìn)行改進(jìn)，使用Kalman濾波對(duì)接收機(jī)原子鐘進(jìn)行鐘差建模，研究了Kalman濾波中系統(tǒng)過(guò)程噪聲協(xié)方差矩陣、初始狀態(tài)向量及協(xié)方差的確定方法，并重點(diǎn)分析了附加原子鐘物理模型對(duì)單站時(shí)間傳遞及站間時(shí)間同步結(jié)果的影響。試驗(yàn)結(jié)果表明：附加原子鐘物理模型后，解決了傳統(tǒng)算法中收斂過(guò)程的耗時(shí)問(wèn)題，顯著提高單站時(shí)間傳遞和站間時(shí)間同步的精度，可將單站時(shí)間傳遞精度平均提高58%，尤其是WAB2和PTBB站精度提高水平高達(dá)71%和73%，同時(shí)可將站間時(shí)間同步精度平均提高51%。研究附加原子鐘物理模型的PPP時(shí)間傳遞算法，對(duì)于國(guó)際時(shí)間實(shí)驗(yàn)室實(shí)現(xiàn)高精度的時(shí)間傳遞具有一定的參考意義，對(duì)于拓展GNSS在時(shí)頻傳遞領(lǐng)域的應(yīng)用具有重要意義。本文算法并非適用于所有的GNSS接收機(jī)，其前提是要判定接收機(jī)鐘本身的物理特性。它要求接收機(jī)配備有高精度和高穩(wěn)定的原子鐘，否則加入鐘差物理模型后不但不能提高時(shí)間傳遞精度，反而會(huì)引起比較大的誤差。

[1] 李紅濤. 基于GPS和GLONASS的單站授時(shí)和時(shí)差監(jiān)測(cè)研究[D]. 西安: 長(zhǎng)安大學(xué), 2012.

LI Hongtao. Research on the Single Station Time Service and Time Offset Based on the GPS and GLONASS Data[D]. Xi’an: Chang’an University, 2012.

[2] DEFRAIGNE P, AERTS W, HARMEGNIES A, et al. Advances in Multi-GNSS Time Transfer[C]∥2013 Joint European Frequency and Time Forum & International Frequency Control Symposium (EFTF/IFC). Prague, Czech Republic: IEEE, 2013: 508-512.

[3] DEFRAIGNE P, AERTS W, POTTIAUX E. Monitoring of UTC(k)’s Using PPP and IGS Real-time Products[J]. GPS Solutions, 2015, 19(1): 165-172.

[4] 張小紅, 蔡詩(shī)響, 李星星, 等. 利用GPS精密單點(diǎn)定位進(jìn)行時(shí)間傳遞精度分析[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2010, 35(3): 274-278.

ZHANG Xiaohong, CAI Shixiang, LI Xingxing, et al. Accuracy Analysis of Time and Frequency Transfer Based on Precise Point Positioning[J]. Geomatics and Information Science of Wuhan University, 2010, 35(3): 274-278.

[5] 黃觀文. GNSS星載原子鐘質(zhì)量評(píng)價(jià)及精密鐘差算法研究[D]. 西安: 長(zhǎng)安大學(xué), 2012.

HUANG Guanwen. Research on Algorithms of Precise Clock Offset and Quality Evaluation of GNSS Satellite Clock[D]. Xi’an: Chang’an University, 2012.

[6] 閆偉, 袁運(yùn)斌, 歐吉坤, 等. 非組合精密單點(diǎn)定位算法精密授時(shí)的可行性研究[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2011, 36(6): 648-651.

YAN Wei, YUAN Yunbin, OU Jikun, et al. Feasibility of Precise Timing with Uncombined PPP[J]. Geomatics and Information Science of Wuhan University, 2011, 36(6): 648-651.

[7] 張寶成, 歐吉坤, 袁運(yùn)斌, 等. 基于GPS雙頻原始觀測(cè)值的精密單點(diǎn)定位算法及應(yīng)用[J]. 測(cè)繪學(xué)報(bào), 2010, 39(5): 478-483.

ZHANG Baocheng, OU Jikun, YUAN Yunbin, et al. Precise Point Positioning Algorithm Based on Original Dual-frequency GPS Code and Carrier-phase Observations and Its Application[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(5): 478-483.

[8] 張小紅, 陳興漢, 郭斐. 高性能原子鐘鐘差建模及其在精密單點(diǎn)定位中的應(yīng)用[J]. 測(cè)繪學(xué)報(bào), 2015, 44(4): 392-398. DOI: 10.11947/j.AGCS.2015.20140287.ZHANG Xiaohong, CHEN Xinghan, GUO Fei. High-performance Atomic Clock Modeling and Its Application in Precise Point Positioning[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(4): 392-398. DOI: 10.11947/j.AGCS.2015.20140287.

[9] WANG Kan, ROTHACHER M. Stochastic Modeling of High-stability Ground Clocks in GPS Analysis[J]. Journal of Geodesy, 2013, 87(5): 427-437.

[10] LICHTEN S M, BORDER J S. Strategies for High-precision Global Positioning System Orbit Determination[J]. Journal of Geophysical Research: Solid Earth, 1987, 92(B12): 12751-12762.

[11] JONES R H, TRYON P V. Continuous Time Series Models for Unequally Spaced Data Applied to Modeling Atomic Clocks[J]. SIAM Journal on Scientific and Statistical Computing, 1987, 8(1): 71-81.

[12] 林旭, 羅志才. 一種新的衛(wèi)星鐘差Kalman濾波噪聲協(xié)方差估計(jì)方法[J]. 物理學(xué)報(bào), 2015, 64(8): 080201.

LIN Xu, LUO Zhicai. A New Noise Covariance Matrix Estimation Method of Kalman Filter for Satellite Clock Errors[J]. Acta Physica Sinica, 2015, 64(8): 080201.

[13] 郭海榮, 楊元喜, 何海波, 等. 導(dǎo)航衛(wèi)星原子鐘Kalman濾波中噪聲方差-協(xié)方差的確定[J]. 測(cè)繪學(xué)報(bào), 2010, 39(2): 146-150.

GUO Hairong, YANG Yuanxi, HE Haibo, et al. Determination of Covariance Matrix of Kalman Filter Used for Time Prediction of Atomic Clocks of Navigation Satellites[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(2): 146-150.

[14] SU W, FILLER R L. Application of Kalman Filtering Techniques to the Precision Clock with Non-constant Aging[C]∥Proceedings of the 46th Frequency Control Symposium. Hershey, PA: IEEE, 1992: 231-237.

[15] STEIN S R, FILLER R L. Kalman Filter Analysis for Real Time Applications of Clocks and Oscillators[C]∥Proceedings of the 42nd Annual Frequency Control Symposium. Baltimore, MD: IEEE, 1988: 447-452.

[16] 張清華, 隋立芬, 賈小林. 應(yīng)用Jones-Tryon Kalman濾波器對(duì)在軌GPS Rb鐘進(jìn)行狀態(tài)監(jiān)測(cè)[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2012, 37(4): 436-440.

ZHANG Qinghua, SUI Lifen, JIA Xiaolin. Monitor State of GPS Rb Clock Using Jones-Tryon Kalman Filter[J]. Geomatics and Information Science of Wuhan University, 2012, 37(4): 436-440.

[17] YANG Yang, YUE Xiaokui, YUAN Jianping, et al. Enhancing the Kinematic Precise Orbit Determination of Low Earth Orbiters Using GPS Receiver Clock Modelling[J]. Advances in Space Research, 2014, 54(9): 1901-1912.

[18] WEINBACH U, SCH?N S. Improved GRACE Kinematic Orbit Determination Using GPS Receiver Clock Modeling[J]. GPS Solutions, 2013, 17(4): 511-520.

[19] WEINBACH U, SCHON S. Improved GPS Receiver Clock Modeling for Kinematic Orbit Determination of the GRACE Satellites[C]∥Proceedings of European Frequency and Time Forum. Gothenburg, Sweden: IEEE, 2012: 157-160.

[20] 肖國(guó)銳, 隋立芬, 陳泉余, 等. 利用接收機(jī)鐘差建模提升PPP收斂速度及精度[J]. 測(cè)繪科學(xué)技術(shù)學(xué)報(bào), 2015, 32(6): 555-558, 564.XIAO Guorui, SUI Lifen, CHEN Quanyu, et al. Improving PPP Convergence and Accuracy Using Receiver Clock Modeling[J]. Journal of Geomatics Science and Technology, 2015, 32(6): 555-558, 564

[21] WEINBACH U, SCH?N S. GNSS Receiver Clock Modeling When Using High-precision Oscillators and Its Impact on PPP[J]. Advances in Space Research, 2011, 47(2): 229-238.

[22] FILHO E A M, KUGA H K, LOPES R V F. Real Time Estimation of GPS Receiver Clock Offset by the Kalman Filter[J]. Personal Communication, 2003.

[23] 于合理, 郝金明, 劉偉平, 等. 一種衛(wèi)星鐘差異常實(shí)時(shí)監(jiān)測(cè)算法[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2016, 41(1): 106-110.

YU Heli, HAO Jinming, LIU Weiping, et al. A Real-time Anomaly Monitoring Algorithm for Satellite Clock[J]. Geomatics and Information Science of Wuhan University, 2016, 41(1): 106-110.

[24] RILEY W J. Handbook of Frequency Stability Analysis[R]. NIST Special Publication 1065, 2007: 1-123.

[25] MANNING D M. AF/NGA GPS Monitor Station High-performance Cesium Frequency Standard Stability 2005/2006: from NGA Kalman Filter Clock Estimates[C]∥Proceedings of the 38th Annual Precise Time and Time Interval Systems and Applications Meeting. Reston, Virginia: [s.n.], 2006: 137-152.

[26] CLKLOG. A Summary File of the Deployment History for GPS Receiver, Antenna, Frequency Standards, and Other Equipment at IGS Stations[EB/OL]. (2010-01-20) [2016-06-02]. ftp:∥igscb.jpl.nasa.gov/igscb/station/general/loghist.txt.

[27] 于合理, 郝金明, 謝建濤, 等. 硬件延遲偏差對(duì)衛(wèi)星鐘差解算的影響[J]. 導(dǎo)航定位學(xué)報(bào), 2015, 3(1): 71-73.

YU Heli, HAO Jinming, XIE Jiantao, et al. Impact of DCB on the Estimate of Satellite Clock Bias[J]. Journal of Navigation and Positioning, 2015, 3(1): 71-73.

(責(zé)任編輯：叢樹平)

A Time Transfer Algorithm of Precise Point Positioning with Additional Atomic Clock Physical Model

YU Heli1,2，HAO Jinming1,2，LIU Weiping1,2，TIAN Yingguo1,2，DENG Ke1

1. School of Navigation and Aerospace Engineering, Information Engineering University, Zhengzhou 450001, China; 2. BeiDou Navigation Technology Collaborative Innovation Center of Henan, Zhengzhou 450001, China

In the traditional time transfer algorithms of precise point positioning, the receiver clock offset is estimated epoch by epoch as independent white noise and the short-term correlation between the epoch of the atomic clock offset is ignored. In order to solve this problem, a time transfer algorithm of precise point positioning with additional atomic clock physical model is proposed. The receiver clock parameters are expanded in the traditional model of precise point positioning by using Kalman filter to model the clock errors of high stability atomic clock, and the method of determining the covariance matrices of process noise and the initial state vector of the Kalman filter is presented. Experimental results show that the proposed algorithm can effectively avoid the convergence process of the traditional precise point positioning algorithm, and the result of the algorithm is more consistent with the physical properties of the atomic clock, this algorithm can significantly improve the accuracy and stability of the time transfer results, the accuracy of time transfer results of single station can be increased on average by 58%, and the accuracy of time synchronization between stations can be improved on average by 51%.

atomic clock; physical model; precise point positioning; time transfer

Specific Project of Chinese Second-generation Satellite Navigation System (No. GFZX0301040308); The Open Research Foundation of State Key Laboratory of Geo-information Engineering (No. SKLGIE2015-M-1-6)

YU Heli(1989—)，male，PhD candidate，majors in GNSS time and frequency transfer.

于合理，郝金明，劉偉平，等.附加原子鐘物理模型的PPP時(shí)間傳遞算法[J].測(cè)繪學(xué)報(bào)，2016,45(11)：1285-1292.

10.11947/j.AGCS.2016.20160217.

YU Heli，HAO Jinming，LIU Weiping，et al.A Time Transfer Algorithm of Precise Point Positioning with Additional Atomic Clock Physical Model[J]. Acta Geodaetica et Cartographica Sinica,2016,45(11):1285-1292. DOI:10.11947/j.AGCS.2016.20160217.

P228

A

1001-1595(2016)11-1285-08

中國(guó)第二代衛(wèi)星導(dǎo)航系統(tǒng)重大專項(xiàng)(GFZX0301040308)；地理信息工程國(guó)家重點(diǎn)實(shí)驗(yàn)室開放研究基金(SKLGIE2015-M-1-6)

2016-05-05

修回日期： 2016-09-05

于合理(1989—)，男，博士生，研究方向?yàn)镚NSS時(shí)頻傳遞方法。

E-mail： yuhl1989@163.com