压缩感知中OMP算法的C/C++实现

压缩感知中OMP算法的C/C++实现

• 背景介绍
• 算法实现部分
• 总结

---

阅读之前注意：

本文阅读建议用时：30min
本文阅读结构如下表：

项目	下属项目	测试用例数量
背景介绍	无	0
算法实现部分	无	1
总结	无	0

背景介绍

国内关于压缩感知的研究正在迅速发展，这里简要介绍一下压缩感知：所谓的压缩感知，就是把采样和压缩结合起来同时进行，这突破了传统的奈奎斯特采样定理（也叫香农定理），是信号采集理论的一次重大突破。压缩感知对于信号采集端的要求很低，并且只需要采集少量的信息，但采集信息后，在后台进行信号恢复时，需要进行大量的运算，因此信号重构算法和重构算法的GPU加速也是一大研究方向。
更多关于压缩感知，以及压缩感知的信号恢复重构算法，可以参考这两个博主：彬彬有礼的压缩感知专栏和Rachel Zhang的专栏，文章都写得很棒。所谓见贤思齐焉，我还得多多向他们学习！

算法实现部分

这一部分的算法实现（C/C++），改写自博客中的Matlab程序。
以下代码在配置Eigen库后，即可直接运行。

//本程序需要配置Eigen库
//OMP算法参考博客：https://blog.csdn.net/jbb0523/article/details/45130793
//程序改写自博客中的Matlab程序
//OMP_Version_1.0 2018.03.28_16:35 by Cooper Liu
//Questions ? Contact me : angelpoint@foxmail.com
#include"iostream"
#include <ctime>  
#include <Eigen/Eigenvalues>
using namespace std;
using namespace Eigen;

void swap(int *a, int i, int j)
{
   
	int tmp = a[i];
	a[i] = a[j];
	a[j] = tmp;
}

void randperm(int *a,int N)
{
   
	int i = 0;
	for (i = 0; i < N; i++)
		a[i