function[p]=cm_sm41(snr_in_dB)
%[p]=cm_sm41(snr_in_dB)
%        CM_SM41求出64 QAM以dB為單位的給定信噪比的
%        符號(hào)誤碼率
N=10000;
d=1;                       %符號(hào)間的最小距離
Eav=10*d^2;                %每符號(hào)能量
snr=10^(snr_in_dB/10);     %信噪比
sgma=sqrt(42/(16*snr));    %噪聲標(biāo)準(zhǔn)方差
M=64;
%數(shù)據(jù)源產(chǎn)生如下
for i=1:N,
	temp=rand;                    %在區(qū)間（0，1 )間的一個(gè)均勻隨機(jī)變量
    dsource(i)=1+floor(M*temp);   %在1到64間的一個(gè)數(shù)，均勻的
end;
%隨后為信號(hào)星座的映射
mapping=[ -7*d 7*d ; -5*d 7*d;  -3*d 7*d;  -1*d  7*d;  1*d 7*d;  3*d 7*d;  5*d 7*d;  7*d 7*d;       
          -7*d 5*d ; -5*d 5*d;  -3*d 5*d;  -1*d  5*d;  1*d 5*d;  3*d 5*d;  5*d 5*d;  7*d 5*d;         
          -7*d 3*d ; -5*d 3*d;  -3*d 3*d;  -1*d  3*d;  1*d 3*d;  3*d 3*d;  5*d 3*d;  7*d 3*d;         
          -7*d 1*d ; -5*d 1*d;  -3*d 1*d;  -1*d  1*d;  1*d 1*d;  3*d 1*d;  5*d 1*d;  7*d 1*d;          
          -7*d -1*d; -5*d -1*d; -3*d -1*d; -1*d  -1*d; 1*d -1*d; 3*d -1*d; 5*d -1*d; 7*d -1*d;
          -7*d -3*d; -5*d -3*d; -3*d -3*d; -1*d  -3*d; 1*d -3*d; 3*d -3*d; 5*d -3*d; 7*d -3*d;
          -7*d -5*d; -5*d -5*d; -3*d -5*d; -1*d  -5*d; 1*d -5*d; 3*d -5*d; 5*d -5*d; 7*d -5*d;
          -7*d -7*d; -5*d -7*d; -3*d -7*d; -1*d  -7*d; 1*d -7*d; 3*d -7*d; 5*d -7*d; 7*d -7*d];
for i=1:N,
    qam_sig(i,:)=mapping(dsource(i),:);
end;
%接收信號(hào)
for i=1:N,
    [n(1) n(2)]=gngauss(sgma);
    r(i,:)=qam_sig(i,:)+n;
end;
%判決、錯(cuò)誤概率的計(jì)算
numoferr=0;
for i=1:N,
    %隨后為量度計(jì)算
    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;