高等数字信号处理吴正国尹为民第6章3课件教学.pptVIP

第三节 小波级数的快速算法 一、二进尺度的小波级数计算 二、小波级数计算中的边界延拓 ⒉周期延拓及矩阵表示(以简单周期延拓为例) 利用周期信号的小波变換特性來研究有限长度信 号小波级数计算中的周期延拓问题。令周期序列为 记 的长度皆为 L ，则对于 N 奌输入 信号的 Mallat 算法公式变为 若令 对于 0≤ j ≤J 成立，则上式 中对每一个 k 仅有 L 次非零乘法，因此上式可记为如 下矩阵形式 考虑小波级数的周期性，将相同项合并后得： 用矢量和矩阵形式可记为 令 ，它与 皆为 维方阵 三、基于Mellin变换的快速算法 利用Mellin变换的性质，对于某给定时刻，一次计 算N个(按等比数列变化的)尺度下的小波级数。 ① Mellin变换的定义及性质 定义： 可证明逆变换为： 性质： 对变量c的富氏变换 时域的尺度变換等效变换域的线性相移(m不一定为整数) 乘积定理 ②算法的理论依据： 富氏变换实现 变量离散化 取样： ③算法的讨论 对给定的平移参数，每次计算N亇尺度的小波级数 尺度： 要求信号按几何级数增长的非均匀取样，第n个取样 时刻为： ，取样间隔为： 。对于给定的信 号均匀取样，应通过插值法获得非均匀取样。 由于取样间隔随n的增加而不断增大，从而产生信号 混迭失真。因此限制q不能太大，即所计算的尺度跨度 不能太大。 该算法的优奌是可直接利用FFT算法，速度较快。 在线教务辅导网： 更多课程配套课件资源请访问在线教务辅导网 * ■ ▲ ■ * 假定信号长度为N，扩張函数长度为 K ( j )，M为 两者中最大者 ，随着 j 的增加，K ( j ) 将迅速增加。 从而使 M 迅速增加，即使用 FFT 算法，计算所有 j 尺 度的小波级数的计算量也是很大的。 对于二进尺度的小波级数计算可使用 mallat 算法， 由于二倍再取样的结果，后面将证明，用该算法计算 所有 二进尺度的小波级数其乘法次数只需2NL ( L为双 尺度序列长度，一般 L远小于N )。 Mallat 算法只解决了 DWT 的计算问题，为了计 算连续函数 f ( t ) 的小波级数，还需计算 ⒈ Shensa 算法:利用Mallat算法计算小波级数 ①算法推导 时域取样间 隔 滿足取样定 理条件 可以由 f ( n ) 通过滤波器 q ( n ) 而得到 ②算法框图 ③ q ( n )与对偶尺度函数 的关系 滤波器的频率响应为 ( -π，π)内对偶尺度函数 频谱周期延拓的结果。(ω为模拟频率) 但滿足取样 定理条件 为数字频率 若尺度函数频谱的最大非零频率 则 若尺度函数频谱的低通带宽充分大(或b0充分大)，使得频谱在积分范围内近似为1，则有 即只有在这种情况下，可用信号的取样值直接进行小 波级数计算。在保证 f ( t ) 尽可能小的取样失真条件 下，取样间隔不可能太大。但是，对于与小波函数对应的 尺度函数而言，由于其时窗充分的小，其低通帶宽是可能 充分大，从而使上式成立。因此，对于实际应用的小波变 换而言，通常令信号的取样值(取样率滿足取样定理要求) 为 由于T充分小，对于充分大的ω仍有下式成立： ⒉ A Trous 算法 ①基本思想：通过使用特殊的尺度函数，达到直接使 用 f ( n ) 进行小波级数计算的目的。 使用具有无限低通带宽的尺度函数的特殊尺度函数， 即有下述算法 若尺度函数φ( t ) 滿足 则 ②公式推导 为多孔 ( a trous )滤波器，h(n) 的偶数点除零外皆为零。 定义序列g(n)满足下式(注意此式不是双尺度方程)： 据Mallat算法 由于小波函数不满足双尺度方程，需重新推导 的递推式。 据小波变换定义得： 引用g(n)的定义式及h(n)的双尺度方程可导出下式： 类推得 实施变量代换l=k-m ③计算框图 g (n )是给定小波函数的取样，为帶通滤波器，而不是 高通滤波器； h ( n )为无限序列，它的偶数奌 ( 除零外 )固定为零；但奇数点可选 由 f ( n ) 直接求小波级数 g ( n )与h ( n )不满足正交关係 φ( t –k ) 巳不是正交基 ⒊算法使用中应注意的问题 ①计算复杂度:以 Mallat 算法为例 假定 的长度皆为L，总分解Ｊ个尺度， 则对于每个输入奌的实数乘法次数为 输入序列长度为 ，则乘法次数为 ②时间间隔加密 在Mallat算法中，由于二倍再取样运算，使小波 级数的时间间隔 随尺度 j