过完备非正交基的典型理论

过完备非正交基的典型理论

引言：从正交到过完备的数学演进

在现代信号处理和数学分析中，过完备非正交基系统的出现标志着从经典正交分解向更灵活表示方法的重要转变。与传统的正交基不同，过完备系统通过引入冗余性获得了更强的表示能力和鲁棒性。这一理论框架的核心在于框架理论（Frame Theory），它为信号的稳定分解与重构提供了严格的数学保证。

框架理论的数学基础

设$\mathcal{H}$为可分希尔伯特空间，序列{ek}k∈N\{e_k\}_{k\in\mathbb{N}}{ek​}k∈N​称为$\mathcal{H}$的框架，当且仅当存在常数$0<A\le B<\infty$，使得对所有$v\in \mathcal{H}$满足框架条件：

$$A\Vert v{\Vert }^2\le \sum _k|⟨v,e_k⟩{|}^2\le B\Vert v{\Vert }^2$$

这个基本不等式确保了信号能量在变换域中的有界性。常数$A$和$B$分别称为下框架界和上框架界，它们刻画了框架的稳定性。当$A=B$时，框架称为紧框架，此时重构过程特别简单。

框架算子$S:\mathcal{H}\to \mathcal{H}$定义为$Sv={\sum }_k⟨v,e_k⟩e_k$，它是正定自伴算子，其谱包含在区间$[A,B]$内。框架算子的可逆性保证了信号的完美重构。对于任意信号$v$，重构公式为：

$$v=\sum _k⟨v,e_k⟩{\tilde{e}}_k=\sum _k⟨v,{\tilde{e}}_k⟩e_k$$

其中{e~k}={S−1ek}\{\tilde{e}_k\} = \{S^{-1}e_k\}{e~k​}={S−1ek​}构成标准对偶框架。这种双向重构公式体现了过完备系统的灵活性——同一信号可以通过不同方式表示。

小波框架

小波框架将经典正交小波理论推广到过完备情形，提供了更灵活的时频分析工具。与正交小波基${\psi }_{j,k}(x)=2^{j/2}\psi (2^{j}x-k)$不同，小波框架系统包含多个生成函数${\psi }^{(m)},m=1,...,M$，形成：

{ψj,k(m)(x)=2j/2ψ(m)(2jx−k):j,k∈Z,m=1,...,M}\{\psi_{j,k}^{(m)}(x) = 2^{j/2}\psi^{(m)}(2^j x - k) : j,k \in \mathbb{Z}, m = 1,...,M\}{ψj,k(m)​(x)=2j/2ψ(m)(2jx−k):j,k∈Z,m=1,...,M}

这种多生成元结构允许设计具有对称性、平滑性等理想性质的小波，而这些性质在正交约束下往往无法同时实现。

多分辨分析的过完备扩展：

传统MRA基于嵌套逼近空间$V_j\subset V_{j+1}$和正交细节空间$W_j\perp V_j$，满足$V_{j+1}=V_j\oplus W_j$。在框架情形，我们引入过完备多分辨分析（OMRA），其中直和关系被松弛为：

$$V_{j+1}=V_j+W_j^{(1)}+\cdots +W_j^{(M)}$$

这里$W_j^{(m)}$是由第$m$个小波生成的细节空间。重要的是，这些空间不必正交，但需满足稳定性条件：

$$A\sum _{m=1}^M\Vert P_{W_j^{(m)}}f{\Vert }^2\le \Vert P_{V_{j+1}}f-P_{V_j}f{\Vert }^2\le B\sum _{m=1}^M\Vert P_{W_j^{(m)}}f{\Vert }^2$$

其中$P_V$表示到空间$V$的正交投影。

完美重构条件的Fourier域刻画：

离散小波框架通过过采样滤波器组实现。设$\varphi$为尺度函数，${\psi }^{(m)}$为小波函数，它们满足细化方程：

$$\varphi (x)=\sqrt{2}\sum _k{h}_0[k]\varphi (2x-k)$$
$${\psi }^{(m)}(x)=\sqrt{2}\sum _k{h}_m[k]\varphi (2x-k)$$

在Fourier域，这些关系表示为：

$$\hat{\varphi }(2\omega )=\frac{1}{\sqrt{2}}{H}_0(\omega )\hat{\varphi }(\omega )$$
$${\hat{\psi }}^{(m)}(2\omega )=\frac{1}{\sqrt{2}}{H}_m(\omega )\hat{\varphi }(\omega )$$

其中$H_m(\omega )={\sum }_k{h}_m[k]e^{-ik\omega }$是滤波器的频率响应。

完美重构的必要充分条件可通过Zak变换表述。定义Zak变换：

$$Z_f(\omega ,t)=\sum _{k\in \mathbb{Z}}f(t+k)e^{-2\pi ik\omega }$$

则小波框架构成紧框架当且仅当：

$$\sum _{m=0}^M|Z_{H_m}(\omega ,t){|}^2=M+1\text{a.e. }(\omega ,t)\in [0,1{)}^2$$

这个条件的深刻之处在于它将时域的框架性质转化为频域的相对简单条件。

非下采样小波变换的算子理论：

非抽取小波变换（Undecimated Wavelet Transform, UWT）定义为线性算子U:L2(R)→ℓ2(Z2×{0,1,...,M})\mathcal{U}: L^2(\mathbb{R}) \to \ell^2(\mathbb{Z}^2 \times \{0,1,...,M\})U:L2(R)→ℓ2(Z2×{0,1,...,M})：

Uf={(⟨f,ϕj,k⟩,⟨f,ψj,k(m)⟩):j,k∈Z,m∈{1,...,M}}\mathcal{U}f = \{(\langle f, \phi_{j,k} \rangle, \langle f, \psi_{j,k}^{(m)} \rangle) : j,k \in \mathbb{Z}, m \in \{1,...,M\}\}Uf={(⟨f,ϕj,k​⟩,⟨f,ψj,k(m)​⟩):j,k∈Z,m∈{1,...,M}}

该变换具有平移协变性：对于平移算子$T_{a}f(x)=f(x-a)$，

$$\mathcal{U}(T_{a}f)=T_a\mathcal{U}f$$

其中右侧的$T_a$是在系数空间上的相应平移。

通过在每个尺度保持原始采样率，UWT实现了框架算子：

$$Sf=\sum _{j,k}⟨f,{\varphi }_{j,k}⟩{\varphi }_{j,k}+\sum _{j,k,m}⟨f,{\psi }_{j,k}^{(m)}⟩{\psi }_{j,k}^{(m)}$$

当滤波器满足紧框架条件时，$S=(M+1)I$，重构简化为：

$$f=\frac{1}{M+1}\sum _{j,k,m}⟨f,{\psi }_{j,k}^{(m)}⟩{\psi }_{j,k}^{(m)}$$

稳态小波变换的数学结构：

稳态小波变换（Stationary Wavelet Transform, SWT）采用"算法à trous"，避免了传统小波变换中的下采样。在第$j$层，滤波器被零插值：

$${\tilde{h}}_j[k]=\{\begin{array}{ll}h[k/2^j]& \text{if }k=m\cdot 2^j\\ 0& \text{otherwise}\end{array}$$

递归关系为：
$$c_j[n]=\sum _k{\tilde{h}}_j[k]c_{j-1}[n-k]=\sum _{k}h[k]c_{j-1}[n-k\cdot 2^j]$$
$$d_j[n]=\sum _k{\tilde{g}}_j[k]c_{j-1}[n-k]=\sum _{k}g[k]c_{j-1}[n-k\cdot 2^j]$$

这种结构保持了每层的采样率，产生冗余度为$J+1$的表示（$J$为分解层数）。

重构通过cycle spinning技术实现：

$$f=\frac{1}{2^J}\sum _{s=0}^{2^J-1}{T}_{-s}\text{ISWT}(T_s\text{SWT}(f))$$

其中$T_s$表示循环移位算子。这种平均过程有效抑制了伪吉布斯现象，同时保持了平移不变性。

Gabor框架

Gabor框架通过对窗函数$g\in L^2(\mathbb{R})$进行时频平移构造，形成系统：

G(g,a,b)={e2πimbxg(x−na):m,n∈Z}\mathcal{G}(g, a, b) = \{e^{2\pi imbx} g(x - na) : m, n \in \mathbb{Z}\}G(g,a,b)={e2πimbxg(x−na):m,n∈Z}

其中$a>0$是时间采样参数，$b>0$是频率采样参数。当该系统满足框架条件时，称为Gabor框架。

密度定理与Nyquist临界：

Gabor框架的一个基本约束是密度定理：为使$\mathcal{G}(g,a,b)$构成$L^2(\mathbb{R})$的框架，必要条件是$ab\le 1$。当$ab=1$时，系统处于临界密度；当$ab<1$时，系统过完备。这种过完备性允许良好的时频局部化。密度定理的证明基于测度论论证。设$\Lambda =a\mathbb{Z}\times b\mathbb{Z}$为时频格点集，其在基本周期$[0,a)\times [0,b)$中的密度为$1/(ab)$。Landau定理表明，对于$L^2(\mathbb{R})$中的任何框架，必须有：

$$\text{lower density}\le 1\le \text{upper density}$$

因此$1/(ab)\le 1$，即$ab\le 1$。

Balian-Low定理与局部化的根本限制：

Balian-Low定理揭示了正交性与局部化之间的根本矛盾：若$\mathcal{G}(g,a,b)$在$ab=1$时构成$L^2(\mathbb{R})$的正交基，则：

$$\left({\int }_{-\infty }^{\infty }{t}^2|g(t){|}^{2}dt\right)\left({\int }_{-\infty }^{\infty }{\xi }^2|\hat{g}(\xi ){|}^{2}d\xi \right)=\infty$$

这意味着不存在同时具有良好时频局部化的正交Gabor基。然而，过完备Gabor框架（$ab<1$）可以绕过这一限制。

定理的证明依赖于量子调和振子的谱理论。定义算子：

$$\hat{H}=\frac{1}{2}(-\frac{d^2}{d{x}^2}+x^2)-\frac{1}{2}$$

其特征函数为Hermite函数$h_n$，特征值为$n$（$n=0,1,2,...$）。

对于Gabor正交基$g_{m,n}(x)=e^{2\pi imx/a}g(x-na)$，假设$\int t^2|g(t){|}^{2}dt<\infty$且$\int {\xi }^2|\hat{g}(\xi ){|}^{2}d\xi <\infty$。利用量子调和振子算子的谱性质，可以证明这与正交性条件$⟨g_{m,n},g_{k,l}⟩={\delta }_{m,k}{\delta }_{n,l}$矛盾。

Zak变换与框架界的精确计算：

Zak变换为Gabor分析提供了强大的数学工具。对于$g\in L^2(\mathbb{R})$，其Zak变换定义为：

$$(Zg)(t,\omega )=\sum _{k\in \mathbb{Z}}g(t+ka)e^{-2\pi ika\omega /b}$$

Zak变换具有准周期性：

$$Zg(t+a,\omega )=e^{2\pi ia\omega /b}Zg(t,\omega )$$
$$Zg(t,\omega +b/a)=Zg(t,\omega )$$

这使得Zak变换本质上定义在环面${\mathbb{T}}^2=[0,a)\times [0,b/a)$上。

通过Zak变换，框架界可精确表示为：

$$A={\text{ess inf}}_{(t,\omega )\in [0,a)\times [0,b/a)}|Zg(t,\omega ){|}^2$$
$$B={\text{ess sup}}_{(t,\omega )\in [0,a)\times [0,b/a)}|Zg(t,\omega ){|}^2$$

特别地，当$ab=1$时，Gabor系统构成紧框架当且仅当$|Zg(t,\omega )|=\text{const}$ a.e.

Wexler-Raz恒等式与对偶窗计算：

Wexler-Raz恒等式建立了原始窗$g$和对偶窗$\gamma$的Zak变换之间的关系：

$$Zg(t,\omega )\overline{Z\gamma (t,\omega )}=ab$$

这个恒等式的证明基于Poisson求和公式。对于Gabor框架的分析算子$C_{g}f=\{⟨f,g_{m,n}⟩{\}}_{m,n}$和合成算子$D_{\gamma }c={\sum }_{m,n}{c}_{m,n}{\gamma }_{m,n}$，完美重构条件$D_{\gamma }{C}_g=I$在Zak域转化为点乘关系。

当$ab<1$时，标准对偶窗通过下式计算：

$$Z\tilde{g}(t,\omega )=\frac{Zg(t,\omega )}{{\sum }_{k,l}|Zg(t+ka,\omega +lb/a){|}^2}$$

其中求和在所有使分母非零的$(k,l)$上进行。

Gabor框架的Janssen表示：

对于$ab=1$的情形，Janssen给出了Gabor系统构成框架的几何条件。设$g$的Zak变换在单位正方形$[0,1{)}^2$上的零点集为$Z=\{(t,\omega ):Zg(t,\omega )=0\}$。

则$\mathcal{G}(g,1,1)$构成框架当且仅当零点集$Z$在测度意义下"稀疏"：

$$\mu (Z\cap \{(t,\omega ):|Zg(t,\omega )|<ϵ\})=O({ϵ}^{\alpha })$$

对某个$\alpha >0$成立。这个深刻结果将框架性质与复分析中的零点分布理论联系起来。

标准对偶框架{g~m,n}\{\tilde{g}_{m,n}\}{g~​m,n​}产生对偶Gabor框架，重构公式为：

$$f=\sum _{m,n\in \mathbb{Z}}⟨f,g_{m,n}⟩{\tilde{g}}_{m,n}=\sum _{m,n\in \mathbb{Z}}⟨f,{\tilde{g}}_{m,n}⟩g_{m,n}$$

其中对偶窗$\tilde{g}$满足双正交关系：

$$\sum _{k,l\in \mathbb{Z}}\tilde{g}(t-ka)g(t-ka)e^{2\pi iklb}=\frac{1}{a}{\delta }_l^0$$

紧框架与特殊构造

紧框架因其简单的重构性质而具有特殊重要性。当框架界满足$A=B$时，框架算子$S=AI$，重构简化为：

$$v=\frac{1}{A}\sum _k⟨v,e_k⟩e_k$$

特别地，当$A=B=1$时，框架称为Parseval框架，此时框架算子为恒等算子，框架本身即为其对偶框架。每个正交规范基都是Parseval框架，但反之不然——Parseval框架可以是过完备的。

有限维空间中，对于${\mathbb{R}}^N$中的$K$个单位范数向量构成的框架，归一化冗余度满足：

$$0<A\le \frac{1}{N}\sum _{k=1}^K|⟨e_k,e_k⟩{|}^2=\frac{K}{N}\le B<\infty$$

冗余比$K/N$直接反映了系统的过完备程度。