function[p]=cm_sm41(snr_in_db)
%[p]=cm_sm41(snr_in_db)
% 求出以DB為單位的給定信噪比的比特誤碼率和符號誤碼率
N=10000;
d=1; %符號間最小距離
eav=10*d^2; %每符號能量
snr=10^(snr_in_db/10); %每比特信噪比（給定）
sgma=sqrt(eav/(8*snr)); %噪聲方差
M=16;
%產(chǎn)生原始數(shù)據(jù)序列
for i=1:N,
    temp=rand; %在區(qū)間（0，1）間的均勻瑞利變量
    dsource(i)=1+floor(M*temp); %在1到16間的一個數(shù)，均勻的
end;
%為信號星座的映射
mapping=[-3*d 3*d; -d 3*d; d 3*d; 3*d 3*d;
         -3*d d; -d d; d d; 3*d d;
         -3*d -d; -d -d; d -d; 3*d -d;
         -3*d -3*d; -d -3*d; d -3*d; 3*d -3*d];
     for i=1:N,
         qam_sig(i,:)=mapping(dsource(i),:);
     end;
     %接收信號
     for i=1:N,
         [n(1) n(2)]=gngauss(sgma);
         r(i,:)=qam_sig(i,:)+n;
     end;
     %判決、誤碼率的計算
     numoferr=0;
     for i=1:N,
         %量度計算
         for j=1:M,
             metrics(j)=(r(i,1)-mapping(j,1))^2+(r(i,2)-mapping(j,2))^2;
         end;
     [min_metric decis]=min(metrics);
     if (decis~=dsource(i)),
         numoferr=numoferr+1;
     end;
     end;
     p=numoferr/(N);