MRA小波分析.docVIP

3.1 haar小波分析 3.1.4 haar小波重构算法 由如下信号开始： 应用等式： (1) (2) (3) (4) 利用(1), 用代替,则 可化为: 类似的, 利用(2), 用代替,则 可化为: 这里， 结合和有： 这里 接下来，把用相同的方式加到中： 得到： 这里 注意：系数和决定了，和决定了……. 定理（Haar重构）设 这里 那么 这里 3.1.4滤波器实现 分解和重构算法可用离散滤波器和作用于序列的简单算子来表述，同时算法还可以方块图形象地表示。 首先考虑分解算法。定义两个离散滤波器（卷积算子）和，其冲击响应分别为和： 若，那么 假如只保留偶数下标项（丢弃奇数项系数，下采样），则： 应用以上表示法，由阶尺度系数得到阶尺度和小波系数。有： 图(分解算法)：下采样算子用更常用的符号“”代替 再考虑重构算法。定义两个离散滤波器和，其冲击响应分别为和： 对，有 若、是两个奇部为0的序列，那么 把以上两式相加，得到： 与重构公式对比 令，可得到： 序列和称为对序列和进行了上取样，则重构公式可表示为： ,,,分别为低通、高通、低通重构、高通重构滤波器 4 多分辨率分析MRA Haar小波方案（分解和重构）依赖于尺度函数和小波函数。 优点：易描述，算法简单； 缺点：两个函数均不连续，仅为连续函数提供了一个粗略近似。 Mallat提出了MRA理论，其框架与Haar小波方案类似，而构筑块和是连续的 多分辨率分析MRA的一般定义： 令为中的函数子空间序列。若以下条件成立，则空间集合称为依尺度函数的多分辨率分析： 1.（嵌套性） 2.（稠密性） 3.（分立性） 4.（尺度性） 当且仅当 5.（标准正交基）函数，且是的标准正交基。 定义注解： 注1：表明分辨率时的逼近包含计算较低分辨率时的逼近应有的全部信息； 注2：“稠密性”意味着当时，分辨率越来越高，包含越来越多信息，对的逼近程度越高； 注3：称为近似空间，根据近似空间的不同，有多种选择，从而对应不同的多分辨率分析； 注4：haar尺度函数是紧支的，与DB小波函数相关的尺度函数不仅紧支，而且连续，相应的分解和重构算法速度快，对分析和重构信号特别有利。 例：线性样条MRA：线性样条是连续的、逐段线性的函数。 令是所有能量有限信号所组成的空间，是连续且分段线性的，其拐点只可能在二进点处（如图）。这些近似空间满足多分辨分析的定义。尺度函数是“跳转函数”： 易证平移系是的基，但不正交，但可以通过一定的方法使其标准正交。（见《小波与傅立叶分析基础》5.3节） 易见： MRA性质1：设是一个依尺度函数的多分辨率分析。那么对任意，函数集是的一个标准正交基。 证明：为证明张成，必须先证明对任一可以写成的线性组合。设任意，应用尺度特性（MRA中条件4），函数，因此是的一个线性组合。以取代，即可发现是的一个线性组合。 为证明是标准正交的，只需证明 或 进行变量替换，则有 根据MRA中条件5给出的标准正交基性质，上式等于 MRA 性质2：多分辨率分析的中心公式——尺度关系式 定理1 设是一个依尺度函数的多分辨率分析，那么下列尺度关系式成立： （1） 而且还有 （2） 或者，等价于 这里 注1：设为实函数可以两边与做内积得到： 注2：(2)式只需将代为即可； 注释： 1：上式表示了与的关系，称为尺度关系式或双尺度关系； 2：当足够大时，的支撑将在的支撑之外，当的支撑为紧支撑时，则仅有有限个非零，所以定理中的和式是有限的。 3：haar系中为： 4：Db2中相应的为： 5：特别注意与无关。 MRA 性质3：（尺度关系中满足的关系）：设为一个具有尺度函数的正交多分辨率分析，则下列等式成立： 1 2 3 4 证（1）：(简单起见，可视为实数) 证明参考《小波与傅立叶分析基础》。 小波和小波空间的 构造： 一旦确定下来，利用尺度关系即可构建产生的函数. 定理2 设是一个多分辨率分析，相应的尺度函数为 （为定理1中的系数）。令是由张成，这里 那么是中的正交补。而且是的一个标准正交基。 多分辨率分析MRA的后续内容（分解和重构算法）请大家参考《小波与傅立叶分析基础》P177--P191。 总结： 分解公式：