卡爾曼濾波在光纖授時時延測量中的應用

王 超，楊君剛，王 凱（西安通信學院信息服務系，西安710106）

?

中文核心期刊

卡爾曼濾波在光纖授時時延測量中的應用

王 超，楊君剛，王 凱
（西安通信學院信息服務系，西安710106）

摘要：針對Round-t ri p時延測量方法中時延測量值受光功率波動和時鐘頻率漂移所引起的時延波動范圍大、測量值準確性低的問題，提出利用卡爾曼濾波對時延測量值進行處理的方法，合理設計了濾波器相關參數。通過M at l ab數值仿真表明，濾波后的時延估計值與原始測量值相比精度約提高8倍，證明了應用所提方法可以得到準確的時延值。在此基礎上，將該方法應用于光纖高精度授時設備中，實測結果表明該方法比統計平均的時延處理方法獲得的系統授時精度提高約130ps。

關鍵詞：光纖授時；授時精度；時延測量；卡爾曼濾波

0　引言

高精度的時間同步網絡作為國家最基礎、最關鍵的基礎設施之一，在確保國家安全方面具有重要的支撐作用，已經成為一種重要的戰略資源。目前，我國高精度授時仍然以GPS為主，北斗衛星導航系統的授時服務尚待進一步完善，而且衛星授時服務系統受其技術體制限制，存在一定的可靠性和安全性隱患，所以有必要對基于其它手段的高精度授時服務系統進行研究，以提高我國整個授時服務體系的可靠性、安全性和完備性。分布廣泛的光傳輸網憑借其低成本、高穩定性和抗干擾性等優勢被用來承載授時業務具有明顯的優勢，可以作為衛星授時服務系統的有效補充和備份手段［1］。由文獻［2］可知，在光纖授時中，要實現主從授時節點之間的時間同步，關鍵在于對時間信息的傳遞時延進行精確的測量和準確補償。本文就針對此問題，通過對經典的Round-trip時延測量法中造成時延測量值波動的主要原因進行分析，設計了卡爾曼濾波算法對時延測量值進行濾波處理。

1　Round-trip時延測量方法

目前對于主從站間傳輸時延的精確測量主要采用雙向法，包括雙向比對法和Round-trip法［3］。其中在雙向比對法中近似認為兩個方向的傳輸時延相等，相較于Round-trip法，雙向對比法精度較低，因此下面主要對Round-trip時延測量方法進行介紹。

圖1　Round-trip時延測量方法

1.1Round-trip時延測量原理

如圖1所示，Round-trip時延測量方法在文獻［4］中被首次提出，主從站通過波分設備及光纖進行連接。時延測量基本過程如下：首先，主站的基準時鐘產生參考時間1PPS信號，經光發送端調制分別以波長λ1和λ2發送到從站，同時將調制前的1PPS上升沿作為時間間隔測量芯片（TIC1）的時延測量啟動信號；接收并恢復出從站以波長λ2傳回的1PPS信號，將其上升沿作為TIC1的測量停止信號，即可完成主從站間往返時延τsum的測量。在從站，以波長λ1和λ2傳輸過來的信號通過光接收端解調恢復出1PPS信號。兩個波長的傳輸群時延不同，因此存在一定的時延差。假定λ1＞λ2，則波長λ1通道恢復出的1PPS信號被作為TIC2的測量啟動信號，同時將該信號經光發送端以波長λ2回傳主站；以波長λ2通道恢復出1PPS信號被作為TIC2的測量停止信號，即可完成時延差τdiff的測量。

假設波長λ1和波長 λ2的單程時延分別為 τ1和τ2，則通過上述過程測量得到的往返時延τsum和時延差τdiff即可求得τ1。然后通過將主站的時間信息和求得的單程傳輸時延τ1傳至從站，通過對從站進行精確的延時控制即可實現主從站的時間同步。具體延時控制方法在此不再贅述。對于時延測量中用到的TIC芯片，目前主要以采用TDC-GP系列芯片為主，測量分辨率從百皮秒到幾十皮秒不等。因此在本方法中，在傳輸距離一定的條件下，需要選擇合適的波長，以確保時延差能夠被TIC精確測量［5］。

1.2時延測量值分析

由上面的時延測量原理可知，只要TIC芯片測量得到的時延值足夠準確，就可以實現主從站的時間同步，保證授時系統的授時精度，然而，從實際測試結果可知，TIC芯片測量得到的時延值往往存在較大的波動［6］。分析其產生的原因主要有兩個方面：一是由于光纖傳輸鏈路環境（溫度、應力等）的變化造成傳輸時延的波動，主要影響時延的長期穩定性；二是由于授時設備中的光收發模塊的發送/接收光功率的波動和本地時鐘頻率漂移所帶來的時延波動，主要影響時延的短期穩定性。由于前者對時延的短期穩定性影響甚微，后者在一定條件下可能會帶來納秒級的時延波動。因此為了克服由此所帶來的時延波動，降低設備不穩定性對時延測量值的影響，需要對時延測量值進行必要的處理。目前，相關的處理算法［7，8］已經被用于類似問題的研究。本文采用卡爾曼濾波對時延測量值進行濾波處理，以獲得準確的時延值，然后將此時延值按照一定算法處理作為從站的補償值進行精確校準，實現高精度的主從站時間同步。

2　卡爾曼濾波在時延測量值處理中的應用

卡爾曼濾波是解決以最小均方誤差為準則的最佳線性濾波問題［9］，其算法迭代過程中只需要根據前一次的估計值和當前的觀測值就可以對觀測變量進行估計，適用于對時延測量值的處理。

2.1卡爾曼濾波原理

卡爾曼濾波的信號模型是由狀態方程和觀測方程組成的。設在k時刻觀測變量的真值為sk，對應的觀測值為xk，則系統的信號模型為：

其中，Ak|k-1為狀態轉移因子，表示相鄰時刻觀測變量的相關程度，wk為系統獨立白噪聲，Hk為觀測參數，nk為觀測白噪聲 （由于測量或其它原因引起的隨機誤差）。wk和nk的統計特性為：

利用卡爾曼濾波對觀測變量進行估計時，需要前一時刻的估計值和當前時刻的測量值作為先驗知識，然后進行遞推估計。算法具體實現步驟如圖2所示。

2.2卡爾曼濾波器設計

圖2　卡爾曼濾波算法實現步驟

在上面介紹的Round-trip時延測量方法中，涉及到的測量數據包括往返時延和時延差。在此以往返時延的測量為例，對卡爾曼濾波器的設計進行說明。在主站，我們利用TIC芯片對主站產生的1PPS信號與從站返回的1PPS信號的時間間隔（即往返時延）進行連續測量。假設sk是k時刻往返時延的真值，xk為包含系統噪聲和測量噪聲的往返時延測量值，系統的維數為1。

在保持連接配置關系不變的情況下，即便考慮到溫度、應力等因素對往返時延的影響，其測量值應該滿足短期穩定性。因此，狀態方程中的狀態轉移因子Ak|k-1和觀測方程中的觀測參數Hk均設置為1。在本設計中，受溫度、應力等因素影響的系統獨立白噪聲wk對時延值的影響較小，同時為了避免造成時延估計值整體偏大或偏小，一般依據經驗對協方差R取值為1。對于由光功率變化、時鐘頻率漂移以及TIC芯片本身測量分辨率引起的觀測噪聲nk的協方差Q的取值，一般通過實驗獲得，取值的依據是要使濾波后的時延值波動范圍與實際相符［10］。

在開始卡爾曼濾波迭代運算前，需要給定狀態的初始估計值s?0和誤差協方差M0。一般情況下，通過預先采集的一組時延測量值的統計量得到這兩個初始值，以加快濾波器的收斂速度。

全部相關參數確定后，就意味著卡爾曼濾波器的設計已經完成。在接下來的處理中，只需按照上述卡爾曼濾波過程進行迭代，就可以得到濾波后的時延估計值。我們將處理后得到的準確時延值，應用到從站的時間補償中就可以降低同步的誤差、提高授時的精度。

3　數值仿真與實際應用

3.1數值仿真與分析

為了驗證卡爾曼濾波算法及濾波器參數的正確性，在實驗室條件下，本文以兩臺光纖高精度授時設備作為主從站，搭建了200km傳輸距離的光纖高精度授時系統。按照Round-trip時延測量方法，在主站通過SR620對往返時延值進行測量，測量時間約30分鐘，測量周期為1次/秒。本文對SR620測量所得往返時延值利用上述設計的卡爾曼濾波器對進行濾波處理，得到的數值仿真結果如圖3、圖4所示。

圖3　200　km傳輸距離下往返時延測量值

由圖3可知，SR620直接測量的往返時延值數據波動范圍較大，經統計，其峰峰值和標準偏差分別為8.32ns和0.89ns。假如將此時延測量值直接用于從站的補償校準，必然會帶來較大的同步誤差，降低授時的精度。如圖4所示，經過濾波處理后的往返時延估計值波動范圍明顯減小，往返時延估計值在800點后基本收斂，收斂后的峰峰值和標準偏差分別為0.66ns 和0.11ns，逐步趨近于常量1，979，218ns。

通過以上數值仿真及分析可知，經過卡爾曼濾波處理可以有效消除由于設備不穩定性和測量精度所造成的時延抖動，獲得準確的時延值用于后期從站的時延補償。

3.2實際應用與結果

為了進一步驗證對時延測量值進行卡爾曼濾波處理有助于提升光纖高精度授時設備的同步性能，我們將某光纖高精度授時設備中的時延處理方法（統計平均方法）改用卡爾曼濾波方法。系統配置采用100km光纖進行連接，授時系統正常運行1小時以上。由實際測量結果對比發現，使用卡爾曼濾波對時延測量值進行處理后，主從站同步性能得到了有效改善，同步誤差均值、標準偏差分別減小242.5ps、133.1ps，進而提高了光纖授時系統的授時精度。

圖4　卡爾曼濾波處理后往返時延估計值

4　結束語

本文通過對光纖高精度授時設備在采用卡爾曼濾波時延處理方法前后同步性能測試結果表明，本文所提出的利用卡爾曼濾波算法消除光纖授時中時延波動的方法適用于實際的光纖授時系統，可以有效提升光纖高精度授時設備的同步性能，提高光纖授時系統的授時精度。

參考文獻：

［1］李恩，羅青松，劉志強，等.光纜網高精度授時技術現狀及應用［J］.中國電子科學研究院學報，2013，8（4）：344-349.

［2］高超，王波.基于光纖鏈路的高精度時間頻率傳輸與同步［J］.科技導報，2014，32（34）：41-46.

［3］丁小玉，盧麟，張寶富，等.光纖Round-trip法授時的誤差分析［J］.激光與光電子學進展，2010，47（4）：040603-1-040603-5.

［4］ATSUSHI IMAOKA，MASAMI KIHARA.Accutate time/frequency transfer method using bidirectional WDM transmission［J］.IEEE Trans.Instrum.Meas，1998，47（2）：537-542.

［5］付永杰，才瀅.光纖雙向時間傳遞的誤差分析［J］.計量學報，2014，35 （3）：276-280.

［6］TSENG WH，LIN SY.A survey of time transfer via a bidirectional fiber link for precise calibration services［J］.NCSLI Measure J.Meas.Sci，2013，8 （2）：70-77.

［7］邢開亮，尹義蓉，黃永華，等.提高網絡授時精度的卡爾曼濾波方法［J］.儀器儀表學報，2011，32（4）：949-954.

［8］張妍，邱天爽，任福全.分布式光纖中的穩健時間延遲估計方法［J］.應用光學，2012，33（4）：816-820.

［9］趙樹杰，趙建勛.信號檢測與估計理論［M］.北京：清華大學出版社，2005.

［10］杜文建.GPS與恒溫晶振互補的高穩定時鐘研究［D］.南京：南京大學，2012：49-50.

中圖分類號：TN929.11

文獻標識碼：A

文章編號：1002-5561（2016）06-0034-04

DOI：10.13921/j.cnki.issn1002-5561.2016.06.010

收稿日期：2015-12-30。

作者簡介：王超（1990-），男，碩士研究生，主要從事光纖時間傳遞技術的研究。

Application of Kalman filter on
propagation delay measurement in fiber time service

WANG Chao，YANG Jun-gang，WANG Kai
（Department of Information service，Xi'an Communications Institute，Xi'an 710106，China）

Abstract：For the problems of fluctuation and low-accuracy of measurements caused by optical power fluctuations and clock frequency drift in Round-trip propagation delay measurement method，a method have been proposed to reduce the effects of fluctuation and improve accuracy of the measurements，presented the reasonable design of the filter parameters.Numerical simulation result by Matlab shows that the precision of propagation delay estimates were improved eight times compared with the original measurements，proved that the application of the proposed method can get the exact propagation delay value；Then the method was applied to a high precision timing devices based on optical fiber，test results show that the system timing precision increased more than 130ps compared with the processing method for propagation delay based on statistical average.

Key words：fiber timing，timing accuracy，propagation delay measurement，Kalman filter