noise1.m

一个基于MATLAB去噪的程序

load leleccum1;
indx=1:1024;
x=leleccum(indx);
%产生含噪信号
init=2055615866;
randn('seed',init);
nx=x+18*randn(size(x));
%获取消噪阈值
[thr,sorh,keepapp]=ddencmp('den','wv',nx);
%对噪音进行消噪
xd=wdencmp('gbl',nx,'haar',2,thr,sorh,keepapp);
subplot(2,2,1);
plot(x);
title('原始信号');
subplot(2,2,2);
plot(nx);
title('含噪信号');
subplot(2,2,3);
plot(xd);