wf.m

该程序基于小波变换

%设置信噪比
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('傅立叶消噪后的信号');