wavelet5.m

来自「小波在信号处理中的应用」· M 代码 · 共 21 行

data=textread('data.txt');
indx=1:1024;
x=data(indx);
% x=noisbump;
% load noisbump;
wname = 'db6'; lev = 1;

[c,l]=wavedec(x,lev,wname);
sigma=wnoisest(c,l,1);
alpha=2;
thr=wbmpen(c',l',sigma,alpha)
% thr =2.7681
keepapp=1;
xd = wdencmp('gbl',c,l,wname,lev,thr,'s',keepapp);
% xd = wdencmp('gbl',c,l,wname,lev,thr,'s',keepapp);
figure(1);
subplot(211);
plot(x);
subplot(212);
plot(xd);