正交小波滤波器系数和的合理性研究

时间:2022-09-16 05:55:54

正交小波滤波器系数和的合理性研究

【摘 要】正交小波滤波器系数是小波分析过程中十分重要的参数。然而,很多书藉及文献中对这个重要参数的总和却推导出了两种不同的结论。这两种结论必然有一个是不合理的。在对正交小波的滤波器系数总和的获得过程进行严格的、完整的数学推导基础上,得到了一个明显合理的结论。这一推导过程可以补充很多文献中关于正交小波分析的内容。并研究了有关文献中产生不合理结论的过程,分析了产生不合理结论的原因。

【关键词】正交小波变换;滤波器系数;Mallat算法

0 引言

多分辨分析理论提供了构造正交小波的一个基本框架。其基本方法是:先构造一个满足多分辨率的嵌套闭子空间序列{Vj}j∈z,并使它构成整个L2(R)空间。在V0空间找一个函数φ(t),它的平移函数序列{φ(t-k)}k∈z构成V0空间的Riesz基。特别指出,Riesz基仅能保证序列{φ(t-k)}k∈z线性无关,但并不能保证其构成V0空间的一个标准正交基。一般的,我们可以对一个非正交的函数φ(t)通过正交化方法使其正交。则当?准(t)∈V0时,序列{?准(t-k)}k∈z构成V0空间的一个标准正交基,?准j,k(t)=2j/2?准(2jt-k),j,k∈Z构成Vj的一个标准正交基。于是,我们说?准(t)生成正交多分辨分析。这样就可以通过?准(t)的二尺度方程求出相应的滤波器系数{hk},并根据?准(t)和{hk}采用Mallat算法对正交小波进行分解与重构。

通过上述分析,滤波器系数{hk}是正交小波分析过程中十分重要的参数。然而,很多书藉及文献中对这个重要参数的总和却推导出了不同的结论。如文献[1-2]中推导出滤波器系数总和hk=2;而文献[3-5]中却推导出hk=。这两种结论必然有一个是不合理的。这里仅用“不合理”而没有用“错误”来说明这一结论,是因为即使在得到不合理结论的文献中,并没有因为这个不合理结论而造成更多其它理论上的错误,原因在于这两个结论之间仅仅相差了一个的常倍数。但是,两种截然不同的结论往往给读者的阅读造成困惑,给学习者特别是初学者的理解带来不少困难。本文将对正交小波的滤波器系数总和这一结论的获得过程重新进行严格、完整的数学推导,从而得到一个合理的结论,以助学习者理解。之所以要重新完成该推导过程,是因为诸多文献中在涉及到这块内容时鲜有完整的推导和证明过程,也许正是这个原因疏忽了这个不合理结论的存在。

1 正交小波滤波器系数总和的合理结论推导过程

正交小波滤波器系数的总和需由正交小波分析的二尺度方程得到,所以推导过程需由二尺度方程的有关定义开始推导。本文中所涉及的有关多尺度逼近、二尺度方程及多分辨分析的相关定义及内容请参考相关文献。

由多分辨分析概念得知,如果?准(t)为尺度空间V0的一个标准正交基函数,其平移序列{?准(t-k)}k∈z构成V0空间的一个标准正交基,而对任意的j,k∈Z,?准j,k(t)=2j/2?准(2jt-k)构成Vj的一个标准正交基,则称?准(t)为尺度函数。又由于V0?奂V1,所以?准(t)也必然属于V1空间,也即?准(t)必然可用V1空间的正交基?准1,k(t)线性展开(这里的?准(t)可认为是尺度函数在0尺度空间下,平移量j=0时所得到的函数,故可用?准0,0(t)方式来表示其含义):

其中展开系数hk=,也称低通滤波器系数。由于式(1)描述的是相邻两个尺度空间基函数之间的关系,所次称此式为二尺度方程。

又由于Vj?奂Vj+1,?准j,0(t)∈Vj,所以?准j,0(t)也必然属于Vj+1空间,也即?准j,0(t)可用Vj+1空间的正交基?准j+1,k(t)线性展开,如下:

其中展开系数hk′=,且有

式(3)说明,正交小波的滤波器系数hk不随尺度j的变化而变化。即二尺度关系存在于任意相邻尺度j,j+1之间,那么式(1)与式(2)的表达式可统一为:

对式(4)两边同时对t积分得:

将式(6)代入式(5),则得hk=(7)

在式(1)两里边同时作Fourier变换,可得到二尺度方程的频域表示:

由式(8)得:

令ω=0,由式(9)可得,将此式代入式(10),可得:

式(1)到式(11)的数学推导过程是严密而完整的,因而得到的式(7)和式(11)的结论也应该是完全合理的。式(7)是从时域方面,式(11)是从频域方面,说明了正交小波的滤波器系数总和应为,在使用过程中应注意不要用错。

2 正交小波滤波器系数总和的不合理结论的产生原因

第1节中的数学推导过程分别从时域和频域两个方面得到了正交小波滤波器系数总和的结论,这也是各种文献得到该结论最常用的两种方法。而有些文献得到的不合理结论的原因也正是在这两种推导方法中产生的。

2.1 从时域方面得到不合理结论的过程

在多尺度逼近中,通常取?准j,k(t)=?准(2jt-k),特别是当?准j,k(t)具有内插基函数特点时,对构造各尺度空间Vj的近似函数f j十分便利。有的文献中直接采用该?准j,k(t)定义出二尺度方程形式:

注意,式(12)表示的二尺度方程与式(4)比较差了一个,而这种二尺度方程格式在不少文献中均有采用。式(12)既表示了尺度空间V0和V1的关系,也表示了任意尺度空间Vj和Vj+1之间的传递关系,即可由式(12)得到更通用的二尺度方程:

根据式(1)至式(7)的推导过程可知,对式(13)两边同时积分可得到滤波器系数总和应该为2,即hk=2。

2.2 从频域方面得到不合理结论的过程

从式(8)到式(12)的数学推导可知,对式(12)所表示的二尺度方程两边作Fourier变换,也能得到正交小波滤波器系数总和为2的结论,即hk=2。原因也在于式(12)的表达式与式(1)比较相差一个的倍数。

上一篇:高校就业指导工作微博客定位与日常建设的思考 下一篇:天津经济技术开发区“科技小巨人”发展现状及...