function?hat_x=cs_omp（yT_Matm）
%?y=T_Mat*x?T_Mat?is?n-by-m
%?y?-?measurements
%?T_Mat?-?combination?of?random?matrix?and?sparse?representation?basis
%?m?-?size?of?the?original?signal
%?the?sparsity?is?length（y）/4

n=length（y）;
s=floor（n/4）;?????????????????????????????????????%??測量值維數(shù)
hat_x=zeros（1m）;?????????????????????????????????%??待重構(gòu)的譜域（變換域）向量?????????????????????
Aug_t=[];?????????????????????????????????????????%??增量矩陣（初始值為空矩陣）
r_n=y;????????????????????????????????????????????%??殘差值?

for?times=1:s;??????????????????????????????????%??迭代次數(shù)（稀疏度是測量的1/4）

????product=abs（T_Mat‘*r_n）;

????[valpos]=max（product）;???????????????????????%??最大投影系數(shù)對應(yīng)的位置
????Aug_t=[Aug_tT_Mat（:pos）];???????????????????%??矩陣擴充
????T_Mat（:pos）=zeros（n1）;??????????????????????%??選中的列置零（實質(zhì)上應(yīng)該去掉，為了簡單將其置零）
????aug_x=（Aug_t‘*Aug_t）^（-1）*Aug_t‘*y;???????????%??最小二乘使殘差最小
????r_n=y-Aug_t*aug_x;????????????????????????????%??殘差
????pos_array（times）=pos;?????????????????????????%??紀(jì)錄最大投影系數(shù)的位置

end
hat_x（pos_array）=aug_x;???????????????????????????%??重構(gòu)的向量?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-10-21?16:25??OMP\
?????文件??????262308??2015-10-20?09:52??OMP\barbara.pgm
?????文件????????1186??2015-10-20?09:52??OMP\cs_omp.m
?????文件?????????164??2015-10-20?09:52??OMP\dct_basis.m
?????文件????????1931??2015-10-21?19:56??OMP\Demo_CS_OMP.m
?????文件?????????599??2015-10-20?09:52??OMP\GenerateProjection.m
?????文件??????262299??2015-10-20?09:52??OMP\goldhill.pgm
?????文件??????262302??2015-10-20?09:52??OMP\lenna.pgm
?????文件??????262218??2015-10-20?09:52??OMP\mandrill.pgm
?????文件?????????857??2015-10-20?09:52??OMP\omp.m
?????文件??????262230??2015-10-20?09:52??OMP\peppers.pgm
?????文件?????2014489??2015-10-20?09:52??OMP\projections.1.mat
?????文件????????1329??2015-10-20?09:52??OMP\PSNR.m
?????文件????????1021??2015-10-20?09:52??OMP\test_omp.m
?????文件????????1473??2015-10-20?09:52??OMP\Untitled3.m