cm_sm32.m

基于MATLAB实现信号的调制解调源代码的简单事例.

function[pb,ps]=cm_sm32(snr_in_dB)
N=10000;
E=1;
snr_in_dB=5;
snr=10^(snr_in_dB);
sgma=sqrt(E/snr)/2;
%信号映射，采用pi/M系统
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)=1;
    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
    [n(1),n(2)]=gngauss(sgma);              %产生高斯噪声
    %相干解调
    if((dsource1(i))==0)&(dsource2(i)==0),
        r=s00+i;
    elseif((dsource1(i))==0)&(dsource2(i)==1),
        r=s01+i;
    elseif((dsource1(i))==1)&(dsource2(i)==0),
        r=s10+i;
    else
        r=s11+i;
    end;
    %判决
    c00=dot(r,s00);
    c01=dot(r,s01);
    c10=dot(r,s10);
    c11=dot(r,s11);
    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;
    else
        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;
    end;
end;
%误码结果
ps=numofsymbolerror/N;
pb=numofbiterror/(2*N);