Loading [MathJax]/jax/output/HTML-CSS/jax.js
高级搜索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于采样值随机压缩矩阵核空间的亚奈奎斯特采样重构算法

盖建新 杜昊辰 刘琦 童子权

盖建新, 杜昊辰, 刘琦, 童子权. 基于采样值随机压缩矩阵核空间的亚奈奎斯特采样重构算法[J]. 电子与信息学报, 2019, 41(2): 484-491. doi: 10.11999/JEIT180323
引用本文: 盖建新, 杜昊辰, 刘琦, 童子权. 基于采样值随机压缩矩阵核空间的亚奈奎斯特采样重构算法[J]. 电子与信息学报, 2019, 41(2): 484-491. doi: 10.11999/JEIT180323
Jianxin GAI, Haochen DU, Qi LIU, Ziquan TONG. Sub-Nyquist Sampling Recovery Algorithm Based on Kernel Space of the Random-compression Sampling Value Matrix[J]. Journal of Electronics & Information Technology, 2019, 41(2): 484-491. doi: 10.11999/JEIT180323
Citation: Jianxin GAI, Haochen DU, Qi LIU, Ziquan TONG. Sub-Nyquist Sampling Recovery Algorithm Based on Kernel Space of the Random-compression Sampling Value Matrix[J]. Journal of Electronics & Information Technology, 2019, 41(2): 484-491. doi: 10.11999/JEIT180323

基于采样值随机压缩矩阵核空间的亚奈奎斯特采样重构算法

doi: 10.11999/JEIT180323
基金项目: 国家自然科学基金(61501150),黑龙江省自然科学基金(QC2014C074)
详细信息
    作者简介:

    盖建新:男,1980年生,博士,副教授,研究方向为压缩感知、亚奈奎斯特采样理论、频谱感知技术等

    杜昊辰:男,1991年生,硕士生,研究方向为电子与通信工程

    刘琦:男,1994年生,硕士生,研究方向为仪器仪表工程

    童子权:男,1968年生,教授,研究方向为现代电子测量仪器与系统、信号处理等

    通讯作者:

    盖建新 gjx800608@126.com

  • 中图分类号: TP391

Sub-Nyquist Sampling Recovery Algorithm Based on Kernel Space of the Random-compression Sampling Value Matrix

Funds: The National Natural Science Foundation of China (61501150), The Natural Science Foundation of Heilongjiang Province (QC2014C074)
  • 摘要:

    针对现有调制宽带转换器亚奈奎斯特采样重构算法性能不高问题,该文提出一种基于采样值核空间的支撑重构算法和随机压缩降秩方法,将两者结合得到一种高性能采样重构算法。首先利用随机压缩变换在不改变未知矩阵稀疏特性的前提下将采样方程转化为多个新的多测量向量问题,然后利用采样值矩阵核空间与采样矩阵支撑正交的关系获取联合稀疏支撑集,最后通过伪逆完成重构。从理论和实验两个方面对所提方法进行了分析和验证。数值实验表明,与传统重构算法相比,所提算法提高了重构成功率、降低了高概率重构所需的通道数,而且重构性能总体上随压缩次数增加而提高。

  • 以奈奎斯特采样理论为基础的传统信息处理技术,目前在雷达电子监视、认知无线电频谱感知等应用中正面临着采样率过高的压力[13]。近年来,压缩感知(Compressive Sensing, CS)理论[4,5]提出了一种充满前景的信息获取方法。该理论指出,对于稀疏或可压缩的信号,当测量矩阵满足约束等距特性时,可以从少量低速采样中恢复出原信号。CS理论掀起了亚奈奎斯特信息处理技术的发展新高潮。近年来,基于CS理论,学者们提出了随机解调[6]、随机采样[7]、多陪集采样[8]、调制宽带转换[9,10]等一系列亚奈奎斯特采样方法。研究的焦点主要集中在采样处理和重构算法两个方面。在这些方法中,调制宽带转换器(Modulated Wideband Converter, MWC)亚奈奎斯特采样,以其所适用的稀疏宽带信号模型应用广泛、重构框架运算负担轻、可采用现有商用元器件实现等优势,受到了雷达通信、医疗成像、认知无线电等多个领域学者的广泛关注[11]

    MWC的信息获取过程由采样和重构两个阶段组成。在采样阶段采用多个通道的结构来获取稀疏模拟信号,通道数量决定着硬件的复杂程度,并影响重构成功率。如何采用更少的通道,以更高成功率从低速数据流中重构原始信号,是MWC的核心问题[12]。MWC的重构模型可归结为多测量向量(Multiple Measurement Vectors, MMV)问题,属于CS理论的一个分支。MMV重构算法大致可以分为贪婪算法[13]和凸松弛算法[14]两大类。前者运算耗时较长难以达到测量系统对实时性的要求,后者具有理论保障且执行效率较高,近年来发展快、应用多。在MWC提出之初,贪婪类的正交匹配追踪(Orthogonal Matching Pursuit, OMP)算法被用于MWC重构。在一定条件下,以较高的概率实现了精确重构。文献[15]对迭代硬门限(IHT)、正则化IHT(NIHT)、硬门限追踪(IHP)、压缩感知匹配追踪(CoSaMP)等4种贪婪类单测量向量(Single Measurement Vector, SMV)算法进行了推广,从理论保障方面进行了推导论证,得出了与之对应的MMV重构算法:SIHT, SNIHT, SIHP, SCoSaMP,其中SCoSaMP算法与其他贪婪算法相比重构性能有较大提高。文献[16]在剖析CS与阵列信号处理之间关系的基础上,将多信号分类(MUSIC)思想在压缩感知框架下进行了改进,提出了可压缩MUSIC(CSMUSIC)算法,该算法充分利用CS与MUSIC算法各自的优点,即使源信号相干时,只要能够通过压缩感知算法估计出特定的部分支撑集,便可以确定性地得到剩余支撑集。但鉴于MWC重构问题的具体特点(如随机频谱分割、采样矩阵与未知矩阵均共轭对称),并不是所有高性能MMV重构算法都可以在MWC重构过程中发挥出优势。上述算法当用于MWC重构时其性能距离理论上限仍存在较大的改进空间。文献[17]针对一般的MMV问题提出了一种降维多测量向量与性能提升(Reduce MMV and Boost, ReMBo)重构算法,该算法将MMV问题随机地组合成SMV问题进行求解,通过多次尝试的方法来提升重构性能,使得ReMBo算法的重构性能在无噪声条件下呈现出了优势,然而在有噪声条件下用于MWC重构时,该算法立即失效。

    本文为提高MWC重构算法性能,将核空间思想引入MWC重构算法。为满足核空间(Kernel Space, KS)条件,针对采样值矩阵的秩较大的问题,提出K列随机压缩(Randem Compression, RC)降低维度的方法,在不改变未知矩阵的稀疏性的前提下,实现了降秩功能。将两种思想相结合最终得到一种名为RCKS的高性能重构算法。理论分析和实验结果验证了本文算法的有效性。

    设信号x(t)为实值连续时间信号,带限于一个较宽的频率范围F=[fNYQ/2,fNYQ/2]内(fNYQ为奈奎斯特率),F内仅存在N个非零子频带,子频带可分布在F中的任意位置,子频带最大宽度小于等于B Hz。该类信号被称为稀疏宽带信号。图1给出了典型稀疏宽带信号频谱示意图。

    图 1  稀疏宽带信号频谱示意图

    MWC采样前端由m个通道组成,输入信号x(t)同时进入m个通道,在第i通道信号x(t)被周期为Tp(fp=1/Tp)的伪随机信号pi(t)混频。混频后的信号在频域被截止频率为1/2T的低通滤波器h(t)截断,最后利用采样率为fs=1/T的低速ADC采样,获取m组低速采样值序列yi[n],i{1,2,...,m}。如图2所示。

    图 2  MWC系统框图

    从频域分析MWC第i个通道可得

    DTFT(yi[n])=Yi(ej2πfT)=λ0n=λ0ci,nX(fnfp),fFS (1)

    式中,DTFT表示离散时间傅里叶变换,ci,npi(t)的傅里叶级数的系数,X(f)x(t)的连续时间傅里叶变换,λ0=[(fNYQ+fs)/2fp]1, []表示向上取整,FS[fs/2,fs/2]

    从式(1)可以看出,MWC采样相当于原信号频谱X(f)fp为步进进行移位,并由滤波器截断成宽度为fs的频谱片段,最后以ci,n,n[λ0,λ0]为系数进行线性组合。综合考虑m个通道,令

    y(f)=[Y1(ej2πfT),···,Yi(ej2πfT),···,Ym(ej2πfT)]T, z(f)=[X(f+λ0fp),···,X(f),···,X(fλ0fp)]T,则式(1)可表示为

    y(f)=Φz(f),fFS (2)

    式中,Φi,n=ci,nλ01,文中称其为采样矩阵,维数为m×λ(λ=2λ0+1),通常m<λ

    式(2)两边同时右乘傅里叶反变换矩阵,可得到对应的时域形式

    Y(n)=ΦZ(n),n[1,L] (3)

    其中,Y(n)=[y1[n],y2[n],···,ym[n]]T, Z(n)=DTFT1(z(f)), L为采样记录长度。

    假设图1为输入信号x(t)的频谱,fs=fpB,则MWC的时域采样方程如图3所示,图中DTFT–1表示对z(f)的每个行向量取DTFT逆运算。可见只要由采样方程得到z(f),通过频谱搬移等处理方法立即恢复出图1所示的原始信号。

    图 3  MWC采样方程示意图

    MWC重构是指从采样值矩阵Y(n)通过特定算法恢复出未知矩阵z(f)Z(n)的过程。由于m<λ,式(2)和式(3)均为欠定性方程组,常规方法有无穷多解。考虑z(f)Z(n)具有联合稀疏特性,即只有少量不全为零的行,当满足定理1条件时,存在唯一的稀疏解。

    定理 1[9] 设多带信号x(t)N个子频带构成,各子频带最大带宽为B,按照图2所示的MWC结构进行采样,如果以下条件成立:

    (1) fsfpB并且fs/fp数值不是很大;

    (2) 一个周期内序列pi(t)的符号(±1)数M2[fNYQ2fp+12]1

    (3)m2N

    (4) 矩阵Φ的任意2N列线性无关。

    则对于fFS, z(f)是式(2)唯一的N-稀疏解。

    定理1满足时,MWC重构过程可分成两个步骤进行:首先通过稀疏优化等算法获得支撑集Ω=supp(Z(n)),然后通过式(4)伪逆运算计算未知矩阵。

    ZΩ(n)=ΦΩY(n)Zi(n)=0,iΩ} (4)

    其中,ΦΩ=(ΦHΩΦΩ)1ΦHΩΦΩ的左伪逆矩阵。MWC提出时将OMP作为MMV问题的支撑重构算法,在一定条件下实现了精确重构。然而OMP算法在重构成功率等性能方面距离理论上限还存在较大的改进空间。针对该问题本文将核空间思想引入到MWC重构过程中。为了满足核空间思想的重构条件,提出一种采样值K列随机压缩降秩变换思想。

    前已述及,MWC重构的关键在于求解MMV问题未知矩阵的联合稀疏支撑集Ω。下面提出一种利用采样值矩阵的核空间矩阵来判断支撑集的方法。

    Y=AX为某一MMV问题,其中,采样值矩阵YRm×K, Y的秩rank(Y)=K,采样矩阵ACm×λ, A的Kruskal秩σ(A)K+1,即A的任意K+1列均线性无关,X的联合稀疏支撑集定义为Ω={i|Xi0}|Ω|=K,其中,Xi表示以i为索引X的行子集,|Ω|表示集合Ω的势(cardinality)。

    考虑支撑集的定义,采样值矩阵Y可以表示为Y=AΩXΩ,其中AΩ表示以Ω为索引A的列子集。设QY各列所张成的空间的核空间矩阵,即有QTY=0成立。于是可得QTY=QTAΩXΩ=0。根据前提假设易知K=rank(Y)min[rank(AΩ),rank(XΩ)],即有rank(XΩ)K,结合X的联合稀疏度为K,因此rank(XΩ)=K,即XΩ可逆,于是有QTAΩ=0成立,即对于支撑集中的任意元素jΩ均有QTAj=0成立。由此可见,核空间思想为我们提供了一种可能的支撑集判决条件。问题是满足QTAj=0的所有A的列索引都是支撑集的元素吗?

    假设除了jΩ外,还存在j0Ω使QTAj0=0成立,则有QTA[Ω,j0]=0成立。然而由QTY=0可知Q的秩rank(QT)=mK,于是由QTA[Ω,j0]=0可推出rank(A[Ω,j0])=K,这与前提条件σ(A)K+1A的任意K+1列均线性无关相矛盾,因此综上所述,当且仅当jΩQTAj=0成立。

    需要注意的是,核空间思想假设采样值矩阵的秩与未知信号的稀疏度相同,然而在MWC实际系统中,一般采样长度较长,即Y包含的列向量较多,导致rank(Y)K图4给出了不同条件下采样值矩阵的秩rank(Y)的统计结果。如图4所示,信号的频带数N=6时(此时稀疏度为12),在采样长度分别为L=13,20,40,60,80时,通道数m以3为步进在区间[13,78]内变化,在不同信噪比(Signal to Noise Ratio, SNR)条件下,只要m13,采样值矩阵的秩尽管随着m有波动,但均大于未知矩阵的稀疏度12。

    图 4  不同条件下采样值矩阵的秩随通道数的变化情况

    前面分析中可知若想采用核空间思想,必须保证Y的秩等于K,若rank(Y)>Krank(Q)<mK,造成集合{j|QTAj=0}的势大于K以至于最终按照核空间方法求出的支撑元素数量会增加,即|Ω|K,导致重构失败。为此,需要寻找一种矩阵降秩方法,在不影响未知矩阵稀疏性的前提下达到压缩维度的目的。

    文献[17]为了获取支撑集,提出了一种名为ReMBo的随机线性组合方法将MMV问题变成SMV问题来求解,同时证明了随机组合前后MMV问题的未知矩阵与SMV问题的未知向量的(联合)支撑集不变。受此启发,本文提出一种K列随机压缩变换方法,采用符合特定分布的K列“窄”随机矩阵将采样值矩阵进行线性变换,进而达到降秩的目的。本节引入K列随机压缩变换方法有两个目的:一是降低采样值矩阵的秩以满足核空间思想的条件,进而利用采样值矩阵的核空间判断支撑集合;二是将MWC采样方程压缩变换成与之具有相同支撑的多个新的MMV问题,利用新的MMV问题间接地求解原MWC采样方程的支撑集,并在有可能的条件下提高重构性能。ReMBo方法中,线性组合相当于采样值矩阵及未知矩阵各行向量在随机向量上进行投影,若采用K个随机向量构成的矩阵作为变换矩阵,可将采样值矩阵压缩变换成m×K的新矩阵。需要注意的是压缩变换后MMV问题的唯一解的稀疏性不应发生改变,下面采用定理的形式描述如下:

    定理 2 设ˆX为MMV问题Y=AX的唯一K稀疏解,其中YRm×L, σ(A)2K。令T={T1,T2,···,TK}K(K<m)个长度为L的单位范数随机列向量,其元素符合某一特定连续分布。令ˉY=YT, ˉX=^XT,则ˉX为新MMV问题ˉY=AU的唯一K稀疏解,且supp(ˉX)=supp(ˆX)成立的概率为1。

    证明 令ˉyn=YTn, ˉxn=ˆXTn,其中n{1,2,···,K}。由文献[17]可知,ˉxn为SMV问题ˉyn=Au的唯一K稀疏解,且P[supp(ˆX)=supp(ˉxn)]= 1成立,即supp(ˆX)=supp(ˉx1)=supp(ˉx2)···=supp(ˉxK)成立的概率为1。由ˉY=YT={ˉy1,ˉy2,···,ˉyK}, ˉX=^XT={ˉx1,ˉx2,···,ˉxK},因此ˉY=AˉX成立。结合联合支撑集定义知:

    supp(ˉX)=supp{ˉx1,ˉx2,···ˉxK}=n{1,2,···,K}supp(ˉxn)=supp(ˉxn) (5)

    于是有supp(ˉX)=supp(ˉxn)=supp(ˆX)成立概率为1,且|supp(ˉX)|=|supp(ˆX)|=K,即ˉX联合稀疏度为K。考虑ˉY=AˉX, |supp(ˉX)|=K, σ(A)2K,由文献[14]中定理2.2知ˉXˉY=AU的唯一K稀疏解。 证毕

    定理2保证了K列随机压缩后新MMV问题稀疏解的支撑集与原MMV问题稀疏解的支撑集相同。因此可以通过求解压缩变换后的MMV问题间接求解原MMV问题。

    ˉY=YT知,理论上rank(ˉY)rank(T)=K。在MWC实际采样中,由于信号频谱分割的随机性及频谱泄露效应,大量实验表明实际采样值矩阵在K列随机压缩后的秩恰好为K图5给出了在不同的通道数m、不同采样长度L及不同信噪比条件下,当信号的频带数为6时(稀疏度为12)随机压缩后采样值矩阵ˉY的秩的情况。从图5中可以看出,由于随机矩阵的满秩特性及MWC处理信号的实际情况,随机压缩变换保证了变换后采样值矩阵的秩与未知矩阵的稀疏度相同(K=12)。可见通过K列随机压缩的方法可以达到应用核空间方法的基本条件,K列随机压缩方法为核空间思想应用于MWC重构奠定了基础。

    图 5  不同条件下随机压缩后采样值矩阵秩的统计结果

    本文提出的随机压缩方法,采用K列随机矩阵作为变换矩阵,鉴于随机性,每次重构时可以利用多个随机矩阵进行多次压缩变换,得到多个不同的MMV问题,进而增加MMV问题的多样性。问题是,每次重构时的支撑集判据如何确定。由核空间思想可知无噪声时采样值矩阵的核空间与以支撑集为索引的采样矩阵列子集正交,对于每次压缩得到的MMV问题可以采用2范数的方法通过计算||QTAj||2是否为0来判断j是否属于支撑集。当然有噪声时可以通过判断||QTAj||2的大小的方法,通过排序来确定支撑集合。一次重构过程进行r次压缩变换时,考虑到各MMV问题解的联合稀疏性特征,可以将r个支撑集判据综合起来,将和式ri=1||Q(i)TAj||2, j{1,2,···,λ}作为最终排序依据。由于该式反映了r个MMV问题采样值矩阵与采样矩阵之间的正交关系,即反映了r个未知矩阵之间的公共信息(联合稀疏信息),当个别随机压缩过程不利于重构时,由于综合条件ri=1||Q(i)TAj||2保留了MMV问题的多样性,仍然有希望完成最终重构,进而提高判决的成功概率。

    本文提出的MWC重构算法结合了K列RC与KS思想,简称为RCKS算法,具体算法描述如下。

    输入:采样值矩阵Y(n)Rm×L,采样矩阵ΦCm×λ,频带数量N,随机压缩次数r

    输出:支撑集估计ˆΩ,未知矩阵估计ˆX(n)

    步骤 1 取r个随机矩阵作为压缩变换矩阵T(i)=randn(L,2N)i{1,2,···,r},实现r次压缩变换ˉY(i)=YT(i),其中T(i)表示第i个随机矩阵,ˉY(i)表示第i个压缩后的采样值矩阵;

    步骤 2 求解ˉY(i)的核空间矩阵Q(i), i{1,2,···r}

    步骤 3 计算Pj=ri=1||Q(i)TΦj||2, j{1,2,···,λ}

    步骤 4 从向量P=[P1P2···Pλ]中找到2N个最小值所对应的索引值,即支撑集ˆΩ=min(P,2N)

    步骤 5 按照式(4)伪逆运算完成重构,返回未知矩阵的估计ˆX(n)

    为了验证所提出的RCKS算法的有效性,本节设计了3个实验:首先考察K列随机压缩的次数r对重构性能的影响;其次在相同条件下对比RCKS算法及当前典型MMV算法用于MWC重构时的成功率情况;最后给出RCKS算法用于频谱感知时的重构效果。为了便于说明问题,实验中稀疏宽带信号由学者们普遍采用的通用表达式产生[9,12]

    x(t)=Nn=1EnBnsinc(Bn(tτn))cos(2πfn(tτn))+n(t) (6)

    其中,参数N, En, Bn, fn, τn分别代表子频带数、第n个子频带的能量系数、子带宽度、载波频率和延迟时间,n(t)为高斯白噪声。每个实验中,以下过程重复500次,将成功次数的百分率作为成功概率:

    (1) 在pi(t)的一个周期内,等概率随机产生±1作为调制波形;

    (2) 在区间[–fNYQ/2, fNYQ/2]内按均匀分布随机产生各子频带的载波频率fn

    (3) 用重构算法估计支撑集,当估计支撑集与实际支撑集相同即ˆΩ=Ω时,则称重构成功,如果ˆΩΩΦˆΩ列满秩,也作为成功标志[9]

    首先,考察RCKS算法中随机压缩次数r对重构性能的影响。前已述及,本文提出的RCKS算法采用随机压缩实现降秩功能,对于同一个重构问题,该过程可以随机进行多次,而且次数r对算法重构性能有影响。为此,本实验检验在不同的r值时,RCKS的重构成功率情况。不失一般性,以具有6个(对称的3对)子频带的稀疏宽带信号为例。信号参数设置如下:SNR=30 dB; En={1, 2, 3}; Bn={50, 50, 50} MHz; τn={0.394, 0.985, 1.379} μs;载波频率fn随机地分布在[–5, 5] GHz内。MWC采样参数设置如下:λ=2λ0+1=195; fs=fp=fNYQ/λ=51.28MHz。图6给出了RCKS算法的重构成功率随通道数m的变化情况,图中分别给出了压缩变换次数r=1, 2, 3, 4, 6, 8, 10时所对应的曲线。如图6所示,当12<m<20时,总体上随着r增大,RCKS的重构成功率具有逐渐提高的趋势。如当m=14, r=1, 2, 3, 4, 6, 8, 10时,成功率分别为36.0%, 60.6%, 66.2%, 74.4%, 80.2%, 82.8%, 82.6%。为了清楚地看到该特点,图7统计了m=14, 15, 16, 17时重构成功率随压缩次数的变化曲线。从图7可以看出,当r数值较小时(如r4时),随着r增大,重构性能改进较大,当r较大时,重构成功率并不是无限度的提高,当r5时,曲线开始上下波动,此时重构性能提高较小或没有提高。可见采用RCKS算法实现MWC重构时,随机压缩的次数并不需要太多,后文取r=4

    图 6  不同压缩次数时RCKS重构性能随通道数的变化
    图 7  RCKS重构成功率随压缩次数的变化

    其次,比较在相同条件下RCKS(r=4)与典型的MMV算法(CSMUSIC, SCoSaMP, OMPMMV)用于MWC重构时的重构性能。本实验中,被测信号的信噪比分别设置为无噪声,40 dB, 30 dB和20 dB,信号其它参数及MWC采样参数与前面实验设置相同。如图8所示,在各种信噪比下,当通道数m<25时在重构成功率方面RCKS(r=4)与其它算法相比均表现出了较大的优势。如当m=14,信号中无噪声及SNR=40 dB, 30 dB, 20 dB时,RCKS(r=4)的重构成功率分别比CSMUSIC高37.4%, 35.8%, 28.6%, 7.2%,比SCoSaMP 高75.0%, 69.8%, 57.4%, 18.2%,比OMPMMV高86.6%, 81.0%, 61.4%, 21.8%。此外从图中还可以看出,信号无噪声及SNR=40 dB, 30 dB, 20 dB时,RCKS(r=4)高概率(大于99%)重构所需要的最小通道数分别比OMPMMV少9, 5, 4, 2。由于通道数与硬件开销及总体采样率有直接关联,因此RCKS算法可降低MWC系统的硬件复杂度。

    图 8  不同信噪比时RCKS(r = 4)与CSMUSIC, SCoSaMP, OMPMMV重构成功率比较

    最后,为了进一步说明本文算法有效性,图9展示了RCKS(r=4)算法用于亚奈奎斯特频谱感知时的重构效果。本实验中N=20,fn={2.30, 2.50, 2.90, 4.00, 4.04, 4.08, 3.00, 3.10, 3.30, 4.58} GHz, En={5, 4, 4, 2, 2, 2, 4, 3, 4, 4}, Bn={40, 37, 48, 40, 50, 40, 15, 6, 6, 50} MHz, τn={ 0.319, 0.798, 0.957, 0.128, 0.160, 0.479, 0.798, 1.117, 1.197, 1.037} μs, SNR = 25 dB, m=50,MWC其它采样参数与前面实验相同。从图9可以看出,尽管原信号的奈奎斯特率高达9.16 GHz,而MWC等效采样率仅为mfNYQ/λ=2.3487 GHz,采用RCKS(r=4)算法仍然能够实现重构,在子频带位置、频带宽度及幅度等方面均实现了精确恢复,进一步说明了本文算法的有效性。

    图 9  RCKS(r = 4)算法重构效果

    针对现有MWC重构算法性能不高的问题,本文结合随机压缩与核空间思想,提出一种基于采样值随机压缩矩阵核空间的重构算法。该算法首先对采样值矩阵进行多次K列随机压缩,在不改变未知量稀疏特性的前提下实现降秩变换,同时获得多个MMV问题,然后将多个基于核空间思想的支撑集判据综合起来,实现最终支撑重构。实验结果表明,与现有CSMUSIC, SCoSaMP, OMPMMV等典型重构算法相比,本文提出的算法在相同条件下具有重构成功率高、高概率重构所需通道数少等优点。

  • 图  1  稀疏宽带信号频谱示意图

    图  2  MWC系统框图

    图  3  MWC采样方程示意图

    图  4  不同条件下采样值矩阵的秩随通道数的变化情况

    图  5  不同条件下随机压缩后采样值矩阵秩的统计结果

    图  6  不同压缩次数时RCKS重构性能随通道数的变化

    图  7  RCKS重构成功率随压缩次数的变化

    图  8  不同信噪比时RCKS(r = 4)与CSMUSIC, SCoSaMP, OMPMMV重构成功率比较

    图  9  RCKS(r = 4)算法重构效果

  • PARK J, JANG J, IM S, et al. A sub-nyquist radar electronic surveillance system[J]. IEEE Access, 2018, 6: 10080–10091. doi: 10.1109/ACCESS.2018.2799304
    QIN Zhijin, GAO Yue, PLUMBLEY M D, et al. Wideband spectrum sensing on real-time signals at sub-Nyquist sampling rates in single and cooperative multiple nodes[J]. IEEE Transactions on Signal Processing, 2016, 64(12): 3106–3117. doi: 10.1109/TSP.2015.2512562
    郑仕链, 杨小牛. 用于调制宽带转换器压缩频谱感知的重构失败判定方法[J]. 电子与信息学报, 2015, 37(1): 236–240. doi: 10.11999/JEIT140127

    ZHENG Shilian and YANG Xiaoniu. A reconstruction failure detection scheme for modulated wideband converter based compressed spectrum sensing[J]. Journal of Electronics &Information Technology, 2015, 37(1): 236–240. doi: 10.11999/JEIT140127
    DONOHO D L. Compressed sensing[J]. IEEE Transaction on Information Theory, 2006, 52(4): 1289–1306. doi: 10.1109/TIT.2006.871582
    BARANIUK R G. More is less: Signal processing and the data deluge[J]. Science, 2011, 331(6018): 717–719. doi: 10.1126/science.1197448
    LASKA J N, KIROLOS S, DUARTE M F, et al. Theory and implementation of an analog-to-information converter using random demodulation[C]. IEEE International Symposium on Circuits and Systems, New Orleans, USA, 2007: 1959–1962.
    FARDAD M, SAYEDI S M, and YAZDIAN E. Hardware implementation of iterative method with adaptive thresholding for random sampling recovery of sparse signals[J]. IEEE Transactions on Very Large Scale Integration (VLSI) Systems, 2018, 26(5): 867–877. doi: 10.1109/TVLSI.2018.2791351
    REN Shiyu, ZENG Zhimin, GUO Caili, et al. A low complexity sensing algorithm for wideband sparse spectra[J]. IEEE Communications Letters, 2017, 21(1): 92–95. doi: 10.1109/LCOMM.2016.2616471
    MISHALI M and ELDAR Y C. From theory to practice: Sub-Nyquist sampling of sparse wideband analog signals[J]. IEEE Journal of Selected Topics in Signal Processing, 2010, 4(2): 375–391. doi: 10.1109/JSTSP.2010.2042414
    ADAMS D, ELDAR Y C, and MURMANN B. A mixer front end for a four channel modulated wideband converter with 62-dB blocker rejection[J]. IEEE Journal of Solid-State Circuits, 2017, 52(5): 1286–1294. doi: 10.1109/JSSC.2017.2647941
    LEXA M A, DAVIES M E, and THOMPSON J S. Reconciling compressive sampling systems for spectrally sparse continuous-time signals[J]. IEEE Transactions on Signal Processing, 2012, 60(1): 155–171. doi: 10.1109/TSP.2011.2169408
    ZHANG Ruoyu, ZHAO Honglin, JIA Shaobo, et al. Sparse multi-band signal recovery based on support refining for modulated wideband converter[C]. IEEE 13th International Conference on Signal Processing (ICSP), Chengdu, China, 2016: 304–309.
    COTTER F S, RAO D B, ENGAN K, et al. Sparse solutions to linear inverse problems with multiple measurement vectors[J]. IEEE Transactions on Signal Processing, 2005, 53(7): 2477–2488. doi: 10.1109/TSP.2005.849172
    CHEN Jie and HUO Xiaoming. Theoretical results on sparse representations of multiple-measurement vectors[J]. IEEE Transactions on Signal Processing, 2006, 54(12): 4634–4643. doi: 10.1109/TSP.2006.881263
    BLANCHARD J D, CERMAK M, HANLE D, et al. Greedy algorithms for joint sparse recovery[J]. IEEE Transactions on Signal Processing, 2014, 62(7): 1694–1704. doi: 10.1109/TSP.2014.2301980
    KIM J M, LEE O K, and YE J C. Compressive MUSIC: Revisiting the link between compressive sensing and array signal processing[J]. IEEE Transactions on Information Theory, 2012, 58(1): 278–301. doi: 10.1109/TIT.2011.2171529
    MISHALI M and ELDAR Y C. Reduce and boost: Recovering arbitrary sets of jointly sparse vectors[J]. IEEE Transactions on Signal Processing, 2008, 56(10): 4692–4702. doi: 10.1109/TSP.2008.927802
  • 期刊类型引用(2)

    1. 罗沅,党娇娇,宋祖勋,王保平. 认知无线电中实现自适应压缩频谱感知. 系统工程与电子技术. 2020(01): 15-22 . 百度学术
    2. 王兵锐,张新刚,马晓普. 通信教学研究中带通采样的统一推导与验证. 现代信息科技. 2019(24): 55-57 . 百度学术

    其他类型引用(5)

  • 加载中
图(9)
计量
  • 文章访问数:  2169
  • HTML全文浏览量:  888
  • PDF下载量:  57
  • 被引次数: 7
出版历程
  • 收稿日期:  2018-04-11
  • 修回日期:  2018-10-29
  • 网络出版日期:  2018-11-08
  • 刊出日期:  2019-02-01

目录

/

返回文章
返回