郭曉彤，楊秋翔

（中北大學軟件學院，山西太原030051）

光纖光柵（Fiber Bragg Grating，F(xiàn)BG）傳感器是近幾年來使用頻率高且應用廣泛的一種新型光纖傳感器。由于它具有耐腐蝕，耐高溫，可靠性高，安全性強等諸多優(yōu)越特性，很容易被制作成檢測氣體濃度、溫度以及壓力等物理量信息的傳感器[1-2]。通過感測外界物理量的變化，與光譜解調(diào)儀采集得到的反射譜中心波長的移位量進行判斷，實現(xiàn)對大型設(shè)備結(jié)構(gòu)狀態(tài)的實時監(jiān)測[3-5]。由于中心波長的移位量對應的是反射譜峰值位置的變化，所以快速準確的尋找峰值點是檢測的關(guān)鍵。在FBG數(shù)據(jù)采集過程中，光纖光柵波長信號會存在較大噪聲的干擾，造成了在峰值點附近會出現(xiàn)波動現(xiàn)象，導致了反射光譜的不對稱[6-7]。而使用傳統(tǒng)的尋峰算法[8-11]，如直接尋峰法、高斯擬合法等解決此類問題時，需要忽略波動現(xiàn)象，從而會產(chǎn)生較大的誤差，導致精確度下降；人工智能算法解決此類問題時，會導致時間運行成本過高。

針對上述問題，Wang等提出了基于Steger尋峰算法[12]，Liu等提出了自適應半峰檢測尋峰算法[13]，其它等提出的相關(guān)算法可參見文獻[14-16]?；赟teger尋峰算法，引入超高斯模型完成不對稱尋峰，但超高斯函數(shù)參數(shù)選擇使此方法精度受到影響；自適應半峰檢測尋峰算法，通過引入實驗經(jīng)驗值，對非對稱峰值進行補償，但此經(jīng)驗值的確定存在誤差，導致尋峰精度受到影響。為了彌補上述算法的不足，本文提出一種基于FPGA狀態(tài)機的非對稱高斯擬合尋峰算法。本算法采用FPGA超強的可編程并行計算來對光纖光柵波長信號進行濾波尋峰處理，狀態(tài)機為FPGA中的靈魂核心部分，基于狀態(tài)機的轉(zhuǎn)換過程中對數(shù)據(jù)處理完畢后再通過對波形的判斷及峰值的修正完成尋峰過程，避免了存儲空間的浪費，同時在速度和精度上，都得到了提升。

圖1 原始反射光譜圖

圖2 濾波后反射光譜圖

1 算法概述及原理

1.1 算法概述

由于外界環(huán)境及信號采集過程中所夾雜的噪聲影響，使得反射譜的峰值點發(fā)生偏移，很大程度上影響了尋峰的精度，所以要先對采集得到的離散數(shù)據(jù)進行濾波處理后再尋峰。算法過程為：對采集得到的數(shù)字信號濾波，完成預處理部分；然后通過狀態(tài)機的轉(zhuǎn)換，將去噪后的波形數(shù)據(jù)調(diào)節(jié)閾值進行分割，粗略估計峰值所在位置范圍；對此范圍的波形數(shù)字進行判斷，若為對稱波形，峰值則直接為高斯擬合后的峰值；反之則進行非對稱高斯擬合，對高斯擬合所得到的峰值進行補償修正。此算法基于FPGA狀態(tài)機對預處理后的數(shù)據(jù)再次進行處理，然后尋峰。不僅提高了速度并且彌補了以往常見尋峰算法對于非對稱波形的忽略問題，提高了算法的精度。

1.2 尋峰算法

1.2.1 濾波預處理

為了消除設(shè)備自身以及外界環(huán)境所產(chǎn)生的噪聲影響，首先采用巴特沃斯低通濾波對采集的數(shù)字信號做預處理。如圖1所示，分別為原始波形圖和濾波處理后的波形圖。

1.2.2 峰值區(qū)域分割采集模塊

本模塊的作用為采集濾波后的有效數(shù)據(jù)。由圖1可見，通過巴特沃斯濾波處理后的數(shù)據(jù)可以有效的剔除“毛刺”等噪聲干擾，得到一條較為平滑的曲線（見圖2），此時若直接進行尋峰計算，數(shù)據(jù)量大，且其中有大部分數(shù)據(jù)均為冗余數(shù)據(jù)，會對數(shù)據(jù)的處理以及系統(tǒng)運行的速度帶來不必要的工作量，占用資源，降低速度。所以進行數(shù)據(jù)篩選是非常必要的一步，此處可設(shè)定閾值，小于該點則不予考慮，完成初步的分峰截幅。通過大量的實驗比對調(diào)整后得出，此處設(shè)置的閾值為信號中心波長的75%處值，可有效濾除掉少量雜波跟噪聲。

1.2.3 峰值數(shù)據(jù)檢測模塊

本模塊的作用為采集峰值處的數(shù)據(jù)，為此后的尋峰算法提供更為可靠有效的數(shù)據(jù)信息。具體過程為上升、下降、錯誤三態(tài)之間進行相互轉(zhuǎn)化，并在FPGA系統(tǒng)內(nèi)分別設(shè)計了用于記錄保持穩(wěn)定同狀態(tài)的數(shù)據(jù)，上升態(tài)的起始位置數(shù)據(jù)以及下降態(tài)終止位置數(shù)據(jù)的寄存器，將其存入數(shù)組a[n]中。

具體過程為：

首先與預處理后的數(shù)據(jù)逐個進行對比，若當前數(shù)據(jù)大于設(shè)置的閾值且前一數(shù)據(jù)小于閾值，將此數(shù)據(jù)計入起始位置寄存器中，轉(zhuǎn)入下一數(shù)據(jù)；若當前數(shù)據(jù)大于閾值且大于前一數(shù)據(jù)，則開始進入波形的上升態(tài)，將當前數(shù)據(jù)與上一數(shù)據(jù)進行差值比較，差值范圍設(shè)定為K（K值為大量實驗數(shù)據(jù)對比后所確定的值），若差值小于等于K，則此時保持同狀態(tài)；若差值大于K，則認為此時波形發(fā)生波動，進行下一步判斷：當差值為正時，若前一狀態(tài)為下降態(tài)，則下一狀態(tài)為錯誤態(tài)，否則仍為上升態(tài)；若差值為負，則轉(zhuǎn)換為下降態(tài)。若當前數(shù)據(jù)小于閾值且前一數(shù)據(jù)大于閾值，判斷兩者之間的差值，若小于K，則轉(zhuǎn)入下一數(shù)據(jù)；若大于K，則過程結(jié)束，將上一數(shù)據(jù)計入結(jié)束位置寄存器中。

數(shù)據(jù)篩選過程結(jié)束后，得到的離散數(shù)據(jù)進入尋峰算法階段。

具體過程為：

1）在狀態(tài)機轉(zhuǎn)換處理后，得到了一組離散數(shù)據(jù)，將此數(shù)據(jù)進行高斯擬合，得到高斯函數(shù)的峰值點，記為A′。

2）非對稱高斯擬合。由于受到多因素影響，采集的FBG信號波長在實際情況下為非對稱，所以此處將非對稱高斯模型（ft）引入，對高斯函數(shù)擬合得到的峰值進行校驗補償。公式如下：

非對稱高斯模型 （ft）為：

式中μ為高斯擬合函數(shù)所得峰值點所對應的時間點；σ1為高斯擬合函數(shù)左側(cè)部分的時間點采樣數(shù)；σ2為高斯擬合函數(shù)右側(cè)部分的時間點采樣數(shù)。

通過方差的判斷并且由式1可推導出補償后的峰值為：

式中A為校驗補償后的峰值。

2 實驗結(jié)果與分析

本實驗采用CCD分光儀檢測法進行波長信號的解調(diào)，光譜解調(diào)儀波長檢測范圍為1 510～1 590 nm，解調(diào)分辨率為0.5 pm，精確度為1～3 pm，工作溫度為0～50℃。在每個周期內(nèi)，取其中的FBG3，F(xiàn)BG4反射譜分析，算法分別選取直接比較法，高斯擬合法，自適應半峰尋峰法，Steger尋峰算法以及本文提出的基于狀態(tài)機的非對稱高斯算法做對比實驗，通過對比誤差和速度，驗證本文提出算法的有效性。

2.1 恒溫尋峰誤差比較

將恒溫試驗箱的溫度設(shè)定為25℃，保證噪聲等其他參數(shù)的穩(wěn)定，光譜解調(diào)儀采集得到的光纖光柵理論中心波長值分別為：1 538.95 nm、1 545.02 nm。為了降低算法測量峰值的外在因素的影響，將所需算法分別采集10次進行平均，所求得的平均值與理論峰值進行比較，得出平均誤差值，如表1所示。

表1 25℃下各尋峰算法的平均誤差

通過表1，可分析得出，當外界環(huán)境穩(wěn)定時，直接比較法誤差較大，精度不高；高斯擬合法比直接比較法精度有所提高，但由于FBG反射光譜并非為對稱光譜，所以與理論峰值比較起來，誤差明顯；自適應半峰檢測尋峰法，雖然考慮到了非對稱譜型的情況，通過補償量進行峰值的修正，誤差有所降低，但是由于補償量為實驗統(tǒng)計值，精度有待提高；基于Steger尋峰算法，結(jié)合超高斯模型，精度比自適應半峰尋峰算法高；本文所提出的算法，通過狀態(tài)機的轉(zhuǎn)換進行區(qū)域分割并考慮到峰值的偏移量問題，引入非對稱高斯模型進行修正，誤差降低，精確度得到提高。

2.2 變溫尋峰誤差比較

將傳感器放入0～50℃的溫控試驗箱中，并且保證與上述恒溫實驗中有相同的參數(shù)，光譜解調(diào)儀分別5℃，10℃，15℃，20℃，25℃，30℃，35℃，40℃，45℃溫度下采集10次數(shù)據(jù)取其平均值，并用上述算法同樣在不同溫度下采集10次得到平均峰值數(shù)據(jù)進行比較，如圖3、圖4所示。

圖3 FBG3在變溫下的各算法誤差圖

圖4 FBG4在變溫下的各算法誤差

表2 變溫下各尋峰算法的平均誤差

通過圖3，圖4以及表2的分析可以得出，直接比較法無論從誤差還是穩(wěn)定性來說，效果都不是很理想；高斯擬合法由于未考慮不對稱峰譜的情況，所以受到影響，穩(wěn)定性和尋峰精度有待提升；自適應半峰檢測法尋峰和基于Steger尋峰算法精度都優(yōu)于高斯擬合法且受溫度影響??；本文提出的算法，精度與其他算法相比，平均誤差為0.40～0.55 pm之間，且穩(wěn)定性好。

2.3 噪聲尋峰誤差比較

保持溫度恒定在25℃，加入高斯白噪聲，使得加入的噪音幅度占信號幅度的范圍為0.001～0.1之間，測量10次取其平均值，計算誤差，如圖5、6所示。

圖5 FBG3在不同噪音下的各算法誤差圖

圖6 FBG4在不同噪音下的各算法誤差圖

由圖5、6分析可知，隨著噪音幅度的加大，算法的精度都會受到影響。直接比較法影響最大，高斯擬合法次之，自適應半峰尋峰法平均誤差在分別為4.21 pm、3.94 pm，基于Steger尋峰算法和本文所提出的算法受到影響產(chǎn)生的誤差都在1.8 pm以下，抗干擾能力較強。

2.4 運行時間分析

通過上述實驗，統(tǒng)計各算法的平均運行時間，可知直接比較法運行時間最短，為5 ms；自適應半峰檢測尋峰算法次之，為55 ms；本文算法時間為90 ms；基于Steger尋峰算法為122 ms；高斯擬合法運行時間最長，為180 ms。以上算法在時間上都可以滿足實際應用中波長解調(diào)的實時性要求。

3 結(jié)束語

文中首先在巴特沃斯低通濾波預處理的基礎(chǔ)上，通過狀態(tài)機之間的相互轉(zhuǎn)換完成對數(shù)據(jù)的進一步處理，引入非對稱高斯模型，對峰值進行補償修正。通過實驗，分析比較了在恒溫環(huán)境下以及不同溫度和噪音的干擾下對尋峰算法的影響。實驗結(jié)果表明，本文所提算法相較于直接比較尋峰法和高斯擬合算法在精度和穩(wěn)定性上都有顯著優(yōu)勢，而對于同樣考慮到非對稱問題的自適應半峰檢測尋峰法和Steger算法，則在精度和時間上優(yōu)于此兩種算法。綜上所述，本算法在FPGA超強的可編程并行計算對波長信號進行濾波尋峰處理的基礎(chǔ)上，既彌補了傳統(tǒng)算法對于非對稱忽略的問題，又避免了存儲空間的浪費，同時在速度和精度上，都得到了提升。

[1]Wada A，Tanaka S，Takahashi N.Optical fiber vibration sensor using FBG Fabry- Perot interferometer with wavelength scanning and Fourier analysis[J].IEEE Sensor Journal，2012，12（1）:225-229.

[2]Bodendorfer T，Muller M S，Hirth F，et al.Comparison ofdifferentpeak algorithmswith regards to spectrumetric fiber Bragg grating interrogation systems[C]//International Symposium on Optomechatronic Technologies，2009：122-126.

[3]蔡能宏，閆連山，陳娟子，等.基于可調(diào)諧F-P濾波器的光纖光柵傳感解調(diào)系統(tǒng)尋峰算法對比分析[J].儀表技術(shù)與傳感器，2013，33（4）:100-106.

[4]薛雙喜，馬艷，王向輝，等.多通道光纖布喇格光柵傳感器同步解調(diào)系統(tǒng)[J].電子設(shè)計工程，2012，20（7）:173-175.

[5]陳志軍，白劍，吳祖堂，等.光纖布喇格光柵反射譜尋峰算法優(yōu)化及比較[J].光子學報，2015，44（11）:83-88.

[6]Frederic C J，Cicero F M.Peak detection algorithm for fiber Bragg grating sensors[C]//Latin America Optics and Photonics Conference，2012:LM2A.8.

[7]王擁軍，劉永超，張靖濤，等.高速高分辨率光纖布拉格光柵傳感系統(tǒng)的解調(diào)技術(shù)[J].中國激光，2013，40（2）:151-156.

[8]沈漫，董雷，于本化，等.一種用于光纖布拉格光柵傳感的自適應尋峰算法[J]，電子測試，2014，9（6）:27-28.

[9]朱梅，張淼，胡立章，等.光纖布喇格光柵中心波長檢測中的尋峰算法[J].光通信研究，2011，167（5）:60-63．

[10]胡遼林，王亞萍，王斌.光纖光柵傳感信號尋峰誤差研究[J].西安理工大學學報，2015，31（4）:409-413.

[11]余有龍，王雪微，王浩.不同采樣方式下光纖布喇格光柵反射譜尋峰算法的分析[J].光子學報，2012，41（11）:1274-1278.

[12]王巧妮，楊遠洪.基于Steger圖像算法的光纖布拉格光柵尋峰技術(shù)[J].光學學報，2014，34（8）:133-138.

[13]劉泉，蔡林均，李政穎等.高速度高精度光纖布拉格光柵解調(diào)的尋峰算法研究[J].光電子·激光，2012，23（7）:1233-1239.

[14]鄭?？?，李根，李靜，等.基于非對稱廣義高斯模型的渡越時間確定方法[J].大慶石油地質(zhì)與開發(fā)，2015，34（4）:143-147.

[15]陳勇，楊凱，劉煥淋.多峰光纖布拉格光柵傳感信號的自適應尋峰處理[J].中國激光，2015，42（8）:192-197.

[16]白志慧，周志光，楊瑞飛，等．非對稱高斯函數(shù)的時變體數(shù)據(jù)特征跟蹤及可視化[J]．中國圖象圖形學報，2015，20（4）:576-584．