小波包分解算法-Readppt课件.ppt
资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值第10章 小波包变换及其应用 简介简介小波包的定义与性质小波包的定义与性质 小波空间的精细分割小波空间的精细分割 小波包滤波器组小波包滤波器组 最佳小波包基的选取最佳小波包基的选取 小波包变换的应用小波包变换的应用 小波分析及其工程应用小波分析及其工程应用-清华大学计算机系清华大学计算机系-孙延奎孙延奎-2005春春资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值简介由于由于正交小波变换正交小波变换只对信号的低频部分做进一步只对信号的低频部分做进一步分解，而对高频部分也即信号的细节部分不再继分解，而对高频部分也即信号的细节部分不再继续分解，所以小波变换能够很好地表征一大类以续分解，所以小波变换能够很好地表征一大类以低频信息为主要成分的信号，但它不能很好地分低频信息为主要成分的信号，但它不能很好地分解和表示包含大量细节信息（细小边缘或纹理）解和表示包含大量细节信息（细小边缘或纹理）的信号，如非平稳机械振动信号、遥感图象、地的信号，如非平稳机械振动信号、遥感图象、地震信号和生物医学信号等。与之不同的是，震信号和生物医学信号等。与之不同的是，小波小波包变换可以对高频部分提供更精细的分解包变换可以对高频部分提供更精细的分解，而且，而且这种分解既无冗余，也无疏漏，所以对包含大量这种分解既无冗余，也无疏漏，所以对包含大量中、高频信息的信号能够进行更好的时频局部化中、高频信息的信号能够进行更好的时频局部化分析。分析。资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值小波包的定义正交小波包正交小波包 的一般解释的一般解释:本章仅考虑实系数滤波器本章仅考虑实系数滤波器.为便于表示小波包函数，本章引入以下新的记号：为便于表示小波包函数，本章引入以下新的记号：资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值小波包的定义通过通过 在固定尺度下可定义一组称为小波包的函数。在固定尺度下可定义一组称为小波包的函数。由由 递归定义的函数递归定义的函数 称为由称为由正交尺度函数正交尺度函数 确定的确定的小波包小波包。资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值小波包的性质(习题10.1)的傅立叶变换可以由的傅立叶变换可以由 性质性质10.1 表示。表示。性质性质10.2 具有平移正交性，即具有平移正交性，即 性质性质10.3 与与 之间具有正交性，即之间具有正交性，即 小波包的平移系小波包的平移系 性质性质10.4 构成构成 的一组规范正交基的一组规范正交基 资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值小波空间的精细分割 小波分析存在的不足：小波分析存在的不足：随着随着j的增大，的增大，相应小波基函数相应小波基函数 分辨率越高，而其频谱的局部性变得越差即频谱分辨率越粗。分辨率越高，而其频谱的局部性变得越差即频谱分辨率越粗。的空间局部性越好的空间局部性越好即空间应对措施：应对措施：对小波空间对小波空间Wj做进一步分解做进一步分解小波空间的分解：小波空间的分解：令令 表示由小波包表示由小波包 的二进伸缩和平移的二进伸缩和平移 的线性组合生成的的线性组合生成的 的闭子空间，则的闭子空间，则 资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值小波空间的精细分割 小波空间的分解：小波空间的分解：对于每个对于每个 且对给定的且对给定的，及 函数系函数系 是空间是空间 的一个规范正交基。的一个规范正交基。频带划分性质频带划分性质:小波包具有划分较高频率频带的能力，可得到比较好小波包具有划分较高频率频带的能力，可得到比较好的频率局部化。的频率局部化。资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值一个逼近空间的小波分解及小波包分解 小波分解小波分解 小波包分解小波包分解 小波包正交基小波包正交基:在图在图10.2b的二分树上取一组子空间集合，如果其直和恰的二分树上取一组子空间集合，如果其直和恰能将能将 V3空间覆盖，相互间又不重叠，则这组空间集合的正交规范基便组空间覆盖，相互间又不重叠，则这组空间集合的正交规范基便组成一个小波包正交基。成一个小波包正交基。资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值小波包滤波器组 已知已知:长度为长度为离散输入信号离散输入信号 的均匀采样的的均匀采样的，首先将，首先将 与在尺度与在尺度 下的一个逼近函数下的一个逼近函数 联系起来，联系起来，它的分解系数它的分解系数 满足：满足：问题问题:已知已知f(t)在子空间在子空间 的子空间的子空间 上的小波包系数，计算出上的小波包系数，计算出在在 的两个子空间的两个子空间 和和 上的小波包系数上的小波包系数.?资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值小波包滤波器组 小波包分解算法：小波包分解算法：小波包重构：小波包重构：小波包多级分解与多级重构：小波包多级分解与多级重构：资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值最佳小波包基的选取 称函数族称函数族 尺度函数尺度函数 为正交为正交导出的导出的小波库小波库。1.Mallat正交小波基正交小波基 2.3.结论结论:小波库中包含许多规范正交基即小波包基小波库中包含许多规范正交基即小波包基.问题问题:什么是最佳小波包基什么是最佳小波包基?如何从小波库中快速选取如何从小波库中快速选取?资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值最佳小波包基的选取 信息代价函数信息代价函数把信号把信号 在一个正交小波包基下展开，使得它与一个小波包系数在一个正交小波包基下展开，使得它与一个小波包系数 序列序列对应对应,我们在该序列上定义一个信息代价函数我们在该序列上定义一个信息代价函数 M,它满足它满足 如下两个条件：如下两个条件：（1）可加性条件）可加性条件（2）代价函数）代价函数M的取值应该反映系数的集中程度的取值应该反映系数的集中程度.最佳小波包基最佳小波包基对于一个给定信息代价函数对于一个给定信息代价函数M,小波包基小波包基 B称为信号称为信号f(t)相对于该代价函数相对于该代价函数 的最佳基，如果在的最佳基，如果在 的所有小波包基中，的所有小波包基中，f(t)在小波包基在小波包基 B下对应的对应的小波包系数序列具有最小的信息代价值小波包系数序列具有最小的信息代价值 资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值最佳小波包基的选取 常用的一些信息代价函数：常用的一些信息代价函数：（1）幅值大于某阈值的系数个数）幅值大于某阈值的系数个数（2）范数的集中度（范数的集中度（concentration）（3）对数熵）对数熵，约定，约定（4）信息熵）信息熵资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值最佳小波包基的选取 在一般情况下，具有最小代价函数值的序列不易计算出来。所幸的是，在一般情况下，具有最小代价函数值的序列不易计算出来。所幸的是，正如正如10.3节所谈到的，在实际应用中我们通常考虑的是节所谈到的，在实际应用中我们通常考虑的是 L2(R)的一个子的一个子空间的小波包分解空间的小波包分解,这种分解可以用一个小波包二叉树表示这种分解可以用一个小波包二叉树表示.我们可以我们可以采采取自底向顶的快速搜索法发现最佳小波包基。取自底向顶的快速搜索法发现最佳小波包基。资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值最佳小波包基的选取(1)(2)(3)(4)资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值最佳小波包基的选取 从以上讨论可知，从以上讨论可知，最佳基搜索算法主要由两步组成：最佳基搜索算法主要由两步组成：1）搜索构成）搜索构成最佳基的节点；最佳基的节点；2）抽取离树根最近的最佳基节点中的小波包系数。）抽取离树根最近的最佳基节点中的小波包系数。顺便指出，如果小波包分解采用深度优先顺序（顺便指出，如果小波包分解采用深度优先顺序（depth-first order），则最佳基节点的标记过程可以在计算节点中小波包系数），则最佳基节点的标记过程可以在计算节点中小波包系数的同时完成。由于小波包树具有有限深度，所以以深度优先的搜索的同时完成。由于小波包树具有有限深度，所以以深度优先的搜索算法可在有限步终止。算法可在有限步终止。如果原信号的长度为如果原信号的长度为，则最佳基算法的计算复杂度为，则最佳基算法的计算复杂度为 资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值小波包变换的应用 几种不同变换对应的时频平面铺砌几种不同变换对应的时频平面铺砌 资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值小波包变换的应用 信号小波包分析的基本实现步骤如下：信号小波包分析的基本实现步骤如下：1）选择适当的小波滤波器，对给定的采样信号进行小波包变换，）选择适当的小波滤波器，对给定的采样信号进行小波包变换，获得树形结构的小波包系数。获得树形结构的小波包系数。2）选择信息代价函数，利用最佳小波包基选取算法选取最佳基。）选择信息代价函数，利用最佳小波包基选取算法选取最佳基。3）对最佳正交小波包基对应的小波包系数进行处理。）对最佳正交小波包基对应的小波包系数进行处理。4）对处理后的小波包系数采用小波包重构算法得到重构信号。）对处理后的小波包系数采用小波包重构算法得到重构信号。资金是运动的价值，资金的价值是随时间变化而变化的，是时间的函数，随时间的推移而增值，其增值的这部分资金就是原有资金的时间价值小波包变换的应用 小波包在信号去噪、滤波等方面的应用原理和方法小波包在信号去噪、滤波等方面的应用原理和方法 （1）滤波与去噪）滤波与去噪（2）非平稳机械振动信号的故障诊断）非平稳机械振动信号的故障诊断（3）特征提取）特征提取注意注意:习题习题10.1与与10.2可以作为作业题选做可以作为作业题选做.
