基于matlab的16QAM调制解调仿真

% 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]')

        data_mod(:,k)=[-1 1]';

    elseif (data_re(:,k)==[0 0 1 1]')

        data_mod(:,k)=[-1 3]';

    elseif (data_re(:,k)=