%三层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