clc;
%clear;
% 裝載噪聲chirp信號
%load noischir; x = noischir;
% 使用sym6對信號進行5層小波包分解
x=s;
wname = 'db3'; lev = 3;
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('降噪后信號')