function[p]=cm_sm8(snr_in_dB)
%[p]=cm_sm8(snr_in_dB)
%     CM_SM8求出以dB為單位的給定信噪比的比特誤碼率和符號誤碼率
N=1000;
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;