小波变换

一．作业内容

（1） 任意给定一离散的1D数字信号（或实际采集1D信号），编成实现Harr小波的一尺度快速分解，画出原始信号及分解后近似分量及细节分量图。 

（2）利用上述分解得到的近似分量和细节分量，进行Harr小波反变换，即重构信号。画出重构信号的曲线。 

（3）计算重构信号与原始信号的残差，画出残差曲线示意图，并对残差曲线作进一步分析。

二．作业分析

1.任意给定一离散的1D数字信号（或实际采集1D信号），编成实现Harr小波的一尺度快速分解，画出原始信号及分解后近似分量及细节分量图。

[plain]  view plain  copy

1. %装载离散的1D数字信号  
2. load leleccum;  
3. s = leleccum(1:1000);%1000个采样点  
4. %用小波函数db1对信号进行单尺度小波分解  
5. [cA1,cD1]=dwt(s,'haar');  
6. subplot(3,2,1); plot(s);   
7. title('原始信号');  
8. subplot(3,2,2); plot(cA1);  
9. title('近似分量');  
10. subplot(3,2,3); plot(cD1);  
11. title('细节分量');  

2.利用上述分解得到的近似分量和细节分量，进行Harr小波反变换，即重构信号。画出重构信号的曲线。

[plain]  view plain  copy

1. s_1 = idwt(cA1,cD1,'haar');%重构信号  
2. subplot(3,2,4); plot(s_1);  
3. title('重构信号');  

3.计算重构信号与原始信号的残差，画出残差曲线示意图，并对残差曲线作进一步分析。

[plain]  view plain  copy

1. d = s_1 - s;%残差  
2. subplot(3,2,5); plot(d);  
3. title('残差');  

分析：Haar小波分解与重构的残差很小，证明重构很有效果。