%裝入變換放大器輸入輸出數(shù)據(jù)

%bf_150ms.dat為正常系統(tǒng)輸出信號

%bf_160ms.dat為故障系統(tǒng)輸出信號

load bf_150ms.dat;

load bf_160ms.dat;

s1＝bf_150ms(1:1000);%s1為正常信號

s2＝bf_160ms(1:1000);%s2為故障信號

%畫出正常信號與故障信號的原始波形

tittle(“原始信號’);

Ylabel('s1');

subplot(922); plot(s2);

title('故障信號');

Ylabel('s2');

%============================================

%用dbl小波包對正常信號s1進(jìn)行三層分解

[t,d]＝wpdec(sl,3,'db','shannon');

%plontree(t)%畫小波包樹結(jié)構(gòu)的圖形

%下面對正常信號第三層各系數(shù)進(jìn)行重構(gòu)

%s130是指信號sl的[3,0]結(jié)點的重構(gòu)系數(shù);其他依次類推

sl30＝wprcoef(t,d,[3,0]);

s13l＝wprcoef(t,d,[3,1]);

s132＝wprcoef(t,d,[3,2]);

sl33＝wprcoef(t,d,[3,3]);

sl34＝wprcoef(t,d,[3,4]);

s135＝wprcoef(t,d,[3,5]);

s136＝wprcoef(t,d,[3,6]);

s137＝wprcoef(t,d,[3,7]);

%畫出至構(gòu)系數(shù)的波形

subplot(9,2,3); plot(s130);

Ylabel('S130');

subpolt(9,2,5); plot(s131);

Ylabel('S13l');

subplot(9,2,7); plot(s132);

Ylabel('S132');

subplot(9,2,9); plot(s133);

Ylabel('S133');

subplot(9,2,11);plot(s134);

Ylabel('S134');

subplot(9,2,13);plot(s135);

Ylabel('S135');

subplot(9,2,15);plot(s136);

Ylabel('S136');

subplot(9,2,17);plot(s137);

Ylabel('S137');

%--------------------------------------

%計算正常信號各重構(gòu)系數(shù)的方差

%s10是指s130的方差，其他依此類推

s10＝norm(sl30);

sll＝norm(s131);

s12＝norm(sl32);

s13＝norm(sl33);

sl4＝norm(s134);

s15＝norm(s135);

s16＝norm(sl36);

s17＝norm(sl37);

%向量ssl是針對信號s1構(gòu)造的向量

disp=('正常信號的輸出向量')

ssl＝[sl0,s11,sl2,sl3,s14,s15,sl6,s17]

%＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

%用db1小波包對故障信號s2進(jìn)行三層分解

[t,d]=wpdec(s2,3,'db1','shannon');

%plottree(t)%畫小波包樹結(jié)構(gòu)的圖形

%s230是指信號S2的[3,0]結(jié)點的重構(gòu)系數(shù),其他以此類推

s230=wprcoef(t,d,[3,0]);

s231=wprcoef(t,d,[3,1]);

s232=wprcoef(t,d,[3,2]);

s233=wprcoef(t,d,[3,3]);

s234=wprcoef(t,d,[3,4]);

s235=wprcoef(t,d,[3,5]);

s236=wprcoef(t,d,[3,6]);

s237=wprcoef(t,d,[3,7]);

%畫出重構(gòu)系數(shù)的波形

subplot(9,2,4);plot(s230);

Ylabel('S230');

subplot(9,2,6);plot(s231);

Ylabel('S231');

subplot(9,2,8);plot(s232);

Ylabel('S232');

subplot(9,2,10);plot(s233);

Ylabel('S233');

subplot(9,2,12);plot(s234);

Ylabel('S234');

subplot(9,2,14);plot(s235);

Ylabel('S235');

subplot(9,2,16);plot(s236);

Ylabel('S236');

subplot(9,2,18);plot(s237);

Ylabel('S237');

%----------------------------------------------------------

%計算故障信號各重構(gòu)系數(shù)的方差

%s20是指s230的方差，其他依次類推

s20＝norm(s230);

s21＝norm(s231);

s22＝norm(s232);

s23＝norm(s233);

s24＝norm(s234);

s25＝norm(s235);

s26＝norm(s236);

s27＝norm(s237);

%向量ss2是針對信號S1構(gòu)造的向量

disp('故障信

總結(jié)

以上是生活随笔為你收集整理的matlab小波分析特征提取,基于小波分析的特征信号提取的matlab程序的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。