有限等距性质RIP理解
有限等距性质(Restricted Isometry Property, RIP)
记录当下的一些个人理解,以后学多了再改改。
疑问:
为什么需要RIP准则?
为什么介绍K阶RIP准则后恢复K稀疏信号又要用2K阶RIP的性质?
解惑:
在考虑稀疏信号的可恢复性时,我们先考虑简单的,零空间特性(NSP)和spark(Φ)spark\left( \varPhi \right)spark(Φ)(spark = sparse + rank)适用于没有噪声的情况。一个矩阵的spark定义为这个矩阵列向量中最少线性相关列向量的个数,所以我们有:
定理 当且仅当spark(Φ)>2Kspark\left( \varPhi \right) >2Kspark(Φ)>2K,对于任何矢量y∈RMy\in \mathbb{R}^My∈RM,至多存在一个信号x∈ΣKx\in \varSigma _Kx∈ΣK使得y=Φxy=\varPhi xy=Φx。
从这个定理的反证法中,我们可以看出一些问题,当spark(Φ)⩽2Kspark\left (\varPhi\right)\leqslant2Kspark(Φ)⩽2K时,矩阵Φ\varPhiΦ中至多存在2K个线性相关的列向量,那么一定存在一个向量h=x−x′h=x-x'h=x−x′属于Φ\varPhiΦ的零空间且满足Φh=Φ(x−x′)=0\varPhi h=\varPhi\left(x-x'\right)=0Φh=Φ(x−x′)=0,因而Φx=Φx′\varPhi x=\varPhi x'Φx=Φx′。因此我们可以理解为,要想恢复出K稀疏的信号,矩阵Φ\varPhiΦ中最少线性相关的列向量个数一定要大于2K。这个更好理解所以先放在前面。
问 1:那么当测量值中包含噪声时,能否唯一的恢复出K稀疏信号,就要对矩阵Φ\varPhiΦ施加更严格的限制,因此引入了RIP。注意,别看下面的公式是Φ\varPhiΦ与x相乘,实际上这仅仅是Φ\varPhiΦ的性质,你可以把Φ\varPhiΦ与x相乘看作是从Φ\varPhiΦ中取了K个列出来按一定系数加在了一起。
定义 如果存在δK∈(0,1)\delta_K\in \left(0,1\right)δK∈(0,1),使得
(1−δK)∥x∥22⩽∥Φx∥22⩽(1+δK)∥x∥22\left( 1-\delta _K \right) \lVert x \rVert _{2}^{2}\leqslant \lVert \varPhi x \rVert _{2}^{2}\leqslant \left( 1+\delta _K \right) \lVert x \rVert _{2}^{2} (1−δK)∥x∥22⩽∥Φx∥22⩽(1+δK)∥x∥22对所有x∈ΣK{x:∥x∥2⩽K}x\in\varSigma_K \{x:\lVert x\rVert_2\leqslant K\}x∈ΣK{x:∥x∥2⩽K}都成立,其中如果针对所有K阶稀疏矢量x均满足上式的最小常数δK\delta_KδK,则矩阵Φ\varPhiΦ满足K阶约束等距特性,δK\delta_KδK称为约束等距常数。
我们先来看看这个定义,可以用能量来解释,向量的2范数的平方一般表示信号能量,在上式中可以看出的是,向量x在乘上一个矩阵Φ\varPhiΦ后能量将会在一个范围内波动,这个范围由参数δK\delta_KδK控制。如当δK\delta_KδK为0时,信号x在乘上一个矩阵后能量不变,那么就可以认为矩阵Φ\varPhiΦ是一个正交阵,他对向量x做了正交变换。然而实际上我们不会有这种随便取K列就完全正交的矩阵(况且通常这个K列的子集就不是个方阵),这个参数δK\delta_KδK的范围为(0,1)(0,1)(0,1),当参数δK\delta_KδK很小的时候,可以解释为K个列是“近似正交”的。那么RIP可以看作是一个矩阵中K列组成的子集与正交阵的相似程度。
问 2:前面介绍的是K阶RIP,它就是一个矩阵的性质,一个定义,与我们是否恢复多少稀疏的信号并不直接相关。为了证明压缩感知的方法能恢复出一个K阶信号,我们需要的是一个满足2K阶RIP特性的矩阵。对于这样的矩阵Φ\varPhiΦ来说,我们可以看作从矩阵Φ\varPhiΦ中任意拿出2K个列向量构成的子集应该都是“近似正交的”,那么可以理解为这2K个列都是不相关的,这个矩阵的最少线性相关的列向量个数至少为2K+1,那它也就相当于满足了spark条件,可以唯一的恢复出一个K稀疏的信号。具体的证明过程在很多书和论文中都有,与spark的证明蛮相似的,主要关心概念理解就不抄一遍证明了。后续研究有结论1是当δK<1\delta_K<1δK<1时零范数可以恢复出唯一的K稀疏信号,结论2是当δK<2−1\delta_K<\sqrt{2}-1δK<2−1时可以保证一范数也可以恢复出唯一的K稀疏信号,结论2就衍生出了很多凸优化相关的解法。
总结一下性质:传感矩阵Φ\varPhiΦ只有满足了2K阶RIP才能保证一个K阶稀疏信号对观测信号y的映射是唯一的,这个传感矩阵不能有2K列线性相关。 然后还有对于一个感知矩阵而言,一定阶数的 RIP 常数越小越好,其物理意义在于RIP 常数越小说明该矩阵内的任意小于该阶数的列向量之间的正交性越好,利用其对稀疏信号进行感知采样时,感知矩阵各列“采集”到的信息差异性越大,越有利于进行信号的重构和恢复。对于稀疏测向而言,这个传感矩阵都是给定的,我们对这个矩阵的优化也主要集中在如何是其列间相关性变得更弱,相关性越弱测向结果的分辨力应该就越好,测角结果也越好。
有限等距性质RIP理解相关推荐
- 压缩感知测量矩阵之有限等距性质(Restricted Isometry Property, RIP)
题目:压缩感知测量矩阵之有限等距性质(Restricted Isometry Property,RIP) 阅读压缩感知的文献,RIP绝对是一个抬头不见低头见的英文简写,也就是有限等距性质(Restri ...
- 【转】压缩感知测量矩阵之有限等距性质(Restricted Isometry Property, RIP)
出处:http://blog.csdn.net/jbb0523/article/details/44565647 题目:压缩感知测量矩阵之有限等距性质(Restricted Isometry Prop ...
- 【压缩感知合集7】压缩感知RIP有限等距性:定义解析,理解说明,数学原理推导
0 前情提要 0.1 数学模型和总体框图如下 给定输入信号X∈RN×1\boldsymbol{X} \in \mathbb{R}^{N\times1}X∈RN×1,最终想要得到压缩信号A∈RM×1\b ...
- 压缩感知高斯测量矩阵matlab,浅谈压缩感知(十七):测量矩阵之有限等距常数RIC的计算...
有限等距常数(RestrictedIsometry Constant, RIC)是与有限等距性质(Restricted IsometryProperty, RIP)紧密结合在一起的一个参数. 一.RI ...
- python分析数据的相关性质_理解数据的性质_Python数据分析实战应用_数据挖掘与分析视频-51CTO学院...
为什么学Python: 重要:数据分析是职业技能必备,Python是大数据分析** 趋势:Python是目前非常火的编程语言,使用人多 好学:学习简单,容易上手,使用灵活,可扩展强 **:会Pytho ...
- 压缩感知(compressed sensing)的通俗解释
作者:咚懂咚懂咚 来源: 知乎 链接:https://zhuanlan.zhihu.com/p/22445302 ------------------------------------------- ...
- 形象易懂讲解算法II——压缩感知
作者:咚懂咚懂咚 链接:https://zhuanlan.zhihu.com/p/22445302 来源:知乎 著作权归作者所有.商业转载请联系作者获得授权,非商业转载请注明出处. 之前曾经写过一篇关 ...
- 浅谈压缩感知(十六):感知矩阵之RIP
在压缩感知中,总是看到"矩阵满足RIP"之类的字眼,没错,这是一个压缩感知绕不开的术语,有限等距性质(Restricted Isometry Property, RIP). 注意: ...
- 【压缩感知合集6】压缩感知为什么可以恢复信号;为什么需要满足稀疏性条件、RIP条件、矩阵不相关等限制条件才可以恢复信号的逻辑分析
0 压缩感知的理论依据前言 主要想讲清楚的问题是: 为什么压缩感知在随机采样的情况下可以对信号进行恢复? 其实这个问题也可以换一个方式理解: 在满足什么条件的情况下,信号可以通过压缩感知进行压缩并恢复 ...
最新文章
- Python基础教程学习笔记:第一章 基础知识
- linux:交叉编译后程序的下载方式
- Cache一致性与2种基本写策略(1)
- 原生js实现简单的焦点图效果
- 网络基础 2-1(应用层,HTTP三点注意,HTTP协议格式, 最简单的HTTP服务器)
- 让IE6-IE8 支持HTML5
- 修改hostname不重启机器并生效
- Pyinstaller打包附带DLL、图标和压缩EXE方法
- python给函数设置超时时间已到_在 Linux/Mac 下为Python函数添加超时时间
- android连接打印机打印pdf文件,在android中打印pdf
- PS零基础自学笔记:常见操作方法记录(去水印、抠图、调色调)
- 幼儿园数学目标_大班数学领域目标
- html标题字体文本居中字体颜色,HTML字体标签:如何控制字体大小与颜色
- linux中iso源码解压_linux下解压iso镜像文件方法
- android自定义壁纸制作,Android 自定义View实现画背景和前景(ViewGroup篇)
- Android5.0后外置SDCard读写
- 解决filebeat 报错 Failed to publish events
- 电话薄程序java程序实现
- iperf3.1源码解析笔记(1)-网络压力测试工具iperf的使用
- 学Java可以不学算法吗-学java有必要学算法吗?