盧 敏，李小霞，2

(1.西南科技大學(xué) 信息工程學(xué)院，四川 綿陽(yáng) 621010； 2.特殊環(huán)境機(jī)器人技術(shù)四川省重點(diǎn)實(shí)驗(yàn)室，四川 綿陽(yáng) 621010)

基于EMD的反射式太赫茲光譜信號(hào)分解與處理*

盧 敏1，李小霞1，2

(1.西南科技大學(xué) 信息工程學(xué)院，四川 綿陽(yáng) 621010； 2.特殊環(huán)境機(jī)器人技術(shù)四川省重點(diǎn)實(shí)驗(yàn)室，四川 綿陽(yáng) 621010)

針對(duì)反射式太赫茲時(shí)域光譜系統(tǒng)所測(cè)時(shí)域信號(hào)存在反射峰、頻譜出現(xiàn)振蕩、頻域信噪比低等問(wèn)題，提出基于經(jīng)驗(yàn)?zāi)B(tài)分解處理太赫茲時(shí)域反射峰的方法，同時(shí)消除頻域振蕩的影響。通過(guò)獲取本征模函數(shù)并代替反射峰，增加時(shí)域信號(hào)有效長(zhǎng)度、提高太赫茲頻域分辨率，消除頻域振蕩影響。古代壁畫(huà)樣本的太赫茲反射譜實(shí)驗(yàn)結(jié)果表明，該方法能自適應(yīng)地處理反射峰，對(duì)太赫茲時(shí)域信號(hào)修復(fù)效果良好，頻域振蕩抑制效果良好。

太赫茲反射譜；經(jīng)驗(yàn)?zāi)B(tài)分解；反射峰處理

0 引言

太赫茲(THz)波是指頻率在0.3～3 THz范圍內(nèi)的電磁波，對(duì)應(yīng)的波長(zhǎng)范圍為100 μm～1 mm[1]。由于太赫茲波有許多優(yōu)越的特性，其在食品安全[2]、無(wú)損檢測(cè)[3]、生物分子探測(cè)[4]和毒品檢測(cè)[5]等方面都具有非常重要的科學(xué)價(jià)值。

太赫茲時(shí)域光譜(THz-TDS)技術(shù)是20世紀(jì)80年代由IBM公司和AT&T Bell實(shí)驗(yàn)室聯(lián)合發(fā)展起來(lái)的一種有效的相干探測(cè)技術(shù)[6]，該技術(shù)不僅能夠獲得THz脈沖的振幅，同時(shí)也能獲取相位信息。THz-TDS系統(tǒng)主要有反射式THz-TDS系統(tǒng)和透射式THz-TDS系統(tǒng)，如果需要測(cè)量的樣品是較厚的介質(zhì)或樣品吸收特別強(qiáng)，那么需要使用反射THz-TDS系統(tǒng)來(lái)對(duì)其進(jìn)行測(cè)量。因此，針對(duì)古代壁畫(huà)樣本，本文采用反射式THz-TDS系統(tǒng)對(duì)樣品進(jìn)行測(cè)量。

經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition,EMD)是20世紀(jì)90年代后期由華裔美國(guó)院士黃鄂提出的自適應(yīng)時(shí)頻分析方法[7]。相比于傳統(tǒng)的時(shí)域相關(guān)分析、時(shí)頻域小波變換等方法，EMD[8]在處理非平穩(wěn)、非線性數(shù)據(jù)上具有明顯的優(yōu)勢(shì)。

1 算法原理

1.1 EMD算法簡(jiǎn)介

EMD的核心思想是能將復(fù)雜信號(hào)分解為若干個(gè)本征模函數(shù)(Intrinsic Mode Function，IMF)之和。 黃鄂認(rèn)為，IMF分量必須滿足下面兩個(gè)條件：(1)在全部時(shí)間范圍內(nèi)，局部極值點(diǎn)數(shù)和過(guò)零點(diǎn)數(shù)必須相等，或最多相差一個(gè)；(2)在任意時(shí)刻，局部最大值的包絡(luò)與局部最小值的包絡(luò)的平均值必須為零。本文利用EMD將復(fù)雜的原始信號(hào)分解成瞬時(shí)頻率不同的若干個(gè)IMF的特性去除THz反射峰。

1.2 THz時(shí)域反射譜的EMD分解

下面采用EMD方法將復(fù)雜的原始測(cè)量信號(hào)分解成瞬時(shí)頻率不同的若干個(gè)IMF：

(1)找出原始測(cè)量信號(hào)中反射信號(hào)x(t)的所有局部最大值點(diǎn)和最小值點(diǎn)；

(2)分別得到(1)中的上包絡(luò)線xmax(t)和下包絡(luò)線xmin(t)：

xmin(t)≤x(t)≤xmax(t)

(1)

(3)計(jì)算上包絡(luò)線和下包絡(luò)線的均值：

(2)

(4)得到x(t)去掉均值的新數(shù)據(jù)序列h11(t)：

h11(t)=x(t)-m11(t)

(3)

(5)檢查h11(t)是否滿足本征模函數(shù)的條件，對(duì)于復(fù)雜信號(hào)，一般情況下，需要對(duì)h11(t)重復(fù)上述處理過(guò)程，直到滿足本征模函數(shù)的定義要求為止，則：

h1k(t)=h1(k-1)(t)-m1k(t)

(4)

(5)

令：r1(t)=x(t)-c1(t)

(6)

(6)以r1(t)為反射信號(hào)，重復(fù)上述處理步驟，則可以得到其他本征模函數(shù)：c2(t)、c3(t),...,cn(t)，即

cn(t)=rn-1(t)-rn(t)

(7)

其中rn(t)稱為余項(xiàng)。至此，原始反射信號(hào)的EMD分解結(jié)束，原信號(hào)可以表示為：

(8)

分解出的本征模函數(shù)按頻率由高到低依次排列，本文選擇最后一個(gè)本征模函數(shù)，并用它替代原始的反射峰信號(hào)，獲得處理后的THz時(shí)域光譜信號(hào)。

2 實(shí)驗(yàn)結(jié)果與討論

采用1.2節(jié)的方法處理THz時(shí)域反射峰，然后通過(guò)快速傅里葉變換(FFT)得到時(shí)域反射峰被處理后的頻域譜。

圖1和圖2分別為參考和測(cè)量信號(hào)EMD處理前后THz時(shí)域反射譜對(duì)比圖。通過(guò)比較圖1、圖2可以明顯看出，兩個(gè)信號(hào)的反射峰的處理效果良好。

圖1 處理前的THz時(shí)域反射譜

圖2 處理后的THz時(shí)域反射譜

圖3 信號(hào)時(shí)域反射峰EMD處理前的頻域譜

圖3和圖4分別為參考和測(cè)量信號(hào)EMD處理前后THz頻域反射譜對(duì)比圖。通過(guò)比較可以明顯地看出，參考信號(hào)和測(cè)量信號(hào)的時(shí)域反射峰被處理掉后，其頻域譜的振蕩被很好地抑制了，可以在頻率為0.15～1.1 THz范圍內(nèi)觀察到很明顯的3條水的吸收峰，提高了頻域信噪比，更有利于研究樣本的吸收特性。

圖4 信號(hào)時(shí)域反射峰EMD處理后的頻域譜

3 結(jié)論

針對(duì)太赫茲反射譜中時(shí)域存在反射峰與頻域出現(xiàn)振蕩的問(wèn)題，本文首次提出利用EMD自適應(yīng)分解和處理THz時(shí)域反射峰的解決方法。以古代壁畫(huà)樣本的THz反射譜為例，比較了EMD處理時(shí)域反射峰前后的時(shí)域譜，反射峰消除效果良好；通過(guò)比較經(jīng)EMD處理前后的頻域譜，發(fā)現(xiàn)頻域振蕩影響被很好地抑制，水的吸收峰明顯，提高了頻域信噪比，有利于進(jìn)一步研究樣本的吸收特性。由于反射峰的去除可增加時(shí)域有效數(shù)據(jù)，因此也提高了頻率分辨率。

[1] SIEGEL P H. Terahertz technology[J]. IEEE Transactions on Microwave Theory and Techniques, 2002, 50(3): 910-928.

[2] MASSAOUTI M, CHRISTINA D, GORODETSKY A, et al. Detection of harmful residues in honey using terahertz time-domain spectroscopy[J].Applied Spectroscopy, 2013,67(11):1264-1269.

[3]STOIK C, BOHN M, BLACKSHIRE J. Nondestructive evaluation of aircraft composites using reflective terahertz time domain spectroscopy[J]. NDT&E International 43, 2010, 43(2):106-115.

[4] VALERIA N C, MARTINA H. New insights into the role of water in biological function: studying solvated biomolecules using terahertz absorption spectroscopy in conjunction with molecular dynamics simulations [J]. Journal of the American Chemical Society, 2014, 136(37): 12800-12807.

[5] He Ting, Shen Jingling. Applications of terahertz spectroscopy in illicit drugs detection [J]. Spectroscopy and Spectral Analysis, 2013, 33(9):2348-2353.

[6] 趙國(guó)忠. 太赫茲科學(xué)技術(shù)研究的新進(jìn)展[J]. 國(guó)外電子測(cè)量技術(shù),2014,33(2):1-6,20.

[7] HUANG N E,Shen Zheng, STEVEN R L, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences . 1998,454:903-995.

[8] 寇艷廷,范濤濤,劉晨,等.EMD過(guò)程中數(shù)據(jù)擬合的算法改進(jìn)與實(shí)現(xiàn)[J].微型機(jī)與應(yīng)用,2013,32(5):66-68.

Reflective terahertz spectrum signal decompose and remove based on the empirical mode decomposition

Lu Min1， Li Xiaoxia1,2

(1. School of Information Engineering, Southwest University of Science and Technology, Mianyang 621010, China;2. Robot Technology Used for Special Environment Key Laboratory of Sichuan Province, Mianyang 621010, China)

When utilizing the reflective terahertz time domain spectroscopy system to measure the time domain signals, there will be problems like reflection peak, spectrum oscillation and low signal to noise rate in frequency domain. To address these issues, this paper proposes a method based on the empirical mode decomposition to tackle THz reflection peak in time domain as well as to eliminate the effects caused by frequency domain oscillation. We first obtain the intrinsic mode function to substitute the reflection peak, thus making time domain signal's effective length longer and THz frequency domain resolution higher in order to lessen remove the impact of frequency domain oscillation. The result of the THz reflection spectrum experiment on ancient murals proves that this approach can self-adaptively handle the reflection peak and acquire a good reparation of THz time domain signal and good inhibition of frequency domain oscillation.

terahertz reflection spectrum; empirical mode decomposition; reflection peak reparation

四川省教育廳資助科研項(xiàng)目(11ZB106)

O433.4

A

1674-7720(2016)02-0014-03

盧敏，李小霞. 基于EMD的反射式太赫茲光譜信號(hào)分解與處理[J] .微型機(jī)與應(yīng)用，2016,35(2)：14-16.

2015-08-26)

盧敏(1990-)，女，碩士，主要研究方向：太赫茲光譜信號(hào)檢測(cè)及處理。

李小霞(1976-)，通信作者，女，博士，教授，主要研究方向：模式識(shí)別、光譜信號(hào)處理。E-mail：664368504@qq.com。