基于Wiener過(guò)程的GPS校準(zhǔn)晶振型頻率源守頻方法

楊少塵，胡昌華，李紅增

（火箭軍工程大學(xué)控制工程系，陜西 西安 710025）

基于Wiener過(guò)程的GPS校準(zhǔn)晶振型頻率源守頻方法

楊少塵，胡昌華，李紅增

（火箭軍工程大學(xué)控制工程系，陜西 西安 710025）

GPS校準(zhǔn)晶振型頻率源系統(tǒng)接收不到GPS信號(hào)時(shí)，僅能依靠晶振自身保持較高頻率標(biāo)準(zhǔn)輸出。該文針對(duì)短期GPS信號(hào)丟失后晶振輸出頻率漂移問(wèn)題，提出一種基于Wiener過(guò)程的頻率漂移校正補(bǔ)償方法。該方法通過(guò)建立晶振輸出頻率的老化特性退化模型，將溫度作為加速應(yīng)力引入模型參數(shù)。利用GPS信號(hào)丟失前的系統(tǒng)歷史校準(zhǔn)數(shù)據(jù)，采用極大似然估計(jì)法估計(jì)模型擴(kuò)散參數(shù)，采用Kalman濾波方法估計(jì)模型漂移參數(shù)，進(jìn)而預(yù)測(cè)晶振失鎖后的頻率漂移量，完成對(duì)系統(tǒng)輸出頻率的補(bǔ)償校正。將該方法應(yīng)用到某GPS校準(zhǔn)晶振型頻率源系統(tǒng)中，結(jié)果表明在GPS信號(hào)丟失后的短期內(nèi)，系統(tǒng)的輸出頻率準(zhǔn)確度能夠保持在9.7×10-11，滿(mǎn)足短期內(nèi)的計(jì)量需求。

GPS校頻；晶振；頻率漂移；Wiener過(guò)程

0　引　言

隨著航空航天測(cè)控技術(shù)的不斷發(fā)展，對(duì)頻率測(cè)試設(shè)備的計(jì)量需要更高準(zhǔn)確度的時(shí)頻標(biāo)準(zhǔn)。目前高準(zhǔn)確度頻率標(biāo)準(zhǔn)源主要有銣鐘、銫鐘等原子鐘及高準(zhǔn)確度晶振。其中氫、銫等原子鐘的頻率準(zhǔn)確度可達(dá)10-12～10-15量級(jí)，長(zhǎng)期穩(wěn)定度好，但其使用條件嚴(yán)格，價(jià)格昂貴。晶振的頻率準(zhǔn)確度一般為10-7～10-9量級(jí)，短期穩(wěn)定性好，價(jià)格低廉，但其精度易受溫度、老化等因素影響，長(zhǎng)期穩(wěn)定性較差。

GPS系統(tǒng)每顆衛(wèi)星上都攜帶有原子鐘，因此GPS信號(hào)具有原子頻標(biāo)的長(zhǎng)期穩(wěn)定度。利用GPS信號(hào)校準(zhǔn)地面晶振的方法可獲得短期和長(zhǎng)期穩(wěn)定度都很高的頻率標(biāo)準(zhǔn)源，性?xún)r(jià)比高，實(shí)用性強(qiáng)［1-2］。但在地鐵、隧道等場(chǎng)所使用時(shí)，GPS接收機(jī)短期內(nèi)可能接收不到GPS信號(hào)。此時(shí)系統(tǒng)只能依靠晶振自身實(shí)現(xiàn)高精度頻率輸出。但晶振頻率準(zhǔn)確度易受各種因素的影響，其中溫度變化和自身元器件老化是主要的影響因素。為了使系統(tǒng)在丟失GPS信號(hào)后繼續(xù)保持相對(duì)較高的頻率準(zhǔn)確度，需要對(duì)晶振輸出頻率的漂移量進(jìn)行實(shí)時(shí)補(bǔ)償。對(duì)完成航空航天設(shè)備在復(fù)雜環(huán)境下的便捷高效伴隨計(jì)量具有重要意義。

研究晶振頻率漂移預(yù)測(cè)補(bǔ)償?shù)乃悸肥峭ㄟ^(guò)建立晶振性能退化模型，利用GPS信號(hào)正常時(shí)系統(tǒng)的歷史校準(zhǔn)數(shù)據(jù)得到晶振的性能退化數(shù)據(jù)，估計(jì)出模型相關(guān)參數(shù)，預(yù)測(cè)晶振輸出頻率的漂移量，從而對(duì)系統(tǒng)進(jìn)行補(bǔ)償校正，保證系統(tǒng)輸出頻率的精度，直至GPS信號(hào)恢復(fù)后重新鎖定晶振。

建立退化預(yù)測(cè)模型時(shí)，文獻(xiàn)［3］提出采用多項(xiàng)式法表示晶振相位偏差，利用歷史數(shù)據(jù)計(jì)算出多項(xiàng)式參數(shù)進(jìn)行預(yù)測(cè)，該方法簡(jiǎn)單但短期內(nèi)精度不高；文獻(xiàn)［4］提出采用Kalman濾波算法進(jìn)行預(yù)測(cè)，但是僅考慮了老化因素對(duì)頻率漂移的影響，沒(méi)有考慮溫度變化的影響；文獻(xiàn)［5］提出通過(guò)數(shù)字濾波器分離出溫度、老化對(duì)頻率漂移的影響，而后采用Kalman濾波算法分別對(duì)溫度和老化的影響進(jìn)行預(yù)測(cè)后再相加，該方法參數(shù)選擇復(fù)雜，模型噪聲難確定，效果不太明顯；文獻(xiàn)［6］提出采用神經(jīng)網(wǎng)絡(luò)算法進(jìn)行建模，該方法收斂速度較慢，效率較低，不能很好地滿(mǎn)足實(shí)時(shí)性要求。

針對(duì)上述問(wèn)題，本文提出采用基于Wiener過(guò)程的可靠性建模方法，將短期內(nèi)的晶振退化問(wèn)題近似為線性問(wèn)題，建立晶振老化與頻率漂移之間的線性模型，并將溫度變化作為系統(tǒng)的加速應(yīng)力融合到模型參數(shù)中。首先利用極大似然估計(jì)和Kalman濾波算法分別對(duì)模型的擴(kuò)散參數(shù)和漂移參數(shù)進(jìn)行估計(jì)，然后實(shí)時(shí)預(yù)測(cè)晶振的頻率漂移量，對(duì)系統(tǒng)輸出頻率進(jìn)行補(bǔ)償，確保系統(tǒng)保持高精度頻率輸出。

1　問(wèn)題描述

GPS信號(hào)校準(zhǔn)晶振型頻率源的基本原理如圖1所示。當(dāng)GPS信號(hào)正常時(shí)，系統(tǒng)將GPS接收機(jī)輸出的秒信號(hào)與地面晶振分頻后的信號(hào)進(jìn)行比對(duì)，獲得相對(duì)頻差，再將頻差轉(zhuǎn)換成控制修正電壓，不斷調(diào)節(jié)晶振的振蕩頻率，最終實(shí)現(xiàn)晶振高準(zhǔn)確度頻率輸出。當(dāng)GPS信號(hào)丟失后，系統(tǒng)通過(guò)晶振頻率漂移預(yù)測(cè)模塊實(shí)時(shí)預(yù)測(cè)晶振輸出頻率的漂移量，再將其轉(zhuǎn)換成控制修正電壓對(duì)晶振進(jìn)行校正補(bǔ)償，實(shí)現(xiàn)相對(duì)高準(zhǔn)確度的頻率輸出［7］。

圖1　地面晶振馴服保持技術(shù)原理圖

其中，晶振頻率漂移預(yù)測(cè)模塊的工作原理如圖2所示。首先建立晶振基于溫度、老化特性的頻率漂移模型，然后利用壓控電壓值將GPS信號(hào)正常時(shí)的頻差信號(hào)量等效轉(zhuǎn)換成相應(yīng)的晶振頻率漂移量，利用對(duì)應(yīng)的溫度值和時(shí)刻值估計(jì)模型參數(shù)。當(dāng)GPS信號(hào)丟失后利用該模型，通過(guò)當(dāng)前時(shí)刻值和溫度值預(yù)測(cè)出相應(yīng)的晶振輸出頻率漂移量。

圖2　頻率漂移預(yù)測(cè)模塊原理圖

1.1晶振老化預(yù)測(cè)模型

一般情況下，晶振的老化模型是非線性的，頻率溫度變化模型是線性的。但老化率的非線性是對(duì)于較長(zhǎng)時(shí)間而言的，如一月或者一年。在較短時(shí)間內(nèi)，老化率可以認(rèn)為是線性的，如一天。由于GPS信號(hào)丟失的時(shí)間一般都很短，并且頻率源系統(tǒng)要求的精度保持時(shí)間基本都是一天。因此在保證精度的前提下，晶振的老化模型可以簡(jiǎn)化為線性，從而大大簡(jiǎn)化了算法處理。

基于Wiener過(guò)程的退化模型是一類(lèi)典型的描述設(shè)備隨機(jī)退化過(guò)程的線性隨機(jī)模型。到目前為止，這類(lèi)模型已經(jīng)在設(shè)備腐蝕、機(jī)械振動(dòng)、液晶顯示器劣化等退化中得到了廣泛應(yīng)用［8］。本文將短期內(nèi)的晶振退化過(guò)程描述為一元Wiener過(guò)程：

式中：X（t）——t時(shí)刻的退化量；

μ、σ——漂移參數(shù)和擴(kuò)散參數(shù)；

B（t）——標(biāo)準(zhǔn)布朗運(yùn)動(dòng)，且有B（t）～N（0，t）。

由于一元Wiener過(guò)程是齊次馬爾科夫過(guò)程，則其均值和方差分別為

可以看出X（t）的均值和方差隨時(shí)間線性增加。其變異系數(shù)為

1.2基于溫度應(yīng)力下的晶振老化預(yù)測(cè)模型

在某些Wiener過(guò)程的可靠性建模中，為了提高產(chǎn)品的性能退化過(guò)程，以便在短期內(nèi)獲得較多性能退化數(shù)據(jù)，通常需要提高試驗(yàn)應(yīng)力加快其退化過(guò)程，如改變溫度、壓力、濕度等。若設(shè)某應(yīng)力變化量為L(zhǎng)，在不同應(yīng)力下一元Wiener過(guò)程的參數(shù)是應(yīng)力的函數(shù)，即：

其中θμ是未知參數(shù)向量。

若考慮將溫度看作加速應(yīng)力作用于晶振的退化過(guò)程，由于頻率溫度變化模型可近似為線性［9］，那么基于溫度應(yīng)力下的退化模型為

其中a、b為模型未知參數(shù)。

由于本文所研究的GPS校準(zhǔn)晶振型頻率源系統(tǒng)使用的是恒溫晶振，其頻率漂移受外界溫度變化影響較小，在研究中可忽略不計(jì)。下面重點(diǎn)討論晶振的老化對(duì)頻率漂移的影響。

2　基于Wiener過(guò)程的GPS信號(hào)校準(zhǔn)晶振型頻率源守頻方法

2.1模型參數(shù)估計(jì)

2.1.1估計(jì)擴(kuò)散參數(shù)σ2

利用前文所述頻率源系統(tǒng)，設(shè)頻率源在測(cè)量時(shí)刻t1，t2，…，tn的頻率漂移量分別為x1，x2，…，xn，令0＜i＜n，其中t0=0，n≥1，x0=0。Δxi=xi-xi-1是時(shí)刻ti-1，ti之間的性能退化量。由Wiener過(guò)程的性質(zhì)［10-11］得

其中，Δti=ti-ti-1，i=1，2，…，n。由式（8）得到似然函數(shù)為

那么利用極大似然估計(jì)方法可求得σ2的估計(jì)值為

2.1.2估計(jì)漂移參數(shù)μ

對(duì)于參數(shù)μ的估計(jì)，若采用極大似然估計(jì)法有：

由式（11）可知，若采用極大似然估計(jì)方法估計(jì)參數(shù)μ時(shí)，使用的狀態(tài)值為當(dāng)前時(shí)刻值和初始時(shí)刻值的差值。如果當(dāng)前時(shí)刻值存在較大噪聲，就會(huì)對(duì)μ的估計(jì)值產(chǎn)生較大影響。因此，為了降低噪聲因素的影響，充分利用先前的歷史校準(zhǔn)數(shù)據(jù)，提高μ的估計(jì)準(zhǔn)確度，本文采用Kalman濾波算法對(duì)μ進(jìn)行預(yù)測(cè)估計(jì)。

Kalman濾波適合解決最優(yōu)線性濾波和估計(jì)問(wèn)題，它以最小均方誤差為準(zhǔn)則，根據(jù)前一次估計(jì)值和最近一次觀測(cè)數(shù)據(jù)估計(jì)信號(hào)的當(dāng)前值，適用于實(shí)時(shí)處理［12］。通過(guò)之前建立的Wiener過(guò)程模型推導(dǎo)，可得關(guān)于參數(shù)μ的狀態(tài)方程和觀測(cè)方程為

其中，yk=Δxk=xk-xk-1，Δtk=tk-tk-1，k=1，2，…，n；Wk為系統(tǒng)噪聲，具有零均值和已知協(xié)方差陣Qk；Vk為觀測(cè)噪聲，且服從N（0，σ2）。

Kalman濾波過(guò)程如下：

1）由公式Pk-=APk-1AT+Qk，計(jì)算向前推算誤差協(xié)方差Pk-。

2）根據(jù) Kk=Pk-HkT（HkPk-HkT+Rk）-1，計(jì)算 Kalman增益Kk。

5）由Pk=（I-KkHk）Pk-，更新誤差協(xié)方差Pk。

其中，A=1，Hk=Δtk，k-為 k時(shí)刻的先驗(yàn)狀態(tài)估計(jì)為k-1時(shí)刻的后驗(yàn)狀態(tài)估計(jì)，Pk-為k時(shí)刻的先驗(yàn)估計(jì)誤差的協(xié)方差。

實(shí)際中，需給定初值μ0估計(jì)誤差協(xié)方差矩陣初值P0。模型噪聲協(xié)方差陣Qk采用經(jīng)驗(yàn)估計(jì)，按照使濾波效果最優(yōu)的原則由實(shí)驗(yàn)確定。觀測(cè)噪聲協(xié)方差矩陣

2.2頻率漂移量預(yù)測(cè)

由于一元Wiener過(guò)程是齊次馬爾科夫過(guò)程，根據(jù)性質(zhì)可知其狀態(tài)增量是平穩(wěn)相互獨(dú)立的。當(dāng)x1，x2，…，xn已知時(shí)，設(shè)t＞n，對(duì)于t時(shí)刻的晶振輸出頻率漂移量X（t），有：

由上式可得晶振輸出頻率漂移X（t）的預(yù)測(cè)值為

2.3頻率輸出校正補(bǔ)償

預(yù)測(cè)出晶振的頻率漂移量后，將其轉(zhuǎn)換成壓控電壓，實(shí)時(shí)輸入晶振來(lái)調(diào)節(jié)輸出頻率漂移，從而完成晶振的自校正補(bǔ)償，使系統(tǒng)在GPS信號(hào)丟失的短期內(nèi)保持相對(duì)高精度的頻率輸出，如圖3所示。

圖3　頻漂校正補(bǔ)償原理圖

晶振的性能退化是連續(xù)的，頻漂預(yù)測(cè)校正補(bǔ)償是負(fù)反饋環(huán)節(jié)，具有一定的延遲。但是由于晶振在相鄰兩時(shí)刻的頻漂量之差，與頻漂預(yù)測(cè)值相差幾個(gè)數(shù)量級(jí)，可以直接將前一時(shí)刻的頻漂預(yù)測(cè)值作為當(dāng)前時(shí)刻的頻率漂移值用于校正補(bǔ)償。

綜上所述，解決GPS信號(hào)丟失后晶振輸出頻率漂移預(yù)測(cè)補(bǔ)償問(wèn)題的步驟如下：

1）分析老化、溫度等因素對(duì)晶振頻率漂移的影響關(guān)系，建立頻率漂移預(yù)測(cè)的一元Wiener過(guò)程模型。

2）根據(jù)頻率漂移預(yù)測(cè)模型推導(dǎo)出狀態(tài)似然函數(shù)和狀態(tài)方程，利用極大似然估計(jì)方法估計(jì)擴(kuò)散參數(shù)σ2，利用Kalman濾波算法估計(jì)漂移參數(shù)μ。

3）根據(jù)推導(dǎo)出的頻率漂移量預(yù)測(cè)公式，計(jì)算在GPS信號(hào)丟失情況下的晶振輸出頻率漂移量。

4）將預(yù)測(cè)出的晶振輸出頻率漂移量輸入D/A轉(zhuǎn)換模塊轉(zhuǎn)換為控制電壓，實(shí)時(shí)調(diào)節(jié)晶振輸出頻率，實(shí)現(xiàn)系統(tǒng)高精度頻率輸出。

3　實(shí)例驗(yàn)證

針對(duì)某GPS校準(zhǔn)晶振型頻率源系統(tǒng)，在GPS信號(hào)正常鎖定晶振的情況下，設(shè)置40s為采樣周期，采集頻差信號(hào)值及對(duì)應(yīng)的晶振壓控端輸入電壓。在具體測(cè)試驗(yàn)證中，使用的晶振為恒溫晶振，性能如表1所示［13］。其中采集的頻差信號(hào)值反映了系統(tǒng)的頻率輸出準(zhǔn)確度，如圖4所示。

表1　VCOCXO的老化率和頻率溫度系數(shù)

圖4　晶振鎖定后系統(tǒng)的頻率準(zhǔn)確度（采樣周期40s）

由圖4可以看出，GPS信號(hào)鎖定晶振信號(hào)后，系統(tǒng)輸出頻率值在晶振標(biāo)稱(chēng)頻率上下起伏，頻率準(zhǔn)確度最大起伏約為8.7×10-11。通過(guò)計(jì)算，圖中所顯示的頻率平均準(zhǔn)確度可達(dá)2.17×10-11，相對(duì)于表1所示晶振的老化率約為5×10-10/d有明顯改進(jìn)。由此可以說(shuō)明利用GPS信號(hào)校準(zhǔn)晶振后，系統(tǒng)的輸出頻率準(zhǔn)確度有明顯提高，GPS衛(wèi)星時(shí)鐘的長(zhǎng)期準(zhǔn)確度指標(biāo)在鎖定后的晶振上得到了復(fù)現(xiàn)，系統(tǒng)滿(mǎn)足大多數(shù)情況下的計(jì)量標(biāo)準(zhǔn)要求［14］。

根據(jù)控制電壓值將頻差信號(hào)進(jìn)行等效轉(zhuǎn)換，得到晶振頻率漂移量，如圖5所示。

由圖可以看出，晶振的頻率準(zhǔn)確度隨時(shí)間不斷降低，即頻率漂移量越來(lái)越大，整體上可近似為以老化率為斜率的線性退化。圖4與圖5比較可知，正是由于GPS信號(hào)對(duì)晶振信號(hào)進(jìn)行了校準(zhǔn)，才使得系統(tǒng)能夠保持較高的輸出精度。

當(dāng)晶振被GPS信號(hào)鎖定一段時(shí)間以后，斷開(kāi)GPS校準(zhǔn)信號(hào)，系統(tǒng)進(jìn)入晶振自校準(zhǔn)補(bǔ)償模式。圖6為使用保持方法后系統(tǒng)輸出頻率準(zhǔn)確度變化曲線。

圖6　使用保持方法后的頻率準(zhǔn)確度（采樣周期40s）

由圖可以看出，當(dāng)采用了保持方法之后，晶振的頻率漂移問(wèn)題得到了很好的修正，輸出頻率準(zhǔn)確度在GPS信號(hào)丟失后的11h內(nèi)仍保持在4.0×10-10之上，保持期間的頻率平均準(zhǔn)確度為9.7×10-11。但從圖中的曲線波動(dòng)趨勢(shì)發(fā)現(xiàn)，隨著時(shí)間的延長(zhǎng)，系統(tǒng)的頻率準(zhǔn)確度呈下降趨勢(shì)，說(shuō)明了該頻率保持方法可以使恒溫晶振在短期內(nèi)保持較高的輸出精度，但由于恒溫晶振的老化和系統(tǒng)計(jì)算時(shí)迭代誤差的影響，使得系統(tǒng)的輸出頻率準(zhǔn)確度不斷降低?？傮w來(lái)講，該方法能夠使系統(tǒng)在短期內(nèi)保持相對(duì)較高的頻率準(zhǔn)確度，滿(mǎn)足短期內(nèi)的計(jì)量需求。

4　結(jié)束語(yǔ)

本文研究了一種基于Wiener過(guò)程的頻率源守頻保持方法。首先建立了晶振性能退化的一元Wiener過(guò)程模型，利用極大似然估計(jì)和Kalman濾波算法估計(jì)出模型參數(shù)，計(jì)算得出晶振輸出頻率漂移量，當(dāng)GPS信號(hào)丟失后，實(shí)現(xiàn)晶振輸出頻率漂移的自校正補(bǔ)償。將該方法用于某GPS校準(zhǔn)晶振型頻率源系統(tǒng)，能夠滿(mǎn)足系統(tǒng)在GPS信號(hào)丟失的短期內(nèi)，仍保持高準(zhǔn)確度頻率輸出的要求。

［1］SANLEMD， DRISCOLLOC，LACHAPELLEG. Methodology for comparing two carrier phase tracking techniques［J］.GPS Solutions，2012，16（2）：197-207.

［2］HAN H Z，WANG J，WANG J，et al.Reliable partial ambiguity resolution for single-frequency GPS/BDS and INS integration［J］.GPS Solutions，2016，20（14）：1-14.

［3］DIAOYH，LIM，WANGXJ，etal.Satellite disciplined crystal oscillator system based on Kalman filter and PI algorithm［C］∥The 9th IEEE Conference on Industrial Electronics and Applications，2014：624-628.

［4］曾亮，孟慶杰，徐偉.利用GPS馴服校頻技術(shù)提高晶振性能［J］.計(jì)量技術(shù)，2008，42（5）：6-8.

［5］NICHOLLE C W T，CARLETON G C.Adaptive OCXO driftcorrectionalgorithm［C］∥IEEE International UFFC Joint 50th Anniversary Conference，2004：509-517.

［6］柳麗，陳之純，曾元峰.石英晶振老化的建模與模型求解［J］.上海航天，2004，21（3）：1-6.

［7］SONG J，AN Q，LIU S B.A high-resolution time to digital converter implemented in field programmable gate arrays［J］.IEEE Trans on Nuclear Science，2006，53（1）：236-241.

［8］司小勝，胡昌華，周東華.基于退化建模的備件需求預(yù)測(cè)［C］∥25th Chinese Control and Decision Conference，2013：4817-4820.

［9］HU C H，LEE M Y，TANG J.Optimum step-stress accelerateddegradationtestforWienerdegradation processunderconstraints［J］.EuropeanJournalof Operational Research，2015，241（2）：412-421.

［10］彭寶華.基于Wiener過(guò)程的可靠性建模方法研究［D］.長(zhǎng)沙：國(guó)防科學(xué)技術(shù)大學(xué)，2010.

［11］陳逸成，何華鋒，鄧會(huì)選，等.基于Kalman濾波和Wiener過(guò)程的系統(tǒng)性能退化數(shù)據(jù)建模方法［J］.電光與控制，2015，22（7）：107-110.

［12］王宇譜，呂志平，陳正生，等.一種新的鐘差預(yù)處理方法及在WNN鐘差中長(zhǎng)期預(yù)報(bào)中的應(yīng)用［J］.武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2016，41（3）：373-379.

［13］倪媛媛，胡永輝，何在民.北斗衛(wèi)星校準(zhǔn)銣鐘單元的設(shè)計(jì)與實(shí)現(xiàn)［J］.測(cè)試技術(shù)學(xué)報(bào)，2011，15（2）：56-59.

［14］蔣陸萍，曾祥君，李澤文，等.基于GPS實(shí)現(xiàn)電力系統(tǒng)高準(zhǔn)確度同步時(shí)鐘［J］.電網(wǎng)技術(shù)，2011，35（2）：201-206.

（編輯：徐柳）

Timekeeping method for frequency sources calibrated by GPS based on Wiener progress

YANG Shaochen，HU Changhua，LI Hongzeng
（Dept.of Control Engineering，the Second Artillery Engineering University，Xi’an 710025，China）

The GPS calibrating crystal oscillator frequency source system can only rely on the crystal oscillator to keep the output at a high standard after GPS signals are lost.In this paper，an algorithm of maintaining frequency source punctuality based on Wiener process is proposed to correct and compensate for the output frequency drift of the crystal oscillator.Through analysis of historical calibrating data，an output frequency drift model was established with temperature and agingcharacteristics.Maximumlikelihoodestimationwasappliedtoestimatethediffusion parameters of the model and Kalman filtering method was used to estimate its drift parameters. Furthermore，the output frequency drift in real time was predicted and the frequency of the system was compensated.The proposed method was applied to a GPS calibrating crystal oscillator frequency source system，and the results show that it can maintain the output precision of frequency at a degree as high as 9.7×10-11after GPS signals are lost.

GPS calibrating frequency；crystal oscillator；frequency drift；Wiener process

A

1674-5124（2016）06-0014-05

10.11857/j.issn.1674-5124.2016.06.004

2016-01-22；

2016-03-07

國(guó)家杰出青年基金（61025014）國(guó)家自然科學(xué)基金（61573365）

楊少塵（1991-），男，北京市人，碩士研究生，專(zhuān)業(yè)方向?yàn)楣收显\斷與預(yù)測(cè)。