小波去噪:
load leleccum;
index = 1:1024;
x = leleccum(index);
%产生噪声信号
init = 2055615866;
randn(‘seed’,init);
nx = x + 18*randn(size(x));
%获取消噪的阈值
[thr,sorh,keepapp] = ddencmp(‘den’,‘wv’,nx);
%对信号进行消噪
xd1 = wdencmp(‘gbl’,nx,‘db4’,1,thr,sorh,keepapp);
xd2 = wdencmp(‘gbl’,nx,‘db4’,2,thr,sorh,keepapp);
xd3 = wdencmp(‘gbl’,nx,‘db4’,3,thr,sorh,keepapp);
subplot(231);
plot(x);
title(‘原始信号’);
subplot(232);
plot(nx);
title(‘含噪信号’);
subplot(233);
plot(xd1);
title(‘阈值1消噪后的信号’);
subplot(234);
plot(xd2);
title(‘阈值2消噪后的信号’);
subplot(235);
plot(xd3);
title(‘阈值3消噪后的信号’);
均值滤波:
load leleccum;
index = 1:1024;
x = leleccum(index);
%产生噪声信号
init = 2055615866;
randn(‘seed’,init);
nx = x + 18*randn(size(x));
xd1=smooth(nx,3);
xd2=smooth(nx,5);
xd3=smooth(nx,7);
subplot(231);
plot(x);
title(‘原始信号’);
subplot(232);
plot(nx);
title(‘含噪信号’);
subplot(233);
plot(xd1);
title(‘窗口为3消噪后的信号’);
subplot(234);
plot(xd2);
title(‘窗口为5消噪后的信号’);
subplot(235);
plot(xd3);
title(‘窗口为7消噪后的信号’);
中值滤波:
load leleccum;
index = 1:1024;
x = leleccum(index);
%产生噪声信号
init = 2055615866;
randn(‘seed’,init);
nx = x + 18*randn(size(x));
xd1=medfilt1(nx,3);
xd2=medfilt1(nx,5);
xd3=medfilt1(nx,7);
subplot(231);
plot(x);
title(‘原始信号’);
subplot(232);
plot(nx);
title(‘含噪信号’);
subplot(233);
plot(xd1);
title(‘窗口为3消噪后的信号’);
subplot(234);
plot(xd2);
title(‘窗口为5消噪后的信号’);
subplot(235);
plot(xd3);
title(‘窗口为7消噪后的信号’);