noise2.m

一个基于MATLAB去噪的程序

load leleccum;
indx=1:1024;
x=leleccum(indx);
%产生含噪信号
init=2055615866;
randn('seed',init);
nx=x+18*randn(size(x));
%用小波函数‘db6’对信号进行3层分解
[c,l]=wavedec(nx,3,'db6');
%估计尺度1的噪声标准偏差
sigma=wnoisest(c,l,1);
alpha=2;
%获取消噪过程中的阈值
thr=wbmpen(c,l,sigma,alpha);
keepapp=1;
%将信号nx用小波函数sym5分解到第5层
%用minimaxi阈值选者对系数进行处理，消除噪声信号
lev=5;
xd=wden(nx,'rigrsure','s','sln',lev,'sym5');
subplot(2,2,1);
plot(x);

subplot(2,2,2);
plot(nx);

subplot(2,2,3);
plot(xd);