load leleccum;
indx = 1:1024; 
x = leleccum(indx);
% 產生含噪信號
init=2055615866; 
randn('seed',init); 
nx = x + 18*randn(size(x)); 
[c,l]=wavedec(nx,3,'db5');
%設置尺度向量
n=[1,2,3];
%設置閥值向量
p=[120,110,100];
%對高頻系數進行閥值處理
nc=wthcoef('d',c,l,n,p);
%對修正后的小波分解結構進行重構
rx = waverec(nc,l,'db5');
subplot(2,2,1);
plot(x);
title('原始信號');
subplot(2,2,2);
plot(nx);
title('含噪信號');
subplot(2,2,3)
plot(rx);
title('消噪后的信號');