QPSK和16QAM调制matlab编码

function [outputinput] = data_produce（NcarrierN_symPerframeModulation）
    %子载波数、每帧中的ofdm符号数、调制方式  
    % 1--BPSK调制 2--QPSK调制， 4--16QAM调制
 Es = 1;                 % 在16QAM 64QAM调制方式下符号能量都被归一化 

 N_data_frame = Ncarrier * N_symPerframe;  % number of data per OFDM frame （实际产生数据部分的比特数，子载波数*一帧中的ofdm的符号数）
 
 data = randint（ModulationN_data_frame[01]）;
 input = zeros（1Modulation*N_data_frame）;
 output = zeros（NcarrierN_symPerframe）;
 
 input = reshape（data1Modulation*N_data_frame）;
 
 switch Modulation
     case 2                  %QPSK调制
         kmod = 1/sqrt（2）;
         datatmp = data.‘;
         datatmptmp = bi2de（datatmp‘left-msb‘）;
         % 00->1+i10->-1+i01->1-i11->-1-i
         QAMTable = [1+i1-i-1+i-1-i];
         qpsk_data = QAMTable