Laplace小波學(xué)習(xí)整理.doc_第1頁(yè)
Laplace小波學(xué)習(xí)整理.doc_第2頁(yè)
Laplace小波學(xué)習(xí)整理.doc_第3頁(yè)
Laplace小波學(xué)習(xí)整理.doc_第4頁(yè)
Laplace小波學(xué)習(xí)整理.doc_第5頁(yè)
已閱讀5頁(yè),還剩2頁(yè)未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

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

文檔簡(jiǎn)介

Laplace小波特征波混合基分解的精髓:使用與信號(hào)波形最匹配的基函數(shù)對(duì)信號(hào)進(jìn)行分解、提取隱含故障特征。Laplace小波定義: 該小波為單邊衰減的復(fù)指數(shù)小波,定義如下圖所示:在MATLAB中程序:Fs=1e4;t=1/Fs:1/Fs:0.1;f=500;w=2*pi*f;A=0.08; %歸一化小波函數(shù)ep=0.03; %粘滯阻尼比tal=0.1; %時(shí)間參數(shù)y=A.*exp(-ep/(sqrt(1-ep2).*w.*(t-tal).*exp(-j.*w.*(t-tal);其中,y為復(fù)指數(shù)小波,歐拉公式知y包含實(shí)部和虛部?jī)刹糠?。Laplace小波特性: 小波的波形:對(duì)應(yīng)的實(shí)部波形圖及頻譜圖: 實(shí)部Re(y)實(shí)際上是一個(gè)高通濾波器。Laplace小波基函數(shù)庫(kù)濾波方法(小波檢測(cè)方法):思想:信號(hào)x(t)是某個(gè)系統(tǒng)S的輸出,信號(hào)x(t)與x(t)的Laplce小波原子的內(nèi)積,可以估計(jì)出它們之間的相似性,從而得到系統(tǒng)S的模態(tài)參數(shù),與的頻率、阻尼特性的對(duì)應(yīng)關(guān)系。內(nèi)積:若x(t)與完全相關(guān),那么它們之間的夾角。可以定義一個(gè)相關(guān)系數(shù)來(lái)量化與之間的夾角。其中為L(zhǎng)aplace小波基函數(shù)庫(kù)中的一個(gè),一個(gè)對(duì)應(yīng)一個(gè)。補(bǔ)充知識(shí):內(nèi)積定義: 其中二范數(shù)定義:那么應(yīng)用例子:通過(guò)相關(guān)濾波法寫出MATLAB程序%-模擬信號(hào)檢驗(yàn),用Laplace小波檢驗(yàn)自由度響應(yīng)的脈沖響應(yīng)信號(hào)%-Laplace小波濾波例子%clear;clc;%An=0.01;%噪聲幅度z0=0.05; %阻尼比t0=0.5; %固有時(shí)間參數(shù)f0=10; %固定頻率Fs=200; %采樣率t=-5:1/Fs:5; %時(shí)間范圍L=length(t); %采樣點(diǎn)數(shù)s=zeros(1,L); %構(gòu)造信號(hào)lap=zeros(1,L); %構(gòu)造小波原子函數(shù)noise =zeros(1,L); %隨機(jī)噪聲%-構(gòu)造脈沖響應(yīng)信號(hào),噪聲noise=An*random(Normal,0,1,1,L);%-自由度系統(tǒng)的脈沖響應(yīng)信號(hào)for i=1:L if t(i)=t0 s(i)=exp(-z0/(sqrt(1-z02)*2*pi*f0.*(t(i)-t0)*(sin(2*pi*f0.*(t(i)-t0)+noise(i); else s(i)=noise(i); endend%-構(gòu)造Laplace小波特征波庫(kù)F=5:0.5:20; %頻率組Z=0.005:0.005:0.2,0.3:0.1:0.9; %阻尼比組T=-5:0.1:5; %時(shí)間參數(shù)組for m=1:length(F) for n=1:length(Z) for p=1:length(T) for i=1:L if t(i)=T(p) lap(i)=exp(-Z(n)/(sqrt(1-Z(n)2)*2*pi*F(m)*(t(i)-T(p). *(sin(2*pi*F(m)*(t(i)-T(p); else lap(i)=1e-5; end end norm_lap(m,n,p)=norm(lap,2); norm_s(n,m,p)=norm(s,2); kr(m,n,p)=sqrt(2)*sum(lap.*s)/(norm(s,2)*norm(lap,2);% plot3(F(m),Z(n),kr(m,n,p);hold on; end endend%norm求歐幾里得范數(shù)%-求信號(hào)s的二范數(shù),小波原子的二范數(shù)% 求 kr最大值并找出最大值對(duì)應(yīng)的坐標(biāo)size_kr=size(kr);max_kr=max(kr(:);i_max=find(kr=max_kr);c,k,g=ind2sub(size_kr,i_max);gy_f=F(c);gy_z=Z(k);gy_t=T(g);gy_s=zeros(1,L); %固有參數(shù)對(duì)應(yīng)的結(jié)果% 根據(jù)得到的固有參數(shù)得到的結(jié)果for i=1:L if t(i)=gy_t gy_s(i)=exp(-gy_z/(sqrt(1-gy_z2)*2*pi*f0*(t(i)-gy_t)*(sin(2*pi*gy_f*(t(i)-gy_t); else gy_s(i)=0; endendfigure(1)subplot(1,2,1)plot(t,s,-r);grid on;subplot(1,2,2)plot(t,gy_s,-k);grid on;% 每個(gè)T時(shí)刻對(duì)應(yīng)的相關(guān)系數(shù)峰值kr(T)for p=1:length(T)-1 krT(p)=max(max(kr(:,:,p); row,clum=find(kr(:,:,p)=krT(p); f_lap(p)=F(row); z_lap(p)=Z(clum);endfigure(2)subplot(3,1,1)plot(T(1,1:100),krT); %每個(gè)時(shí)刻T的相關(guān)系數(shù)kr的峰值,對(duì)應(yīng)的圖像title(每個(gè)時(shí)刻T的相關(guān)系數(shù)kr峰值);grid on;xlabel(T,時(shí)間);ylabel(kr峰值);subplot(3,1,2)plot(T(1,1:100),f_lap);title(每個(gè)時(shí)刻T的相關(guān)系數(shù)頻率值);grid on;xlabel(T,時(shí)

溫馨提示

  • 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ù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 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)論