qammod2.m

信道编码采用卷积码

% qammod2.m
% 16QAM调制程序

function [iout,qout]=qammod2(paradata,para,nd,ml)

%****************** variables *************************
% paradata : 输入数据 (para-by-nd matrix)
% iout     : Ich路输出数据
% qout     : Qch路输出数据
% 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=ml
       if jj<=m2
           isi=isi+2.^(m2-jj).*paradata((l:para),count2+jj);
       else
           isq=isq+2.^(ml-jj).*paradata((l:para),count2+jj);
       end
   end
   iout((l:para),ii)=reshape(iv(isi+1)./k,para,1);
   qout((l:para),ii)=reshape(iv(isq+1)./k,para,1);
   count2=count2+ml;
end
%********************end of file********************