% 當前延拓模式是補零
% 裝載信號
load noisdopp; x = noisdopp;
figure(1); subplot(211); 
plot(x); title('原始信號');
% 使用db1小波包對x進行3層分解，使用shannon熵
wpt = wpdec(x,3,'db1');
% 畫出小波包樹
plot(wpt)
% 讀取小波包(2,1)的系數
cfs = wpcoef(wpt,[2 1]);
figure(1); subplot(212); 
plot(cfs); title('小波包(2,1)的系數');