%設置信噪比
snr=4;
%設置隨機數初值
init=2055615866;
%產生原始信號xref和含白噪聲信號x
[xref,x]=wnoise(4,11,snr,init);
xref=xref(1:2000);
x=x(1:2000);
%用sym8小波進行三層分解并用heursure軟閾值進行小波系數閾值量化并消噪
xd=wden(x,'heursure','s','one',10,'sym8');
subplot(3,2,1);
plot(xref);
title('原始信號');
subplot(3,2,2);
plot(x);
title('含噪信號');
%用傅立葉變換進行消噪,
%對原始信號進行傅立葉變換
xxref=fft(xref);
xxref=abs(xxref);
%對含噪信號進行傅立葉變換
xx=fft(x);
absxx=abs(xx);
%畫出傅立葉變換后的頻譜圖
subplot(3,2,3);
plot(xxref);
title('原始信號的頻譜圖');
subplot(3,2,4);
plot(absxx);
title('含噪信號的頻譜圖');
%進行低通濾波
indd2=200:1800;
xx(indd2)=zeros(size(indd2));
%進行傅立葉逆變換
xden=ifft(xx);
%提取逆變換后的實數部分
xden=real(xden);
xden=abs(xden);
subplot(3,2,5);
plot(xd);
title('小波消噪后的信號');
subplot(3,2,6);
plot(xden);
title('傅立葉消噪后的信號');