《基于小波分析的特征信号提取的matlab程序.docx》由会员分享，可在线阅读，更多相关《基于小波分析的特征信号提取的matlab程序.docx（6页珍藏版）》请在三一办公上搜索。
1、基于小波分析的特征信号提取的matlab程序%装入变换放大器输入输出数据%bf_150ms.dat为正常系统输出信号%bf_160ms.dat为故障系统输出信号load bf_150ms.dat; load bf_160ms.dat;s1bf_150ms(1:1000);%s1为正常信号s2bf_160ms(1:1000);%s2为故障信号%画出正常信号与故障信号的原始波形tittle(“原始信号);Ylabel(s1);subplot(922); plot(s2);title(故障信号);Ylabel(s2);%=%用dbl小波包对正常信号s1进行三层分解t,dwpdec(sl,3,db,s
2、hannon);%plontree(t)%画小波包树结构的图形%下面对正常信号第三层各系数进行重构%s130是指信号sl的3,0结点的重构系数;其他依次类推sl30wprcoef(t,d,3,0);s13lwprcoef(t,d,3,1);s132wprcoef(t,d,3,2);sl33wprcoef(t,d,3,3);sl34wprcoef(t,d,3,4);s135wprcoef(t,d,3,5);s136wprcoef(t,d,3,6);s137wprcoef(t,d,3,7);%画出至构系数的波形subplot(9,2,3); plot(s130);Ylabel(S130);subp
3、olt(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是指s
4、130的方差，其他依此类推s10norm(sl30);sllnorm(s131);s12norm(sl32);s13norm(sl33);sl4norm(s134);s15norm(s135);s16norm(sl36);s17norm(sl37);%向量ssl是针对信号s1构造的向量disp=(正常信号的输出向量)sslsl0,s11,sl2,sl3,s14,s15,sl6,s17%用db1小波包对故障信号s2进行三层分解t,d=wpdec(s2,3,db1,shannon);%plottree(t)%画小波包树结构的图形%s230是指信号S2的3,0结点的重构系数,其他以此类推s230=w
5、prcoef(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);sub
6、plot(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的方差，其他依次类推s20norm(s230);s21norm(s231);s22norm(s232);s23norm(s233);s24norm(s234);s25norm(s235);s26norm(s236);s27norm(s237);%向量ss2是针对信号S1构造的向量disp(故障信号的输出向量)ss2=s20,s21,s22,s23,s24,s25,s26,s27
