%装入变换放大器输入输出数据

%bf_150ms.dat为正常系统输出信号

%bf_160ms.dat为故障系统输出信号

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进行三层分解

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

%plontree(t)%画小波包树结构的图形

%下面对正常信号第三层各系数进行重构

%s130是指信号sl的[3,0]结点的重构系数;其他依次类推

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]);

%画出至构系数的波形

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');

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

%计算正常信号各重构系数的方差

%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构造的向量

disp=('正常信号的输出向量')

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

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

%用db1小波包对故障信号s2进行三层分解

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

%plottree(t)%画小波包树结构的图形

%s230是指信号S2的[3,0]结点的重构系数,其他以此类推

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]);

%画出重构系数的波形

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');

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

%计算故障信号各重构系数的方差

%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构造的向量

disp('故障信