%装入变换放大器输入输出数据
%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('故障信
号的输出向量')
ss2=[s20,s21,s22,s23,s24,s25,s26,s27]