The code that implements the Two Stage Subspace Pursuit + Re-weighted Lp reconstruction algorithm featured on Monday has been cleaned up and is now here. It will be cleaned up further but this version is a good one to start with.
function [x1]=SP_lp_re(K,phi,y,p)
%
% first call a new version of SP (cleaner interface)
%
[xfinal,That]=SPnew(K, phi, y);
%
% second call to re-weighted l_p
%
[x1] = lp_re1(phi,y,p,xfinal);
No comments:
Post a Comment