小波故障診斷系統(tǒng)

1、小波分析故障診斷系統(tǒng)小波分析故障診斷系統(tǒng)目錄 1、小波定義 2、小波分析原理 3、小波診斷故障原理 4、編程芯片、計算機(jī)中如何使用 （快速算法） 5、matlab仿真小波定義： 一種特殊的函數(shù) （小波函數(shù)） “小”：有限寬度 “波”：有震蕩變化，均值為零 自己可不可設(shè)計一個小波？小波分析原理：（畫圖） 傅里葉變換與小波變換比較：公式 Wf:相似系數(shù) 圖解說明：頻率與放大倍數(shù)（尺度）故障位置、類型 傅里葉變換對位置、突變無能為力 故障類型：1、噪聲干擾 2、頻率突變 3、奇點(diǎn) 快速算法-mallat算法互補(bǔ)濾波器高頻：細(xì)節(jié)低頻：近似、逼近層數(shù): 改變小波放大倍數(shù) 即尺度(篩選例子)程序：類似數(shù)字

2、濾波器（卷積）Matlab仿真噪聲Matlab-頻率突變db5小波細(xì)節(jié)高頻有用、與零比較判斷位置Haar小波Matlab奇點(diǎn)在相似系數(shù)中找最值、確定突變位置Haar小波clear;load cuspamax; %載入含有突變點(diǎn)的原始信號whos; %列出當(dāng)前工作空間中所有變量figure(1)plot(cuspamax)xlabel(時間);ylabel(幅值); %以時間為橫軸，幅值為縱軸構(gòu)造圖形title(突變信號);figure(2) %對原始信號使用Daudechies6小波在尺度2,4,8,16,32上進(jìn)行小波變換c,l=wavedec(cuspamax,5,db6); cfd=ze

3、ros(5,1024);for k=1:5 d=detcoef(c,l,k); d=d(ones(1,2k),:); cfd(k,:)=wkeep(d(:),1024)endcfd=cfd(:);I=find(abs(cfd)sqrt(eps);cfd(I)=zeros(size(I);cfd=reshape(cfd,5,1024);colormap(pink(64);img=image(flipud(wcodemat(cfd,64,row);set(get(img,parent),YtickLabel,);title(離散小波變換系數(shù)的絕對值)ylabel(層數(shù))clear;load nea

4、rbrk; %載入頻率突變的近似信號whos;figure(1); %以時間為橫軸，幅? 值為縱軸構(gòu)建圖形plot(nearbrk)xlabel(時間);ylabel(幅值);title(頻率突變信號);figure(2); %采用Daudechies5小波對信號進(jìn)行處理d,a=wavedec(nearbrk,3,db5);a3=wrcoef(a,d,a,db5,3);d3=wrcoef(d,d,a,db5,3);d2=wrcoef(d,d,a,db5,2);d1=wrcoef(d,d,a,db5,1);subplot(411);plot(a3);ylabel(近似信號 a3);title(小波分解后示意圖);subplot(412);plot(d3);ylabel(細(xì)節(jié)信號 d3);subplot(413);plot(d2)