单尺度一维离散小波变换DWT的Mallat算法C++和MATLAB实现.docxVIP

小波学习之一（单尺度一维离散小波变换DWT的Mallat算法C++和MATLAB实现） 来源：/v_hyx/article/details/8557071Mallat算法小波变换实现过程CC小波变换结果序列长度获取滤波器组信号边界延拓小波分解小波重构c实现结果小波变换实现MATLAB小结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小波)为例，其变换与重构滤波器组的结果如下：[cpp]view plaincopyprint?//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 //matlab输入获取命令 [Lo_D,Hi_D,Lo_R,Hi_R] = wfilters(db1)//获取的结果Lo_D = 0.7071 0.7071Hi_D = -0.7071 0.7071Lo_R = 0.7071 0.7071Hi_R = 0.7071 -0.70712.3 信号边界延拓在Mallat算法中，假定输入序列是无限长的，而实际应用中输入的信号是有限的采样序列，这就会出现信号边界处理问题。对于边界信号的延拓一般有3种方法，即零延拓、对称延拓和周期延拓。3种延拓方法比较情况如下：对于正交小波变换来说，前两种延拓方法实现起来比较简单，但重建时会产生边界效应，而且分解的层数越多，产生的边界效应越显著。零延拓方法给人一种跳跃的感觉。至于对称性延拓，由于正交小波滤波器一般都是非对称性的（Haar小波基虽然是正交的，但它是非连续的），重建图象给人一种错位的感觉。相比较而言，只有最后一种延拓方式可以得到比较精确的重建结果，它不仅能保证分解与重建正确计算，而且恢复的质量也好。不过，周期性延拓方法虽然是常用的三种方法中比较好的方法，但会导致信号边缘的非连续性，从而会使得较高频率（子带）层的小波系数很大，即使信号本身相当平滑。从信号压缩的角度看，大的系数是希望避免的。信号的对称延拓可避免边缘的非连续性问题。然而，对称延拓只能和对称的小波滤波器一起适用。如果降低正交性要求，选择双正交小波变换，对称性延拓不失为一种好的方法。周期延拓可适用于任何小波变换，但可能导致输入序列边缘的不连续，使得高频系数较大。而对称延拓则避免了输入序列边界的不连续，是当前广泛采用的一种延拓方法。下式中给出了最常用的对称延拓表达式。当原序列长sLEN为偶数时延拓后的序列长为sLEN+2*(filterLEN)，而原序列长为奇数时则需要在右端再延拓一个元素。注：在Matlab中默认使用了对称延拓。2.4 小波分解在db1小波函数下，离散序列的Mallat算法分解公式展开如下：其它的db小波，不再详述。小波分解C++源码如下。[cpp]view plaincopyprint?/** * @brief 小波变换之分解* @paramsourceData源数据* @paramdataLen源数据长度* @paramdb过滤器类型* @paramcA分解后的近似部分序列-低频部分* @paramcD分解后的细节部分序列-高频部分* @return * 正常则返回分解后序列的数据长度，错误则返回-1 */ int Wavelet::Dwt(double *sourceData, intdataLen, Filter db, double *cA, double *cD) { if(dataLen 2) return -1; if((NULL == sourceData)||(NULL ==