qam_4.m

来自「数字调制解调的MATLAB仿」· M 代码 · 共 37 行

function [p]=QAM_4(snr_in_dB)
N=10000;
d=1;				  	% 符号间最小距离
Eav=2*d^2;		 	  	% 符号能量
snr=10^(snr_in_dB/10);	 	  	% 每比特符号的信噪比
sgma=sqrt(Eav/(4*snr));	  	  	% 噪声方差
M=4;
% 产生数据
for i=1:N,	
  temp=rand;			  
  dsource(i)=1+floor(M*temp);	  	% 生成[1，4]的随机数
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);