基于信息算子正交三角分解的测量矩阵优化方法技术

本发明专利技术涉及一种基于信息算子正交三角分解的测量矩阵优化方法，属于压缩感知测量矩阵优化领域。本发明专利技术的过程是：首先，根据生成的初始随机测量矩阵Φ与选定的信号稀疏基（或稀疏字典）Ψ求得信息算子D，然后对信息算子的转置进行正交三角分解，再将分解后得到的上三角阵的非对角线元素全部置零后更新信息算子，该方法可降低信息算子各列向量之间的相关性，最后根据此优化信息算子求得的最终的优化测量矩阵。应用优化测量矩阵与初始测量矩阵对相同的信号进行压缩感知采样再重构后，应用优化测量矩阵时重构信号的质量更好。

Optimal Method of Measuring Matrix Based on Orthogonal Triangular Decomposition of Information Operator

【技术实现步骤摘要】
基于信息算子正交三角分解的测量矩阵优化方法
本专利技术涉及压缩感知信号处理领域，主要涉及压缩感知中测量矩阵的优化设计。
技术介绍
应用压缩感知理论，可以根据信号的少量线性采样值精确或近似精确恢复原信号。从信号采样频率的角度理解，压缩感知的信号采样频率可以远低于奈奎斯特采样频率，现阶段压缩感知已成功应用于医学上的核磁共振等领域。测量矩阵是压缩感知理论的信号采样矩阵。压缩感知理论诞生之初，测量矩阵多采用正交变换矩阵的随机行抽取矩阵，如傅里叶变换、哈达玛变换的随机行抽取矩阵。另一类应用与研究较广泛的测量矩阵是便于硬件实现的二值矩阵。而等距约束性质指出，测量矩阵与信号的稀疏基的相关性越小越好。高斯随机矩阵被证明与信号的稀疏基之间能够以很大概率满足等距约束性质的要求，所以高斯随机测量矩阵在测量矩阵优化设计的研究中也备受关注。测量矩阵的优化方法之一是降低测量矩阵与信号稀疏基之间的相关性。根据此二者与信息算子之间的关系，该相关性可用信息算子各个列向量之间的相关性表示。所以降低信息算子各列向量之间的相关性可以实现测量矩阵的优化设计。本专利技术根据上述原理，通过将信息算子正交三角分解的上三角阵的非对角线元素置零的方法降低其列向量之间的相关性的方法，来实现测量矩阵的优化设计。
技术实现思路
针对降低测量矩阵与信号稀疏基之间相关性的问题，本专利技术公开了一种基于信息算子正交三角分解的测量矩阵优化方法。该方法能降低信息算子各列向量之间的相关性，实现测量矩阵的优化设计。本专利技术提出的基于信息算子正交三角分解的测量矩阵优化方法，实现步骤为：步骤一、产生随机测量矩阵Φ∈Rm×l(m<;l)，选定信号稀疏基Ψ∈Rl×n(l≤n)；步骤二、计算信息算子D=ΦΨ；步骤三、对信息算子D的转置DT进行正交三角分解DT=QR，其中Q∈Rn×n为正交阵，R∈Rn×m为上三角阵；步骤四、将上三角阵R的非对角线元素置零求得矩阵；步骤五、计算新的信息算子=；步骤六、根据更新后的信息算子求得优化后的测量矩阵。本专利技术优点：本专利技术通过将信息算子的正交三角分解的上三角阵的非对角线元素置零的方法降低信息算子各列向量之间的相关性，实现测量矩阵的优化设计，算法实现步骤简单，不需迭代，运算量小。附图说明图1是本专利技术的实现流程图。图2是使用完备正交稀疏基时，根据初始测量矩阵求得的信息算子各列向量之间的相关性分布与用本专利技术方法对初始测量矩阵的优化结果求得的信息算子各列向量之间的相关性分布比较。图3是使用过完备稀疏基时，根据初始测量矩阵求得的信息算子各列向量之间的相关性分布与用本专利技术的方法对初始测量矩阵的优化结果求得的信息算子各列向量之间的相关性分布比较。图4是使用初始测量矩阵的压缩感知图像重构结果。图5是使用本专利技术的优化测量矩阵的压缩感知图像重构结果。具体实施方式为了更好的描述本专利技术提出的基于信息算子正交三角分解的测量矩阵优化方法，根据图1所示的本专利技术的测量矩阵实现方式，对本专利技术的具体的技术方案予以详细说明。如图1所示，本专利技术的具体步骤如下。步骤一、设置测量矩阵的行数m、列数l，产生初始测量矩阵Φ。选择信号的稀疏基矩阵，确定其行、列数l、n，生成稀疏基矩阵Ψ。以图2的实验为例，初始测量矩阵Φ是元素均值为0，方差为1的高斯随机矩阵，其行数m=100、列数l=256。稀疏基矩阵Ψ也采用元素均值为0，方差为1的高斯随机矩阵，其行、列数分别为l=256，n=256。步骤二、计算信息算子D=ΦΨ。步骤三、应用施密特正交化方法对信息算子D的转置DT进行正交三角分解DT=QR，其中Q∈Rn×n为正交阵，R∈Rn×m为上三角阵。步骤四、将上三角阵R的非对角线元素置零得矩阵，式中diag(Rm)表示由R的前m行m列的对角线元素构成的对角阵。步骤五、求得更新后的信息算子=。步骤六、根据更新后的信息算子求得优化后的测量矩阵。图3的初始测量矩阵Φ是元素均值为0，方差为1的高斯随机矩阵，其行数m=45、列数l=100。稀疏字典Ψ采用离散余弦变换矩阵，其行、列数分别为l=100，n=256。然后采用步骤二到步骤六求解优化后的测量矩阵。图4使用初始测量矩阵Φ为元素均值为0，方差为1的高斯随机矩阵，其行数m=180、列数l=256。稀疏基矩阵Ψ采用Symlets小波变换矩阵，其行、列数分别为l=256、n=256。该图是直接使用初始测量矩阵对图像进行压缩感知采样再重构的结果。图5的初始测量矩阵、稀疏基都与图4的相同，Φ为元素均值为0，方差为1的高斯随机矩阵，其行数m=180、列数l=256。稀疏基矩阵Ψ采用Symlets小波变换矩阵，其行、列数分别为l=256、n=256。对初始测量矩阵应用本专利技术的方法优化后的测量矩阵作为图5的测量矩阵。本文档来自技高网...

【技术保护点】
1.基于信息算子正交三角分解的测量矩阵优化方法，包括以下步骤：步骤一、产生随机测量矩阵

【技术特征摘要】
1.基于信息算子正交三角分解的测量矩阵优化方法，包括以下步骤：步骤一、产生随机测量矩阵Φ∈Rm×l(m<l)，选定信号稀疏基（或稀疏字典矩阵）Ψ∈Rl×n(l≤n)，步骤二、计算信息算子D=ΦΨ；步骤三、对信息算子D的转置DT进行正交三角分解DT=QR，其中Q∈Rn×n为正交阵，R∈Rn×m为上三角阵；步骤四、将上三角阵R的非对角线元素置零得矩阵；步骤五、求得更新后的信息算子=；步骤六、根据更新后的...

【专利技术属性】
技术研发人员：潘金凤，申晋，刘发英，毛帅，马立修，尹丽菊，
申请(专利权)人：山东理工大学，
类型：发明
国别省市：山东,37