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