于 帆，陳 偉

(西安工業(yè)大學 計算機科學與工程學院，陜西 西安 710021)

GPS/BD雙模融合的高精度時間同步方法研究

于 帆，陳 偉

(西安工業(yè)大學 計算機科學與工程學院，陜西 西安 710021)

為滿足電力和測控通訊等領(lǐng)域?qū)Ω呔葧r鐘的要求，結(jié)合GPS/BD時鐘無累積誤差和晶振時鐘無隨機誤差的特性，提出了一種新型的基于GPS/BD雙模接收機的時鐘馴服的軟、硬件設(shè)計及其實現(xiàn)方法。該方法通過構(gòu)建權(quán)重系數(shù)將GPS和BD信號濾波得到的鐘差估計值和鐘差速度估計值進行融合，實現(xiàn)了基于GPS/BD雙模融合授時的時鐘馴服系統(tǒng)。在系統(tǒng)軟件設(shè)計中，采用滑動中位數(shù)方法剔除測量數(shù)據(jù)中的異常點，應(yīng)用無偏FIR濾波算法對鐘差數(shù)據(jù)進行濾波，以提高時鐘偏差的濾波精度；在ARM的控制下，應(yīng)用數(shù)字PID控制算法產(chǎn)生控制電壓，由D/A轉(zhuǎn)換器調(diào)節(jié)晶振輸出頻率，實現(xiàn)快速馴服本地時鐘。實驗結(jié)果表明，所提出時鐘馴服方法的時間同步精度達到了50 ns，優(yōu)于基于單模GPS和單模北斗的時間同步系統(tǒng)。

GPS/BD雙模；時間同步；鎖相環(huán)技術(shù)；無偏FIR濾波

0 引 言

預計2020年全球?qū)⑦M入5G時代，在5G時代人們可以享受千倍提速的網(wǎng)絡(luò)、通信等服務(wù)，這些便利的服務(wù)要求時鐘系統(tǒng)具有極高的準確性和穩(wěn)定性，對時鐘精度要求甚至達到納秒級別，并且各個系統(tǒng)都要求達到嚴格的時間同步。

衛(wèi)星授時是目前主流的時間同步技術(shù)，其中美國的全球衛(wèi)星導航系統(tǒng)(GPS)技術(shù)最為成熟，憑借覆蓋面廣、精度高等特點成為了衛(wèi)星授時的首選[1]。但是GPS歸美國政府所有，由美國軍方開發(fā)和控制，存在著故意降低精度的可能，甚至在戰(zhàn)爭等不確定因素下可能導致中國等其他地區(qū)不能使用GPS服務(wù)，對國內(nèi)各種GPS應(yīng)用造成了潛在隱患。

北斗衛(wèi)星導航系統(tǒng)(BD)是中國自主建設(shè)、獨立運行的衛(wèi)星導航系統(tǒng)，全面建立完成后將為全球用戶提供全天候、全天時、高精度的定位、導航和授時服務(wù)，精度也將更高?？梢酝耆毩⒐ぷ鳎瑸閲鴥?nèi)電力、導航、勘測等領(lǐng)域進行服務(wù)[2]。

現(xiàn)有的GPS/BD雙模授時技術(shù)是將GPS授時信號作為優(yōu)先選擇，BD授時信號作為備份，GPS信號不能正常工作時再使用BD授時信號，這種基于單模接收機的時間同步方案[3-5]，精度通常較低，系統(tǒng)的性能指標還有進一步提升的空間。為此，提出了一種新型的基于GPS/BD雙模融合授時的時間同步方法，通過數(shù)據(jù)融合方法將GPS鐘差與BD鐘差相融合，提高了鐘差的估計精度，應(yīng)用無偏FIR濾波算法對鐘差數(shù)據(jù)進行濾波，提高了鐘差的濾波精度。

1 系統(tǒng)總體設(shè)計

1.1 時間同步方法原理

衛(wèi)星授時信號具備較高的長期穩(wěn)定性和高準確度等特點，利用衛(wèi)星授時技術(shù)產(chǎn)生一個穩(wěn)定的秒脈沖信號(1 PPS)周期性地對本地時鐘源進行馴服，使本地時鐘源輸出高穩(wěn)定的頻率信號且濾波后的誤差明顯減小[6]，系統(tǒng)可以獲得一個短期準確度和長期穩(wěn)定性都較高的時間頻率標準[7-9]。所使用恒溫晶體震蕩器作為本地時鐘源。

系統(tǒng)的基本調(diào)控原理基于鎖相環(huán)的調(diào)節(jié)機制[10-11]，使用PID控制算法調(diào)節(jié)恒溫晶振的輸入信號，使輸出信號和輸入信號保持一個固定關(guān)系。系統(tǒng)由衛(wèi)星接收模塊、時間間隔測量模塊、數(shù)據(jù)處理模塊和恒溫晶振等構(gòu)成，其中時間間隔測量模塊在FPGA中編程實現(xiàn)，數(shù)據(jù)處理模塊由ARM編程實現(xiàn)。時間同步方法原理圖如圖1所示。

圖1 系統(tǒng)原理圖

衛(wèi)星接收模塊輸出1PPS作為FPGA中時間間隔測量模塊的輸入，恒溫晶振經(jīng)分頻后產(chǎn)生頻率為1 Hz的周期信號作為時間間隔測量模塊的另一輸入，時間間隔測量模塊的功能包括脈沖計數(shù)、分頻和測頻等，測量出恒溫晶振分頻1 Hz信號與GPS的1PPS之間的鐘差；數(shù)據(jù)處理模塊使用無偏FIR濾波算法對鐘差數(shù)據(jù)進行處理，可以提高鐘差的濾波精度，PID控制算法計算出恒溫晶振相對于GPS的頻率偏差，將偏差值送入到D/A轉(zhuǎn)換器轉(zhuǎn)換成模擬電壓輸入到恒溫晶振的壓控輸入端，調(diào)整恒溫晶振的輸出頻率。

1.2 雙模融合

為進一步提高授時精度，提出了一種新型的基于GPS/BD雙模融合授時的時間同步方法。通過構(gòu)建權(quán)重系數(shù)將GPS和BD信號濾波得到的鐘差估計值和鐘差速度估計值進行融合，實現(xiàn)基于GPS/BD雙模融合授時的時鐘馴服系統(tǒng)。使用無偏FIR濾波算法對鐘差數(shù)據(jù)進行濾波，消除了衛(wèi)星授時信號不穩(wěn)定以及測量引起的誤差，提高了鐘差的濾波精度，獲得了高精度的鐘差估計值和鐘差速度估計值。融合計算公式如下：

(1)

其中，x為融合后的鐘差估計值；xGPS為基于GPS的鐘差估計值；xBD為基于BD的鐘差估計值。

(2)

GPS信號和BD信號都含有不同的噪聲成分，通過提出方法可以相互抵消一部分噪聲，提高精度。

2 軟件設(shè)計及算法

2.1 軟件設(shè)計

軟件部分包括FPGA軟件和ARM軟件，時間間隔測量在FPGA中進行，精確測量出1 PPS和恒溫晶振分頻產(chǎn)生的1 Hz信號上升沿的時間間隔。FPGA中采用Verilog硬件描述語言并結(jié)合自身的IP核相對較容易實現(xiàn)脈沖計數(shù)、分頻和測頻等功能。在此不考慮FPGA軟件程序設(shè)計。

數(shù)據(jù)處理模塊由ARM編程實現(xiàn)，執(zhí)行數(shù)據(jù)處理程序計算、分析和判斷時間間隔測量模塊送入的鐘差數(shù)據(jù)，并計算出頻率偏差值，根據(jù)頻率偏差值通過D/A轉(zhuǎn)換器對恒溫晶振進行頻率修正。系統(tǒng)的工作過程可以抽象為4個狀態(tài)，包括初始化狀態(tài)、跟蹤狀態(tài)、穩(wěn)定狀態(tài)和守時狀態(tài)，主要是根據(jù)接收機接收信號狀態(tài)和濾波得到的鐘差估計值大小來區(qū)分。系統(tǒng)狀態(tài)切換圖如圖2所示。

圖2 系統(tǒng)狀態(tài)切換圖

初始化狀態(tài)：加電后初始化整個系統(tǒng)，包括ARM初始化、下載配置FPGA等。

跟蹤狀態(tài)：快速調(diào)節(jié)本地晶振輸出頻率達到與參考頻率高度同步的狀態(tài)。

穩(wěn)定狀態(tài)：區(qū)別于跟蹤狀態(tài)每個控制周期對恒溫晶振進行調(diào)節(jié)，穩(wěn)定狀態(tài)下當且僅當鐘差值超過設(shè)定閾值時才進行調(diào)節(jié)。

守時狀態(tài)：衛(wèi)星信號不能正常工作時系統(tǒng)進入守時狀態(tài)，通過鐘差預測模型對恒溫晶振進行校正，保持恒溫晶振秒脈沖和頻率輸出的精度。

2.2 算法設(shè)計

無偏FIR濾波算法：有限長單位沖激響應(yīng)濾波器，又稱為非遞歸型濾波器，是數(shù)字信號處理系統(tǒng)中最基本的元件，它可以在保證任意幅頻特性的同時具有嚴格的線性相頻特性，同時其單位抽樣響應(yīng)是有限長的，因而濾波器是穩(wěn)定的系統(tǒng)。

時鐘馴服系統(tǒng)中實現(xiàn)對恒溫晶振的精密控制，關(guān)鍵一部分在于使用濾波算法提高鐘差的濾波精度。無偏FIR濾波算法在接收機鐘差估計的應(yīng)用，首先由國外Yuriy S. Shmaliy教授提出，在其文獻中說明了相應(yīng)的公式推導過程，并通過實驗比對得出在接收機鐘差估計中無偏FIR濾波比卡爾曼濾波效果更好的結(jié)論[12-13]，與Kalman濾波相比，其Allan方差要稍微好些。在軟件實現(xiàn)中將無偏FIR濾波算法改進成迭代形式，不需要存儲大量數(shù)據(jù)，便于算法移植。

無偏FIR濾波的思想是，通過采集一定數(shù)量的觀測值再利用卷積公式來計算要求的狀態(tài)估計值，計算公式如下：

(3)

其中，λ(n)為鐘差觀測值；Nk為采樣周期長度；W(i)為卷積系數(shù)。

(4)

將無偏FIR濾波算法應(yīng)用到系統(tǒng)中，采用的基本思想是通過低階的狀態(tài)估計量計算高階觀測量，從而利用無偏FIR濾波得到控制算法需要的鐘差估計值和鐘差速度估計值。時間間隔測量模塊送入數(shù)據(jù)處理模塊的是鐘差觀測量，鐘差估計值的計算方法如下：

(5)

將一階的觀測量λ1(i)在采樣周期N1中與卷積系數(shù)ξ1(i)進行卷積計算，可以得到一階的狀態(tài)估計值ω1(n)，即鐘差估計值。通過一階狀態(tài)估計值ω1(n)計算得到二階觀測值λ2(i)，計算公式如下：

(6)

計算得到N0個二階觀測值λ2(i)后，與卷積系數(shù)ξ0(i)進行卷積計算，得到二階狀態(tài)估計值ω2(n)，即系統(tǒng)中的鐘差速度，計算公式如下：

(7)

系統(tǒng)中使用數(shù)字PID控制算法計算電壓控制字，數(shù)字PID控制算法則是根據(jù)PID控制原理應(yīng)用在計算機或硬件控制中的程序算法，該算法具有穩(wěn)定、快速和準確的特性[14-15]。在設(shè)計中，為了簡化編程，減少了對存儲單元的占用。通過無偏FIR濾波得到控制算法需要的鐘差量和鐘差速度量。系統(tǒng)中的時鐘模型如下：

(8)

因為恒溫晶振具有較好的短期穩(wěn)定性，假設(shè)一個控制周期內(nèi)鐘差速度保持不變，即鐘差是速度的累積，則下一個控制節(jié)點的鐘差為：

(9)

其中，n為控制周期。

在一個周期n后需要將鐘差消除為零，因此需要施加的控制量u(i)為：

(10)

其中，kd為微分系數(shù)；kp為比例系數(shù)。

鐘差調(diào)節(jié)示意圖如圖3所示。

圖3 控制量計算示意圖

將控制量u(i)輸入D/A轉(zhuǎn)換器轉(zhuǎn)換成模擬電壓輸入恒溫晶振壓控輸入端，調(diào)整恒溫晶振輸出穩(wěn)定。

3 實驗結(jié)果

對所提出的一種新型的基于GPS/BD雙模融合授時的時間同步方法進行驗證，分別使用單模BD授時信號、單模GPS授時信號和GPS/BD雙模融合授時對恒溫晶振進行馴服，計算馴服后的恒溫晶振分頻產(chǎn)生的秒脈沖信號與GPS的秒脈沖信號之間的差值，采樣周期150s，采集100組數(shù)據(jù)，三種模式下的秒脈沖精度如圖4所示。

圖4 秒脈沖精度

由圖4可知，在系統(tǒng)中衛(wèi)星授時信號正常工作情況下，系統(tǒng)進入鎖定狀態(tài)后恒溫晶振分頻產(chǎn)生的秒脈沖信號與GPS的秒脈沖信號相位跟蹤誤差保持在一定范圍內(nèi)，基于單模接收機的時間同步方案，其中基于單模GPS的時間同步精度優(yōu)于基于單模BD時間同步系統(tǒng)，但是其時間同步精度仍然大于60ns，并且存在較大波動；基于GPS/BD雙模融合授時的時間同步方法提高了鐘差的估計精度，時間同步精度始終保持在50ns以內(nèi)。所提出的新型基于GPS/BD雙模接收機的時鐘馴服方法的秒脈沖精度明顯優(yōu)于基于單模GPS或單模BD的授時方案。

4 結(jié)束語

針對電力和測控通訊等領(lǐng)域?qū)Ω呔葧r鐘的要求，提出了一種新型的基于GPS/BD雙模接收機的時鐘馴服方法?；陔p模融合的設(shè)計原理、軟件設(shè)計流程和算法設(shè)計等要點，設(shè)計實現(xiàn)了ARM和FPGA相結(jié)合的時鐘馴服系統(tǒng)。實驗結(jié)果表明，采用數(shù)據(jù)融合方法提高了鐘差的估計精度，時間同步精度也進一步得到了提高，分別優(yōu)于基于單模GPS和單模北斗的時間同步系統(tǒng)，且能夠輸出高精度的頻率基準和時間同步信號，可為各領(lǐng)域提供精度更高、更準確的時間頻率信號。

[1] 單慶曉，楊 俊.衛(wèi)星馴服時鐘系統(tǒng)的新進展[J].測試技術(shù)學報，2009,23(5)：396-401.

[2] 郭 彬.基于北斗/GPS雙模授時的電力系統(tǒng)時間同步技術(shù)研究[D].長沙：湖南大學，2010.

[3] 周大勇，劉 鵬，劉樹昌，等.BD2/GPS高精度同步時鐘裝置的設(shè)計與應(yīng)用[J].吉林大學學報:信息科學版，2014,32(3)：262-266.

[4] 黨曉園，單曉慶，肖昌炎，等.基于GPS與北斗雙模授時的壓控晶振校頻系統(tǒng)的研究與設(shè)計[J].計算機測量與控制，2009,17(11)：2246-2248.

[5] 郭 彬，單慶曉，肖昌炎,等.電網(wǎng)時鐘系統(tǒng)的北斗/GPS雙模同步技術(shù)研究[J].計算機測量與控制，2011，19(1)：139-141.

[6]MosaviMR.UseofaccurateGPStimingbasedonradialbasisprobabilisticneuralnetworkinelectricsystems[C]//Internationalconferenceonelectricalandcontrolengineering.[s.l.]:[s.n.],2010：2572-2575.

[7] 張國琴，吳玉蓉.基于GPS校準晶振的高精度時鐘的設(shè)計[J].儀表技術(shù)，2010(4)：23-24.

[8] 謝 強，錢光弟.基于授時GPS的高精度頻率源設(shè)計與實現(xiàn)[J].工業(yè)控制計算機，2007,20(3):15-16.

[9] 于躍海，張道農(nóng)，胡永輝，等.電力系統(tǒng)時間同步方案[J].電力系統(tǒng)自動化，2008,32(7)：82-86.

[10]KimB.Dual-loopDPLLgear-shiftingalgorithmforfastsynchronization[J].IEEETransactionsonCircuitsandSystems-II:AnalogandDigitalSignalProcessing,1997,44(7):577-586.

[11] 侯興勃,崔保健,趙海鷹.全數(shù)字低頻鎖相環(huán)的研究與應(yīng)用[J].計算機測量與控制，2010，18(9)：2181-2182.

[12]ShmaliyYS.AnunbiasedFIRfilterforTIEmodelofalocalclockinapplicationstoGPS-basedtimekeeping[J].IEEETransactionsonUltrasonics,Ferroelectrics,andFrequencyControl,2006,53(5):862-870.

[13]ShmaliySY,Munoz-DiazJ,Arceo-MiquelL,etal.OptimumtimestepandmemoryforGPS-basedunbiasedFIRestimatesofthelocalclockTIEmodel[C]//IEEEinternationalfrequencycontrolsymposiumjointlywiththe21stEuropeanfrequencyandtimeforum.[s.l.]:IEEE,2007:291-296.

[14] 楊劍青，楊曉琴，謝 亮，等.基于ARM與CPLD的可馴頻率基準的設(shè)計與實現(xiàn)[J].電子測量技術(shù),2014,37(2)：75-77.

[15] 楊前利.基于數(shù)字PID的閉環(huán)溫度控制系統(tǒng)的設(shè)計[J].計算機與數(shù)字工程,2013,41(12):2004-2008.

Research on High-precision Time Synchronization Method with GPS/BD Dual Mode Fusion

YU Fan，CHEN Wei

(School of Computer Science and Engineering,Xi’an Technological University,Xi’an 710021,China)

In order to meet the request of high precision clock for the electric power and control communication fields,a new method of software and hardware design and implementation for clock disciplining with GPS/BD dual-mode receiver has been proposed combined with the characteristics of the GPS/BD clock without random error and oscillator without accumulation error,which could improve the estimation accuracy of clock bias by fusing the GPS clock bias data and the BD clock bias data.Furthermore,the sliding median method has been employed in the software design to eliminate the abnormal data in the measurement data and the accuracy of the clock filter has been improved with unbiased FIR filtering algorithm to filter the clock bias data.Under the control of ARM processor,the control voltage for oscillator has been acquired by adopting the digital PID control algorithm and the oscillator output frequency has been adjusted with D/A converter,which aimed at achieving a fast tamable local oscillator.According to the experimental results,the time synchronization precision of the proposed method has reached 50 ns,which is higher than single mode GPS and single mode BD time synchronization system.

GPS/BD dual-mode;time synchronization;phase-locked loop technology;unbiased FIR filter

2016-08-09

2016-11-09 網(wǎng)絡(luò)出版時間：2017-03-13

國家自然科學基金資助項目(60904089)

于 帆(1963-)，男，教授，研究方向為計算機控制技術(shù)與人工智能；陳 偉(1991-)，男，碩士研究生，研究方向為嵌入式控制。

http://kns.cnki.net/kcms/detail/61.1450.tp.20170313.1547.084.html

TP31

A

1673-629X(2017)05-0201-04

10.3969/j.issn.1673-629X.2017.05.042