Optimised_OMP 压缩感知信号恢复算法

一种压缩感知信号恢复算法，针对贪婪迭代类算法中的正交匹配追踪（OMP)算法的改进。OMP在每次迭代过程中选择出的原子并不是最优的，无法使本轮迭代中残差的减少最大化。本例程实现了改进的最优OMP算法，即Optimise

%% Input
%% T: the measuring matrix sized M*N
%% s: the measured signal whose length is M
%% maxCycle: limit of iteration if k is known it should be set as k+1
%% k: sparse of original signal
%% Output
%% hat_y: the estimated value of original signal
%% normR: the L2 of residual in each interation
%% pos_array: selected atom set
function [ hat_y normR pos_array] = OOMP（ TsmaxCyclek ）%
M=size（T1）;
N=size（T2）;
r_n=s;
Aug_t=[];
hat_y=zeros（N1）;
normR=zeros（1maxCycle）;

%Initially set
k=1;
gama=T; % N collumns after orthogonalize
d=ones（N1）;% 分母上
b=zeros（N1）;% 分子上
e=zeros（N1）;% 用于选择原子的标准参数

product=abs（T‘*r_n）;
[varpos]=max（product）;
kesi（:1）=T（:pos）;
pos_array（1）=pos;
Aug_t=[Aug_tT（:pos）];
normR（1）=sqrt（s‘*s-（abs（T（:pos）‘*s））^2）;
normR（1）=normR（1）;
var（

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1604  2012-11-28 10:10  Optimised_OMP.m