【小波分析】四、正交多分辨分析

【小波分析】四、正交多分辨分析

内容回顾

形而上的理解

上一次我们引入了多分辨分析，及其与正交小波的一个关系。

多分辨分析 $V_j$ 有这样一些特征。$V_j$ 包含于 $V_{j+1}$ 里面的，它可以通过在时间上压缩一倍直接得到。而每个 $V_j$ 空间，都有自己的一组标准正交基，$V_{j+1}$ 的标准正交基，并不是通过 $V_j$ 的标准正交基直接扩充得到，两组基之间，可能完全变了样。它最小的是一个只包含零元素的平凡空间，$V_j$ 中最大的一个，刚好充满整个 $L^2$ 空间。这个过程就像盖房子，你要盖更大的房子，你就要把原来的房子推倒重来。

而正交小波有这样一些特征。$W_{j+1}$ 也是可以通过 $W_j$ 在时间上压缩一倍得到，更神奇的是，$W_{j+1}$ 和 $W_j$ 是正交的，这就好比，一个东西，你把在时间上压了一压，就变得不是自己了，反而是和自己正交的东西。每个 $W_j$ 都有一组标准正交基，基和基之间也是正交的，$W_{j+1}$ 的基可以通过 $W_j$的基时间伸缩和单位化得到。把这些所有的 $W_j$ 基放在一块，刚好可以张成整个的 $L^2$ 空间。这个过程就好比是贴瓷砖，每次都在原来的瓷砖周围贴上一层新的瓷砖，一直到铺满一面墙。

$V_j$ 和 $W_j$ 之间有这样一个关系。

$$V_{j+1}=V_j\oplus W_j$$

所以，从这里来看，如果给定了一个正交多辨分析，只要找到 $V_{j+1}$ 在 $V_j$ 上的补空间的一组标准正交基，使得它和前面所有的 $W_j$ 是正交的，那么，我们其实就找到了正交小波。我们可以把这个过程写得更加具体而可操作，那么这就是正交多辨分析的主要内容。

多分辨分析与正交小波

回顾一下多分辨的定义：

定义 设 { V j ; j ∈ Z } \left\{V_{j} ; j \in Z\right\} {Vj​;j∈Z} 是$L(R)$上的一列闭子空间, $\varphi (x)$ 是 $L^2(R)$ 中 的一个函数，如果它们满足如下的五个条件，即
(1) 单调性:
$$V_j\subset V_{j+1},\forall j\in Z$$
(2) 唯一性:
⋂ j ∈ Z V j = { 0 } \bigcap_{j \in Z} V_{j}=\{0\} j∈Z⋂​Vj​={0}
(3) 稠密性:
$$\overline{\left(\bigcup _{j\in Z}{V}_j\right)}=L^2(R)$$
(4) 伸缩性:
$$f(t)\in V_j\iff f(2t)\in V_{j+1}\forall j\in Z$$
(5) 可构造性:
{ ϕ ( t − n ) ; n ∈ Z } \{\phi(t-n) ; n \in Z\} {ϕ(t−n);n∈Z}
构成子空间 $V_0$ 的标准正交基。

那么，称 { { V j ; j ∈ Z } ; ϕ ( x ) } \left\{\left\{V_{j} ; j \in Z\right\} ; \phi(x)\right\} {{Vj​;j∈Z};ϕ(x)} 是 $L^2(R)$ 上的一个正交多分辨分析(MRA，Multi-Resolution Analysis)。

由如上的定义，我们容易知道：
$$\left\{{\varphi }_{j,n}(t)=2^{\frac{j}{2}}\varphi \left(2^{j}t-n\right);n\in Z\right\}$$
构成了$V_j$ 空间的一组标准正交基。

仿照 Shannon 小波的构造方法，对 $\forall j\in Z$, 定义如下的子空间 $W_j$，
$$W_j\perp V_j,V_{j+1}=W_j\oplus V_j$$

容易验证，子空间序列 { W j ; j ∈ Z } \left\{W_{j} ; j \in Z\right\} {Wj​;j∈Z} 具有下述性质：

• $\forall j\ne l,W_j\perp W_l$;
• $L^2(R)=\underset{l\in Z}{\oplus }{W}_l$;
• $\forall j\in Z,g(t)\in W_j\iff g(2t)\in W_{j+1}$

我们要构造正交小波，我们只要找到一组标准正交基，要得到一组标准正交基，由第一条和第二条我们知道，只要找到所有 $W_j$ 的一组标准正交基就可以了，再由第三条，我们其实只要找到 $W_0$ 的一组标准正交基 $\psi (t-k)$ 就可以了。这是因为如果 $\psi (t-k)$ 是 $V_0$ 的标准正交基，那么

$$\left\{2^{\frac{j}{2}}\psi \left(2^{j}t-k\right);(j,k)\in Z\times Z\right\}$$

必然是$W_j$ 的标准正交基。

所以，我们要构造正交小波，要做的无非就是找到 $W_0$ 的一组标准正交基。

尺度方程和小波方程

尺度方程

由于 $\varphi (x)\in V_0\subseteq V_1$, $V_1$ 有标准正交基 { 2 ϕ ( 2 t − n ) ; n ∈ Z } \{\sqrt{2} \phi(2 t-n) ; n \in Z\} {2 ​ϕ(2t−n);n∈Z}，必存在唯一的系数序列 { h n ; n ∈ Z } ∈ l 2 ( Z ) \left\{h_{n} ; n \in Z\right\} \in l^{2}(Z) {hn​;n∈Z}∈l2(Z), 使得
$$\varphi (t)=\sqrt{2}\sum _{n\in Z}{h}_n\varphi (2t-n)$$
这个方程叫做尺度方程。系数计算的方式为，
$$h_n=⟨\varphi (t),\sqrt{2}\varphi (2t-n)⟩=\sqrt{2}{\int }_R\varphi (t)\bar{\varphi }(2t-n)\mathrm{d}x$$
这个系数叫做低通滤波器系数。对尺度方程做傅里叶变换可以得到，
$$\hat{\varphi }(\omega )=\mathrm{H}\left(\frac{\omega }{2}\right)\hat{\varphi }\left(\frac{\omega }{2}\right)$$
其中，
$$\mathrm{H}(\omega )=\frac{1}{\sqrt{2}}\sum _{n\in Z}{h}_n{\mathrm{e}}^{-\mathrm{i}n\omega }$$
这个称为低通滤波器。

小波方程

类比尺度方程，对于$\psi (x)\in W_0\subseteq V_1$，存在 { g n ; n ∈ Z } ∈ l 2 \left\{g_{n} ; n \in Z\right\} \in l^2 {gn​;n∈Z}∈l2，使得，
$$\psi (t)=\sqrt{2}\sum _{n\in Z}{g}_n\varphi (2t-n)$$

做傅里叶变换可以得到频域形式，
$$\hat{\psi }(\omega )=G\left(\frac{\omega }{2}\right)\hat{\varphi }\left(\frac{\omega }{2}\right)$$
其中，
$$G(\omega )=\frac{1}{\sqrt{2}}\sum _{n\in Z}{g}_n{\mathrm{e}}^{-\mathrm{i}n\omega }$$
这个称为带通滤波器。$g_n$ 称为脉冲响应系数。

标准正交系的频域表示

引理 设函数 $s(x)\in L^2(R)$, 那么 { s ( t − n ) ; n ∈ Z } \{s(t-n) ; n \in Z\} {s(t−n);n∈Z} 构成 $L^2(R)$ 的 标准正交系，即
$$⟨s(t-n),s(t-l)⟩=\delta (n-l)$$
的充分必要条件是
$$\sum _{k\in Z}|\hat{s}(\omega +2k\pi ){|}^2=1\text{ a.e. }\omega \in R$$

证明：
$$\begin{array}{rl}⟨s(t-n),s(t-l)⟩& =\frac{1}{2\pi }{\int }_R\hat{s}(\omega ){\mathrm{e}}^{-\mathrm{i}n\omega }\left(\overline{\hat{s}(\omega ){\mathrm{e}}^{-\mathrm{i}l\omega }}\right)\mathrm{d}\omega =\\ & \frac{1}{2\pi }{\int }_0^{2\pi }\sum _{k\in Z}|\hat{s}(\omega +2k\pi ){|}^2{\mathrm{e}}^{-\mathrm{i}(n-l)\omega }\mathrm{d}\omega \end{array}$$

低通滤波器和带通滤波器性质

$$|\mathrm{H}(\omega ){|}^2+|\mathrm{H}(\omega +\pi ){|}^2=1,\text{ a.e. }\omega \in R$$

满足这种共轭条件的叫共轭滤波器。

$$|G(\omega ){|}^2+|G(\omega +\pi ){|}^2=1,\text{ a.e. }\omega \in R$$

$$\mathrm{H}(\omega )\bar{G}(\omega )+\mathrm{H}(\omega +\pi )\bar{G}(\omega +\pi )=0,\text{ a.e. }\omega \in R$$

从证明这三条的证明过程中可以看出，这三条的相对应的充要条件分别是：$\varphi (t-n)$ 构成了 $L^2(\mathbb{R})$ 的标准正交系；$\psi (t-n)$ 构成了 $L^2(\mathbb{R})$ 的标准正交系；$\varphi (t-n)$ 张成的空间和 $L^2(\mathbb{R})$ 张成的空间正交。

正交小波的构造

正交小波的充要条件

首先，引入矩阵记号 $\mathbf{M}(\omega )$
$$M(\omega )=\left(\begin{array}{cc}\mathrm{H}(\omega )& G(\omega )\\ \mathrm{H}(\omega +\pi )& G(\omega +\pi )\end{array}\right)$$

构造定理
如果按照前述形式构造小波函数 $\psi (x)$ ，那么，函数族 { ψ ( x − k ) ; k ∈ Z } \{\psi(x-k) ; k \in Z\} {ψ(x−k);k∈Z} 构成 $W_0$ 的标准正交基即 $\psi (x)$ 成为正交小波的充要条件是，矩阵 $\mathbf{M}(\omega )$ 是西矩阵，即
$${\mathbf{M}}^{\ast }(\omega )\mathbf{M}(\omega )=\mathrm{I},\text{ a.e. }\omega \in R$$

正交小波的构造

我们可以选择，
$$G(\omega )={\mathrm{e}}^{-\mathrm{i}\omega }\overline{\mathrm{H}}(\omega +\pi )$$

容易证明，我们这么选 $G$ ，必有 $M$ 是酉矩阵。这是，

$$\hat{\psi }(\omega )={\mathrm{e}}^{-\mathrm{i}\omega /2}\overline{\mathrm{H}}(\pi +\omega /2)\hat{\varphi }(\omega /2)$$

这时候的系数关系是，
$$g_n=(-1{)}^{n-1}{\bar{h}}_{1-n}n\in Z$$

从而，小波函数 $\psi (x)$ 的时域形式为
$$\psi (t)=\sqrt{2}\sum (-1{)}^{n-1}{\bar{h}}_{1-n}\varphi (2t-n)$$

特例：Haar 多分辨分析

定义函数 $\varphi (t)$，
$$\varphi (t)=\{\begin{array}{ll}1& t\in [0,1)\\ 0& t\in \text{others}\end{array}$$

构造，
$$V_m=\text{ Closespan }\left\{2^{\frac{m}{2}}\varphi \left(2^{m}t-n\right);n\in Z\right\}$$

容易证明 ( { V m ; m ∈ Z } ; φ ( t ) ) \left(\left\{V_{m} ; m \in Z\right\} ; \varphi(t)\right) ({Vm​;m∈Z};φ(t)) 是一个 $L^2(T)$ 上的一个正交多分辨分析。

容易计算得到 $h_0=h_1=1/\sqrt{2}$，对应的尺度方程，
$$\phi (t)=\sqrt{2}\left(\frac{1}{\sqrt{2}}\phi (2t)+\frac{1}{\sqrt{2}}\phi (2t-1)\right)$$

这是，可以得到，$g_0=-1/\sqrt{2},g_1=1/\sqrt{2}$，由此得到小波方程，

$$\psi (t)=\sqrt{2}\left(-\frac{1}{\sqrt{2}}\phi (2t)+\frac{1}{\sqrt{2}}\phi (2t-1)\right)$$

这个为，
$$h(t)=\{\begin{array}{cc}-1& 0<\text{ t}<0.5\\ 1& 0.5\le t<1\\ 0& 其他\end{array}$$

和之前的 Haar 小波函数，只差一个符号。