wavelet4.m

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

data=textread('data.txt');
% load leleccum;
indx=1:2000;
x=data(indx);
% init=20556158; 
% randn('seed',init); 
% x = x+ 18*randn(size(x));
figure(1);
subplot(2,1,1);
plot(x);
[thr,sorh,keepapp]=ddencmp('den','wv',data);
xd=wdencmp('gbl',x,'db4',2,thr,sorh,keepapp);

title('1');
subplot(2,1,2);
plot(xd);
title('2');