%qammod2.m
%16QAM調制程序
function[iout,qout] = qammod2(paradata,para,nd,m1);
%*************************variables********************************
%paradata:輸入數據
%iout
%qout
%para:并行信道數
%nd:OFDM數據數
%m1:調制數(4)      QPSK為2；16QAM為4
%*******************************************************************
%the constellation power
%m1 = 4;nd = 2;para = 4;paradata = [1 1 1 1 0 0 0 0;0 1 0 1 0 1 0 1;0 0 0 0 1 1 1 1;1 0 1 0 1 0 1 0];
%測試數據
k = sqrt(10);
iv = [-3 -1 3 1];

m2 = m1/2;
count2 = 0;
for ii = 1:nd
    
    isi = zeros(para,1);
    isq = zeros(para,1);
    
    for jj = 1:m1
        if jj <= m2;
            isi = isi + 2.^(m2 - jj).*paradata((1:para),count2 + jj);
        else
            isq = isq + 2.^(m1 - jj).*paradata((1:para),count2 + jj);
        end
    end
    
    iout((1:para),ii) = reshape(iv(isi + 1)./k,para,1);
    qout((1:para),ii) = reshape(iv(isi +1 )./k,para,1);
    
    count2 = count2 + m1;
end
%end of file**********************************************************************