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个比特