% 裝載噪聲chirp信號
load noischir; x = noischir;
% 使用sym6對信號進行5層小波包分解
wname = 'sym6'; lev = 5;
tree = wpdec(x,lev,wname);
% 由第1層高頻系數估計噪聲的標準差，結點索引號為2
det1 = wpcoef(tree,2);
sigma = median(abs(det1))/0.6745;
% 使用wpbmpen選擇閾值
alpha = 2;
thr = wpbmpen(tree,sigma,alpha)
% 上面的閾值，軟閾值處理，保存低頻，使用wpdencmp進行信號降噪
keepapp = 1;
xd = wpdencmp(tree,'s','nobest',thr,keepapp);
% 畫出原始信號和降噪后的信號
figure(1)
subplot(211), plot(x),
title('原始信號')
subplot(212), plot(xd)
title('降噪后信號')