DWT小波变换根据系数重建原始信号

%三层DWT分解，而后根据第三层近似分量重构将噪信号

data = importdata('C:\Users\DELL\Desktop\spot.txt');

data = data(:,4);

len = length(data);

[A1,D1] = dwt(data,'Haar'); %第一次DWT

[A2,D2] = dwt(A1,'Haar');     %第二次DWT

[A3,D3] = dwt(A2,'Haar');     %第三次DWT

figure
plot(A1);
title('一层近似分量');

figure
plot(D1);
title('一层细节分量');

figure
plot(A2);
title('二层近似分量');

figure
plot(D2);
title('二层细节分量');

figure
plot(A3);
title('三层近似分量');

figure
plot(D3);
title('三层细节分量');

%重构
Denoise = upcoef('a', A3, 'haar',3,len); % 根据第三层的近似分量重构降噪信号
figure

plot(Denoise);
title('降噪后的信号');

结果如下：

                   原始信号

最后一张为根据第三层近似分量重构的降噪信号。

另附一个小波工具箱博客：http://blog.sina.com.cn/s/blog_6163bdeb0102dw7a.html