基于小波分析的特征信号提取的matlab程序

%装入变换放大器输入输出数据%bf_150ms.dat为正常系统输出信号%bf_160ms.dat为故障系统输出信号loadbf_150ms.dat;loadbf_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]