vblast.m

v-blast4x4_dast程序

function u=VBLAST(r_in,M_in,m,s)
%产生接收信号和信道矩阵；
y=[real(r_in);imag(r_in)];
M1=real(M_in);
M2=imag(M_in);
M=[M1,-M2;M2,M1];
n=2*m;
e=eye(n);
p=zeros(n,1);
H=M;
%u = zeros(1,4);
%G=inv(H);
G = pinv(H);
   for i=1:n
       x=G*y;
       k(i)=argmin(G);
       u(k(i))=zuijing(x(k(i)),s);
       e(k(i),k(i))=0;
       p(k(i))=1;
       y=y-H*p*u(k(i));
       p(k(i))=0;
       H=H*e;
       G=pinv(H);     
   end;
   %for i=1:n
    %   for p=1:n
     %      if k(p)==i
      %         v(i)=u(p);
       %        break;
       %   end;
       %end;
       %end;
   u=reshape(u,m,2);
   u=u(:,1)+j.*u(:,2);
   u=reshape(u,1,size(u,1));