decode.m

ofdm中互补序列降低par的ber程序

function dec=decode(in)
n=length(in);
H=[1 1 1 2;3 3 0 3;3 0 3 3;1 0 0 0;0 3 3 3;0 1 0 0;0 0 1 0;0 0 0 1];
b=[0 0 0 2 0 0 2 0];
dec=zeros(1,n/2);
for j=1:n/8
    tt=in(j*8-7:j*8)-b;
    tt=mod(tt,4);
    dd=tt*H;
    dd=mod(dd,4);
    if dd(1)==1&&dd(2)==1&&dd(3)==1&&dd(4)==2
        tt(1)=mod(tt(1)-1,4);
    elseif dd(1)==2&&dd(2)==2&&dd(3)==2&&dd(4)==0
        tt(1)=mod(tt(1)-2,4);
    elseif dd(1)==3&&dd(2)==3&&dd(3)==3&&dd(4)==2
        tt(1)=mod(tt(1)-3,4);
    elseif dd(1)==3&&dd(2)==3&&dd(3)==0&&dd(4)==3
        tt(2)=mod(tt(2)-1,4);
    elseif dd(1)==2&&dd(2)==2&&dd(3)==0&&dd(4)==2
        tt(2)=mod(tt(2)-2,4);
    elseif dd(1)==1&&dd(2)==1&&dd(3)==0&&dd(4)==1
        tt(2)=mod(tt(2)-3,4);
    elseif dd(1)==3&&dd(2)==0&&dd(3)==3&&dd(4)==3
        tt(3)=mod(tt(3)-1,4);
    elseif dd(1)==2&&dd(2)==0&&dd(3)==2&&dd(4)==2
        tt(3)=mod(tt(3)-2,4);
    elseif dd(1)==1&&dd(2)==0&&dd(3)==1&&dd(4)==1
        tt(3)=mod(tt(3)-3,4);
    elseif dd(1)==0&&dd(2)==3&&dd(3)==3&&dd(4)==3
        tt(5)=mod(tt(5)-1,4);
    elseif dd(1)==0&&dd(2)==2&&dd(3)==2&&dd(4)==2
        tt(5)=mod(tt(5)-2,4);
    elseif dd(1)==0&&dd(2)==1&&dd(3)==1&&dd(4)==1
        tt(5)=mod(tt(5)-3,4);
    end
    dec(4*j-3)=tt(1);
    dec(4*j-2)=tt(2);
    dec(4*j-1)=tt(3);
    dec(4*j)=tt(5);
end