%設(shè)置信噪比和隨機(jī)種子值
snr=4;
init=2055615866;
%產(chǎn)生原始信號(hào)sref和被高斯白噪聲污染的信號(hào)s
[sref,s]=wnoise(1,11,snr,init);

%用sym8小波對(duì)信號(hào)s進(jìn)行3層分解并對(duì)細(xì)節(jié)系數(shù)
%選用sure閾值模式和尺度噪聲
xd=wden(s,'heursure','s','one',3,'sym8');

%顯示信號(hào)波形
subplot(3,1,1);
plot(sref);
xlabel('樣本序號(hào) n');
ylabel('幅值 A');
title('參考信號(hào)');
subplot(3,1,2);
plot(s);
title('含噪信號(hào)');
xlabel('樣本序號(hào) n');
ylabel('幅值 A');
subplot(3,1,3);
plot(xd);
title('消噪信號(hào)');
xlabel('樣本序號(hào) n');
ylabel('幅值 A');