load leleccum1;
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,'haar',2,thr,sorh,keepapp);
subplot(2,2,1);
plot(x);
title('原始信號');
subplot(2,2,2);
plot(nx);
title('含噪信號');
subplot(2,2,3);
plot(xd);