【误码率仿真】基于matlab 8PSK调制解调误码率仿真【含Matlab源码 3837期】.zip

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);