pskmoto.m

qpsk通信系统的MonteCalo仿真程序

function [pb,ps]=pskmoto(snr_in_dB)
%pskmoto求出以db为单位的给定信噪比的比特误码率和符号误码率
N=1000;%符号数
E=1;
snr=10^(snr_in_dB/10);         %计算信噪比的数值
sgma=sqrt(E/snr)/2;
%信号映射
n=[0 0];
s00=[1  0];   s01=[0  1];   s11=[-1  0];  s10=[0  -1];
for  i=1:N
     temp=rand
     if(temp<0.25)
         dsource1(i)=0;
         dsource2(i)=0;
     elseif(temp<0.5)
         dsource1(i)=0;
         dsource2(i)=0;
     elseif(temp<0.75)
         dsource1(i)=1;
         dsource2(i)=0;
     else
         dsource1(i)=1;
         dsource2(i)=1;
     end;
end;
%判决、误码率的计算
numofsymbolerror=0;   numofbiterror=0;
for i=1:N
    %在判决器的接收端的信号，对于第i个符号为：
   n=gngauss(sgma);
   if((dsource1(i)==0)&(dsource2(i)==0))
       r=s00+n;
   elseif((dsource1(i)==0)&(dsource2(i)==1))
       r=s01+n;
   elseif((dsource1(i)==1)&(dsource2(i)==0))
       r=s10+n;
   else
       r=s11+n;
   end
   %以下计算互相关量度
 c00=dot(r,s00);
 c01=dot(r,s01);
 c10=dot(r,s10);
 c11=dot(r,s11);
 %第i个符号的判决如下进行
 c_max=max([c00  c01  c10  c11]);
 if(c00==c_max)
    decis1=0;  decis2=0;
 elseif(c01==c_max)
     decis1=0; decis2=1;
 elseif(c10==c_max)
     decis1=1; decis2=0;
 elseif(c11==c_max)
     decis1=1; decis2=1;
 end;
 %若判决不对，记错器数加1
 symbolerror=0;
 if(decis1~=dsource1(i))
     numofbiterror=numofbiterror+1;  symbolerror=1;
 end;
  if(decis2~=dsource2(i))
     numofbiterror=numofbiterror+1;  symbolerror=1;
 end;
 if(symbolerror==1)
     numofsymbolerror=numofsymbolerror+1;%误符号数加1
 end;
end;
ps=numofsymbolerror/N;        %总共发出N个符号
pb=numofbiterror/(2*N);       %总共发出2N个比特