%調頻信號 
for i=1:512
    x(i)=sin(500*pi*i/1000+20*sin(20*pi*i/1000));
end
figure(1);
plot(1:length(x),x);
% 3層Haar小波包
wpt = wpdec(x,3,'haar');
%最優基
T=besttree(wpt);
% 小波包結構
plot(wpt);
%Haar小波基下的相平面
wpviewcf(wpt,1);