qam_8.m

数字调制解调的MATLAB仿

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