特征检测小波包matlab程序与数据

小波包原创特征检测Matlab程序,包含试验测出来的数据，可以用于学习和毕业论文，直接下载可用

clear
clc
N=1400; 
% t=5e-9:5e-9:1400*5e-9;
t=0:5e-9:1399*5e-9;
xt1=（load（‘C:\Users\lenovo\Desktop\123.txt‘））‘;%数据文件所在的地址
xt2=（load（‘C:\Users\lenovo\Desktop\223.txt‘））‘;
xt3=（load（‘C:\Users\lenovo\Desktop\323.txt‘））‘;
xt4=（load（‘C:\Users\lenovo\Desktop\423.txt‘））‘;
xt5=（load（‘C:\Users\lenovo\Desktop\523.txt‘））‘;
figure（1）;
plot（txt1‘r‘）;
hold on
plot（txt2‘b‘）;
hold on
plot（txt3‘y‘）;
hold on
plot（txt4‘k‘）;
hold on
plot（txt5‘g‘）;
legend（‘正常‘‘表面一个气孔‘‘表面两个气孔‘‘下方有气孔‘‘表面裂痕‘）
title（‘信号图‘）; xlabel（‘时间（t）‘）; ylabel（‘幅值‘）; 
%***************************************************************
%进行小波包分解
T1=wpdec（xt14‘db4‘‘shannon‘）;
T2=wpdec（xt24‘db4‘‘shannon‘）;
T3=wpdec（xt34‘db4‘‘shannon‘）;  
T4=wpdec（xt44‘db4‘‘shannon‘）; 
T5=wpdec（xt54‘db4‘‘shannon‘）; 
figure（2）  %小波包树
plot（T1）
%最后一层时间频率图
wpviewcf（T11）;  
wpviewcf（T21）;
wpviewcf（T51）;
%重构小波包系数
w1=wprcoef（T1[40]）; 
w2=wprcoef（T1[41]）;
w3=wprcoef（T1[42]）;
w4=wprcoef（T1[43]）;
w5=wprcoef（T1[44]）;
w6=wprcoef（T1[45]）;
w7=wprcoef（T1[46]）;
w8=wprcoef（T1[47]）;
w9=wprcoef（T1[48]）;
w10=wprcoef（T1[49]）;
w11=wprcoef（T1[410]）;
w12=wprcoef（T1[411]）;
w13=wprcoef（T1[412]）;
w14=wprcoef（T1[413]）;
w15=wprcoef（T1[414]）;
w16=wprcoef（T1[415]）;
%***********************
w21=wprcoef（T2[40]）; 
w22=wprcoef（T2[41]）;
w23=wprcoef（T2[42]）;
w24=wprcoef（T2[43]）;
w25=wprcoef（T2[44]）;
w26=wprcoef（T2[45]）;
w27=wprcoef（T2[46]）;
w28=wprcoef（T2[47]）;
w29=wprcoef（T2[48]）;
w210=wprcoef（T2[49]）;
w211=wprcoef（T2[410]）;
w212=wprcoef（T2[411]）;
w213=wprcoef（T2[412]）;
w214=wprcoef（T2[413]）;
w215=wprcoef（T2[414]）;
w216=wprcoef（T2[415]）;
%***********************
w31=wprcoef（T3[40]）; 
w32=wprcoef（T3[41]）;
w33=wprcoef（T3[42]）;
w34=wprcoef（T3[43]）;
w35=wprcoef（T3[44]）;
w36=wprcoef（T3[45]）;
w37=wprcoef（T3[46]）;
w38=wprcoef（T3[47]）;
w39=wprcoef（T3[48]）;
w310=wprcoef（T3[49]）;
w311=wprcoef（T3[410]）;
w312=wprcoef（T3[411]）;
w313=wprcoef（T3[412]）;
w314=wprcoef（T3[413]）;
w315=wprcoef（T3[414]）;
w316=wprcoef（T3[415]）;
%***********************
w41=wprcoef（T4[40]）; 
w42=wprcoef（T4[41]）;
w43=wprcoef（T4[42]）;
w44=wprcoef（T4[43]）;
w45=wprcoef（T4[44]）;
w46=wprcoef（T4[45]）;
w47=wprcoef（T4[46]）;
w48=wprcoef（T4[47]）;
w49=wprcoef（T4[48]）;
w410=wprcoef（T4[49]）;
w411=wprcoef（T4[410]）;
w412=wprcoef（T4[411]）;
w413=wprcoef（T4[412]）;
w414=wprcoef（T4[413]）;
w415=wprcoef（T4[414]）;
w416=wprcoef（T4[415]）;
%***********************
w51=wprcoef（T5[40]）; 
w52=wprcoef（T5[41]）;
w53=wprcoef（T5[42]）;
w54=wprcoef（T5[43]）;
w55=wprcoef（T5[44]）;
w56=wprcoef（T5[45]）;
w57=wprcoef（T5[46]）;
w58=wprcoef（T5[47]）;
w59=wprcoef（T5[48]）;
w510=wprcoef（T5[49]）;
w511=wprcoef（T5[410]）;
w512=wprcoef（T5[411]）;
w513=wprcoef（T5[412]）;
w514=wprcoef（T5[413]）;
w515=wprcoef（T5[414]）;
w516=wprcoef（T5[415]）;
%计算第四层各个频段的能量

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件       14344  2018-02-20 22:34  123.txt
     文件       14326  2018-02-19 15:01  223.txt
     文件       14218  2018-02-27 10:37  323.txt
     文件       14234  2018-02-27 10:42  423.txt
     文件       14222  2018-02-27 10:42  523.txt
     文件        5926  2018-03-06 21:08  WaveletPacketTransform.m