压缩采样学习

传统的采样信号方法必须遵守香浓采样定理。CS理论认为可以通过远少于传统方法次数来恢复图像。依赖两个原理:1 稀疏度2非相干性
稀疏度:连续信号的“信息率”比由其带宽表现的小很多,离散信号由其自由度决定,这比它本身的长度小得多。更精确的说,CS利用了这样的事实:大多数自然信号是稀疏且可压缩的,他们在合适的基下有更紧凑的表示。
非相干性扩展了时间和频率的二元性,表达了如下的观点:有稀疏表示的物体必须在他们获得的域被铺开延展,就像时域Dirac函数在频域里被铺开延展。换句话说,非相干性认为不同于感兴趣的信号,采样的波形在 Ψ \Psi Ψ下有着稠密的表示。
关键的观测是一个人可以设计有效感知或者采样工具来捕获在一个稀疏信号的有效信息并压缩成小的数据量。这些工具是非改进的并且只需要关联信号和固定的波形(和稀疏基非相干)。在不理解信号下这些采样工具让传感器有效捕获稀疏信号的信息。并且,有方法使用数值优化来重构完整信号

感知问题

y k = ⟨ f , φ k ⟩ , k = 1 , … , m y_{k}=\left\langle f, \varphi_{k}\right\rangle, \qquad k=1, \ldots, m yk=f,φk,k=1,,m关联我们物体和波形 φ k ( t ) \varphi_{k}(t) φk(t)(可以是正弦的波形,y即为傅立叶系数)
我们感兴趣的问题为欠采样问题: A f ~ = y A \tilde{f}=y Af~=y对于n维度的f的信号进行远少于n的m次测量。场景诸如:传感器数目有限,测量昂贵,只能测量几次如(MRI)

稀疏度:大多自然信号在某种基下可以稀疏表示。此时可以丢弃小的稀疏而感觉上没什么损失。
f ( t ) = ∑ i = 1 n x i ψ i ( t ) f(t)=\sum_{i=1}^{n} x_{i} \psi_{i}(t) f(t)=i=1nxiψi(t)考虑 x S x_S xS为保留含有S非0元素的信号。 f S : = Ψ x S f_{S} :=\Psi x_{S} fS:=ΨxS由于 Ψ \Psi Ψ是正交基,存在 ∥ f − f S ∥ ℓ 2 = ∥ x − x S ∥ ℓ 2 \left\|f-f_{S}\right\|_{\ell_{2}}=\left\|x-x_{S}\right\|_{\ell_{2}} ffS2=xxS2
若x是稀疏的那么 x i x_i xi强度排序后下降迅速,那么 ∥ f − f S ∥ ℓ 2 \left\|f-f_{S}\right\|_{\ell_{2}} ffS2的误差很小,也就是说我们可以扔掉稀疏而没有多大损失。这也是很多JPEG-2000编码器的原理。
非相干采样
假设我们有一对 ( Φ , Ψ ) (\Phi, \Psi) (Φ,Ψ)正交基。第一组基 Φ \Phi Φ用来感知物体f,第二个用来表示f。两个正交基的相干性可以表达为: μ ( Φ , Ψ ) = n ⋅ max ⁡ 1 ≤ k , j ≤ n ∣ ⟨ φ k , ψ j ⟩ ∣ \mu(\Phi, \Psi)=\sqrt{n} \cdot \max _{1 \leq k, j \leq n}\left|\left\langle\varphi_{k}, \psi_{j}\right\rangle\right| μ(Φ,Ψ)=n max1k,jnφk,ψj相干性测量了两组基下最大的关联的关联性。
压缩采样主要用低的关联对。例如: Φ \Phi Φ φ k ( t ) = δ ( t − k ) \varphi_{k}(t)=\delta(t-k) φk(t)=δ(tk),而 Ψ \Psi Ψ是傅立叶基, ψ j ( t ) = n − 1 / 2 e i 2 π j t / n \psi_{j}(t)=n^{-1 / 2} e^{i 2 \pi j t / n} ψj(t)=n1/2ei2πjt/n。由于 Φ \Phi Φ是感知矩阵,对应经典采样方案,满足 μ ( Φ , Ψ ) = 1 \mu(\Phi, \Psi)=1 μ(Φ,Ψ)=1
第二个例子 Φ \Phi Φ是噪声基,而 Ψ \Psi Ψ是小波基。和小波基的相干性为 2 \sqrt2 2 D4和D8小波基为2.2,2.9
最后,随机矩阵和任何 Ψ \Psi Ψ是非相干的。
欠采样和稀疏信号恢复
y k = ⟨ f , φ k ⟩ , k ∈ M y_{k}=\left\langle f, \varphi_{k}\right\rangle, \quad k \in M yk=f,φk,kM
通过下面式子
min ⁡ x ~ ∈ R n ∥ x ~ ∥ ℓ 1  subject to  y k = ⟨ φ k , Ψ x ~ ⟩ , ∀ k ∈ M \min _{\tilde{x} \in \mathbb{R}^{n}}\|\tilde{x}\|_{\ell_{1}} \quad \text { subject to } \quad y_{k}=\left\langle\varphi_{k}, \Psi \tilde{x}\right\rangle, \quad \forall k \in M minx~Rnx~1 subject to yk=φk,Ψx~,kM
在满足 f ~ = Ψ x ~ \tilde{f}=\Psi \tilde{x} f~=Ψx~下找到最 l 1 l1 l1范数的x
理论
1相干度越小,需要越少的采样,
2
小波域的稀疏性可以用来信号获取后的压缩,广泛使用在jpeg使用图像 8 × 8 8\times8 8×8的小块的离散余弦变换。视频信号mpeg使用一维的DCT和复杂的视频编码格式如H.264
信号获取后使用有损的编码的缺点在于大量采集的信号并未实际参与最后信号的表达。在很多设备中,测量的分辨率(采样时间和测量的空间分辨率)这会严重限制系统的带宽。测量时间也会被限制。
因此,一个更好地获取稀疏表达下信号的方法是直接记录稀疏域的信号。使用一个线性变换 Φ \Phi Φ

你可能感兴趣的:(压缩采样学习)