function [decide_bits] = VBLAST(F,H,in_match,in, P,Ns,NPW);

decide_bits=zeros(P,Ns);
dec_table=(1:1:P).';  % a index table helps make decision in tx_order
for k=1:P
   G=inv(F*F' + NPW*F)*F;
   [ordered_norm,index]=sort(diag(G'*G));
   soft_symbol=G(:,index(1))'*in_match;
   decide = 1/sqrt(2)*(sign(real(soft_symbol)) + j*sign(imag(soft_symbol)));
   decide_bits(dec_table(index(1)),:) = decide;
   in=in-H(:,index(1))*decide_bits(dec_table(index(1)),:);
   dec_table(index(1))=[];
   H(:,index(1))=[];
   in_match=H'*in;
   F=H'*H;
end