function[p]=cm_sm41(snr_in_dB) 
%[p]=cm_sm41(snr_in_dB) 
%        CM_SM41求出16 QAM以dB為單位的給定信噪比的 
%        符號誤碼率 
N=10000; 
d=1;                       %符號間的最小距離 
Eav=10*d^2;                %每符號能量 
snr=10^(snr_in_dB/10);     %信噪比 
sgma=sqrt(Eav/(3*snr));    %噪聲標準方差 
M=4; 
%數據源產生如下 
for i=1:N, 
	temp=rand;                    %在區間（0，1 )間的一個均勻隨機變量 
    dsource(i)=1+floor(M*temp);   %在1到64間的一個數，均勻的 
end; 
%隨后為信號星座的映射 
mapping=[-d d;
         d d;
        -d,-d;
        d,-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;