%qamdemod2.m
%16QAM解調程序
function[demodata]=qamdemod2(idata,qdata,para,nd,ml)
%******************variables*************************
%idata:I信道輸入數據
%qdata:Q信道輸入數據
%demodata:解調出的數據
%para:并行信道數
%nd:數據數
%ml:調制數
%(QPSK->2 16QAM->4)
%*****************************************************
k=sqrt(10);
idata=idata.*k;
qdata=qdata.*k;
demodata=zeros(para,ml*nd);

m2=ml/2;
count2=0;

for ii=1:nd
  for kk=1:para
    a=1; 
    b=1;
    i_lngth=0;
    q_lngth=0;
      for jj=1:m2
        if jj~=1
         if demodata(kk,jj-1+count2)==1
           a=-a;
      end
    if demodata(kk,m2+jj-1+count2)==1
      b=-b;
  end
  i_lngth=i_lngth+i_plrty.*2.^(m2-jj+1);
  q_lngth=q_lngth+q_plrty.*2.^(m2-jj+1);
end
if idata(kk,ii)>=i_lngth
   demodata(kk,jj+count2)=a>=0;
   i_plrty=1;
else
   demodata(kk,jj+count2)=a<=0;
   i_plrty=-1;
end
if qdata(kk,ii)>=q_lngth
   demodata(kk,m2+jj+count2)=b>=0;
   q_plrty=1;
else
   demodata(kk,m2+jj+count2)=b<=0;
   q_plrty=-1;
end
   end%for jj=1:m2
 end%kk=1:para
 count2=count2+ml;
end%for ii=1:nd
%********************end of file***************************