OMP算法代码学习

正交匹配追踪(OMP)算法的MATLAB函数代码并给出单次测试例程代码

测量数M与重构成功概率关系曲线绘制例程代码

信号稀疏度K与重构成功概率关系曲线绘制例程代码

 

参考来源：http://blog.csdn.net/jbb0523/article/details/45130793

参考文献：Joel A. Tropp and Anna C. Gilbert. Signal Recovery From Random Measurements Via Orthogonal Matching Pursuit[J]. IEEETransactions on Information Theory, VOL. 53, NO. 12, DECEMBER 2007.

---

0、符号说明如下

　　 压缩观测y=Φx，其中y为观测所得向量M×1，x为原信号N×1（M<<N）。x一般不是稀疏的，但在某个变换域Ψ是稀疏的，即x=Ψθ，其中θ为K稀疏的，即θ只有K个非零项。此时y=ΦΨθ，令A=ΦΨ，则y=Aθ。

        (1)y为观测所得向量，大小为M×1

 

        (2)x为原信号，大小为N×1

 

        (3)θ为K稀疏的，是信号在x在某变换域的稀疏表示

 

        (4)Φ称为观测矩阵、测量矩阵、测量基，大小为M×N

 

        (5)Ψ称为变换矩阵、变换基、稀疏矩阵、稀疏基、正交基字典矩阵，大小为N×N

 

        (6)A称为测度矩阵、传感矩阵、CS信息算子，大小为M×N

　　 上式中，一般有K<<M<<N，后面三个矩阵各个文献的叫法不一，以后我将Φ称为测量矩阵、将Ψ称为稀疏矩阵、将A称为传感矩阵。

1、OMP重构算法流程

2、正交匹配追踪(OMP)MATLAB代码(CS_OMP.m)

function[theta]=CS_OMP(y,A,t)
%CS_OMP Summary of this function goes here
 %Version: 1.