EMD 包絡(luò)線擬合算法改進(jìn)及在泡沫尺寸趨勢提取中的應(yīng)用

盧 青， 王雅琳 ， 彭 凱， 薛永飛， 謝永芳

(中南大學(xué) 信息科學(xué)與工程學(xué)院，湖南 長沙410083)

經(jīng)驗(yàn)?zāi)B(tài)分解法(Empirical Mode Decomposition，EMD)在非平穩(wěn)及非線性信號的處理及分析上具有明顯優(yōu)勢，被廣泛應(yīng)用于海洋內(nèi)波探測［1］及溫度預(yù)測［2］、大氣時間序列預(yù)測［3］、天體觀測［4］、地震記錄［5］、機(jī)械故障診斷［6］等方面。包絡(luò)線擬合是EMD 算法中非常關(guān)鍵的一步。目前，常用的包絡(luò)線擬合方法是三次樣條插值法［7-10］。但現(xiàn)有三次樣條插值擬合包絡(luò)線的方法存在著插值點(diǎn)選取不當(dāng)、端點(diǎn)發(fā)散［11］，擬合松散［7］等問題，致使擬合誤差的出現(xiàn)，該擬合誤差會隨著篩分過程逐漸向內(nèi)，使得分解的整個數(shù)據(jù)被嚴(yán)重污染，導(dǎo)致分解結(jié)果失真。

目前，在包絡(luò)線擬合方面的研究主要是針對端點(diǎn)發(fā)散問題，在端點(diǎn)效應(yīng)抑制方面國內(nèi)外學(xué)者提出了很多種方法，具體可歸為3 類［12］:波形延拓法、數(shù)據(jù)預(yù)測延拓法和極值延拓法。波形匹配［13］、波形鏡像等方法屬于波形延拓法，神經(jīng)網(wǎng)絡(luò)預(yù)測［14］、ARMA 模型預(yù)測［15］、時間序列預(yù)測等方法是數(shù)據(jù)預(yù)測延拓，多項式擬合、極值對稱延拓法［11］等方法屬于極值延拓法。其中，常用的端點(diǎn)效應(yīng)抑制方法為波形匹配、神經(jīng)網(wǎng)絡(luò)預(yù)測、ARMA 模型預(yù)測、極值對稱延拓。然而波形匹配法計算量大，占用內(nèi)存;神經(jīng)網(wǎng)絡(luò)預(yù)測方法的學(xué)習(xí)過程時間較長，大大影響了EMD 的分解速度，不適合實(shí)時信號處理;ARMA 模型用于處理非線性非平穩(wěn)的信號時參數(shù)估計比較復(fù)雜，且分析效果受限;極值對稱延拓是鏡像延拓的延續(xù)與提升，方法簡單，但處理非周期信號時誤差很大;其余各種方法針對非線性非平穩(wěn)信號延拓結(jié)果因人而異，很難保證廣泛的適用性。因此需要提出一種簡單、容易實(shí)現(xiàn)、適用性強(qiáng)、有效的方法。

文中通過加權(quán)最小二乘多項式擬合方法進(jìn)行極值延拓，根據(jù)延拓的實(shí)際情況對得到的端點(diǎn)值進(jìn)行合理性判斷并加以取舍，然后作為插值擬合曲線的端點(diǎn);同時結(jié)合包絡(luò)線擬合中存在的另外兩個問題，提出了EMD 中包絡(luò)線擬合的改進(jìn)方法。為評價改進(jìn)方法的包絡(luò)線擬合效果，文中定義了包絡(luò)線與信號的方差和以及曲線長度協(xié)同作為評判標(biāo)準(zhǔn)。

鋁土礦泡沫浮選中不同槽間的泡沫尺寸變化趨勢存在關(guān)聯(lián)規(guī)則，要挖掘相關(guān)規(guī)則，需先提取每個槽的泡沫尺寸的變化趨勢，從而獲取泡沫尺寸變化的趨勢圖，為關(guān)聯(lián)性分析打下基礎(chǔ)。文中將所提包絡(luò)線擬合改進(jìn)方法應(yīng)用于粗選槽泡沫尺寸的趨勢提取，獲得了較好的粗選槽泡沫尺寸變化趨勢，驗(yàn)證了所提方法的有效性。

1 EMD 算法及其存在的問題

1.1 EMD 算法簡介

美籍華人HUANG 等于1998 年提出了一種適合于分析非線性、非平穩(wěn)信號序列的經(jīng)驗(yàn)?zāi)B(tài)分解法EMD。該方法可以將復(fù)雜信號分解為有限數(shù)量的含原信號不同時間尺度的局部特征信號的本征模函數(shù)(Intrinsic Mode Function，IMF)［3］，進(jìn)而根據(jù)不同頻率分量的特性進(jìn)行濾波或預(yù)測。

EMD 算法步驟如下:

1)尋找信號X(t)的所有極值點(diǎn)。

2)分別以所有極大值點(diǎn)、極小值點(diǎn)為插值點(diǎn)進(jìn)行三次樣條插值，得到信號X(t)的上、下包絡(luò)線。

3)求出上下包絡(luò)線的平均值m(t)，用以表示信號的平均值曲線。

4)X(t)－m(t)= h1，完成一次篩分過程。第2次篩分過程用h1代替X(t)作為分解信號，重復(fù)k次，直到h1k滿足IMF 條件:①函數(shù)在整個時間范圍內(nèi)，局部極值點(diǎn)和過零點(diǎn)的數(shù)目必須相等，或最多相差一個;②在任意時刻，局部最大值的包絡(luò)(上包絡(luò)線)和局部最小值的包絡(luò)(下包絡(luò)線)平均為零。

5)記c1= h1k，做運(yùn)算X(t)－ c1= r1，再將r1作為分解數(shù)據(jù)，重復(fù)篩分過程，得到c2，r2，…，cn，rn，直到rn或cn小于給定值或rn成為單調(diào)函數(shù)不可再篩分出IMF。則分解完成，有

其中，ci為第i 個IMF 分量，代表了原始信號X(t)中不同特征時間尺度的信號分量;rn為殘余函數(shù)，代表信號的平均趨勢。

1.2 存在問題及原因

EMD 的關(guān)鍵步驟之一是用三次樣條函數(shù)擬合得到信號的上、下包絡(luò)，它將直接影響分解得到的所有IMF 分量的精度。然而在進(jìn)行信號的包絡(luò)線擬合時，由于信號在端點(diǎn)位置往往不是極值點(diǎn)，使得擬合的包絡(luò)線在邊界附近出現(xiàn)嚴(yán)重失真，這種現(xiàn)象被稱為邊界效應(yīng)或端點(diǎn)效應(yīng)。實(shí)際應(yīng)用中信號中存在的噪聲及偏差不可能完全去除，可能導(dǎo)致極值點(diǎn)數(shù)量增加;同時由于邊界效應(yīng)產(chǎn)生的誤差及極值點(diǎn)數(shù)量的增加將引起擬合誤差，且隨著篩分過程的不斷反復(fù)逐漸傳播到信號的內(nèi)部并污染整個信號，最終導(dǎo)致得到不可靠的IMF 分量和殘余分量。

例如，分段信號

圖1 是用三次樣條插值方法在時間段［0，20］內(nèi)擬合信噪比為70 的信號y0(t)的波形圖。

圖1 三次樣條插值法擬合y0(t)的包絡(luò)線Fig.1 Envelope extraction of y0(t)using CSI

包絡(luò)線應(yīng)當(dāng)完整地包絡(luò)信號，且是光滑的曲線，上包絡(luò)線的所有值不小于對應(yīng)點(diǎn)的信號值，下包絡(luò)線的所有值不大于對應(yīng)點(diǎn)的信號值。由圖1 可以看出，三次樣條插值法對于非平穩(wěn)信號的擬合效果不佳。究其原因:

1)所采用的信號不是平滑信號。實(shí)際應(yīng)用中由于采集設(shè)備誤差、記錄誤差及現(xiàn)場環(huán)境影響等信號數(shù)據(jù)本身就不完全平滑，且信號采集過程中無法避免噪聲信號的干擾，濾波和去噪在一定程度上可以凈化信號，但提取包絡(luò)線時信號中必然還殘留一些噪聲，容易出現(xiàn)極值點(diǎn)冗余現(xiàn)象，若把這些極值點(diǎn)全部作為插值點(diǎn)，會出現(xiàn)包絡(luò)線貼合于原信號曲線的現(xiàn)象，導(dǎo)致曲線擬合失真。

2)實(shí)際應(yīng)用中采集到的信號可能會出現(xiàn)某一小段時間內(nèi)幅值不變的情況，由極值點(diǎn)的定義可知這一時間段內(nèi)沒有極值點(diǎn)，若依舊只是把極值點(diǎn)作為插值點(diǎn)就會出現(xiàn)一段包絡(luò)線的空缺。

3)端點(diǎn)無約束，包絡(luò)線在端點(diǎn)處出現(xiàn)發(fā)散現(xiàn)象。

4)相鄰的兩個插值點(diǎn)之間的插值曲線斜率出現(xiàn)過零點(diǎn)，即擬合曲線不單調(diào)，部分包絡(luò)值在臨近極值之外，出現(xiàn)了“凸峰”或“凹谷”現(xiàn)象，導(dǎo)致包絡(luò)線擬合的不“緊”。

2 算法改進(jìn)及效果評判

針對1.2 中提到的包絡(luò)線擬合中插值點(diǎn)的選取、端點(diǎn)效應(yīng)的抑制以及插值曲線擬合松散的問題，經(jīng)綜合考慮，提出有效的解決辦法對EMD 中包絡(luò)線的擬合予以改進(jìn)。

以同時滿足包絡(luò)線與信號的方差和盡可能小且包絡(luò)線長度盡可能小，作為評判包絡(luò)線擬合效果的標(biāo)準(zhǔn)。包絡(luò)線與信號的方差和小說明包絡(luò)線與信號曲線貼合度好，無過沖現(xiàn)象及端點(diǎn)發(fā)散現(xiàn)象，但包絡(luò)線與信號的方差和太小又說明包絡(luò)線和信號有重疊趨勢，包絡(luò)線擬合效果不好，因此以包絡(luò)線與信號的方差和單獨(dú)作為評價標(biāo)準(zhǔn)有一定合理性卻并不科學(xué)。包絡(luò)線長度值小說明無過沖、欠沖、端點(diǎn)發(fā)散現(xiàn)象，但包絡(luò)線長度最短，即以直線連接各插值點(diǎn)時，不滿足包絡(luò)線光滑的要求，因此包絡(luò)線長度最短單獨(dú)作為評價標(biāo)準(zhǔn)也有一定合理性卻并不科學(xué)。所以，以兩項指標(biāo)加權(quán)和e 作為評判包絡(luò)線擬合效果的標(biāo)準(zhǔn)。權(quán)值選擇的原則是使兩個標(biāo)準(zhǔn)計算值的數(shù)量級相同，即es+ex與| ys| +| yx| 的數(shù)量級相同，選取辦法是將數(shù)量級高的一個標(biāo)準(zhǔn)的數(shù)值降級。

記上包絡(luò)線為ys(t)，下包絡(luò)線為yx(t)，上、下包絡(luò)線與信號的方差和分別為es和ex，上下包絡(luò)線的長度分別為| ys| 和| yx|，則由方差和及曲線長度的定義可得

式中:i 為第i 個上包絡(luò)線插值點(diǎn);p 為上包絡(luò)線插值點(diǎn)個數(shù);ysi為第i 與i + 1 個插值點(diǎn)間的插值函數(shù);y'si為第i 與i +1 個插值點(diǎn)間插值函數(shù)的一階導(dǎo)數(shù)，i = 0 時取上包絡(luò)線左端點(diǎn)，i = p 時i +1 = p +1 取上包絡(luò)線右端點(diǎn);j 為第j 個下包絡(luò)線插值點(diǎn);q 為下包絡(luò)線插值點(diǎn)個數(shù);yxj為第j 與j +1 個插值點(diǎn)間的插值函數(shù);y'xj為第j 與j +1 個插值點(diǎn)間插值函數(shù)的一階導(dǎo)數(shù)，j = 0 時取下包絡(luò)線左端點(diǎn)，j = q 時j +1 =q +1 取下包絡(luò)線的右端點(diǎn)。

對于y0(t)，鑒于包絡(luò)線與信號的方差和是包絡(luò)線長度數(shù)量級的100 倍，所以取

圖1 中三次樣條插值法直接擬合曲線y0(t)的誤差見表1。

表1 三次樣條插值法擬合包絡(luò)線的誤差Tab.1 Error of the envelope extraction using CSI

2.1 基于特殊點(diǎn)的三次樣條插值

首先，針對小段值不變的平穩(wěn)數(shù)據(jù)導(dǎo)致的無極值點(diǎn)問題，提出一種基于特殊點(diǎn)的三次樣條插值方法。特殊點(diǎn)既包括傳統(tǒng)意義上的極值點(diǎn)又包括平穩(wěn)數(shù)據(jù)段的所有點(diǎn)。具體選取方法為:上包絡(luò)線插值點(diǎn)x(ti)滿足

下包絡(luò)線插值點(diǎn)x(tj)滿足

其中，M，N 分別為上、下包絡(luò)線插值點(diǎn)的個數(shù)。

其次，針對信號誤差引起的小段時間內(nèi)極值點(diǎn)數(shù)量增多且極值波動很小導(dǎo)致擬合失真的問題，選擇插值點(diǎn)時，在尋找到特殊點(diǎn)后增加一步篩選操作。以下包絡(luò)線的擬合為例闡述篩選原理:將所有極小值點(diǎn)組成新的時間序列，按照一定的時間間隔Δ 對x1(t)進(jìn)行分段，對每一段的所有時間數(shù)據(jù)求均值，該時間均值處的點(diǎn)作為插值點(diǎn)進(jìn)行包絡(luò)線擬合。具體步驟如下:

1)由時間序列的時間特性可知，后一時刻的時間值大于前一時刻的時間值，即ti+1＞ti＞ti－1。將時間序列x1(t)在ti－1與ti之間分段，分界點(diǎn)產(chǎn)生規(guī)則是:①ti與前一點(diǎn)ti－1的時間間隔大于一定的閾值σ，即Δ = ti－ti－1＞σ，σ 由時間序列本身的波動周期(周期序列)及序列長度決定。②若x1(ti)＞x1(ti+1)，則必須成立;否則，若x1(ti)＜x1(ti+1)，則必須成立。閾值R由序列本身決定。

2)對每一段數(shù)據(jù)的時間坐標(biāo)求算數(shù)平均ˉt，用該均值處的點(diǎn)(ˉt，x1(ˉt))近似作為該段的唯一插值點(diǎn)。改進(jìn)后基于特殊點(diǎn)的包絡(luò)線擬合如圖2 所示，擬合誤差見表2。

圖2 基于特殊點(diǎn)的包絡(luò)線擬合Fig.2 Envelope extraction based on the special points

表2 基于特殊點(diǎn)的包絡(luò)線擬合誤差Tab.2 Error of the envelope extraction based on the special points

通過圖2 及表2 與圖1 及表1 的對比可發(fā)現(xiàn)，3個信號擬合包絡(luò)線的長度及包絡(luò)線與信號的方差和都有明顯減小，誤差e 均大幅減小，從圖2 中也可以看出擬合效果有明顯改善。

2.2 端點(diǎn)效應(yīng)的抑制

端點(diǎn)發(fā)散是因?yàn)槎它c(diǎn)無約束，所以解決辦法是增加約束。若獲取的信號量足夠多，采樣時在有用信號的前后分別多采集一段數(shù)據(jù)，以增加極值點(diǎn)的個數(shù)作為約束，這樣得到的包絡(luò)線在端點(diǎn)處不會嚴(yán)重扭曲。若已知曲線方程，可以在信號前后再計算出幾個極值點(diǎn)作為插值點(diǎn)進(jìn)行包絡(luò)線擬合。但是在沒有多余數(shù)據(jù)又不知道曲線方程的前提下，在擬合包絡(luò)線時，需要對已知插值數(shù)據(jù)進(jìn)行延拓，以在數(shù)據(jù)兩端獲得附加的極大值和極小值作為約束抑制端點(diǎn)效應(yīng)。

加權(quán)最小二乘多項式擬合是對所有的點(diǎn)進(jìn)行整體擬合，而三次樣條插值是在相鄰兩個插值點(diǎn)間建立三次多項式進(jìn)行擬合的分段擬合，端點(diǎn)處為自由端點(diǎn)，所以加權(quán)最小二乘多項式擬合得到的端點(diǎn)比自由端點(diǎn)精確得多。為此，提出加權(quán)最小二乘多項式擬合實(shí)現(xiàn)EMD 包絡(luò)線擬合中端點(diǎn)發(fā)散的抑制。

以下包絡(luò)線左端點(diǎn)的確定為例說明采用加權(quán)最小二乘多項式擬合確定約束條件的方法。加權(quán)最小二乘原理是找到一個函數(shù)

可以最佳地擬合于諸觀測點(diǎn)

即應(yīng)當(dāng)滿足

此處選取擬合多項式的最高次數(shù)為3 次，盡管擬合次數(shù)越高，擬合精度越高，擬合效果越好，但4次擬合計算量急劇增大，且易出現(xiàn)放大微小誤差的問題。且3 次擬合與2 次擬合相比，計算量相差不多，擬合精度卻會有所提升。相應(yīng)地，插值點(diǎn)至少選取4 個。但當(dāng)特殊點(diǎn)數(shù)量不足4 個時，選取全部的特殊點(diǎn)，并相應(yīng)地可降低擬合多項式的次數(shù)。此時

需滿足

權(quán)值選取原則是誤差項方差越大則權(quán)值越小，權(quán)值的最佳取值是方差和的倒數(shù)。但是由于方差和不可知，習(xí)慣選取方法是取橫坐標(biāo)x 的倒數(shù)。此處，根據(jù)實(shí)際情況，x 應(yīng)理解為觀測點(diǎn)橫坐標(biāo)xi與左端點(diǎn)橫坐標(biāo)x0的距離，即

若增加的下包絡(luò)線端點(diǎn)值大于臨近的極大值點(diǎn)，顯然不符合實(shí)際，舍棄該端點(diǎn)不用，以靠近端點(diǎn)的第一個下包絡(luò)特殊點(diǎn)和第一個上包絡(luò)特殊點(diǎn)的均值作為端點(diǎn)值。雖然這樣取得的端點(diǎn)值存在誤差，但由于本意并不是獲得精確的端點(diǎn)值，只要可以有效抑制端點(diǎn)發(fā)散問題就是可取的。同理，對上包絡(luò)線添加端點(diǎn)，如圖3 所示的繪出包絡(luò)線，擬合誤差見表3。

圖3 抑制端點(diǎn)效應(yīng)后的包絡(luò)線Fig.3 Extraction with the restraining end effect

表3 抑制端點(diǎn)效應(yīng)后的擬合誤差Tab.3 Error of extraction with the restraining end effect

由圖3 及表3 與圖2 及表2 相比較可以看出，發(fā)現(xiàn)擬合誤差e 減小，曲線擬合效果有一定的改善。

2.3 保形分段三次Hermite 插值

由圖3 可以看出，在抑制端點(diǎn)效應(yīng)之后三次樣條插值法擬合出的包絡(luò)線依然存在與信號曲線有交點(diǎn)(下包絡(luò)線部分點(diǎn)取值大于信號值或上包絡(luò)線部分取值小于信號值)，且相鄰插值點(diǎn)間的下包絡(luò)線有過沖現(xiàn)象，擬合的包絡(luò)線松散。原因是相鄰兩插值點(diǎn)間存在單調(diào)性不一致，包絡(luò)線缺少約束。

針對上述問題，首先要增加一組對包絡(luò)線擬合的約束條件:上包絡(luò)線值始終不小于相同坐標(biāo)處的信號值;同時下包絡(luò)線的所有值都不大于相同坐標(biāo)處的信號值。即對于所有時間點(diǎn)t，有ys(t)≥y(t)≥yx(t)。然后采用分段三次埃爾米特多項式插值［16］(Cubic HermiteInterpolation，CHI)進(jìn)行包絡(luò)線擬合，可以有效改善這一缺陷。

在區(qū)間［a，b］上構(gòu)造三次埃爾米特插值時，首先取a = x0＜x1＜x2＜… ＜xn= b，記xi處函數(shù)值yi，導(dǎo)數(shù)值y'i;然后在每個子區(qū)間［xi，xi+1］，(i =0，1，…，n － 1)上構(gòu)造三次插值多項式Hi(x)，Hi(x)需滿足Hi(xi)= yi，H'i(xi)= y'i。但一階導(dǎo)數(shù)y'i通常不可知，所以y'i的確定就成為CHI 的關(guān)鍵。

具體步驟為:

1)借鑒Lagrange 插值基函數(shù)法，引入4 個三次多項式α0(x)，α1(x)，β0(x)，β1(x)作為基函數(shù)，令Hi(x) = yiα0(x) + yi+1α1(x) + y'iβ0(x) +y'i+1β1(x)，并求解基函數(shù)。

2)于是，x ∈［xi，xi+1］時，兩點(diǎn)Hermite 插值

所以，［a，b］上的插值函數(shù)為

3)一階導(dǎo)數(shù)的確定:(xi，yi)(i = 1，2，…，n －1)前后子區(qū)間(xi－1，xi)，(xi，xi+1)的斜率分別為

所以(xi，yi)點(diǎn)處斜率di為

其中

CHI 插值具有在每個子區(qū)間單調(diào)的特性，且公式簡單、計算量小、穩(wěn)定性好［17］。基于CHI 插值的包絡(luò)線擬合如圖4 所示，擬合誤差見表4。

圖4 CHI 擬合包絡(luò)線Fig.4 Envelope extraction by using CHI

表4 CHI 擬合后的誤差Tab.4 Error of extraction using CHI

與圖3 及表3 相比，誤差e 減小，同時可以看出包絡(luò)線在兩插值點(diǎn)間擬合曲線的過沖、欠沖，即擬合松散問題得到改善。由此證明了文中所提出的改進(jìn)方法有效。

2.4 仿真驗(yàn)證

圖5(a)、(b)分別為信噪比為30 的非線性非平穩(wěn)信號

直接用三次樣條插值法擬合及改進(jìn)方法擬合的包絡(luò)線，擬合誤差見表5。

圖6(a)、(b)分別為信噪比為30 的信號y2(t)= sin(8πt)+5sin(0.4πt)+2sin(2πt +π/3)，直接用三次樣條插值法擬合及改進(jìn)方法擬合的包絡(luò)線，擬合誤差見表6。

圖5 y1(t)擬合包絡(luò)線Fig.5 Envelope extraction of y1(t)

表5 y1(t)擬合誤差對比Tab.5 Comparison of extraction error of y1(t)

圖6 y2(t)擬合包絡(luò)線Fig.6 Envelope extraction of y2(t)

表5 和表6 均取

通過上述兩例的對比可以發(fā)現(xiàn)，e 減小，同時可以看出改進(jìn)后的包絡(luò)線擬合方法擬合出的包絡(luò)線包的更緊密，端點(diǎn)發(fā)散現(xiàn)象得到抑制。由此驗(yàn)證了所提改進(jìn)方法在包絡(luò)線擬合上的有效性。

表6 y2(t)擬合誤差對比Tab.6 Comparison of extraction error of y2(t)

3 改進(jìn)的包絡(luò)線提取方法在泡沫尺寸變化趨勢提取上的應(yīng)用

典型的鋁土礦泡沫浮選操作流程，包括粗選、精選、粗掃和精掃作業(yè)。實(shí)際生產(chǎn)中，操作工人通過觀察泡沫進(jìn)行生產(chǎn)操作。工業(yè)上通過分布機(jī)器視覺系統(tǒng)在每個工序獲取泡沫圖像并從中提取出顏色、尺寸、紋理等信息。鑒于各工序泡沫存在關(guān)聯(lián)信息，對不同工序之間的泡沫信息進(jìn)行關(guān)聯(lián)性分析得到關(guān)聯(lián)規(guī)則，可以用于指導(dǎo)生產(chǎn)，對泡沫尺寸的包絡(luò)線分析可以獲取泡沫尺寸變化的規(guī)律，便于進(jìn)行關(guān)聯(lián)性分析。

選取基于分布機(jī)器視覺的鋁土礦泡沫浮選過程在粗選槽獲取的圖像中提取的一班泡沫尺寸數(shù)據(jù)，時間間隔為5 min，剔除異常值后作為原始數(shù)據(jù)。針對這組數(shù)據(jù)，直接用三次樣條插值法擬合包絡(luò)線與用改進(jìn)后的方法擬合的包絡(luò)線如圖7 所示，擬合誤差見表7。

對比圖7(a)和圖7(b)及表7 中所列擬合誤差，可以發(fā)現(xiàn)包絡(luò)線與信號的方差和以及曲線長度同步減小，誤差e 減小;同時可以看出改進(jìn)后的包絡(luò)線擬合方法擬合出的包絡(luò)線包的更緊密，不存在插值點(diǎn)冗余及端點(diǎn)發(fā)散問題。再次驗(yàn)證了所提改進(jìn)方法在包絡(luò)線擬合上的有效性。

表7 包絡(luò)線擬合方法改進(jìn)前后的擬合誤差Tab.7 Extraction error before and after improvement

4 結(jié) 語

文中分析了EMD 算法中基于三次樣條插值的包絡(luò)線擬合方法，列出了存在的問題并分析了原因，提出了基于特殊點(diǎn)的EMD 包絡(luò)線擬合改進(jìn)算法。在極值點(diǎn)的基礎(chǔ)上增加了轉(zhuǎn)折點(diǎn)與值不變的平穩(wěn)點(diǎn)組成特殊點(diǎn)，經(jīng)過分段篩選刪除冗余插值點(diǎn)。針對端點(diǎn)發(fā)散問題，通過加權(quán)最小二乘多項式擬合對所有極值點(diǎn)進(jìn)行多項式擬合計算出端點(diǎn)值。三次埃爾米特插值基于其在每個子區(qū)間單調(diào)的特性及對上下包絡(luò)線數(shù)值的約束可以有效改善三次樣條插值包絡(luò)線擬合中存在的包絡(luò)線松散且與信號有交叉的情況。

在有針對性的提出相應(yīng)改進(jìn)方法的同時，通過定義的擬合優(yōu)良性評判標(biāo)準(zhǔn)，及逐步繪出改進(jìn)后的包絡(luò)圖，構(gòu)成清晰明確的客觀數(shù)據(jù)及主觀視覺的對比，表明了所提改進(jìn)方法在包絡(luò)線擬合上的良好效果。最后將改進(jìn)方法應(yīng)用于鋁土礦泡沫浮選粗選槽獲取的泡沫尺寸數(shù)據(jù)，提取其變化趨勢，再次證明了所提方法對包絡(luò)線擬合的有效性。

［1］宋海斌，拜陽，董崇志，等. 南海東北部內(nèi)波特征——經(jīng)驗(yàn)?zāi)B(tài)分解方法應(yīng)用初探［J］. 地球物理學(xué)報，2010，53(2):393-400.

SONG Haibin，BAI Yang，DONG Chongzhi，et al.A preliminary study of application of empirical mode decomposition method in understanding the features of internal waves in the northeastern south china sea［J］.Chinese Journal of Geophysics，2010，53(2):393-400.(in Chinese)

［2］Bhuiyan S M A，Carey S，Khan J F.Trend analysis of sea surface temperature and sea level pressure employing a pseudo-EMD method［C］//2013 Proceedings of IEEE Southeastcon.Jacksonville，F(xiàn)L:IEEE，2013:1-4.

［3］玄兆燕，楊公訓(xùn).經(jīng)驗(yàn)?zāi)B(tài)分解法在大氣時間序列預(yù)測中的應(yīng)用［J］.自動化學(xué)報，2008，34(1):97-101.

XUAN Zhaoyan，YANG Gongxun.Application of EMD in the atmosphere time series prediction［J］.Acta Automatic Sinica，2008，34(1):97-101.(in Chinese)

［4］XU T，WU J，WU Z S，et al.Long-term sunspot number prediction based on EMD analysis and AR model［J］.Chinese Journal of Astronomy and Astrophysics，2008，8(3):337-342.

［5］錢昌松，劉代志，劉志剛，等. 基于遞歸高通濾波的經(jīng)驗(yàn)?zāi)B(tài)分解及其在地震信號分析中的應(yīng)用［J］. 地球物理學(xué)報，2010，53(5):1215-1225.

QIAN Changsong，LIU Daizhi，LIU Zhigang，et al.EMD based on recursive high-pass filter and its application on seismic signal analysis［J］.Chinese Journal of Geophysics，2010，53(5):1215-1216，1224-1225.(in Chinese)

［6］程軍圣，于德介，楊宇.基于內(nèi)稟模態(tài)奇異值分解和支持向量機(jī)的故障診斷方法［J］.自動化學(xué)報，2006，32(3):475-480.

CHENG Junsheng，YU Dejie，YANG Yu. Fault diagnosis approach based on intrinsic mode singular valuedecomposition and support vector machines［J］.Acta Automatic Sinica，2006，32(3):475-480.(in Chinese)

［7］朱偉芳，趙鶴鳴，陳小平.一種最小長度約束的EMD 包絡(luò)擬合方法［J］.電子學(xué)報，2012，40(9):1909-1912.

ZHU Weifang，ZHAO Heming，CHEN Xiaoping. Aleast length constrained envelope approach for EMD［J］. Acta Electronica Sinica，2012，40(9):1909-1912.(in Chinese)

［8］張勇，陳濱.黃變換的穿越篩分法［J］.電子科技大學(xué)學(xué)報，2007，36(1):8-10.

ZHANG Yong，CHEN Bin.Crossing sifting method based on huangtransform［J］.Journal of University of Electronic Science and Technology of China，2007，36(1):8-10.(in Chinese)

［9］JIA X F，AN H Q，ZHANG S G.Cubic spline interpolation method for theenvelope tracking of middle and low frequency voltage flicker［J］.Thermal Power and Electrical Engineering III，2014(960/961):704-709.

［10］QIU Mianhao，LIU Jing，CONG Hua.Research of direct sifting EMD based on cubic B-spline interpolation curve and its application in processing mechanical vibrating signals［J］.Journal of Academy of Armored Force Engineering，2007，21(3):29-33.

［11］GUO Yanwei，ZHOU Tao. The new method based on improvement EMD end effect for power quality analysis［J］. Electrical Engineering，2013，12:13-16.

［12］曹端超，康建設(shè)，趙勁松，等.EMD 端點(diǎn)效應(yīng)抑制方法仿真比較與實(shí)例分析研究［J］.機(jī)械傳動，2013，37(3):83-87.

CAO Duanchao，KANG Jianshe，ZHAO Jingsong，et al. Research on the comparison with methods of restraining endingeffect of EMD and its application in fault diagnosis［J］. Mechanical Transmission，2013，37(3):83-87.(in Chinese)

［13］李合龍，王龍，李明建，等.基于平均包絡(luò)線匹配算法的EMD 端點(diǎn)效應(yīng)分析及在股價趨勢分解中的應(yīng)用［J］.系統(tǒng)工程理論與實(shí)踐，2013，33(8):2072-2079.

LI Helong，WANG Long，LI Mingjian，et al.The end effect of EMD based on matching mean envelope and itsapplications in trend decomposition of stock price［J］.Systems Engineering-Theory and Practice，2013，33(8):2072-2079.(in Chinese)

［14］孟宗，顧海燕，李姍姍.基于神經(jīng)網(wǎng)絡(luò)集成的B 樣條經(jīng)驗(yàn)?zāi)B(tài)分解端點(diǎn)效應(yīng)抑制方法［J］.機(jī)械工程學(xué)報，2013，49(9):106-112.

MENG Zong，GU Haiyan，LI Shanshan. Restraining method for end effect of B-spline empirical modedecomposition based on neural network ensemble［J］.Journal of Mechanical Engineering，2013，49(9):106-112.(in Chinese)

［15］程軍圣，于德介，楊宇.Hilbert-Huang 變換端點(diǎn)效應(yīng)問題的探討［J］.振動與沖擊，2005，24(6):40-42，47.

CHENG Junsheng，YU Dejie，YANG Yu.Discussion of the end effects in Hilbert-Huang transform［J］.Journal of Vibration and Shock，2005，24(6):40-42，47.(in Chinese)

［16］鄧林峰，趙榮珍.基于CHI-LMD 方法的轉(zhuǎn)子振動信號分析［J］.振動與沖擊，2014，35(15):58-64.

DENG Linfeng，ZHAO Rongzhen. Rotor vibration signal analysis based on the CHI-LMD method［J］. Journal of Vibration and Shock，2014，35(15):58-64.(in Chinese)

［17］ZHU W，ZHAO H，CHEN X. Improving empirical mode decomposition with an optimized piecewise cubic hermiteinterpolation method［C］//2012 International Conference on Systems and Informatics (ICSAI). Yantai:IEEE，2012:1698-1701.