一种基于改进StOMP的压缩感知重构算法制造技术

本发明专利技术提出一种基于改进StOMP的压缩感知重构算法。首先，基于StOMP算法更新阈值，并对阈值大小进行判断，由判断结果选择下一步执行回溯原子选择算法或ROMP算法，即当阈值大于判断标准，执行回溯原子选择算法；当阈值小于判断标准，执行ROMP算法。其次，以相邻两次重构信号能量差作为算法转换条件和停止条件；当回溯原子选择算法迭代至满足算法转换条件时执行ROMP算法，执行ROMP至满足停止条件时输出重构信号。本发明专利技术利用StOMP算法收敛速度快，回溯思想原子可信赖性高，以及ROMP重构精度高的优点，将回溯原子选择算法和ROMP算法嵌入StOMP中，有效实现了重构精度及收敛速度的折中。

【技术实现步骤摘要】
一种基于改进StOMP的压缩感知重构算法
：本专利技术涉及无线通信领域，尤其涉及压缩感知领域中一种基于StOMP的改进压缩感知重构算法。
技术介绍
：压缩感知理论是一种新兴的信号压缩采样技术。该理论中的核心问题是信号的重构问题，常用的重构方法主要有贪婪追踪算法、凸松弛算法和组合算法这三大类。贪婪追踪算法由于算法结构简单、计算量小等优点备受关注，而正交匹配追踪(OrthogonalMatchingPursuit，OMP)类算法为其主流，更成为研究者研究的重点，国内外许多研究者对此类算法进行了研究与改进。贪婪追踪算法是通过贪婪迭代的方法来更新支撑集，逐步逼近原始解，目前常用的贪婪算法有OMP算法、分段正交匹配追踪(StagewiseOrthogonalMatchingPursuit，StOMP)算法、正则化正交匹配追踪(RegularizedOrthogonalMatchingPursuit，ROMP)算法、压缩采样匹配追踪(CompressiveSamplingMatchingPursuit，CoSaMP)算法、子空间追踪(SubpuistPursuit，SP)算法、稀疏度自适应匹配追踪(SparsityAdaptiveMatchingPursuit，SAMP)算法等。StOMP算法因每次迭代根据阈值选择得到的是多个匹配原子而不是单个原子，减少了迭代次数，相对于OMP算法进行了一定程度的简化，使运行时间为常用匹配算法里较少的一种，且它不需要知道稀疏度。由于并没有运用稀疏度这一先验信息，且原子的选择跟阈值设置紧密相关，固其重构精度不够理想。SP算法在OMP的基础上，引入回溯思想，为了提高算法的收敛速度和算法效率，通过回溯的思想从原子库里选择多个相关原子同时剔除部分不相关原子，在每次迭代中保证了原子的可信赖度；但其对稀疏度依赖性大，若错误估计了稀疏度的值，算法精确重构的能力会下降。ROMP算法在OMP的基础上改进了原子的选择标准，通过正则化过程对原子进行了二次选择。它结合了贪婪算法的速度和凸优化方法的强有力的理论保证，给出了不同终止准则下的重构误差的上限。ROMP算法的运行时间与OMP算法的运行时间在理论上是相当的，且需要估计稀疏度，但它具有较高的重构精度。综上，对于现在已经出现的压缩感知重构算法，均存在重构精度低或复杂度高的问题，这就亟需本领域技术人员解决相应的技术问题。
技术实现思路
：本专利技术旨在至少解决现有技术中存在的技术问题，特别创新地提出了一种基于改进StOMP的压缩感知重构算法。为了实现本专利技术的上述目的，本专利技术提供了一种基于改进StOMP的压缩感知重构算法，其特征在于，包括：S1，基于分段正交匹配追踪算法，根据阈值大小选择执行回溯原子选择算法或ROMP算法。S2，引入子空间正交匹配追踪算法的回溯思想，即回溯原子选择算法，使每次选择的原子个数等于分段正交匹配追踪算法中选出的列数。S3，将回溯原子选择算法和ROMP嵌入StOMP中，根据相邻两次重构信号能量差判断是否进行算法转换或者迭代停止。所述的一种基于改进StOMP的压缩感知重构算法，其特征在于，所述S1包括：考虑原信号且为可稀疏信号，根据压缩感知理论，将其投影到稀疏基空间上可表示为：X＝ψθ其中，为稀疏系数，含有k个非零元素，则称原信号X是k稀疏的；为信号X的稀疏基。选择随机高斯矩阵作为观测矩阵Φ，其内元素均满足高斯分布N(0,1/n)，则N维信号可由M个观测值表示：Y＝ΦX＝Φψθ＝Aθ其中，是信号的观测值，为传感矩阵。令r0＝Y，并以h＝β×max(th)作为判断标准，确定下一步执行回溯原子选择算法或ROMP算法，其中th＝tsσs，根据经验值取ts＝2.5，1/n为高斯观测矩阵中元素所服从的高斯分布的方差。则根据StOMP，以rt表示残差，t表示迭代次数，表示空集，Jo表示每次迭代找到的列序号，Λt表示t次迭代的列序号集合，aj表示矩阵A的第j列，At表示按索引Λt选出的传感矩阵A的列集合。首先计算u＝<rt-1,aj>,1≤j≤N，选u中大于阈值th＝tsσs的c列，并记对应的A列序号j构成集合Jo；然后，令Λt＝Λt-1∪Jo，At＝At-1∪aj(j∈Jo)。再则令ts＝2.5，判断th＝tsσs与h的大小关系后，选择执行回溯原子选择算法或ROMP算法。当th＞h时，将StOMP选出的At执行回溯原子选择算法；当th＜h时，将StOMP选出的At执行ROMP算法。所述的一种基于改进StOMP的压缩感知重构算法，其特征在于，所述S2包括：当th＞h时，将StOMP选出的At执行回溯原子选择算法，引入回溯思想，确保原子可靠性，选出的稀疏系数估计值个数等于分段正交匹配选出的原子列数c。首先，求Y＝Atθt的最小二乘解：为第t次迭代估计得到的稀疏系数；再从中选出绝对值最大的c项，记为对应的At中的c列记为Atc，对应的列序号记为Λtc，更新集合Λt＝Λtc；其次，更新残差利用稀疏矩阵可得重构信号：所述的一种基于改进StOMP的压缩感知重构算法，其特征在于，所述S3包括：将回溯原子选择算法和ROMP嵌入StOMP中，根据相邻两次重构信号能量差的大小确定算法的转换或迭代的停止。在StOMP中，阈值th是残差的函数，随迭代过程中残差的变小而变小，则以阈值th反映重构的进度。当残差值较大，满足th＞h时，迭代得到的重构信号离目标值较远，此时执行回溯原子选择算法，迭代一次得到c个有效原子，迭代至相邻两次重构信号的能量差小于ε1时停止执行回溯原子选择算法，转而执行ROMP算法。当残差值较小，满足th＜h时，迭代得到的重构信号较靠近目标值，此时选择ROMP算法，迭代一次得到不大于c个有效原子，迭代至相邻两次重构信号的能量差小于ε2时，停止迭代，并输出重构信号作为最终结果。当th＞h时，执行回溯原子选择算法；迭代至相邻两次重构得到的信号能量差满足：时，停止执行回溯原子选择算法转而执行ROMP算法。当th＜h时，执行ROMP算法，对原子进行正则化；首先，构建子空间，使其内的元素满足：选择所有满足要求的子集J中具有最大能量者，即选择其次，令Λt＝Λt-1∪Jo，At＝At-1∪aj(j∈Jo)；然后，更新残差则可利用稀疏矩阵可得重构信号：进入ROMP算法后，迭代直到相邻两次重构信号能量差满足时，停止迭代，输出重构信号。本专利技术通过改进StOMP压缩感知重构算法，有效地利用阈值判断标准保证算法收敛速度和精度，同时将相邻两次重构信号的能量差作为算法转化依据和停止迭代的依据，结合了回溯原子选择算法收敛速度快和ROMP算法重构精度高的优点，实现重构精度和收敛速度的折中。附图说明本专利技术的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：图1是本专利技术总体流程图。具体实施方式下面详细描述本专利技术的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本专利技术，而不能理解为对本专利技术的限制。在本专利技术的描述中，需要理解的是，术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底本文档来自技高网...

【技术保护点】
一种基于改进StOMP的压缩感知重构算法，其特征在于，包括：S1，基于分段正交匹配追踪算法，根据阈值大小选择执行回溯原子选择算法或ROMP算法；S2，引入子空间正交匹配追踪算法的回溯思想，即回溯原子选择算法，使每次选择的原子个数等于分段正交匹配追踪算法中选出的列数；S3，将回溯原子选择算法和ROMP嵌入StOMP中，根据相邻两次重构信号能量差判断是否进行算法转换或者迭代停止。

【技术特征摘要】
1.一种基于改进StOMP的压缩感知重构算法，其特征在于，包括：S1，基于分段正交匹配追踪算法，根据阈值大小选择执行回溯原子选择算法或ROMP算法；S2，引入子空间正交匹配追踪算法的回溯思想，即回溯原子选择算法，使每次选择的原子个数等于分段正交匹配追踪算法中选出的列数；S3，将回溯原子选择算法和ROMP嵌入StOMP中，根据相邻两次重构信号能量差判断是否进行算法转换或者迭代停止。2.根据权利要求1所述的一种基于改进StOMP的压缩感知重构算法，其特征在于，所述S1包括：考虑原信号且为可稀疏信号，根据压缩感知理论，将其投影到稀疏基空间上可表示为：X＝ψθ其中，为稀疏系数，含有k个非零元素，则称原信号X是k稀疏的；为信号X的稀疏基；选择随机高斯矩阵作为观测矩阵Φ，其内元素均满足高斯分布N(0,1/n)，则N维信号可由M个观测值表示：Y＝ΦX＝Φψθ＝Aθ其中，是信号的观测值，为传感矩阵；令r0＝Y，并以h＝β×max(th)作为判断标准，确定下一步执行回溯原子选择算法或ROMP算法，其中th＝tsσs，根据经验值取ts＝2.5，1/n为高斯观测矩阵中元素所服从的高斯分布的方差；则根据StOMP，以rt表示残差，t表示迭代次数，表示空集，Jo表示每次迭代找到的列序号，Λt表示t次迭代的列序号集合，aj表示矩阵A的第j列，At表示按索引Λt选出的传感矩阵A的列集合。首先计算u＝<rt-1,aj>,1≤j≤N，选u中大于阈值th＝tsσs的c列，并记对应的A列序号j构成集合Jo；然后令Λt＝Λt-1∪Jo，At＝At-1∪aj(j∈Jo)；再则令ts＝2.5...

【专利技术属性】
技术研发人员：廖勇，谢金豆，周昭坤，陈玲，姜大，
申请(专利权)人：重庆大学，
类型：发明
国别省市：重庆,50