% 當前延拓模式是補零
% 載入信號
load noisdopp; x = noisdopp;
figure(1); subplot(211); 
plot(x); title('原始信號');
% 利用db1小波包在第3層對x進行分解，使用Shannon熵
t = wpdec(x,3,'db1','shannon');
% 畫出小波包樹
plot(t)
% 重構小波包結點(2,1)
rcfs = wprcoef(t,[2 1]);
figure(1); subplot(212); 
plot(rcfs); title('重構小波包結點(2,1)');