t2.txt

8psk的源程序

8psk的源程序： 
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);