MMA14.zip

MMA14.zip

%   感谢亲亲使用此代码，此代码解决您的问题了吗~（@^_^@）~
%   没解决的话告诉亲亲一个好消息，登录淘宝店铺“大成软件工作室”，可以下载（????）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 < 101 
[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
 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           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