齊 震,程浪浪,劉 強,王 昕,趙 凡,高玉平,馮 平
(1. 中國科學院國家授時中心,西安 710600; 2. 封開低頻時碼授時臺,廣東 肇慶 526500; 3. 中國科學院大學,北京 100049)
低頻時碼授時技術(shù)屬于長波授時技術(shù)的一種擴展,是國際電信聯(lián)盟(International Telecommunication Union, ITU)一直提倡的一種授時技術(shù)[1]。目前,已有多國建有低頻時碼授時系統(tǒng),如中國BPC、美國WWVB及日本JJY等[2]。低頻時碼授時的優(yōu)勢在于借助低頻載波通過天波和地波進行時間信號傳輸,具有較高的相位穩(wěn)定性,傳播損耗小,可以遠距離傳輸[1]。我國于2007年建成的商丘低頻時碼授時臺,采用幅度調(diào)制體制,調(diào)制速率為1 bit/s,通過4種不同脈寬傳遞時間信息,在每秒開始時刻使幅度下跌90%標志UTC(NTSC)秒的起始時刻[3-4]。其中,低頻時碼授時系統(tǒng)的系統(tǒng)時間與UTC(NTSC)的時間同步通過衛(wèi)星共視實現(xiàn)。
隨著數(shù)字電子技術(shù)的不斷發(fā)展,低頻時碼接收機也朝著數(shù)字化、小型化的趨勢發(fā)展,所采用的定時手段也趨向數(shù)字化解決方案[5]。低頻時碼授時信號從發(fā)播臺至用戶接收終端的過程中,會受到各種環(huán)境的干擾以及信號傳播鏈路時延的影響,因此接收端往往接收到的低頻時碼秒脈沖會存在隨機抖動現(xiàn)象,在一定程度上影響定時精度[6]。為提高定時精度,需對秒脈沖的抖動進行平滑。本文將利用接收到的我國商丘低頻時碼發(fā)播臺發(fā)播的BPC授時信號,分別對最小均方誤差(least mean square, LMS)自適應濾波和卡爾曼濾波這兩種數(shù)字濾波方法在低頻時碼秒脈沖抖動平滑處理中的平滑性能進行仿真驗證分析。
BPC授時信號到達接收終端后通過解調(diào)得到的1PPS與國家授時中心標準時間UTC(NTSC) 1PPS之間存在一個由于傳播路徑及設備時延等因素產(chǎn)生的固定相位差,該相位差可通過理論模型計算和接收終端設備標定獲得。同時,接收終端的內(nèi)部工作時鐘源產(chǎn)生的本地1PPS與BPC 1PPS之間存在一個隨機的相位差,該相位差是由于設備上電瞬間本地1PPS初始相位隨機產(chǎn)生導致,可通過接收終端內(nèi)部測量獲得,每次斷電重啟后會發(fā)生變化。
而在實際BPC信號接收時,受接收環(huán)境干擾等影響,解調(diào)得到的BPC 1PPS存在抖動現(xiàn)象。因此BPC 1PPS與UTC(NTSC) 1PPS之間的相位差會對應地出現(xiàn)幅度抖動,同理,BPC 1PPS與本地1PPS之間的相位差會出現(xiàn)同樣幅度變化的抖動。
所以,接收終端解調(diào)得到的BPC 1PPS的抖動現(xiàn)象可以通過接收終端內(nèi)部測量BPC 1PPS與本地1PPS的相位差波動情況進行量化反映。對該抖動的相位差進行濾波平滑處理后,結(jié)合已標定的時延對本地1PPS進行修正,使其與UTC(NTSC) 1PPS保持同步,接收終端即可輸出標準穩(wěn)定的1PPS。原理如圖1所示。
圖1 BPC接收終端秒脈沖同步修正原理圖Fig.1 BPC receiver 1PPS synchronization correction schematic diagram
BPC接收終端秒信號修正原理如圖2所示。天線接收BPC授時信號并解調(diào),然后得到BPC 1PPS(one pulse per second),同時本地高穩(wěn)晶振提供的工作頻率信號經(jīng)分頻倍頻得到本地1PPS,再對上述兩個1PPS的脈沖邊沿進行捕獲,利用時鐘計數(shù)器得到兩個脈沖之間的時差值,該時差值即為BPC 1PPS和本地1PPS之間的相位差。對獲得的時差值進行數(shù)字濾波,將濾波結(jié)果作為本地1PPS進行第一次移相修正,然后再附加時延修正,即可輸出與UTC(NTSC) 1PPS保持同步的1PPS。
圖2 BPC接收終端秒信號修正原理圖Fig.2 BPC receiver second signal correction schematic diagram
1.2.1 LMS自適應濾波
本文采用基于最小均方誤差(least mean square, LMS)算法的自適應濾波方法,是一種改進的最陡下降算法,不需要對信號提取先驗統(tǒng)計,具有計算量小、易于實現(xiàn)的特點[7]。如圖3所示自適應濾波器原理,輸入信號x(n)經(jīng)過可調(diào)節(jié)參數(shù)濾波器后得到輸出信號y(n),再與期望信號d(n)對比后得到誤差信號e(n),然后通過自適應算法對濾波器參數(shù)進行調(diào)整,最終使e(n)的均方值最小[8]。自適應濾波利用前一時刻已得的濾波器參數(shù)的結(jié)果,自動調(diào)節(jié)當前時刻的濾波器參數(shù),以適應信號和噪聲隨時間變化的統(tǒng)計特性,從而實現(xiàn)最優(yōu)濾波[9]。
圖3 LMS自適應濾波器原理圖Fig.3 Principle of LMS adaptive filter
所以,LMS自適應濾波算法主要分為兩個部分:濾波和和自適應調(diào)整。其算法過程如下[7]
(1)
e(n)=d(n)-y(n)
(2)
w(n+1)=w(n)+μe(n)x(n)
(3)
式中,N為濾波器的階數(shù);μ為步長因子,取值需要滿足公式(4),即0<μ<1/λmax,λmax為輸入信號自相關矩陣的最大特征值。
|1-2μλmax|<1
(4)
根據(jù)LMS自適應濾波原理,在對秒脈沖進行濾波處理時需設置4個輸入信號,分別為時差信號、時差期望信號、步長因子信號和濾波器階數(shù)信號。時差信號為BPC 1PPS和本地1PPS之間的相位差;時差期望信號為一定時間內(nèi)時差信號的均值,同時利用該時間段內(nèi)的時差信號計算自相關矩陣的最大特征值,從而確定步長因子;濾波器階數(shù)可根據(jù)硬件資源情況確定,一般為8階,而本文通過仿真進一步分析了不同階數(shù)對濾波性能的影響。
1.2.2 卡爾曼濾波
卡爾曼濾波利用遞歸算法解決時域離散數(shù)據(jù)的濾波問題,是基于最小誤差準則的最優(yōu)線性估計[10-11]??柭鼮V波可對平穩(wěn)和非平穩(wěn)隨機信號進行較高精度的估計,由于系統(tǒng)的當前狀態(tài)只與前一時刻的系統(tǒng)狀態(tài)以及當前測量值有關,所以無需過去所有的測量值,只需前一時刻的系統(tǒng)狀態(tài)與當前時刻的測量值就可以對當前狀態(tài)實現(xiàn)最優(yōu)估計[12-13]??柭鼮V波是通過預測和更新兩個階段實現(xiàn)對系統(tǒng)的最優(yōu)估計,基本模型如下[14]。
首先構(gòu)建一個離散線性系統(tǒng)模型,見式(5)、式(6)[15]
xk=A×xk-1+B×uk+ωk-1
(5)
zk=H×xk+vk
(6)
其中,xk為k時刻系統(tǒng)的狀態(tài)矩陣;xk-1為k-1時刻系統(tǒng)的狀態(tài)矩陣;A為狀態(tài)轉(zhuǎn)移矩陣;B為控制輸入矩陣;uk為系統(tǒng)控制輸入;ωk-1為過程噪聲;zk為k時刻系統(tǒng)的觀測量;H為狀態(tài)觀測矩陣;vk為系統(tǒng)噪聲。
基于該離散系統(tǒng)模型,假設過程噪聲與系統(tǒng)噪聲均為高斯白噪聲,即P(ω)∈N(0,Q),P(v)∈N(0,R),其中Q和R分別為ω和v的協(xié)方差??傻妙A測方程為[15]
(7)
(8)
更新方程為
(9)
(10)
(11)
在秒脈沖卡爾曼濾波處理中,由于信號在高斯信道中傳輸,當前時刻系統(tǒng)的預測值為上一時刻系統(tǒng)的輸出值,所以令狀態(tài)轉(zhuǎn)移矩陣A設為1,控制輸入矩陣B設為0,狀態(tài)觀測矩陣H設為1,單位矩陣E設為1。故式(7)~(11)可更新為
(12)
(13)
(14)
(15)
(16)
本文利用已有的BPC接收設備所采集的BPC 1PPS和本地1PPS時差作為濾波模塊的輸入信號進行MATLAB仿真分析。該接收設備的本地1PPS由JKOC36A-10 MHz型號的恒溫晶振提供,其頻率準確度為±0.03×10-6,日老化率為±0.5×10-9。數(shù)據(jù)選用23點后采集的1 h數(shù)據(jù),共3 600個采樣點。
LMS自適應濾波的仿真結(jié)果如圖4~圖7所示。為研究濾波器不同階數(shù)對濾波結(jié)果的影響,本文分別從8階、16階、32階和64階分析LMS自適應濾波的性能。根據(jù)仿真結(jié)果發(fā)現(xiàn),LMS自適應濾波算法可以有效平滑秒脈沖的抖動現(xiàn)象,并且隨著濾波器階數(shù)的增加,平滑效果越好。對濾波結(jié)果一次差分如圖8~圖11所示,在濾波器階數(shù)為8階時,相鄰采樣點之間的波動在1 ms左右,甚至在部分時間段波動超出1 ms。當階數(shù)增加至16階時,平滑效果稍有改善,隨著階數(shù)繼續(xù)增加,濾波結(jié)果的穩(wěn)定性也明顯改善。分別統(tǒng)計8階、16階、32階和64階LMS自適應濾波器濾波結(jié)果標準差以及原始數(shù)據(jù)標準差如表1所示,原始數(shù)據(jù)標準差結(jié)果為2.894 ms,經(jīng)LMS自適應濾波處理后,標準差可以提升到0.569 ms。如圖12所示的LMS自適應濾波器的階數(shù)同濾波結(jié)果標準差之間的變化趨勢可以發(fā)現(xiàn),標準差的變化趨勢隨著階數(shù)變大而慢慢趨于平穩(wěn),最終穩(wěn)定在0.5 ms附近。統(tǒng)計濾波前后均值可以發(fā)現(xiàn),與原始數(shù)據(jù)均值相差最大的是采用8階濾波,相差0.201 ms,與原始數(shù)據(jù)均值相差最小的是采用32階濾波,相差僅0.002 ms。所以,在實際應用中采用LMS自適應濾波應當選擇合適的濾波階數(shù),以達到較好的濾波效果。雖然隨著濾波階數(shù)的增加,會獲得更好的濾波效果,但是在電路實現(xiàn)中,多階數(shù)濾波對應消耗更多的電路邏輯資源。綜上所述,LMS自適應濾波可以有效地平滑低頻時碼秒脈沖的抖動,但是在實際應用中應根據(jù)濾波性能和資源使用情況合理選擇濾波階數(shù)。
表1 LMS自適應濾波結(jié)果的標準差及均值統(tǒng)計Tab.1 Standard deviation and mean statistics of LMS adaptive filtering results ms
圖4 8階LMS自適應濾波仿真結(jié)果Fig.4 Simulation results of 8-order LMS adaptive filtering
圖5 16階LMS自適應濾波仿真結(jié)果Fig.5 Simulation results of 16-order LMS adaptive filtering
圖6 32階LMS自適應濾波仿真結(jié)果Fig.6 Simulation results of 32-order LMS adaptive filtering
圖7 64階LMS自適應濾波仿真結(jié)果Fig.7 Simulation results of 64-order LMS adaptive filtering
圖8 8階LMS自適應濾波仿真結(jié)果一階差分Fig.8 First-order difference of 8-order LMS adaptive filtering simulation results
圖9 16階LMS自適應濾波仿真結(jié)果一階差分Fig.9 First-order difference of 16-order LMS adaptive filtering simulation results
圖10 32階LMS自適應濾波仿真結(jié)果一階差分Fig.10 First-order difference of 32-order LMS adaptive filtering simulation results
圖11 64階LMS自適應濾波仿真結(jié)果一階差分Fig.11 First-order difference of 32-order LMS adaptive filtering simulation results
圖12 LMS自適應濾波器階數(shù)與濾波結(jié)果標準差關系圖Fig.12 The relationship between the order of LMS adaptive filter and the standard deviation of filtering results
卡爾曼濾波仿真結(jié)果如圖13所示,可以發(fā)現(xiàn)卡爾曼濾波在平滑秒脈沖抖動應用上具有非常明顯的效果。圖中紅線為濾波結(jié)果,收斂后的結(jié)果接近于一條直線。卡爾曼濾波結(jié)果的一階差分結(jié)果如圖14所示,在濾波的初始階段大約有近700 s的收斂時間,但收斂后的結(jié)果表現(xiàn)十分穩(wěn)定。計算已收斂時間段內(nèi)的均值和標準差如表2所示,標準差由2.894 ms提升至0.084 ms,均值稍有0.342 ms的偏差。
圖13 卡爾曼濾波仿真結(jié)果Fig.13 Kalman filter simulation results
圖14 卡爾曼濾波仿真結(jié)果一階差分Fig.14 First-order difference results of Kalman filter simulation
對比LMS自適應濾波和卡爾曼濾波在秒脈沖平滑處理的結(jié)果,發(fā)現(xiàn)卡爾曼濾波的效果要優(yōu)于LMS自適應濾波,但是卡爾曼濾波存在一定時間的收斂過程,所以濾波響應時間上,LMS自適應濾波要優(yōu)于卡爾曼濾波。如圖15所示,利用Allan方差對比兩種濾波結(jié)果的頻率穩(wěn)定度可以發(fā)現(xiàn),卡爾曼濾波的結(jié)果明顯優(yōu)于LMS自適應濾波。所以在實際應用中,可根據(jù)需求合理選擇濾波策略進行秒脈沖抖動處理,若對定時穩(wěn)定性及精度要求較高可以選擇卡爾曼濾波,若對系統(tǒng)響應時間要求較高,而對定時精度沒有過高需求可以選擇LMS自適應濾波。
圖15 原始數(shù)據(jù)與不同濾波方法的Allan方差Fig.15 Allan variance of original data and different filtering methods
本文通過研究BPC 1PPS抖動平滑原理,并利用BPC接收機實際采集的BPC 1PPS和本地1PPS相位差數(shù)據(jù),對LMS自適應濾波和卡爾曼濾波秒脈沖抖動平滑仿真分析,得到以下結(jié)論:
1)BPC信號在傳播過程中受環(huán)境等因素的影響而發(fā)生的秒脈沖抖動可以通過BPC 1PPS和本地1PPS的相位差反映,對該相位差抖動平滑后用于校準本地1PPS可以有效提高BPC定時精度和穩(wěn)定性。
2)卡爾曼濾波和LMS自適應濾波對于低頻時碼秒脈沖的抖動平滑都具有明顯效果,卡爾曼濾波的平滑效果和頻率穩(wěn)定度表現(xiàn)更優(yōu),收斂后濾波結(jié)果標準差為0.084 μs。LMS自適應濾波隨著濾波器階數(shù)的增加,濾波效果也不斷改善,最終濾波結(jié)果標準差穩(wěn)定在0.5 ms附近。
3)上述兩種濾波方法有各自的優(yōu)勢和缺陷,LMS自適應濾波的濾波結(jié)果響應時間快,但濾波性能受階數(shù)的影響,而高階濾波器在電路實現(xiàn)中會消耗更多資源;卡爾曼濾波具有較好的濾波性能,但需要一定的時間進行收斂。所以,在秒脈沖抖動平滑的實際應用中,應根據(jù)實際需求合理選擇濾波器。