码迷,mamicode.com
首页 > 其他好文 > 详细

小波变换教程(五)

时间:2019-11-15 12:25:37      阅读:63      评论:0      收藏:0      [点我收藏+]

标签:等等   频率分辨率   image   一个   解决   net   pre   tco   正是   

译文转:https://blog.csdn.net/alihouzi/article/details/45190303

原文转:http://users.rowan.edu/~polikar/WTpart1.html

五、终极解决方案:小波变换

        小波变换是这样一种变换:它提供了信号的时频表示(还有一些变换可以提供这些信息,如短时傅立叶变化,魏格纳分布等等)。

        在任一刻出现的特殊的频谱分量都有特殊的意义。这种情况下,如果知道了这些特殊的频谱分量出现的时间会比较有益。举例来说,在脑电图中,一个事件相关电位的潜伏期是主要关注点(事件相关电位是指大脑对某一特定刺激的反应,例如闪灯,这种反应的潜伏期是从刺激的开始到反应发生这段时间)。

        小波变换可以同时提供时间和频率信息,因此给出了信号的一种时频表示。

        但是小波变换到底是如何变换的仍然是一个不同的有趣故事,需要在理解了短时傅立叶变换(STFT)之后再解释。小波变换是被用来替代短时傅立叶变换(STFT)的。我们将在后面详细阐述STFT。现在可以说,一些在STFT中遇到的有关分辨率的问题,可以用小波变换解决。

        为了长话短说,我们略过一些时域信号的高通和低通滤波处理,这些滤波器用来滤除信号中的低频和高频部分分量。这个过程是重复进行的,每一刻都可以从信号中滤除一些频率分量。

        这里解释一下滤波过程是如何工作的:假定我们有一个信号,其中最高频分量为1000Hz。第一步,我们通过高通和低通滤波器把信号分成两部分(滤波器必须满足某些特定的条件,即允许条件),结果产生了同一信号的两个不同版本,0-500Hz的信号(低通)和500-1000Hz的信号(高通)。

        然后,我们可以拿任意一部分(通常是低频部分)或者二者来做相同的处理。这个过程叫做分解。

        假设我们拿低频部分做了处理,现在我们就有了3列数据,分别为0-250Hz,250-500Hz和500-1000Hz。

        然后再对低通滤波过的信号做高通和低通滤波处理,现在我们就有了4列数据,分别为0-125Hz,125-250Hz,250-500Hz和500-1000Hz。我们持续进行这个处理过程,直到将信号分解到一个预定义的级别。这样我们就有了一系列信号,这些信号实际上表示相同的信号,但是每一个序列都有不同的频带。我们知道哪些信号对应哪个频段,如果我们将这些信号放在一起画出三维图,一个轴表示时间,频率在另外一个轴上,幅度在第三个轴上。这幅图会告诉我们在某时刻出现的是什么频率的信号(这里有一个问题,叫做“不确定性原理”,即我们不能确切的知道哪个频率出现在哪个时间点上,我们仅仅知道那个频段出现在哪个时间段内,后文中将有更多关于此的介绍)。

        不过,我仍然想用一个简明扼要的方式解释它:

        最初是由海森堡发现并阐述了不确定性原理(测不准原理),这个原理是这样的,移动粒子的动量和位置不能同时确定,对我们这个课题则是这样的:

        在时频平面内的一个确切的点上,信号的频率和发生时间不能同时确定。换句话说:在任意一个时间点,我们不能确定哪个频谱分量存在。我们能做到的是在一个给定的时间段内确定哪个频谱分量存在。这是一个设计到分辨率的问题,正是因为这个原因,才使得研究者们从快速傅立叶变换(STFT)转到了小波变换(WT)上。快速傅立叶变换的分辨率是固定的,而小波变换则能给出不同的分辨率:

        高频信号在时域内可以得到很好的解决,低频信号则可以在频域内得到很好的解决,这意味着,相对于低频分量,高频分量更容易在时域内定位(因为有更小的相对误差)。反而言之,低频分量更容易在频域内定位。看下面这张图:

  f ^
          |*******************************************         continuous
          |*  *  *  *  *  *  *  *  *  *  *  *  *  *  *         wavelet transform
          |*     *     *     *     *     *     *
          |*           *           *           *
          |*                       *
          --------------------------------------------> time

        对上图的解释是:最上面一行是对高频信号的多个采样,其采样间隔时间也较短。就是说高频信号更容易在时域内处理。最下面的一行是对低频信号的采样,特征点较少,因此,低频信号在时域内并不容易处理。

       ^frequency
        |
        |
        |
        | *******************************************************
        |
        |
        |
        | *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *   discrete time
        |                                                           wavelet transform
        | *     *     *     *     *     *     *     *     *     *
        |
        | *           *           *           *           *
        | *                       *                       *
        |----------------------------------------------------------> time
      

        在离散时间域内,信号的时间分辨率如上图所示,但是现在,频率信息的分辨率在每一个阶段都不相同。注意,低频信号更容易在频域内处理,高频则正好相反。注意到图中后续频率分量的间隔是如何随频率增高(这里应为降低)而增大的。

        下图是一些连续小波变换的例子。让我们以一个正弦波为例,这个正弦波包含了出现在不同时刻的两个频率分量。前半部分是低频信号,后面是高频信号。

                        技术图片技术图片

                                                                                        图1.10

对上图的连续小波变换如下:

                        技术图片技术图片

 

 

                                                                         图1.11

        注意到上图中,以“尺度”为标签的轴代表频率。“尺度”这个概念将会在后续章节进行阐述,但是这时需要注意的是尺度正好与频率成反比,即:尺度越大频率越低,尺度越小频率越高。因此,图中的小尖峰反映了信号中的高频分量,大的尖峰则反映了信号中的低频分量(在时域内低频信号是先出现的)。

        你可能被图中的频率分辨率搞迷惑了,因为高频信号也得到了很好的频率分辨率。但是需要注意的是,高频信号的尺度分辨率很高,高尺度分辨率意味着低频率分辨率,反之则反。更多关于此的介绍在文章的后续部分。

小波变换教程(五)

标签:等等   频率分辨率   image   一个   解决   net   pre   tco   正是   

原文地址:https://www.cnblogs.com/sggggr/p/11865196.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!