一种维数任意的部分哈达玛测量矩阵构造方法及基于其的信号处理方法
【专利摘要】本发明公开了一种维数任意的部分哈达玛测量矩阵构造方法及基于其的信号处理方法，测量矩阵构造方法先取得待测稀疏或可压缩信号x的维数N，然后判断N的大小，若N＝2k，则测量矩阵ΦM×N(M<N)由N阶哈达玛矩阵HN的前M行构成；否则，将l个阶数分别为2k1，2k2，…，2kl(k1、k2、…、kl为互不相等的非负整数，且2k1+2k2+…+2kl＝N)的哈达玛矩阵子块对角排列成矩阵ΦN×N，并取其前M行作为测量矩阵ΦM×N(M<N)；最后利用ΦM×N对信号x进行测量得到M维测量值，再用匹配追踪类重建算法获得重建信号。本发明适用范围广，占用存储空间小，降低了硬件实现成本。
【专利说明】
一种维数任意的部分哈达玛测量矩阵构造方法及基于其的信 号处理方法
技术领域
[0001] 本发明属于信息技术领域中的信号压缩处理方法，特别涉及信号压缩感知中一种 维数任意的部分哈达玛测量矩阵构造方法及基于其的信号处理方法。
【背景技术】
[0002] 随着信息技术的快速发展，人们对信息的需求量也以惊人的速度增长，以奈奎斯 特采样定理为支撑的传统信息获取、处理、传输和存储方式面临着前所未有的挑战，这主要 是因为：1)奈奎斯特米样定理要求的米样频率过尚，有效彳目息的提取效率低下；2)先米样后 压缩的方法造成了大量的数据冗余和存储资源浪费。因此，寻找高效且能降低存储和传输 成本的信息处理方法是必然的。
[0003] 2004年，Donoho，Cand6s和Tao等学者提出了压缩感知理论，该理论指出，如果信号 在一组正交基下是稀疏或可压缩的，那么就可以用一个与稀疏矩阵不相关的矩形测量矩阵 对其进行观测，直接获取信号的压缩测量值，然后通过求解非线性优化问题高概率重建出 原信号。
[0004] 测量矩阵是压缩感知理论中实现信号采样的关键。2006年，Cand6s，Romberg和Tao 等学者在研究欠定的核磁共振成像时发现，当部分傅里叶矩阵作为测量矩阵时，只需要0 (K ? log(N))的测量值数目便能将稀疏度为K的N维信号精确重建。但是，部分傅里叶矩阵只 适用于时域或频域稀疏的信号，而大部分的自然信号都不是时域或频域稀疏的，所以这种 矩阵不具有普遍适用性。随后，Cand6s和Tao等学者证明，独立同分布的高斯随机矩阵可以 观测任意维数的可压缩信号，并能获得较高的重建质量。为了指导构造更多的测量矩阵， 2007年Cand6s和Tao等提出了著名的有限等距性质(RIP)，即对于任意长度为N的K稀疏信号 x和常数S，如果测量矩阵〇满足：
[0005] (1-^)||(D|' <|〇x|g <(1 + ^)|〇>|' , ^€ (0,1)
[0006] 则称它满足K阶有限等距性质。理论上讲，矩阵〇只要满足2K阶有限等距性质，就 能保证把任意一个K稀疏信号x映射到唯一的采样集合y中。但是，要证明一个矩阵是否满足 有限等距性质并不容易。为了解决该问题，Baraniuk给出了有限等距性质的等价条件一一 不相干性。测量矩阵〇和稀疏矩阵W的不相干性通过互相关系数y来定义，其数学表达式 为：
[0008]式中，ai为矩阵A第i列的向量，这里A= 〇 W，称为感知矩阵。互相关系数y越小，测 量值里包含的原始信号的信息越多，信号的重建效果就越好。除此之外，Donoho在论文 Compressedsensing中给出压缩感知概念的同时，也定性和定量地给出了测量矩阵要满足 的三个特征:①由测量矩阵的列向量组成的子矩阵的最小奇异值必须大于一定的常数;② 测量矩阵的列向量体现某种类似噪声的独立随机性;③满足稀疏度的解是满足1范数最小 的向量。
[0009] 在这些理论的指导下，研究者们先后构造出了三类测量矩阵：
[0010] 1)随机测量矩阵，如高斯随机矩阵、贝努利随机矩阵等，相应地，矩阵元素独立地 服从高斯分布、贝努利分布。这类矩阵与绝大多数的稀疏矩阵均不相关，且能很好地满足有 限等距性质。但是，它们也有着明显的不足，如矩阵元素无规律，随机性太强，计算复杂度 高，硬件实现困难，在仿真实验中尚且需要多次测量再取平均值。
[0011] 2)确定性测量矩阵，如多项式测量矩阵、伪随机测量矩阵等。多项式测量矩阵是最 早提出的一种确定性测量矩阵，定义在有限域上，解决了随机测量矩阵稳定性差、硬件实现 繁琐的问题，缺点是测量值数目不能任意选取，矩阵构造时间较长，重建效果较差;伪随机 测量矩阵由伪随机序列生成，这种矩阵的重建效果介于随机测量矩阵和多项式测量矩阵之 间，缺点是矩阵元素均为浮点数，需要占用较大的存储空间。
[0012] 3)结构随机测量矩阵，如托普利兹测量矩阵、部分哈达玛测量矩阵等。托普利兹测 量矩阵的特点是主对角线上的元素相等，平行于主对角线的各线上的元素也相等。实际中， 这种矩阵只用一个简单的移位寄存器便可实现，但缺点是重建效果差。部分哈达玛测量矩 阵由N阶哈达玛矩阵的任意M行构成，不仅构造简单，在现有测量矩阵中其重建效果也是最 好的，但由于哈达吗矩阵的维数~=2 1{仏=1，2,3，...），这极大地限制了其使用范围，因为 图像处理中经常需要选择感兴趣区域(R0I)进行分析，这种情况下的图像维数具有任意性。 此外，该矩阵元素均为非零值，需要较大的存储空间，硬件实现的成本相对高。
[0013] 因此，有必要设计一种适用范围广、占用空间小的哈达玛测量矩阵。
【发明内容】
[0014] 本发明的目的就是要解决将部分哈达玛矩阵作为压缩感知测量矩阵时，信号不能 为任意维数且测量阵占用内存空间大的问题，提供一种维数任意的部分哈达玛测量矩阵的 构造方法及基于其的信号处理方法，适用范围广，占用空间小。
[0015] 本发明的具体实现步骤为：
[0016] 步骤一:通过直接输入或者计数或者使用测量函数等方式，取得待测的原始稀疏 或可压缩信号x的维数，记为N;
[0017] 步骤二:若N = 2k，则构造 N阶哈达玛矩阵Hn，然后选取Hn的前M行构成测量矩阵〇MXN (1〈《;否则，将1(1为大于等于2的整数)个阶数分别为21{1，21{2，~，2 1{1仏1、1^、一、1^1为互不 相等的非负整数，且2kl+2k2+~+2 kl = N)的哈达玛矩阵作为子块序列，并将此子块序列沿主 对角线排列构成块对角矩阵Onxn，然后选取?NXN的前M行作为测量矩阵C>mxn(M〈N);顺序选 取矩阵Hn前M行而不是随机选取任意M行，更好地保证了测量矩阵的行正交性以及测量结果 的稳定性;M的取值由采样率决定；
[0018] 步骤二中1个阶数不同的哈达玛矩阵子块的构造步骤如下：
[0019] a:按哈达玛矩阵构成原理构造哈达玛矩阵子块序列中的第一个子块出，其阶数2kl (kl=0，l，2,3....)为2的整数次幂中小于信号维数N的最大值；
[0020] b:按哈达玛矩阵构成原理构造哈达玛矩阵子块序列中的第二个子块出，其阶数为 2k2(k2 = 0，l，2,3....);若N-2kl是2的整数次幂，即令2k2 = N-2k1，结束哈达玛矩阵子块的构 造;否则，2k2取2的整数次幂中小于N-2kl的最大值；
[0021] c:按哈达玛矩阵构成原理构造哈达玛矩阵子块序列中的第三个子块H3，其阶数为 2k3(k3 = 0，l，2,3....);若N-2kl-2k2是2的整数次幂，则令2k3 = N-2kl-2k2,结束哈达玛矩阵子 块的构造;否则，2k3取2的整数次幂中小于N-2kl-2k2的最大值；
[0022] d:按步骤b和c的方式，逐个构造哈达玛矩阵子块序列的后续低阶子块，直至第1个 哈达玛矩阵子块Hi(l = 2,3,4,...)的阶数2kl = N-2kl-2k2--2k(1_1)〇
[0023] 上述1个阶数不同的哈达玛矩阵子块序列出，出，…，Hr沿主对角线排列形成块对角 矩阵C>nxn，Onxn如下式所不： 丨 ~
[0024] ^ 。 _
[0025] 所述步骤二中测量矩阵〇 MXN为Emxn^nxn，其中，Emxn；^N|^单位矩阵的前赠丁〇
[0026]阶数为n的哈达玛矩阵Hn是以1和-1为元素构成的并且满足HnHn = nIn的方阵，其 中，InSn阶单位矩阵，Hn是其中的任意两行或两列均相互正交的矩阵。
[0027]基于上述的维数可控的部分哈达玛测量矩阵的信号处理方法，利用步骤二得出的 测量矩阵〇MXN对步骤一所述信号x进行测量，即将测量矩阵〇MXN与待测的N维原始稀疏或 可压缩的信号x相乘，得到M维测量值;然后通过匹配追踪类重建算法（如正交匹配追踪算 法），获得重建信号。
[0028] 与现有技术相比，本发明的有益效果是：
[0029] 本发明提出的方法突破了以部分哈达玛矩阵作为压缩感知测量矩阵时，信号维数 不能为任意值的限制，扩大了其适用范围，同时，因采用分块对角结构，避免了稠密矩阵占 用存储空间大的缺点，降低了硬件实现成本。
[0030] -、针对目前重建效果最好的部分哈达玛测量矩阵适用范围受限而且窄的缺点， 本发明提出了一种维数任意的部分哈达玛测量矩阵的构造方法，使得按该方法构造的测量 矩阵可以观测任意维数的稀疏或可压缩信号，并获得较高的重建质量；
[0031] 二、针对稠密矩阵占用存储空间大的缺点，本发明在构造过程中采用了块对角阵 结构，增加了矩阵的稀疏性，降低了内存需求和硬件实现成本；
[0032]三、针对随机矩阵稳定性差的缺点，本发明顺序选取矩阵Onxn的前M行而不是随机 选取任意M行，使得压缩测量的结果更加稳定。
【附图说明】
[0033]图1是本发明提出的一种维数任意的部分哈达玛测量矩阵构造方法及基于其的信 号处理方法的实现流程图；
[0034]图2是采样率为0.5时，采用本发明提出的测量矩阵对大小为120X 120的Lena、 Cameraman、Fruits和Boat原始图像进行压缩测量后用相同重建算法得到的重建图像；图2 (a)、（b)分别为Lena原始图像和重建图像;图2(c)、（d)分别为Cameraman原始图像和重建图 像；图2(e)、（f)分别为Fruits原始图像和重建图像；图2(g)、（h)分别为Boat原始图像和重 建图像；
[0035]图3是采样率为0.5时，采用本发明提出的测量矩阵对大小为200X200的Lena、 Cameraman、Fruits和Boat原始图像进行压缩测量后用相同重建算法得到的重建图像；图3 (a)、（b)分别为Lena原始图像和重建图像;图3(c)、（d)分别为Cameraman原始图像和重建图 像；图3(e)、（f)分别为Fruits原始图像和重建图像；图3(g)、（h)分别为Boat原始图像和重 建图像；
[0036]图4是分别采用高斯随机矩阵、Logistic混沌测量矩阵、托普利兹测量矩阵和本发 明提出的测量矩阵，对大小为120X120的Lena图像进行压缩测量后以相同重建算法进行重 建所得采样率和重建图像的峰值信噪比(PSNR)的关系曲线图。
【具体实施方式】
[0037] 下面结合附图对本发明的【具体实施方式】做进一步的描述。
[0038] 如附图1所示，本发明提出的一种维数任意的部分哈达玛测量矩阵构造方法，具体 实现步骤为：
[0039] 步骤一、使用测量函数取得待测原始稀疏或可压缩信号x的维数，记为N;
[0040] 步骤二、判断N的大小，若N=2k，贝IJ构造 N阶哈达玛矩阵Hn，然后选取Hn的前M行构成 测量矩阵cDmxn(M〈N);否则，将1(1为大于等于2的整数)个阶数分别为2 kl，2k2,…，2kl(kl、 k2、…、kl为互不相等的非负整数，且2kl+2k2+…+2kl = N)的哈达玛矩阵作为子块序列，并将 此子块序列沿主对角线排列形成块对角矩阵〇nxn，然后选取〇NXN的前M行作为测量矩阵 OmxnCKN)；
[0041] 基于上述的维数可控的部分哈达玛测量矩阵的信号处理方法，利用步骤二所得测 量矩阵〇 mxn(M〈N)对步骤一中的信号x进行测量，得到M维测量值，然后通过匹配追踪类重建 算法，获得重建彳目号；
[0042] 步骤二中2k阶哈达玛矩阵的构造方式为： r//..., //,, 1 ? ^
[0043] H0k =Tj-~ .左=其中，Hi= (1)。
[0044] 步骤二中1个阶数不同的哈达玛矩阵子块的构造方法如下：
[0045] a:按哈达玛矩阵构成原理，构造哈达玛矩阵子块序列中的第一个子块出，其阶数 2kl(kl = 0，l，2,3....)为2的整数次幂中小于信号维数N的最大值；
[0046] b:按哈达玛矩阵构成原理，构造哈达玛矩阵子块序列中的第二个子块H2，其阶数 为2k2(k2 = 0，1，2，3....);若N-2kl是2的整数次幂，即令2k2 = N-2kl，结束哈达玛矩阵子块的 构造;否则，2k2取2的整数次幂中小于N-2 kl的最大值；
[0047] c :按哈达玛矩阵构成原理，构造哈达玛矩阵子块序列中的第三个子块H3，其阶数 为2k3(k3 = 0，l，2,3....);若N-2kl-2k2是2的整数次幂，则令2k3 = N-2kl-2k2,结束哈达玛矩阵 子块的构造;否则，2k3取2的整数次幂中小于N-2 kl-2k2的最大值；
[0048] d:按步骤b和c的方式，逐个构造哈达玛矩阵子块序列的后续低阶子块，直至第1个 哈达玛矩阵子块Hi(l = 2,3,4,...)的阶数2kl = N-2kl-2k2--2k(1_1)〇
[0049] 步骤二中矩阵Onxn由哈达玛矩阵子块序列(出，出，…，H〇沿主对角线排列形成块 对角矩阵：
[0050]边關= ' " 杜.2 -
[0051] 步骤二中顺序选取矩阵〇nxn前M行而不是随机选取任意M行是为了更好的保证测 量矩阵的行正交性以及测量结果的稳定性。
[0052]下面通过具体实施例来验证本发明所提方法的有效性。需要指出的是，该实施例 只是示例性的，并不是要限制本发明的适用范围。
[0053] 实施例1:
[0054] 选取大小为120 X 120的Lena、Cameraman、Fruits和Boat图像，并对其进行压缩重 建。具体过程为：
[0055] 一、采用离散余弦变换(DCT)矩阵分别对上述三幅图像进行稀疏化；
[0056] 二、构造测量矩阵，按照本发明提供的一种维数可控的部分哈达玛测量矩阵构造 方法，120可以表示为26+2 5+24+23 ,所以各哈达玛矩阵子块的阶数分别为64,32，16,8，按从 小到大顺序沿主对角线排列得到：
[0057] ，編= 好32 _
[0058] 然后选取?12QX12Q的前M行作为最后的测量矩阵，其中，M的取值由采样率决定； [0059] 三、在采样率M/N为0.5时，对上述四幅图像进行压缩测量，然后采用正交匹配追踪 算法对其恢复重建。
[0060]附图2给出了分辨率为120x120的原始图像和通过压缩感知重建的对应图像，从中 可以看出，采样率为0.5时，本发明提出的测量矩阵构造方法可以实现对120 X 120图像的压 缩重建。
[0061 ] 实施例2:
[0062] 选取大小为200 X 200的Lena、Cameraman、Fruits和Boat图像，并对其进行压缩重 建。其具体实现过程同实施例1，其中因200可表示为27+26+23 ,所以对应的哈达玛矩阵子块 的阶数分别为128，64，8。附图3给出了分辨率为200 X 200原始图像和通过压缩感知重建的 对应图像，从中可以看出，采样率为0.5时，本发明提出的测量矩阵构造方法同样可以实现 对200X200图像的压缩重建。另外，该方法还可以推广到对任意维信号的压缩重建。需要补 充说明的是，图3显示的重建图像质量优于图2显示的重建图像质量，是因为图3的原始图像 的分辨率高于图2的原始图像。
[0063] 实施例3:
[0064]分别采用高斯随机矩阵、Logistic混沌测量矩阵、托普利兹测量矩阵和本发明提 出的测量矩阵对大小为120X120的Lena图像进行压缩测量的性能比较，这里是将峰值信噪 比作为性能比较的评价指标。若峰值信噪比大，则认为重建效果好。从图4中对应不同压缩 测量矩阵的采样率与重建信号的峰值信噪比的关系曲线可以看出，采样率不超过0.8时，采 用本发明提出的测量矩阵进行压缩感知的图像重建效果明显优于其他测量矩阵的。对于采 样率大于0.8以后的情况，采用其他测量矩阵测量后重建的效果会略好，但高采样率有违采 用压缩感知方法的本意，因此，在要采用压缩感知方法的前提下以高采样率来改进重建性 能没有实际意义，故本例不考虑采样率大于0.8的情况。从实际应用的角度来说，本发明构 造的测量矩阵进行压缩感知的图像重建效果优于其他测量矩阵。
【主权项】
1. 一种维数任意的部分哈达玛测量矩阵构造方法，其特征在于，该方法包括以下步骤： 步骤一:取得待测的原始稀疏或可压缩的信号X的维数，记为N; 步骤二:若N = 2k，则构造 N阶哈达玛矩阵Hn，然后选取Hn的前M行构成测量矩阵Φ MXN，M〈 N;否则，将1个阶数分别为2kl，2k2，…，2kl的哈达玛矩阵作为子块序列，并将此子块序列按对 角排列构成矩阵Φνχν，然后选取Φνχν的前M行作为测量矩阵Φ ΜΧΝ，Μ〈Ν;其中，1为大于等于2 的整数，1^1、1^2、~、1^1为互不相等的非负整数，且21{1+2 1{2+~+21{1 = 12. 根据权利要求1所述的维数任意的部分哈达玛测量矩阵构造方法，其特征在于，所述 步骤二中1个阶数不同的哈达玛矩阵子块序列的构造方法如下： 步骤a:按哈达玛矩阵构成原理，构造哈达玛矩阵子块序列中的第一个子块出，其阶数 2kl为2的整数次幂中小于信号维数N的最大值;kl = 0，l，2,3....为非负整数； 步骤b:按哈达玛矩阵构成原理，构造哈达玛矩阵子块序列中的第二个子块出，其阶数为 2k2，k2 = O，1，2，3....为非负整数;若N-2kl是2的整数次幂，即令2k2 = N-2kl，结束哈达玛矩阵 子块序列的构造;否则，2k2取2的整数次幂中小于N-2 kl的最大值； 步骤c:按哈达玛矩阵构成原理，构造哈达玛矩阵子块序列中的第三个子块H3，其阶数为 2k3，k3 = O，1，2，3....为非负整数;若N-2kl-2k2是2的整数次幂，则令2k3 = N-2kl-2k2，结束哈 达玛矩阵子块序列的构造;否则，2k3取2的整数次幂中小于N-2 kl-2k2的最大值； 步骤d:按步骤b和c的方式，逐个构造哈达玛矩阵子块序列的后续低阶子块，直至第1个 哈达玛矩阵子块Hi，其阶数为2kl = N-2kl-2k2--3. 根据权利要求2所述的维数任意的部分哈达玛测量矩阵构造方法，其特征在于，1个 阶数不同的哈达玛矩阵子块序列H1，H 2，…，H1沿主对角线排列形成块对角矩阵ΦNXN，ΦNXN 如下式所示：4. 根据权利要求1~3中任一项所述的维数任意的部分哈达玛测量矩阵构造方法，其特 征在于，通过直接输入、计数或者使用测量函数的方式，取得待测的原始稀疏或可压缩的信 号X的维数，记为N。5. -种基于权利要求1~3中任一项所述的维数任意的部分哈达玛测量矩阵的信号处 理方法，其特征在于，利用步骤二得出的测量矩阵Φμχν对步骤一所述信号X进行测量，即将 测量矩阵Φ ΜΧΝ与待测的N维原始稀疏或可压缩的信号X相乘，得到M维测量值;然后通过匹配 追踪类重建算法，获得重建信号。
【文档编号】G06T11/00GK105931282SQ201610232004
【公开日】2016年9月7日
【申请日】2016年4月14日
【发明人】刘少强, 闫晓艳, 樊晓平, 徐文, 李飞
【申请人】中南大学