diff --git a/VanillaID/tac_reconstruction.m b/VanillaID/tac_reconstruction.m index dce3b21..ca6577a 100644 --- a/VanillaID/tac_reconstruction.m +++ b/VanillaID/tac_reconstruction.m @@ -36,7 +36,7 @@ UU=zeros(N, MAXITER); w_estimate=zeros(N, MAXITER); WWW=ones(N, MAXITER); -fprintf(1, 'Finding a sparse feasible point using l1-norm heuristic ...'); +fprintf(1, 'Finding a sparse feasible point using l1-norm heuristic ...\n'); for iter=1:1:MAXITER @@ -52,7 +52,7 @@ variable W(N) w_estimate(:,iter)=W; WWW(:,iter)=W; - Gamma(:,iter)=U(:,iter).^-1.*W; + Gamma(:,iter)=U(:,iter).^-1.*abs(W); Dic0=lambda*eye(M)+Dic*diag(Gamma(:,iter))*Dic'; UU(:, iter)=diag(Dic'*(Dic0\Dic)); U(:,iter+1)=abs(sqrt(UU(:, iter))); @@ -64,4 +64,3 @@ variable W(N) end end -