quasidecode88.m

包括八天线准正交空时码的编码和解码函数

function  out=quasidecode88(H,Nr,r,n,M)

m1=[];m2=[];m3=[];g1=[];g2=[];g3=[];out=[];F=[];
 s=qammod(0:M-1,M);

     if n==2
     for t=1:Nr
        m1=[conj(H(t,2)),-conj(H(t,1)),-conj(H(t,4)),conj(H(t,3)),-conj(H(t,6)),conj(H(t,5)),conj(H(t,8)),-conj(H(t,7))];
        g1(2*t-1:2*t,:)=[H(t,:);m1];
     end
     r11=[r(:,1) conj(r(:,2))];
     r_1=reshape(r11.',[],1);
     elseif n==3
        for t=1:Nr 
        m1=[conj(H(t,2)),-conj(H(t,1)),-conj(H(t,4)),conj(H(t,3)),-conj(H(t,6)),conj(H(t,5)),conj(H(t,8)),-conj(H(t,7))];
        m2=[conj(H(t,3)),conj(H(t,4)),-conj(H(t,1)),-conj(H(t,2)),-conj(H(t,7)),-conj(H(t,8)),conj(H(t,5)),conj(H(t,6))];
        g1(3*t-2:3*t,:)=[H(t,:);m1;m2];
        end
    r11=[r(:,1) conj(r(:,2)) conj(r(:,3))];
    r_1=reshape(r11.',[],1);  
     elseif n==4
         for t=1:Nr
        m1=[conj(H(t,2)),-conj(H(t,1)),-conj(H(t,4)),conj(H(t,3)),-conj(H(t,6)),conj(H(t,5)),conj(H(t,8)),-conj(H(t,7))];
        m2=[conj(H(t,3)),conj(H(t,4)),-conj(H(t,1)),-conj(H(t,2)),-conj(H(t,7)),-conj(H(t,8)),conj(H(t,5)),conj(H(t,6))];
       m3=[H(t,4),-H(t,3),H(t,2),-H(t,1),-H(t,8),H(t,7),-H(t,6),H(t,5)];
       g1(4*t-3:4*t,:)=[H(t,:);m1;m2;m3];
         end
    r11=[r(:,1) conj(r(:,2)) conj(r(:,3)) r(:,4)];
    r_1=reshape(r11.',[],1);
     elseif n==5
         for t=1:Nr
       m1=[conj(H(t,2)),-conj(H(t,1)),-conj(H(t,4)),conj(H(t,3)),-conj(H(t,6)),conj(H(t,5)),conj(H(t,8)),-conj(H(t,7))];
       m2=[conj(H(t,3)),conj(H(t,4)),-conj(H(t,1)),-conj(H(t,2)),-conj(H(t,7)),-conj(H(t,8)),conj(H(t,5)),conj(H(t,6))];
       m3=[H(t,4),-H(t,3),H(t,2),-H(t,1),-H(t,8),H(t,7),-H(t,6),H(t,5)];
       m4=[conj(H(t,5)),conj(H(t,6)),conj(H(t,7)),conj(H(t,8)),-conj(H(t,1)),-conj(H(t,2)),-conj(H(t,3)),-conj(H(t,4))];
       g1(5*t-4:5*t,:)=[H(t,:);m1;m2;m3;m4];  
         end
         r11=[r(:,1) conj(r(:,2)) conj(r(:,3)) r(:,4) conj(r(:,5))];
    r_1=reshape(r11.',[],1);
     elseif n==6
         for t=1:Nr
       m1=[conj(H(t,2)),-conj(H(t,1)),-conj(H(t,4)),conj(H(t,3)),-conj(H(t,6)),conj(H(t,5)),conj(H(t,8)),-conj(H(t,7))];
       m2=[conj(H(t,3)),conj(H(t,4)),-conj(H(t,1)),-conj(H(t,2)),-conj(H(t,7)),-conj(H(t,8)),conj(H(t,5)),conj(H(t,6))];
       m3=[H(t,4),-H(t,3),H(t,2),-H(t,1),-H(t,8),H(t,7),-H(t,6),H(t,5)];
       m4=[conj(H(t,5)),conj(H(t,6)),conj(H(t,7)),conj(H(t,8)),-conj(H(t,1)),-conj(H(t,2)),-conj(H(t,3)),-conj(H(t,4))];
       m5=[H(t,6),-H(t,5),H(t,8),-H(t,7),H(t,2),-H(t,1),H(t,4),-H(t,3)];
       g1(6*t-5:6*t,:)=[H(t,:);m1;m2;m3;m4;m5];  
         end
         r11=[r(:,1) conj(r(:,2)) conj(r(:,3)) r(:,4) conj(r(:,5)) r(:,6)];
    r_1=reshape(r11.',[],1); 
     elseif n==7
        for t=1:Nr
       m1=[conj(H(t,2)),-conj(H(t,1)),-conj(H(t,4)),conj(H(t,3)),-conj(H(t,6)),conj(H(t,5)),conj(H(t,8)),-conj(H(t,7))];
       m2=[conj(H(t,3)),conj(H(t,4)),-conj(H(t,1)),-conj(H(t,2)),-conj(H(t,7)),-conj(H(t,8)),conj(H(t,5)),conj(H(t,6))];
       m3=[H(t,4),-H(t,3),H(t,2),-H(t,1),-H(t,8),H(t,7),-H(t,6),H(t,5)];
       m4=[conj(H(t,5)),conj(H(t,6)),conj(H(t,7)),conj(H(t,8)),-conj(H(t,1)),-conj(H(t,2)),-conj(H(t,3)),-conj(H(t,4))];
       m5=[H(t,6),-H(t,5),H(t,8),-H(t,7),H(t,2),-H(t,1),H(t,4),-H(t,3)];
       m6=[H(t,7),-H(t,8),-H(t,5),H(t,6),H(t,3),-H(t,4),-H(t,1),H(t,2)];
       g1(7*t-6:7*t,:)=[H(t,:);m1;m2;m3;m4;m5;m6];  
         end
    r11=[r(:,1) conj(r(:,2)) conj(r(:,3)) r(:,4) conj(r(:,5)) r(:,6) r(:,7)];
    r_1=reshape(r11.',[],1); 
     else
       for t=1:Nr
       m1=[conj(H(t,2)),-conj(H(t,1)),-conj(H(t,4)),conj(H(t,3)),-conj(H(t,6)),conj(H(t,5)),conj(H(t,8)),-conj(H(t,7))];
       m2=[conj(H(t,3)),conj(H(t,4)),-conj(H(t,1)),-conj(H(t,2)),-conj(H(t,7)),-conj(H(t,8)),conj(H(t,5)),conj(H(t,6))];
       m3=[H(t,4),-H(t,3),H(t,2),-H(t,1),-H(t,8),H(t,7),-H(t,6),H(t,5)];
       m4=[conj(H(t,5)),conj(H(t,6)),conj(H(t,7)),conj(H(t,8)),-conj(H(t,1)),-conj(H(t,2)),-conj(H(t,3)),-conj(H(t,4))];
       m5=[H(t,6),-H(t,5),H(t,8),-H(t,7),H(t,2),-H(t,1),H(t,4),-H(t,3)];
       m6=[H(t,7),-H(t,8),-H(t,5),H(t,6),H(t,3),-H(t,4),-H(t,1),H(t,2)];
       m7=[conj(H(t,8)),conj(H(t,7)),-conj(H(t,6)),-conj(H(t,5)),conj(H(t,4)),conj(H(t,3)),-conj(H(t,2)),-conj(H(t,1))];
       g1(8*t-7:8*t,:)=[H(t,:);m1;m2;m3;m4;m5;m6;m7];  
         end
    r11=[r(:,1) conj(r(:,2)) conj(r(:,3)) r(:,4) conj(r(:,5)) r(:,6) r(:,7) conj(r(:,8))];
    r_1=reshape(r11.',[],1);   
     end
  out1=pinv(g1)*r_1;

  out=qamdemod(out1,M);