%設置信噪比和隨機種子值
snr=4;
init=2055615866;
%=============================
%產生原始信號sref和被高斯白噪聲污染的信號s
[sref,s]=wnoise(1,11,snr,init);
sref=sref(1:2500);
s=s(1:1000)
%=============================
%用sym8小波對信號s進行3層分解并對細節系數
%選用SURE閾值模式和尺度噪聲
xd=wden(s,'heursure','s','one',3,'sym8');
%=============================
%對上述信號進行圖示
subplot(3,1,1);
plot(sref);
title('參考信號');
subplot(3,1,2);
plot(s);
title('染噪信號');
subplot(3,1,3);
plot(xd);
title('消噪信號');