《小波变换》PPT课件.ppt

《《小波变换》PPT课件.ppt》由会员分享，可在线阅读，更多相关《《小波变换》PPT课件.ppt（122页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第章小波变换第章小波变换v8.1 连续小波变换的基本概念和性质连续小波变换的基本概念和性质v8.2 常用的小波函数常用的小波函数v8.3 尺度因子离散化的小波变换及小波标架尺度因子离散化的小波变换及小波标架v8.4 离散小波变换的多分辨率分析离散小波变换的多分辨率分析v8.5 Mallat算法及实现算法及实现v8.6 小波变换小结小波变换小结v 第章小波变换第章小波变换v 自从1822年傅里叶(Fourier)发表“热传导解析理论”以来，傅里叶变换一直是信号处理领域中最完美、应用最广泛、效果最好的一种分析手段，但傅里叶变换只是一种纯频域的分析方法，它在频域的定位性是完全准确的(即频域分辨率最高

2、)，而在时域无任何定位性(或分辨能力)，也即傅里叶变换所反映的是整个信号全部时间下的整体频域特征，而不能提供任何局部时间段上的频率信息。相反，当一个函数用函数展开时，它在时间域的定位性是完全准确的，而在频域却无任何定位性(或分辨能力)，也即函数分析所反映的只是信号在全部频率上的整体时域特征，而不能提供任何频率段所对应的时间信息。实际中，对于一些常见的非平稳信号，如音乐信号，在不同时间演奏不同音符；语音信号，在不同时间对应不同音节；地震信号，在目标出现的位置对应一个回波信号等，它们的频域特性都随时间而变化，因此也可称它们为时变信号。对这一类时变信号进行分析，v ，通常需要提取某一时间段(或瞬间)

3、的频域信息或某一频率段所对应的时间信息。因此，寻求一种介于傅里叶分析和分析之间的，并具有一定的时间和频率分辨率的基函数来分析时变信号，一直是信号处理界及数学界人士长期以来努力的目标。v 为了研究信号在局部时间范围的频域特征，1946年Gabor提出了著名的Gabor变换，之后又进一步发展为短时傅里叶变换(Short Time Fourier Transform，简记为STFT，又称为加窗傅里叶变换)。目前，STFT已在许多领域获得了广泛的应用，但由于STFT的定义决定了其窗函数的大小和形状均与时间和频率无关而保持固定不变，这对于分析时变信号来说是不利的。高频信号一般持续时间很短，而低频信号持续

4、时间较长，因此，我们期望对于高频信号采用小时间窗，对于低频信号则采用大时间窗进行分析。在进行信号分析时，这种变时间窗的要求同STFT的固定时窗(窗不随频率而变化)的特性是相矛盾的，这表明STFT在处理这一类问题时已无能为力了。此外，在进行数值计算时，人们希望将基函数离散v离散化，以节约计算时间及存储量，但Gabor基无论怎样离散，都不能构成一组正交基，因而给数值计算带来不便，这些是Gabor变换的不足之处，但恰恰是小波变换的特长所在。小波变换不仅继承和发展了STFT的局部化的思想，而且克服了窗口大小不随频率变化，缺乏离散正交基的缺点，是一种比较理想的进行信号处理的数学工具。v8.1 8.1 连

5、续小波变换的基本概念和性质连续小波变换的基本概念和性质v8.1.1 8.1.1 小波变换的定义小波变换的定义v 给定一个基本函数，令v v （8.18.1）v式中 均为常数，且 。显然，是基本函数先作移位再作伸缩以后得到的。若不断地变化，我们可得到一族函数 。给定二次方可积的信号 ，即 ，则 的小波变换（WT，Wavelet Transform）定义为v v （8.2）v式中 和 均是连续变量，是 的共轭函数，因此该式又称为连续小波变换（CWT）。如无特别说明，式中及以后各式中的积分都是从 到 。信号 的小波变换 是 和 的函数，是时移，是尺度因子。基本小波函数 又称为基本小波，或母小波。是母

6、小波经移位和伸缩所产生的一族函数，我们称之为小波基函数，或简称小波基。这样，式（8.2）中的WT又可解释为信号 和一族小波基的内积。v 基本小波可以是实函数，也可以是复函数。若 是实信号，则 是实函数，也是实函数，反之，为复函数。v在式（8.1）中，b的作用是确定对 分析的时间位置，也即时间中心。尺度因子 的作用是把基本小波 作伸缩。我们知道，由 变成 ，当 时，若 越大，则 的时域支撑范围（即时域宽度）较之 变得越大；反之，当 时，若 越小，则 的宽度越窄。这样，和 联合确定了对 分析的中心位置及分析的时间宽度，如图8.1所示。v 这样，式（8.2）的WT可理解为用一族分析宽度不断变化的基函

7、数对 作分析，由下一节的讨论可知，这一变化正好适应了我们对信号分析时在不同频率范围所需要不同的分辨率这一基本要求。式（8.1）中的因子 是为了保证在不同的尺度 时，始终能和基本小波 有着相同的能量，即v a)v b)v c)d)v 图8.1 基本小波的伸缩及参数和对分析范围的控制va)基本小波，b），c)不变，d)分析范围v令 ，则 ，这样，上式的积分即等于 。v令 的傅里叶变换为 ，的傅里叶变换为 ，由傅里叶变换的性质，的傅里叶变换为：v （8.3）v由Parsevals定理，式（8.2）可重新表示为：v （8.4）v此式即为小波变换的频域表达式。v8.1.2 8.1.2 小波变换的特点小波

8、变换的特点v 下面，我们从小波变换的恒Q性质、时域及频率分辨率以及和其它变换方法的对比来讨论小波变换的特点，以便对v小波变换有更深入的理解。v 比较式（8.2）和式（8.4），对小波变换的两个定义可以看出，如果 在时域是有限支撑的，那么它和 作内积后将保证 在时域也是有限支撑的，从而实现我们所希望的时域定位功能，也即使 反映的是 在b附近的性质。同样，若 具有带通性质，即 围绕着中心频率是有限支撑的，那么 和 作内积后也将反映 在中心频率处的局部性质，从而实现好的频率定位性质。显然，这些性能正是我们所希望的。问题是如何找到这样的母小波 ，使其在时域和频域都是有限支撑的。v 若 的时间中心是 ，

9、时宽是 ，的频率中心是v ，带宽是 ，那么 的时间中心仍是 ，但时宽变成 ，的频谱 的频率中心变为 ，带宽变成 。这样，的时宽带宽积仍是 ，与 无关。v这一方面说明小波变换的时频关系也受到不定原理的制约，但另一方面，也即更主要的是揭示了小波变换的一个性质，也即恒Q性质。定义:v =带宽/中心频率 (8.5)v为母小波 的品质因数，对 ，其带宽/中心频率=v因此，不论 为何值 ，始终保持了和 具有性同的品质因数。恒Q性质是小波变换的一个重要性质，也是区别于其它类型的变换且被广泛应用的一个重要原因。图8.2说明了和的带宽及中心频率随变化的情况。v(a)(b)(c)v 图8.2 随 变化的说明v将图

10、8.1和图8.2结合起来，我们可看到小波变换在对信号分析时有如下特点：当 变小时，对 的时域观察范围变窄，但对 在频率观察的范围变宽，且观察的中心频率向高频处移动，如图8.2.c所示。反之，当 变大时，对 v 的时域观察范围变宽，频域的观察范围变窄，且分析的中心频率向低频处移动，如图8.2b所示。将图8.1和8.2所反映的时频关系结合在一起，我们可得到在不同尺度下小波变换所分析的时宽、带宽、时间中心和频率中心的关系，如图8.3所示。v由于小波变换的恒Q性质，因此在不同尺度下，图8.3中三个时、频分析区间（即三个矩形）的面积保持不变。由此我们看到，小波变换为我们提供了一个在时、频平面上可调的分析

11、窗口。该分析窗口在高频端（图中 处）的频率分辨率不好（矩形窗的频率边变长），但时域的分辨率变好（矩形的时间边变短）；反之，在低频端（图中 处），频率分辨率变好，而时域分辨率变差。但在不同的值下，图8.3中分析窗的面积保持不变，也即时、频分辨率可以随分析任务的需要作出调整。v 图8.3 a取不同值时小波变换对信号分析的时频区间v众所周知，信号中的高频成份往往对应时域中的快变时则要求时域分辨率要好以适应快变成份成份，如陡峭的前沿、后沿、尖脉冲等。对这一类信号分析间隔短的需要，对频域的分辨率则可以放宽，当然，时、频分析窗也应处在高频端的位置。与此相反，低频信号往往是信号中的慢变成份，对这类信号分析时

12、一般希望频率的分辨率要好，而时间的分辨率可以放宽，同时分析的中心频率也应移到低频处。显然，小波变换的特点可以自动满足这些客观实际的需要。v总结上述小波变换的特点可知，当我们用较小的对信号作高频分析时，我们实际上是用低频小波对信号作概貌观察,小波变换的这一特点即既符合对信号作实际分析时的规律，也符合人们的视觉特点。综上所述，由于小波变换信号作概貌观察,小波变换的这一特点即既符合对信号作实际分析时的规律，也符合人们的视觉特点。综上所述，由于小波变换具有恒Q性质及自动调节对信号分析的时宽/带宽等一系列突出优点，因此被人们称为信号分析的“数学显微镜”。v8.1.3 8.1.3 连续小波变换的基本性质连

13、续小波变换的基本性质v1时移性质v若 的CWT是 ，那么 的CWT是 。该结论极易证明。记 ，则v （8.6）2尺度转换性质 如果 的CWT是 ，令 ，则 （8.7）v证明：，令 ，v则 v该性质指出，当信号的时间轴按 作伸缩时，其小波变换在 和 两个轴上同时要作相同比例的伸缩，但小波变换的波形不变。这是小波变换优点的又一体现。v3微分性质 如果 的CWT是 ，令 ，则 （8.8）v证明：v由式（8.6）的移位性质，有v即 v4两个信号卷积的CWT 令 的CWT是 及 ，并令 ，则有 v （8.9）v式中符号 表示对变量 作卷积。v证明：v由式（8.6）的移位性质，有v同理，v于是式（8.9）

14、得证。v5两个信号和的CWTv 令 的CWT分别是 ，且v ，则v （8.10）v同理，如果 ，则v （8.11）v式（8.10），（8.11）说明两个信号和的CWT等于各自CWT的和，也即小波变换满足叠加原理。v6小波变换的内积定理v 设 和 ，的小波变换分别是v 和 ，则v （8.12）v式中 ，为 的傅里叶变换。v证明：由式（8.4）关于小波变换的频域定义，式（8.12）的左边有：v假定积分v存在，再由Parseval定理，上述的推导最后为v于是定理得证。v式（8.12）实际上可看作是小波变换的Parseval定理。该式又可写成更简单的形式,即v (8.13)v进一步，如果令 ，有v (

15、8.14)v该式更清楚地说明，小波变换的幅平方在尺度位移平面上的加权积分等于信号在时域的总能量，因此，小波变换的幅平方可看作是信号能量时频分布的一种表示形式。v式（8.12）和式（8.13）中对 的积分是从 ，这是因为我们假定 总为正值。这两个式子中出现的 是由于定义小波变换时在分母中出现了 ，而式中又要对 作积分所引入的。v读者都熟知傅里叶变换中的Parseval定理，即时域中的能量等于频域中的能量。但小波变换的Parseval定理稍为复杂，它不但要有常数加权，而且以 的存在为条件。v8.1.48.1.4小波反变换及小波容许条件小波反变换及小波容许条件v 下面给出连续小波反变换的公式及反变换

16、存在的条件。v设 ，记 为 的傅里叶变换，若v则 可由其小波变换 来恢复，即v (8.15)v证明：设 ，,则v将它们分别代入式（8.12）的两边，再令 ，于是有v v于是定理得证。v在式(8.12)、式(8.15)中，结论的成立都是以 为前提条件的，这又称为“容许条件（admissibility condition）。该容许条件含有多层的意思：v 1.并不是时域的任一函数 都可以充当小波。其可以作为小波的必要条件是其傅里叶变换满足该容许条件；v2.若 ，则必有 ，否则 必趋于无穷。这等效地告诉我们，小波函数 必然是带通函数；v3.由于 ，因此必有v (8.16)v这一结论指出，的取值必然是有

17、正有负，也即它是振荡的。v以上三条给我们勾画出了作为小波的函数所应具有的大致特征，即 是一带通函数，它的时域波形应是振荡的。此外，从时频定位的角度，我们总希望 是有限支撑的，因此它应是快速衰减的。这样，时域有限长且是振荡的这一类函数即是被称作小波（wavelet）的原因。由上述讨论，自然应和一般的窗函数一样满足：v （8.17）v并且由后面的讨论可知，尺度因子 常按 来离散化，v。由式（8.3），对应的傅里叶变换 ，由于需要在不同的尺度下对信号进行分析，同时也需要在该尺度下由 来重建 ，因此要求 是有界的，当 由 时，应有v (8.18)v式中，。该式称为小波变换的稳定性条件，它是在频域对小波

18、函数提出的又一要求。满足式（8.18）的小波称作“二进（dyadic）”小波。v8.1.5小波变换的充要条件v我们在上一节指出，并不是时域任一函数都可以用作小波 。可以作为小波的函数至少要满足容许条件。与此结论相类似，并不是 平面上的任一二维函数 都v对应某一函数的小波变换。如果是某一时域信号，如v 的小波变换，它应满足一定的条件。v设 是 平面上的任一点，上的二维函数 欲是某一函数的小波变换的充要条件是它必须满足如下的重建核方程，即v （8.19）v式中 是 在 处的值，v （8.20）v称为重建核。v证明：由式（8.2）小波变换的定义，有v将式（8.15）代入该式，有v式（8.19）的重建

19、核方程和式（8.20）的重建核公式说明，若 是 的小波变换，那么在 平面上某一点 v 处小波变换的值 可由半平面上 的值 来表示，也即，是半平面上 的总贡献。v既然 平面上各点的 可由式（8.19）互相表示，因此这些点上的值是相关的，也即式（8.15）对 的重建是存在信息冗余的。这一结论告诉我们可以用 平面上离散栅格上的 来重建 ，以消除重建过程中的信息冗余。v我们知道，当用 的短时傅里叶变换 来重建v 时，平面上的信息也是有冗余的，即 平面上各点的 是相关的，因此引出了离散栅格上的STFT，进一步的发展即是信号的Gabor展开与Gabor变换。由此可以得出，将一个一维的函数映射为一个二维函数

20、后，在二维平面上往往会存在信息的冗余，由此引出了二维函数的离散化问题及标架理论。有关离散小波变换及小波标架的内容将在后面予以讨论。v重建核 是小波 和 处的小波 的内积，因此 反映了 和 的相关性。v若 ，即两个小波重合时，取最大值；若 远离 ，则 将迅速减小。若能保证 v，则 平面上各点小波变换的值将互不相关。这等效地要求对任意的尺度 及位移 ，由母小波 形成的一族v是两两正交的。可以想象，若 连续取值，要想找到这样的母小波 使 两两正交，那将是非常困难地。因此，连续小波变换 必然存在信息冗余。然而，当 离散取值时，则有可能得到一族正交小波基 。v8.2 8.2 常用的小波函数常用的小波函数

21、v由前面叙述可知，作为一个小波的函数 ，它一定要满足容许条件，在时域一定要是有限支撑的，同时，也希望在频域也是有限支撑的，当然，若时域越窄，其频域必然是越宽，反之亦然。在时域和频域的有限支撑方面我们往往只能取一个折中。此外，希望由母小波 形成的 是两两正v交的或是双正交的,进一步，希望 有高阶的消失矩，希望与 相关的滤波器具有线性相位等,可以根据上述要求对现已提出的大量的小波函数作一粗略地分类。在下面的分类中，第一类是所谓地“经典类小波”，在MATLAB中把它们称作“原始（Crude）小波”。这是一批在小波发展历史上比较有名的小波；第二类是Daubecheis构造的正交小波，第三类是由Cohe

22、n，Daubechies构造的双正交小波。v8.2.18.2.1经典类小波经典类小波v8.2.1.1 Haar小波v Haar小波来自于数学家Haar于1910年提出的Haar正交函数集，其定义是：v （8.21）v其波形如图8.4（a）所示。v 的傅里叶变换是：v （8.22）vHaar小波有很多优点，如：vHaar小波在时域是紧支撑的，即其非零区间为（0，1）；v若取 ，那么Haar小波不但在其整数位移处是正交的，即 ，而且在 取不同值时也是两两正交的，即 ，如图8.3(b)和(c)所示。所以Haar小波属正交小波；vHaar波是对称的。我们知道，系统的单位抽样响应若具有对称性，则该系统具

23、有线性相位，这对于去除相位失真是非常有利的。Haar小波是目前唯一一个既具有对称性又是有限支撑的正交小波；vHaar小波仅取1和1，因此计算简单。但Haar小波是不连续小波，由于 ，因此 在 处只有一阶零v点，这就使得Haar小波在实际的信号分析与处理中受到了限制。但由于Haar小波有上述的多个优点，因此在教科书与论文中常被用作范例来讨论。v 图8.4 Harr小波v8.2.1.2 Morlet小波vMorlet小波定义为v (8.23)v其傅里叶变换v （8.24）v它是一个具有高斯包络的单频率复正弦函数。考虑到待分析的信号一般是实信号，所以在MATLAB中将式（8.23）改造为：v （8.

24、25）v并取 。该小波不是紧支撑的，理论上讲 可取v 。但是当 ，或再取更大的值时，和 在时域和频域都具有很好的集中，如图8.5所示。vMorlet小波不是正交的，也不是双正交的，可用于连续小波变换。但该小波是对称的，是应用较为广泛的一种小波。v (a)时域波形 (b)频谱v 图8.5 Morlet小波v8.2.1.3 Mexican hat小波v该小波的中文名字为“墨西哥草帽”小波，又称Marr小波。它定义为 (8.26)v式中 ，其傅里叶变换为v (8.27)v 该小波是由一高斯函数的二阶导数所得到的，它沿着中心轴旋转一周所得到的三维图形犹如一顶草帽，故由此而得名。其波形和频谱如图8.6所

25、示。v 该小波不是紧支撑的，不是正交的，也不是双正交的，但它是对称的，可用于连续小波变换。由于该小波在 处有二阶零点，因此它满足容许条件，且该小波比较接近人眼视觉的空间响应特征，因此它在1983年即被用于计算机视觉中的图像边缘检测。v 图8.6 墨西哥草帽小波，(a)时域波形，(b)频谱v8.2.1.4 Gaussian小波v高斯小波是由一基本高斯函数分别求导而得到的，定义为：v ，(8.28)v式中定标常数是保证 。v该小波不是正交的，也不是双正交的，也不是紧支撑的。当 取偶数时 v 正对称，当 取奇数时，反对称。图8.7给出了 时 的时域波形及对应的频谱。v 图8.7 高斯小波（取 ）(a

26、)时域波形，(b)频谱v8.2.28.2.2正交小波正交小波v 目前提出的正交小波大致可分为四种，即Daubechies小波，对称小波，Coiflets小波和Meyer小波。这些正交小波和前面所讨论的“经典小波”不同，它们一般不能由一个简洁的表达式给出，而是通过一个叫做“尺度函数（Scalling function）”的 的加权组合来产生的。尺度函数是小波变换的又一个重要概念。由下节讨论可知，小波函数 ，尺度函数 同时和一个低通滤波器 及高通滤波器 相关连，和 可构成一个两通道的分析滤波器组。这些内容构成了小波变换的多分辨率分析的理论基础。因此，在讨论正交小波时，同时涉及到尺度函数 ，分析滤波

27、器组 ，及综合滤波器组 ，。MATLAB中的Wavelet Toolbox中有相关的软件来产生各类正交小波及其相应的滤波器。v8.2.2.1 Daubechies小波v Daubechies小波简称db小波。它是由法国学者Ingrid Dauechies于90年代初提出并构造的。Daubechies对小波变v换的理论做出了突出的贡献，特别是在尺度取2的整数次幂时的小波理论及正交小波的构造方面进行了深入的研究，其代表作Ten Lectures on Wavelet深受同行们的欢迎。vdbN中的表示db小波的阶次，。当时，db1即是Haar小波。因此，前述的Haar小波应归于“正交小波”类。Dau

28、bechies计算出了 时的 及 。在MATLAB5.3中，N的阶次还可以扩展。db小波是正交小波，当然也是双正交小波，并是紧支撑的。的支撑范围在v ，的支撑范围在 。小波 具有N阶消失矩，在 处具有阶零点。但db小波是非对称的，其相应的滤波器组属共轭正交镜像滤波器组（CQMFB）。图8.8给出了 时，及 ，的波形。图8.8 时db小波(a)，(b)，(c)，(d)v8.2.2.2 对称小波v对称小波简记为symN，它是db小波的改进，也是由Daubechies提出并构造的。它除了有db小波的特点外，主要是 是接近对称的，因此，所用的滤波器可接近于线性相位。图8.9是 时的对称小波。v (a)

29、(b)v 图8.9 时的对称小波(a)，(b)v8.2.2.3 Coiflets小波v该小波简记为coifN，.在db小波中，Daubechies小波仅考虑了使小波函数 具有消失矩（N阶），而没考虑尺度函数 。R.Coifman于1989年向Daubechies提出建议，希望能构造出使 也具有高阶消失矩的正交紧支撑小波。Daubechies接受了这一建议，构造出了这一类小波，并以Coifman的名字命名。vcoifN是紧支撑正交、双正交小波，支撑范围为 ，也是接近对称的。的消失矩是2N，的消失矩是2N-1。图8.10是N=4时的coif4小波。v 图8.10 N=4时v 的Coiflets小v

30、 波，a)，v (b)v8.2.2.4 Meyer小波vMeyer小波简记为meyr，它是由Meyer于1986年提出的。该小波无时域表达式，它是由一对共轭正交镜像滤波器组的频谱来定义的。vMeyer小波是正交、双正交的，但不是有限支撑的，但其有效的支撑范围在 到 之间。该小波是对称的，且有着非常好的规则性。图8.11给出了Meyer小波的尺度函数 和小波函数 。v 图8.11 v Meyer小波v a)，v (b)v8.2.3 8.2.3 双正交小波双正交小波v我们知道，两通道正交镜像滤波器组的分析滤波器 和v 是功率对称的，且 和 之间有着正交性，再者 ，和 有着同样的长度，都不是线性相位

31、的。为了取得线性相位的滤波器组，需放弃 的功率互补性质。这也就放弃了 和 之间的正交性，代之的是双正交关系。v 由于离散小波变换最后是由两通道滤波器组来实现。因此，正交小波条件下的 ，和 与 都不具有线性相位（Haar小波除外）。为此，Daubechies和Cohen提出并构造了双正交小波，其目的是在放宽小波正交性的条件下得到线性相位的小波及相应的滤波器组。v双正交滤波器组简称biorNr,Nd，其中Nr是低通重建滤波器的阶次，Nd是低通分解滤波器的阶次。在MATLAB中，Nr和Nd的可能组合是：v Nr=1,Nd=1,3,5v Nr=2,Nd=2,4,6,8v Nr=3,Nd=1,3,5,7

32、,9v Nr=4,Nd=4v Nr=5,Nd=5v Nr=6,Nd=8v这一类小波自然不是正交的，但它们是双正交的，是紧支撑的，更主要的是它们是对称的，因此具有线性相位。分解小波 的消失矩为Nr-1。图8.12给出的bior3.7的分解小波、尺度函数及重建小波和尺度函数。v 图8.12 双正交小波bior3.7 v (a)分解尺度函数 (b)分解小波 v (c)重建尺度函数 (d)重建小波v8.2.48.2.4连续小波变换的计算连续小波变换的计算v在式（8.2）关于小波变换的定义中，变量 ，和 都是连续的，当我们在计算机上实现一个信号的小波变换时，和 均应离散化。对 离散化最常用的方法是取 ，

33、如取 ，这样 。对于 按2的整次幂取值所得到的小波习惯上称之为“二进（dyadic）”小波。对这一类小波的小波变换，我们可用有关离散小波变换的方法来实现。然而取 ，在实际工作中有时显得尺度跳跃太大。当希望 任意取值 ，也即在 的范围内任意取值时，这时的小波变换即是连续小波变换。v计算式（8.2）的最简单的方法是用数值积分的方法，即令v (8.29)v由于在 的区间内，所以上式又可写为：v (8.30)v由该式可以看出，小波变换 可看作是 和v的卷积后的累加所得到的结果，卷积的中间变量是 ，卷积后的变量为 及 。MATLAB中的cwt.m即是按此思路来实现的。具体过程大致如下：v1.先由指定的小

34、波名称得到母小波 及其时间轴上的刻度，假定刻度长为 ；v2.从时间轴坐标的起点开始求积分 ，v3.由尺度因子 确定对上述积分值选择的步长，越大，上述积分值被选中的越多；v4求 和所选中的积分值序列的卷积，然后再作差分，即完成式（8.30）。v本方法的不足之处是在 变化时，式（8.30）中括号内的积分、差分后的点数不同，也即和 卷积后的点数不同。解决的方法是在不同的尺度下对 作插值，使其在不同的尺度下，在其有效支撑范围内的点数始终相同。有关CWT快速计算的方法还可借助于CZT及梅林变换等方法，此处不再讨论。v例8.1令 为一正弦加噪声信号，它取自MATLAB中的noissin.mat。对该信号作

35、CWT，分别等于2和128，时，小波变换的结果对应信号中的高频成份，时，小波变换对应信号中的低频成份。其原始信号及变换结果见图8.13(a)，(b)和（c）。v例8.2 仍然使用例8.1的信号“noissin”，对其作CWT时 分别取10，30，60，90，120及150。所得到的图8.14是在各个尺度下的小波系数的灰度图。颜色越深，说明在该尺度及该位移（水平轴）处的小波系数越大。此例旨在说明对小波变换的结果具有不同的表示方式。v图8.13 信号“noissin”的小波变换，(a)原信号，v (b)，(c)v 图8.14 多尺度下小波变换的灰度表示v8.3 8.3 尺度离散化的小波变换及小波标

36、架尺度离散化的小波变换及小波标架v我们在式（8.2）定义了信号 的连续小波变换，式中，v 和 都是连续变量。为了在计算机上有效地实现小波变换，自然应取离散值，和 也应取离散值。从减少信息冗余的角度，和 也没有必要连续取值。v 和 形成了一个二维的“尺度位移”平面。前已述及，v 越大，对应的频率越低，反之，对应的频率越高。因此，平面也可视为“时频平面”。对同一个信号 ，我们已给出过不同的表示形式，如STFT，Gabor变换，WVD及本章的小波变换。v现重写几个有关的公式，即v (8.31)v (8.32)v (8.33)v (8.34)v其中式（8.32）是用时频平面离散栅格 上的点来表示，即G

37、abor展开，式（8.33）是具有双线性变换的表示形式，它和其它三种表示形式有较大的区别。式（8.31）和式（8.34）说明同一信号在时频平面上具有不同的表示形式。式（8.31）的反变换是有信息冗余的，即不需要 的所有的值就可恢复。同理，式（8.34）的小波变换也存在着信息冗余。在这两个式子中，我们只需取时频平面上的离散栅格处的点即可。问题的关键是如何决定 和 抽样的步长以保证对 的准确重建。下面，我们首先考虑尺度因子 的离散化，然后再考虑 和 的同时离散化。v8.3.1 8.3.1 尺度离散化的小波变换尺度离散化的小波变换v目前通用的对 离散化的方法是按幂级数的形式逐步加大 ，即令 。若取

38、，则v (8.35)v称为“半离散化二进小波”，而v (8.36)v称为二进小波变换。v设母小波 的中心频率为 ，带宽为 ，当 时，v 的中心频率变为 ，带宽 v。若 时，的中心频率和带宽分别是：v 和 。从对信号作频域分析的v角度，我们希望当 由 变成 时，和 在频域对应的分析窗为 v和 能够相连接。这样，当 由0变至无穷时，的傅里叶变换可以覆盖整个 轴。显然，若令母小波 的 ，则上面两个频域窗首尾相连，即 和 首尾相连。通过对母小波作合适的调制，可以方便地做到 。v现在，我们来讨论如何由式（8.36）的 来恢复 ，设 是 的对偶小波，并令 和 取类似的形式，即 (8.37)v这样，通过对偶

39、小波，我们希望能重建 ：v (8.38)v为了寻找 和 应满足的关系，现对上式作如下改变：v式中F代表求傅里叶变换。由式（8.3）和式（8.4），有v (8.39)v显然，若v (8.40)v则式（8.39）的右边变成的傅里叶反变换，自然就是 。v对于满足容许条件的小波 ，当 时，其二进制小波 对应的傅里叶变换应满足式（8.18）的稳定性条件。这样，结合式（8.18）和式（8.40），我们可由下式得到对偶小波 :v (8.41)v由于式（8.41）的分母满足式（8.18），因此有v (8.42)v这样，对偶小波 也满足稳定性条件，也即，总可以找到一个“稳定的”对偶小波 由式（8.38）重建出

40、。下面定理更完整地回答了在半离散二进小波变换情况下的 重建问题。v定理8.1 如果存在常数 ，使得v (8.43)v则v (8.44)v如果 满足v (8.45)v则v (8.46)v该定理指出，若 的傅里叶变换满足稳定性条件，则 在v 上的小波变换的幅平方的和是有界的。进而，和 的傅里叶变换若满足式（8.45）（也即式（8.40），则可由式（8.46）重建。v总之，若满足容许条件，且再满足稳定性条件，由二进小波变换总可以重建，也即一个满足稳定性条件的对偶小波总是存在的。但是，满足稳定性条件的对偶小波不一定是唯一的。如何构造“好”的小波及得到唯一的对偶小波是小波理论中的重要内容。v若式（8.4

41、3）的稳定性条件满足，则 的容许条件必定满足，且v (8.47)v从而，由连续小波变换 总可以恢复 ，也即式（8.15）总是成立。v以上讨论的是仅对 作二进制离散化的情况，现在考虑 和 同时离散化的相应理论问题。v8.3.28.3.2离散栅格上的小波变换离散栅格上的小波变换v令 ，我们可实现对 的离散化。若 ，则 。欲对 离散化，最简单的方法是将 均匀抽样，如令 ，的选择应保证能由 来恢复出 。当 时，将 由 变成 时，即是将 扩大了 倍，这时小波 的中心频率比 的中心频率下降了 倍，带宽也下降了 倍。因此，这时对v 抽样的间隔也可相应地扩大 倍。由此可以看出，当尺度 分别取 ，对 的抽样间隔

42、可以取v ，这样，对 和 离散化后的结果是：v (8.48)v对给定的信号 ，式（8.2）的连续小波变换可变成如下离散栅格上的小波变换，即v (8.49)v此式称为“离散小波变换（Discrete Wavelet Transform，DWT）”，注意式中 仍是连续变量。v记 ，我们可以仿照傅里叶级数和Gabor展开那样来重建 ，即v (8.50)v该式称为小波级数，称为小波系数，是 的对偶函数，或对偶小波。v我们知道，对任一周期信号 ，若周期为T，且v ，则 可展成傅里叶级数，即v (8.51)v式中 是 的傅里叶系数，它由下式求出：v (8.52)v小波级数和傅里叶级数形式上类似，但其物理概

43、念却有着明显的不同：v(1)傅里叶级数的基函数 ，是一组正交基，即v 。而小波级数所用的一族函数 不一定是正交基，甚至不一定是一组“基”；v(2)对傅里叶级数来说，基函数是固定的，且分析和重建的基函数是一样的，即都是 （差一负号）；对小波级数来说，分析所用的函数 是可变的，且分析和重建所用的函数是不相同的，即分析时是 ，而重建时是 ；v(3)在傅里叶级数中，时域和频域的分辨率是固定不变的，而小波级数在 轴上的离散化是不等距的，这正体现了小波变换“变焦”和“恒Q”性的特点。v将式（8.2）的连续小波变换改变成式（8.49）的离散小波变换，人们自然会问：v(1)一族小波函数 ，在空间 上是否是完备

44、的？所谓完备，是指对任一 ，它都可以由这一组函数（即 ）来表示；v(2)如果 是完备的，那么 对 的表示是否有信息的冗余？v(3)如果 是完备的，那么对 和 的抽样间隔如何选取才能保证对 的表示不存在信息的冗余？vDaubechies对上述问题进行了深入的研究，给出了“小波标架”的理论，现介绍一下其中主要的结论。v8.3.38.3.3小波标架理论小波标架理论v标架的基本理论要点是：v(1)若 是Hilbert空间中的一组向量，对给定的v ，若存在常数 ，满足v (8.53)v则 构成了一个标架；v(2)若 则称 为紧标架，若 ，则 构成一正交基；（3）定义标架算子为v (8.54)v则 (8.

45、55)v记 为 的对偶函数族，则 也构成一个标架，标架界分别 为和 ；v（4）用标架来表征一个信号x，也即对x作分解时，标架可给出完备的且是稳定的表示，但这种表示是冗余的，即v 之间是线性相关的，因此 不是唯一的。对信号的冗余表示有时并不一定是坏事，它在表示的稳定性、对噪声的鲁棒性（robustness）方面都优于正交基；v（5）标界边界B和A之比值，即B/A称为冗余比。在实际工作中，总希望B/A接近于1，即 为紧标架。当A=B 时，我们有v (8.56)v将以上要点内容用于小波变换，即得小波标架。在式（8.48）中，令 ，我们从而得到了一族在尺度和位移上均是离散的小波 。能否由离散小波变换

46、来重建 ，显然取决于 。v 和 越小，重建越容易，当然冗余度也越大，对不同的 是线性相关的，这时将有无数的 存在。当然 ，过大，准确重建将不会可能。v下面两个定理给出了小波标架的主要内容。v定理8.2 如果 构成 中的一个标架，且标架边界分别为A和B，则母小波须满足：v (8.57)v该定理又称 构成标架的必要条件。这一条件实际上即是连续小波变换中的容许条件。当仅对 取二进制离散化，v 保持连续时，该必要条件也就是充分条件。v若 构成紧标架，即A=B，那么其标架边界v (8.58)v若 构成 中正交基，则v v (8.59)v定理8.3 定义v (8.60)v及 v (8.61)v如果 和 的

47、选取保证v (8.62)v及v (8.63)v则 是 中的一个标架。、分别是标架界A和B的下界与上界。v总之，以上的标架理论及边界A、B值的计算给我们一个大致估计 选取的原则，即两者的选取要保持离散化后的v 至少要构成一个标架，以保证对信号稳定、完备的表示。但在一般情况下，标架并不是正交基，除非A=B=1。v若 是由母小波 通过伸缩与移位生成的v 上的“稠密”的二维函数族，并且存在常数A和B，使得v (8.64)v对于所有满足平方和的序列 成立，式中v则称 是 上的一个Riesz基，常数 vA和B分别称为Riesz基的下界和上界。v上述定义中“稠密”的含义是指 中的任一函数都可由二维序列 的线

48、性组合来表示。其实，该定义可简单地解释为如下：v（1）首先，是一个标架；v（2）对任意的 ，之间是线性无关的。v这样，Riesz基可以比标架最大限度地去除冗余度。此外，生成Riesz基 的母小波 称为Riesz函数。v可以证明，Riesz基的对偶函数序列 也是一个Riesz基，因此 对任意的 是线性无关的，对给定的 ，其对称基 是唯一的。这样，有v （8-56）v下面，我们在前面关于小波分类的基础上再给出几个有关小波的定义：v1正交小波v若Riesz基 满足v (8.66)v则称生成 的母小波 为正交小波。式中v (8.67)v式（8.66）指出，在同一尺度j下，不同移位之间的 是正交的。同时

49、，在同一位移k下，不同尺度j之间的 也是正交的。v2半正交小波v若 满足v 对 （8.68）v该式的含义是，若 ，则 。这时，对不同的位移k之间 不是正交的。因此，生成 的 称为半正交小波。v3.双正交小波v若 和其对偶小波 之间满足v (8.69)v则称生成 的 为双正交小波。v半正交小波不是正交小波，双正交小波指的是 和其对偶v 之间的关系，因此也不是正交小波。但一个正交小波必定是半正交的，也是双正交的。下面给出正交小波、半正交小波及双正交小波之间的关系。v令 是一个半正交小波，其傅里叶变换为 ，定义v v (8.70)v并记 的傅里叶反变换为 ，则由 和 分别作二进制伸缩和移位生成的 和

50、 之间是双正交的，即它们满足式（8.69）。v若 为正交小波，则式（8.70）的分母为1，这样，v ，也即 。这正是我们以前所指出的，即正交基和其对偶基是一样的。因此，令v (8.71)v并记 为 的傅里叶反变换。由上述可知，的对偶小波 应由下式给出：v （8.72)v可以证明，即 和其对偶函数 是自对偶的，因此，即是正交小波。v8.4 8.4 离散小波变换的多分辨率分析离散小波变换的多分辨率分析v前面我们介绍了连续小波变换的定义与性质，给出了在 平面上离散栅格上小波变换的定义及与其有关的标架问题。在这两种情况下，时间t仍是连续的。在实际应用中，特别是在计算机上实现小波变换时，信号总要取成离散