function?[Sestd]=cosamp2（PhiuutrueK）

%?Cosamp?algorithm
%???Input
%???????K?:?sparsity?of?Sest
%???????Phi?:?measurement?matrix
%???????u:?measured?vector
%???????tol1?:?tolerance?for?approximation?between?successive?solutions.?
%???Output
%???????Sest:?Solution?found?by?the?algorithm
%???????d???:?success?index?（d=1?is?success?d?=?0?is?no?convergence）
%
%?Algorithm?as?described?in?“CoSaMP:?Iterative?signal?recovery?from?
%?incomplete?and?inaccurate?samples“?by?Deanna?Needell?and?Joel?Tropp.
%?
%?This?implementation?was?written?by?David?Mary
%
%?This?script/program?is?released?under?the?Commons?Creative?Licence
%?with?Attribution?Non-commercial?Share?Alike?（by-nc-sa）
%?http://creativecommons.org/licenses/by-nc-sa/3.0/
%?Short?Disclaimer:?this?script?is?for?educational?purpose?only.
%?Longer?Disclaimer?see??http://igorcarron.googlepages.com/disclaimer

%?Initialization
Sest=zeros（size（utrue））;%不明白不知道utrue是什么
v=u;
t=1;?T2=[];
while?t?[kz]=sort（abs（Phi‘*v））;k=flipud（k）;z=flipud（z）;%
Omega=z（1:2*K）;%取了前2K個。
T=sort（union（OmegaT2））;phit=Phi（:T）;%T為支撐集合并后的集合，PHIT為觀察矩陣phi只取下標為T的列后所組成的矩陣
%keyboard
b=abs（pinv（phit）*u）;%通過最小二乘法進行信號估計，pinv表示偽逆
[k3z3]=sort（（b））;k3=flipud（k3）;z3=flipud（z3）;%把得到的信號的新的估計值進行從大到小的排序
Sest=zeros（size（utrue））;
Sest（T（z3（1:K）））=abs（b（z3（1:K）））;??%信號估計后進行刪減，只取了前K個，并取絕對值后賦給sestsest此時為信號估計
[k2z2]=sort（abs（Sest））;k2=flipud（k2）;z2=flipud（z2）;
T2=z2（1:K）;%T2為信號估計值的最大K個下標集合。
v=u-Phi*Sest;%殘差進行更新
d=0;n2=norm（abs（Sest-utrue）‘inf‘）;
if?n2?<1e-3
????d=1;t=1e10;
????disp（‘CoSaMP:?success‘）;
end
?t=t+1;
?end
?if?d==0
?????disp（‘CoSaMP:?failed‘）
?end
?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1780??2009-12-17?21:01??OMP_SP\cosamp2.m

?????文件???????1369??2008-04-07?06:59??OMP_SP\IRLSregcomp.m

?????文件????????767??2008-04-07?06:59??OMP_SP\krored.m

?????文件???????1105??2013-07-04?10:06??OMP_SP\OMP.m

?????文件???????2142??2008-04-07?06:58??OMP_SP\OMPmod3.m

?????文件???????1256??2009-12-22?11:27??OMP_SP\SPA.m

?????文件???????4932??2013-07-04?09:30??OMP_SP\test.m

?????目錄??????????0??2013-07-04?10:23??OMP_SP

-----------?---------??----------?-----??----

????????????????13351????????????????????8