小波变换教程(五)

小波变换网文精粹：小波变换教程(五)

原文：ROBI POLIKAR. THE ENGINEER'S ULTIMATE GUIDE TO WAVELET ANALYSIS：The Wavelet Tutorial

网址：http://users.rowan.edu/~polikar/WAVELETS/WTtutorial.html

译文转自：http://blog.163.com/renfengyuee@126/blog/static/359431362009109102040700/

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

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

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

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

        但是小波变换到底是如何变换的仍然是一个不同的有趣故事，需要在理解了短时傅立叶变换(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.7

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

                        

                                                                         图1.8

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

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