cm_sm32.m

本课程提供的示范程序包括三个方面

function[pb,ps]=cm_sm32(snr_in_dB)

N=10000;
E=1;
snr=10^(snr_in_dB/10);
sgma=sqrt(E/(3*2*snr));

s000=[1 0];
s001=[cos(pi/4) sin(pi/4)];
s011=[0 1];
s010=[cos(3*pi/4) sin(3*pi/4)];
s110=[-1 0];
s111=[cos(5*pi/4) sin(5*pi/4)];
s101=[0 -1];
s100=[cos(7*pi/4) sin(7*pi/4)];

for i=1:N,
   temp=rand;
   
   if(temp<0.125),
      dsource1(i)=0;
      dsource2(i)=0;
      dsource3(i)=0;
   elseif(temp<0.25),
      dsource1(i)=0;
      dsource2(i)=0;
      dsource3(i)=1;
   elseif(temp<0.375),
      dsource1(i)=0;
      dsource2(i)=1;
      dsource3(i)=0;
   elseif(temp<0.5),
      dsource1(i)=0;
      dsource2(i)=1;
      dsource3(i)=1;
   elseif(temp<0.625),      
      dsource1(i)=1;
      dsource2(i)=0;
      dsource3(i)=0;
   elseif(temp<0.75),
      dsource1(i)=1;
      dsource2(i)=0;
      dsource3(i)=1;
   elseif(temp<0.875),
      dsource1(i)=1;
      dsource2(i)=1;
      dsource3(i)=0;
   else
      dsource1(i)=1;
      dsource2(i)=1;
      dsource3(i)=1;
   end;
end;

numofsymbolerror=0;
numofbiterror=0;

for i=1:N,
   n(1)=gngauss(sgma);
   n(2)=gngauss(sgma);
   
   if((dsource1(i)==0)&(dsource2(i)==0)&(dsource3(i)==0)),
      r=s000+n;
   elseif((dsource1(i)==0)&(dsource2(i)==0)&(dsource3(i)==1)),
      r=s001+n;
   elseif((dsource1(i)==0)&(dsource2(i)==1)&(dsource3(i)==0)),
      r=s010+n;
   elseif((dsource1(i)==0)&(dsource2(i)==1)&(dsource3(i)==1)),
      r=s011+n;
   elseif((dsource1(i)==1)&(dsource2(i)==0)&(dsource3(i)==0)),
      r=s100+n;
   elseif((dsource1(i)==1)&(dsource2(i)==0)&(dsource3(i)==1)),
      r=s101+n;
   elseif((dsource1(i)==1)&(dsource2(i)==1)&(dsource3(i)==0)),
      r=s110+n;
   else
      r=s111+n;
   end;
   
   c000=dot(r,s000);
   c001=dot(r,s001);
   c010=dot(r,s010);
   c011=dot(r,s011);
   c100=dot(r,s100);
   c101=dot(r,s101);
   c110=dot(r,s110);
   c111=dot(r,s111);
   
   c_max=max([c000 c001 c010 c011 c100 c101 c110 c111]);
   
   if(c000==c_max),
      decis1=0;decis2=0;decis3=0;
   elseif(c001==c_max),
      decis1=0;decis2=0;decis3=1;
   elseif(c010==c_max),
      decis1=0;decis2=1;decis3=0;
   elseif(c011==c_max),
      decis1=0;decis2=1;decis3=1;
   elseif(c100==c_max),
      decis1=1;decis2=0;decis3=0;
   elseif(c101==c_max),
      decis1=1;decis2=0;decis3=1;
   elseif(c110==c_max),
      decis1=1;decis2=1;decis3=0;
   else
      decis1=1;decis2=1;decis3=1;
   end;
   
   symbolerror=0;
   if(decis1~=dsource1(i)),
      numofbiterror=numofbiterror+1;
      symbolerror=1;
   end;
   
   if(decis2~=dsource2(i)),
      numofbiterror=numofbiterror+1;
      symbolerror=1;
   end;
   
   if(decis3~=dsource3(i)),
      numofbiterror=numofbiterror+1;
      symbolerror=1;
   end;
   
   if(symbolerror==1),
      numofsymbolerror=numofsymbolerror+1;
   end;
end;

ps=numofsymbolerror/N;
pb=numofbiterror/(3*N);