load leleccum;
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,'db4',2,thr,sorh,keepapp);
subplot(2,2,1);
plot(x)
subplot(2,2,2);
plot(nx);
subplot(2,2,3);
plot(xd);
title('消噪后的信號')