Hi everybody,

I am currently working in "low-rankness channel matrix as a sparsity approach to use compressive sensing for channel estimation in massive MIMO systems.", and I am trying to recreate a model and perform simulations related to this paper:

"Vanidevi, M., & Selvaganesan, N. (2017). Channel Estimation for Finite Scatterers Massive Multi-User MIMO System. Circuits, Systems, and Signal Processing, 1-17."

In fact, I developed a code  attached below to recreate figure 4 in this paper.

Concerning the algorithm(singular value thresholding+landweber) solving the nuclear norm minimization I developped this matlab code:

**************************************************

function [H_optimal,numiter_NN]=NN_estimator(lambda,Phi_matrix,y_rec_per_antenna_mat_with_noise,beta,M,K,L)

%L = 10;%pilot length

%M=antennas number

%initialisation

H=zeros(M,K);

psi_matrix=kron(Phi_matrix',eye(M));

maxiter = 1000;

dHmax=1e-2;

niter = 0;

dH= 1000000;

thr=lambda; %thresholding_coef

%(niter