高等数字信号处理 吴正国 尹为民 第6章4新.ppt

第四节 小波包 2、双正交小波包 ②性质 ③平方可积空间分解 算子的性质： 注意：是空间”合”的关系，而非”直和”。 ④双正交小波包算法 重构式 分解式 注意：重构与分解应分别用互为对偶的函数 3、最优小波包基的选择 按不同信号分析的要求，如何从小波包库中 选择一组最优的小波包基。 ①最优代价函数(以数据压缩为例) 代价函数为关于序列 x (n ) 的实函数M ( x ( n )) 信息熵： 对数熵： 具有可加性： 集中度最大为优，对于上述代价函数则是最 小化为优。 对于信号分类，则可选择描述两序列相对信 息含量的相对熵；对于子带分解，可选择各子 帶间的互相关函数为代价函数等。 ②快速搜索算法-----单树算法 在各子空间内，以小波包系数为元素计算M ③时变小波包分解算法------双树算法 对于时变信号，在不同的时间段取不同的 小波包分解方案。 实现步骤：假定待分解信号按时变性质可分 为A、B、C、D四段，则可首先按单树算法求最优子空间分解方案，然后计算分为不同段组合的最优子空间分解方案，再在此基础上形成第二个 树---时间树，最后在时间树上搜索最优小波包基 注意在每个时 间分段的最优 基搜索中，信 号长度的二进 制缩短是子空 间分解(高、低 通滤波)的结果。 第五节 小波提升方案 (lifting scheme) 第一代小波：构造方法依赖富里叶变换，时域 上依赖 伸缩平移变换；算法相对复杂，且有有限字长效应问题。 第二代小波：(90年代中期提出) 不依赖富里叶变换，直接从时域上实现小波构造(尤其是按予定要求构造双正交小波)；可提供一种比mallat算法更快速的小波级数分解和重构算法；实現整数小波变换，并能实現同址处理。 一、提升策略的基本原理 1、Vetterli-Herley定理：对于所给定的双正交尺 度函数 和 ，其对应的双尺度序列为 和 ，那么所有满足如下条件的序列 与 也将满足双正交关系。反之亦然。(下式中 为一个三角函数多项式。 * ■ ▲ ■ * 为提高高频段频率分辨率，而对小波空间再进行 多尺度剖分 一、正交小波包的定义与性质 1、定义：给定正交尺度函数φ( t ) 及正交小波 函数ψ( t ) ，若令 则双尺度方程可记为 定义函数序列： 称此函数序列为由φ( t ) 生成的小波包。 2、正交小波包的频潽 若令 对s0位的二进制整数n可记为 3、正交小波包的性质(为简便计，假定实函数) 用数学归纳法可证： 对于s0+1位二进制整数n，不超过0.5n的最大整数可记为 正交性对于n=0或1时皆成立，所以使用数学归纳法。即令上式对于二进制位数不超过s0位的所有整数成立，需证明该式对于二进制位数不超过(s0+1)位的整数n成立 n相同时正交小波包函数整数平移正交 n为相邻奇偶数时正交小波包函数整数平移正交 二、空间的正交小波包分解 1、空间算子 完全重构(c=1,k=0) 单位算子 2、 j 尺度空间 定理：如果上述空间序列对n大于等于零、j为全体整数成立，则有 证明：第一式是小波包性质2的结果，仅需证第二式， 即证合空间的基为两子空间基的线性组合。 3、平方可积空间的正交小波包分解 kj=0时Wj不分解，kj 表示分解的程度 随着尺度的增加，频窗半径減小，频率分辨 率增加，但时间分辨率下降。 4、小波库与小波包基 称为由φ( t ) 所导出的小波包库。 n 为振荡参数 j 为尺度参数 l 为平移参数 从小波库中抽取的、能构成平方可积空间 的一组正交基称之小波包基。 选择不同的 组合，将给出 不同的小波包基。 三、小波包变换及快速算法 1、正交小波包变换 从小波库中选择一组最佳小波包基，则有 ①正交小波包变换 为 f ( t ) 在 子空间的投影系数， f ( t )在上式各子空间投影系数的集合，称之 f ( t ) 的正交小波包变换。 ②正交小波包变换算法 定理：如果令 则有 ③正交小波包逆变换 ①定义