fs=5000;
nfft=10240;
load a.txt;
xdata= a;
figure(1);
plot(xdata);
fbh=cwt(xdata,1:2:64,'cgau4','plot');
%采樣頻率
s=imag(fbh)
%5層小波包分解
T=wpdec(s,5,'db2');
%重構低頻信號
y10=wprcoef(T,[1,0]);
y11=wprcoef(T,[1,1]);
%e10 = wentropy(y10,'shannon')
%e11 = wentropy(y11,'shannon')
figure(2);
subplot(2,2,1);
plot(y10);
xlabel('第一層低頻系數時間 t/s');
ylabel('第一層低頻系數幅值 A/W');
Y10=abs(fft(y10,nfft));
subplot(2,2,2);
plot((0:nfft/2-1)/nfft*fs,Y10(1:nfft/2));
xlabel('第一層低頻系數頻率 f/Hz');
ylabel('第一層低頻系數功率譜 P/W');
subplot(2,2,3);
plot(y11);
xlabel('第一層高頻系數時間 t/s');
ylabel('第一層高頻系數幅值 A/W');
Y11=abs(fft(y11,nfft));
subplot(2,2,4);
plot((0:nfft/2-1)/nfft*fs,Y11(1:nfft/2));
xlabel('第一層高頻系數頻率 f/Hz');
ylabel('第一層高頻系數功率譜 P/W');
%------------------------------------------------------------------

%-----------------------------------------------------------------
%重構第二層各系數
y20=wprcoef(T,[2,0]);
y21=wprcoef(T,[2,1]);
y22=wprcoef(T,[2,2]);
y23=wprcoef(T,[2,3]);
figure(3);
subplot(4,2,1);
plot(y20);
xlabel('Y10分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y20=abs(fft(y20,nfft));
subplot(4,2,2);
plot((0:nfft/2-1)/nfft*fs,Y20(1:nfft/2));
xlabel('Y10分解的低頻系數頻率f/Hz');
ylabel('功率譜P/W');
subplot(4,2,3);
plot(y21);
xlabel('Y10分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y21=abs(fft(y21,nfft));
subplot(4,2,4);
plot((0:nfft/2-1)/nfft*fs,Y21(1:nfft/2));
xlabel('Y10分解的高頻系數頻率f/Hz');
ylabel('功率譜P/W');
subplot(4,2,5);
plot(y22);
xlabel('Y11分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y22=abs(fft(y22,nfft));
subplot(4,2,6);
plot((0:nfft/2-1)/nfft*fs,Y22(1:nfft/2));
xlabel('Y11分解的低頻系數頻率f/Hz');
ylabel('功率譜P/W');
subplot(4,2,7);
plot(y23);
xlabel('Y11分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y23=abs(fft(y23,nfft));
subplot(4,2,8);
plot((0:nfft/2-1)/nfft*fs,Y23(1:nfft/2));
xlabel('Y11分解的高頻系數頻率f/Hz');
ylabel('功率譜P/W');
%重構第三層各系數
y30=wprcoef(T,[3,0]);
y31=wprcoef(T,[3,1]);
y32=wprcoef(T,[3,2]);
y33=wprcoef(T,[3,3]);
y34=wprcoef(T,[3,4]);
y35=wprcoef(T,[3,5]);
y36=wprcoef(T,[3,6]);
y37=wprcoef(T,[3,7]);
figure(4);
subplot(9,2,1);
plot(y30);
xlabel('Y20分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y30=abs(fft(y30,nfft));
subplot(9,2,2);
plot((0:nfft/2-1)/nfft*fs,Y30(1:nfft/2));
xlabel('Y20分解的低頻系數頻率f/Hz');
ylabel('功率譜P/W');
subplot(9,2,3);
plot(y31);
xlabel('Y20分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y31=abs(fft(y31,nfft));
subplot(9,2,4);
plot((0:nfft/2-1)/nfft*fs,Y31(1:nfft/2));
xlabel('Y20分解的高頻系數頻率f/Hz');
ylabel('功率譜P/W');
subplot(9,2,5);
plot(y32);
xlabel('Y21分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y32=abs(fft(y32,nfft));
subplot(9,2,6);
plot((0:nfft/2-1)/nfft*fs,Y32(1:nfft/2));
xlabel('Y21分解的低頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,7);
plot(y33);
xlabel('Y21分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y33=abs(fft(y33,nfft));
subplot(9,2,8);
plot((0:nfft/2-1)/nfft*fs,Y33(1:nfft/2));
xlabel('Y21分解的高頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,9);
plot(y34);
xlabel('Y22分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y34=abs(fft(y34,nfft));
subplot(9,2,10);
plot((0:nfft/2-1)/nfft*fs,Y34(1:nfft/2));
xlabel('Y22分解的低頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,11);
plot(y35);
xlabel('Y22分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y35=abs(fft(y35,nfft));
subplot(9,2,12);
plot((0:nfft/2-1)/nfft*fs,Y35(1:nfft/2));
xlabel('Y22分解的高頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,13);
plot(y36);
xlabel('Y23分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y36=abs(fft(y36,nfft));
subplot(9,2,14);
plot((0:nfft/2-1)/nfft*fs,Y36(1:nfft/2));
xlabel('Y23分解的低頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,15);
plot(y37);
xlabel('Y23分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y37=abs(fft(y37,nfft));
subplot(9,2,16);
plot((0:nfft/2-1)/nfft*fs,Y37(1:nfft/2));
xlabel('Y23分解的低頻系數時間t/s');
ylabel('功率譜P/W');
y30=wprcoef(T,[3,0]);
y31=wprcoef(T,[3,1]);
y32=wprcoef(T,[3,2]);
y33=wprcoef(T,[3,3]);
y34=wprcoef(T,[3,4]);
y35=wprcoef(T,[3,5]);
y36=wprcoef(T,[3,6]);
y37=wprcoef(T,[3,7]);
figure(4);
subplot(9,2,1);
plot(y30);
xlabel('Y20分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y30=abs(fft(y30,nfft));
subplot(9,2,2);
plot((0:nfft/2-1)/nfft*fs,Y30(1:nfft/2));
xlabel('Y20分解的低頻系數頻率f/Hz');
ylabel('功率譜P/W');
subplot(9,2,3);
plot(y31);
xlabel('Y20分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y31=abs(fft(y31,nfft));
subplot(9,2,4);
plot((0:nfft/2-1)/nfft*fs,Y31(1:nfft/2));
xlabel('Y20分解的高頻系數頻率f/Hz');
ylabel('功率譜P/W');
subplot(9,2,5);
plot(y32);
xlabel('Y21分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y32=abs(fft(y32,nfft));
subplot(9,2,6);
plot((0:nfft/2-1)/nfft*fs,Y32(1:nfft/2));
xlabel('Y21分解的低頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,7);
plot(y33);
xlabel('Y21分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y33=abs(fft(y33,nfft));
subplot(9,2,8);
plot((0:nfft/2-1)/nfft*fs,Y33(1:nfft/2));
xlabel('Y21分解的高頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,9);
plot(y34);
xlabel('Y22分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y34=abs(fft(y34,nfft));
subplot(9,2,10);
plot((0:nfft/2-1)/nfft*fs,Y34(1:nfft/2));
xlabel('Y22分解的低頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,11);
plot(y35);
xlabel('Y22分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y35=abs(fft(y35,nfft));
subplot(9,2,12);
plot((0:nfft/2-1)/nfft*fs,Y35(1:nfft/2));
xlabel('Y22分解的高頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,13);
plot(y36);
xlabel('Y23分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y36=abs(fft(y36,nfft));
subplot(9,2,14);
plot((0:nfft/2-1)/nfft*fs,Y36(1:nfft/2));
xlabel('Y23分解的低頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,15);
plot(y37);
xlabel('Y23分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y37=abs(fft(y37,nfft));
subplot(9,2,16);
plot((0:nfft/2-1)/nfft*fs,Y37(1:nfft/2));
xlabel('Y23分解的低頻系數時間t/s');
ylabel('功率譜P/W');
y40=wprcoef(T,[4,0]);
y41=wprcoef(T,[4,1]);
y42=wprcoef(T,[4,2]);
y43=wprcoef(T,[4,3]);
y44=wprcoef(T,[4,4]);
y45=wprcoef(T,[4,5]);
y46=wprcoef(T,[4,6]);
y47=wprcoef(T,[4,7]);
y48=wprcoef(T,[4,8]);
y49=wprcoef(T,[4,9]);
y410=wprcoef(T,[4,10]);
y411=wprcoef(T,[4,11]);
y412=wprcoef(T,[4,12]);
y413=wprcoef(T,[4,13]);
y414=wprcoef(T,[4,14]);
y415=wprcoef(T,[4,15]);
figure(5);
subplot(9,2,1);
plot(y40);
xlabel('Y30分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y40=abs(fft(y40,nfft));
subplot(9,2,2);
plot((0:nfft/2-1)/nfft*fs,Y40(1:nfft/2));
xlabel('Y30分解的低頻系數頻率f/Hz');
ylabel('功率譜P/W');
subplot(9,2,3);
plot(y41);
xlabel('Y30分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y41=abs(fft(y41,nfft));
subplot(9,2,4);
plot((0:nfft/2-1)/nfft*fs,Y41(1:nfft/2));
xlabel('Y30分解的高頻系數頻率f/Hz');
ylabel('功率譜P/W');
subplot(9,2,5);
plot(y42);
xlabel('Y31分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y42=abs(fft(y42,nfft));
subplot(9,2,6);
plot((0:nfft/2-1)/nfft*fs,Y42(1:nfft/2));
xlabel('Y31分解的低頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,7);
plot(y43);
xlabel('Y31分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y43=abs(fft(y43,nfft));
subplot(9,2,8);
plot((0:nfft/2-1)/nfft*fs,Y43(1:nfft/2));
xlabel('Y31分解的高頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,9);
plot(y43);
xlabel('Y32分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y44=abs(fft(y44,nfft));
subplot(9,2,10);
plot((0:nfft/2-1)/nfft*fs,Y44(1:nfft/2));
xlabel('Y32分解的低頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,11);
plot(y45);
xlabel('Y32分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y45=abs(fft(y45,nfft));
subplot(9,2,12);
plot((0:nfft/2-1)/nfft*fs,Y45(1:nfft/2));
xlabel('Y32分解的高頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,13);
plot(y46);
xlabel('Y33分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y46=abs(fft(y46,nfft));
subplot(9,2,14);
plot((0:nfft/2-1)/nfft*fs,Y46(1:nfft/2));
xlabel('Y33分解的低頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,15);
plot(y47);
xlabel('Y33分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y47=abs(fft(y47,nfft));
subplot(9,2,16);
plot((0:nfft/2-1)/nfft*fs,Y47(1:nfft/2));
xlabel('Y33分解的低頻系數時間t/s');
ylabel('功率譜P/W');
figure(6);
subplot(9,2,1);
plot(y48);
xlabel('Y34分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y48=abs(fft(y48,nfft));
subplot(9,2,2);
plot((0:nfft/2-1)/nfft*fs,Y48(1:nfft/2));
xlabel('Y30分解的低頻系數頻率f/Hz');
ylabel('功率譜P/W');
subplot(9,2,3);
plot(y49);
xlabel('Y34分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y49=abs(fft(y49,nfft));
subplot(9,2,4);
plot((0:nfft/2-1)/nfft*fs,Y49(1:nfft/2));
xlabel('Y34分解的高頻系數頻率f/Hz');
ylabel('功率譜P/W');
subplot(9,2,5);
plot(y410);
xlabel('Y35分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y410=abs(fft(y410,nfft));
subplot(9,2,6);
plot((0:nfft/2-1)/nfft*fs,Y410(1:nfft/2));
xlabel('Y35分解的低頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,7);
plot(y411);
xlabel('Y35分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y411=abs(fft(y411,nfft));
subplot(9,2,8);
plot((0:nfft/2-1)/nfft*fs,Y411(1:nfft/2));
xlabel('Y36分解的高頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,9);
plot(y412);
xlabel('Y36分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y412=abs(fft(y44,nfft));
subplot(9,2,10);
plot((0:nfft/2-1)/nfft*fs,Y412(1:nfft/2));
xlabel('Y36分解的低頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,11);
plot(y413);
xlabel('Y36分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y413=abs(fft(y413,nfft));
subplot(9,2,12);
plot((0:nfft/2-1)/nfft*fs,Y413(1:nfft/2));
xlabel('Y36分解的高頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,13);
plot(y414);
xlabel('Y37分解的低頻系數時間t/s');
ylabel('幅值A/W');
Y414=abs(fft(y46,nfft));
subplot(9,2,14);
plot((0:nfft/2-1)/nfft*fs,Y414(1:nfft/2));
xlabel('Y37分解的低頻系數時間t/s');
ylabel('功率譜P/W');
subplot(9,2,15);
plot(y415);
xlabel('Y37分解的高頻系數時間t/s');
ylabel('幅值A/W');
Y415=abs(fft(y415,nfft));
subplot(9,2,16);
plot((0:nfft/2-1)/nfft*fs,Y415(1:nfft/2));
xlabel('Y37分解的低頻系數時間t/s');
ylabel('功率譜P/W');