陸加海,張同雙,陶小紅,陳浩謙,李春海

(1.上海交通大學(xué),上海200030;2.中國(guó)衛(wèi)星海上測(cè)控部,江蘇江陰 214431)

1 引 言

測(cè)控系統(tǒng)要求有高精度的時(shí)間頻率系統(tǒng)。為了保持站點(diǎn)間時(shí)間同步,一般采用銣原子鐘或恒溫晶振(Oven-Controlled Crystal Oscillator,OCXO)來(lái)提供時(shí)間頻率信號(hào)。前者的成本較高,后者頻率準(zhǔn)確度較差,不能滿足測(cè)控系統(tǒng)的指標(biāo)要求。

為解決上述問(wèn)題,本文采用GPS對(duì)恒溫晶振進(jìn)行馴服,綜合了GPS時(shí)間信號(hào)長(zhǎng)穩(wěn)指標(biāo)高的優(yōu)點(diǎn)和恒溫晶振短穩(wěn)指標(biāo)高的優(yōu)點(diǎn),從而獲得了高穩(wěn)定度和高準(zhǔn)確度的頻率信號(hào)。

為了計(jì)算恒溫晶振相對(duì)于GPS的頻率偏差,需要長(zhǎng)時(shí)間測(cè)量GPS秒信號(hào)與恒溫晶振分頻信號(hào)之間的時(shí)差,通過(guò)時(shí)差計(jì)算得到頻率偏差,再根據(jù)頻率偏差來(lái)校準(zhǔn)恒溫晶振,所以時(shí)差測(cè)量的精度直接影響到頻率馴服的精度。文獻(xiàn)[1]采用FGPA延遲線內(nèi)插法測(cè)時(shí)差,測(cè)量精度達(dá)到121 ps,但測(cè)量線性化較差、成本高。文獻(xiàn)[2]介紹了模擬內(nèi)插法測(cè)時(shí)差,測(cè)量精度達(dá)到400 ps,但易受噪聲干擾、實(shí)現(xiàn)困難;

本文采用時(shí)間數(shù)字轉(zhuǎn)換器(Time to Digital Converter,TDC)測(cè)量時(shí)差,精度達(dá)到250 ps,極大減小了測(cè)量誤差,且設(shè)計(jì)簡(jiǎn)單、成本低,有效實(shí)現(xiàn)了GPS馴服恒溫晶振系統(tǒng)的設(shè)計(jì)。

2 系統(tǒng)設(shè)計(jì)

2.1 系統(tǒng)硬件組成

系統(tǒng)組成框圖如圖1所示,主要由GPS模塊、TDC時(shí)間數(shù)字轉(zhuǎn)換器、分頻模塊、數(shù)據(jù)處理單元、D/A轉(zhuǎn)換模塊和恒溫晶振等組成。GPS模塊用于輸出標(biāo)準(zhǔn)的GPS秒脈沖信號(hào),分頻模塊將恒溫晶振頻率信號(hào)分頻成1 PPS信號(hào)送TDC模塊,TDC模塊用于測(cè)量恒溫晶振分頻1 PPS信號(hào)與GPS秒信號(hào)之間的時(shí)差,數(shù)據(jù)處理單元對(duì)時(shí)差數(shù)據(jù)進(jìn)行處理后得到恒溫晶振相對(duì)于GPS的頻率偏差,并控制D/A轉(zhuǎn)換模塊對(duì)恒溫晶振進(jìn)行頻率修正。

圖1 系統(tǒng)組成框圖Fig.1 Block diagram of the system

2.2 時(shí)差測(cè)量電路設(shè)計(jì)

TDC是時(shí)間數(shù)字轉(zhuǎn)換器的英文縮寫。這類芯片把時(shí)差直接轉(zhuǎn)化為高精度的數(shù)字值。它們被描述為與ADC(模擬數(shù)字轉(zhuǎn)換器)功能相似的器件。通常TDC是用在分辨率小于1 ns的測(cè)量場(chǎng)合。本設(shè)計(jì)使用TDC-GP1芯片進(jìn)行測(cè)量,其測(cè)量電路框圖如圖2所示。將待測(cè)的GPS秒信號(hào)和晶振1 PPS信號(hào)送TDC模塊;START信號(hào)由分頻模塊產(chǎn)生并送TDC模塊,用于控制測(cè)量時(shí)序。TDC模塊將時(shí)差轉(zhuǎn)化為數(shù)字值送數(shù)據(jù)處理單元。

圖2 時(shí)差測(cè)量電路框圖Fig.2 Block diagram of time difference measurement circuit

時(shí)差在TDC模塊的測(cè)量時(shí)序如圖 3所示,將GPS秒信號(hào)與恒溫晶振分頻得到的1 PPS信號(hào)進(jìn)行比對(duì),得到的時(shí)間間隔即是待測(cè)的時(shí)差。

圖3 時(shí)差測(cè)量時(shí)序圖Fig.3 Time difference measurement timing diagram

圖3中START信號(hào)為測(cè)試開(kāi)始信號(hào),GPS秒信號(hào)與晶振1 PPS信號(hào)之間的時(shí)間間隔T為待測(cè)的時(shí)差。信號(hào)脈沖寬度要求大于2.5 ns,開(kāi)始信號(hào)與被測(cè)信號(hào)之間時(shí)間間隔要求大于3 ns、小于7.6 μ s。

為了保證每次被測(cè)信號(hào)落在允許的測(cè)量范圍內(nèi)。測(cè)量開(kāi)始之前應(yīng)進(jìn)行初步的時(shí)間同步,確保GPS秒與晶振1 PPS之間的時(shí)差小于3 μ s。START信號(hào)通過(guò)將GPS秒延時(shí)約999 996 μ s后產(chǎn)生,確保START信號(hào)與GPS秒信號(hào)之間的時(shí)差保持在4 μ s左右。

2.3 分頻模塊的設(shè)計(jì)

分頻模塊采用CPLD設(shè)計(jì)實(shí)現(xiàn),通過(guò)CPLD分頻計(jì)數(shù)將恒溫晶振的10 MHz頻率信號(hào)分頻為1 PPS信號(hào),同時(shí)將 GPS秒延時(shí)999 996 μ s后生成 START信號(hào),確保START信號(hào)超前GPS秒信號(hào)4 μ s。當(dāng)晶振1 PPS與GPS秒之間的時(shí)差大于3 μ s時(shí)通過(guò)計(jì)數(shù)器清零的方法進(jìn)行同步,確保兩者之間的時(shí)差小于3 μ s。其時(shí)序如圖4所示。

圖4 分頻模塊時(shí)序圖Fig.4 Timing diagram of frequency division module

3 數(shù)據(jù)處理

根據(jù)頻率和周期互為倒數(shù)的關(guān)系,可采用測(cè)時(shí)差的方法[4](以GPS的秒信號(hào)為參考)來(lái)測(cè)量本地頻率源的頻率準(zhǔn)確度(Δf/f),作為校頻的依據(jù)。計(jì)算公式如下:

按照公式(1)可以計(jì)算出恒溫晶振相對(duì)于GPS的頻率偏差[5],其中 Δt2、Δt1分別為 t2、t1時(shí)刻測(cè)得的晶振1 PPS與GPS秒之間的時(shí)差值。根據(jù)頻率偏差的大小,將需要調(diào)整的頻率偏差量轉(zhuǎn)換為模擬電壓施加在壓控恒溫晶振的電壓控制端,便可以實(shí)現(xiàn)對(duì)恒溫晶振進(jìn)行頻率修正。

但GPS信號(hào)在傳輸過(guò)程中存在各種干擾[3],這些干擾造成GPS輸出的秒信號(hào)存在一定的抖動(dòng),其短穩(wěn)較差。通過(guò)對(duì)時(shí)差值進(jìn)行遞推平均濾波可以有效消除GPS秒信號(hào)抖動(dòng)的影響[6]。取測(cè)量得到的時(shí)差數(shù)據(jù) Δt1,Δt2,Δt3,…,Δtn,…,以 n 作為窗口長(zhǎng)度對(duì)其中的n個(gè)測(cè)量值取平均,通過(guò)移動(dòng)窗口得到新的時(shí)差數(shù)據(jù)序列[1]ΔT1,ΔT2,ΔT3,…,ΔTn,…,如公式(2)所示:

為了使頻率準(zhǔn)確度 Δf/f測(cè)量更準(zhǔn)確,對(duì)遞推平均濾波算法進(jìn)行優(yōu)化,取 ΔTj,ΔTn+j,ΔT2n+j,ΔT3n+j,…,ΔTmn+j,…,作為實(shí)際測(cè)量序列進(jìn)行計(jì)算,如公式(3)所示:

將式(2)代入式(3)得到

圖5是使用本系統(tǒng)測(cè)量得到的恒溫晶振相對(duì)于GPS系統(tǒng)的時(shí)差曲線,該曲線由于GPS秒信號(hào)抖動(dòng)而存在一定的波動(dòng),使用這些數(shù)據(jù)計(jì)算得到的頻率偏差也會(huì)受到影響而發(fā)生波動(dòng),所以不能直接使用。圖6是添加遞推平均濾波處理后得到的時(shí)差數(shù)據(jù),其濾波窗口n設(shè)置為120。濾波后時(shí)差數(shù)據(jù)的抖動(dòng)減小約兩個(gè)數(shù)量級(jí)。

圖5 時(shí)差數(shù)據(jù)濾波前曲線Fig.5 Curve before time difference data filtering

圖6 時(shí)差數(shù)據(jù)濾波后曲線Fig.6 Curve after time difference data filtering

頻率調(diào)節(jié)通過(guò)16位DAC(數(shù)字模擬轉(zhuǎn)換器)實(shí)現(xiàn),將需要調(diào)整的頻率偏差值轉(zhuǎn)換為模擬電壓,施加在恒溫晶振的控制電壓輸入端,實(shí)現(xiàn)對(duì)恒溫晶振頻率的調(diào)控。

4 試 驗(yàn)

本設(shè)計(jì)中采用MV180 OCXO雙層恒溫晶振,其頻率牽引范圍為±3×10-7,外部控制電壓范圍為0～+5 V。GPS模塊采用Motolora M12+,在位置保持模式下,輸出秒信號(hào)的抖動(dòng)方差 σ<20 ns。時(shí)間數(shù)字轉(zhuǎn)換器使用TDC-GP1芯片,其測(cè)量精度達(dá)到250 ps,測(cè)量誤差遠(yuǎn)小于GPS秒信號(hào)抖動(dòng)引入的誤差。16位數(shù)模轉(zhuǎn)換芯片采用DAC8571,其輸出電壓對(duì)恒溫晶振的理論頻率調(diào)諧靈敏度優(yōu)于1×10-11。圖7所示為根據(jù)濾波后的時(shí)差數(shù)據(jù)計(jì)算得到的頻率準(zhǔn)確度曲線。不采用濾波算法頻率馴服時(shí)間為數(shù)個(gè)小時(shí),采用濾波算法后頻率馴服時(shí)間縮短為約1 200 s。實(shí)驗(yàn)表明,馴服后恒溫晶振的頻率準(zhǔn)確度優(yōu)于3×10-11。

圖7 馴服后頻率準(zhǔn)確度曲線Fig.7 Tamed frequency accuracy curve

5 小 結(jié)

本文通過(guò)采用時(shí)間數(shù)字轉(zhuǎn)換器測(cè)量時(shí)差,簡(jiǎn)化了GPS馴服恒溫晶振系統(tǒng)的硬件設(shè)計(jì),測(cè)量精度達(dá)到250 ps,極大減小了測(cè)量誤差。設(shè)計(jì)了優(yōu)化的遞推平均濾波算法對(duì)時(shí)差數(shù)據(jù)進(jìn)行濾波處理,消除了GPS秒信號(hào)抖動(dòng)引入的干擾,縮短了頻率馴服時(shí)間。頻率馴服時(shí)間約為1 200 s,與國(guó)內(nèi)同類應(yīng)用相比,縮短約一個(gè)數(shù)量級(jí);頻率馴服精度優(yōu)于3×10-11。

[1]盧祥弘,陳儒軍,何展翔.基于FPGA的恒溫晶振頻率校準(zhǔn)系統(tǒng)的設(shè)計(jì)[J].電子技術(shù)應(yīng)用,2010(7):104-107.LU Xiang-hong,CHEN Ru-jun,HE Zhan-xiang.The designof frequency calibrationsystem for OCXO based onFPGA[J].Measurement Control Technology and Instruments,2010(7):104-107.(in Chinese)

[2]孫杰,潘繼飛.高精度時(shí)間間隔測(cè)量方法綜述[J].計(jì)算機(jī)測(cè)量與控制,2007,15(2):145-148.SUN Jie,PAN Ji-fei.Methods of High Precision Time-Interval Measurement[J].Computer Measurement&Control,2007,15(2):145-148.(in Chinese)

[3]單慶曉,楊俊.衛(wèi)星馴服時(shí)鐘系統(tǒng)的新進(jìn)展[J].測(cè)試技術(shù)學(xué)報(bào),2009,23(5):396-401.SHAN Qing-xiao,YANG Jun.Review of Satellite Disciplined Clock System[J].Journal of Test and Measurement Technology,2009,23(5):396-401.(in Chinese)

[4]楊旭海,翟惠生,胡永輝,等.基于新校頻算法的GPS可馴銣鐘系統(tǒng)研究[J].儀器儀表學(xué)報(bào),2005,26(1):41-44.YANG Xu-hai,ZHAI Hui-sheng,HU Yong-hui,et al.Study on GPS Disciplined Rb Clock Based on New Frequency Accuracy Measurement Algorithm[J].Chinese Journal of Scientific Instrument,2005,26(1):41-44.(in Chinese)

[5]ValdesM D,Villares IFarina J,Moure M J.A FPGA-based frequency measurement system for high—accuracy QCM sensors[J].Industrial Electronics,2008,34(11):1707-1712.

[6]Cadenas O,Megson G.Pipelining considerations for an FPGA case[C]//Proceedings of 2001IEEE International Symposium on Digital Systems Design.Warsaw,Poland:IEEE,2001:276-283.