%qammod2.m
%16QAM調制程序
function[iout,qout]=qammod2(paradata,para,nd,ml)
%******************variables*************************
%paradata:輸入數據
%iout:I信道輸出數據
%qout:Q信道輸出數據
%para:并行信道數
%nd:數據數
%ml:調制數
%(QPSK->2 16QAM->4)
%*****************************************************
%The constellation power
k=sqrt(10);
iv=[-3-1 3 1];
m2=ml/2;
count2=0;
for ii=1:nd
isi=zeros(para,1);
isq=zeros(para,1);
for jj=1:ml
if jj<=m2
  isi=isi+2.^(m2-jj).*paradata((1:para),count2+jj);
else
  isq=isq+2.^(ml-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(isq+1)./k,para,1);
count2=count2+ml;
end
%********************end of file***************************