小波学习之一（单层一维离散小波变换DWT的Mallat算法C++和MATLAB实现）

最新推荐文章于 2024-08-27 13:33:51 发布

v_hyx

最新推荐文章于 2024-08-27 13:33:51 发布

阅读量2.3w

点赞数 10

分类专栏：数字信号处理文章标签： c++ DWT Mallat算法 matlab 小波变换

本文链接：https://blog.csdn.net/share_happy_1984/article/details/8557071

版权

本文详细介绍了Mallat算法，包括离散序列的Mallat算法分解和重构公式，并在C++和MATLAB中实现了一维离散小波变换。内容涵盖了小波变换结果序列长度、滤波器组获取、信号边界延拓方法的比较，以及C++和MATLAB的代码实现。通过对db1小波函数的实例分析，验证了C++实现与MATLAB效果的一致性，强调了C++实现的效率优势。

摘要由CSDN通过智能技术生成

1 Mallat算法

离散序列的Mallat算法分解公式如下：

其中，H(n)、G(n)分别表示所选取的小波函数对应的低通和高通滤波器的抽头系数序列。

从Mallat算法的分解原理可知，分解后的序列就是原序列与滤波器序列的卷积再进行隔点抽取而来。

离散序列的Mallat算法重构公式如下：

其中，h(n)、g(n)分别表示所选取的小波函数对应的低通和高通滤波器的抽头系数序列。

2 小波变换实现过程(C/C++)

2.1 小波变换结果序列长度

小波的Mallat算法分解后的序列长度由原序列长SoureLen和滤波器长FilterLen决定。从Mallat算法的分解原理可知，分解后的序列就是原序列与滤波器序列的卷积再进行隔点抽取而来。即分解抽取的结果长度为(SoureLen+FilterLen-1)/2。

2.2 获取滤波器组

对于一些通用的小波函数，简单起见，可以通过Matlab的wfilters(‘wavename’)获取4个滤波器；特殊的小波函数需要自行构造获得。

下面以db1小波函数(Haar小波)为例，其变换与重构滤波器组的结果如下：

//matlab输入获取命令
>> [Lo_D,Hi_D,Lo_R,Hi_R] = wfilters('db1')

//获取的结果
Lo_D =
    0.7071    0.7071
Hi_D =
   -0.7071    0.7071
Lo_R =
    0.7071    0.7071
Hi_R =
    0.7071   -0.7071