日韩欧美在,日本在线免费视频,中文字幕夜夜嗨

您現(xiàn)在的位置：特種加工機(jī)床網(wǎng)> 技術(shù)前沿>基于EMD分解在電火花數(shù)據(jù)分析的應(yīng)用

基于EMD分解在電火花數(shù)據(jù)分析的應(yīng)用

2016-12-13 來(lái)源：合肥工業(yè)大學(xué)電氣與自動(dòng)化工程學(xué) 作者：劉春杜雲(yún)

摘要：經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode dccomposition，EMD)是處理非平穩(wěn)、非線(xiàn)性信號(hào)一種有效的新方法。運(yùn)用EMD分解法將脈沖典型負(fù)荷一電火花信號(hào)中本征模態(tài)函數(shù)(intrinsic mode function，IMF)分量逐級(jí)分離出來(lái)，得到原信號(hào)的多尺度振蕩特性；分析各個(gè)分量與原序列的顯著性檢驗(yàn)和相關(guān)系數(shù)以及各個(gè)分量自身的周期性。用IMF分量和趨勢(shì)項(xiàng)合成原序列并作誤差分析，得到信號(hào)的總誤差率通過(guò)理論計(jì)算為4．3％。利用bior4．4和db2兩種小波基在3層分解系數(shù)基礎(chǔ)上對(duì)電火花數(shù)據(jù)展開(kāi)小波包變換。最后借助MATLAB平臺(tái)對(duì)EMD和小波包分解方法進(jìn)行仿真并做理論對(duì)比分析。結(jié)果表明，EMD分解法在提取機(jī)床運(yùn)行放電特性及不同工況下信號(hào)更具優(yōu)越性和有效性，提取結(jié)果更能滿(mǎn)足實(shí)際工程應(yīng)用需求。

關(guān)鍵詞：脈沖電能；電火花；經(jīng)驗(yàn)?zāi)B(tài)分解；小波包變換

1.引言

近幾年，采用電火花加工技術(shù)在加工復(fù)雜型面、深槽和窄縫等應(yīng)用中所占比重越來(lái)越大。脈沖電能負(fù)荷不同的工況應(yīng)用條件，造成脈沖電能質(zhì)量擾動(dòng)因素種類(lèi)繁多。放電狀態(tài)又是對(duì)電火花加工過(guò)程實(shí)施準(zhǔn)確控制的主要依據(jù)，直接決定了加工質(zhì)量、效率和電極損耗等工藝指標(biāo)的優(yōu)劣。然而，傳統(tǒng)電火花加工工藝效果受多種工藝參數(shù)(尤其是電參數(shù))影響，加工工藝效果難以實(shí)現(xiàn)自動(dòng)化生產(chǎn)。因此，研究電火花加工電參數(shù)的優(yōu)化選擇，對(duì)工藝效果控制及自動(dòng)化加工的實(shí)現(xiàn)具有重要意義¨。。從而，進(jìn)一步提高加工效率、精度及加工過(guò)程的穩(wěn)定性。

為了更好地提高電火花加工表面質(zhì)量，國(guó)內(nèi)外學(xué)者對(duì)其加工表面微觀形貌做了大量的研究。研究比較多的有高斯濾波和小波濾波。Raja等人”1將快速傅立葉變換引入高斯濾波。Hara等人H。用兩階Butterwonh高通濾波器模擬高斯濾波的轉(zhuǎn)換特性，提高了運(yùn)算速度卻導(dǎo)致了相移，降低了運(yùn)算精度。但是，傅里葉變換由于信號(hào)的復(fù)雜性導(dǎo)致頻率變化劇烈，且它的窗口的大小固定不變，導(dǎo)致分辨率固定。針對(duì)該缺陷，陳慶虎等人∞1選擇Shannon小波，采用Mallat快速分解算法，提出基于小波分析的理論和方法，小波濾波產(chǎn)生的三維基準(zhǔn)

面，光滑自然，能精確地把表面粗糙度和其他表面結(jié)構(gòu)成份分離。小波變換繼承了窗口變換的局部化思想并改進(jìn)了傅里葉變換的弊端，是一種白適應(yīng)的時(shí)頻分析方法；但是小波變換在小波基選擇、小波基函數(shù)固定、閥值確定對(duì)于不同信號(hào)很難做出最佳選擇。為徹底解決傳統(tǒng)分析方法對(duì)電火花信號(hào)的的局限性∞{j，將不需要信號(hào)先驗(yàn)知識(shí)和無(wú)需定義基函數(shù)的經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)法引入到電火花信號(hào)分析中，把電火花信號(hào)內(nèi)在的放電特性通過(guò)本征模態(tài)函數(shù)(IMF)分量展現(xiàn)出來(lái)。

本文運(yùn)用EMD分解法對(duì)機(jī)床放電特性的電火花數(shù)據(jù)進(jìn)行逐級(jí)分解，并綜合已有的小波分析方法作對(duì)比，重點(diǎn)分析了IMF分量與原數(shù)據(jù)顯著性檢驗(yàn)和相關(guān)系數(shù)，采用自相關(guān)函數(shù)求功率譜尋找每個(gè)分量周期值；將IMF分量和趨勢(shì)分量逆變換疊加還原數(shù)據(jù)，對(duì)比原數(shù)據(jù)做誤差分析并考量原序列的主成分。本文最后部分，將實(shí)測(cè)電火花的電流信號(hào)分段做EMD生產(chǎn)實(shí)際分析。EMD分解法是分析電火花機(jī)床的一種有效方法，IMF分量可體現(xiàn)電火花機(jī)床的內(nèi)在多尺度振蕩特性～放電特性，為電火花機(jī)床的運(yùn)行狀態(tài)進(jìn)行實(shí)時(shí)監(jiān)控提供指導(dǎo)。

2.理論基礎(chǔ)

2.1 IMF

EMD將原序列分解成有限個(gè)IMF分量，并且分量的瞬時(shí)頻率各不相同；使得IMF分量有物理意義的必要定義是：函數(shù)關(guān)于零均值線(xiàn)是局部對(duì)稱(chēng)，并且穿零點(diǎn)和極值點(diǎn)的數(shù)目是一樣的。鑒于此，給出本征模態(tài)函數(shù)的確切定義：1)在信號(hào)序列里，極值點(diǎn)的數(shù)目與穿零點(diǎn)的數(shù)目必須相等或者最多相差一個(gè)；2)由局部極大值所構(gòu)成的包絡(luò)線(xiàn)以及由局部極小值所構(gòu)成的包絡(luò)線(xiàn)的平均值為零。對(duì)于非平穩(wěn)信號(hào)而言，定義一要對(duì)“局部平均”一段一段求，然而局部時(shí)間長(zhǎng)度很難知道；定義二是一個(gè)新思路，用上下包絡(luò)線(xiàn)的局部平均值為零來(lái)表征局部對(duì)稱(chēng)性，巧妙的避開(kāi)了局部時(shí)間尺度的問(wèn)題。

zMF表征了信號(hào)序列內(nèi)在的，固有的振蕩狀態(tài)；“狀態(tài)”一詞，泛指漲落、起伏、變化，不一定具有周期性；但是如果一個(gè)振蕩經(jīng)過(guò)顯著性檢驗(yàn)具有周期性，則稱(chēng)之為周期性振蕩。

2．2 經(jīng)驗(yàn)?zāi)B(tài)函數(shù)分量mF的篩選過(guò)程

EMD的關(guān)鍵性問(wèn)題就是如何把一個(gè)非平穩(wěn)非線(xiàn)性的信號(hào)分解成有限個(gè)IMF分量和一個(gè)趨勢(shì)項(xiàng)。采用合適的篩選方法，篩選出符合定義的IMF分量。給定原信號(hào)，把所有局部極大值用3階樣條曲線(xiàn)連接起來(lái)得到上包絡(luò)線(xiàn)，運(yùn)用同樣方法把局部極小值連接起來(lái)得到下包絡(luò)線(xiàn)；上下包絡(luò)線(xiàn)分別涵蓋信號(hào)全部極大值和極小值；再取上下包絡(luò)線(xiàn)的均值線(xiàn)m，，將原序列x(￡)與均值線(xiàn)m。(￡)做差得到：

條連接后把某些極大值夸大了，出現(xiàn)了一些過(guò)高峰和過(guò)低谷；2)信號(hào)復(fù)雜，在峰或谷上有很多騎行波，就是零均值線(xiàn)的上下方起伏的彎曲，騎行波不滿(mǎn)足IMF要求且一次篩選不能完全消除騎行波；3)非線(xiàn)性信號(hào)來(lái)說(shuō)，包絡(luò)的均值線(xiàn)和真正的局部平均線(xiàn)不一樣，某些非對(duì)稱(chēng)波形依然存在，一次篩選無(wú)法消除非對(duì)稱(chēng)波形。鑒于上述原因，必須進(jìn)行重復(fù)多次篩選。

sD的門(mén)限值要根據(jù)具體的信號(hào)來(lái)設(shè)定，門(mén)限的參考數(shù)值可取O．2一O．3。前后2次篩選結(jié)果的如值達(dá)到預(yù)先設(shè)定值，則停止篩選過(guò)程。第l步篩選結(jié)果得到的C，應(yīng)是信號(hào)的特征時(shí)間尺度最小的高頻IMF分量。接下來(lái)，就把C．從原序列中分離出去，得到一個(gè)剩余序列r1(t)：

2.3 小波包理論及濾波算法過(guò)程

3.仿真結(jié)果及分析

本文以一組實(shí)測(cè)機(jī)床放電電流全過(guò)程電火花信號(hào)作為擾動(dòng)信號(hào)樣本進(jìn)行EMD分解，得到5個(gè)IMF分量和一起趨勢(shì)分量R如圖1所示。根據(jù)實(shí)際經(jīng)驗(yàn)門(mén)限值通常去O．2～0．3。本次實(shí)驗(yàn)的EMD門(mén)限值取0．2。

3．1蹦F分量的周期性分析

為了分析IMF分量的周期性特征，本文采用基于自柜關(guān)函數(shù)法求本征模態(tài)函數(shù)的功率譜進(jìn)而尋找周期性‘13‘14。。在MATIAB中對(duì)各個(gè)分量進(jìn)行功率譜分析，將獲得的最大譜對(duì)應(yīng)的周期如表1所示。

圖1 本征模態(tài)函數(shù)(IMF分量和趨勢(shì)分量)

表1 自相關(guān)函數(shù)求功率譜尋找電流時(shí)間序列周期

表1結(jié)果表明：本組數(shù)據(jù)EMD分解后的每個(gè)細(xì)致分量IMF都表現(xiàn)出一定的周期性。在最大落后步長(zhǎng)肘=10，30，50不同情況下，IMF分量表現(xiàn)出不同的周期性，這是因?yàn)樽畲舐浜蟛介L(zhǎng)越小導(dǎo)致時(shí)間序列越短，因而在基于自相關(guān)函數(shù)的功率譜最大譜上面對(duì)應(yīng)的周期越小；但在低頻分量IMF4和IMF5，對(duì)于最大落后步長(zhǎng)對(duì)其周期性的影響減弱。深入分析IMF2和IM乃分量周期在不同M值的情況下周期跳躍性較大，結(jié)合圖1中IMF2和IMF3的幅值是其余分量的10和100倍，假定IMF2和IMF3可能是原始信號(hào)的主成分。

3．2肼F分量的顯著性和相關(guān)系數(shù)分析

為了更加準(zhǔn)確地表征每個(gè)IMF分量與原電流序列統(tǒng)有無(wú)顯著性差別和相關(guān)性程度¨5|，在MATLAB中對(duì)每個(gè)分量進(jìn)行顯著性檢驗(yàn)和相關(guān)系數(shù)測(cè)定，結(jié)果如表2所示。

表2 m伍分量與原電流時(shí)間序列的顯著性檢驗(yàn)和相關(guān)系數(shù)

表2結(jié)果表明：在每個(gè)IMF分量對(duì)原序列做顯著性檢驗(yàn)時(shí)，發(fā)現(xiàn)顯著性水平日值都為0，說(shuō)明所有的IMF分量都通過(guò)了的顯著性檢驗(yàn)，結(jié)合上述的周期性分析，理論上認(rèn)為本組數(shù)據(jù)是周期性振蕩的。然而，每個(gè)分量的顯著性概率卻不相同，表明每個(gè)分量對(duì)于原序列的影響有大小之分；結(jié)合相關(guān)系數(shù)的數(shù)據(jù)特點(diǎn)，有理由證實(shí)IMF2和IMF3是原序列的主成分，并且其攜帶了原信號(hào)大部分信息。

3．3 mF分量的逆變換和誤差分析

圖2 用IMF分量和趨勢(shì)項(xiàng)合成原電流時(shí)間序列和誤差分析

表3誤差分析

圖2和表3聯(lián)合結(jié)果表明：IMF2+IMF3分量雖然有誤差但在可接受范圍內(nèi)，能較好好擬合原電流時(shí)間序列，說(shuō)明IMF2和IMF3是原信號(hào)的主要分量，實(shí)驗(yàn)驗(yàn)證了前面理論分析的正確性。當(dāng)所有的IMF分量相加后得到的信號(hào)與原信號(hào)的總誤差率通過(guò)理論計(jì)算為4．3％，說(shuō)明EMD分解法不僅可以表現(xiàn)信號(hào)的本質(zhì)性質(zhì)，同時(shí)不會(huì)減少信號(hào)的完備性；也實(shí)現(xiàn)了EMD分解法的逆變換。

4.應(yīng)用分析

為了更清晰地分析EMD分解法的有效性，將同時(shí)利用小波分析和EMD分解法對(duì)電火花數(shù)據(jù)進(jìn)行分析。

4．1 小波包的電火花應(yīng)用分析

以采樣時(shí)間為20個(gè)工頻周期，采樣點(diǎn)Ⅳ=1 000，采樣頻率，=1 kHz，取信號(hào)額定頻率為50 Hz。本文小波基函數(shù)分解是在3層分解系數(shù)基礎(chǔ)上展開(kāi)。利用MATLAB的強(qiáng)大信號(hào)處理功能，實(shí)現(xiàn)上述過(guò)程，結(jié)果如圖3和4所示。

圖3 電火花50 A信號(hào)biol4．4小波分解

4．2 EMD分解的電火花應(yīng)用分析

利用EMD分解法對(duì)實(shí)測(cè)的機(jī)床電流放電過(guò)程進(jìn)行分段分割分析，并分析每個(gè)分量的頻譜特性，找到各分量的物理意義，從而更好地理解機(jī)床電流放電的內(nèi)部過(guò)程。EMD門(mén)限值取0．2，利用MATLAB繪制結(jié)果如圖5和6所示。

圖4 電火花50 A信號(hào)db2小波分解

4．3小結(jié)

從圖3和4可以看到：db2小波在對(duì)電火花50A信號(hào)的濾波效果不如bior4．4好，然而在頻率保持方面db2小波變現(xiàn)的更為出色。bioI’4．4小波在不同閥值對(duì)電火花去噪影響情況下，默認(rèn)閥值的波形沒(méi)有調(diào)節(jié)閥值的好，而默認(rèn)閥值的卻濾除了其余諧波的影響。

圖5和6結(jié)果表明：對(duì)于不同的機(jī)床運(yùn)行過(guò)程，從IMFl分量到IMF5分量，都可以逐級(jí)清晰地分解出原電流時(shí)間序列不同階段的高頻和低頻分量。結(jié)合圖1的趨勢(shì)項(xiàng)以及上一節(jié)的理論分析，可以看出經(jīng)過(guò)EMD分解后的電流序列個(gè)分量隨著時(shí)間變化呈現(xiàn)一些周期性變化。對(duì)于各分量可以很容易看出：IMFl包絡(luò)了信號(hào)的5次、7次和9次諧波含量，IMF2包絡(luò)了號(hào)的3次、5次和7次諧波含量；IMF3包絡(luò)了信號(hào)的基波含量；IMF4和IMF5包絡(luò)了信號(hào)的低頻分量。

圖5 未放電電流時(shí)間序列EMD分解和頻譜

圖6放電過(guò)程中電流時(shí)間序列EMD分解和頻譜

對(duì)實(shí)測(cè)信號(hào)的試驗(yàn)表明，采用需選擇合適的小波包算法，并給定合適的閾值才能分析出點(diǎn)火花信號(hào)的部分特性。與小波分析法相比，一方面，EMD分解不僅無(wú)需信號(hào)的先驗(yàn)知識(shí)和固定的基函數(shù)；且可獲得電火花信號(hào)的內(nèi)在放電特性。另一方面，EMD分解法在提取機(jī)床運(yùn)行的早期信號(hào)、微弱信號(hào)和突變信號(hào)等各方面信號(hào)更具優(yōu)越性和高效性，提取結(jié)果更能滿(mǎn)足實(shí)際工程應(yīng)用需求。

5.結(jié)論

EMD分解法是一種新的序列分解方法，本文以一組機(jī)床放電電流特性數(shù)據(jù)為基礎(chǔ)，借助MA7ⅡAB平臺(tái)對(duì)EMD分解過(guò)程和其分量進(jìn)行細(xì)致的理論和實(shí)驗(yàn)對(duì)比論證，得出其在保證信號(hào)完備性的基礎(chǔ)上，具備處理間歇性、非平穩(wěn)和非線(xiàn)性序列的良好效果，同時(shí)具有直觀、簡(jiǎn)潔和適應(yīng)性強(qiáng)等優(yōu)點(diǎn)。因此EMD分解方法特別適用于電火花這類(lèi)非平穩(wěn)與非線(xiàn)性信號(hào)的處理給研究電火花加工工藝參數(shù)和加工工藝效果之間的關(guān)系的研究人員提供了一種新思路。從而實(shí)現(xiàn)電火花加工向高效、高精、自動(dòng)化方向發(fā)展。EMD分解法剛剛起步，研究領(lǐng)域眾多，本文缺乏將EMD分解法與經(jīng)典信號(hào)分解法的對(duì)比性分析，將在以后的研究中改進(jìn)。

投稿箱：
如果您有機(jī)床行業(yè)、企業(yè)相關(guān)新聞稿件發(fā)表，或進(jìn)行資訊合作，歡迎聯(lián)系本網(wǎng)編輯部，郵箱：skjcsc@vip.sina.com

更多相關(guān)信息

名企推薦

業(yè)界視點(diǎn)

| 更多

行業(yè)數(shù)據(jù)

| 更多

博文選萃

| 更多

数控机床市场网-专业的数控机床网站

數(shù)控機(jī)床