浅谈压缩感知(十五):测量矩阵之spark常数

在压缩感知中,有一些用来评价测量矩阵的指标,如常见的RIP等,除了RIP之外,spark常数也能够用来衡量能否成为合适的测量矩阵。

1、零空间条件NULL Space Condition

在介绍spark之前,先考虑一下测量矩阵的零空间。

这里从矩阵的零空间来考虑测量矩阵需满足的条件:对于K稀疏的信号x,当且仅当测量矩阵的零空间与2K个基向量张成的线性空间没有交集,或者说零空间中的向量不在2K个基向量张成的线性空间中。

上述描述的性质似乎有点难懂,那么与之等价的表述就是spark常数。

2、spark常数定义

简单来说就是,矩阵列线性相关向量组的最小数目

3、spark评价定理

当且仅当spark(Φ)>2k时,可以通过最小0范数优化问题得到k-稀疏信号x的精确近似。

4、线性相关

定义:

定理:

性质:

5、定理证明

反证法。

第一步证明:对于任意的向量y,存在至多一个K稀疏的信号x,使得,则

证明:假设,即根据定义测量矩阵的线性相关列小于或等于2K,从线性相关的定义出发,

存在某个向量,即,使得且h不等于0。

由于,则h可以表示为:,因此

从而得到

但是我们的条件中说明了至多只有一个K稀疏的信号x,因此与原条件矛盾,故假设不成立,原命题成立。

第二步证明:对于任意的向量y,满足的K稀疏信号x至多有一个。

证明:假设K稀疏信号x至少有两个,设为,则

因为,即根据定义测量矩阵的线性相关列大于2K,从线性相关的定义出发,

存在某个向量,使得,且h不等于0。而测量矩阵的零空间应该大于2K维,而假设中的h所在子空间小于或等于2K维,要满足,当且仅当h=0的时候,即,与原假设矛盾,因此假设不成立,原命题成立。

时间: 2024-10-22 08:22:19

浅谈压缩感知(十五):测量矩阵之spark常数的相关文章

浅谈压缩感知(五):单像素相机

前面介绍了关于压缩感知的一些理论知识,这里介绍压缩感知最简单最开始的应用--单像素相机Single Pixel Camera. 1.单像素相机的模型与结构: 如下图所示: PD是光感器件(即单像素),对应公式中的yi: 场景图像对应公式中的f: DMD是数字微镜阵列,用来生成测量矩阵,每一次对应公式中的Φ其中一行. 实际的实验场景: 2.单像素相机的重建结果

浅谈压缩感知(二十四):压缩感知重构算法之子空间追踪(SP)

主要内容: SP的算法流程 SP的MATLAB实现 一维信号的实验与结果 测量数M与重构成功概率关系的实验与结果 SP与CoSaMP的性能比较 一.SP的算法流程 压缩采样匹配追踪(CoSaMP)与子空间追踪(SP)几乎完全一样,因此算法流程也基本一致. SP与CoSaMP主要区别在于"Ineach iteration, in the SP algorithm, only K new candidates are added, while theCoSAMP algorithm adds 2K

浅谈压缩感知(十六):感知矩阵之RIP

在压缩感知中,总是看到"矩阵满足RIP"之类的字眼,没错,这是一个压缩感知绕不开的术语,有限等距性质(Restricted Isometry Property, RIP). 注意:RIP性质针对的同样是感知矩阵而非测量矩阵. 0.相关概念与符号 1.RIP定义 中文版: 英文版: 概括: (RIP)矩阵满足2K阶RIP保证了能够把任意一个K稀疏信号θK映射为唯一的y,也就是说要想通过压缩观测y恢复K稀疏信号θK,必须保证传感矩阵满足2K阶RIP,满足2K阶RIP的矩阵任意2K列线性无关

浅谈压缩感知(七):常见测量矩阵的MATLAB实现

1.随机高斯测量矩阵 function [ Phi ] = GaussMtx( M,N ) %GaussMtx Summary of this function goes here % Generate Bernoulli matrix % M -- RowNumber % N -- ColumnNumber % Phi -- The Gauss matrix %% Generate Gauss matrix Phi = randn(M,N); %Phi = Phi/sqrt(M); end 2

浅谈压缩感知(二十六):压缩感知重构算法之分段弱正交匹配追踪(SWOMP)

主要内容: SWOMP的算法流程 SWOMP的MATLAB实现 一维信号的实验与结果 门限参数a.测量数M与重构成功概率关系的实验与结果 SWOMP与StOMP性能比较 一.SWOMP的算法流程 分段弱正交匹配追踪(Stagewise Weak OMP)可以说是StOMP的一种修改算法,它们的唯一不同是选择原子时的门限设置,这可以降低对测量矩阵的要求.我们称这里的原子选择方式为"弱选择"(Weak Selection),StOMP的门限设置由残差决定,这对测量矩阵(原子选择)提出了要求

浅谈压缩感知(二十八):压缩感知重构算法之广义正交匹配追踪(gOMP)

主要内容: gOMP的算法流程 gOMP的MATLAB实现 一维信号的实验与结果 稀疏度K与重构成功概率关系的实验与结果 一.gOMP的算法流程 广义正交匹配追踪(Generalized OMP, gOMP)算法可以看作为OMP算法的一种推广.OMP每次只选择与残差相关最大的一个,而gOMP则是简单地选择最大的S个.之所以这里表述为"简单地选择"是相比于ROMP之类算法的,不进行任何其它处理,只是选择最大的S个而已. gOMP的算法流程: 二.gOMP的MATLAB实现(CS_gOMP

浅谈压缩感知(三十):压缩感知重构算法之L1最小二乘

主要内容: l1_ls的算法流程 l1_ls的MATLAB实现 一维信号的实验与结果 前言 前面所介绍的算法都是在匹配追踪算法MP基础上延伸的贪心算法,从本节开始,介绍基于凸优化的压缩感知重构算法. 约束的凸优化问题: 去约束的凸优化问题: 在压缩感知中,J函数和H函数的选择: 那么,后面要解决的问题就是如何通过最优化方法来求出x. 一.l1_ls的算法 l1_ls,全称?1-regularized least squares,基于L1正则的最小二乘算法,在标准内点法的基础上,在truncate

浅谈压缩感知(十三):压缩感知与传统压缩

导言: 压缩感知,顾名思义,就是感知压缩,这里包含两层意思,1.感知,即采集或采样,在传统的信号采集中,为了不失真,必须满足Nyquist采样定理,在上一篇博文已经介绍了压缩感知在采样上与传统信号采集的联系与区别,参照http://www.cnblogs.com/AndyJee/p/5050321.html : 2.压缩,即数据压缩,这里是有区别与传统压缩的一个概念,因此本文主要从压缩的角度出发,介绍压缩感知与传统压缩的联系与区别. 传统压缩: 传统的信号获取与处理过程包括:采样.压缩.传输.解

浅谈压缩感知(二):理论基础

主要内容: 信号的稀疏表示 编码测量(采样过程) 恢复算法(非线性) 一.信号与图像的稀疏表示 在DSP(数字信号处理)中,有个很重要的概念:变换域(某个线性空间:一组基函数支撑起来的空间) 一般而言,我们的信号都是在时域或空域中来表示,其实我们可以在其他变换域中通过某些正交基函数的线性组合来表示信号.如:sinusoids, wavelets, curvelets, Gabor functions,. . . 对于某个变换域或空间,其基函数是确定的,只要得到系数α的这一组值,即可通过该系数向量