% 當前的延拓模式是補零
% 裝載信號
load noisdopp; 
x = noisdopp;
% 使用默認熵（shannon熵），用db1小波對x進行3層分解
wpt = wpdec(x,3,'db1'); 
% 分解小波包[3 0]
wpt = wpsplt(wpt,[3 0]);
% 畫出小波包樹
plot(wpt)
% 計算最優小波包樹
blt = bestlevt(wpt);
% 畫出最優小波包樹
plot(blt)