%設置信噪比和隨機種子值
snr=4;
init=2055615866;
%產生原始信號sref和被高斯白噪聲污染的信號s
[sref,s]=wnoise(1,11,snr,init);
%用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('消噪信號');