[注意:本文中所有的傅里叶变换和反变换均含对称因子$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}$