多量测压缩感知的感知矩阵构建方法、系统及存储介质与流程

本发明涉及信号处理技术领域，特别涉及一种多量测压缩感知的感知矩阵构建方法、系统及存储介质。

背景技术

压缩感知是一种全新的信号处理方法，其核心思想是通过对信号非自适应、不完全的量测，恢复出原始的稀疏信号。由于压缩感知可以突破奈奎斯特采样定理的限制，因此，广泛应用于数据压缩、图像处理、医学信号处理、信号参数估计等相关领域。

传统的压缩感知在信号的稀疏量测阶段和稀疏重构阶段均采用同一量测矩阵，并通过恢复算法实现对信号的恢复。而所述量测矩阵是一个冗余矩阵，其列向量间存在较强的相关性，当用贪婪类算法对信号进行稀疏重构时会使得信号支撑集的恢复产生错误，造成信号的稀疏重构效果差，信号成功恢复概率低，影响信号重构的性能。

因此，现有技术还有待改进和提高。

技术实现要素：

本发明有必要为了解决现有技术中未利用量测数据以高概率重构原始信号，且在信号稀疏重构时信号支撑集不准确而影响信号重构效果问题，提供一种多量测压缩感知的感知矩阵构建方法、系统及存储介质，旨在使得感知矩阵和量测矩阵对应的列向量间强相关的约束条件下，获得较小的局部累计互相关(lccc，localcumulativecross-coherence)，降低稀疏信号支撑集恢复错误产生率，提升压缩感知中贪婪类恢复算法的性能，从而准确还原原始信号。

本发明解决上述技术问题所采用的技术方案如下：

本发明提供一种多量测压缩感知的感知矩阵构建方法，所述多量测压缩感知的感知矩阵构建方法包括：

获取原始信号的采样数据；

生成随机矩阵作为量测矩阵对采样数据进行稀疏量测，得到量测数据；

根据量测的次数与所述量测矩阵中行数的大小关系，通过待构建的感知矩阵对所述量测数据进行重构，输出对原始信号的估计以实现原始信号的恢复。

所述的多量测压缩感知的感知矩阵构建方法，其中，所述获取原始信号的采样数据具体包括：

获取原始信号；

对所述原始信号进行采样得到采样数据。

所述的多量测压缩感知的感知矩阵构建方法，其中，所述生成随机矩阵作为量测矩阵对采样数据进行稀疏量测，得到量测数据具体包括：

通过软件生成一个随机矩阵作为量测矩阵；

设置量测的次数；

根据量测的次数，通过所述量测矩阵对所述采样数据进行稀疏量测，构建所有量测次数对应的量测数据。

所述的多量测压缩感知的感知矩阵构建方法，其中，所述根据量测的次数与所述量测矩阵中行数的大小关系，通过待构建的感知矩阵对所述量测数据进行重构，输出对原始信号的估计以实现原始信号的恢复具体包括：

获取量测的次数和所述量测矩阵的行数；

比较量测次数与行数的大小；

根据比较结果，构建感知矩阵；

通过所述构建的感知矩阵对所述量测数据进行重构，输出对原始信号的估计。

所述的多量测压缩感知的感知矩阵构建方法，其中，所述根据比较结果，构建感知矩阵具体包括：

当所述行数不大于所述量测次数时，构建第一感知矩阵；

当所述行数大于所述量测次数时，构建第二感知矩阵。

所述的多量测压缩感知的感知矩阵构建方法，其中，所述当所述行数不大于所述量测次数时，构建第一感知矩阵具体包括：

当所述行数不大于所述量测次数时，构建所述量测数据的第一重构模型；

将所述第一重构模型进行优化处理得到第一感知矩阵。

所述的多量测压缩感知的感知矩阵构建方法，其中，所述当所述行数大于所述量测次数时，构建第二感知矩阵具体包括：

当所述行数大于所述量测次数时，构建正则项；

根据所述正则项，构建所述量测数据的第二重构模型；

将所述第二重构模型进行优化处理得到第二感知矩阵。

所述的多量测压缩感知的感知矩阵构建方法，其中，所述优化处理是指求解一定约束条件下第一重构模型或第二重构模型的最优解。

本发明还提供一种系统，所述系统包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行所述多量测压缩感知的感知矩阵构建程序，所述多量测压缩感知的感知矩阵构建程序被所述处理器执行时实现如上述所述的多量测压缩感知的感知矩阵构建方法的步骤。

本发明还提供一种存储介质，所述存储介质存储有多量测压缩感知的感知矩阵构建的程序，所述多量测压缩感知的感知矩阵构建程序被处理器执行时实现上述所述多量测压缩感知的感知矩阵构建方法的步骤。

有益效果：

1.充分利用量测数据，在信号重构阶段，通过构建的感知矩阵替换传统的量测矩阵，减少信号支撑集的恢复的错误产生，提高原始信号估计的准确性。

2.基于量测次数与随机生成的量测矩阵行数关系，灵活选择对应的感知矩阵，使得信号压缩感知过程更具调节性和人为控制，最大程度地还原数据如原始图像。

附图说明

图1是本发明实施例一多量测压缩感知的感知矩阵构建方法的流程图；

图2是本发明实施例m≤l时的稀疏信号支撑集成功恢复概率与稀疏度的仿真结果；

图3是本发明实施例m≤l时恢复信号的均方根误差与信号稀疏度的仿真结果；

图4是本发明实施例m＞l时稀疏信号支撑集成功恢复概率与稀疏度的仿真结果；

图5是本发明实施例m＞l时恢复信号的均方根误差与稀疏度的仿真结果。

图6是本发明系统的较佳实施例的运行环境示意图。

具体实施方式

为使本发明的目的、技术方案及优点更加清楚、明确，以下参照附图并举实施例对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

为使本发明的目的、技术方案及优点更加清楚、明确，以下参照附图并举实施例对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

需要说明的是，本发明是基于压缩感知理论，其处理过程包括三个阶段，分别为信号的稀疏表示、信号的稀疏量测以及信号的稀疏重构，以实现本发明。

本发明提供了一种多量测压缩感知的感知矩阵构建方法，如图1所示，所述多量测压缩感知的感知矩阵构建方法包括：

s10,获取原始信号的采样数据。

具体地，预先进行采样，即获取原始信号，并进行采样得到采样数据，其中，所述原始信号指的是源端向终端发送数据时，通过对应信号相互传递消息，当接收到相应信号时才能知道对方所要表达的消息。例如用户a需要向用户b发送图像，则用户a向用户b发送图像信号(即为原始信号)，用户b接收到该图像信号即开始接收该图像，并向用户a反馈接收图像的信号，从而完成一次完整的数据传输。再如，医生需要探查患者患病部位，通过医学仪器扫描探测的光子转换为电子，形成电脉冲信号(即为原始信号)，经信号分析、数模转换及数据处理等成像。

s20,生成随机矩阵作为量测矩阵对采样数据进行稀疏量测，得到量测数据。

即步骤s20具体包括：

s21,通过软件生成一个随机矩阵作为量测矩阵；

s22,设置量测的次数；

s23,根据量测的次数，通过所述量测矩阵对所述采样数据进行稀疏量测，构建所有量测次数对应的量测数据。

本实施例中，如对获取的图像信号采样后，通过软件生成一个随机矩阵作为量测矩阵对该采样后的图像信号进行稀疏量测，其中，所述量测矩阵用φ表示，(m表示量测矩阵行的个数，n表示量测矩阵列的个数，m和n具体的值由实际的工程问题确定)，所述随机矩阵服从高斯分布，此时将量测矩阵φ对所述采样数据(如图像信号)进行稀疏量测得到若干个对应的量测数据，经过预设的量测次数l后，将得到l个量测数据，然后基于每一次量测数据，构建采样信号稀疏量测后的多向量量测模型(mmv，multiplemeasurementvectors)，如式(1)所示：

y＝φx+n(1)

其中，y＝[y1y2…yl]表示量测数据矩阵，表示第l个量测向量，x＝[x1x2…xl]表示每一次量测时对应的采样信号构成的集合，简称联合稀疏信号，其中，该x中只有某些行的元素为非零值而其它行的元素均为零，表示x中非零行序号构成的集合表示稀疏信号的支撑集，n表示量测噪声，表示量测矩阵，m表示量测矩阵行的个数，n表示量测矩阵列的个数，并且行数远远小于列数，即m＜＜n；l＝1,2,…,l表示对联合稀疏信号x的量测次数，在第l次量测时，第l次稀疏信号的支撑集(即采样信号)为xl，量测数据为yl，l次量测后集合得到量测数据矩阵y。

s30,根据量测的次数与所述量测矩阵中行数的大小关系，通过待构建的感知矩阵对所述量测数据进行重构，输出对原始信号的估计以实现原始信号的恢复。

即步骤s30具体包括：

s31,获取量测的次数和所述量测矩阵的行数；

s32,比较量测次数与行数的大小；

s33,根据比较结果，构建感知矩阵；

s34,通过所述构建的感知矩阵对所述量测数据进行重构，输出对原始信号的估计。

进一步的，本实施例中，步骤s33具体包括：

s331,当所述行数不大于所述量测次数时，构建第一感知矩阵；

s332,当所述行数大于所述量测次数时，构建第二感知矩阵。

本发明中，在压缩感知信号的重构阶段，采用感知矩阵ψ替换传统压缩感知的量测矩阵φ来重构信号，用于恢复原始信号。

具体地，通过量测次数l与上述步骤s20生成的量测矩阵φ中行数m大小关系来构建对应的感知矩阵。

1)当m≤l时，定义待构建的第一感知矩阵ψ1，用待构建的第一感知矩阵ψ1对量测数据进行重构，以得到原始信号的估计值即得到所述量测数据的第一重构模型是如式(2)所示：

其中，ψ1·i表示第一感知矩阵的第i列形式，表示的第i行。(·)^t表示矩阵的转置操作。

为了使得量测数据重构后能逼近原始信号，得到准确的原始信号的估计值，需要保证式(2)中(其中i＝1,2,…,n)无畸变输出，即要求同时为了抑制或降低其它干扰和噪声的影响，要求将公式(2)等式右面最小化，即式子(2)右边转换为因此，将所述第一重构模型进行优化处理就得到第一感知矩阵ψ1。而所述优化处理指的是求解一定约束条件下第一重构模型的最优解，即取最小值。本实施例中，将该优化处理转化为最小二范数优化问题求解，即求解在的约束条件下，目标函数的最小值，如式(3)所示，这样也就可以凭借量测数据y通过求解优化问题而精确恢复出原始信号，进而求解得到第一感知矩阵ψ1。

其中，φ·i表示量测矩阵的第i列，i＝1,2,…,n，(·)^t表示矩阵的转置操作，||·||2表示向量的二范数，min(·)表示取最小值操作，表示约束条件为第一感知矩阵ψ1转置后的列向量与量测矩阵φ对应的列向量乘积为1，也就是说，保证了待构建的第一感知矩阵ψ1与所述量测矩阵φ的局部累计相关(lccc，localcumulativecross-coherence)最小化，降低了第一感知矩阵ψ1和量测矩阵φ对应的列向量间强相关的影响。

公式(3)用于表示在第一感知矩阵ψ1转置后的列向量与量测矩阵φ对应的列向量乘积为1时目标函数即求取最小值，求解出的最小值即为待构建的第i列的形式的第一感知矩阵ψ1，即m≤l下如式(4)所示：

其中，r1＝yy^t表示量测数据矩阵y的第一协方差矩阵，i＝1,2,…,n；上标-1表示矩阵的取逆操作。

此时，通过迭代计算式(4)，即i取不同值构成的集合以得到完整的第一感知矩阵ψ1。

最后，终端通过已知的第一感知矩阵ψ1对量测数据进行重新构建，恢复出原始信号，即还原输出原始数据，如原始图像。

如图2和图3所示，图2和图3示例了当m≤l时的仿真结果，以验证本发明技术方案，其仿真条件如下：信噪比snr＝20db；为量测矩阵，计算机仿真中设定量测矩阵行的个数m为128，列的个数n为256，其中的元素服从均值为零方差为一的高斯分布，即高斯随机矩阵；为了获得统计性能，每次实验独立重复500次，即l＝500；采样信号的稀疏度(sparsityofsignal，记为k)从5至100逐渐递增；为了方便对比，同时给出了传统压缩感知(ψ＝φ)、交替投影法apm、重加权算法rwa的仿真结果；本发明所采用的恢复算法是联合正交匹配追踪(simultaneousorthogonalmatchingpursuit，somp)算法。其中，proposedmethod对应本发明所采用的算法结果。

图2示例了在snr＝20db，l＝500下稀疏信号支撑集成功恢复概率随稀疏度(k)的变化情况。图2中，横坐标表示稀疏度k，纵坐标表示稀疏信号支撑集成功恢复概率。随着稀疏度k的增加，四种算法对稀疏信号支撑集恢复的成功率都呈下降趋势。传统压缩感知(ψ＝φ)、交替投影法apm、重加权算法rwa随着k的增加相继失效，当k＝70时，本发明提出的算法依然能够以100％的概率恢复出稀疏信号的支撑集，说明了本发明提出方法更具有有效性，更准确还原出原始信号，重构效果更佳。

图3示例了在snr＝20db，l＝500下恢复信号(即原始信号，也即联合稀疏信号重构后的信号)的均方根误差随稀疏度的变化情况。图3中，横坐标表示稀疏度k，纵坐标表示恢复信号的均方根误差。随着稀疏度k的增加，传统压缩感知(ψ＝φ)、交替投影法apm、重加权算法rwa以及本发明提出方法重构出的稀疏信号，均方根误差都相继升高，但是由本发明所提算法重构出的信号均方根误差最小，说明了在同等条件下，本发明所提方法更保证了还原原始数据的准确性与完整性，重构效果更佳。

2)当m＞l时，定义待构建的第二感知矩阵ψ2，采用与m≤l时第一感知矩阵ψ1对量测数据进行重构的同样原理，用待构建的第二感知矩阵ψ2对量测数据进行重构，即构建所述量测数据的第二重构模型，得到原始信号数据的估计值也就得到第i列的形式的第二感知矩阵其中，r2＝yy^t表示量测数据矩阵y的第二协方差矩阵。

但是，当m＞l时，yy^t为欠定矩阵，由于逆矩阵不存在，因此，无法基于1)式(3)中同样原理求解出第二感知矩阵ψ2，此时构建正则项，即用于解决当m＞l时出现的欠定矩阵yy^t的逆矩阵不存在的问题。其中，diag(·)表示将向量构造为一个对角矩阵，λ表示正则化参数。

进一步地，基于1)中转换为求解最小二范数优化问题相同的原理，即可以凭借量测数据y和正则项通过求解优化问题，如式(5)所示而精确恢复出原始信号，进而求解得到第二感知矩阵ψ2。

其中，diag(·)表示将向量构造为一个对角矩阵，λ表示正则化参数，λ的取值位于区间(0,1)中，用户可根据实际情况进行调整其值；||·||2表示向量的二范数，min(·)表示取最小值操作，s.t.表示约束条件为第二感知矩阵ψ2转置后的列向量与量测矩阵φ对应的列向量乘积为1，也就是说，保证了待构建的第二感知矩阵ψ2与所述量测矩阵φ的局部累计相关lccc最小化，降低了第二感知矩阵ψ2和量测矩阵φ对应的列向量间强相关的影响。

公式(5)表示第二感知矩阵ψ2转置后的列向量与量测矩阵φ对应的列向量乘积为1，并且使得目标函数即取最小值，求解得到m＞l下第二感知矩阵ψ2的第i列的形式，如式(6)所示：

其中，r2＝yy^t+λφww^tφ^t表示正则化后量测数据矩阵y的第二协方差矩阵，矩阵r2中yy^t表示欠定矩阵，λφww^tφ^t表示用于解决矩阵yy^t的欠定问题的构建的第二正则项，上标-1表示取逆操作。

此时，通过迭代计算式(6)即i取不同值构成的集合以得到完整的第二感知矩阵ψ2。

最后，终端通过已知的第二感知矩阵ψ2对量测数据进行重新构建，恢复出原始信号，即还原输出原始数据，如原始图像。

如图4和图5所示，图4和图5示例了m＞l下的仿真结果，其仿真条件如下：信噪比snr＝20db；为量测矩阵，计算机仿真中设定量测矩阵行的个数m为128，列的个数n为256，其中的元素服从均值为零方差为一的高斯分布，即高斯随机矩阵；为了获得统计性能，每次实验独立重复500次，即l＝500；信号的稀疏度(sparsityofsignal，记为k)从5至100逐渐递增；为了方便对比，同时给出了传统压缩感知(ψ＝φ)、交替投影法apm、重加权算法rwa的仿真结果；本发明所采用的恢复算法是联合正交匹配追踪(simultaneousorthogonalmatchingpursuit，somp)算法。其中，proposedmethod对应本发明所采用的算法结果。

图4示例了在snr＝20db，l＝50，λ＝0.8下稀疏信号支撑集成功恢复概率随稀疏度的变化情况。图4中，横坐标表示稀疏度k，纵坐标表示稀疏信号支撑集成功恢复概率。随着稀疏度的增加，四种方法对稀疏信号支撑集恢复的成功率都呈下降趋势。传统压缩感知(ψ＝φ)、交替投影法apm、重加权算法rwa随着k的增加相继失效，当k＝40的时候本发明提出的算法依然能够以100％的概率恢复出稀疏信号的支撑集，说明了在同等条件下，本发明所提方法更保证了还原原始数据的准确性与完整性，重构效果更佳。

图5示例了在snr＝20db，l＝50，λ＝0.8下恢复信号(即原始信号，也即联合稀疏信号)的均方根误差随稀疏度的变化情况。图5中，横坐标表示稀疏度k，纵坐标表示恢复信号的均方根误差。随着稀疏度的增加，传统压缩感知(ψ＝φ)、交替投影法apm、重加权算法rwa以及本发明提出方法重构出的稀疏信号，均方根误差都相继升高，但是k≤60时由本发明所提算法重构出的信号均方根误差最小，说明了本发明所提方法的有效性，更能降低稀疏信号支撑集恢复错误，重构效果更优。当k＞60时，本发明所提方法的均方根误差大于rwa算法的均方根误差。但是通过图5可知当k＞60时四种算法均不能正确恢复稀疏信号的支撑集，因此考虑k＞60时的均方根误差没有意义。

实施例二

进一步地，如图6所示，基于上述多量测压缩感知的感知矩阵构建方法，本发明还相应提供了一种系统，所述系统包括处理器10、存储器20及显示器30。图6仅示出了系统的部分组件，但是应理解的是，并不要求实施所有示出的组件，可以替代的实施更多或者更少的组件。

所述存储器20在一些实施例中可以是所述系统的内部存储单元，例如系统的硬盘或内存。所述存储器20在另一些实施例中也可以是所述系统的外部存储设备，例如所述系统上配备的插接式硬盘，智能存储卡(smartmediacard,smc)，安全数字(securedigital,sd)卡，闪存卡(flashcard)等。进一步地，所述存储器20还可以既包括所系统的内部存储单元也包括外部存储设备。所述存储器20用于存储安装于所述系统的应用软件及各类数据，例如所述安装系统的多量测压缩感知的感知矩阵构建程序代码等。所述存储器20还可以用于暂时地存储已经输出或者将要输出的数据。在一实施例中，存储器20上存储有多量测压缩感知的感知矩阵构建程序40，该多量测压缩感知的感知矩阵构建程序40可被处理器10所执行，从而实现多量测压缩感知的感知矩阵构建方法。

所述处理器10在一些实施例中可以是一中央处理器(centralprocessingunit,cpu)，微处理器或其他数据处理芯片，用于运行所述存储器20中存储的多量测压缩感知的感知矩阵构建程序代码或处理数据，例如执行所述多量测压缩感知的感知矩阵构建方法等。

所述显示器30在一些实施例中可以是led显示器、液晶显示器、触控式液晶显示器以及oled(organiclight-emittingdiode，有机发光二极管)触摸器等。所述显示器30用于显示在所述系统的信息以及用于显示可视化的用户界面。所述系统的部件10-30通过系统总线相互通信。

在一实施例中，当处理器10执行所述存储器20中多量测压缩感知的感知矩阵构建程序40时实现以下步骤：

获取原始信号的采样数据；

生成随机矩阵作为量测矩阵对采样数据进行稀疏量测，得到量测数据；

根据量测的次数与所述量测矩阵中行数的大小关系，通过待构建的感知矩阵对所述量测数据进行重构，输出对原始信号的估计以实现原始信号的恢复，具体如上述s10-s30所述。

实施例三

本发明还提供一种存储介质，所述存储介质存储有多量测压缩感知的感知矩阵构建程序，所述多量测压缩感知的感知矩阵构建程序被处理器10执行时实现上述所述多量测压缩感知的感知矩阵构建方法的步骤，具体如上所述。

综上所述，本发明公开了一种多量测压缩感知的感知矩阵构建方法、系统及存储介质，所述多量测压缩感知的感知矩阵构建方法包括：获取原始信号的采样数据；生成随机矩阵作为量测矩阵对采样数据进行稀疏量测，得到量测数据；根据量测的次数与所述量测矩阵中行数的大小关系，通过待构建的感知矩阵对所述量测数据进行重构，输出对原始信号的估计以实现原始信号的恢复。本发明通过将感知矩阵设计问题转化为带有约束条件的最小二范数优化问题，使得感知矩阵和量测矩阵对应的列向量间强相关以实现二者的局部累积互相关最小化，从而使得接收后的数据在重构后能够还原输出准确且完整的原始数据，提高了数据恢复的成功率和准确性。

当然，本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关硬件(如处理器，控制器等)来完成，所述的程序可存储于一计算机可读取的存储介质中，所述程序在执行时可包括如上述各方法实施例的流程。其中所述的存储介质可为存储器、磁碟、光盘等。

应当理解的是，本发明的应用不限于上述的举例，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，所有这些改进和变换都应属于本发明所附权利要求的保护范围。