谐波小波转换(Harmonic Wavelet Transform)为学者大卫‧纽兰德(David E. Newland)于1993年所提出,是一个以小波为基底的线性转换,得以将讯号变换至时频域(Time-Frequency Domain)上。谐波小波转换结合了短时距傅立叶变换和连续小波转换两者之优点的讯号分析工具,而其离散版本则可以用快速傅立叶变换做有效率的运算。
定义与性质
基础推理
考量一个偶对称的实数函数
,其傅立叶变换定义为:
![{\displaystyle W_{e}(\omega )={\begin{cases}{\frac {1}{4\pi }}&{\mbox{for }}-4\pi \leq \omega <-2\pi \\{\frac {1}{4\pi }}&{\mbox{for }}\quad 2\pi \leq \omega <4\pi \\0&{\mbox{elsewhere}}\end{cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/87e4314a7b3e3e7990ae8def903ab941a33d6f4a)
则透过反傅立叶变换,我们可以得到该函数
为:
![{\displaystyle w_{e}(x)=\int _{-\infty }^{\infty }W_{e}(\omega )e^{i\omega x}dw={\frac {sin4\pi x-sin2\pi x}{2\pi x}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/721c79fe80c22d5799cf9c5455f3f854d59ef76e)
而考量另一奇对称的函数
,若定义其傅立叶变换为:
![{\displaystyle W_{e}(\omega )={\begin{cases}{\frac {i}{4\pi }}&{\mbox{for }}-4\pi \leq \omega <-2\pi \\{\frac {-i}{4\pi }}&{\mbox{for }}2\pi \leq \omega <4\pi \\0&{\mbox{elsewhere}}\end{cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c463f6d4b51faf4810051e1d6d993c76f23d798a)
则其反傅立叶变换会得到
为:
![{\displaystyle w_{o}(x)=\int _{-\infty }^{\infty }W_{o}(\omega )e^{i\omega x}dw={\frac {-(cos4\pi x-cos2\pi x)}{2\pi x}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4aed04f8a667511b05102f4c69ea68af4fa330af)
假如结合
和
,透过
的关系,我们会得到一复数函数,并定义它为谐波小波(Harmonic Wavelet)。本谐波小波将为以下数学形式:
![{\displaystyle w(x)={\frac {e^{i4\pi x}-e^{i2\pi x}}{i2\pi x}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0f135703dc0877a57841776e7a6720eec81f2f68)
也由于傅立叶转换的特性和
和
的定义,谐波小波的傅立叶转换对为:
![{\displaystyle W(\omega )=W_{e}(\omega )+iW_{o}(\omega )={\begin{cases}{\frac {1}{2\pi }}&{\mbox{for }}2\pi \leq \omega <4\pi \\0&{\mbox{elsewhere}}\end{cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ec7eab60c372619dbc375cc928e708e2d94cc82c)
一系列的谐波小波
接着,考量到小波转换中的精神--母小波的缩放(Dilation)和平移,透过伸张方程式(Dilation Equation)我们可以写出一系列的谐波小波(其中
和
皆为整数):
![{\displaystyle w(x)\Rightarrow w(2^{j}x-k)={\frac {e^{i4\pi (2^{j}x-k)}-e^{i2\pi (2^{j}x-k)}}{i2\pi (2^{j}-k)}}=v(x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cfe2d96e79117d81d3581e60084d0a36fac853d6)
根据前文对
的定义,或是透过直接计算傅立叶转换对,我们也可以得到缩放和平移后的一系列谐波小波在频域上的表示法:
![{\displaystyle V(\omega )={\frac {1}{2^{j}}}e^{\frac {-i\omega k}{2^{j}}}W({\frac {\omega }{2^{j}}})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/603dcd764ecf0f2da8aa4e410ea5a6befad34cf4)
而若我们将不同的正整数
带入上式,例如
和
,我们会发现后者的振幅会是前者的一半,然而其频带宽会是前者的两倍。这样的特性使得每一阶(Level,对应到不同的
)的谐波小波,其频域将随着阶数越高而越宽,由是达到多分辨率的效果。
低频频带(Zero-frequency band)
随着
的阶数比0越来越小,频带的振幅将越来越高、越来越窄,一路向频率为0的位置延伸。而根据多分辨率分析的理论,我们可以将这些阶数小于0的频带全部收为一个频带,并定义为-1阶(
)。它涵盖了DC到
的频带范围。以小波转换的术语来说,这样具低通滤波性质的函式,被称之为缩放函数(Scaling Function),又称为父小波(Father Wavelet)。谐波小波的缩放函数定义为:
,其频域特性将是一个介于
的方波,振幅为![{\displaystyle {\frac {1}{2\pi }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/129204d50704b07e6a4223870954242b21170354)
若要证明谐波小波有正交的特性,必须分两个层次讨论,
(不同阶的谐波小波)和
(不同位移量)。首先讨论不同阶的谐波小波。根据傅立叶理论,若两任意阶数的谐波小波正交,它将有下列关系(参考David Newland,1993):
![{\displaystyle \int _{-\infty }^{\infty }w(x)v(x)dx=\int _{-\infty }^{infty}W(\omega )V(-\omega )dx}](https://wikimedia.org/api/rest_v1/media/math/render/svg/84998e5077a9bab1e7788712650b4c1601b905f5)
因为任意阶数之谐波小波其频谱皆分布在正频率轴,故
永远为0。我们还必须证明下式也成立:
![{\displaystyle \int _{-\infty }^{infty}w(x)v^{*}(x)dx=\int _{-\infty }^{infty}W(\omega )V^{*}(\omega )dx}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0a848cc9174d17434d16d0787cb14a3e989e949f)
而因为不同阶数之谐波小波其频带不相交,故上式的右式也为0,由是证明不同阶数谐波小波的正交特性。至于同阶数、不同位移量的谐波小波,因为傅立叶变换的特性,在时域的位移相当于在频域的讯号必须乘上一个线性相位,因此对位移之谐波小波来说,必须满足下式:
![{\displaystyle \int _{2\pi }^{4\pi }e^{i\omega k}d\omega =0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3cdbf49cb0eec41663217f776c5151dde9a54c6a)
当k不为0的时候,上式将会成立。换言之,当具有位移存在时,谐波小波正交的特质成立。最后,我们也可以用相似的证明方式,证明谐波小波之父小波也具有正交特性。
谐波小波转换
如同傅立叶级数一般,由于父小波和母小波皆具有正交的特性,我们可以用它们做为基底来对一个函数做展开:
![{\displaystyle f(t)=\sum _{k=-\infty }^{\infty }\left[a_{k}\phi (t-k)+{\tilde {a}}_{k}\phi ^{*}(t-k)\right]+\sum _{j=0}^{\infty }\sum _{k=-\infty }^{\infty }\left[a_{j,k}w(2^{j}t-k)+{\tilde {a}}_{j,k}w^{*}(2^{j}t-k)\right].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/555700c6c4cfaa9449789125bc50d5f227ca4f59)
而它们的系数,根据投影,分别为:
![{\displaystyle {\begin{aligned}a_{j,k}&{}=2^{j}\int _{-\infty }^{\infty }f(t)\cdot w^{*}(2^{j}t-k)\,dt\\{\tilde {a}}_{j,k}&{}=2^{j}\int _{-\infty }^{\infty }f(t)\cdot w(2^{j}t-k)\,dt\\a_{k}&{}=\int _{-\infty }^{\infty }f(t)\cdot \phi ^{*}(t-k)\,dt\\{\tilde {a}}_{k}&{}=\int _{-\infty }^{\infty }f(t)\cdot \phi (t-k)\,dt.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/94c8d69edffd0d2c2c039068ffbb3f9a95a68646)
而对于实数函数
,系数将会有这样的关系:
、
。
参见
参考资料
- Newland, David. "Harmonic Wavelet Analysis". Proceedings of the Royal Society of London, Series A (Mathematical and Physical Sciences). Oct 1993, 443 (1917): 203 – 225.
- Lokenath Debnath. Wavelet Transforms and Their Applications. Boston: Birkhäuser. 2002: 475-490 [2013-01-17]. ISBN 0817642048.