noise3.m

来自「一个基于MATLAB去噪的程序」· M 代码 · 共 27 行

load leleccum1;
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,3);
alpha=2;
%获取小噪过程中的阈值
thr=wbmpen(c,l,sigma,alpha);
keepapp=1;
%对信号进行消噪
xd=wdencmp('gbl',c,l,'db6',3,thr,'s',keepapp);
subplot(2,2,1);
plot(x);
title('原始信号');
subplot(2,2,2);
plot(nx);
title('含噪信号');
subplot(2,2,3);
plot(xd);
title('消噪后的信号')