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