Featured image of post 压缩感知

压缩感知

压缩感知是基于信号稀疏性先验,通过压缩测量和优化重建实现信号处理的技术。

利用信号的稀疏性作为先验知识,完成对信号的压缩。

基本公式

  • $x$ 为信号
  • $s$ 为稀疏系数,$x=\Psi s$
  • $y$ 为观测值,$y = \Phi\Psi s = Rs$
  • $\Phi$ 为观测矩阵
  • $\Psi$ 为稀疏矩阵

如何重建

$L_0$最小范数(理论解)

$$ \min_s \Vert s\Vert_0\quad s.t. \quad y=Rs $$

$L_1$最小范数

求解$L_0$最小范数是NP-hard问题,但如果信号足够稀疏($M \gg O(Klog(N/K))$),可以使用$L_1$最小范数得到相同结果。

$$ \min_s \Vert s\Vert_1\quad s.t. \quad y=Rs $$

存在噪声的情况下

没有K的先验知识

$$ \hat s = \arg\min_s{\{\Vert y-Rs\Vert_2^2+\lambda_K\Vert s\Vert_1\}} $$

$\lambda_K$为拉格朗日乘子,取决于采样点数$N$和噪声水平$\sigma_\varepsilon$

有K的先验知识

$$ \hat s= \min_s \Vert s\Vert_1\quad s.t. \quad \Vert y-Rs\Vert_2 \lt \sigma_\varepsilon $$

性质

非相关性

通过最大化每次测量的信息效率分散信号特征到所有采样点中。

$$ \mu(\mathrm{\bf{R}})=\mu(\bf{\Phi},\bf{\Psi})=\sqrt{N}\operatorname*{max}_{1\leq k,j\leq N}\frac{\vert\langle\varphi_{k},\psi_{j}\rangle\vert}{\Vert\varphi_{k}\Vert_2\Vert\psi_{j}\Vert_2} $$

典型测量基$\Phi$ 和稀疏基$\Psi$组合 1. 傅里叶矩阵+单位矩阵:当信号自身具有稀疏性时 2. 噪声基+小波基 3. 随机矩阵+任意固定基

限制等距条件(RIP)

保证测量矩阵对稀疏信号的几何结构保持稳定

$$ \left(1-\delta_{s}\right)\Vert \nu\Vert_{2}^{2}\ \leq\ \Vert\mathbf{R}\,\nu\Vert_{2}^{2}\ \leq\ \left(1+\delta_{s}\right)\Vert\nu\Vert_{2}^{2} $$

典型重建方法

基于松弛理论

凸优化

凸优化算法利用$L_1$范数最小化来求解$L_0$问题等效的凸优化问题问题,问题如下

$$ \hat s = \arg\min_s{\lVert s\rVert_1} \qquad \text{s.t.} \quad \lVert Rs-y\rVert_2\le\varepsilon. $$

该类算法可以分为内点法(Interior Point Method, IPM)和一阶算法。 内点法的典型算法有基追踪(BP)、$L_1$正则化最小二乘($L_1$-LS)、$L_1$-magic等。它们属于重构精度较高的高阶算法,但是它们对高阶问题的求解效率很低。 一阶算法的代表有梯度投影(GPSR)、软阈值迭代(IST, TwIST, FIST)、Bregman、不动点延拓(FPC, FPC-AC)、基于Nesterov理论的算法(NESTA)、一阶锥求解器(TFOCS)、交替方向法(ADM/ADMM)等。这类方法对高维问题的求解非常有效。 $L_1$范数最小化算法可以用来重构<u>稀疏性较差的信号</u>,对加性噪声具有良好的鲁棒性。

非凸优化

非凸优化与之类似,主要的研究方向是基于迭代重加权原理和基于$L_q(0<q<1)$正则化理论。 基于迭代重加权原理的算法由优化最小法(Majorization Minimization)框架推导而来。基本思想是利用与未知信号相关的加权矩阵求解一个优化问题序列。典型的算法有迭代重加权最小二乘(IRLS)和FOCUSS等。 基于$L_q(0<q<1)$正则化理论的算法所要求解的问题如下

$$ \hat s = \arg\min_s \{{\lVert Rs-y\rVert_2}+\lambda\lVert s\rVert^q_q\}. $$

典型的算法有$L_{1/2}$阈值迭代算法和$L_{2/3}$阈值迭代算法。 相比于凸优化算法,非凸优化算法通常对稀疏信号的<u>重构精度更高</u>,然而若参数选择不合适<u>极易收敛到局部极小值</u>,导致信号重构结果出现偏差。

基于贪婪追踪

贪婪追踪算法是一类搜索算法,此类算法主要包含两个步骤,即支集选择和系数更新。现有各种贪婪追踪算法主要由匹配追踪(MP)和正交匹配追踪(OMP)发展而来。典型的算法有OMP、IHT等,此类算法优势是<u>在稀疏性较强时,能迅速重建所需信号。当信号稀疏性较弱或测量存在噪声时,运算效率低、重建结果与真实值之间会存在很大偏差</u>

基于贝叶斯

贝叶斯重构算法是通过将信号稀疏性与先验信息关联,从概率论角度求解稀疏信号的一类算法。典型的有基于稀疏性贝叶斯学习和基于图模型的信息传递算法。

使用 Hugo 构建
主题 StackJimmy 设计