qam_64.m

64QAM 基于MATLAB调制解调仿真

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