基于双极性混沌序列的托普利兹块状感知矩阵*
2021-02-06干红平张涛花燚舒君何立军
干红平 张涛 花燚 舒君 何立军
1) (西北工业大学软件学院,西安 710072)
2) (清华大学电子工程系,北京 100084)
3) (西北工业大学航空学院,西安 710072)
4) (四川大学电子信息学院,成都 610065)
感知矩阵的构造是压缩感知从理论走向工程应用的关键技术之一.由于托普利兹感知矩阵能够支持快速算法且与离散卷积运算相对应,因此具有重要的研究意义.然而常用的随机托普利兹感知矩阵因其元素的不确定性,使得它在实际应用中受到了诸多约束,例如内存消耗较高和不易于硬件加载.基于此,本文结合双极性混沌序列的内在确定性和托普利兹矩阵的优点,提出了基于双极性混沌序列的托普利兹块状感知矩阵.具体地,首先介绍了双极性混沌序列的产生并分析了它的统计特性.其次,构造了双极性托普利兹块状混沌感知矩阵,从相关性方面证明了新建的感知矩阵具有近乎最优的理论保证,并同时证实了它满足约束等距条件.最后,研究了该感知矩阵针对一维信号和图像的压缩测量效果,并与典型感知矩阵进行了对比.结果表明,提出的感知矩阵对这些测试信号具有更好的测量效果,而且它在内存开销、计算复杂度和硬件实现等方面均具有明显的优势.特别地,该感知矩阵非常适用于多输入-单输出线性时不变系统的压缩感知测量问题.
1 引 言
近年来,Candès 等[1]提出的压缩感知(compressed sensing,CS)理论为数据获取与重建提供了一种全新的策略.该理论创造性地将待采样信号的频域带宽特性扩展到更为广泛的变换域稀疏特性,以此利用满足一定特性的感知矩阵对信号进行不相关测量,从而突破了香农-奈奎斯特采样理论的极限.换句话说,压缩感知打破了从信号到信息域的壁垒,实现了对信息直接感知和处理.因此,它一经提出便成为当今世界最为热门的课题之一,并被广泛地应用到众多研究领域,如压缩成像[2]、计算机科学[3]和数据安全[4]等.
为准确描述CS 的数学模型,令x={xi}ni ∈Rn(i=1,2,··· ,n)表示待测量的目标信号,并假定该信号满足以下两个条件之一: 1)x是稀疏度为s(s 2)x可以等价表示为Ψη且η是稀疏度为s的稀疏信号,也就是说,x=Ψη,其中Ψ代表某个变换域或稀疏域.那么可以通过特殊的信息感知算子A:Rn →Rm实现对该信号降维信息感知与重构.数学上,CS 的线性测量过程可表示为 式中y ∈Rm表示测量样本数据,A∈Rm×n又被称为感知矩阵.记ϖ=m/n(m ≪n) 为采样率.由于存在数据降维过程 Rn →Rm,因此从y中恢复信号x似乎是不可能的.但是,只要A能够保证在低维空间 Rm中有且只有一个y与信号x建立起y=Ax的关系,那么便可借助一定的恢复算法将x精确重构出来. 直观上,重建信号xℓ0最简单的方法就是优化范数问题: 但遗憾的是,ℓ0(·) 是一个极其特殊的函数,导致直接求解(1)式是一个NP-hard 问题.为解决该问题,CS 研究者提出了不同的近似策略来平滑ℓ0函数,如高斯-牛顿平滑ℓ0算法[5].除此之外,贪婪算法也常被用来稀疏逼近ℓ0问题的最优解从而得到重构信号xˆ ,如正交匹配追踪(orthogonal matching pursuit,OMP)[6].值得注意的是,当A能够做到降维信息保真时,那么可以将ℓ0问题转化为ℓ1范数稀疏求解问题, 显然,上述优化问题的优化目标ℓ1范数是一个凸函数.因此,可通过已有的优化算法对其进行高效求解,例如基追踪(basis pursuit,BP)[7]. 从上述分析可知,亟需解决以下重要问题: 如何设计A才能做到降维信息保真,即A满足何种性质时才能保证信号x从高维空间 Rn降维到Rm时使y=Ax建立起唯一的对应关系.Candès[8]证明了当感知矩阵满足约束等距性(restricted isometry property,RIP)时,A就能够保证所有稀疏度为s的信号x∈Σs降维后的唯一性[8],其中,Σs表示稀疏度为s的信号的集合. 定义1[8]若矩阵A能够使 成立,则称A(s,δs)满足 -RIP. 满足(3)式的最小常数δs被称为A的s阶约束等距常数(restricted isometry constant,RIC).当RIC 越小时,那么矩阵的约束等距性越好.实际上,(s,δs) -RIP 能够确保A的所有大小为m×s的子矩阵均类似于正交矩阵,从而使得x的“长度”得以保持,以此实现信息的保真,即∥Ax∥22=∥x∥22.RIP 为信号的降维信息保真提供了优雅的几何解释.然而,验证A是否满足RIP 是一个NP-hard 问题,即必须检验A中的所有子矩阵都满足(3)式.因此,有必要去寻求更易于计算且能够保证唯一性的其他标准.其中,矩阵的相关性(coherence)就是这样的一个标准. 定义2[9]给定一个矩阵A,它的自相关系数µ(A)等于A中两两不相等列的最大内积: 其中ai,ak分别表示矩阵A 的第i,k列. 引理1假定矩阵A具有单位ℓ2范数的列向量,其自相关系数µ=µ(A) ,那么A满足 (s,δs) -RIP 且其约束等距常数为δs≤(s −1)µ. 在CS 理论诞生之初,由于随机矩阵(例: 高斯或伯努利随机矩阵)易满足RIP 且具有普适性,因此常被用作感知矩阵.但随机感知矩阵不仅需要巨量的内存空间、高复杂度的运算,而且也难以加载到硬件电路中,这潜在地束缚了它们的应用.为克服这些缺点,CS 研究人员朝着下述三大方向做了尝试.第一,通过优化策略降低感知矩阵对硬件电路实现的要求,其中最常用的方式是将矩阵元素二元化,例如,基于二分图(周长大于4)的双极性感知矩阵[10].但应当指出的是,当信号的维数n增长时,这类矩阵的相关性可能会随之增加,进而影响到它们的感知效率.第二,设计存储要求低且对硬件友好的确定性感知矩阵,例如,Ansari 正则单位感知矩阵[11].应当注意的是,尽管确定性感知矩阵能够在特定的场合良好工作,但它们通常不具有对感知场景的普适性,这大大制约了它们的实际应用.第三,构建与实际系统的测量机制相对应的感知矩阵.这当中最吸引人的是托普利兹(Toeplitz)和循环(Circulant)感知矩阵,例如托普利兹随机感知矩阵[12].该类矩阵被广泛应用于解决线性时不变(linear time invariant,LTI)系统的辨识问题,如稀疏信道估计问题.但令人遗憾的是,托普利兹随机感知矩阵依旧没有摆脱随机元素,这会限制它们的推广和拓展. 显然,为推动CS 从理论走向更多实际应用场景,有必要融合以上3 个方向的优势,并规避它们的弊端.一方面,从数据获取的角度来说,人们希望构建出低内存消耗、低计算复杂度,具有普适性且测量效率优良的感知矩阵.另一方面,从硬件实现的角度来看,构造的感知矩阵应当对应于可行的确定性硬件架构,如LTI 系统.基于这两方面的考虑,可探索一些新的机制来生成确定性感知矩阵,以便减轻在数据感知和硬件实现等方面的负担.实际上,用混沌系统生成的混沌序列来构建感知矩阵可以很好地解决上述两个维度的问题.作为最早的启发性工作,Yu 等[13]提出利用Logistic 混沌序列来构造感知矩阵,并通过一个简单的组合问题证明了该混沌算子满足RIP,这也从根本上保证了它的采样效率.随后,Gan 等[14]设计了Chebysehv 混沌感知矩阵,并通过Johnson-Lindenstrauss 引理证明了该混沌感知算子具有与随机矩阵类似的感知性能.特别地,郭静波等结合混沌序列和循环矩阵的优点,开发了Cat 混沌循环感知矩阵[15],并将其用于二进制信号的信息感知与重构[16].构造这些感知矩阵的一个关键操作是利用抽样操作将实值混沌序列变成独立同分布序列,以保证这些实值混沌序列的独立性,进而使其构造的混沌感知矩阵具有良好的感知效率.然而,抽样操作需要较为复杂的运算,这在一定程度上限制了它们的实际应用. 本文将结合双极性混沌序列(元素仅有–1 和1)的内在确定性和托普利兹矩阵的优点,提出基于双极性混沌序列的托普利兹块状感知矩阵.它不仅直接保留了托普利兹块状感知矩阵的优势,而且可以继承双极性混沌序列的天然优点.特别地,本文提出的构造方法只需要生成双极性混沌序列,毋须对实值混沌序列进行抽样操作.另外,使用托普利兹块状感知矩阵测量数据实际上是对LTI 系统的多输入信号做离散卷积,能够解决众多与卷积运算相关的理论与应用问题,而且可以支持快速算法.因此,新建的托普利兹块状感知矩阵特别适用于多输入-单输出LTI 系统的压缩感知测量问题,例如,多输入-单输出有限脉冲响应系统的参数与时滞估计问题[17].换句话说,它具有易于数字电路实现的潜在优势.本文的主要贡献有: 1) 提出了一种基于对称阈值的双极性混沌序列的生成方法,并分析了该双极性混沌序列的统计性质; 2) 利用双极性混沌序列构建了双极性托普利兹块状感知矩阵,并通过相关性和约束等距条件证明了构造的感知矩阵具有优良的感知效率; 3) 提出的构建框架可推广至不同混沌系统(如Logistic 或Cat 混沌系统),而且也可以派生出Hankel 块、循环块以及堆积块矩阵,进而推动将压缩感知应用到更多LTI 系统的压缩测量问题中. 此外,通过对一维稀疏信号和图像进行数值仿真测试,验证了新建的托普利兹块状混沌感知矩阵与经典的感知矩阵相比在采样效率和恢复效果等方面具有更好的性能. 本文的内容安排如下: 第2 节给出双极性混沌序列的生成方法,并分析其统计特性; 随后,在第3 节中介绍双极性托普利兹块状混沌感知矩阵的构造和性能分析; 第4 节通过实验仿真与其他传统感知矩阵进行比较; 在第5 节得出相关结论. 回顾一个简单α阶一维映射函数,其定义为 式 中zj=τj(z0)∈I=[−1,1] ,z0表示初始值.当1<α ∈N+且−1 ≤z0≤1 时,(5)式就是著名的切比雪夫(Chebyshev)混沌系统,其连续不变测度为.通过对(5)式反复迭代,就可以产生一组切比雪夫实值混沌序列即.随后,为该混沌实值序列引入切比雪夫对称混沌阈值函数,即 其中T(z) 的互补函数为 结合(5)式和(6)式,便可得到一组双极性切比雪夫混沌序列{T(τj(z0))}∞j=0,也就是{T(zj)}∞j=0.在实际应用中,通过一个比较器和一些简单的移位寄存器就可以很容易地硬件生成双极性混沌序列[18].值得注意的是,双极性混沌序列的统计均值可通过下式求得: 代入双极性切比雪夫混沌系统的连续不变测度可得〈T〉τ=1/2 . 引 理2[19]Perron-Frobenius 算 子Pτ作 用 于关于τ(z) 的有界变差函数H(z) 时满足 其中gi(z)=τi−1(z) . 值得注意的是,连续不变测度ρ∗(z)dz满足等式Pτρ∗(z)=ρ∗(z) ,该式又被称为Perron-Frobenius 等式.结合(6)式的混沌阈值函数T(z) 和引理2,可以很容易得到 其 中,li≥0(1 ≤i≤k −1) .因 此,对 于 一 组,它的高阶相关系数可以定义为 引理3对于一组切比雪夫混沌阈值函数,当lj−1≥1 时,可以得到 证明1根据引理2,将Pτ算子作用于每一个切比雪夫混沌阈值函数Tj(·) 的高阶相关系数,并结合(7)式得到重复以上步骤,得证引理3. 设Wk=W0W1···Wk−1代表一个由k位双极性变量组成的任意数字串,其中Wj ∈{1,−1}(j=0,1,··· ,k −1) .显然,Wk有2k种可能性.令表示第h个数字串,其中.对于切比雪夫混沌阈值函数,可以引进一个二元变量 定理1对于一组切比雪夫混沌阈值函数,其生成的双极性混沌序列满足 证明2根据引理3 可知 式中β代表中元素1 的总个数.将〈T〉τ=1/2 代入(9)式中,可以得到 A(i)(i=2,··· ,b) 的定义与(10)式类似,即A(i)的元素为双极性混沌序列确定.这里将A(i)的元素简记为.显然,只需要长度为b(m+d −1) 的双极性切比雪夫混沌序列便可完全确定Bi-TpCM.(10)式中标量用于归一化A的列,以保证降维过程 Rn →Rm中原信号x与测量样本数据y的能量保持一致. 值得注意的是,通过形如(10)式的构造方式不仅可以构造出托普利兹块状矩阵,还能派生出Hankel 块、循环块以及堆积块矩阵,进而推动将压缩感知应用到更多LTI 系统的压缩测量问题中.由于分析这些矩阵性能的方法与托普利兹块状矩阵类似且后者对应于多输入-单输出LTI 系统的压缩感知测量问题,故本文仅分析托普利兹块状感知矩阵. 定理2当b≥2 时,大小为m×n(n=bd) 的双极性托普利兹块状混沌感知矩阵A的自相关系数µ(A) 将以超过 ( 1−4/n) 的概率满足 证明3由定义2 可知,Bi-TpCM 的自相关系数是由其托普利兹子块矩阵决定的,也就是说 式中µ(A(i)) 表示A(i)的自相关系数,∥·∥max代表矩阵的最大范数.根据文献[12]的定理4,可以推出双极性托普利兹子块矩阵A(i)满足 式中ξ ∈(0,1).因此可以推出 结 合Qi,k有d2个 不 同 的 元 素,且存在个组合方式,故有 根据(13)式和(15)式且b≥2 ,可知 当b=1 时,Bi-TpCM 将完全退化为一个双极性托普利兹感知矩阵,文献[12]已分析了这类感知矩阵的性能,这里不再赘述.从定理2 可知,Bi-TpCM 的自相关系数µ(A) 的上限为这与感知矩阵自相关系数的Welch 界只差一个标量 l ogn.定理2 也证实了Bi-TpCM 在自相关性方面接近于最优理论保证.结合引理1 和定理2,可以推出Bi-TpCM 满足约束等距条件,其约束等距常数δs不超过.除此之外,由于双极性托普利兹块状混沌感知矩阵本质上是一个Rademacher 矩阵,因此它必然满足度量集中(concentration of measure)不等式.结合Richard 等[21]关于RIP 的工作,可以有以下推论. 推论1对于一个大小为m×n的双极性托普利兹块状混沌感知矩阵A,当m≥c0slog(n/s) 时,它将以超过 1−2 exp−c1m的概率满足RIP,其中c0和c1为两个常数. 推论1 证实了提出的Bi-TpCM 在约束等距条件下均具有良好的理论保证.为节约篇幅,这里省略了推论1 的详细证明过程.读者可以根据文献[21]的工作,推出相应的证明步骤.在此基础上,可以将以上结论推广至所有混沌系统. 推论2任意能够产生混沌Rademacher 序列的混沌系统均能够构造出Bi-TpCM,并且构建的Bi-TpCM 的自相关系数为,且它能够以接近于 1 的概率满足约束等距条件. 根据以上分析可知,本文提出的构造框架不仅适用于切比雪夫混沌系统,还可以推广至不同的混沌系统,如Logistic 或Cat 混沌系统,并且它们对应的Bi-TpCM 的采样效率能够接近于最优.实际上,将Geršhgorin 圆盘定理应用于Bi-TpCM 的格拉姆(Gram)矩阵,可以得出以下推论. 推论3对于一个大小为m×n的Bi-TpCM,令G=ATA,其 中AT表 示A的 转 置,即G为A的格拉姆矩阵.若存在两个正常数δd和δo,δd+δo=δs ∈(0,1) ,使 得G的 每 一 个 对 角 元 素Gi,i满足|Gi,i −1|<δd,且非负对角元素满足Gi,k(i ′=k)满足|Gi,k|<δo/s,那么该矩阵满足 (s,δs) -RIP. 推论3 表明了可以通过数值实验来验证构造的Bi-TpCM 是否满足约束等距条件.在仿真分析的小节中,会证实Bi-TpCM 的确具有RIP. 为了更好地评估Bi-TpCM 性能,表1 将Bi-TpCM 和其他常用感知矩阵做了一个全面比较.传统的常用感知矩阵包括Den-RgM,Den-Bol,Den-CbM,Top-Rad 和Cir-CaM.其中,Den-RgM和Den-Bol 分别代表高斯和伯努利感知矩阵;Den-CbM 是Gan 等[14]利用Chebyshev 实值混沌序列直接构造出的混沌感知矩阵; Top-Rad 代表随机双极性序列构建的双极性托普利兹随机感知矩阵[12]; Cir-CaM 是Guo 等[15]开发的Cat 混沌循环感知矩阵. 由表1 可知,Bi-TpCM 不仅具有随机感知矩阵的信息感知能力和普适性,而且也继承了托普利兹矩阵支持快速计算和对硬件友好的优势,同时还引入了双极性混沌序列易于存储、计算复杂度低等天然优点.特别地,相比于已存在的混沌感知矩阵,本文提出的构造方法只需要利用生成的双极性混沌序列直接构造Bi-TpCM,毋须对实值混沌序列进行抽样操作.除此之外,构造的托普利兹块状感知矩阵特别适用于多输入-单输出LTI 系统的压缩感知测量问题. 表1 不同感知矩阵的性能比较Table 1.Performance comparisons of different sensing matrices. 本节将通过仿真实验来验证Bi-TpCM 的性能.首先采用数值实验来观察Bi-TpCM 的约束等距现象.随后,通过对一维稀疏信号和图像分别进行压缩采样,来验证Bi-TpCM 的采样性能.为了对比分析,传统的感知矩阵也被加入了同样的稀疏感知场合,涉及的矩阵包括Den-RgM,Den-Bol,Den-CbM,Top-Rad 和Cir-CaM.其中,Den-CbM是由α=6 ,初始值z0=0.17 和抽样步长为10 的实值Chebyshev 混沌序列构造而成.Cir-CaM 则是由初始值z0=0.09 和抽样步长为200 的实值Cat 混沌序列派生而成.Bi-TpCM 是由本文介绍的双极性Chebyshev 混沌序列构造而成.本节的所有实验均在台式电脑,英特尔至强银4110 CPU@2×2.10 GHz (双处理器),64 GB 内存下的Matlab-R2018b 的环境下进行. 首先,根据3.1 节构造出一个Bi-TpCM∈R100×256,构 造 参 数 为α=6 ,z0=0.17 ,b=4 及d=64 .忽略标量,图1(a)给出了构造的Bi-TpCM 的元素分布.可以看出,双极性托普利兹块状混沌感知矩阵的元素–1 与+1 的数量基本相同,这也从侧面验证了用于构造该感知矩阵的双极性混沌序列是Rademacher 序列. 图1(b)和图1(c)分别给出了Bi-TpCM 对应的格拉姆矩阵.可以看出,Bi-TpCM 的Gram 矩阵的对角元素Gi,i均收敛于1,而非对角元素均在0附近浮动.根据推论3 可知,Bi-TpCM 的确具有约束等距性,正如标准的随机感知矩阵一样. 先生成一个稀疏度为s的一维信号x∈R256作为待采样信号,其非零元素服从高斯分布N(0,1/2) .为恢复信号和评估重建质量,欠定重构算法将使用文献[7]提出的基追踪算法; 重构质量标准则使用信噪比(signal-to-noise ratio,SNR),其定义为 图1 忽 略 标 量 的Bi-TpCM及其对应的Gram 矩阵展示图 (a) Bi-TpCM; (b) Gram 矩阵的三维渲染图; (c) Gram 矩阵的等高图Fig.1.Bi-TpCM without the factor 1 / and its Gram matrix: (a) Bi-TpCM; (b) three dimensional rendering of its Gram matrix; (c) contour map of its Gram matrix. 式中表示重构信号.特别地,对于一个时域稀疏信号x,若重构时 S NR(x)≥100 dB ,那么称该信号被完美重构(perfect recovery). 在这些参数设置下,图2 给出了利用Bi-TpCM∈R100×256分别对稀疏度为20 和30 的一维稀疏信号x∈R256进行压缩测量的具体表现.具体地,图2 第一列对应的恢复误差(定义为)和重构SNR 分别是3.58×10–13和124.45 dB,第二列的恢复误差和重构SNR 分别是1.26×10–12和118.99 dB.由图2 可知,一维稀疏信号x能够以极小的误差(恢复误差的数量级是10–12)被重构出来,这也说明了本文新建的Bi-TpCM 具有优良的测量效率. 图2 Bi-TpCM 压缩测量一维信号的重构 (a) s = 20,条形图; (b) s = 20,细节图; (c) s = 30,条形图; (d) s = 30,细节图Fig.2.Reconstructions of one-dimensional signal using Bi-TpCM: (a) s = 20,stem rendering; (b) s = 20,detailed drawing; (c) s =30,stem rendering; (d) s = 30,detailed drawing. 接下来,利用Den-RgM,Den-Bol,Den-CbM,Top-Rad,Cir-CaM 和Bi-TpCM (均∈R100×256)分别去采样稀疏度s变化的一维稀疏信号x∈R256.随后,使用BP 算法去重构原信号以得到,并且针对每一个稀疏度s的x重复实验100 次.最后,通过求平均的方法便可获得这些感知矩阵的平均恢复结果.图3(a)—(c)分别比较了这些感知矩阵的恢复误差、重构SNR 和完美重构概率.由图3可知,采用的典型感知矩阵(包括Den-RgM,Den-Bol,Den-CbM,Top-Rad 和Cir-CaM)整体呈现出了近乎一致的稀疏感知能力.但在细节上,实值混沌感知矩阵(Den-CbM)要略优于随机感知矩阵(Den-RgM 和Den-Bol,值得强调的是,人们普遍认为Den-RgM 具有良好的稀疏感知效率),而它们都优于结构性感知矩阵(Top-Rad 和Cir-CaM). 特别地,图3 清晰地展示出了Bi-TpCM 的性能要明显优于这些传统的常用感知矩阵,这也从侧面说明了Bi-TpCM 比这些感知矩阵更具有接近于理论最优的采样效率.Bi-TpCM 强大的稀疏感知能力为它在CS 中的实际应用提供了根本性保障. 本部分将比较不同感知矩阵对图像信号进行压缩采样时的性能,使用的感知矩阵包括Den-RgM,Den-Bol,Den-CbM,Top-Rad,Cir-CaM和Bi-TpCM.这里选取了两张自然图像作为测试信号x∈R256×256,即“Lena”和“Lin”.这两张图像分别展示在图4(a)和图4(e)中,它们在离散小波变换域中均是可压缩的.本文将采用经典的OMP算法作为图像重建算法,并使用峰值信噪比(peaksignal-to-noise ratio,PSNR)作为重建图像的质量评价指标,其定义为 式中代表重构后的图像.应当指出的是,这里在图像的压缩测量过程中引入了高斯噪声来模拟实际应用场景中的加性噪声和乘性噪声. 图3 分别使用不同的感知矩阵对稀疏度变化的x 进行压缩测量时的重建性能比较 (a) 重建误差; (b) 信噪比 (dB); (c) 完美重建的概率Fig.3.Performance comparisons for recovering x with different sparsity using various sensing matrices,respectively: (a) Recovery error; (b) SNR (dB); (c) perfect recovery probability. 图4 原始图像和Bi-TpCM 在不同采样率下的恢复图像,其中第一行是(a) 原 始“Lena”,(b) ϖ =0.3 ,(c) ϖ =0.6 ,(d) ϖ =0.8 ; 第二行是(e) 原始“Lin”,(f) ϖ =0.3 ,(g) ϖ =0.6 ,ϖ =0.8Fig.4.Original and reconstructed images using Bi-TpCM at different sampling rates.The first row: (a) Original “ Lena” ;(b) ϖ =0.3 ; (c) ϖ =0.6 ; (d) ϖ =0.8 .The second row: (e) Original “Lin”; (f) ϖ =0.3 ; (g) ϖ =0.6 ; (h) ϖ =0.8 . 首先,图像“Lena”和“Lin”分别被Den-RgM,Den-Bol,Den-CbM,Top-Rad,Cir-CaM 和 Bi-TpCM 在不同的采样率下压缩测量和重建,即采样率随着测量样本数据的维度增加而增加.在重构算法OMP 的作用下,可以获得这些感知矩阵对应的重构图像和重建PSNR.值得注意的是,为增加随机感知矩阵(即Den-RgM,Den-Bol 和Top-Rad)性能的稳定性,这里让它们重复执行该实验100 次,后取其平均结果.由于Den-CbM,Cir-CaM 和Bi-TpCM 是确定性感知矩阵,因此它们的实验结果是稳定的,不需要重复测试.为了更好地视觉展示,图4 给出了使用提出的Bi-TpCM 在采样率分别为 0 .3 ,0 .6 和 0 .8 时的恢复图像,其中图像“Lena”对应的重建PSNR 分别是22.52,29.19和32.36 dB; “Lin”的恢复PSNR 分别是23.69,30.22 和33.59 dB.由图4 可知,随着采样率(测量样本数据)的增加,重建的图像也越来越清晰.并且当≥0.5 时,便可得到一幅较为清晰的恢复图像.为节约篇幅,这里省略了其他感知矩阵的恢复图像. 图5 比较了在不同采样率下,Bi-TpCM 和传统感知矩阵(包括Den-RgM,Den-Bol,Den-CbM,Top-Rad 和Cir-CaM)重建图像“Lena”和“Lin”时的PSNR 值.从图5 可知,在低采样率(≤ 0 .35 )下,提出的Bi-TpCM 的性能要明显优于其他感知矩阵,而在较高的采样率(ϖ≥0.4 )时,这些感知矩阵的采样性能相差不大.但整体上,新建的Bi-TpCM 要比传统的感知矩阵做得更好. 图5 在不同采样率下利用不同的感知矩阵对图像进行压缩测量时的重建PSNR 比较 (a) “Lena”; (b) “Lin”Fig.5.Reconstructed PSNR comparisons for image compressed sensing using different sensing matrices at various sampling rates,respectively: (a) “Lena”; (b) “Lin”. 以上对一维信号和图像的测试表明,本文的Bi-TpCM 与经典感知矩阵相比在采样效率和恢复效果等方面具有明显的优势.结合Bi-TpCM 的确定性托普利兹块状结构,同时具有双极性混沌序列的天然优点,因此,Bi-TpCM 具有更好的应用潜力. 本文结合双极性混沌序列的内在确定性和托普利兹矩阵的优点,构造了基于双极性混沌序列的托普利兹块状感知矩阵.与其他混沌感知矩阵生成方法不同的是,构造Bi-TpCM 只需要双极性混沌序列,无须对实值混沌序列进行抽样操作.Bi-TpCM 不仅直接继承了托普利兹块状感知矩阵的优势,而且还引入了双极性混沌序列的天然优点.理论分析表明,构造的Bi-TpCM 具有约束等距性,且在自相关性方面接近于最优的采样保证,此外,它在内存开销、计算复杂度和硬件实现等方面具有明显优势.通过对一维信号和图像进行压缩采样,验证了本文构造的Bi-TpCM 相比于其他典型感知矩阵具有更好的性能. 一方面,本文提出的构建Bi-TpCM 框架可推广至不同的混沌系统,如Logistic 和Cat 混沌系统,这样就可以建造出一大族托普利兹块状感知矩阵.另一方面,Bi-TpCM 的构造方式还能派生出Hankel 块、循环块以及堆积块感知矩阵等.读者可以利用类似本文的方法来分析这些混沌感知矩阵的性能.除此之外,可以将Bi-TpCM 和其衍生出的感知矩阵应用于各种各样的多输入-单输出LTI系统的压缩感知测量问题,这也是该方向未来工作的重点. 附录A 定理3 (Geršhgorin 圆盘定理)[22]设矩阵G∈Rn×n,其元素记为Gi,k,1 ≤i,k≤n,那么矩阵G的特征值存在于n个 圆 盘di=di(ci,uk),1 ≤i≤n的 结 合 处,其 中ci=Gi,i表示di的中心,为di的半径.
2 双极性混沌序列及其统计特性
2.1 双极性混沌序列的产生
2.2 双极性混沌序列的统计特性
3 双极性托普利兹块状混沌感知矩阵
3.1 Bi-TpCM 的构造
3.2 Bi-TpCM 的性能分析
3.3 Bi-TpCM vs.传统感知矩阵
4 仿真结果及分析
4.1 约束等距现象
4.2 一维稀疏信号
4.3 图像信号
5 结 论
