% 產(chǎn)生heavy sine信號(hào)和噪聲
init = 2055615866;
[xref,x] = wnoise(5,11,7,init);
xref=xref(1:1000);
x=x(1:1000);
% 畫出原始信號(hào)和噪聲信號(hào)
subplot(221), plot(x),
title('原始信號(hào)')
subplot(222), plot(xref),
title('噪聲信號(hào)')
% 使用wpdencmp函數(shù)進(jìn)行信號(hào)降噪
n = length(x); 
thr = sqrt(2*log(n*log(n)/log(2))); 
xwpd = wpdencmp(x,'s',4,'sym4','sure',thr,1);
% 畫出小波包降噪信號(hào)
subplot(223), plot(xwpd),
title('小波包降噪信號(hào)')
% 比較小波基降噪結(jié)果
xwd = wden(x,'rigrsure','s','one',4,'sym4');
% 畫出小波降噪信號(hào)
subplot(224), plot(xwd),
title('小波降噪信號(hào)')