紅外光譜預處理去噪研究范文

本站小編為你精心準備了紅外光譜預處理去噪研究參考范文，愿這些范文能點燃您思維的火花，激發您的寫作靈感。歡迎深入閱讀并收藏。

《海軍航空工程學院學報》2016年第5期

摘要：

采用小波包變換軟閾值去噪方法對潤滑油的紅外光譜進行去噪預處理，結合偏最小二乘法進行定量預測，并與常用的平滑去噪方法進行對比。結果表明，小波包變換能有效地去除紅外光譜的噪聲，以此為基礎建立的潤滑油酸值模型的預測精度高于常用的平均平滑法和Savitzky-Golay卷積平滑法這2種平滑去噪方法的預測精度。

關鍵詞：

小波包變換；紅外光譜；平滑；偏最小二乘法

紅外光譜分析技術由被測樣品的紅外光譜主導，由紅外光譜儀得到的光譜信號中不僅含有樣品的信息，還包含了噪聲和各種外界干擾因素[1]。為了使定量分析模型更加穩健和精確，需對光譜進行去噪預處理[2-3]。苗福生[4]、吳海云[5]等對常用的平均平滑方法和Savitzky-Golay卷積平滑法這2種平滑方法做出了比較研究。本文以某特種車的潤滑油為研究對象，研究了小波包變換[6]對紅外光譜數據的去噪預處理，對比了兩種平滑方法下利用偏最小二乘法（PLS）[7]對酸值進行建模的預測精確度。

1實驗與方法

1.1樣品的光譜測定

采用BRUKERRT-DLaTGS紅外光譜儀[8]。收集潤滑油樣本21個，采用固定光程為100mm的ZnSe樣品池，以空氣為參比進行光譜掃描，掃描10次取平均值，減少光譜采集過程中隨機因素干擾，保證光譜精度。潤滑油樣品的原始紅外光譜圖如圖1所示。

1.2樣品的酸值測定

潤滑油的酸值是指中和1g油液中全部酸性組分所需要的堿量，以mgKOH/g表示。潤滑油的酸值測定選擇顏色指示劑法[9]，采用GB/T264-83石油產品酸值測定法對所選的潤滑油樣品進行酸值的測定。實驗時，對同一樣品進行2次實驗取平均值。實驗樣品的酸值分布見圖2。其中最大值為0.060mgKOH/g，最小值為0.034mgKOH/g。均符合要求的范圍，它們的平均值為0.042mgKOH/g，標準差為0.0074。

2平滑預處理

2.1平均平滑法

根據紅外光譜圖中的吸收峰是否消失或采用偏最小二乘法（PLS）建模后得出的均方根誤差（RMSE）確定出預處理方法的最佳參數[10]。平均平滑的數據點數通常從20以內的奇數中選擇，平滑的點數越高，光譜越平滑，去噪效果越好，但所得光譜的分辨率越低，光譜的有些肩峰會消失[11]。考察范圍為從5點開始。間隔為2點，分別進行平滑，結果如圖3~7所示。當平滑點數為11點時，波數為967處的肩峰開始消失，當平滑點數為13點時，波數為967處的肩峰完全消失，因而樣品的平均平滑點數選擇9點。

2.2Savitzky-Golay平滑法

Savitzky-Golay卷積平滑法是1964年Savitzky和Golay提出并廣泛運用于信號濾波處理的平滑方法[12]，該方法與移動平均平滑法的基本思想是類似的，只是沒有進行簡單的平均，而是通過多項式來對移動窗口內的數據進行多項式最小二乘擬合，其本質是一種加權平均法，強調中心點的中心作用[13]。采用Savitzky-Golay平滑，平滑的窗寬值是一個重要參數。如果窗寬值選擇較小，平滑去噪效果不好；如果窗寬值選擇較大，則會造成光譜信號失真。因此，要選擇合適的窗寬值，考察范圍為7~23，間隔為2，在選擇不同的窗寬值后發現紅外光譜圖的特征峰均還存在，因而需要用偏最小二乘法對不同的窗寬值的平滑結果建立模型后比較均方根誤差（RMSE）的值來確定出最佳的窗寬值，結果如圖8所示。當窗寬值為11時，均方根誤差（RMSE）最小，所以樣品的Savitzky-Golay平滑窗寬值選擇11。

3小波包變換預處理

小波包變換是基于小波變換的進一步發展，能夠提供比小波變換更高的分辨率。小波包分解與小波分解相比，是一種更精細的分解方法[14]。在多分辨分析中，L2(R)=⊕j∈ZWj，表明多分辨分析是按照不同的尺度因子j把Hilbert空間L2(R)分解為所有子空間Wj（j∈Z）的正交和，其中，Wj為小波函數φ(t)的閉包（小波子空間）。小波包分析就是進一步對小波子空間Wj按照二進制分式進行頻率的細分，以達到提高頻率分辨率的目的。

3.1小波包降噪的步驟

小波包分析的一般步驟[15]如下。1）信號的小波包分解。選擇一個小波并確定小波分解的層次N，然后對信號進行小波包分解。2）確定最優小波包基。對于一個給定的熵標準，計算最優樹。3）小波包分解系數的閾值量化。對于每一個小波包分解系數，選擇一個恰當的閾值并對系數進行閾值量化。4）信號的小波包重構。根據最底層的小波包分解系數和經過量化處理的系數進行小波包重構。

3.2小波包變換參數選擇

本文運用小波包變換對紅外光譜圖進行去噪處理，采用了軟閾值和“db3”小波包基[16]。但信號尺度的分解層數是影響去噪效果的一個重要因素，一般情況下分解層數較少，去噪效果不理想；但分解層數較多，導致運算量增大，且會造成信息的丟失。通常分解層數在3層或3層以上，所以選擇從考察3層開始考察，分別進行小波包變換去噪，結果如圖9~11所示。在分解層數為4層時，波數為1422處的吸收峰已經基本消失，在分解層數為5層時，波數為1422處的吸收峰已經完全消失，所以信號尺度的分解層數選擇3層。

4預處理方法的對比分析

4.1建模及預測

通過上述3種方法分別對獲得的紅外光譜進行預處理，再采用偏最小二乘法（PLS）對這3種預處理后的紅外光譜進行總酸值建模預測，并以相對誤差值（RE）作為評價標準，結果如圖12和表1所示。圖12中，“○”線、“□”線、“+”線分別是是平均平滑法（9點）、Savitzky-Golay平滑法（窗寬值11）和小波包變換預處理后建模的預測相對誤差。

4.2對比分析

樣本的預測相對誤差分布在1.1%~5.5%之間，大部分樣本在利用小波包變換進行紅外光譜去噪預處理后建模的相對誤差較其余兩種平滑去噪方法小，樣本的相對誤差超過4.0%的有2個，大部分建模樣本的相對誤差在3.0%以內，符合模型預測要求。

5結論

由于紅外光譜分析中測量得到的光譜信息可能受到來自各方面因素的影響，在建立定量分析模型時，對光譜進行去噪預處理是必要的。本文運用小波包變換去噪方法對紅外光譜進行了去噪預處理，對比了常用的平均平滑法和Savitzky-Golay平滑法2種平滑去噪方法。結果表明，在對潤滑油樣品的紅外光譜進行去噪預處理時，小波包變換不僅能夠有效地去除紅外光譜的噪聲，而且結合偏最小二乘法所建立的酸值模型預測精度高，相對誤差在1.1%~5.2%之間。因此，針對潤滑油的酸值定量預測，小波包變換去噪預處理方法的優勢明顯。

參考文獻：

[1]高榮強，范世福，嚴衍祿，等.近紅外光譜的數據預處理研究[J].光譜學與光譜分析，2004，24（12）：1563-1564.

[2]張銀，周孟然.近紅外光譜分析技術的數據處理方法[J].紅外技術，2007，29（6）：345-348.

[4]苗福生，馬毅，汪西原，等.不同預處理方法對PLS模型檢測鮮長棗糖度的影響[J].寧夏大學學報，2011，32（2）：130-133.

[5]吳海云，劉洋，左月明，等.近紅外光譜數據分析方法的研究進展[J].農產品加工，2010（3）：76-79.

[6]郎文杰.基于小波包分析的圖像融合去噪方法[J].河北理工大學學報，2010，32（2）：42-44.

[7]盧小泉，陳晶，周喜賓，等.化學計量學方法[M].北京：科學出版社，2013：15-20.

[8]翁詩甫.傅里葉變換紅外光譜儀[M].北京：化學工業出版社，2005：35-40.

[10]夏俊芳，李培武，李小昱，等.不同預處理對近紅外光譜檢測臍橙VC含量的影響[J].農業機械學報，2007，38（6）：107-111.

[11]徐廣通，袁洪福，陸婉珍.CCD近紅外光譜譜圖預處理方法研究[J].光譜學與光譜分析，2000，20（5）：619-622.

[12]趙杰文，張海東，劉木華.簡化蘋果糖度預測模型的近紅外光譜預處理方法[J].光學學報，2006，26（1）：137-141.

[13]孫輝，楊其俊，裴峻峰.小波包變換在往復泵活塞狀態監測中的應用[J].石油大學學報：自然科學版，1999，23（3）：60-62.

[14]李自國，郝偉，李凌均.基于小波包分解和支持向量數據描述的故障診斷方法[J].機械強度，2007，29（3）：365-369.

[15]段禮祥，張來斌，王朝暉，等.柴油機振動信號的小波包奇異值降噪[J].中國石油大學學報：自然科學版，2006，30（1）：93-100.

[16]宦克為，劉小溪，王欣，等.基于連續小波變換參數選擇的小麥近紅外光譜模型優化方法研究[J].長春理工大學學報，2014，37（5）：146-150.

作者:史令飛 瞿軍 單位:海軍航空工程學院