%???感謝親親使用此代碼，此代碼解決您的問題了嗎~（@^_^@）~
%???沒解決的話告訴親親一個好消息，登錄淘寶店鋪“大成軟件工作室”，可以下載（????）1分錢成品代碼（′▽‘〃）哦~
%???是的，親親真的沒有看錯，撓破頭皮的問題真的1分錢就可以解決了（づ??????）づ
%???小的這就把傳送門給您，記得要收藏好哦（づ￣3￣）づ╭?～
%???傳送門：https://item.taobao.com/item.htm?spm=a1z10.1-c.w4004-15151018122.5.uwGoq5&id=538759553146
%???如果傳送門失效，親親可以來店鋪討要，客服MM等親親來騷擾哦~（*/ω╲*）
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只取下標(biāo)為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個下標(biāo)集合。
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
?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2016-09-22?16:54??MMA14\
?????文件????????1787??2009-12-17?21:00??MMA14\cosamp2.asv
?????文件????????2254??2016-09-22?15:47??MMA14\cosamp2.m
?????文件????????1346??2016-09-22?15:47??MMA14\disclaimer.m
?????文件????????1837??2016-09-22?15:47??MMA14\IRLSregcomp.m
?????文件????????1235??2016-09-22?15:47??MMA14\krored.m
?????文件????????1548??2016-09-22?15:47??MMA14\OMP.m
?????文件????????2610??2016-09-22?15:48??MMA14\OMPmod3.m
?????文件????????4920??2009-12-22?11:22??MMA14\scripttf2dcosamp.asv
?????文件????????5394??2016-09-22?15:49??MMA14\scripttf2dcosamp.m
?????文件????????1256??2009-12-22?11:22??MMA14\SPA.asv
?????文件????????1730??2016-09-22?15:49??MMA14\SPA.m