cm_sm34.m

通原实验课的作业

function[p]=cm_sm34(snr_in_dB)
%[p]=cm_sm8(snr_in_dB)
%     CM_SM8求出以dB为单位的给定信噪比的比特误码率和符号误码率
N=5000;
E=1;                      %每符号能量
snr=10^(snr_in_dB/10);    %信噪比
sgma=sqrt(E/(4*snr));     %噪声方差
%随后为数据源的产生
for i=1:3*N,
    temp=rand;            %在区间(0,1)间的一个均匀随机变量
    if (temp<0.5),
       dsource(i)=0;
    else
       dsource(i)=1; 
    end;
end;
%随后为数据源的差分编码
mapping=[0 1 3 2 7 6 4 5];
M=8;
[diff_enc_output]=cm_dpske(E,M,mapping,dsource);
%接收信号
for i=1:N
   [n(1) n(2)]=gngauss(sgma);
   r(i,:)=diff_enc_output(i,:)+n;
end;
%判决、误码率的计算
numoferr=0;
prev_theta=0;
for i=1:N,
  theta=angle(r(i,1)+j*r(i,2));
  delta_theta=mod(theta-prev_theta,2*pi);
  if ((delta_theta<pi/8) | (delta_theta>15*pi/8)),
    deris=[0 0 0];
  elseif (delta_theta<3*pi/8),
    deris=[0 0 1];
  elseif (delta_theta<5*pi/8),
    deris=[0 1 1];
  elseif (delta_theta<7*pi/8),
     deris=[0 1 0];
  elseif (delta_theta<9*pi/8),
    deris=[1 1 0];
  elseif (delta_theta<11*pi/8),
    deris=[1 1 1];
  elseif (delta_theta<13*pi/8),
     deris=[1 0 1];
  else 
     deris=[1 0 0];     
  end;
  prev_theta=theta;
  %如果判决不对，计错器加一
  if ((deris(1)~=dsource(3*i-2)) | (deris(2)~=dsource(3*i-1)) |(deris(3)~=dsource(3*i))),
     numoferr=numoferr+1;
  end;
end;
p=numoferr/N;