%装入变换放大器输入输出数据
%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]);