% File name:16QAM.m
clear all;
clc;
%/**************产生信源***************/
N=10000;
data=randint(1,4*N);
data_re=zeros(4,N);
Eav=10;
data_re=reshape(data,4,N);
%/**************符号映射***************/
% I_data=zeros(1,2*N);
% Q_data=zeros(1,2*N);
% I_data(1,:)=data(1:2:end);
% Q_data(1,:)=data(2:2:end);
data_mod=zeros(2,N);
for k=1:N
if (data_re(:,k)==[0 0 0 0]')
data_mod(:,k)=[-1 -1]';
elseif (data_re(:,k)==[0 0 0 1]')
data_mod(:,k)=[-1 -3]';
elseif (data_re(:,k)==[0 0 1 0