matlab学习之降噪平滑算法

平滑降噪测试，代码如下

% 平滑降噪

% FFT变换和小波变换

clc

clf

clear

length_of_sig=128;

x=linspace(0,2*pi,length_of_sig);

% signal=5*sin(x)+2*sin(5*x)+randn(x);书上的这个是错的,随机数里的参数要求是整数

raw=5*sin(x)+2*sin(5*x);

signal=5*sin(x)+2*sin(5*x)+randn(1,length_of_sig);

jiequ=16;

transf=fft(signal);

filter_transf(1:jiequ)=transf(1:jiequ);

filter_transf(length_of_sig-jiequ:length_of_sig)=transf(length_of_sig-jiequ:length_of_sig);

filter_signal=ifft(filter_transf);

% before_hist=real(transf);

% after_hist=hist(before_hist);

% after_hist=histogram(before_hist);

[xd1,cxd1,lxd1] = wden(signal,'sqtwolog','h','sln',3,'db3');

[xd2,cxd2,lxd2] = wden(signal,'sqtwolog','s','one',2,'db3');

plot(x,signal,'k--')

hold on;

% plot(x,raw,'r')'无噪声信号',

plot(x,filter_signal,'b')

plot(x,xd1);

plot(x,xd2);

legend('无噪声信号+噪声','FFT恢复信号',...

    '小波1恢复信号','小波2恢复信号')

% 从结果看来小波2和FFT效果都不错。

　　

参考：

1、matlab教程，例1.6 傅里叶变换与信号分析，之后发现改教程是基于matlab5.2的。

2、七中滤波方法测试matlab实现（转）