谐波小波转换

✍ dations ◷ 2025-06-20 20:14:38 #信号处理

谐波小波转换(Harmonic Wavelet Transform)为学者大卫‧纽兰德(David E. Newland)于1993年所提出,是一个以小波为基底的线性转换,得以将讯号变换至时频域(Time-Frequency Domain)上。谐波小波转换结合了短时距傅立叶变换和连续小波转换两者之优点的讯号分析工具,而其离散版本则可以用快速傅立叶变换做有效率的运算。

考量一个偶对称的实数函数 w e ( x ) {\displaystyle w_{e}(x)} ,其傅立叶变换定义为:

则透过反傅立叶变换,我们可以得到该函数 w e ( x ) {\displaystyle w_{e}(x)} 为:

而考量另一奇对称的函数 w o ( x ) {\displaystyle w_{o}(x)} ,若定义其傅立叶变换为:

则其反傅立叶变换会得到 w o ( x ) {\displaystyle w_{o}(x)} 为:

假如结合 w e ( x ) {\displaystyle w_{e}(x)} w o ( x ) {\displaystyle w_{o}(x)} ,透过 w ( x ) = w e ( x ) + i w o ( x ) {\displaystyle w(x)=w_{e}(x)+iw_{o}(x)} 的关系,我们会得到一复数函数,并定义它为谐波小波(Harmonic Wavelet)。本谐波小波将为以下数学形式:

也由于傅立叶转换的特性和 W e ( ω ) {\displaystyle W_{e}(\omega )} W o ( ω ) {\displaystyle W_{o}(\omega )} 的定义,谐波小波的傅立叶转换对为:

接着,考量到小波转换中的精神--母小波的缩放(Dilation)和平移,透过伸张方程式(Dilation Equation)我们可以写出一系列的谐波小波(其中 j {\displaystyle j} k {\displaystyle k} 皆为整数):

根据前文对 W ( ω ) {\displaystyle W(\omega )} 的定义,或是透过直接计算傅立叶转换对,我们也可以得到缩放和平移后的一系列谐波小波在频域上的表示法:

而若我们将不同的正整数 j {\displaystyle j} 带入上式,例如 j = 3 {\displaystyle j=3} j = 4 {\displaystyle j=4} ,我们会发现后者的振幅会是前者的一半,然而其频带宽会是前者的两倍。这样的特性使得每一阶(Level,对应到不同的 j {\displaystyle j} )的谐波小波,其频域将随着阶数越高而越宽,由是达到多分辨率的效果。

随着 j {\displaystyle j} 的阶数比0越来越小,频带的振幅将越来越高、越来越窄,一路向频率为0的位置延伸。而根据多分辨率分析的理论,我们可以将这些阶数小于0的频带全部收为一个频带,并定义为-1阶( j = 1 {\displaystyle j=-1} )。它涵盖了DC到 2 π {\displaystyle 2\pi } 的频带范围。以小波转换的术语来说,这样具低通滤波性质的函式,被称之为缩放函数(Scaling Function),又称为父小波(Father Wavelet)。谐波小波的缩放函数定义为:

若要证明谐波小波有正交的特性,必须分两个层次讨论, j {\displaystyle j} (不同阶的谐波小波)和 k {\displaystyle k} (不同位移量)。首先讨论不同阶的谐波小波。根据傅立叶理论,若两任意阶数的谐波小波正交,它将有下列关系(参考David Newland,1993):

因为任意阶数之谐波小波其频谱皆分布在正频率轴,故 W ( ω ) V ( ω ) {\displaystyle W(\omega )V(-\omega )} 永远为0。我们还必须证明下式也成立:

而因为不同阶数之谐波小波其频带不相交,故上式的右式也为0,由是证明不同阶数谐波小波的正交特性。至于同阶数、不同位移量的谐波小波,因为傅立叶变换的特性,在时域的位移相当于在频域的讯号必须乘上一个线性相位,因此对位移之谐波小波来说,必须满足下式:

当k不为0的时候,上式将会成立。换言之,当具有位移存在时,谐波小波正交的特质成立。最后,我们也可以用相似的证明方式,证明谐波小波之父小波也具有正交特性。

如同傅立叶级数一般,由于父小波和母小波皆具有正交的特性,我们可以用它们做为基底来对一个函数做展开:

而它们的系数,根据投影,分别为:

而对于实数函数 f ( t ) {\displaystyle f(t)} ,系数将会有这样的关系: a ~ j , k = a j , k {\displaystyle {\tilde {a}}_{j,k}=a_{j,k}^{*}} a ~ k = a k {\displaystyle {\tilde {a}}_{k}=a_{k}^{*}}

相关

  • 抗原呈递交叉呈递是特定抗原呈现细胞吞噬并利用MHC I呈现外来抗原给细胞毒性T细胞的能力。交叉致敏,是交叉呈递后的结果,其描述的是透过交成呈递而使初始T细胞(英语:Naive T cell)变成活
  • 红外光谱法红外光谱学是光谱学中研究电磁光谱红外部分的分支。它包括了许多技术,到目前为止最常用的是吸收光谱学。同所有的分光镜技术一样,它可以被用来鉴别一种化合物和研究样品的成分
  • 福州船政局福州船政局,又名马尾船政局,1866年(清同治五年)8月19日由左宗棠在福州马尾设立。9月左宗棠改任陕甘总督,推荐由沈葆桢任总理船政大臣执掌船政局事务。初期委任法国人为监督,并且雇
  • 北京体育大学坐标:40°01′48″N 116°19′09″E / 40.0299048°N 116.3192935°E / 40.0299048; 116.3192935北京体育大学是位于中华人民共和国北京市的一所公立体育大学,主管部门为国家
  • 埃德蒙·马斯基埃德蒙·西克斯图斯·马斯基(英语:Edmund Sixtus Muskie,1914年3月28日-1996年3月26日),美国政治家,曾任缅因州州长和第58任美国国务卿。
  • 核酸序列核酸序列(英语:Nucleic acid sequence,亦称为核酸的一级结构)使用一串字母表示的真实的或者假设的携带基因信息的DNA分子的一级结构。每个字母代表一种核碱基,两个碱基形成一个碱
  • 分子机器分子机器(molecular machine)也称为纳米机器(nanomachine)是由少量的分子所组成,可以对特定的刺激(输入)产生准机械运动(输出)的物体。分子机器一词也常用来表示模仿巨观世界机器功能
  • 国际文字设计协会国际文字设计协会(法语:Association Typographique Internationale,简称:ATypI)是字体排印和字体设计行业的非营利性国际组织,于1957年由法国Deberny & Peignot(英语:Deberny & Peig
  • 前卫通用光盘FVD(Forward Versatile Disc),中文名称为前瞻多功能光碟或红光高清影音光碟(依据FVD官方网站之定名),是由台湾的工业技术研究院(工研院)所主导开发的光碟储存格式。其目标是在现有DV
  • 吉姆·亚伯特詹姆斯·安东尼·亚伯特(英语:James Anthony Abbott,1967年9月19日-),小名为吉姆·亚伯特(Jim Abbott),生于美国密歇根州弗林特,前美国职棒大联盟投手。他的右手掌先天残障,只能以左手