第六章连续小波变换
6.1概述
连续小波变换(CWT)通常被用于信号分析,即科学研究类。小波变换现在被大量不同的应用领域采纳,有时甚至会取代了傅里叶变换的位置,在许多领域都有这样的转变。例如很多物理学的领域亦经历了这个转变,包括分子动力学,从头计算(abinitiocalculations),天文物理学,密度矩阵局部化,地球物理学,光学,湍流,和量子力学。其他经历了这种变化的学科有图像处理,血压,心率和心电图分析,DNA分析,蛋白质分析,气象学,通用信号处理,语言识别,计算机图形学,和多分形分析。
小波分析的产生是由于,最初用于处理信息的技术,FT傅里叶变换,仅可在忽略时频分量的情况下进行,但是相对于实际应用情形,绝大多数信息的时频分量是一个不可忽视的因素,为此对FT进行一定程度的优化,得到STFT即短时傅里叶变换。短时距傅立叶变换是傅立叶变换的一种变形用于决定随时间变化的信号局部部分的正弦频率和相位。实际上,计算短时距傅立叶变换(STFT)的过程是将长时间信号分成数个较短的等长信号,然后再分别计算每个较短段的傅立叶转换。通常拿来描绘频域与时域上的变化,为时频分析中其中一个重要的工具。但是在实际应用的过程中我们又遇见了新的问题,人们无法知道信号的确切时频表示,即人们无法获知在何种时间实例中存在何种频谱分量。人们可知晓的是某些频段存在的时间间隔,(其根源可以追溯到海森堡不确定性原理,但在这里我们不做详细论述。)这是一个分辨率问题。
6.2例证
我们使用的窗口函数只是一个高斯函数,形式为:
e^{-a\left(\frac{t^2}{2}\right)}
其中a确定窗口的长度,t是时间。下图显示了由a的值确定的不同支持区域的四个窗口函数。请忽略a的数值,因为计算此函数的时间间隔也决定了函数。只需记下每个窗口的长度即可。上面给出的示例是使用第二个值a=0.计算得出的。现在,我们将显示上面给出的与其他窗口计算的相同信号的STFT。
首先,让我们看一下第一个最窄的窗口。我们预计STFT具有非常好的时间分辨率,但频率分辨率相对较差:
下图显示了此STFT。该图从顶部鸟瞰图显示,并带有一个角度,以便更好地解释。请注意,这四个峰值在时间上彼此之间有很好的分离。另请注意,在频域中,每个峰值都覆盖一个频率范围,而不是单个频率值。现在,让我们将窗口变宽,并查看第三个窗口(第二个窗口已在第一个示例中显示)。
请注意,与之前的情况不同,峰值在时间上彼此之间没有很好的分离,但是,在频域中,分辨率会进一步提升,我们进一步增加窗口的宽度。左图就明显的体现出了当窗口过宽时分辨率比较差。
鉴于上述问题CWT(连续小波变换)应运而生。
6.3STFT与CWT之间的两个主要区别
不采用窗口信号的傅里叶变换,因此将看到对应于正弦的单个峰值,即不计算负频率。
当为每个光谱分量计算变换时,窗口的宽度会发生变化,这可能是小波变换最重要的特征。
公式表达:
CWT_x^\psi(\tau,s)=\Psi_x^\psi(\tau,s)=\frac{1}{\sqrt{
s
}}\intx(t)\psi^*\left(\frac{t-\tau}{s}\right)dt
如上面的等式所示,变换后的信号是两个变量的函数,tau和s,分别是平移和标度参数。\psi(t)是变换函数,称为母小波,母小波是用于生成其他窗口函数的原型。
如图所示,信号由30Hz、20Hz、10Hz和5Hz的四个频率分量组成。
请注意,轴是平移和缩放,而不是时间和频率。然而,平移与时间严格相关,因为它指示母小波的位置。母小波的平移可以被认为是自t=0以来经过的时间。
较小的刻度对应于较高的频率,即频率随着刻度的增加而降低,因此,比例约为零的图形部分实际上对应于分析中的最高频率,而具有高尺度的标度对应于最低频率。
信号首先具有30Hz(最高频率)分量,这在0到30的转换时以最低刻度显示。然后是20Hz分量,第二高频率,依此类推。5Hz分量出现在平移轴的末端(如预期的那样),并按预期再次以较高的比例(较低频率)显示。
未完待续
前期回顾
第一章熵编码器/第二章BWT算法
第三章基于BWT的改进算法
第四章Lempel-ZivParsing
第五章双标准数据压缩
下期预告
第七章:OpenHarmony内核子系统之文件系统和压缩器参考文献
[1]DeorowiczS.UniversalLosslessDataCompressionAlgorithms[J].PhilosophyDissertationThesis,.