第十二章 离散小波变换.ppt

理想的小波包时—频空间分解 最优小波基 针对每一类待分析信号，选定小波包，信号由不同的子空间组合构成。 正交小波包变换 第j级的第n个子空间 第j级的第n个子空间奇偶分解 当n=0,j=1 例，HARR小波包的空间分解 第0层只有一个子空间 第1层有2个子空间 第2层有4个子空间 正交小波包变换的滤波器实现结构 小波包的空间分解的滤波器实现结构 小波包的空间综合的滤波器实现结构 [例12-9-1]已知信号由幅值等于1，频率分别为0.2 KHz、0.6 KHz、1.0 KHz、1.4 KHz、1.8 KHz、2.2 KHz、2.6 KHz、3.0 KHz的正弦波构成，采样频率等于6.4 KHz，采样长度等于0.08秒，试用db20小波对其进行三层小波包分解，绘出第三层各子空间系数恢复出的信号波形，以及对应的幅频特性，并说明第三层各子空间对应的频谱区间 [例12-9-1]小波包分解第三层的8个子空间系数的恢复信号及其频谱 小波包分解结果表明小波包频谱划分并不是按子空间序号顺序递增 第三层子空间d0至d7的频谱分布为： [0，1/8]，[1/8，2/8]，[3/8，4/8]，[2/8，3/8]， [7/8，8/8]，[6/8，7/8]，[4/8，5/8]，[5/8，6/8] 第三层子空间对应的频谱剖分排列为：0\1\3\2\7\6\4\5 小波包分解的频谱特性 小波包消噪和压缩算法过程 据信号特性和分析要求选择合适的小波包函数确定分解层数； 对采样信号进行小波包分解； 根据给定的熵标准，计算最优树 采用一定的算法，对最优树中的每个小波系数进行阀值处理 应用阀值处理后最优树中的每个小波系数重构信号，得到消噪结果处理 对阀值处理后最优树中的每个小波系数进行量化编码，得到信号的压缩结果 [例12-9-2]采用相同的小波函数、分解层数，以及相同的阀值算法，应用小波包变换对例12-8-2的信号进行消噪处理，并比较两种方法的消噪结果 [XD,TREED,PERF0,PERFL2] = wpdencmp(X,SORH,N,wname,CRIT,PAR,KEEPAPP) V3 V0 V1 V2 多分辨率分析的空间示意图 小波空间和小波函数 尺度空间之间相互包含不具有正交性 定义尺度空间的正交补空间 称 是小波空间 小波空间示意图 V1 W2 W3 V3 V0 V2 W1 正交基存在性： 小波空间的特性 正交性 逼近性 尺度伸缩规则 W1空间的一组正交基 一个多分辨率分析对应一个小波函数 信号的多分辨率分析 对 空间按如下方式分解 计算信号x(t)在各空间下的投影 当分解尺度J趋向无穷大，信号全部分解在小波空间 小波级数分析时框架常数A=B=1时信号重建公式，与上式完全相同。MRA分析从空间二进剖分角度得到了信号的正交小波分解 尺度函数和小波函数的性质 不同位移的两个尺度函数正交 小波函数的伸缩平移函数对尺度和平移参数正交 Wj和Vj空间正交 12．6 Mallat算法和双通道滤波器组 MRA分析中的双尺度方程 V0空间的一次分解 V0空间正交基 V1空间的尺度函数 在V0空间展开 W1空间的小波函数 在V0空间展开 MRA分析的双尺度方程 Vj空间的分解 MRA分析的表示方法 MRA分析的表示方法 信号在MRA分解空间上的投影 设信号在Vj-1空间的投影是: Vj空间的投影与Vj-1空间投影的关系 Wj空间的投影与Vj-1空间投影的关系 MRA投影过程 MRA分解的滤波器实现 Vj空间的投影与Vj-1空间投影关系的滤波器结构 Wj空间的投影与Vj-1空间投影关系的滤波器结构 MRA分析分解过程的完整滤波器实现结构 信号空间分解的重建过程 设信号在空间Vj的投影序列是: 设信号在空间Wj的投影序列是: 设信号在空间Vj-1的投影序列是: 空间Vj的投影和Wj空间投影重建Vj-1空间投影 空间Vj的投影和Wj空间投影重建Vj-1空间投影的滤波结构 序列 经过 滤波输出 序列 经过 滤波输出 空间Vj的投影和Wj空间投影重建Vj-1空间投影的完整滤波结构 信号空间分解的滤波器实现结构---Mallta算法 12．7正交小波和小波滤波器 MRA分析的表示方式 尺度函数的小波函数(适合CWT和二进WT) 高通滤波器g(n)和低通滤波器h(n),适合离散正交小波变换 通过滤波器组的精确重建理论寻找小波函数 MRA分析的一级空间分解和一级空间综合 一级空间分解和一级空间综合的滤