zoukankan      html  css  js  c++  java
  • Multiresolution Analysis(多分辨率分析)

    [注意:本文中所有的傅里叶变换和反变换均含对称因子$frac{1}{sqrt{2pi}}$,且$z=e^{-ikomega}$]

    1. 多分辨率分析

    1.1 概念

    多分辨率分析指的是一系列$L^2(R)$的子空间$V_j$,每个子空间$V_{j+1}$都是它“上一级”子空间$V_j$的“精细化”:

    (i) $V_jsubset V_{j+1}$    # 子空间逐级嵌套

    (ii) $overline{U_{j in Z}V_j}=L^2(R)$  # 所有子空间的并“构成”$L^2(R)$。注意此处上横线不是求补的意思(待补充:)

    (iii) $cap_{j in Z} V_j={0}$    # 全部子空间的交为空集。可这样理解:交的结果趋向于最“不精细”的那个子空间$V_{-infty}$

    (iv) $f(t) in V_0 iff f(2^jt)in V_j$ # 子空间之间存在“缩放”关系

    对子空间$V_0$,存在一组尺度函数${ phi (t-k) }$作为规范正交基。

    我们假定$phi(t)$是实值且归一化的函数,那么下式成立:

    $int_Rphi(t)dt=sqrt{2pi}Phi(0)=1$   # 傅里叶变换式取$omega=0$即得

    根据子空间之间的“包含”和“缩放”关系,可以写出空间$V_j$的规范正交基:

    $phi_{j,k}(t)=2^{j/2}phi(2^jt-k)$

    系数$2^{j/2}$使得时间尺度变换后的基的范数仍为1。

    如果函数$f(t) in sum_{k in Z}c_kphi_{j,k}(t)$,那么可以写出$f(t)$在$V_j$的规范正交基上的投影:

    $c_k=langle f(t),phi_{j,k}(t) angle$

    因为各级子空间存在包含与被包含的关系,所以子空间$V_j$的基总可以用$V_{j+1}$的基的线性组合表示。例如:

    $phi (t) = sqrt{2} sum_{k in Z} h_k phi (2t-k)$

    上式中,$h_k=langle phi(t), phi_{1,k}(t) angle$

    对于更一般的$phi_j(t)$:

    $phi_j(t)=sum_{k in Z} h_k phi_{j+1,k}(t)$

    注意对上式中的各级子空间$V_0, V_1, ..., V_j, V_{j+1}, ...$,$h_k$总是相同的。

    对于更更一般的$phi_{j,k}(t)$,将$h_k$进行移位即可:

    $phi_{j,l}(t)=sum_{k in Z} h_{k-2l} phi_{j+1, k}(t)$

    考虑到各级子空间之间的尺度关系,$h$的下标为$k-2l$是自然的事。

    系数$h_k$称为尺度滤波器,具有如下性质:

    (i) $sum_{k in Z} h_k=sqrt{2}$    # 为了保持归一化,$sum h_k$必须等于$sqrt{2}$。回想下HAAR的尺度滤波器系数

    (ii) $sum_{k in Z} h_k h_{k-2l} = delta(l)$  # 因为$phi(t)$是正交基,所以距离为$2l$的权序列不相关

    (iii) $sum_{k in Z} h_k^2=1$  # 为了保持尺度变换前后(尺度滤波器内的部分)能量不变

    如果将$L^2(R)$上的任意函数$f(t)$投影到$V_j$:

    $P_{f,j}(t)=sum_{k in Z}langle f(t), phi_{j,k}(t) anglephi_{j,k}(t)$

    如果我们有一个函数$f_{j+1}in V_{j+1}$,那么有:

    $f_{j+1}(t)=sum_{k in Z}a_kphi_{j+1,k}(t)$

    现在将$f_{j+1}(t)$投影到$V_j$:

    $f_j(t)=sum_{lin Z}b_lphi_{j,l}(t)$

    为了求出$b_l$的表达式,我们可以这样考虑:

    想象一下将$phi_{j+1}$投影到$V_j$,每一个$phi_{j+1}$都要用$phi_j$的线性组合来表示。所以$phi_{j,l}$前的系数是以$a_k$为权的$h_k$的线性组合(加上偏移量$2l$):

    $f_j(t)=sum_{lin Z}b_lphi_{j,l}(t)=sum_{lin Z} left(  sum_{k in Z}a_kh_{k-2l}   ight) phi_{j,l}(t)$

    1.2 小波函数

    $psi(t)=sqrt{2}sum_{k in Z}g_kphi(2t-k)$

    $psi_{j,k}(t)=2^{j/2}psi(2^jt-k)$

    $g_k=(-1)^kh_{1-k}$

    $psi_{j,l}(t)=sum_{k in Z}g_{k-2l}phi_{j+1,k}(t)=sum_{k in Z}(-1)^k h_{1+2l-k}phi_{j+1,k}(t)$

    1.3 $phi(t)$的symbol

    上面我们将$phi(t)$表示为

    $phi(t)=sqrt{2}h_kphi(2t-k)$

    若将$h_k$记作$h(k)$,上面的式子就是一个卷积式子。于是经过一些简单变换就得到:

    $Phi(omega)=H(frac{omega}{2})Phi(frac{omega}{2})$

    $H(omega)=frac{1}{sqrt{2}}sum_{k in Z}h_ke^{-ikomega}$

    $phi(t)$具有如下性质:

    (i) $||Phi||=1$  # 单位长度

    (ii) $Phi(0)=frac{1}{sqrt{2pi}}$ 

    (iii) $H(0)=1$ # 直流增益1

    (iv) $H(omega)$周期为$2pi$  # $h$是离散的,所以这条显然

    1.4 The Stability Function

    根据空间$V_j$的基在时域的正交性质,可以推导出下式(The Stability Function):

    $mathcal{A}(omega)=sum_{l in Z} left|  Phi(omega + 2pi l) ight| ^2=frac{1}{2pi}$

    以HAAR小波为例,其频域为sinc函数,按$2pi$移位平方累加,在整个频域为常数。

    另一方面,在时域上有下式成立:

    $sum_{k in Z}phi(t-k)=1$

    The Stability Function只需$phi(t)$为正交基即可,不需要尺度条件。如果加上尺度条件,还有下面的式子:

    $|H(omega)|^2+|H(omega+pi)|^2=1$

    如果一个函数的频谱满足上式及以下条件,那么该函数有与之对应的尺度函数:

    (i) $H(0)=1$    # $h_k$过直流

    (ii) 满足$H(z)=(frac{1+z}{2})^NS(z)$,其中$max_{|z|=1}|S(z)|le 2^{N-1}$ 

        # 1. 系数$(frac{1+z}{2})^N$使得$H(z)$为低通形式;

        # 2. [1, 1]的z变换为$frac{1+z}{2}$(考虑了归一化系数),所以$(frac{1+z}{2})^N$可看作N组[1, 1]的卷积的z变换。

    1.5 $g(t)$的symbol$G(omega)$

    $G(omega)=frac{1}{sqrt{2}} sum_{k in Z} g_k e^{-ikomega}=frac{1}{sqrt{2}}sum_{k in Z} (-1)^kh_{1-k}e^{-ikomega}$

    $G(omega)$和$Psi(omega)$的性质

    (i) $Psi(omega)=G(frac{omega}{2})Phi(frac{omega}{2})$

    (ii) $G(omega)=-e^{-iomega}overline{H(omega+pi)}$

    (iii) $G(0)=0$  # 高通效果

    (iv) $sum_{k in Z}g_k=0$  # 直流增益0

    (v) $sum_{k in Z}h_{2k}=sum_{k in Z}h_{2k+1}=frac{sqrt{2}}{2}$

  • 相关阅读:
    面向对象设计原则
    简单工厂模式和策略模式结合使用php
    lua string
    js.ajax优缺点,工作流程
    深入理解JavaScript是如何实现继承的
    js中哈希表的几种用法总结
    js,indexOf()、lastIndexOf()
    js获取字符串字节数方法小结
    JS(JavaScript)插入节点的方法appendChild与insertBefore
    js中的this关键字详解
  • 原文地址:https://www.cnblogs.com/byeyear/p/7375929.html
Copyright © 2011-2022 走看看