sim_select_qstbc_ml_n1.txt

QPSK误码率仿真程序

snrindB=0;
snr=10^(snrindB/10);
M=4;
m=log2(M);
sgma=sqrt(0.5);
Eb=snr*sgma^2; 
Es=m*Eb;
es=sqrt(Es);
c=2*pi/M;
s00=es*exp(j*0*c);
s01=es*exp(j*1*c);
s10=es*exp(j*2*c);
s11=es*exp(j*3*c);

SER=0;
for i=1:20000
         temp1=4*rand;temp2=4*rand;temp3=4*rand;temp4=4*rand;
         if(temp1<1)
         s1=s00;
         elseif(temp1<2)
         s1=s01;
         elseif(temp1<3)
         s1=s10;
         else
         s1=s11;
         end
         if(temp2<1)
         s2=s00;
         elseif(temp2<2)
         s2=s01;
         elseif(temp2<3)
         s2=s10;
         else
         s2=s11;
         end
         if(temp3<1)
         s3=s00;
         elseif(temp3<2)
         s3=s01;
         elseif(temp3<3)
         s3=s10;
         else
         s3=s11;
         end
         if(temp4<1)
         s4=s00;
         elseif(temp4<2)
         s4=s01;
         elseif(temp4<3)
         s4=s10;
         else
         s4=s11;
         end
     for jj=1:4
         h(jj)=rayleigh(sqrt(0.5))*exp(j*2*pi*rand);
     end
    h1=h(1);h2=h(2);h3=h(3);h4=h(4);
     
         

      n1=gngauss;n2=gngauss;n3=gngauss;n4=gngauss;
      r1=h1*s1+h2*s2+h3*s3+h4*s4+n1;
      r2=-h1*conj(s2)+h2*conj(s1)-h3*conj(s4)+h4*conj(s3)+n2;
      r3=-h1*conj(s3)-h2*conj(s4)+h3*conj(s1)+h4*conj(s2)+n3;
      r4=h1*s4-h2*s3-h3*s2+h4*s1+n4;

      for m1=1:4
           for m2=1:4
              S1=es*exp(j*(m1-1)*c);
              S4=es*exp(j*(m2-1)*c);
              e14(m1,m2)=f14(h1,h2,h3,h4,S1,S4,r1,r2,r3,r4);
            end
      end
      mine14=min(min(e14));
      for m1=1:4
           for m2=1:4
              S1=es*exp(j*(m1-1)*c);
              S4=es*exp(j*(m2-1)*c);
              if f14(h1,h2,h3,h4,S1,S4,r1,r2,r3,r4)==mine14
              s11=S1;s14=S4;
              end
            end
      end
      
      if s11~=s1
          SER=SER+1;
      end
      if s14~=s4
          SER=SER+1;
      end

      for m1=1:4
           for m2=1:4
              S2=es*exp(j*(m1-1)*c);
              S3=es*exp(j*(m2-1)*c);
              e23(m1,m2)=f23(h1,h2,h3,h4,S2,S3,r1,r2,r3,r4);
            end
      end
      mine23=min(min(e23));
      for m1=1:4
           for m2=1:4
              S2=es*exp(j*(m1-1)*c);
              S3=es*exp(j*(m2-1)*c);
              if f23(h1,h2,h3,h4,S2,S3,r1,r2,r3,r4)==mine23
              s12=S2;s13=S3;
              end
            end
      end      
              
      if s12~=s2
          SER=SER+1;
      end
      if s13~=s3
          SER=SER+1; 
      end
end
      Pser=SER/80000