影響局部均值分解算法分解精度的若干因素.doc_第1頁(yè)
影響局部均值分解算法分解精度的若干因素.doc_第2頁(yè)
影響局部均值分解算法分解精度的若干因素.doc_第3頁(yè)
影響局部均值分解算法分解精度的若干因素.doc_第4頁(yè)
影響局部均值分解算法分解精度的若干因素.doc_第5頁(yè)
免費(fèi)預(yù)覽已結(jié)束,剩余1頁(yè)可下載查看

下載本文檔

版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

影響局部均值分解算法分解精度的若干因素局部均值分解LMD(Local mean decomposition)是一種新的自適應(yīng)時(shí)頻分析方法。LMD可以自適應(yīng)地將任何一個(gè)復(fù)雜信號(hào)分解為若干個(gè)瞬時(shí)頻率具有物理意義的PF(Product function)分量之和,從而獲得原始信號(hào)完整的時(shí)頻分布,其本質(zhì)上是將多分量的信號(hào)自適應(yīng)地分解為若干個(gè)單分量的調(diào)幅調(diào)頻信號(hào)(包絡(luò)信號(hào)和純調(diào)頻信號(hào))之和,非常適合處理非平穩(wěn)、非線性信號(hào),特別是多分量的調(diào)幅調(diào)頻信號(hào)。局部均值分解的基本理論和算法很容易理解,在許多文獻(xiàn)上都給出了詳細(xì)的介紹,在這里不做過(guò)多的闡述,重點(diǎn)討論影響局部均值分解算法的一些問(wèn)題。局部均值分解作為一種新提出的時(shí)頻分析方法,有許多問(wèn)題有待于進(jìn)一步研究,如原始信號(hào)的采樣頻率、求局部均值函數(shù)和包絡(luò)函數(shù)用到的滑動(dòng)平均算法以及端點(diǎn)效應(yīng)。采樣頻率不同會(huì)影響LMD分解精度,尤其對(duì)瞬時(shí)頻率影響較大,采用不同的滑動(dòng)跨度和端點(diǎn)處理方法就會(huì)得到不同的分解結(jié)果,有時(shí)甚至?xí)斐伤惴ǖ牟皇諗俊O旅鎸⒁跃唧w實(shí)例分析二者對(duì)滑動(dòng)平均算法的影響。具體信號(hào)如下:1. 采樣頻率對(duì)LMD分解精度的影響 上面兩個(gè)圖是采樣間隔為0.2s時(shí)得到的PF分量、包絡(luò)信號(hào)分量和對(duì)應(yīng)的瞬時(shí)頻率分量下面兩個(gè)圖是采樣間隔為1s時(shí)得到的PF分量、包絡(luò)信號(hào)分量和對(duì)應(yīng)的瞬時(shí)頻率分量2. 滑動(dòng)平均跨度選擇對(duì)LMD分解精度的影響滑動(dòng)平均的跨度不僅關(guān)系到LMD分解精度,如果滑動(dòng)平均跨度選擇不合理,有可能造成LMD算法不收斂。這是因?yàn)榫钟虬j(luò)函數(shù)容易受到滑動(dòng)平均算法的影響,滑動(dòng)平均的跨度不同,局域包絡(luò)函數(shù)就會(huì)不一樣。在求解每個(gè)PF分量時(shí)采用的是純調(diào)頻信號(hào)的判據(jù),即要求局部包絡(luò)函數(shù)滿足,因此局域包絡(luò)函數(shù)不同,循環(huán)迭代次數(shù)也就不同。其中增減量的取值范圍,需要根據(jù)不同的信號(hào)和不同的精度要求來(lái)設(shè)定,的值越小,計(jì)算量就越大,LMD分解的精度越高,根據(jù)經(jīng)驗(yàn)一般取0.0010.01。以上述信號(hào)為例,在消除端點(diǎn)影響的前提下,研究滑動(dòng)平均跨度對(duì)LMD分解的影響。分別采用不同的滑動(dòng)平均跨度方法求取PF1,觀察其循環(huán)次數(shù)和滑動(dòng)平均跨度的關(guān)系,如下如圖所示:滑動(dòng)平均跨度(Moving Average Span)和循環(huán)次數(shù)(Loop Times)關(guān)系圖當(dāng)滑動(dòng)跨度為5點(diǎn)時(shí),循環(huán)7次結(jié)束,當(dāng)滑動(dòng)跨度超過(guò)5時(shí),循環(huán)次數(shù)逐步增加,直到滑動(dòng)跨度超過(guò)11點(diǎn)(循環(huán)16次)時(shí),循環(huán)次數(shù)顯著增加。當(dāng)滑動(dòng)跨度為13點(diǎn)時(shí),循環(huán)需要63次才能結(jié)束。當(dāng)滑動(dòng)跨度為15點(diǎn)時(shí),循環(huán)超過(guò)100次時(shí)仍未求出純調(diào)頻信號(hào),由程序強(qiáng)制結(jié)束??梢?jiàn),滑動(dòng)平均跨度的選擇會(huì)對(duì)循環(huán)產(chǎn)生一定的影響,當(dāng)選擇不合理時(shí)會(huì)增加計(jì)算量,甚至使算法不收斂。但是由圖也可以看出,對(duì)具體信號(hào)而言,在某一范圍內(nèi),滑動(dòng)平均跨度對(duì)PF1分量影響并不是很大,僅僅是影響局部包絡(luò)函數(shù)的個(gè)數(shù)(循環(huán)次數(shù))。從另一方面講,在合理的取值范圍內(nèi),LMD分解精度是與計(jì)算量成正比的。所以滑動(dòng)平均跨度的選擇要綜合考慮分解精度和計(jì)算量的影響。(對(duì)該信號(hào)如果采樣時(shí)間間隔為0.2s時(shí),跨度選擇在25點(diǎn)左右時(shí)較為合理,不會(huì)對(duì)分解結(jié)果產(chǎn)生大的影響,當(dāng)跨度超過(guò)33點(diǎn)時(shí)波形會(huì)產(chǎn)生明顯的誤差)。根據(jù)經(jīng)驗(yàn)滑動(dòng)平均跨度的選擇一般采用兩種方法,一是取相鄰極值點(diǎn)最長(zhǎng)距離的三分之一;二是設(shè)定為相鄰極值點(diǎn)的最短距離。兩種方法對(duì)不同的振動(dòng)信號(hào)分解精度也不相同,對(duì)具體信號(hào)應(yīng)采用何種方法來(lái)選擇滑動(dòng)平均跨度一般很難確定。但是無(wú)論采用哪一種方法分解精度不會(huì)相差太大。對(duì)于任何一個(gè)信號(hào)序列,當(dāng)采樣頻率和采樣點(diǎn)數(shù)確定以后,其相應(yīng)的滑動(dòng)平均跨度的最合理取值就可以通過(guò)提出的兩種方法計(jì)算得到。經(jīng)過(guò)大量實(shí)驗(yàn)的證明,我們發(fā)現(xiàn)在這個(gè)最合理的跨度值附近取值就不會(huì)對(duì)LMD算法造成大的影響。圖a 七點(diǎn)滑動(dòng)平均 圖b 九點(diǎn)滑動(dòng)平均 圖c 十一點(diǎn)滑動(dòng)平均以上面提到的信號(hào)為例。由于我們已經(jīng)對(duì)端點(diǎn)做了鏡像處理,所以端點(diǎn)處數(shù)據(jù)不會(huì)產(chǎn)生大的變形,對(duì)于該信號(hào),橫坐標(biāo)為500左右的數(shù)據(jù)極值點(diǎn)相對(duì)其它地方較稀疏,信號(hào)在處理時(shí)會(huì)發(fā)生變化,可以選擇該段作為觀察重點(diǎn)部分。其最合理的跨度取值為7點(diǎn)(采樣間隔較長(zhǎng)、每周期采樣點(diǎn)數(shù)較少,跨度就相對(duì)較?。捎?點(diǎn)滑動(dòng)平均得到的結(jié)果如圖a所示。將滑動(dòng)跨度改為3點(diǎn)、5點(diǎn)、9點(diǎn)我們可以發(fā)現(xiàn)圖形并沒(méi)有大的變化,即在這些值之間選擇滑動(dòng)平均跨度均不會(huì)對(duì)LMD分解產(chǎn)生大的誤差,在這里我們選擇跨度為9點(diǎn)時(shí)的圖形(圖b)與之比較。將跨度值進(jìn)一步增加,取跨度值為11點(diǎn),圖形如圖c所示??缍瘸隽似浜侠砣≈捣秶?,分解結(jié)果產(chǎn)生了一定的誤差。當(dāng)繼續(xù)增加跨度值時(shí),誤差會(huì)更明顯,使信號(hào)嚴(yán)重失真。LMD算法在采用滑動(dòng)平均平滑數(shù)據(jù)時(shí),不能事先定義滑動(dòng)跨度的值。對(duì)不同的信號(hào)或者同一信號(hào)不同的采樣頻率其最合理的滑動(dòng)跨度值是不同的。這個(gè)值可以用上述提到的兩種方法求得。一般來(lái)說(shuō)對(duì)同一信號(hào),采樣時(shí)間間隔越短,即一個(gè)周期采樣點(diǎn)數(shù)越多,其跨度值就會(huì)越大。3. 端點(diǎn)效應(yīng)對(duì)LMD分解精度的影響由于實(shí)際分析的信號(hào)序列其長(zhǎng)度總是有限的,無(wú)法知道端點(diǎn)以外的信號(hào),所以端點(diǎn)附近的局域均值函數(shù)和局域包絡(luò)函數(shù)只能根據(jù)已知信號(hào)推測(cè),這樣就不可避免的會(huì)產(chǎn)生誤差。LMD是通過(guò)相鄰極值點(diǎn)(一個(gè)極大值和一個(gè)極小值),計(jì)算得到局域均值函數(shù),再經(jīng)過(guò)若干次滑動(dòng)平均得到光滑的局域均值函數(shù)。由求解局部均值函數(shù)公式:可以看出在信號(hào)左端點(diǎn)至第一個(gè)極值點(diǎn)(n1)之間,局域均值函數(shù)是一段未知信號(hào),如下圖箭頭至左端點(diǎn)部分所示(為突出端點(diǎn),截取了原始信號(hào)端點(diǎn)附近的部分信號(hào))。同理,右端點(diǎn)部分也存在同樣的問(wèn)題。LMD端點(diǎn)附近未知包絡(luò)信號(hào)下面就端點(diǎn)效應(yīng)對(duì)LMD分解精度的影響程度做進(jìn)一步討論:LMD分解會(huì)受到端點(diǎn)效應(yīng)的影響,由于局域均值包絡(luò)函數(shù)在端點(diǎn)處均存在一段未知的信號(hào),若對(duì)端點(diǎn)不進(jìn)行處理,在程序運(yùn)行時(shí),會(huì)自動(dòng)給這部分信號(hào)添加一些虛假信息,從而對(duì)LMD分解產(chǎn)生影響。LMD端點(diǎn)效應(yīng)首先發(fā)生在端點(diǎn)附近,然后在迭代過(guò)程中不斷向內(nèi)部擴(kuò)散,迭代次數(shù)越多端點(diǎn)效應(yīng)污染整個(gè)數(shù)據(jù)段的程度就越嚴(yán)重。端點(diǎn)效應(yīng)會(huì)使分解得到的各分量在端點(diǎn)附近產(chǎn)生一些變形,從而使結(jié)果不容易滿足循環(huán)終止條件,增加了循環(huán)次數(shù),嚴(yán)重的時(shí)候會(huì)使數(shù)據(jù)產(chǎn)生嚴(yán)重失真。為減小端點(diǎn)效應(yīng)對(duì)算法的影響,在分解前要對(duì)端點(diǎn)經(jīng)行一定的處理,應(yīng)用最多的處理方法是鏡像延拓算法,鏡像延拓是在端點(diǎn)以外延拓一段信號(hào)。實(shí)際處理的信號(hào)兩端點(diǎn)一般不是極值點(diǎn),這時(shí)候可以采用鏡像延拓的方法進(jìn)行拓展,在LMD算法中只延拓一個(gè)極值點(diǎn)就可以很好的消除端點(diǎn)效應(yīng)的影響。延拓方法如下圖所示,以離端點(diǎn)最近的一個(gè)極值點(diǎn)為對(duì)稱軸(如圖中虛線所示),將離端點(diǎn)次近的極值點(diǎn)向外延拓。在求解局部均值函數(shù)和局部包絡(luò)函數(shù)的時(shí)候?qū)⑦@個(gè)延拓的極值點(diǎn)代入,即可求得完整的局部均值函數(shù)和局部包絡(luò)函數(shù)。極值點(diǎn)延拓示意圖為加深端點(diǎn)效應(yīng)對(duì)LMD分解精度影響的理解,現(xiàn)采用兩種端點(diǎn)處理方法經(jīng)行比較。方法一:極值點(diǎn)延拓處理端點(diǎn),方法二:將端點(diǎn)值直接作為極值經(jīng)行計(jì)算。分解結(jié)果如圖a、圖b所示。圖a 極值點(diǎn)延拓處理得到的各PF分量 圖b 將端點(diǎn)作為極值點(diǎn)得到的各PF分量從圖中我們可以直觀的看出,采用極值點(diǎn)延拓方法,端點(diǎn)效應(yīng)得到了很好的抑制,兩個(gè)PF分量在端點(diǎn)處沒(méi)有引起大的變形。并且其殘余項(xiàng)的幅值也很小,可以忽略。而采用方法二時(shí),雖然分解的PF分量保持了大致的形狀,可是端點(diǎn)部分卻出現(xiàn)了很明顯的變形,端點(diǎn)處的幅值變得很突出。由于端點(diǎn)效應(yīng)的影響,方法二比方法一多分解出了一個(gè)PF分量(PF3),PF3分量除端點(diǎn)部分外,其它部分幅值均很小,可以看作是端點(diǎn)效應(yīng)所帶來(lái)的誤差。可見(jiàn)端點(diǎn)效應(yīng)對(duì)LMD算法影響很大,端點(diǎn)如

溫馨提示

  • 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

最新文檔

評(píng)論

0/150

提交評(píng)論