cl_complex.asv

v-blast4x4_dast程序

function [u,itera_times]=CL_complex(r_in,M_in,SS,C)
% 解码输入，在实数范围内
% 输入参数：
% r_in为接收到的信号
% M_in为格生成矩阵
% SS为星座集
r=[real(r_in);imag(r_in)];
M1=real(M_in);
M2=imag(M_in);
M=[M1,-M2;M2,M1];
m=size(M,1);
%G=M'*M;
%R=chol(G);   
[Q,R]=qr(M);
maxlen=length(SS);
p=pinv(M)*r;

for ii=1:size(R,1)
    for jj=1:size(R,2)
        if ii==jj
            qq(ii,jj)=R(ii,jj)^2;
        elseif ii<jj
            %qq(ii,jj)=R(ii,jj)*pinv(R(ii,ii));
            qq(ii,jj)=R(ii,jj)/(R(ii,ii));
        else
            qq(ii,jj)=0;
        end;
    end;
end;

d2=C;
T(m)=C;
S(m)=p(m);
ii=m;
itera_times=0;
aa = qq;
while(1>0)
    if qq(ii,ii) == inf
        UU(ii) = Q_latter(S(ii),SS);
        L(ii) = Q_former(S(ii),SS);
    else
        if qq(ii,ii) == 0
            UU(ii)=Q_latter(inf,SS);
            L(ii)=Q_former(-inf,SS);
        else
            UU(ii)=Q_latter(sqrt(T(ii)/qq(ii,ii))+S(ii),SS);%取值上界
            L(ii)=Q_former(-sqrt(T(ii)/qq(ii,ii))+S(ii),SS);%取值下界
        end;
    end;
    N(ii)=length(sorted_enum(L(ii),UU(ii),SS,S(ii)));
    n=maxlen-N(ii);
   extra=linspace(inf,inf,n+1);
    z(ii,:)=[sorted_enum(L(ii),UU(ii),SS,S(ii)) extra];
    x(ii)=0;
     while(1>0)
        while(1>0)
            itera_times=itera_times+1;
            x(ii)=x(ii)+1;
            if(x(ii)>N(ii))
                if (ii==m)
                    u=reshape(u,m/2,2);
                    u=u(:,1)+j.*u(:,2);
                    u=reshape(u,1,size(u,1));%行向量输出
                    dd2=d2;
                    return;
                else
                    ii=ii+1;
                end;
            else
                break;
            end;
        end; %最里层的while
        if (ii>1)
            break;
        else
            d_2=T(m)-T(1)+qq(1,1)*(S(1)-z(1,x(ii)))^2; 
        end;
        if d_2<d2
            d2=d_2;
            T(m)=d_2;
            k=m;
            while(k>=2)
                T(k-1)=T(k)-qq(k,k)*(S(k)-z(k,x(k)))^2;
                k=k-1;
            end;
            for k=1:m
                u(k)=z(k,x(k));%%得到有效格点各维坐标值
                if qq(k,k) == inf
                   UU(k) = Q_latter(S(ii),SS);
                   L(k) = Q_former(S(ii),SS);
                else
                  if qq(k,k) == 0
                     UU(k)=Q_latter(inf,SS);
                     L(k)=Q_former(-inf,SS);
                  else
                    UU(k)=Q_latter(sqrt(T(k)/qq(k,k))+S(k),SS);%取值上界
             L(ii)=Q_former(-sqrt(T(ii)/qq(ii,ii))+S(ii),SS);%取值下界
        end;
    end;
                UU(k)=Q_latter(sqrt(T(k)/qq(k,k))+S(k),SS);%取值上界
                L(k)=Q_former(-sqrt(T(k)/qq(k,k))+S(k),SS);%取值下界
                N(k)=length(Enum(L(k),UU(k),SS));
            end;
        else
            continue;
        end;
    end;  % 第二个while
        E(ii)=p(ii)-z(ii,x(ii));%%  ??  y(ii,x(ii))    if(ii>1)  break;
        T(ii-1)=T(ii)-qq(ii,ii)*(S(ii)-z(ii,x(ii)))^2;%%
        sum=0.0;
        for jj=ii:m
            sum=sum+qq(ii-1,jj)*E(jj);
        end;
        S(ii-1)=p(ii-1)+sum;
        ii=ii-1;
end;% 最外层while