QPSK和4QAM代码(2)

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