64QAM调制解调

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);