function [pb,ps]=pskmoto(snr_in_dB)
%pskmoto求出以db為單位的給定信噪比的比特誤碼率和符號誤碼率
N=1000;%符號數(shù)
E=1;
snr=10^(snr_in_dB/10);         %計算信噪比的數(shù)值
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;
 %若判決不對，記錯器數(shù)加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;%誤符號數(shù)加1
 end;
end;
ps=numofsymbolerror/N;        %總共發(fā)出N個符號
pb=numofbiterror/(2*N);       %總共發(fā)出2N個比特