華科-機(jī)械工程測試信息信號(hào)分析-課件-專題2-小波分析

1、MEASUREMENTINFORMATION SIGNAL ANALYSIS IN MECHANICAL ENGINEERING,機(jī)械工程測試?信息?信號(hào)分析,機(jī)械科學(xué)與工程學(xué)院 機(jī)械電子信息工程系,Wavelet Analysis,小波分析,本周討論內(nèi)容,小波相關(guān)歷史回顧連續(xù)小波變換 (Continuous Wavelet Transform)離散小波變換 (Discrete Wavelet Transform)

2、小波包變換 (Wavelet Packet Transform)應(yīng)用 (Applications),回顧1：Fourier變換,Jean B. Joseph Fourier(1768-1830)兩大貢獻(xiàn)：,不足：對(duì)整個(gè)信號(hào)進(jìn)行分析，不能說明頻率隨時(shí)間變化的關(guān)系。,“周期信號(hào)都可表示為諧波關(guān)系的正弦信號(hào)的加權(quán)和”,“非周期信號(hào)都可用正弦信號(hào)的加權(quán)積分表示”,回顧2：短時(shí)Fourier變換,Dennis Gabor(1900-1979)于

3、1946年采用滑動(dòng)窗口觀察信號(hào)，計(jì)算局部時(shí)間段內(nèi)的FT變換，得局部頻率。,or,?,不足：不利于分析時(shí)變信號(hào)高頻信號(hào)持續(xù)時(shí)間短，而低頻長。希望對(duì)于高頻采用小的時(shí)間窗，低頻使用大時(shí)間窗進(jìn)行分析。,特點(diǎn): (1) 時(shí)頻分辨率固定，取決于窗函數(shù)的時(shí)寬和帶寬 (2) 窗函數(shù)的時(shí)寬△t 和帶寬△f 滿足不確定性原理,不確定性原理(Heisenberg Uncertainty Principle),Cohen, L.. T

4、ime-frequency distributions-a review. Proceedings of the IEEE Vol. 77(7): 941 – 981, Jul 1989.,時(shí)間和頻率的分辨率不能同時(shí)無限高。魚和熊掌不可兼得。選魚還是熊掌呢？,現(xiàn)象:電子位置與動(dòng)量不能同時(shí)準(zhǔn)確的測量,結(jié)論：,模型：,,,回顧3：小波變換,1984年法國的地球物理學(xué)家Jean Morlet在進(jìn)行石油勘探的地震數(shù)據(jù)處理分析時(shí)與法國理論物理學(xué)家

5、A.Grossman一起提出了小波變換(wavelet transform,WT)的概念。1910年Haar給出了矩形基函數(shù)即后來的Haar小波函數(shù)。1989年S.Mallat和Y. Meyer提出了多分辨分析。,基本數(shù)學(xué)概念,向量空間和函數(shù)空間,復(fù)共軛,向量空間,函數(shù)空間,基,向量正交,線性表示,函數(shù)正交,正交基函數(shù),函數(shù)分解,小波的基本思想：函數(shù)的平移與伸縮,,移位,,收縮,,擴(kuò)展,,,收縮 + 平移,擴(kuò)展 + 平移,大尺度看

6、全局,小尺度看細(xì)節(jié),尺度因子a移位因子b,連續(xù)小波變換(CWT),,連續(xù)小波變換(CWT),基本小波函數(shù)的平移與伸縮,基本小波函數(shù),基本小波函數(shù),,滿足允許條件,Haar小波,Mexican Hat小波,Morlet小波,幾種基本小波,,常用于地震波信號(hào)處理,常用于教學(xué),常用于信號(hào)瞬態(tài)成份識(shí)別,基本小波條件 波形能量為1,能量不變性,小波函數(shù)及其頻譜,,收縮,,擴(kuò)展,,FT,,擴(kuò)展,,收縮,CWT的性質(zhì),疊加性質(zhì) 若

7、 的CWT是 ,那么 的CWT是時(shí)移不變性 若 的CWT是 ,那么 的CWT是尺度伸縮共變性 若 的CWT是 ,則 的CWT是Moyal原理 若 的CWT是 ,則 特別地,當(dāng)

8、 則有,,,,,CWT的物理解釋,帶通濾波器解釋 信號(hào)x(t)通過中心分析頻率為f的帶通濾波器 濾波器的帶寬△f與中心分析頻率f有關(guān)，兩者之比為常數(shù),CWT的分辨率,,時(shí)間分辨率由伸縮小波函數(shù)的時(shí)寬△t 所決定 頻率分辨率由伸縮小波函數(shù)的帶寬△f 所決定 時(shí)間分辨率和頻率分辨率都隨中心分析頻率f而變化 對(duì)于高頻信號(hào)，用窄時(shí)窗（寬帶寬），則時(shí)域分辨率高 對(duì)于低頻信號(hào)，用寬時(shí)窗（窄帶寬），則頻率分辨率高 分析頻帶按指

9、數(shù)劃分 (如二進(jìn)劃分) ，每一帶寬對(duì)應(yīng)一個(gè)尺度 (級(jí)),CWT與STFT的比較,,CWT的數(shù)學(xué)解釋,,CWT是原信號(hào)與一系列基函數(shù)的內(nèi)積 基函數(shù)是由基本小波經(jīng)平移 (因子b) 和伸縮 (因子a) 得到的系列函數(shù) CWT表明了原信號(hào)與某個(gè)特定基函數(shù)的“相似”程度,原信號(hào)的重構(gòu)（連續(xù)小波逆變換） 重構(gòu)需要基本小波滿足允許條件 原信號(hào)也可看成是按基函數(shù)的展開 (展開系數(shù)就是CWT),,小波變換示意,CWT實(shí)現(xiàn)過程和結(jié)果,1.取一個(gè)小波

10、與信號(hào)的最前面部分比較;,2.計(jì)算小波系數(shù)C,C代表小波和這段數(shù)據(jù)的相關(guān)性， 即C越大,兩者越相似;,3.移動(dòng)小波,重復(fù)步驟1和2,一直遍歷整個(gè)數(shù)據(jù);,4.對(duì)小波進(jìn)行縮放,重復(fù)步驟1到3;,5.在所有小波尺度下,重復(fù)上述步驟.,離散小波變換(DWT),,基本小波函數(shù)的離散平移與伸縮,,小波“級(jí)”：j 小波級(jí)大，則尺度小,在每個(gè)可能的縮放因子和平移參數(shù)下計(jì)算小波系數(shù)，其計(jì)算量相當(dāng)大，將產(chǎn)生驚人的數(shù)據(jù)量，而且有許多數(shù)據(jù)是無用的。如果

11、縮放因子和平移參數(shù)都選擇為2j（j>0且為整數(shù)）的倍數(shù)， 即只選擇部分縮放因子和平移參數(shù)來進(jìn)行計(jì)算， 就會(huì)使分析的數(shù)據(jù)量大大減少。,離散小波變換(DWT),（二進(jìn)）離散小波變換(DWT),原信號(hào)的重構(gòu)（二進(jìn)離散小波逆變換）,這實(shí)際上是對(duì)CWT的離散化，并不是真正意義上的離散變換,DWT要求,能提供足夠信息分析信號(hào)和重構(gòu)信號(hào)；高的計(jì)算效率；容易實(shí)現(xiàn)；能在不同頻帶內(nèi)以不同的分辨率分析信號(hào)；能將信號(hào)分解成粗糙部分和細(xì)節(jié)部分。,基

12、本尺度函數(shù),,基本尺度函數(shù),基本小波函數(shù),,平移+伸縮,兩尺度方程,{hk} 和{gk}稱為兩尺度序列，且均與j無關(guān)， 只與尺度函數(shù)和小波函數(shù)相關(guān) {hk} 和{gk}都是由有限個(gè)數(shù)組成的序列 如果小波函數(shù)簇全正交，則{hk} 和{gk}滿足,尺度函數(shù)構(gòu)建小波函數(shù),,DWT的快速算法（S. Mallat 1989）,遞推分解與重構(gòu)公式（金字塔算法）,只對(duì) cj 進(jìn)行繼續(xù)分解，分解過程是遞推的。 對(duì)于數(shù)字信號(hào)的分解，可以直

13、接將其定義為c0，則分解過程是完全離散化的，而且分解過程的計(jì)算只需要知道兩尺度序列，甚至不涉及尺度函數(shù)和小波函數(shù)的具體形式。,快速分解,快速重構(gòu),DWT快速算法的物理意義,遞推分解公式的頻域表示（共軛濾波器）,待分析信號(hào)通過H(f)和G(f)分別進(jìn)行低通和高通濾波 濾波的結(jié)果是得到一組低頻分量和一組高頻分量 每次獲得的低頻分量和高頻分量的長度都是原信號(hào)長度的一半 可以看作是濾波后進(jìn)行了“隔點(diǎn)采樣”,DWT的分解結(jié)果,二進(jìn)頻帶劃分,

14、,,級(jí),時(shí),將高頻成份不斷剝離，可得低頻部分常用于降噪,尺度函數(shù)與小波函數(shù)的關(guān)系,,正交補(bǔ)空間,尺度函數(shù)空間的嵌套性,,,,,注意j前面的正負(fù)號(hào),DWT的快速算法降噪實(shí)例過程,第2次變換,第3次變換,生成信號(hào),定義濾波器db(2)參數(shù),信號(hào)濾波,隔點(diǎn)采樣,第1次變換,512采樣點(diǎn),注意H0和H1不變,小波包變換,基本尺度函數(shù),基本小波函數(shù),小波包,快速分解與重構(gòu)遞推算法,改進(jìn)Mallat算法以提高高頻區(qū)域的分辨率,Mallat算法,生

15、成小波包基,,小波包變換的物理解釋,小波包基實(shí)例,3級(jí)Haar小波包基,小波（包）分析的工程應(yīng)用,應(yīng)用基礎(chǔ) 具有頻帶分離功能 可以在緩變信號(hào)中檢測出突變信號(hào) 可以在高頻噪聲中檢測出全局緩變信號(hào) 小波變換用時(shí)級(jí)圖表示，可以得到二進(jìn)劃分頻帶內(nèi)的局部能量 小波包變換用時(shí)級(jí)圖表示，可以得到線性劃分頻帶內(nèi)的局部能量,部分工程應(yīng)用 鋼絲繩斷絲檢測 切削顫振識(shí)別 齒輪故障檢測 銑刀破損監(jiān)測 發(fā)動(dòng)機(jī)故障振動(dòng)監(jiān)測,分析高頻的好處，提取

16、信號(hào)特征有個(gè)性：有棱有角。無個(gè)性：泯然眾人矣,基于小波分析的齒輪故障檢測(1),齒輪故障,,,,Pinion,Broken Tooth,基于小波分析的齒輪故障檢測(2),Abnormal,,去除嚙合頻率分量后振動(dòng)信號(hào)的時(shí)域同步平均,基于小波分析的齒輪故障檢測(3),D4小波分解的時(shí)級(jí)均方圖,Abnormal,,基于小波分析的齒輪故障檢測(4),去除嚙合頻率分量后振動(dòng)信號(hào)的時(shí)域同步平均,基于小波分析的齒輪故障檢測(5),D4小波分解的時(shí)

17、級(jí)均方圖,Abnormal,,基于小波分析的齒輪故障檢測(6),調(diào)諧小波分解的時(shí)級(jí)均方圖,Abnormal,,基于小波分析的齒輪故障檢測(7),高斯小波分解的時(shí)級(jí)均方圖,發(fā)動(dòng)機(jī)故障診斷系統(tǒng)--信號(hào)采集系統(tǒng),信號(hào)采集系統(tǒng)構(gòu)成及傳感器安裝位置,對(duì)象 四沖程柴油機(jī)，工作轉(zhuǎn)速1500轉(zhuǎn)/分鐘,信號(hào)采集方法,采樣方式 同步采樣（測量空間720度）,采樣頻率 每個(gè)通道采樣頻率約25.6kHz,采樣點(diǎn)數(shù) 每個(gè)通道每次循環(huán)采樣2

18、048個(gè)點(diǎn),觸發(fā)同步信號(hào)用的電渦流傳感器，安裝在柴油機(jī)曲軸上,加速度傳感器，安裝在柴油機(jī)缸蓋上（進(jìn)氣門和排氣門）,傳感器,信號(hào)調(diào)理模塊,電荷放大器和濾波器等,數(shù)據(jù)采集卡,美國RTD公司AD3110卡,12位采樣精度，16個(gè)A/D采樣通道，最高采樣頻率200kHz,基于小波包分析的發(fā)動(dòng)機(jī)故障診斷(1),研究背景 利用缸蓋表面振動(dòng)信號(hào)進(jìn)行柴油機(jī)潛在故障的診斷，故障分類 缸蓋振動(dòng)信號(hào)具有局部沖擊特性，是非平穩(wěn)時(shí)變信號(hào) 希望從整循環(huán)振動(dòng)

19、信號(hào)中直接提取特征參數(shù),Acceleration (m/s2),Time (ms),Cylinder 1,Cylinder 2,EVC IVC Combustion EVO IVO EVC,IVC Combustion EVO IVO EVC IVC,TDC BDC TDC BDC

20、 TDC BDC (of cylinder 1),Timing information important.,,,,基于小波包分析的發(fā)動(dòng)機(jī)故障診斷(2),診斷過程 形成特征向量,基于小波包分析的發(fā)動(dòng)機(jī)故障診斷(3),Pattern 1 - Healthy Pattern 2Pattern 3Pattern 4Pattern 5Pattern 6,實(shí)際測量波形,氣門間隙偏小,氣門間隙偏大,

21、氣門漏氣,進(jìn)氣門外彈簧斷裂,氣門內(nèi)彈簧斷裂,基于小波包分析的發(fā)動(dòng)機(jī)故障診斷(4),Pattern 2,實(shí)測波形的小波包分解,Pattern 1 - Healthy,基于小波包分析的發(fā)動(dòng)機(jī)故障診斷(5),Pattern 1 - Healthy Pattern 2Pattern 3Pattern 4Pattern 5Pattern 6,整循環(huán)特征向量,基于小波包分析的發(fā)動(dòng)機(jī)故障診斷(6),Original Algorithm,Wa

22、velet Packet Coefficients (Scale/Frequency Values),基于小波包分析的發(fā)動(dòng)機(jī)故障診斷(7),Original Algorithm with Re-coding,Wavelet Packet Coefficients (Scale/Frequency Values),基于小波包分析的發(fā)動(dòng)機(jī)故障診斷(8),Improved Algorithm,Wavelet Packet Coefficien

23、ts (Scale/Frequency Values),小結(jié)（一些人和事）,一個(gè)原理一個(gè)思想一些概念和方法,入門參考書,Mallat S. A wavelet tour of signal process(3rd edition).機(jī)械工業(yè)出版社，2010年C. S. Burrus, R. A. Gopinath, H.T. Guo. Introduction to wavelets and wavelet transforms.

24、 小波與小波變換導(dǎo)論.(程正興譯）.機(jī)械工業(yè)出版社，2008年Robert X. Gao, R.Q. Yan. Wavelets: theory and applications for manufacturing. Springer,2011.K.P. Soman, K.I. Ramachandran. Insights into wavelets from theory to practice(2nd edition). Pre